CN102313537B - 基于单片InSAR正射影像的调绘方法 - Google Patents
基于单片InSAR正射影像的调绘方法 Download PDFInfo
- Publication number
- CN102313537B CN102313537B CN 201110209941 CN201110209941A CN102313537B CN 102313537 B CN102313537 B CN 102313537B CN 201110209941 CN201110209941 CN 201110209941 CN 201110209941 A CN201110209941 A CN 201110209941A CN 102313537 B CN102313537 B CN 102313537B
- Authority
- CN
- China
- Prior art keywords
- atural object
- insar
- read
- monolithic
- atural
- 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.)
- Expired - Fee Related
Links
Images
Landscapes
- Image Analysis (AREA)
- Length Measuring Devices By Optical Means (AREA)
Abstract
本发明公开了一种基于单片InSAR正射影像的调绘方法,包括步骤:一、单片InSAR正射影像获取;二、基于单片InSAR影像的调绘,其调绘过程如下:201、建立InSAR影像地物特征库;202、室内判读:根据InSAR影像地物特征库对地物进行定性和定位;203、野外核查;204、阴影区和叠掩区的调绘与补测;204、新增地物与隐蔽地物的调绘与补测;205、室内数据处理;三、接边处理。本发明设计合理、方法步骤简单、实现方便且使用效果好、调绘精度高,能解决现有全数字航空摄影测量系统及调绘方法无法满足单片InSAR影像的调绘需要、单片InSAR影像调绘效果差、不能充分发挥InSAR技术优势等实际问题。
Description
技术领域
本发明属于机载合成孔径雷达干涉测量技术领域,尤其是涉及一种基于单片InSAR正射影像的调绘方法。
背景技术
InSAR(Interferometric Synthetic Aperture Radar;简称:干涉测量合成孔径雷达或合成孔径雷达干涉测量)技术源于美国,在欧美发达国家得到不断完善与成熟,其应用领域也得到不断推广。对于InSAR技术发达的国家,如美国和德国的空间信息化产业技术公司,已将实用化的机载高分辨率InSAR技术作为一种新的、先进的技术手段,用于地形测绘、森林测量、资源调查和环境制图、地质环境和灾害监测等方面,并且随着这一技术的快速发展,新的应用领域还在逐步拓宽。
合成孔径雷达干涉测量(InSAR)技术,经过近二十年的研究,其理论日臻成熟,实际使用过程中,雷达干涉测量具有以下优点:第一、不依赖于太阳光,而是利用自身发射的电磁波进行测量,因此可以全天时工作;第二、除了能穿云破雾之外,还不受天气因素的影响,因此可以全天候工作;第三、雷达干涉测量可以直接获取地形的高程信息。因而,现如今合成孔径雷达干涉测量(InSAR)技术的应用领域也得到不断推广。许多欧美国家已将实用化的机载高分辨率InSAR技术作为一种新的、先进的技术手段,用于地形测绘、森林测量、资源调查和环境制图、地质环境和灾害监测等方面。
我国也在近几年开展了机载高分辨率InSAR技术,在地形测绘方面的应用研究。合成孔径雷达干涉测量(InSAR)技术,特别是机载高分辨率InSAR技术作为一种新的、先进的技术手段,已经逐渐应用于地形测绘。在我国它处于起步阶段,主要在机载InSAR系统的研制阶段,地形测绘则刚刚起步,为此迫切需要开展机载高分辨率InSAR技术在地形测绘方面的应用研究。
由于基于单片InSAR影像的测图不同于立体测图,不能直接在影像上量测垂直地物的高度或比高,需要借助三维DEM数据或野外人工量测,这样一来,工作量比较大、作业过程比较繁琐。
如图1所示,实际测量时,由于单片InSAR正射影像特殊的侧视成像机理,对于高出地表的垂直地物,会产生阴影区和叠掩区。其中,阴影区即雷达波束不能到达之处,因此地面上该部分没有回波返回到雷达天线,从而在InSAR图像上形成阴影区,阴影的影像呈现为黑色。阴影区的长度Lg与垂直地物高度h和侧视角θ之间满足公式h=Lgsinθcosθ。叠掩区,即雷达波束到垂直地物顶部A的时间比到垂直地物底部B的时间短,因而垂直地物的顶部影像a先被记录,底部影像b后被记录,因此顶部影像a和底部影像b的位置在单片InSAR影像中不是同一位置(但实际情况是正射投影A和B两点应为同一点),因而在单片InSAR影像中产生了Lg1的投影宽度,而且Lg1与垂直地物的高度h存在如下关系:h=Lg1tanθ。总之,垂直地物的高度h与其产生的阴影宽度Lg或叠掩宽度Lg1之间都存在着一定的关系,这就是基于单片InSAR影像测量垂直地物高度系统的基本原理。
另外,结合图2,由于单片InSAR正射影像特殊的侧视成像机理,对于高出地表的垂直地物,则产生顶底位移现象,即雷达波束到顶部A的时间比到底部B的时间短,因而顶部影像a先被记录,底部影像b后被记录,因此顶部影像a和底部影像b的位置在单片InSAR影像中不是在同一位置上(但实际情况是正射投影A和B两点应为同一点),则在单片InSAR影像中产生了量为Δr的定位误差,因而单片InSAR影像中对垂直地物进行测量时,存在定位误差。其次,由于在单片InSAR影像中,垂直地物顶部的影像往往因未遮挡而多呈现为突出的高亮点,而底部的影像由于遮挡而不易定位。因此,基于单片InSAR影像垂直地物的采集,可以通过采集其顶部位置,经过误差改正,达到底部定位的目的。
在InSAR系统成像过程中,地物目标的位置在方位向上按平台飞行时序记录成像,在距离向上按地物反射信息先后来记录成像,其构象几何学属于斜距投影,这一特点决定了InSAR影像固有的几何特性和特殊的辐射特性,InSAR影像中必然存在叠掩区、阴影区以及其特有的投影差等现象,而且它不同于光学航空摄影所获取的影像,无法采用原有的调绘方法进行地物的判读和定位。因此基于这样的技术背景,需对基于单片InSAR正射影像的调绘方法进行进一步研究。
目前,我国采用的数据采集及处理系统主要是全数字航空摄影测量系统,该系统主要适应于立体像对的数据采集及处理,不适合基于单片InSAR影像的采集及处理,而且对于垂直地物的高度无法从单片InSAR影像中直接获取,因而现有的全数字航空摄影测量系统及调绘方法,无法满足单片InSAR影像调绘的需要。
发明内容
本发明所要解决的技术问题在于针对上述现有技术中的不足,提供一种基于单片InSAR正射影像的调绘方法,其设计合理、方法步骤简单、实现方便且使用效果好、调绘精度高,能有效解决现有全数字航空摄影测量系统及调绘方法无法满足单片InSAR影像的调绘需要、单片InSAR影像调绘效果差、不能充分发挥InSAR技术的优势等问题。
为解决上述技术问题,本发明采用的技术方案是:一种基于单片InSAR正射影像的调绘方法,其特征在于该方法包括以下步骤:
步骤一、单片InSAR正射影像获取:采用机载InSAR系统且按常规的机载合成孔径雷达干涉测量方法,获取测区的多张单片InSAR影像,并将所获取的多张所述单片InSAR影像均存储至与数据处理器相接的存储单元内;多张所述单片InSAR影像均为正射影像;
步骤二、基于单片InSAR正射影像进行调绘:对步骤一中存储至存储单元内的多张所述单片InSAR影像分别进行调绘,且对多张所述单片InSAR影像的调绘方法均相同;对于其中任一张单片InSAR影像而言,其调绘方法包括以下步骤:
201、室内判读,其室内判读过程如下:
2011、影像放大处理及同步显示:通过数据处理器对被判读单片InSAR影像进行放大,并通过与数据处理器相接的显示单元对放大后的被判读单片InSAR影像进行同步显示;
2012、阴影区、叠掩区和可判读区域确定:人为对通过显示单元进行同步显示的被判读单片InSAR影像进行视觉观察,找出被判读单片InSAR影像上存在的阴影区和叠掩区;同时,通过数据处理器对找出的阴影区和叠掩区进行标记;相应地,被判读单片InSAR影像上除被标记的阴影区和叠掩区之外的区域,为可判读区域;
2013、室内判读:利用数据处理器且根据预先在存储单元内所建立的InSAR影像地物特征库,对步骤2012中所述可判读区域内按照成图比例尺的技术规范要求需进行判读的所有地物分别进行判读;所述InSAR影像地物特征库内存储有多种不同类型地物的参考特征参数,且每一种类型地物的参考特征参数均包括地物类型、地物性质、地物在InSAR影像上的形状和大小尺寸以及地物在InSAR影像上的色调和灰度值;实际进行室内判读时,对所述可判读区域内需进行判读的所有地物的判读方法均相同;对于所述可判读区域内需进行判读的任一个地物来说,其室内判读过程包括以下步骤:
20131、地物性质判断:采用处理器自动判读方法或人为判读方法对被判读地物的性质进行室内判读;
当采用处理器自动判读方法对被判读地物的性质进行室内判读时,其室内判读过程如下:
Ⅰ、被判读地物的特征参数获取:所述数据处理器调用尺寸量测模块和灰度值读取模块,在被判读单片InSAR影像上量测出被判读地物的大小尺寸且同步读取出被判读地物的灰度值;同时,根据被判读地物的视觉观察结果,人为通过与数据处理器相接的参数输入单元输入被判读地物在InSAR影像上的形状和在InSAR影像上的色调信息,此时数据处理器获得被判读地物的特征参数并同步存储至存储单元内;被判读地物的特征参数包括被判读地物在InSAR影像上的形状、大小尺寸以及色调信息和灰度值;
Ⅱ、特征匹配并确定被判读地物性质:所述数据处理器调取地物特征参数匹配模块,将被判读地物的特征参数与所述InSAR影像地物特征库进行特征匹配:当匹配出与被判读地物的特征参数一致的参考特征参数时,匹配得出的参考特征参数中的地物类型和地物性质便为被判读地物的类型和性质,此时数据处理器调用标注模块在被判读单片InSAR影像上对被判读地物的类型和性质进行标注;当未匹配出与被判读地物的特征参数一致的参考特征参数时,所述数据处理器调用标注模块将被判读地物的性质标注为未知地物;
当采用人为判读方法对被判读地物的性质进行室内判读时,其室内判读过程如下:
i、InSAR影像地物特征库学习:作业人员通过对所述InSAR影像地物特征库内所存储的多种不同类型地物的参考特征参数进行学习,了解并掌握各种不同类型地物在InSAR影像上的目视判读特征,所述目视判读特征包括地物在InSAR影像上的形状、大小尺寸以及色调信息和灰度值;
ⅱ、人为目视判读:作业人员结合步骤i中了解并掌握的各种不同类型地物在InSAR影像上的目视判读特征,对被判读地物进行目视判读,并相应判读出被判读地物的类型和性质;当目视判读得出被判读地物的类型和性质时,作业人员通过参数输入单元输入目视判读结果,之后数据处理器调用标注模块在被判读单片InSAR影像上对被判读地物的类型和性质进行标注;
当不能目视判读得出被判读地物的类型和性质时,作业人员通过参数输入单元输入目视判读结果且此时目视判读结果为未知地物,之后数据处理器调用标注模块在被判读单片InSAR影像上将被判读地物的性质标注为未知地物;
20132、地物定位:首先根据被判读地物在被判读单片InSAR影像上的阴影方向,判断出被判读地物的顶部位置,被判读地物的顶部位置与其阴影方向相反;其次,所述数据处理器调用位置标定模块且按照常规数字地图标定方法对被判读地物的顶部位置进行标定,标定的同时数据处理器调用平面坐标确定模块且根据预先在被判读单片InSAR影像上所选定参照物的平面坐标,确定被判读地物顶部位置的平面坐标并同步存储至存储单元内;此时,存储单元内所存储被判读地物顶部位置的平面坐标为被判读地物的初步定位结果;
202、野外核查:到实地对步骤20131中判读得出的可判读区域内各地物的地物性质进行核查确认,当核查出步骤20131中判读得出的地物性质错误时,作业人员根据野外核查结果对在被判读单片InSAR影像上所标注的错误地物性质进行修改;同时,对步骤20131中标注为未知地物的地物类型和性质进行实地确定,并根据实地确定结果在被判读单片InSAR影像上对所述未知地物的类型和性质进行标注;
203、阴影区和叠掩区的调绘与补测:采用常规定位测量设备且按照常规的测量和定位方法,对步骤2012中标记的阴影区和叠掩区分别进行调绘与补测,并将调绘与补测结果存储至可移动存储设备内;
204、新增地物与隐蔽地物的调绘与补测:按照常规外业像片调绘方法,对新增地物和隐蔽地物进行调绘与补测,并将调绘与补测结果存储至所述可移动存储设备内;所述新增地物为步骤一中所述机载合成孔径雷达干涉测量结束后且对所述测区的调绘工作结束前,在所述测区内新增设的地物;所述隐蔽地物为步骤一中进行机载合成孔径雷达干涉测量时因与周围相邻地物的InSAR影像混为一体且使得该地物无法辨识的地物,或者步骤一中进行机载合成孔径雷达干涉测量时被周围相邻的高大地物遮挡而使得该地物没有回波返回到机载InSAR系统雷达天线的地物;
205、室内数据处理,其数据处理过程包括以下步骤:
2051、野外调绘补测数据处理:将步骤203中和步骤204中所述可移动存储设备内存储的调绘与补测结果传送至数据处理器,且数据处理器调用数字展绘模块或数字转绘模块,将所述调绘与补测结果展绘或转绘至被判读单片InSAR影像上;
2052、垂直地物高度测定:采用基于单片InSAR影像的室内高度测量方法或基于野外高度测量的室内数据处理方法,对步骤20132中标注的所有垂直地物的高度进行测定,且对所有垂直地物的高度测定方法均相同;所述垂直地物为在测区地面的垂直方向上有高度的地物;
对于任一个垂直地物来说,当采用基于单片InSAR影像的室内高度测量方法对被测量垂直地物高度进行测定时,数据处理器调用垂直地物高度测量模块且根据公式h=Lgsinθcosθ或h=Lg1tanθ,计算得出被测量垂直地物的高度h;式中θ为步骤一中进行机载合成孔径雷达干涉测量时所述机载InSAR系统的侧视角,且侧视角θ通过参数输入单元进行输入;Lg为量测得到的被纠正垂直地物在被判读单片InSAR影像上的阴影宽度,Lg1为量测得到的被纠正垂直地物在被判读单片InSAR影像上的叠掩宽度;之后,所述数据处理器调用标注模块将被测量垂直地物的高度标注为h;
当采用基于野外高度测量的室内数据处理方法对被测量垂直地物高度进行测定时,步骤202中所述的到实地对各地物的地物性质进行核查确认过程中,还需采用常规的定位测量设备对被测量垂直地物的高度进行实地测量,并将被测量垂直地物的高度实地测量结果存储至可移动存储设备内;步骤2052中进行垂直地物高度测定时,根据所述可移动存储设备内所存储的高度实地测量结果,所述数据处理器调用标注模块对被测量垂直地物的高度进行标注;
2053、垂直地物定位误差纠正:采用室内定位误差纠正方法或基于野外定位测量的室内数据纠正方法,对步骤20132中标注的所有垂直地物的定位误差进行纠正,且对所有垂直地物的定位误差纠正方法均相同;
对于任一个垂直地物来说,当采用室内定位误差纠正方法对被纠正垂直地物的定位误差进行纠正时,数据处理器调用误差纠正模块且根据公式Δr=h/tanη,计算得出被纠正垂直地物因高度引起的顶部平面偏移量Δr,式中η为步骤一中进行机载合成孔径雷达干涉测量时所述机载InSAR系统的入射角;之后,数据处理器根据计算得出的顶部平面偏移量Δr,对步骤20132中存储单元内所存储被纠正垂直地物的初步定位结果进行纠正,获得纠正后的定位结果;然后,所述数据处理器调用标注模块将被纠正垂直地物的底部位置平面坐标标注为纠正后的定位结果;
当采用基于野外定位测量的室内数据纠正方法对被纠正垂直地物的定位误差进行纠正时,步骤202中所述的到实地对各地物的地物性质进行核查确认过程中,还需采用常规的定位测量设备对步骤20132中标注的所有垂直地物的底部位置平面坐标进行实地测量,并将所有垂直地物的底部位置平面坐标实地测量结果存储至可移动存储设备内;步骤2053中进行垂直地物定位误差纠正时,所述数据处理器根据所述可移动存储设备内所存储的底部位置平面坐标实地测量结果,对被纠正垂直地物的顶部位置平面坐标进行纠正,并将被纠正的顶部位置平面坐标标注为被纠正垂直地物的底部位置实地测量结果;
206、按照步骤201至步骤205中所述的调绘方法,对多张所述单片InSAR影像分别进行调绘;
步骤三、接边处理:所述数据处理器按照常规数字地图的接边处理方法,对步骤二中调绘完成的多张所述单片InSAR影像进行接边处理,获得所述测区的完整InSAR正射影像调绘图。
上述基于单片InSAR正射影像的调绘方法,其特征是:步骤202中所述的对步骤20131中标注为未知地物的地物类型和性质进行实地确定后,还需获取所述未知地物的参考特征参数,并将所获取未知地物的参考特征参数存入所述InSAR影像地物特征库内,以对所述InSAR影像地物特征库进行完善;所述未知地物的参考特征参数包括经实地确定后获得的地物类型和性质以及所述未知地物在InSAR影像上的形状、大小尺寸、色调和灰度值。
上述基于单片InSAR正射影像的调绘方法,其特征是:步骤202中所述的到实地对各地物的地物性质进行核查确认过程中,还需采用常规定位测量设备对被判读单片InSAR影像上所有地面地物的初步定位结果进行核查确认,并根据核查确认结果对各地面地物的初步定位结果进行修改。
上述基于单片InSAR正射影像的调绘方法,其特征是:步骤2013中预先在存储单元内建立InSAR影像地物特征库时,采用试验区实地调绘法或同类测量资料综合分析法建立所述InSAR影像地物特征库;
当采用所述试验区实地调绘法时,先划定一个内部包括多种类型常见地物的试验区,之后对划定的试验区进行实地调绘,并根据实地调绘结果,对试验区内所包括的多种类型常见地物的参考特征参数同步进行记录;当采用所述同类测量资料综合分析法时,先对所述测区同期或近期内的航空影像、卫星影像和地形图进行收集,然后根据收集到的资料,对所述测区内所存在的各种地物的参考特征参数进行初步确定,之后到实地核查,进一步确认各种地物参考特征参数的正确性。
上述基于单片InSAR正射影像的调绘方法,其特征是:步骤Ⅱ中所述的当未匹配出与被判读地物的特征参数一致的参考特征参数时或者步骤ⅱ中所述的当不能目视判读得出被判读地物的类型和性质时,还需按照常规航空相片调绘中所用的逻辑推理法,且结合前期掌握的所述测区的地形与测区内各地物的信息资料,从被判读地物的空间分布规律及其与相邻地物之间的相互关系角度出发,对被判读地物的性质进行进一步人为判断;当经进一步判断得出被判读地物的性质时,数据处理器根据进一步判断结果且调用标注模块对被判读地物的性质进行标注;否则,数据处理器调用标注模块将被判读地物的性质标注为未知地物。
上述基于单片InSAR正射影像的调绘方法,其特征是:步骤一中所述的数据处理器和存储单元以及步骤二中所述的显示单元和参数输入单元组成一个完整的PC机。
上述基于单片InSAR正射影像的调绘方法,其特征是:步骤2011中对被判读单片InSAR影像进行放大时,放大倍数不小于2倍。
本发明与现有技术相比具有以下优点:
1、调绘方法设计合理、实现方便且调绘效率高,劳动强度低。
2、利用InSAR正射影像地物特征库能简单、快速判读地物的性质,同时能准确定位地物的位置,提高InSAR技术成图的精度和效率。由于InSAR技术的成像机理与光学系统不同,本发明所公开的室内判读方法根据InSAR影像的成像原理,能对InSAR正射影像进行准确、快速判读和定位。
3、调绘方法设计合理,主要包括建立InSAR正射影像地物特征库、室内判读、野外核查、特殊区域(具体指阴影区和叠掩区)的调绘与补测、新增地物与隐蔽地物的调绘与补测、室内数据处理和接边处理多个步骤。其中,InSAR正射影像地物特征库的建立过程是在室内影像判读标绘的基础上,再到实地进行一系列工作的综合过程;建立正射InSAR影像地物特征库的目的是由于InSAR影像不同于光学影像,不易被识别。因此,在进行室内判读之前,需要建立InSAR影像地物特征库对InSAR影像进行了解和掌握。而进行室内判读时,是在学习InSAR影像地物特征库,初步掌握InSAR影像的成像规律后,再根据地物在InSAR影像的直观特征(包括形状、大小、色调、灰度值等)在室内进行地物性质判断和位置标定的过程。实际进行地物性质判断,先调出InSAR影像地物特征库且根据InSAR影像地物的直观特征进行初步判断;再结合了解的有关背景资料,从地物空间分布规律、相互关系进一步判断,从而确定地物性质。进行野外核查是对室内判读的地物性质到实地进行核查确认,对于室内无法判读定性的地物进行野外调绘,并补充完善InSAR影像地物特征库;其次,对室内定位的地物进行野外确认,对室内无法定位的地物采用测量仪器进行野外定位和定性;最后,对需要改正高度引起的定位误差的垂直地物在野外量取地物高度。由于,本发明的调绘方法均是在单片InSAR正射影像上进行,最终提交调绘的数字成果,可直接供内业使用。
4、智能化程度高,且人工实际运算量小。
5、实用价值高,本发明所采用的调绘方法能基于单片InSAR影像快速获取地物的属性信息,并且能准确提取地物的位置坐标信息,因而能进一步推广机载合成孔径雷达干涉测量(InSAR)新技术在地形测绘领域的应用,充分发挥了InSAR新技术的各项优势,从而弥补了传统光学摄影测量的一些不足。由于InSAR系统不受天气影响,可以全天候、全天时的工作,这为灾害监测测绘、应急测绘提供了新的技术手段和方法,尤其在抗震救灾以及灾情监测中的地形测绘中发挥着重要作用。
6、本发明所公开的调绘方法形成基于单片InSAR影像外业调绘的技术要求,促进了InSAR技术的业务化运行,通过对InSAR技术以及外业调绘方法的研究,最终达到形成技术规范,形成业务化运行系统。
综上所述,本发明设计合理、方法步骤简单、实现方便且使用效果好、调绘精度高,能有效解决现有全数字航空摄影测量系统及调绘方法无法满足单片InSAR影像的调绘需要、单片InSAR影像调绘效果差、不能充分发挥InSAR技术的优势等多种缺陷和不足。
下面通过附图和实施例,对本发明的技术方案做进一步的详细描述。
附图说明
图1为本发明采用机载InSAR系统对垂直地物进行实际测量时产生阴影区和叠掩区的侧视成像原理示意图。
图2为本发明采用机载InSAR系统对垂直地物进行实际测量时产生定位误差的侧视成像原理示意图。
图3为本发明进行调绘时所采用调绘系统的电路原理框图。
图4为本发明进行调绘时的方法流程框图。
附图标记说明:
1—数据处理器;2—存储单元;3—显示单元;
4—参数输入单元。
具体实施方式
如图1、图2、图3及图4所示,本发明所述的基于单片InSAR正射影像的调绘方法,包括以下步骤:
步骤一、单片InSAR正射影像获取:采用机载InSAR系统且按常规的机载合成孔径雷达干涉测量方法,获取测区的多张单片InSAR影像,并将所获取的多张所述单片InSAR影像均存储至与数据处理器1相接的存储单元2内;多张所述单片InSAR影像均为正射影像。
步骤二、基于单片InSAR正射影像进行调绘:对步骤一中存储至存储单元2内的多张所述单片InSAR影像分别进行调绘,且对多张所述单片InSAR影像的调绘方法均相同;对于其中任一张单片InSAR影像而言,其调绘方法包括以下步骤:
201、室内判读,其室内判读过程如下:
2011、影像放大处理及同步显示:通过数据处理器1对被判读单片InSAR影像进行放大,并通过与数据处理器1相接的显示单元3对放大后的被判读单片InSAR影像进行同步显示。
本实施例中,对被判读单片InSAR影像进行放大时,放大倍数不小于2倍。
2012、阴影区、叠掩区和可判读区域确定:人为对通过显示单元3进行同步显示的被判读单片InSAR影像进行视觉观察,找出被判读单片InSAR影像上存在的阴影区和叠掩区;同时,通过数据处理器1对找出的阴影区和叠掩区进行标记;相应地,被判读单片InSAR影像上除被标记的阴影区和叠掩区之外的区域,为可判读区域。
2013、室内判读:利用数据处理器1且根据预先在存储单元2内所建立的InSAR影像地物特征库,对步骤2012中所述可判读区域内按照成图比例尺的技术规范要求需进行判读的所有地物分别进行判读;所述InSAR影像地物特征库内存储有多种不同类型地物的参考特征参数,且每一种类型地物的参考特征参数均包括地物类型、地物性质、地物在InSAR影像上的形状和大小尺寸以及地物在InSAR影像上的色调和灰度值;实际进行室内判读时,对所述可判读区域内需进行判读的所有地物的判读方法均相同;对于所述可判读区域内需进行判读的任一个地物来说,其室内判读过程包括以下步骤:
20131、地物性质判断:采用处理器自动判读方法或人为判读方法对被判读地物的性质进行室内判读。
当采用处理器自动判读方法对被判读地物的性质进行室内判读时,其室内判读过程如下:
Ⅰ、被判读地物的特征参数获取:所述数据处理器1调用尺寸量测模块和灰度值读取模块,在被判读单片InSAR影像上量测出被判读地物的大小尺寸且同步读取出被判读地物的灰度值;同时,根据被判读地物的视觉观察结果,人为通过与数据处理器1相接的参数输入单元4输入被判读地物在InSAR影像上的形状和在InSAR影像上的色调信息,此时数据处理器1获得被判读地物的特征参数并同步存储至存储单元2内;被判读地物的特征参数包括被判读地物在InSAR影像上的形状、大小尺寸以及色调信息和灰度值。
Ⅱ、特征匹配并确定被判读地物性质:所述数据处理器1调取地物特征参数匹配模块,将被判读地物的特征参数与所述InSAR影像地物特征库进行特征匹配:当匹配出与被判读地物的特征参数一致的参考特征参数时,匹配得出的参考特征参数中的地物类型和地物性质便为被判读地物的类型和性质,此时数据处理器1调用标注模块在被判读单片InSAR影像上对被判读地物的类型和性质进行标注;当未匹配出与被判读地物的特征参数一致的参考特征参数时,所述数据处理器1调用标注模块将被判读地物的性质标注为未知地物;
当采用人为判读方法对被判读地物的性质进行室内判读时,其室内判读过程如下:
i、InSAR影像地物特征库学习:作业人员通过对所述InSAR影像地物特征库内所存储的多种不同类型地物的参考特征参数进行学习,了解并掌握各种不同类型地物在InSAR影像上的目视判读特征,所述目视判读特征包括地物在InSAR影像上的形状、大小尺寸以及色调信息和灰度值。
实际对InSAR影像地物特征库进行学习时,所述目视判读特征为通过肉眼对进行InSAR影像观察,即可得出的地物在InSAR影像上的形状、大小尺寸以及色调信息和灰度值。
ⅱ、人为目视判读:作业人员结合步骤i中了解并掌握的各种不同类型地物在InSAR影像上的目视判读特征,对被判读地物进行目视判读,并相应判读出被判读地物的类型和性质;当目视判读得出被判读地物的类型和性质时,作业人员通过参数输入单元4输入目视判读结果,之后数据处理器1调用标注模块在被判读单片InSAR影像上对被判读地物的类型和性质进行标注。
当不能目视判读得出被判读地物的类型和性质时,作业人员通过参数输入单元4输入目视判读结果且此时目视判读结果为未知地物,之后数据处理器1调用标注模块在被判读单片InSAR影像上将被判读地物的性质标注为未知地物。
本实施例中,实际对地物性质进行判读时,步骤Ⅱ中所述的当未匹配出与被判读地物的特征参数一致的参考特征参数时或者步骤ⅱ中所述的当不能目视判读得出被判读地物的类型和性质时,还需按照常规航空相片调绘中所用的逻辑推理法,且结合前期掌握的所述测区的地形与测区内各地物的信息资料,从被判读地物的空间分布规律及其与相邻地物之间的相互关系角度出发,对被判读地物的性质进行进一步人为判断;当经进一步判断得出被判读地物的性质时,数据处理器1根据进一步判断结果且调用标注模块对被判读地物的性质进行标注;否则,数据处理器1调用标注模块将被判读地物的性质标注为未知地物。
所述逻辑推理法是指利用各种现象之间的关系,依照专业逻辑推理进行判读。例如河流是地壳垂直运动最敏感的标志,河流向一侧移动,说明另一侧可能有地壳上升,这就是利用地质、地貌专业的知识进行逻辑推理判读出来的。利用判读标志直接从像片上判读出来的地物多半是地面可见物体。对地理工作者来说,不仅要了解个体地物特征,更需要了解地区的综合特点,以及它们的发生发展规律。因此,逻辑推理方法在专业判读中应用相当广泛。所以,判读中应该注意运用地理专业知识,细致观察各种地物之间的相互关系,以取得更好的判读效果。
需注意的是:所述InSAR影像地物特征库为InSAR正射影像地物特征库,且每一种类型地物的参考特征参数均包括地物类型、地物性质、地物在InSAR正射影像上的形状和大小尺寸以及地物在InSAR正射影像上的色调和灰度值。
本实施例中,步骤20131中预先在存储单元2内建立InSAR影像地物特征库时,采用试验区实地调绘法或同类测量资料综合分析法建立所述InSAR影像地物特征库。
当采用所述试验区实地调绘法时,先划定一个内部包括多种类型常见地物的试验区,之后对划定的试验区进行实地调绘,并根据实地调绘结果,对试验区内所包括的多种类型常见地物的参考特征参数同步进行记录;当采用所述同类测量资料综合分析法时,先对所述测区同期或近期内的航空影像、卫星影像和地形图进行收集,然后根据收集到的资料,对所述测区内所存在的各种地物的参考特征参数进行初步确定,之后到实地核查,进一步确认各种地物参考特征参数的正确性。
无论采用上述哪一种建库方法建立InSAR影像地物特征库,均依据测区的具体情况而定。无论哪种建库方法建立InSAR影像地物特征库,都必须经过反复完善修改、补充的过程,不可能一次完成。建立InSAR影像地物特征库的目的是要快速准确识别各种地物,从而达到对地物进行准确定性的目的。
20132、地物定位:首先根据被判读地物在被判读单片InSAR影像上的阴影方向,判断出被判读地物的顶部位置,被判读地物的顶部位置与其阴影方向相反;其次,所述数据处理器1调用位置标定模块且按照常规数字地图标定方法对被判读地物的顶部位置进行标定,标定的同时数据处理器1调用平面坐标确定模块且根据预先在被判读单片InSAR影像上所选定参照物的平面坐标,确定被判读地物顶部位置的平面坐标并同步存储至存储单元2内。此时,存储单元2内所存储被判读地物顶部位置的平面坐标为被判读地物的初步定位结果。
202、野外核查:到实地对步骤20131中判读得出的可判读区域内各地物的地物性质进行核查确认,当核查出步骤20131中判读得出的地物性质错误时,作业人员根据野外核查结果对在被判读单片InSAR影像上所标注的错误地物性质进行修改;同时,对步骤20131中标注为未知地物的地物类型和性质进行实地确定,并根据实地确定结果在被判读单片InSAR影像上对所述未知地物的类型和性质进行标注。
本实施例中,步骤202中所述的对步骤20131中标注为未知地物的地物类型和性质进行实地确定后,还需获取所述未知地物的参考特征参数,并将所获取未知地物的参考特征参数存入所述InSAR影像地物特征库内,以对所述InSAR影像地物特征库进行完善;所述未知地物的参考特征参数包括经实地确定后获得的地物类型和性质以及所述未知地物在InSAR影像上的形状、大小尺寸、色调和灰度值。
203、阴影区和叠掩区的调绘与补测:采用常规定位测量设备且按照常规的测量和定位方法,对步骤2012中标记的阴影区和叠掩区分别进行调绘与补测,并将调绘与补测结果存储至可移动存储设备内。
204、新增地物与隐蔽地物的调绘与补测:按照常规外业像片调绘方法,对新增地物和隐蔽地物进行调绘与补测,并将调绘与补测结果存储至所述可移动存储设备内;所述新增地物为步骤一中所述机载合成孔径雷达干涉测量结束后且对所述测区的调绘工作结束前,在所述测区内新增设的地物;所述隐蔽地物为步骤一中进行机载合成孔径雷达干涉测量时因与周围相邻地物的InSAR影像混为一体且使得该地物无法辨识的地物,或者步骤一中进行机载合成孔径雷达干涉测量时被周围相邻的高大地物遮挡而使得该地物没有回波返回到机载InSAR系统雷达天线的地物。
205、室内数据处理,其数据处理过程包括以下步骤:
2051、野外调绘补测数据处理:将步骤203中和步骤204中所述可移动存储设备内存储的调绘与补测结果传送至数据处理器1,且数据处理器1调用数字展绘模块或数字转绘模块,将所述调绘与补测结果展绘或转绘至被判读单片InSAR影像上。
2052、垂直地物高度测定:采用基于单片InSAR影像的室内高度测量方法或基于野外高度测量的室内数据处理方法,对步骤20132中标注的所有垂直地物的高度进行测定,且对所有垂直地物的高度测定方法均相同。
实际进行调绘工程中,还需人为对所述可判读区域内需进行判读的所有地物是地面地物还是垂直地物进行判断,且判断各地物是地面地物还是垂直地物的方法均为常规判断方法,所述垂直地物为在测区地面的垂直方向上有高度的地物(也就是说,该地物的顶部位置和地面位置不一致),所述地面地物为在测区地面的垂直方向上没有高度的地物(也就是说,该地物的顶部位置和地面位置一致)。
对于任一个垂直地物来说,当采用基于单片InSAR影像的室内高度测量方法对被测量垂直地物高度进行测定时,数据处理器1调用垂直地物高度测量模块且根据公式h=Lgsinθcosθ或h=Lg1tanθ,计算得出被测量垂直地物的高度h;式中θ为步骤一中进行机载合成孔径雷达干涉测量时所述机载InSAR系统的侧视角,且侧视角θ通过参数输入单元4进行输入;Lg为量测得到的被纠正垂直地物在被判读单片InSAR影像上的阴影宽度,Lg1为量测得到的被纠正垂直地物在被判读单片InSAR影像上的叠掩宽度;之后,所述数据处理器1调用标注模块将被测量垂直地物的高度标注为h。
当采用基于野外高度测量的室内数据处理方法对被测量垂直地物高度进行测定时,步骤202中所述的到实地对各地物的地物性质进行核查确认过程中,还需采用常规的定位测量设备对被测量垂直地物的高度进行实地测量,并将被测量垂直地物的高度实地测量结果存储至可移动存储设备内;步骤2052中进行垂直地物高度测定时,根据所述可移动存储设备内所存储的高度实地测量结果,所述数据处理器1调用标注模块对被测量垂直地物的高度进行标注。
一般情况下,由于地面地物的高度为零,则无需对被判读单片InSAR影像上所存在地面地物的高度进行标注。
2053、垂直地物定位误差纠正:采用室内定位误差纠正方法或基于野外定位测量的室内数据纠正方法,对步骤20132中标注的所有垂直地物的定位误差进行纠正,且对所有垂直地物的定位误差纠正方法均相同。
对于任一个垂直地物来说,当采用室内定位误差纠正方法对被纠正垂直地物的定位误差进行纠正时,数据处理器1调用误差纠正模块且根据公式Δr=h/tanη,计算得出被纠正垂直地物因高度引起的顶部平面偏移量Δr,式中η为步骤一中进行机载合成孔径雷达干涉测量时所述机载InSAR系统的入射角;之后,数据处理器1根据计算得出的顶部平面偏移量Δr,对步骤20132中存储单元2内所存储被纠正垂直地物的初步定位结果进行纠正,获得纠正后的定位结果;然后,所述数据处理器1调用标注模块将被纠正垂直地物的底部位置平面坐标标注为纠正后的定位结果。
实际采用室内定位误差纠正方法进行纠正时,根据计算得出顶部平面偏移量Δr,对步骤20132中存储单元2内所存储的初步定位结果(即顶部位置的平面坐标(X1,Y1))进行纠正,获得纠正后的定位结果(即被纠正垂直地物的底部位置平面坐标(X1+Δx,Y1+Δy),其中Δr2=Δx2+Δy2)。
当采用基于野外定位测量的室内数据纠正方法对被纠正垂直地物的定位误差进行纠正时,步骤202中所述的到实地对各地物的地物性质进行核查确认过程中,还需采用常规的定位测量设备对步骤20132中标注的所有垂直地物的底部位置平面坐标进行实地测量,并将所有垂直地物的底部位置平面坐标实地测量结果存储至可移动存储设备内;步骤2053中进行垂直地物定位误差纠正时,所述数据处理器1根据所述可移动存储设备内所存储的底部位置平面坐标实地测量结果,对被纠正垂直地物的顶部位置平面坐标进行纠正,并将被纠正的顶部位置平面坐标标注为被纠正垂直地物的底部位置实地测量结果。
206、按照步骤201至步骤205中所述的调绘方法,对多张所述单片InSAR影像分别进行调绘。
步骤三、接边处理:所述数据处理器1按照常规数字地图的接边处理方法,对步骤二中调绘完成的多张所述单片InSAR影像进行接边处理,获得所述测区的完整InSAR正射影像调绘图。
本实施例中,步骤202中所述的到实地对各地物的地物性质进行核查确认过程中,还需采用常规定位测量设备对被判读单片InSAR影像上所有地面地物的初步定位结果进行核查确认,并根据核查确认结果对各地面地物的初步定位结果进行修改。
当采用本发明所述调绘方法对约16平方公里范围的测区进行调绘后进行了精度统计,同时该测区进行了实地调绘,并对本发明调绘结果与实地调绘结果进行对比,其对比分析结果见表1:
表1本发明调绘结果与实地调绘结果对比分析结果表
中误差:0.87米(规范要求:1:10000比例尺为5米)
由表1可看出,本发明调绘结果与实地调绘结果之间的中误差为0.87米,满足规范要求。
本实施例中,所述数据处理器1和存储单元2以及步骤二中所述的显示单元3和参数输入单元4组成一个完整的PC机,且数据处理器1、存储单元2、显示单元3和参数输入单元4组成调绘系统。
以上所述,仅是本发明的较佳实施例,并非对本发明作任何限制,凡是根据本发明技术实质对以上实施例所作的任何简单修改、变更以及等效结构变化,均仍属于本发明技术方案的保护范围内。
Claims (7)
1.一种基于单片InSAR正射影像的调绘方法,其特征在于该方法包括以下步骤:
步骤一、单片InSAR正射影像获取:采用机载InSAR系统且按常规的机载合成孔径雷达干涉测量方法,获取测区的多张单片InSAR影像,并将所获取的多张所述单片InSAR影像均存储至与数据处理器(1)相接的存储单元(2)内;多张所述单片InSAR影像均为正射影像;
步骤二、基于单片InSAR正射影像进行调绘:对步骤一中存储至存储单元(2)内的多张所述单片InSAR影像分别进行调绘,且对多张所述单片InSAR影像的调绘方法均相同;对于其中任一张单片InSAR影像而言,其调绘方法包括以下步骤:
201、室内判读,其室内判读过程如下:
2011、影像放大处理及同步显示:通过数据处理器(1)对被判读单片InSAR影像进行放大,并通过与数据处理器(1)相接的显示单元(3)对放大后的被判读单片InSAR影像进行同步显示;
2012、阴影区、叠掩区和可判读区域确定:人为对通过显示单元(3)进行同步显示的被判读单片InSAR影像进行视觉观察,找出被判读单片InSAR影像上存在的阴影区和叠掩区;同时,通过数据处理器(1)对找出的阴影区和叠掩区进行标记;相应地,被判读单片InSAR影像上除被标记的阴影区和叠掩区之外的区域,为可判读区域;
2013、室内判读:利用数据处理器(1)且根据预先在存储单元(2)内所建立的InSAR影像地物特征库,对步骤2012中所述可判读区域内按照成图比例尺的技术规范要求需进行判读的所有地物分别进行判读;所述InSAR影像地物特征库内存储有多种不同类型地物的参考特征参数,且每一种类型地物的参考特征参数均包括地物类型、地物性质、地物在InSAR影像上的形状和大小尺寸以及地物在InSAR影像上的色调和灰度值;实际进行室内判读时,对所述可判读区域内需进行判读的所有地物的判读方法均相同;对于所述可判读区域内需进行判读的任一个地物来说,其室内判读过程包括以下步骤:
20131、地物性质判断:采用处理器自动判读方法或人为判读方法对被判读地物的性质进行室内判读;
当采用处理器自动判读方法对被判读地物的性质进行室内判读时,其室内判读过程如下:
Ⅰ、被判读地物的特征参数获取:所述数据处理器(1)调用尺寸量测模块和灰度值读取模块,在被判读单片InSAR影像上量测出被判读地物的大小尺寸且同步读取出被判读地物的灰度值;同时,根据被判读地物的视觉观察结果,人为通过与数据处理器(1)相接的参数输入单元(4)输入被判读地物在InSAR影像上的形状和在InSAR影像上的色调信息,此时数据处理器(1)获得被判读地物的特征参数并同步存储至存储单元(2)内;被判读地物的特征参数包括被判读地物在InSAR影像上的形状、大小尺寸以及色调信息和灰度值;
Ⅱ、特征匹配并确定被判读地物性质:所述数据处理器(1)调取地物特征参数匹配模块,将被判读地物的特征参数与所述InSAR影像地物特征库进行特征匹配:当匹配出与被判读地物的特征参数一致的参考特征参数时,匹配得出的参考特征参数中的地物类型和地物性质便为被判读地物的类型和性质,此时数据处理器(1)调用标注模块在被判读单片InSAR影像上对被判读地物的类型和性质进行标注;当未匹配出与被判读地物的特征参数一致的参考特征参数时,所述数据处理器(1)调用标注模块将被判读地物的性质标注为未知地物;
当采用人为判读方法对被判读地物的性质进行室内判读时,其室内判读过程如下:
i、InSAR影像地物特征库学习:作业人员通过对所述InSAR影像地物特征库内所存储的多种不同类型地物的参考特征参数进行学习,了解并掌握各种不同类型地物在InSAR影像上的目视判读特征,所述目视判读特征包括地物在InSAR影像上的形状、大小尺寸以及色调信息和灰度值;
ⅱ、人为目视判读:作业人员结合步骤i中了解并掌握的各种不同类型地物在InSAR影像上的目视判读特征,对被判读地物进行目视判读,并相应判读出被判读地物的类型和性质;当目视判读得出被判读地物的类型和性质时,作业人员通过参数输入单元(4)输入目视判读结果,之后数据处理器(1)调用标注模块在被判读单片InSAR影像上对被判读地物的类型和性质进行标注;
当不能目视判读得出被判读地物的类型和性质时,作业人员通过参数输入单元(4)输入目视判读结果且此时目视判读结果为未知地物,之后数据处理器(1)调用标注模块在被判读单片InSAR影像上将被判读地物的性质标注为未知地物;
20132、地物定位:首先根据被判读地物在被判读单片InSAR影像上的阴影方向,判断出被判读地物的顶部位置,被判读地物的顶部位置与其阴影方向相反;其次,所述数据处理器(1)调用位置标定模块且按照常规数字地图标定方法对被判读地物的顶部位置进行标定,标定的同时数据处理器(1)调用平面坐标确定模块且根据预先在被判读单片InSAR影像上所选定参照物的平面坐标,确定被判读地物顶部位置的平面坐标并同步存储至存储单元(2)内;此时,存储单元(2)内所存储被判读地物顶部位置的平面坐标为被判读地物的初步定位结果;
202、野外核查:到实地对步骤20131中判读得出的可判读区域内各地物的地物性质进行核查确认,当核查出步骤20131中判读得出的地物性质错误时,作业人员根据野外核查结果对在被判读单片InSAR影像上所标注的错误地物性质进行修改;同时,对步骤20131中标注为未知地物的地物类型和性质进行实地确定,并根据实地确定结果在被判读单片InSAR影像上对所述未知地物的类型和性质进行标注;
203、阴影区和叠掩区的调绘与补测:采用常规定位测量设备且按照常规的测量和定位方法,对步骤2012中标记的阴影区和叠掩区分别进行调绘与补测,并将调绘与补测结果存储至可移动存储设备内;
204、新增地物与隐蔽地物的调绘与补测:按照常规外业像片调绘方法,对新增地物和隐蔽地物进行调绘与补测,并将调绘与补测结果存储至所述可移动存储设备内;所述新增地物为步骤一中所述机载合成孔径雷达干涉测量结束后且对所述测区的调绘工作结束前,在所述测区内新增设的地物;所述隐蔽地物为步骤一中进行机载合成孔径雷达干涉测量时因与周围相邻地物的InSAR影像混为一体且使得该地物无法辨识的地物,或者步骤一中进行机载合成孔径雷达干涉测量时被周围相邻的高大地物遮挡而使得该地物没有回波返回到机载InSAR系统雷达天线的地物;
205、室内数据处理,其数据处理过程包括以下步骤:
2051、野外调绘补测数据处理:将步骤203中和步骤204中所述可移动存储设备内存储的调绘与补测结果传送至数据处理器(1),且数据处理器(1)调用数字展绘模块或数字转绘模块,将所述调绘与补测结果展绘或转绘至被判读单片InSAR影像上;
2052、垂直地物高度测定:采用基于单片InSAR影像的室内高度测量方法或基于野外高度测量的室内数据处理方法,对步骤20132中标注的所有垂直地物的高度进行测定,且对所有垂直地物的高度测定方法均相同;所述垂直地物为在测区地面的垂直方向上有高度的地物;
对于任一个垂直地物来说,当采用基于单片InSAR影像的室内高度测量方法对被测量垂直地物高度进行测定时,数据处理器(1)调用垂直地物高度测量模块且根据公式h=Lgsinθcosθ或h=Lg1tanθ,计算得出被测量垂直地物的高度h;式中θ为步骤一中进行机载合成孔径雷达干涉测量时所述机载InSAR系统的侧视角,且侧视角θ通过参数输入单元(4)进行输入;Lg为量测得到的被纠正垂直地物在被判读单片InSAR影像上的阴影宽度,Lg1为量测得到的被纠正垂直地物在被判读单片InSAR影像上的叠掩宽度;之后,所述数据处理器(1)调用标注模块将被测量垂直地物的高度标注为h;
当采用基于野外高度测量的室内数据处理方法对被测量垂直地物高度进行测定时,步骤202中所述的到实地对各地物的地物性质进行核查确认过程中,还需采用常规的定位测量设备对被测量垂直地物的高度进行实地测量,并将被测量垂直地物的高度实地测量结果存储至可移动存储设备内;步骤2052中进行垂直地物高度测定时,根据所述可移动存储设备内所存储的高度实地测量结果,所述数据处理器(1)调用标注模块对被测量垂直地物的高度进行标注;
2053、垂直地物定位误差纠正:采用室内定位误差纠正方法或基于野外定位测量的室内数据纠正方法,对步骤20132中标注的所有垂直地物的定位误差进行纠正,且对所有垂直地物的定位误差纠正方法均相同;
对于任一个垂直地物来说,当采用室内定位误差纠正方法对被纠正垂直地物的定位误差进行纠正时,数据处理器(1)调用误差纠正模块且根据公式Δr=h/tanη,计算得出被纠正垂直地物因高度引起的顶部平面偏移量Δr,式中η为步骤一中进行机载合成孔径雷达干涉测量时所述机载InSAR系统的入射角,h为被纠正垂直地物的高度;之后,数据处理器(1)根据计算得出的顶部平面偏移量Δr,对步骤20132中存储单元(2)内所存储被纠正垂直地物的初步定位结果进行纠正,获得纠正后的定位结果;然后,所述数据处理器(1)调用标注模块将被纠正垂直地物的底部位置平面坐标标注为纠正后的定位结果;
当采用基于野外定位测量的室内数据纠正方法对被纠正垂直地物的定位误差进行纠正时,步骤202中所述的到实地对各地物的地物性质进行核查确认过程中,还需采用常规的定位测量设备对步骤20132中标注的所有垂直地物的底部位置平面坐标进行实地测量,并将所有垂直地物的底部位置平面坐标实地测量结果存储至可移动存储设备内;步骤2053中进行垂直地物定位误差纠正时,所述数据处理器(1)根据所述可移动存储设备内所存储的底部位置平面坐标实地测量结果,对被纠正垂直地物的顶部位置平面坐标进行纠正,并将被纠正的顶部位置平面坐标标注为被纠正垂直地物的底部位置实地测量结果;
206、按照步骤201至步骤205中所述的调绘方法,对多张所述单片InSAR影像分别进行调绘;
步骤三、接边处理:所述数据处理器(1)按照常规数字地图的接边处理方法,对步骤二中调绘完成的多张所述单片InSAR影像进行接边处理,获得所述测区的完整InSAR正射影像调绘图。
2.按照权利要求1所述的基于单片InSAR正射影像的调绘方法,其特征在于:步骤202中所述的对步骤20131中标注为未知地物的地物类型和性质进行实地确定后,还需获取所述未知地物的参考特征参数,并将所获取未知地物的参考特征参数存入所述InSAR影像地物特征库内,以对所述InSAR影像地物特征库进行完善;所述未知地物的参考特征参数包括经实地确定后获得的地物类型和性质以及所述未知地物在InSAR影像上的形状、大小尺寸、色调和灰度值。
3.按照权利要求1或2所述的基于单片InSAR正射影像的调绘方法,其特征在于:步骤202中所述的到实地对各地物的地物性质进行核查确认过程中,还需采用常规定位测量设备对被判读单片InSAR影像上所有地面地物的初步定位结果进行核查确认,并根据核查确认结果对各地面地物的初步定位结果进行修改。
4.按照权利要求1或2所述的基于单片InSAR正射影像的调绘方法,其特征在于:步骤2013中预先在存储单元(2)内建立InSAR影像地物特征库时,采用试验区实地调绘法或同类测量资料综合分析法建立所述InSAR影像地物特征库;
当采用所述试验区实地调绘法时,先划定一个内部包括多种类型常见地物的试验区,之后对划定的试验区进行实地调绘,并根据实地调绘结果,对试验区内所包括的多种类型常见地物的参考特征参数同步进行记录;当采用所述同类测量资料综合分析法时,先对所述测区同期或近期内的航空影像、卫星影像和地形图进行收集,然后根据收集到的资料,对所述测区内所存在的各种地物的参考特征参数进行初步确定,之后到实地核查,进一步确认各种地物参考特征参数的正确性。
5.按照权利要求1或2所述的基于单片InSAR正射影像的调绘方法,其特征在于:步骤Ⅱ中所述的当未匹配出与被判读地物的特征参数一致的参考特征参数时或者步骤ⅱ中所述的当不能目视判读得出被判读地物的类型和性质时,还需按照常规航空相片调绘中所用的逻辑推理法,且结合前期掌握的所述测区的地形与测区内各地物的信息资料,从被判读地物的空间分布规律及其与相邻地物之间的相互关系角度出发,对被判读地物的性质进行进一步人为判断;当经进一步判断得出被判读地物的性质时,数据处理器(1)根据进一步判断结果且调用标注模块对被判读地物的性质进行标注;否则,数据处理器(1)调用标注模块将被判读地物的性质标注为未知地物。
6.按照权利要求1或2所述的基于单片InSAR正射影像的调绘方法,其特征在于:步骤一中所述的数据处理器(1)和存储单元(2)以及步骤二中所述的显示单元(3)和参数输入单元(4)组成一个完整的PC机。
7.按照权利要求1或2所述的基于单片InSAR正射影像的调绘方法,其特征在于:步骤2011中对被判读单片InSAR影像进行放大时,放大倍数不小于2倍。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201110209941 CN102313537B (zh) | 2011-07-26 | 2011-07-26 | 基于单片InSAR正射影像的调绘方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201110209941 CN102313537B (zh) | 2011-07-26 | 2011-07-26 | 基于单片InSAR正射影像的调绘方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102313537A CN102313537A (zh) | 2012-01-11 |
CN102313537B true CN102313537B (zh) | 2013-01-30 |
Family
ID=45426920
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN 201110209941 Expired - Fee Related CN102313537B (zh) | 2011-07-26 | 2011-07-26 | 基于单片InSAR正射影像的调绘方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102313537B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102707284B (zh) * | 2012-06-25 | 2014-03-26 | 西安煤航信息产业有限公司 | 一种基于机载InSAR的地面控制测量布点方法 |
CN108957454B (zh) * | 2018-08-01 | 2020-10-27 | 中国科学院电子学研究所 | 机载干涉合成孔径雷达的干涉相位仿真方法 |
CN114612631B (zh) * | 2022-03-02 | 2023-06-09 | 自然资源部重庆测绘院 | 一种基于InSAR技术的高精度无漏洞DSM提取方法 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101551455B (zh) * | 2009-05-13 | 2011-10-19 | 西安电子科技大学 | 干涉合成孔径雷达三维地形成像系统及其高程测绘方法 |
CN201837388U (zh) * | 2010-10-26 | 2011-05-18 | 西安煤航信息产业有限公司 | 基于单片InSAR影像的垂直地物高度测量系统 |
CN201837717U (zh) * | 2010-10-26 | 2011-05-18 | 西安煤航信息产业有限公司 | 基于InSAR影像的垂直地物定位误差纠正系统 |
-
2011
- 2011-07-26 CN CN 201110209941 patent/CN102313537B/zh not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
CN102313537A (zh) | 2012-01-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103884321B (zh) | 一种遥感影像成图工艺 | |
CN102034238B (zh) | 基于光学成像测头和视觉图结构的多摄像机系统标定方法 | |
Gong et al. | ICEsat GLAS data for urban environment monitoring | |
CN105865334B (zh) | 地下建(构)筑物测绘方法 | |
CN104463868B (zh) | 一种基于无参数高分辨率影像的建筑高度快速获取方法 | |
CN104200086A (zh) | 宽基线可见光相机位姿估计方法 | |
CN102707284B (zh) | 一种基于机载InSAR的地面控制测量布点方法 | |
CN102980510B (zh) | 一种激光光尺影像测树装置及其测树方法 | |
CN105180890A (zh) | 融合激光点云和数字影像的岩体结构面产状测量方法 | |
CN103186892B (zh) | 利用航拍图像生成等比例实景现场图的方法及系统 | |
CN112833861A (zh) | 一种基于倾斜摄影大比例尺地形图的测绘方法和测绘系统 | |
CN111426303A (zh) | 一种岩溶坡立谷参数测量方法 | |
CN102313537B (zh) | 基于单片InSAR正射影像的调绘方法 | |
CN104729529B (zh) | 地形图测量系统误差判断的方法和系统 | |
CN109631754A (zh) | 一种测量装置坐标系标定的方法及相关装置 | |
Crispel et al. | All-sky photogrammetry techniques to georeference a cloud field | |
Remondino et al. | Surface reconstruction algorithms for detailed close-range object modeling | |
CN103792542A (zh) | 一种基于地面激光雷达技术的危岩落石调查方法 | |
CN109359631A (zh) | 一种基于卷积神经网络的海冰类型遥感分类方法 | |
CN109708570A (zh) | 用于掌子面结构面的信息采集与分析方法及装置 | |
CN106643667B (zh) | 测距方法和装置 | |
CN108195359A (zh) | 空间数据的采集方法及系统 | |
CN101644569B (zh) | 基于gps/ins数码摄影测量像控布点方法 | |
CN102997902B (zh) | 一种用于数码相机数据的快速质量检测方法 | |
CN114973017A (zh) | 融合遥感及gis技术的输变电设备高精度定位方法及系统 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
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: 20130130 Termination date: 20160726 |
|
CF01 | Termination of patent right due to non-payment of annual fee |