State-to-state rotational relaxation of carbon monoxide (CO) has been studied using an ir double resonance technique. Individual rotational lines of the (2-0) vibrational overtone band were pumped by a pulsed tunable ir laser and the subsequent evolution of the v=2 rotational population distribution was monitored by the absorption of a tunable cw ir laser via the (3-2) band transitions. Both the excitation and probe lasers were linearly polarized, with linewidths that were narrower than the CO Doppler width. Consequently, alignment and velocity relaxation effects were observed in these measurements. A data set consisting of 54 time-dependent rotational state population profiles was acquired. The full CO-CO rotational relaxation matrix, which consists of state-to-state rate constants for rotational levels up to J=29, was deduced from computer simulations of the data. Scaling and fitting laws were used to provide parametric representations of the rate constants. The three most common models, modified exponential gap, statistical polynomial-exponential gap (SPEG), and energy corrected sudden with exponential-power gap (ECS-EP) were investigated. We concluded that the SPEG law best reproduced the CO-CO rotational energy transfer data. A propensity to preserve the CO parity in rotational energy transfer was observed for collisions where the amount of energy transferred was small. Hence even DeltaJ processes were favored for transitions between levels with low J values. For near-single collisions events a correlation was noted between the amount of rotational energy transferred and the degree of velocity distribution relaxation. This correlation yielded insights regarding the energy transfer dynamics.