[发明专利]一种同轨激光测高点与三线阵立体影像联合测绘处理方法

权利要求书

1.一种同轨激光测高点与三线阵立体影像联合测绘处理方法，其特征在于，所述方法包括以下步骤：

步骤1以卫星成像摄影轨道为单位，将同一轨多景三线阵立体影像构建一个平差区域网，并采用基于像方仿射变换补偿有理函数模型RFM的区域网平差模型；

步骤2在各三线阵立体影像内部布设连接点，在相邻的三线阵立体影像之间布设公共连接点；

步骤3根据激光测高点的经纬度坐标和高程值，通过正视影像的RFM，计算获取激光测高点在同轨获取的正视影像上的像点坐标；

步骤4分别将前视和后视影像与正视影像开展高精度的同名点匹配，获取正视影像上的激光测高点像点在前视和后视影像上的同名点，即获取了激光测高点在前视影像和后视影像上的像点坐标；

步骤5针对每个激光测高点像点，将影像RFM的像方仿射变换参数和激光测高点像点对应的地面经纬度坐标作为未知数，将激光测高点的高程值作为已知量，逐点构建误差方程；针对每个连接点，将影像RFM的像方仿射变换参数和连接点对应的地面三维坐标作为未知数，逐点构建误差方程；

步骤6对误差方程进行法化，得到法方程；采用最小二乘原理，进行整体平差解算，解算每张影像RFM的仿射变换参数、连接点对应的地面三维坐标改正数和激光测高像点对应的地面经纬度坐标改正数；

步骤7每完成一次平差解算，判断平差是否满足收敛条件，若满足，则平差结束；反之，利用平差解算获取的改正数更新影像RFM、连接点对应地面三维坐标和激光测高像点对应的地面经纬度坐标，重复步骤5-步骤6，再次进行平差迭代解算，不断更新影像RFM的仿射变换参数、连接点对应的地面三维坐标和激光测高像点对应的地面经纬度坐标，直到满足收敛条件为止；最后输出每张影像RFM的仿射变换参数；

步骤8利用最后一次平差迭代计算得到的每张影像RFM的仿射变换参数,对RFM模型进行像方仿射变换补偿，为每张影像生成得到新的RFM；

所述步骤4中，分别匹配正视影像上的激光测高点像点在前视和后视影像上的同名点具体包括以下步骤：

步骤4.1计算正视影像上激光测高点像点坐标位置的尺度不变特征变换SIFT特征，并保存特征描述符；

步骤4.2根据激光测高点的地面经纬度坐标和高程值，利用前视影像的RFM，计算激光测高点在前视影像上的大致像点位置，并以此为中心，向四周外扩一定像素范围作为匹配的搜索范围；计算前视影像搜索范围内每个像素的SIFT特征，并保存其特征描述符；

步骤4.3使用快速最近邻逼近搜索函数库搜索步骤4.2获取的多个特征描述符中与正视影像上激光测高点像点的SIFT特征描述符最近邻值，获取该特征描述符对应的前视影像上的像点坐标，亦即前视影像上激光测高点像点的整像素坐标值；

步骤4.4以前视影像上激光测高点像点的整像素坐标位置为中心，与正视影像上激光测高点像点开展最小二乘匹配，获取前视影像上激光测高点像点的精确坐标值；

步骤4.5重复执行步骤4.2-步骤4.4，同样获取后视影像上激光测高点像点的精确坐标值。