Currently, Mars Express HRSC imagery is an essential data source to derive high accuracy Mars topographic data. In view of the characteristics of Martian surface satellite imagery, a dense image matching scheme for HRSC imagery based on precise point prediction method is proposed. The image matching strategies of our method are elaborated in detail. Based on the proposed method, DEM and DOM of Martian surface are derived and compared with those published by ESA. The experiment results show that the root mean square error in planar direction is about three pixels, while the root mean square error in height direction is about one pixel. Moreover, the mean square error in plane direction show a certain systematic error and the reasons are analysed. Experiment results also demonstrate that the point prediction accuracy for corresponding points is up to 1-3 pixels.