  The propagation of polarized emission in pulsar magnetosphere is investigated
in this paper. The polarized waves are generated through curvature radiation
from the relativistic particles streaming along curved magnetic field lines and
co-rotating with the pulsar magnetosphere. Within the 1/{\deg} emission cone,
the waves can be divided into two natural wave mode components, the ordinary
(O) mode and the extraord nary (X) mode, with comparable intensities. Both
components propagate separately in magnetosphere, and are aligned within the
cone by adiabatic walking. The refraction of O-mode makes the two components
separated and incoherent. The detectable emission at a given height and a given
rotation phase consists of incoherent X-mode and O-mode components coming from
discrete emission regions. For four particle-density models in the form of
uniformity, cone, core and patches, we calculate the intensities for each mode
numerically within the entire pulsar beam. If the co-rotation of relativistic
particles with magnetosphere is not considered, the intensity distributions for
the X-mode and O-mode components are quite similar within the pulsar beam,
which causes serious depolarization. However, if the co-rotation of
relativistic particles is considered, the intensity distributions of the two
modes are very different, and the net polarization of out-coming emission
should be significant. Our numerical results are compared with observations,
and can naturally explain the orthogonal polarization modes of some pulsars.
Strong linear polarizations of some parts of pulsar profile can be reproduced
by curvature radiation and subsequent propagation effect.
