This paper presents bilinear and bicubic interpolation methods tailored for the division of focal plane polarization imaging sensor. The interpolation methods are targeted for a 1-Mega pixel polarization imaging sensor operating in the visible spectrum. The five interpolation methods considered in this paper are: bilinear, weighted bilinear, bicubic spline, an approximated bicubic spline and a bicubic interpolation method. The modulation transfer function analysis is applied to the different interpolation methods, and test images as well as numerical error analyses are also presented. Based on the comparison results, the full frame bicubic spline interpolation achieves the best performance for polarization images.