Imaging polarimetry is a powerful technique to study circumstellar dusty environments. Detailed analyses of polarization disks allow to quantitatively evaluate the disk geometries and the dust properties. We present near-infrared (NIR) polarimetric images and radiative transfer modeling of the high-mass protostar CRL 2136 and the Herbig Be star R Mon. In CRL 2136, the data show a central star feature (IRS 1), a spot-like high polarization region towards IRS 1 (PK = 32%), and a linear polarization vector alignment along the equatorial direction around IRS 1. In our spherical grain model, the optical depth of the disk is estimated to be 210 in the V-band. With such an intermediate optical depth, the linear vector alignment and the visible central star feature can not be reproduced simultaneously. We interpret that the linear vector alignment detected in our observation is probably due to dichroism by aligned nonspherical grains. In R Mon, our NIR polarimetric data detected a butterfly-shaped polarization di...