Abstract: According to the character of the CCD sensors which set up on HJ-1A/1B satellites,a new imaging model which is extended from the collinear equation for linear sensor is built.In this model,quadratic polynomial and Hermite interpolation polynomial are employed for fitting the variation of the exterior orientation elements,and high-order polynomial is used as distortion model,plus engaging iteration method to solve the frequent problem,ill-conditioned system.Furthermore,we also present a method based on bipartite graph matching for correspondence points matching to meet the requirement of the distribution of control-points.The empirical results show that the proposed model compared with existing models such as quadratic polynomial and rational function model is able to achieve a higher accuracy,and requires for smaller condition number,that is sufficient for keeping pace with the demand of the research for global land cover.