专利名称：海岛遥感测图高程精度控制与修正方法及系统的制作方法
技术领域：
本发明涉及大地测量工程技术领域，特别涉及一种海岛遥感测图高程精度控制与修正方法及系统。
背景技术：
我国海域广阔，海岛(礁)分布零散、极不均匀。各类海岛(礁)12000多个，海岛(礁)是国家领土的重要组成部分，其精确位置和基础地理信息是划分领海及其他管辖领域的重要依据；全面准确地掌握我国海岛(礁)位置和地理信息，是维护国家主权、保障国家安全、实施海洋开发的重要保障，具有重要的现实意义和深远的战略意义。随着遥感影像空间和光谱分辨率的大幅提高及雷达影像的丰富，海岛(礁)识别特别是暗礁探测的数据条件得到改善，但我国综合用于海岛(礁)遥感测图的高程精度控制技术手段缺乏，海岛(礁)遥感测图的高程精度控制技术需要进一步的研究。目前，在陆地上主要运用GPS测量技术、CQG2000似大地水准面模型的内插法、水准高程测量等方法提供像片控制点的平面、高程成果，极大的减轻在我国困难地区测图的工作强度，提高的了工作效率，取得了明显的经济效益。由于我国海岛海岛礁分布零散，选择有效的高程控制点存在很大的困难。在海上无法有效的布设GPS控制点，也无法进行水准高程测量。同时，由于CQG2000似大地水准面模型受地形条件的限制，以及地面观测数据的密度和不均匀分布，致使CQG2000似大地水准面的精度在全国范围内是不一致的。总体分辨率为5' X5,，全国平均精度理论上为±0.36米，不足以满足大比例尺的遥感测图需求。除少量较大面积的海岛外，海岛(礁)的高程一般不大，在我国，高程超过IOm的海岛(礁)数量占全部海岛(礁)的比例不到15%。由于海岛被海水包围，其高程有海平面作绝对参考，大的高程误差甚至会导致海平面与海岛的拓扑关系出现错误。因此，在海岛测图中，海岛高程的准确性和精度至关重要。在陆地航空航天遥感测绘中，当地面像控点稀少时，可以采用基于POS和影像纹理匹配的稀少控制摄影测量方法进行立体测图。但值得注意的是，由于目前摄影测量工作站的立体测图软件普遍使用正常高系统作为像控点高程，当地面像控点稀少时，测图的高程会产生形变，其形变程度由像控点分布密度和大地水准面起伏两个因素共同决定。当像控点间距大于200km时，大地水准面起伏在我国陆地部分地区可达到20m，如果直接采用稀少控制的航空航天遥感测图方法，虽然不会引起地形图要素的拓扑关系改变，但可导致绝对高程的最大误差达到10m。与陆地相比，海岛(礁)分布稀疏，无法按常规的密度和点位布测像控点，海域绝大部分被海水覆盖，影像地面特征稀少、纹理少使得航空航天影像匹配的难度大，导致模型连接困难甚至无法实现。采用稀少控制的航空摄影测量方法很可能由于大地水准面起伏大导致海岛地形图高程 不准确。我们知道，与陆地不同，海岛的高程注记不能因为比例尺小而取舍，海岛岸线高程精度应在0.3m左右比较合适。目前陆地地形图中有关高程精度指标普遍较低，用于指导海岛地形测图甚至可能会导致海岛掉到海水面以下的矛盾现象。
发明内容
(一 )要解决的技术问题本发明要解决的技术问题是:如何提高海岛遥感测图所获得的高程精度。(二)技术方案为解决上述技术问题，本发明提供了一种海岛遥感测图高程精度控制与修正方法，包括以下步骤:S1:提取海岛遥感影像中的岸线；S2:通过海岛验潮站的潮位观测数据获得第一平均海面高、以及卫星测高数据获得第二平均海面高，将所述第二平均海面高进行基准转换，以保证与第一平均海面高的基准一致，并将获得的第一平均海面高和基准转换后的第二平均海面高作为标准平均海面闻;S3:通过地面重力观测数据获得重力数据，通过卫星测高数据获得重力异常数据，将重力异常数据进行基准转换，以保证与所述重力数据的基准一致，并根据所述重力数据、以及基准转换后的重力异常数据来计算大地水准面；S4:对所述海岛验潮站的潮位观测数据进行调和分析，以获得分潮调和常数；S5:根据所述分潮调和常数计算获得平均大潮高潮面与深度基准面之间的高度差、以及所述标准平均海面高 与深度基准面之间的高度差，并根据所述平均大潮高潮面与深度基准面之间的高度差、所述标准平均海面高与深度基准面之间的高度差、以及所述标准平均海面高与大地水准面之间的高度差来计算所述平均大潮高潮面的高程；S6:根据所述平均大潮高潮面的高程来确定所述平均大潮高潮面上的岸线，通过所述海岛遥感影像中的岸线、以及所述平均大潮高潮面上的岸线来计算并去除所述海岛遥感影像中的岸线的系统性误差；S7:将去除系统性误差后的所述海岛遥感影像中的岸线与所述平均大潮高潮面上的岸线进行拟合，以实现对所述海岛遥感影像中的岸线的修正。优选地,所述海岛遥感影像包括:RapidEye_one、Alos、Spot5、IK0N0S、QuickBird、WorIdView-one、WorIdView-Two、GeoEye-One>以及低空无人机所获得的海岛遥感数据中的
至少一种。优选地,步骤S2具体包括以下步骤:S21:对海岛验潮站的潮位观测数据计算平均值，以获得第一平均海面高；S22:根据卫星测高数据中的观测海面高度计算平均值，以获得第二平均海面高；S23:对所述第二平均海面高进行基准转换，以保证与第一平均海面高的基准一
致，在所述第二平均海面高进行基准转换时，采用下列公式，
dh = -Wda + β η2φ /其中，a为第一平均海面高对应的第一参考椭球的长半轴的长度，f为所述第一参考椭球的扁率，0为大地讳度，
权利要求
1.一种海岛遥感测图高程精度控制与修正方法，其特征在于，包括以下步骤: 51:提取海岛遥感影像中的岸线； 52:通过海岛验潮站的潮位观测数据获得第一平均海面高、以及卫星测高数据获得第二平均海面高，将所述第二平均海面高进行基准转换，以保证与第一平均海面高的基准一致，并将获得的第一平均海面高和基准转换后的第二平均海面高作为标准平均海面高； 53:通过地面重力观测数据获得重力数据，通过卫星测高数据获得重力异常数据，将重力异常数据进行基准转换，以保证与所述重力数据的基准一致，并根据所述重力数据、以及基准转换后的重力异常数据来计算大地水准面； 54:对所述海岛验潮站的潮位观测数据进行调和分析，以获得分潮调和常数； 55:根据所述分潮调和常数计算获得平均大潮高潮面与深度基准面之间的高度差、以及所述标准平均海面高与深度基准面之间的高度差，并根据所述平均大潮高潮面与深度基准面之间的高度差、所述标准平均海面高与深度基准面之间的高度差、以及所述标准平均海面高与大地水准面之间的高度差来计算所述平均大潮高潮面的高程； 56:根据所述平均大潮高潮面的高程来确定所述平均大潮高潮面上的岸线，通过所述海岛遥感影像中的岸线、以及所述平均大潮高潮面上的岸线来计算并去除所述海岛遥感影像中的岸线的系统性误差； S7:将去除系统性误差后的所述海岛遥感影像中的岸线与所述平均大潮高潮面上的岸线进行拟合，以实现对所述海岛遥感影像中的岸线的修正。
2.如权利要求1所述的方法，其特征在于，所述海岛遥感影像包括:RapidEye_one、Alos、Spot5、IKONOS、QuickBird、WorIdView-one、WorIdView-Two、GeoEye-One、以及低空无人机所获得的海岛遥感数据中的至少一种。·
3.如权利要求1所述的方法，其特征在于，步骤S2具体包括以下步骤: 521:对海岛验潮站的潮位观测数据计算平均值，以获得第一平均海面高； 522:根据卫星测高数据中的观测海面高度计算平均值，以获得第二平均海面高； 523:对所述第二平均海面高进行基准转换，以保证与第一平均海面高的基准一致，在所述第二平均海面高进行基准转换时，采用下列公式，dh = -Wda + β η2φ / 其中，a为第一平均海面高对应的第一参考椭球的长半轴的长度，f为所述第一参考椭球的扁率，P为大地纬度，F = Jl-C2 sin2^，e为所述第一参考椭球的第一偏心率，da =a0-a,df = f0-f, a0为第二平均海面高对应的第二参考椭球的长半轴的长度，为第二参考椭球的扁率，da为长半轴改正值，df为扁率改正值，dh为参考椭球转换引起的高程变化； 524:将获得的第一平均海面高和基准转换后的第二平均海面高作为标准平均海面高。
4.如权利要求1所述的方法，其特征在于，步骤S3具体包括以下步骤: 531:通过地面重力观测数据获得重力数据； 532:通过卫星测高数据获得重力异常数据，将重力异常数据进行基准转换，以保证与所述重力数据的基准一致，在所述重力异常数据进行基准转换时，采用下列公式，dg = -Wda + β η2 φ /其中，a为所述重力数据对应的第一参考椭球的长半轴的长度，f为所述第一参考椭球的扁率，(^为大地纬度
5.如权利要求1所述的方法，其特征在于，步骤S4中，调和分析的公式为，
6.如权利要求1所述的方法，其特征在于，步骤S5中，通过下式来计算所述平均大潮高潮面的高程， L = L3- (Lg-L1) 其中，L为所述平均大潮高潮面的高程，L3为所述平均大潮高潮面与深度基准面之间的高度差，L2为所述标准平均海面高与深度基准面之间的高度差，L1为所述标准平均海面高与大地水准面之间的高度差。
7.如权利要求1所述的方法，其特征在于，步骤S6中，通过下式计算所述系统性误差，
8.如权利要求1所述的方法，其特征在于，步骤S7中，设所述海岛遥感影像中的岸线与所述平均大潮高潮面上的岸线之间任一点的偏差为Z(x，y)且Z(x，y)=aiX+aj+apy+ay+a；^2,其中，a1; a2, a3, a4, a5为常数，x, y分别为横轴坐标和纵轴坐标，对所述偏差按照最小二乘法，获得a1; a2, a3, a4, a5的值,并根据获得了 a1; a2, a3, a4, a5的值的偏差公式实现对所述海岛遥感影像中的岸线的修正。
9.一种海岛遥感测图高程精度控制与修正系统，其特征在于，包括: 岸线提取模块，用于提取海岛遥感影像中的岸线； 平均海面计算模块，用于通过海岛验潮站的潮位观测数据获得第一平均海面高、以及卫星测高数据获得第二平均海面高，将所述第二平均海面高进行基准转换，以保证与第一平均海面高的基准一致，并将获得的第一平均海面高和基准转换后的第二平均海面高作为标准平均海面高；大地水准面计算模块，用于通过地面重力观测数据获得重力数据，通过卫星测高数据获得重力异常数据，将重力异常数据进行基准转换，以保证与所述重力数据的基准一致，并根据所述重力数据、以及基准转换后的重力异常数据来计算大地水准面； 调和分析模块，用于对所述海岛验潮站的潮位观测数据进行调和分析，以获得分潮调和常数； 高潮面计算模块，用于根据所述分潮调和常数计算获得平均大潮高潮面与深度基准面之间的高度差、以及所述标准平均海面高与深度基准面之间的高度差，并根据所述平均大潮高潮面与深度基准面之间的高度差、所述标准平均海面高与深度基准面之间的高度差、以及所述标准平均海面高与大地水准面之间的高度差来计算所述平均大潮高潮面的高程; 系统性误差计算模块，用于根据所述平均大潮高潮面的高程来确定所述平均大潮高潮面上的岸线，通过所述海岛遥感影像中的岸线、以及所述平均大潮高潮面上的岸线来计算并去除所述海岛遥感影像中的岸线的系统性误差； 拟合修正模块，用于将去除系统性误差后的所述海岛遥感影像中的岸线与所述平均大潮高潮面上的岸线进行 拟合，以实现对所述海岛遥感影像中的岸线的修正。
全文摘要
本发明公开了一种海岛遥感测图高程精度控制与修正方法及系统，涉及大地测量工程技术领域，该方法包括S1提取海岛遥感影像中的岸线；S2获得第一平均海面高、以及第二平均海面高，将所述第二平均海面高进行基准转换；S3计算大地水准面；S4进行调和分析获得分潮调和常数；S5计算所述平均大潮高潮面的高程；S6计算并去除所述海岛遥感影像中的岸线的系统性误差；S7将去除系统性误差后的所述海岛遥感影像中的岸线与所述平均大潮高潮面上的岸线进行拟合。本发明通过对海岛遥感影像中的岸线与所述平均大潮高潮面上的岸线进行拟合，来实现对所述海岛遥感影像中的岸线的修正，提高了海岛遥感测图所获得的高程精度。
文档编号G01C5/00GK103245324SQ20121002506
公开日2013年8月14日 申请日期2012年2月6日 优先权日2012年2月6日
发明者柯宝贵, 章传银, 张利明, 李德海 申请人:中国测绘科学研究院