本发明涉及属于遥感影像的地形重建技术方法领域，尤其是一种月球表面高分辨率dem提取方法，目的在于生成高分辨率、高精度月表dem。
背景技术：
月表的三维建模通常以一个数字高程模型(dem)来显示。利用dem进行地形描述有其特有的优势：如容易以多种形式显示地形信息，并且其精度不会损失，在需要对其地理信息进行添加或修改等操作时，容易实现自动化、实时化。月表dem是选择登月着陆点，勘查月球资源，研究月球撞击坑的三维特征等月球探测研究活动不可缺少的基础信息，也是构建基于g/s模式的三维数字月球平台，实现月球三维可视化的重要数据。
随着中国月球探测工程的开展，国内也掀起了月球数字高程模型研究的热潮。国内使用嫦娥工程数据制作全月表dem的方法大概可以分为2类：⑴直接利用嫦娥激光高度计测高数据制作月球表面模型，利用嫦娥一号卫星搭载的激光高度计在第一次正飞阶段所获取的有效数据点，采用最小曲率法、层次多结点样条算法、kriging差值等方法生成月球高程模型。由于卫星本身的不稳定性、激光测高计的系统噪声、月球表面的尘埃干扰等因素的影响，激光测高数据本身存在着粗差点，这些粗差点的存在表现为高程异常，使得构建的数字高程模型分辨率、精度较低。⑵根据嫦娥卫星ccd三线阵相机传回的影像数据与星载激光高度计结合得到全月高程数据。首先使用嫦娥卫星搭载的激光高度计所获得的数据作为月面控制点，在此基础上，采用摄影测量原理进行部分外方位元素值的计算，然后通过最小二乘拟合出所有的外方位元素，实现ccd影像每行外方位元素值的计算。然后根据前方交会的原理计算得到离散的月面点坐标，最后选择内插方法对离散的数据点进行内插得到月球dem。该方法建立的月球dem与直接使用激光高度计数据制作的月球dem相比精度有所提高，但是由于月球影像纹理特征不明显，而且其影像灰度受光照影响，导致提取的控制点有一定的误差，这就直接影响了所生成的dem的精度，其中利用嫦娥一号立体影像提取的月表dem分辨率只能达到百米级，对嫦娥二号传回的立体影像进行处理，理论上可以得到分辨率优于10米的月表dem，但仍不能满足我们对月球进行更细致的观察研究。
技术实现要素：
为了克服已有月球表面dem提取方法的精度较低、分辨率较低的不足，本发明提供一种高精度高分辨率的月球表面dem提取方法。
本发明解决其技术问题所采用的技术方案是：
一种月球表面高分辨率dem提取方法，包括以下步骤：
步骤1)获取月球影像数据，包括ce-2ccd影像数据、ce-1激光高度计lam数据和ce-3全景相机pcam数据；
步骤2)对步骤1)中的ccd影像前视、后视数据和pcam数据进行预处理；
步骤3)提取pcam左右相机的同名像点以及外方位元素；
步骤4)使用摄影测量学中的前方交会方法计算pcam影像同名像点的物方坐标m1；
步骤5)提取ccd影像数据的同名像点与lam数据中的控制点坐标；
步骤6)选取足量的虹湾地区lam数据控制点通过有理函数模型计算影像成像时的外方位元素；
步骤7)使用摄影测量学中前方交会方法计算ccd影像数据虹湾地区的物方坐标m2；
步骤8)统一m1与m2坐标系，计算m1与m2的误差，如果m1和m2的误差大于指定阈值θ，说明有理函数模型系数的定位精度有计算误差，返回步骤6)进行下一次迭代，否则，迭代结束；
步骤9)矫正全月表同名像点的外方位元素，计算ccd影像数据的物方坐标m3；
步骤10)使用m3制作规则的全月dem。
进一步，所述步骤2)中，预处理包括分析选择影像数据以及影像格式转换。
本发明的有益效果主要表现在：生成高分辨率、高精度月表dem。
具体实施方式
下面对本发明作进一步描述。
一种月球表面高分辨率dem提取方法，包括以下步骤：
步骤1)获取月球影像数据，包括ce-2ccd影像数据、ce-1激光高度计lam数据和ce-3全景相机pcam数据；
步骤2)对步骤1)中的ccd影像数据和pcam数据进行预处理(分析选择影像数据以及影像格式转换等)；
步骤3)提取pcam左右相机的同名像点以及外方位元素；
步骤4)使用摄影测量学中的前方交会方法计算pcam影像同名像点的物方坐标m1；
步骤5)提取ccd影像数据的同名像点与lam数据中的控制点坐标；
步骤6)选取足量的虹湾地区lam数据控制点通过有理函数模型计算影像成像时的外方位元素；
步骤7)使用摄影测量学中前方交会方法计算ccd影像数据虹湾地区的物方坐标m2；
步骤8)统一m1与m2坐标系，计算m1与m2的误差，如果m1和m2的误差大于指定阈值θ，说明ccd相机外方位元素有计算误差，返回步骤7)进行下一次迭代，否则，迭代结束；
步骤9)矫正全月表同名像点的外方位元素，计算ccd影像数据的物方坐标m3；
步骤10)使用m3制作规则的全月dem。
建立有理函数模型的前提条件是选取足够数量的月面控制点，对于ccd线阵推扫式成像模型建立三阶有理函数模型，而且为了避免法方程的病态问题选取均匀分布于影像上的点作为控制点，利用相机焦距、像元尺寸、传感器角度这三个基本参数，解算出有理函数模型系数及外方位元素。
通过有理函数模型求解像片外方位元素xs、ys、zs、
其中，摄影基线b的三个坐标分量(bx,by,bz)可由外方位元素计算：
计算像点a(x,y)和a'(x',y')的空间辅助坐标a1(x1,y1,z1)和a2(x2,y2,z2)。
其中，旋转矩阵r1，r2，表示为：
由式(1)、(2)计算投影系数n1、n2：
求解待定点的物方坐标a(xa,ya,za)如下：
所以有效提取ccd影像的控制点，对相机外方位元素解算以及对全月表dem精度都有极大影响。