CN103245324A - 海岛遥感测图高程精度控制与修正方法及系统 - Google Patents
海岛遥感测图高程精度控制与修正方法及系统 Download PDFInfo
- Publication number
- CN103245324A CN103245324A CN2012100250629A CN201210025062A CN103245324A CN 103245324 A CN103245324 A CN 103245324A CN 2012100250629 A CN2012100250629 A CN 2012100250629A CN 201210025062 A CN201210025062 A CN 201210025062A CN 103245324 A CN103245324 A CN 103245324A
- Authority
- CN
- China
- Prior art keywords
- sea level
- mean
- island
- height
- data
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Images
Landscapes
- Image Processing (AREA)
Abstract
本发明公开了一种海岛遥感测图高程精度控制与修正方法及系统,涉及大地测量工程技术领域,该方法包括:S1:提取海岛遥感影像中的岸线;S2:获得第一平均海面高、以及第二平均海面高,将所述第二平均海面高进行基准转换;S3:计算大地水准面;S4:进行调和分析获得分潮调和常数;S5:计算所述平均大潮高潮面的高程;S6:计算并去除所述海岛遥感影像中的岸线的系统性误差;S7:将去除系统性误差后的所述海岛遥感影像中的岸线与所述平均大潮高潮面上的岸线进行拟合。本发明通过对海岛遥感影像中的岸线与所述平均大潮高潮面上的岸线进行拟合,来实现对所述海岛遥感影像中的岸线的修正,提高了海岛遥感测图所获得的高程精度。
Description
技术领域
本发明涉及大地测量工程技术领域,特别涉及一种海岛遥感测图高程精度控制与修正方法及系统。
背景技术
我国海域广阔,海岛(礁)分布零散、极不均匀。各类海岛(礁)12000多个,海岛(礁)是国家领土的重要组成部分,其精确位置和基础地理信息是划分领海及其他管辖领域的重要依据;全面准确地掌握我国海岛(礁)位置和地理信息,是维护国家主权、保障国家安全、实施海洋开发的重要保障,具有重要的现实意义和深远的战略意义。随着遥感影像空间和光谱分辨率的大幅提高及雷达影像的丰富,海岛(礁)识别特别是暗礁探测的数据条件得到改善,但我国综合用于海岛(礁)遥感测图的高程精度控制技术手段缺乏,海岛(礁)遥感测图的高程精度控制技术需要进一步的研究。
目前,在陆地上主要运用GPS测量技术、CQG2000似大地水准面模型的内插法、水准高程测量等方法提供像片控制点的平面、高程成果,极大的减轻在我国困难地区测图的工作强度,提高的了工作效率,取得了明显的经济效益。
由于我国海岛海岛礁分布零散,选择有效的高程控制点存在很大的困难。在海上无法有效的布设GPS控制点,也无法进行水准高程测量。同时,由于CQG2000似大地水准面模型受地形条件的限制,以及地面观测数据的密度和不均匀分布,致使CQG2000似大地水准面的精度在全国范围内是不一致的。总体分辨率为5′×5′,全国平均精度理论上为±0.36米,不足以满足大比例尺的遥感测图需求。
除少量较大面积的海岛外,海岛(礁)的高程一般不大,在我国,高程超过10m的海岛(礁)数量占全部海岛(礁)的比例不到15%。由于海岛被海水包围,其高程有海平面作绝对参考,大的高程误差甚至会导致海平面与海岛的拓扑关系出现错误。因此,在海岛测图中,海岛高程的准确性和精度至关重要。
在陆地航空航天遥感测绘中,当地面像控点稀少时,可以采用基于POS和影像纹理匹配的稀少控制摄影测量方法进行立体测图。但值得注意的是,由于目前摄影测量工作站的立体测图软件普遍使用正常高系统作为像控点高程,当地面像控点稀少时,测图的高程会产生形变,其形变程度由像控点分布密度和大地水准面起伏两个因素共同决定。当像控点间距大于200km时,大地水准面起伏在我国陆地部分地区可达到20m,如果直接采用稀少控制的航空航天遥感测图方法,虽然不会引起地形图要素的拓扑关系改变,但可导致绝对高程的最大误差达到10m。
与陆地相比,海岛(礁)分布稀疏,无法按常规的密度和点位布测像控点,海域绝大部分被海水覆盖,影像地面特征稀少、纹理少使得航空航天影像匹配的难度大,导致模型连接困难甚至无法实现。采用稀少控制的航空摄影测量方法很可能由于大地水准面起伏大导致海岛地形图高程不准确。
我们知道,与陆地不同,海岛的高程注记不能因为比例尺小而取舍,海岛岸线高程精度应在0.3m左右比较合适。目前陆地地形图中有关高程精度指标普遍较低,用于指导海岛地形测图甚至可能会导致海岛掉到海水面以下的矛盾现象。
发明内容
(一)要解决的技术问题
本发明要解决的技术问题是:如何提高海岛遥感测图所获得的高程精度。
(二)技术方案
为解决上述技术问题,本发明提供了一种海岛遥感测图高程精度控制与修正方法,包括以下步骤:
S1:提取海岛遥感影像中的岸线;
S2:通过海岛验潮站的潮位观测数据获得第一平均海面高、以及卫星测高数据获得第二平均海面高,将所述第二平均海面高进行基准转换,以保证与第一平均海面高的基准一致,并将获得的第一平均海面高和基准转换后的第二平均海面高作为标准平均海面高;
S3:通过地面重力观测数据获得重力数据,通过卫星测高数据获得重力异常数据,将重力异常数据进行基准转换,以保证与所述重力数据的基准一致,并根据所述重力数据、以及基准转换后的重力异常数据来计算大地水准面;
S4:对所述海岛验潮站的潮位观测数据进行调和分析,以获得分潮调和常数;
S5:根据所述分潮调和常数计算获得平均大潮高潮面与深度基准面之间的高度差、以及所述标准平均海面高与深度基准面之间的高度差,并根据所述平均大潮高潮面与深度基准面之间的高度差、所述标准平均海面高与深度基准面之间的高度差、以及所述标准平均海面高与大地水准面之间的高度差来计算所述平均大潮高潮面的高程;
S6:根据所述平均大潮高潮面的高程来确定所述平均大潮高潮面上的岸线,通过所述海岛遥感影像中的岸线、以及所述平均大潮高潮面上的岸线来计算并去除所述海岛遥感影像中的岸线的系统性误差;
S7:将去除系统性误差后的所述海岛遥感影像中的岸线与所述平均大潮高潮面上的岸线进行拟合,以实现对所述海岛遥感影像中的岸线的修正。
优选地,所述海岛遥感影像包括:RapidEye-one、Alos、Spot5、IKONOS、QuickBird、WorldView-one、WorldView-Two、GeoEye-One、以及低空无人机所获得的海岛遥感数据中的至少一种。
优选地,步骤S2具体包括以下步骤:
S21:对海岛验潮站的潮位观测数据计算平均值,以获得第一平均海面高;
S22:根据卫星测高数据中的观测海面高度计算平均值,以获得第二平均海面高;
S23:对所述第二平均海面高进行基准转换,以保证与第一平均海面高的基准一致,在所述第二平均海面高进行基准转换时,采用下列公式,
其中,a为第一平均海面高对应的第一参考椭球的长半轴的长度,f为所述第一参考椭球的扁率,为大地纬度,e为所述第一参考椭球的第一偏心率,da=a0-a、df=f0-f,a0为第二平均海面高对应的第二参考椭球的长半轴的长度,f0为第二参考椭球的扁率,da为长半轴改正值,df为扁率改正值,dh为参考椭球转换引起的高程变化;
S24:将获得的第一平均海面高和基准转换后的第二平均海面高作为标准平均海面高。
优选地,步骤S3具体包括以下步骤:
S31:通过地面重力观测数据获得重力数据;
S32:通过卫星测高数据获得重力异常数据,将重力异常数据进行基准转换,以保证与所述重力数据的基准一致,在所述重力异常数据进行基准转换时,采用下列公式,
其中,a为所述重力数据对应的第一参考椭球的长半轴的长度,f为所述第一参考椭球的扁率,为大地纬度,e为所述第一参考椭球的第一偏心率,da=a0-a、df=f0-f,a0为所述重力异常数据对应的第二参考椭球的长半轴的长度,f0为第二参考椭球的扁率,da为长半轴改正值,df为扁率改正值,dg为参考椭球转换引起的重力异常数据变化;
S33:根据所述重力数据、以及基准转换后的重力异常数据来计算大地水准面,大地水准面的计算公式为,
其中,Δg为基准转换后的重力异常数据,σ为单位球面的面积元,S(ψ)为以球面角距ψ为变量的Stokes核函数,R为地球平均半径,γ为所述重力数据,N为大地水准面。
优选地,步骤S4中,调和分析的公式为,
其中,Hi、gi依次为第i个分潮的振幅和迟角,即分潮调和常数,σi为第i个分潮角速率,V0为参考时刻的平衡潮相角,fi、ui为第i个分潮的交点因子和交点订正角,hi为所述海岛验潮站的第i个潮位观测数据,n为分潮的个数。
优选地,步骤S5中,通过下式来计算所述平均大潮高潮面的高程,
L=L3-(L2-L1)
其中,L为所述平均大潮高潮面的高程,L3为所述平均大潮高潮面与深度基准面之间的高度差,L2为所述标准平均海面高与深度基准面之间的高度差,L1为所述标准平均海面高与大地水准面之间的高度差。
优选地,步骤S6中,通过下式计算所述系统性误差,
其中,ΔL为系统性误差,Lc,i为平均大潮高潮面上的岸线第i点高程值,Lo,i为海岛遥感影像中的岸线第i点高程值,Lc,i和Lo,i的平面位置相同,n为岸线上取样点数。
优选地,步骤S7中,设所述海岛遥感影像中的岸线与所述平均大潮高潮面上的岸线之间任一点的偏差为Z(x,y)且Z(x,y)=a1x+a2y+a3xy+a4x2+a5y2,其中,a1,a2,a3,a4,a5为常数,x,y分别为横轴坐标和纵轴坐标,对所述偏差按照最小二乘法,获得a1,a2,a3,a4,a的值,并根据获得了a1,a2,a3,a4,a5的值的偏差公式实现对所述海岛遥感影像中的岸线的修正。
本发明还公开了一种海岛遥感测图高程精度控制与修正系统,包括:
岸线提取模块,用于提取海岛遥感影像中的岸线;
平均海面计算模块,用于通过海岛验潮站的潮位观测数据获得第一平均海面高、以及卫星测高数据获得第二平均海面高,将所述第二平均海面高进行基准转换,以保证与第一平均海面高的基准一致,并将获得的第一平均海面高和基准转换后的第二平均海面高作为标准平均海面高;
大地水准面计算模块,用于通过地面重力观测数据获得重力数据,通过卫星测高数据获得重力异常数据,将重力异常数据进行基准转换,以保证与所述重力数据的基准一致,并根据所述重力数据、以及基准转换后的重力异常数据来计算大地水准面;
调和分析模块,用于对所述海岛验潮站的潮位观测数据进行调和分析,以获得分潮调和常数;
高潮面计算模块,用于根据所述分潮调和常数计算获得平均大潮高潮面与深度基准面之间的高度差、以及所述标准平均海面高与深度基准面之间的高度差,并根据所述平均大潮高潮面与深度基准面之间的高度差、所述标准平均海面高与深度基准面之间的高度差、以及所述标准平均海面高与大地水准面之间的高度差来计算所述平均大潮高潮面的高程;
系统性误差计算模块,用于根据所述平均大潮高潮面的高程来确定所述平均大潮高潮面上的岸线,通过所述海岛遥感影像中的岸线、以及所述平均大潮高潮面上的岸线来计算并去除所述海岛遥感影像中的岸线的系统性误差;
拟合修正模块,用于将去除系统性误差后的所述海岛遥感影像中的岸线与所述平均大潮高潮面上的岸线进行拟合,以实现对所述海岛遥感影像中的岸线的修正。
(三)有益效果
本发明通过对海岛遥感影像中的岸线与所述平均大潮高潮面上的岸线进行拟合,来实现对所述海岛遥感影像中的岸线的修正,提高了海岛遥感测图所获得的高程精度。
附图说明
图1是按照本发明一种实施方式的海岛遥感测图高程精度控制与修正方法的流程图;
图2是平均大潮高潮面与其他参数的示意图。
具体实施方式
下面结合附图和实施例,对本发明的具体实施方式作进一步详细描述。以下实施例用于说明本发明,但不用来限制本发明的范围。
图1是按照本发明一种实施方式的海岛遥感测图高程精度控制与修正方法的流程图,参照图1,本实施方式的方法包括以下步骤:
S1:提取海岛遥感影像中的岸线;优选地,所述海岛遥感影像包括:RapidEye-one、Alos、Spot5、IKONOS、QuickBird、WorldView-one、WorldView-Two、GeoEye-One、以及低空无人机所获得的海岛遥感数据中的至少一种。
S2:通过海岛验潮站的潮位观测数据获得第一平均海面高、以及卫星测高数据获得第二平均海面高,将所述第二平均海面高进行基准转换,以保证与第一平均海面高的基准一致,并将获得的第一平均海面高和基准转换后的第二平均海面高作为标准平均海面高;
S3:通过地面重力观测数据获得重力数据,通过卫星测高数据获得重力异常数据,将重力异常数据进行基准转换,以保证与所述重力数据的基准一致,并根据所述重力数据、以及基准转换后的重力异常数据来计算大地水准面;
S4:对所述海岛验潮站的潮位观测数据进行调和分析,以获得分潮调和常数;
S5:根据所述分潮调和常数计算获得平均大潮高潮面与深度基准面之间的高度差、以及所述标准平均海面高与深度基准面之间的高度差,并根据所述平均大潮高潮面与深度基准面之间的高度差、所述标准平均海面高与深度基准面之间的高度差、以及所述标准平均海面高与大地水准面之间的高度差来计算所述平均大潮高潮面的高程;
S6:根据所述平均大潮高潮面的高程来确定所述平均大潮高潮面上的岸线,通过所述海岛遥感影像中的岸线、以及所述平均大潮高潮面上的岸线来计算并去除所述海岛遥感影像中的岸线的系统性误差;
S7:将去除系统性误差后的所述海岛遥感影像中的岸线与所述平均大潮高潮面上的岸线进行拟合,以实现对所述海岛遥感影像中的岸线的修正。
优选地,步骤S2具体包括以下步骤:
S21:对海岛验潮站的潮位观测数据计算平均值,以获得第一平均海面高;所述第一平均海面高通过如下公式进行计算,
式中,hi为海岛验潮站的第i个潮位观测数据,n为潮位观测数据个数,his为第一平均海面高。
S22:根据卫星测高数据中的观测海面高度计算平均值,以获得第二平均海面高;所述第二平均海面高通过如下公式进行计算,
式中,hk为卫星测高数据中的第i个观测海面高度,m为潮位观测数据个数,hsa为第二平均海面高。
S23:对所述第二平均海面高进行基准转换,以保证与第一平均海面高的基准一致,在所述第二平均海面高进行基准转换时,采用下列公式,
其中,a为第一平均海面高对应的第一参考椭球的长半轴的长度,f为所述第一参考椭球的扁率,为大地纬度,e为所述第一参考椭球的第一偏心率,da=a0-a、df=f0-f,a0为第二平均海面高对应的第二参考椭球的长半轴的长度,f0为第二参考椭球的扁率,da为长半轴改正值,df为扁率改正值,dh为参考椭球转换引起的高程变化;
S24:将获得的第一平均海面高和基准转换后的第二平均海面高作为标准平均海面高。
优选地,步骤S3具体包括以下步骤:
S31:通过地面重力观测数据获得重力数据;
S32:通过卫星测高数据获得重力异常数据,将重力异常数据进行基准转换,以保证与所述重力数据的基准一致,在所述重力异常数据进行基准转换时,采用下列公式,
其中,a为所述重力数据对应的第一参考椭球的长半轴的长度,f为所述第一参考椭球的扁率,为大地纬度,e为所述第一参考椭球的第一偏心率,da=a0-a、df=f0-f,a0为所述重力异常数据对应的第二参考椭球的长半轴的长度,f0为第二参考椭球的扁率,da为长半轴改正值,df为扁率改正值,dg为参考椭球转换引起的重力异常数据变化;
S33:根据所述重力数据、以及基准转换后的重力异常数据来计算大地水准面,大地水准面的计算公式为,
其中,Δg为基准转换后的重力异常数据,σ为单位球面的面积元,S(ψ)为以球面角距ψ为变量的Stokes核函数,R为地球平均半径,γ为所述重力数据,N为大地水准面。
优选地,步骤S4中,调和分析的公式为,
其中,Hi、gi依次为第i个分潮的振幅和迟角,即分潮调和常数,σi为第i个分潮角速率,V0为参考时刻的平衡潮相角,fi、ui为第i个分潮的交点因子和交点订正角,hi为所述海岛验潮站的第i个潮位观测数据,n为分潮的个数,本实施方式中选用M2、S2、N2、K2、01、K1、P1、Q1这8个分潮。
优选地,参照图2,步骤S5中,通过下式来计算所述平均大潮高潮面的高程,
L=L3-(L2-L1)
其中,L为所述平均大潮高潮面的高程,L3为所述平均大潮高潮面与深度基准面之间的高度差,L2为所述标准平均海面高与深度基准面之间的高度差,L1为所述标准平均海面高与大地水准面之间的高度差。
优选地,步骤S6中,通过下式计算所述系统性误差,
其中,ΔL为系统性误差,Lc,i为平均大潮高潮面上的岸线第i点高程值,Lo,i为海岛遥感影像中的岸线第i点高程值,Lc,i和Lo,i的平面位置相同,n为岸线上取样点数。
优选地,步骤S7中,设所述海岛遥感影像中的岸线与所述平均大潮高潮面上的岸线之间任一点的偏差为Z(x,y)且Z(x,y)=a1x+a2y+a3xy+a4x2+a5y2,其中,a1,a2,a3,a4,a5为常数,x,y分别为横轴坐标和纵轴坐标,对所述偏差按照最小二乘法,获得a1,a2,a3,a4,a的值,并根据获得了a1,a2,a3,a4,a5的值的偏差公式实现对所述海岛遥感影像中的岸线的修正。
除少量较大面积的海岛外,海岛(礁)的高程一般不大,在我国,高程超过10m的海岛(礁)数量占全部海岛(礁)的比例不到15%。由于海岛被海水包围,其高程有海平面作绝对参考,大的高程误差甚至会导致海平面与海岛的拓扑关系出现错误。在海岛测图中,海岛高程的准确性和精度至关重要。
在陆地航空航天遥感测图中,当地面像控点稀少时,可以采用基于POS和影像纹理匹配的稀少控制摄影测量方法进行立体测图。像控点分布密度和大地水准面起伏两个因素共同决定在地面像控点稀少时,测图的高程会产生形变程度。
当像控点间距大于200km时,大地水准面起伏在我国陆地部分地区可达到20m,如果直接采用稀少控制的航空航天遥感测图方法,虽然不会引起地形图要素的拓扑关系改变,但可导致绝对高程的最大误差达到10m。
与陆地相比,海岛(礁)分布稀疏,无法按常规的密度和点位布测像控点,海域绝大部分被海水覆盖,影像地面特征稀少、纹理少使得航空航天影像匹配的难度大,导致模型连接困难甚至无法实现。采用稀少控制的航空摄影测量方法很可能由于大地水准面起伏大导致海岛地形图高程不准确。目前陆地地形图中有关高程精度指标普遍较低,用于指导海岛地形测图甚至可能会导致海岛掉到海水面以下的矛盾现象。有效解决这些问题的办法是对海岛遥感测图成果进行高程修正。
本发明的方法计算得到的高程精度仅受计算的误差影响,以目前的技术水平看,高程精度可达到0.2m,能够满足1∶500以小比例尺测图对绝对高程的精度要求。
本发明还公开了一种海岛遥感测图高程精度控制与修正系统,包括:
岸线提取模块,用于提取海岛遥感影像中的岸线;
平均海面计算模块,用于通过海岛验潮站的潮位观测数据获得第一平均海面高、以及卫星测高数据获得第二平均海面高,将所述第二平均海面高进行基准转换,以保证与第一平均海面高的基准一致,并将获得的第一平均海面高和基准转换后的第二平均海面高作为标准平均海面高;
大地水准面计算模块,用于通过地面重力观测数据获得重力数据,通过卫星测高数据获得重力异常数据,将重力异常数据进行基准转换,以保证与所述重力数据的基准一致,并根据所述重力数据、以及基准转换后的重力异常数据来计算大地水准面;
调和分析模块,用于对所述海岛验潮站的潮位观测数据进行调和分析,以获得分潮调和常数;
高潮面计算模块,用于根据所述分潮调和常数计算获得平均大潮高潮面与深度基准面之间的高度差、以及第二平均海面高与深度基准面之间的高度差,并根据所述平均大潮高潮面与深度基准面之间的高度差、所述标准平均海面高与深度基准面之间的高度差、以及所述标准平均海面高与大地水准面之间的高度差来计算所述平均大潮高潮面的高程;
系统性误差计算模块,用于根据所述平均大潮高潮面的高程来确定所述平均大潮高潮面上的岸线,通过所述海岛遥感影像中的岸线、以及所述平均大潮高潮面上的岸线来计算并去除所述海岛遥感影像中的岸线的系统性误差;
拟合修正模块,用于将去除系统性误差后的所述海岛遥感影像中的岸线与所述平均大潮高潮面上的岸线进行拟合,以实现对所述海岛遥感影像中的岸线的修正。
以上实施方式仅用于说明本发明,而并非对本发明的限制,有关技术领域的普通技术人员,在不脱离本发明的精神和范围的情况下,还可以做出各种变化和变型,因此所有等同的技术方案也属于本发明的范畴,本发明的专利保护范围应由权利要求限定。
Claims (9)
1.一种海岛遥感测图高程精度控制与修正方法,其特征在于,包括以下步骤:
S1:提取海岛遥感影像中的岸线;
S2:通过海岛验潮站的潮位观测数据获得第一平均海面高、以及卫星测高数据获得第二平均海面高,将所述第二平均海面高进行基准转换,以保证与第一平均海面高的基准一致,并将获得的第一平均海面高和基准转换后的第二平均海面高作为标准平均海面高;
S3:通过地面重力观测数据获得重力数据,通过卫星测高数据获得重力异常数据,将重力异常数据进行基准转换,以保证与所述重力数据的基准一致,并根据所述重力数据、以及基准转换后的重力异常数据来计算大地水准面;
S4:对所述海岛验潮站的潮位观测数据进行调和分析,以获得分潮调和常数;
S5:根据所述分潮调和常数计算获得平均大潮高潮面与深度基准面之间的高度差、以及所述标准平均海面高与深度基准面之间的高度差,并根据所述平均大潮高潮面与深度基准面之间的高度差、所述标准平均海面高与深度基准面之间的高度差、以及所述标准平均海面高与大地水准面之间的高度差来计算所述平均大潮高潮面的高程;
S6:根据所述平均大潮高潮面的高程来确定所述平均大潮高潮面上的岸线,通过所述海岛遥感影像中的岸线、以及所述平均大潮高潮面上的岸线来计算并去除所述海岛遥感影像中的岸线的系统性误差;
S7:将去除系统性误差后的所述海岛遥感影像中的岸线与所述平均大潮高潮面上的岸线进行拟合,以实现对所述海岛遥感影像中的岸线的修正。
2.如权利要求1所述的方法,其特征在于,所述海岛遥感影像包括:RapidEye-one、Alos、Spot5、IKONOS、QuickBird、WorldView-one、WorldView-Two、GeoEye-One、以及低空无人机所获得的海岛遥感数据中的至少一种。
3.如权利要求1所述的方法,其特征在于,步骤S2具体包括以下步骤:
S21:对海岛验潮站的潮位观测数据计算平均值,以获得第一平均海面高;
S22:根据卫星测高数据中的观测海面高度计算平均值,以获得第二平均海面高;
S23:对所述第二平均海面高进行基准转换,以保证与第一平均海面高的基准一致,在所述第二平均海面高进行基准转换时,采用下列公式,
其中,a为第一平均海面高对应的第一参考椭球的长半轴的长度,f为所述第一参考椭球的扁率,为大地纬度,e为所述第一参考椭球的第一偏心率,da=a0-a、df=f0-f,a0为第二平均海面高对应的第二参考椭球的长半轴的长度,f0为第二参考椭球的扁率,da为长半轴改正值,df为扁率改正值,dh为参考椭球转换引起的高程变化;
S24:将获得的第一平均海面高和基准转换后的第二平均海面高作为标准平均海面高。
4.如权利要求1所述的方法,其特征在于,步骤S3具体包括以下步骤:
S31:通过地面重力观测数据获得重力数据;
S32:通过卫星测高数据获得重力异常数据,将重力异常数据进行基准转换,以保证与所述重力数据的基准一致,在所述重力异常数据进行基准转换时,采用下列公式,
其中,a为所述重力数据对应的第一参考椭球的长半轴的长度,f为所述第一参考椭球的扁率,为大地纬度,e为所述第一参考椭球的第一偏心率,da=a0-a、df=f0-f,a0为所述重力异常数据对应的第二参考椭球的长半轴的长度,f0为第二参考椭球的扁率,da为长半轴改正值,df为扁率改正值,dg为参考椭球转换引起的重力异常数据变化;
S33:根据所述重力数据、以及基准转换后的重力异常数据来计算大地水准面,大地水准面的计算公式为,
其中,Δg为基准转换后的重力异常数据,σ为单位球面的面积元,S(ψ)为以球面角距ψ为变量的Stokes核函数,R为地球平均半径,γ为所述重力数据,N为大地水准面。
5.如权利要求1所述的方法,其特征在于,步骤S4中,调和分析的公式为,
其中,Hi、gi依次为第i个分潮的振幅和迟角,即分潮调和常数,σi为第i个分潮角速率,V0为参考时刻的平衡潮相角,fi、ui为第i个分潮的交点因子和交点订正角,hi为所述海岛验潮站的第i个潮位观测数据,n为分潮的个数。
6.如权利要求1所述的方法,其特征在于,步骤S5中,通过下式来计算所述平均大潮高潮面的高程,
L=L3-(L2-L1)
其中,L为所述平均大潮高潮面的高程,L3为所述平均大潮高潮面与深度基准面之间的高度差,L2为所述标准平均海面高与深度基准面之间的高度差,L1为所述标准平均海面高与大地水准面之间的高度差。
7.如权利要求1所述的方法,其特征在于,步骤S6中,通过下式计算所述系统性误差,
其中,ΔL为系统性误差,Lc,i为平均大潮高潮面上的岸线第i点高程值,Lo,i为海岛遥感影像中的岸线第i点高程值,Lc,i和Lo,i的平面位置相同,n为岸线上取样点数。
8.如权利要求1所述的方法,其特征在于,步骤S7中,设所述海岛遥感影像中的岸线与所述平均大潮高潮面上的岸线之间任一点的偏差为Z(x,y)且Z(x,y)=a1x+a2y+a3xy+a4x2+a5y2,其中,a1,a2,a3,a4,a5为常数,x,y分别为横轴坐标和纵轴坐标,对所述偏差按照最小二乘法,获得a1,a2,a3,a4,a5的值,并根据获得了a1,a2,a3,a4,a5的值的偏差公式实现对所述海岛遥感影像中的岸线的修正。
9.一种海岛遥感测图高程精度控制与修正系统,其特征在于,包括:
岸线提取模块,用于提取海岛遥感影像中的岸线;
平均海面计算模块,用于通过海岛验潮站的潮位观测数据获得第一平均海面高、以及卫星测高数据获得第二平均海面高,将所述第二平均海面高进行基准转换,以保证与第一平均海面高的基准一致,并将获得的第一平均海面高和基准转换后的第二平均海面高作为标准平均海面高;
大地水准面计算模块,用于通过地面重力观测数据获得重力数据,通过卫星测高数据获得重力异常数据,将重力异常数据进行基准转换,以保证与所述重力数据的基准一致,并根据所述重力数据、以及基准转换后的重力异常数据来计算大地水准面;
调和分析模块,用于对所述海岛验潮站的潮位观测数据进行调和分析,以获得分潮调和常数;
高潮面计算模块,用于根据所述分潮调和常数计算获得平均大潮高潮面与深度基准面之间的高度差、以及所述标准平均海面高与深度基准面之间的高度差,并根据所述平均大潮高潮面与深度基准面之间的高度差、所述标准平均海面高与深度基准面之间的高度差、以及所述标准平均海面高与大地水准面之间的高度差来计算所述平均大潮高潮面的高程;
系统性误差计算模块,用于根据所述平均大潮高潮面的高程来确定所述平均大潮高潮面上的岸线,通过所述海岛遥感影像中的岸线、以及所述平均大潮高潮面上的岸线来计算并去除所述海岛遥感影像中的岸线的系统性误差;
拟合修正模块,用于将去除系统性误差后的所述海岛遥感影像中的岸线与所述平均大潮高潮面上的岸线进行拟合,以实现对所述海岛遥感影像中的岸线的修正。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210025062.9A CN103245324B (zh) | 2012-02-06 | 2012-02-06 | 海岛遥感测图高程精度控制与修正方法及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201210025062.9A CN103245324B (zh) | 2012-02-06 | 2012-02-06 | 海岛遥感测图高程精度控制与修正方法及系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103245324A true CN103245324A (zh) | 2013-08-14 |
CN103245324B CN103245324B (zh) | 2015-04-08 |
Family
ID=48925005
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201210025062.9A Expired - Fee Related CN103245324B (zh) | 2012-02-06 | 2012-02-06 | 海岛遥感测图高程精度控制与修正方法及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103245324B (zh) |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104155297A (zh) * | 2014-08-11 | 2014-11-19 | 江苏恒创软件有限公司 | 一种基于无人机的工业水体污染源检测方法 |
CN104567802A (zh) * | 2015-01-06 | 2015-04-29 | 山东科技大学 | 集成船载重力和gnss的测线式陆海高程传递方法 |
CN109341650A (zh) * | 2018-09-06 | 2019-02-15 | 中石油昆仑燃气有限公司 | 一种基于最小检测代价的无人机高程误差双门限修正方法 |
CN111611540A (zh) * | 2020-04-08 | 2020-09-01 | 中国电建集团西北勘测设计研究院有限公司 | 一种基于千寻位置的像控点高程精密云计算的转换方法 |
CN114266796A (zh) * | 2021-11-05 | 2022-04-01 | 广东省国土资源测绘院 | 基于激光点云数据与平均大潮高潮面的自然岸线自动获取方法、设备、介质、产品 |
CN115222909A (zh) * | 2022-09-19 | 2022-10-21 | 深圳市湾测技术有限公司 | 一种三维图像检测模型构建方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5602731A (en) * | 1986-12-23 | 1997-02-11 | Petroscan Ab | Method for mapping sea level undulations with applications to mineral and hydrocarbon prospecting |
US20070050142A1 (en) * | 2005-08-23 | 2007-03-01 | Pamela Posey | System and method for estimating ocean height and current on a personal computer with location adjustment |
CN101957193A (zh) * | 2010-06-12 | 2011-01-26 | 中国测绘科学研究院 | 一种海岛礁高程传递的优化方法 |
CN102175213A (zh) * | 2011-01-30 | 2011-09-07 | 中国测绘科学研究院 | Gps水准高程传递方法 |
-
2012
- 2012-02-06 CN CN201210025062.9A patent/CN103245324B/zh not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5602731A (en) * | 1986-12-23 | 1997-02-11 | Petroscan Ab | Method for mapping sea level undulations with applications to mineral and hydrocarbon prospecting |
US20070050142A1 (en) * | 2005-08-23 | 2007-03-01 | Pamela Posey | System and method for estimating ocean height and current on a personal computer with location adjustment |
CN101957193A (zh) * | 2010-06-12 | 2011-01-26 | 中国测绘科学研究院 | 一种海岛礁高程传递的优化方法 |
CN102175213A (zh) * | 2011-01-30 | 2011-09-07 | 中国测绘科学研究院 | Gps水准高程传递方法 |
Non-Patent Citations (2)
Title |
---|
柯宝贵等: "原理大陆海岛的高程传递", 《测绘通报》 * |
胡建国等: "近海多卫星测高数据联合处理的方法及应用", 《测绘通报》 * |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104155297A (zh) * | 2014-08-11 | 2014-11-19 | 江苏恒创软件有限公司 | 一种基于无人机的工业水体污染源检测方法 |
CN104567802A (zh) * | 2015-01-06 | 2015-04-29 | 山东科技大学 | 集成船载重力和gnss的测线式陆海高程传递方法 |
CN109341650A (zh) * | 2018-09-06 | 2019-02-15 | 中石油昆仑燃气有限公司 | 一种基于最小检测代价的无人机高程误差双门限修正方法 |
CN109341650B (zh) * | 2018-09-06 | 2021-05-25 | 中国石油天然气股份有限公司 | 一种基于最小检测代价的无人机高程误差双门限修正方法 |
CN111611540A (zh) * | 2020-04-08 | 2020-09-01 | 中国电建集团西北勘测设计研究院有限公司 | 一种基于千寻位置的像控点高程精密云计算的转换方法 |
CN111611540B (zh) * | 2020-04-08 | 2024-03-15 | 中国电建集团西北勘测设计研究院有限公司 | 一种基于千寻位置的像控点高程精密云计算的转换方法 |
CN114266796A (zh) * | 2021-11-05 | 2022-04-01 | 广东省国土资源测绘院 | 基于激光点云数据与平均大潮高潮面的自然岸线自动获取方法、设备、介质、产品 |
CN114266796B (zh) * | 2021-11-05 | 2022-09-13 | 广东省国土资源测绘院 | 基于激光点云数据与平均大潮高潮面的自然岸线自动获取方法、设备、介质 |
CN115222909A (zh) * | 2022-09-19 | 2022-10-21 | 深圳市湾测技术有限公司 | 一种三维图像检测模型构建方法 |
CN115222909B (zh) * | 2022-09-19 | 2022-11-29 | 深圳市湾测技术有限公司 | 一种三维图像检测模型构建方法 |
Also Published As
Publication number | Publication date |
---|---|
CN103245324B (zh) | 2015-04-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103245324B (zh) | 海岛遥感测图高程精度控制与修正方法及系统 | |
CN105938193B (zh) | 一种无需地面辅助的升降轨InSAR监测沉降区绝对地表形变方法 | |
Heck | An evaluation of some systematic error sources affecting terrestrial gravity anomalies | |
CN102607506B (zh) | 高填方机场边坡单台全站仪的自由设站变形监测方法 | |
CN103207419B (zh) | 隧道岩层产状三维测量方法 | |
CN106767671B (zh) | 基于三维电子罗盘的地质结构面产状计算方法 | |
CN102175213B (zh) | Gps水准高程传递方法 | |
CN105158760A (zh) | 一种利用InSAR反演地下流体体积变化和三维地表形变的方法 | |
CN102589528B (zh) | 一种多时相影像的海岛岸线量测方法 | |
CN104567802B (zh) | 集成船载重力和gnss的测线式陆海高程传递方法 | |
CN102589517A (zh) | 基于egm2008的区域似大地水准面精化方法 | |
CN108427741A (zh) | 一种基于大量高精度控制点的dem相对误差评价方法 | |
CN102236108A (zh) | 一种磁性地表三维地形改正方法 | |
CN103400416A (zh) | 一种基于概率多层地形的城市环境机器人导航方法 | |
CN109631863A (zh) | 一种空地结合的潮间带一体化测绘方法 | |
CN106568412A (zh) | 一种不规则形状地块的面积测量方法 | |
CN114493929A (zh) | 一种基于空间分析技术的隧道富水区划分方法 | |
CN101957193B (zh) | 一种海岛礁高程传递的优化方法 | |
CN105571598A (zh) | 一种卫星激光高度计足印相机姿态的测定方法 | |
CN102661741B (zh) | 一种小型冲沟沟头形态及侵蚀量的监测方法 | |
Abili | Comparison of vertical accuracy of open-source global digital elevation models: a case study of Adama City, Ethiopia | |
CN113819882B (zh) | 一种跨海高程点间重力位差计算方法 | |
CN113189660B (zh) | 一种阵列式陆地时变重力和梯度场的观测方法和系统 | |
Mirijovsky et al. | UAV photogrammetry in fluvial geomorphology | |
Dai et al. | Application of UAV photogrammetry on ecological restoration of abandoned open-pit mines, Northern Anhui province, China |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C53 | Correction of patent of invention or patent application | ||
CB03 | Change of inventor or designer information |
Inventor after: Ke Baogui Inventor after: Zhang Chuanyin Inventor after: Zhang Liming Inventor after: Li Dehai Inventor before: Ke Baogui Inventor before: Zhang Chuanyin Inventor before: Zhang Liming Inventor before: Li Dehai |
|
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20150408 Termination date: 20180206 |
|
CF01 | Termination of patent right due to non-payment of annual fee |