CN117664087B - 垂轨环扫式卫星影像核线生成方法、系统及设备 - Google Patents
垂轨环扫式卫星影像核线生成方法、系统及设备 Download PDFInfo
- Publication number
- CN117664087B CN117664087B CN202410132864.2A CN202410132864A CN117664087B CN 117664087 B CN117664087 B CN 117664087B CN 202410132864 A CN202410132864 A CN 202410132864A CN 117664087 B CN117664087 B CN 117664087B
- Authority
- CN
- China
- Prior art keywords
- image
- point
- epipolar line
- right image
- satellite
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 48
- 238000010606 normalization Methods 0.000 claims abstract description 16
- 238000012892 rational function Methods 0.000 claims description 21
- 238000010276 construction Methods 0.000 claims description 12
- 238000012545 processing Methods 0.000 abstract description 7
- 238000003384 imaging method Methods 0.000 description 9
- 230000008569 process Effects 0.000 description 9
- 238000004590 computer program Methods 0.000 description 8
- 238000010586 diagram Methods 0.000 description 7
- 230000003287 optical effect Effects 0.000 description 7
- 230000006870 function Effects 0.000 description 5
- 238000004891 communication Methods 0.000 description 4
- 238000011161 development Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 239000004973 liquid crystal related substance Substances 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000000644 propagated effect Effects 0.000 description 2
- 239000004065 semiconductor Substances 0.000 description 2
- 238000006467 substitution reaction Methods 0.000 description 2
- 230000001052 transient effect Effects 0.000 description 2
- 238000013459 approach Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000001444 catalytic combustion detection Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 239000000835 fiber Substances 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 239000013307 optical fiber Substances 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
- G01C11/00—Photogrammetry or videogrammetry, e.g. stereogrammetry; Photographic surveying
- G01C11/04—Interpretation of pictures
- G01C11/06—Interpretation of pictures by comparison of two or more pictures of the same area
- G01C11/08—Interpretation of pictures by comparison of two or more pictures of the same area the pictures not being supported in the same relative position as when they were taken
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T11/00—2D [Two Dimensional] image generation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/70—Determining position or orientation of objects or cameras
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T10/00—Road transport of goods or passengers
- Y02T10/10—Internal combustion engine [ICE] based vehicles
- Y02T10/40—Engine management systems
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Multimedia (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Image Processing (AREA)
Abstract
本发明属于新型卫星载荷影像处理技术领域,具体涉及垂轨环扫式卫星影像核线生成方法及系统,旨在解决目前垂轨环扫式卫星影像核线生成步骤复杂的问题。本发明包括:获取包含待测地区的垂轨环扫式卫星的下视遥感影像和前视遥感影像,作为立体像对的左右影像原始数据,并生成对应的RPC参数文件;基于垂轨环扫成像卫星下定义瞬时影像坐标系,在左右影像内建立瞬时影像坐标系;根据所述左右影像原始数据构建像方和物方关系;据像方和物方的关系,将所述立体像对的左右影像原始数据进行代入,通过联立求解获得每个扫描行对应的右像点归一化坐标;将每个扫描行对应的右像点归一化坐标去归一化并连线,获得左像点的核线,左像的同名点在该核线上。
Description
技术领域
本发明属于新型卫星载荷影像处理技术领域,具体涉及了一种垂轨环扫式卫星影像核线生成方法、系统及设备。
背景技术
根据核线的定义,左像上选取像点的同名点一定位于右像的核线上,据此可将二维的影像匹配转变为沿核线的一维匹配,利用核线模型寻找同一立体像对中的同名点能大大减小搜索范围,提高搜索效率。
随着卫星研制技术、光学传感器姿态控制技术的进步与发展,区别于线阵推扫成像的垂轨环扫成像机制被提出。垂轨环扫卫星上搭载的多条线阵CCD沿卫星飞行方向安置,光学传感器绕飞行方向360°旋转,视场在垂直轨道方向上不断叠加以实现环扫成像。垂轨环扫卫星兼具高分辨率和超大幅宽的特点,在轨道高度500km的情况下,能实现幅宽3000公里,星下点分辨率1m的指标。
卫星搭载双相机,垂直轨道搭载高分辨率相机,成像为下视,多功能相机主光轴沿轨向前倾斜一定角度,成像为前视,两台相机在地球径向方向上夹角180°,在环扫过程中双相机交替成像。这就使得前视和下视之间同名像点的核线模型极其复杂,和线阵推扫卫星核线模型、面阵式中心投影核线模型大不相同。
发明内容
为了解决现有技术中的上述问题,即目前垂轨环扫式卫星影像核线生成步骤复杂的问题,本发明提供了一种垂轨环扫式卫星影像核线生成方法,包括:
步骤S100,获取包含待测地区的垂轨环扫式卫星的下视遥感影像和前视遥感影像,作为立体像对的左右影像原始数据,并生成立体像对的左右影像原始数据对应的RPC参数文件;
步骤S200,基于所述垂轨环扫式卫星的下视遥感影像和前视遥感影像定义瞬时影像坐标系,在立体像对的左右影像内建立瞬时影像坐标系;
由于垂轨环扫式卫星区别于线阵推扫式卫星的成像方式,导致相机姿态角在不断地变化,所以需要重新定义其瞬时影像坐标系。
步骤S300,根据所述立体像对的左右影像原始数据构建像方和物方的关系;
步骤S400,根据像方和物方的关系,将所述立体像对的左右影像原始数据进行归一化后代入所述像方和物方的关系,求解获得每个扫描行对应的右像点归一化坐标;
由于垂轨环扫遥感影像由上万行扫描行组成,需要在所有扫描行上求得所有可能的同名点,以此生成核线。
步骤S500,将每个扫描行对应的右像点归一化坐标去归一化并连线,获得左像点的核线,左像的同名点在该核线上;
所述同名点为同一物方在不同遥感影像上对应的像点;
所述核线为同名点轨迹,左像点在右像上的同名点就在核线上,在核线上搜索同名点,将二维搜索变为了一维搜索。
进一步的,所述立体像对的左右影像原始数据,具体包括:
包括立体像对的左影像原始数据和立体像对的右影像原始数据;
所述立体像对的左影像原始数据和立体像对的右影像原始数据中包括原始控制点数据、传感器的姿态数据和遥感影像的内外方位元素。
进一步的,所述解算RPC参数文件,具体包括:
根据所述立体像对的左右影像原始数据构建虚拟控制格网;
根据构建的虚拟控制格网进行RPC模型拟合,解算出RPC参数。
进一步的,所述瞬时影像坐标系,具体包括:
瞬时影像坐标系原点规定为每个扫描行的中心,垂轨环扫式卫星环扫的方向为瞬时影像坐标系的y轴,垂轨环扫式卫星运行的方向为瞬时影像坐标系的y轴。
进一步的,所述步骤S300,具体包括:
根据所述立体像对的左右影像原始数据建立左像归一化坐标与物方归一化坐标之间的有理函数模型:
;
根据所述立体像对的左右影像原始数据建立右像归一化坐标与物方归一化坐标之间的有理函数模型:
;
其中,为物方点地理坐标的归一化坐标,/>和/>为左像点像方归一化坐标,/>和/>为右像点像方归一化坐标;
多项式、/>、/>和/>均为不超过3次的多项式。
进一步的,所述步骤S400,具体包括:
在左像范围内选取一个像点,将像点坐标进行归一化处理后得到归一化左像点坐标;将归一化左像点坐标代入所述左像归一化坐标与物方归一化坐标之间的有理函数模型中;
遍历右像的每一个扫描行,将右像的每一个扫描行的y坐标进行归一化处理后依次代入所述右像归一化坐标与物方归一化坐标之间的有理函数模型中;
将所述归一化左像点坐标和归一化右像坐标后的有理函数模型和有理函数模型/>联立求解,获得每个扫描行对应的x坐标归一化值;
所述每个扫描行对应的x坐标归一化值为右像每个扫描行对应的右像点归一化坐标。
本发明的另一方面,提出了一种垂轨环扫式卫星影像核线生成系统,包括原始数据获取模块、瞬时影像坐标系构建模块、像方物方关系构建模块、统计输出模块:
所述原始数据获取模块,用于获取包含待测地区的垂轨环扫式卫星的下视遥感影像和前视遥感影像,作为立体像对的左右影像原始数据,并生成立体像对的左右影像原始数据对应的RPC参数文件;
所述瞬时影像坐标系构建模块,用于基于所述垂轨环扫式卫星的下视遥感影像和前视遥感影像定义瞬时影像坐标系,在立体像对的左右影像内建立瞬时影像坐标系;
所述像方物方关系构建模块,用于根据所述立体像对的左右影像原始数据构建像方和物方的关系;
所述统计输出模块,用于根据像方和物方的关系,通过联立求解获得左像点的核线,左像的同名点即在该核线上。
进一步的,所述统计输出模块,具体包括:
联立求解单元,用于根据像方和物方的关系,将所述立体像对的左右影像原始数据进行代入,求解获得每个扫描行对应的右像点归一化坐标;
核线构成单元,用于将每个扫描行对应的右像点归一化坐标去归一化并连线,获得左像点的核线,左像的同名点即在该核线上;
本发明的第三方面,提出了一种电子设备,包括:
至少一个处理器;以及
与至少一个所述处理器通信连接的存储器;其中,
所述存储器存储有可被所述处理器执行的指令,所述指令用于被所述处理器执行以实现上述的一种垂轨环扫式卫星影像核线生成方法。
本发明的第四方面,提出了一种计算机可读存储介质,所述计算机可读存储介质存储有计算机指令,所述计算机指令用于被所述计算机执行以实现上述的一种垂轨环扫式卫星影像核线生成方法。
本发明的有益效果:
(1)本发明通过联立有理函数模型,省去了获取高程的步骤,仅通过左像与物方和右像与物方的有理函数联立便可生成卫星影像核线;
(2)本发明基于有理函数模型,构建物方和像方的关系,此方法求解的核线不仅更接近于真实核线,而且求解模型简单可靠。
附图说明
通过阅读参照以下附图所作的对非限制性实施例所作的详细描述,本申请的其它特征、目的和优点将会变得更明显:
图1是本发明一种垂轨环扫式卫星影像核线生成方法的流程图;
图2是本发明一种垂轨环扫式卫星影像核线生成方法的垂轨环扫式卫星相机环扫示意图;
图3是本发明一种垂轨环扫式卫星影像核线生成方法基于有理函数模型的核线示意图;
图4是于实现本申请方法、系统、设备实施例的服务器的计算机系统的结构示意图。
具体实施方式
下面结合附图和实施例对本申请作进一步的详细说明。可以理解的是,此处所描述的具体实施例仅用于解释相关发明,而非对该发明的限定。另外还需要说明的是,为了便于描述,附图中仅示出了与有关发明相关的部分。
需要说明的是,在不冲突的情况下,本申请中的实施例及实施例中的特征可以相互组合。下面将参考附图并结合实施例来详细说明本申请。
本发明提供一种垂轨环扫式卫星影像核线生成方法,本方法包括:
步骤S100,获取包含待测地区的垂轨环扫式卫星的下视遥感影像和前视遥感影像,作为立体像对的左右影像原始数据,并生成立体像对的左右影像原始数据对应的RPC参数文件;
步骤S200,基于垂轨环扫成像卫星下定义瞬时影像坐标系,在立体像对的左右影像内建立瞬时影像坐标系;
步骤S300,根据所述立体像对的左右影像原始数据构建像方和物方的关系;
步骤S400,根据像方和物方的关系,将所述立体像对的左右影像原始数据进行归一化后代入所述像方和物方的关系,求解获得每个扫描行对应的右像点归一化坐标;
步骤S500,将每个扫描行对应的右像点归一化坐标去归一化并连线,获得左像点的核线,左像的同名点即在该核线上;
所述同名点为同一物方在不同遥感影像上对应的像点;
所述核线为同名点轨迹,左像点在右像上的同名点就在核线上,在核线上搜索同名点,将二维搜索变为了一维搜索。
为了更清晰地对本发明一种垂轨环扫式卫星影像核线生成方法进行说明,下面结合图1对本发明实施例中各步骤展开详述。
本发明第一实施例的一种垂轨环扫式卫星影像核线生成方法,如图1所示,各步骤详细描述如下:
步骤S100,如图2所示,获取包含待测地区的垂轨环扫式卫星的下视遥感影像和前视遥感影像,作为立体像对的左右影像原始数据,并生成立体像对的左右影像原始数据对应的RPC参数文件;
下视和前视遥感影像既可以是同轨影像,也可以是异轨影像,只需要覆盖同一地区即可。卫星下传的数据为相机获取影像时的姿轨数据,需要据此解算出相应的RPC参数文件。
在本实施例中,所述立体像对的左右影像原始数据,具体包括:
立体像对的左影像原始数据、立体像对的右影像原始数据、原始控制点数据、传感器的姿态数据和遥感影像的内外方位元素;
步骤S200,基于垂轨环扫成像卫星下定义瞬时影像坐标系,在立体像对的左右影像内建立瞬时影像坐标系;
在本实施例中,所述瞬时影像坐标系,具体包括:
瞬时影像坐标系原点规定为每个扫描行的中心,垂轨环扫式卫星环扫的方向为瞬时影像坐标系的y轴,垂轨环扫式卫星运行的方向为瞬时影像坐标系的y轴,由于每扫描行姿态均不相同,对于该坐标系,其y值恒为0。
步骤S300,根据所述立体像对的左右影像原始数据构建像方和物方的关系;
在本实施例中,所述步骤S300,具体包括:
通过比值多项式将物方点地理坐标与左像坐标和右像坐标关系起来,构成有理函数模型;
为了增强参数求解的稳定性,将物方点地理坐标与左像坐标和右像坐标归一化处理;
对于左像,根据所述立体像对的左右影像原始数据建立左像归一化坐标与物方归一化坐标之间的有理函数模型:
;
对于右像,根据所述立体像对的左右影像原始数据建立右像归一化坐标与物方归一化坐标之间的有理函数模型:
;
多项式、/>、/>、/>均为不超过3次的多项式,其通用形式为:
;
其中,为物方点地理坐标的归一化坐标,/>和/>为左像点像方归一化坐标,/>和/>为右像点像方归一化坐标;
在本实施例中,所述解算RPC参数文件,具体包括:
根据所述立体像对的左右影像原始数据构建虚拟控制格网;
根据构建的虚拟控制格网进行RPC模型拟合,解算出RPC参数。
步骤S400,根据像方和物方的关系,将所述立体像对的左右影像原始数据进行归一化后代入所述像方和物方的关系,求解获得每个扫描行对应的右像点归一化坐标;
在本实施例中,所述步骤S400,如图3所示,具体包括:
在左像范围内选取一个像点,将像点坐标进行归一化处理后得到归一化左像点坐标;将归一化左像点坐标代入所述左像归一化坐标与物方归一化坐标之间的有理函数模型中;
遍历右像的每一个扫描行,将右像的每一个扫描行的y坐标进行归一化处理后依次代入所述右像归一化坐标与物方归一化坐标之间的有理函数模型中;联立求解非线性方程组,获得每个扫描行对应的x坐标归一化值;
所述每个扫描行对应的x坐标归一化值为右像每个扫描行对应的右像点归一化坐标。
步骤S500,将每个扫描行对应的右像点归一化坐标去归一化并连线,获得左像点的核线,左像的同名点在该核线上;
所述同名点为同一物方在不同遥感影像上对应的像点;
所述核线为同名点轨迹,左像点在右像上的同名点就在核线上,在核线上搜索同名点,将二维搜索变为了一维搜索。
上述实施例中虽然将各个步骤按照上述先后次序的方式进行了描述,但是本领域技术人员可以理解,为了实现本实施例的效果,不同的步骤之间不必按照这样的次序执行,其可以同时(并行)执行或以颠倒的次序执行,这些简单的变化都在本发明的保护范围之内。
本发明第二实施例的一种垂轨环扫式卫星影像核线生成系统,包括原始数据获取模块、瞬时影像坐标系构建模块、像方物方关系构建模块、统计输出模块:
所述原始数据获取模块,用于获取包含待测地区的垂轨环扫式卫星的下视遥感影像和前视遥感影像,作为立体像对的左右影像原始数据,并生成立体像对的左右影像原始数据对应的RPC参数文件;
所述瞬时影像坐标系构建模块,用于基于垂轨环扫成像卫星下定义瞬时影像坐标系,在立体像对的左右影像内建立瞬时影像坐标系;
所述像方物方关系构建模块,用于根据所述立体像对的左右影像原始数据构建像方和物方的关系;
所述统计输出模块,用于根据像方和物方的关系,通过联立求解获得左像点的核线,左像的同名点即在该核线上。
在本实施例中,所述统计输出模块,具体包括:
联立求解单元,用于根据像方和物方的关系,将所述立体像对的左右影像原始数据进行代入,通过联立求解获得每个扫描行对应的右像点归一化坐标;
核线构成单元,用于将每个扫描行对应的右像点归一化坐标去归一化并连线,获得左像点的核线,左像的同名点即在该核线上。
所属技术领域的技术人员可以清楚地了解到,为描述的方便和简洁,上述描述的系统的具体工作过程及有关说明,可以参考前述方法实施例中的对应过程,在此不再赘述。
需要说明的是,上述实施例提供的一种垂轨环扫式卫星影像核线生成系统,仅以上述各功能模块的划分进行举例说明,在实际应用中,可以根据需要而将上述功能分配由不同的功能模块来完成,即将本发明实施例中的模块或者步骤再分解或者组合,例如,上述实施例的模块可以合并为一个模块,也可以进一步拆分成多个子模块,以完成以上描述的全部或者部分功能。对于本发明实施例中涉及的模块、步骤的名称,仅仅是为了区分各个模块或者步骤,不视为对本发明的不当限定。
本发明第三实施例的一种电子设备,包括:
至少一个处理器;以及
与至少一个所述处理器通信连接的存储器;其中,
所述存储器存储有可被所述处理器执行的指令,所述指令用于被所述处理器执行以实现上述的一种垂轨环扫式卫星影像核线生成方法。
本发明第四实施例的一种计算机可读存储介质,所述计算机可读存储介质存储有计算机指令,所述计算机指令用于被所述计算机执行以实现上述的一种垂轨环扫式卫星影像核线生成方法。
所属技术领域的技术人员可以清楚地了解到,为描述的方便和简洁,上述描述的存储装置、处理装置的具体工作过程及有关说明,可以参考前述方法实施例中的对应过程,在此不再赘述。
本领域技术人员应该能够意识到,结合本文中所公开的实施例描述的各示例的模块、方法步骤,能够以电子硬件、计算机软件或者二者的结合来实现,软件模块、方法步骤对应的程序可以置于随机存储器(RAM)、内存、只读存储器(ROM)、电可编程ROM、电可擦除可编程ROM、寄存器、硬盘、可移动磁盘、CD-ROM、或技术领域内所公知的任意其它形式的存储介质中。为了清楚地说明电子硬件和软件的可互换性,在上述说明中已经按照功能一般性地描述了各示例的组成及步骤。这些功能究竟以电子硬件还是软件方式来执行,取决于技术方案的特定应用和设计约束条件。本领域技术人员可以对每个特定的应用来使用不同方法来实现所描述的功能,但是这种实现不应认为超出本发明的范围。
下面参考图4其示出了用于实现本申请方法、系统、设备实施例的服务器的计算机系统的结构示意图。图4出的服务器仅仅是一个示例,不应对本申请实施例的功能和使用范围带来任何限制。
如图4示,计算机系统包括中央处理单元(CPU,Central Processing Unit)601,其可以根据存储在只读存储器(ROM,Read Only Memory)602中的程序或者从存储部分608加载到随机访问存储器(RAM,Random Access Memory)603中的程序而执行各种适当的动作和处理。在RAM 603中,还存储有系统操作所需的各种程序和数据。CPU 601、ROM 602以及RAM603通过总线604彼此相连。输入/输出(I/O,Input/Output)接口605也连接至总线604。
以下部件连接至I/O接口605:包括键盘、鼠标等的输入部分606;包括诸如阴极射线管(CRT,Cathode Ray Tube)、液晶显示器(LCD,Liquid Crystal Display)等以及扬声器等的输出部分607;包括硬盘等的存储部分608;以及包括诸如LAN(局域网,Local AreaNetwork)卡、调制解调器等的网络接口卡的通信部分609。通信部分609经由诸如因特网的网络执行通信处理。驱动器610也根据需要连接至I/O接口605。可拆卸介质611,诸如磁盘、光盘、磁光盘、半导体存储器等等,根据需要安装在驱动器610上,以便于从其上读出的计算机程序根据需要被安装入存储部分608。
特别地,根据本公开的实施例,上文参考流程图描述的过程可以被实现为计算机软件程序。例如,本公开的实施例包括一种计算机程序产品,其包括承载在计算机可读介质上的计算机程序,该计算机程序包含用于执行流程图所示的方法的程序代码。在这样的实施例中,该计算机程序可以通过通信部分609从网络上被下载和安装,和/或从可拆卸介质611被安装。在该计算机程序被中央处理单元(CPU)601执行时,执行本申请的方法中限定的上述功能。需要说明的是,本申请上述的计算机可读介质可以是计算机可读信号介质或者计算机可读存储介质或者是上述两者的任意组合。计算机可读存储介质例如可以是——但不限于——电、磁、光、电磁、红外线、或半导体的系统、装置或器件,或者任意以上的组合。计算机可读存储介质的更具体的例子可以包括但不限于:具有一个或多个导线的电连接、便携式计算机磁盘、硬盘、随机访问存储器(RAM)、只读存储器(ROM)、可擦式可编程只读存储器(EPROM或闪存)、光纤、便携式紧凑磁盘只读存储器(CD-ROM)、光存储器件、磁存储器件、或者上述的任意合适的组合。在本申请中,计算机可读存储介质可以是任何包含或存储程序的有形介质,该程序可以被指令执行系统、装置或者器件使用或者与其结合使用。而在本申请中,计算机可读的信号介质可以包括在基带中或者作为载波一部分传播的数据信号,其中承载了计算机可读的程序代码。这种传播的数据信号可以采用多种形式,包括但不限于电磁信号、光信号或上述的任意合适的组合。计算机可读的信号介质还可以是计算机可读存储介质以外的任何计算机可读介质,该计算机可读介质可以发送、传播或者传输用于由指令执行系统、装置或者器件使用或者与其结合使用的程序。计算机可读介质上包含的程序代码可以用任何适当的介质传输,包括但不限于:无线、电线、光缆、RF等等,或者上述的任意合适的组合。
可以以一种或多种程序设计语言或其组合来编写用于执行本申请的操作的计算机程序代码,上述程序设计语言包括面向对象的程序设计语言—诸如Java、Smalltalk、C++,还包括常规的过程式程序设计语言—诸如“C”语言或类似的程序设计语言。程序代码可以完全地在用户计算机上执行、部分地在用户计算机上执行、作为一个独立的软件包执行、部分在用户计算机上部分在远程计算机上执行、或者完全在远程计算机或服务器上执行。在涉及远程计算机的情形中,远程计算机可以通过任意种类的网络——包括局域网(LAN)或广域网(WAN)—连接到用户计算机,或者,可以连接到外部计算机(例如利用因特网服务提供商来通过因特网连接)。
附图中的流程图和框图,图示了按照本申请各种实施例的系统、方法和计算机程序产品的可能实现的体系架构、功能和操作。在这点上,流程图或框图中的每个方框可以代表一个模块、程序段、或代码的一部分,该模块、程序段、或代码的一部分包含一个或多个用于实现规定的逻辑功能的可执行指令。也应当注意,在有些作为替换的实现中,方框中所标注的功能也可以以不同于附图中所标注的顺序发生。例如,两个接连地表示的方框实际上可以基本并行地执行,它们有时也可以按相反的顺序执行,这依所涉及的功能而定。也要注意的是,框图和/或流程图中的每个方框、以及框图和/或流程图中的方框的组合,可以用执行规定的功能或操作的专用的基于硬件的系统来实现,或者可以用专用硬件与计算机指令的组合来实现。
术语“第一”、“第二”等是用于区别类似的对象,而不是用于描述或表示特定的顺序或先后次序。
术语“包括”或者任何其它类似用语旨在涵盖非排他性的包含,从而使得包括一系列要素的过程、方法、物品或者设备/装置不仅包括那些要素,而且还包括没有明确列出的其它要素,或者还包括这些过程、方法、物品或者设备/装置所固有的要素。
至此,已经结合附图所示的优选实施方式描述了本发明的技术方案,但是,本领域技术人员容易理解的是,本发明的保护范围显然不局限于这些具体实施方式。在不偏离本发明的原理的前提下,本领域技术人员可以对相关技术特征做出等同的更改或替换,这些更改或替换之后的技术方案都将落入本发明的保护范围之内。
Claims (8)
1.一种垂轨环扫式卫星影像核线生成方法,其特征在于,所述方法包括:
步骤S100,获取包含待测地区的垂轨环扫式卫星的下视遥感影像和前视遥感影像,分别作为立体像对的左右影像原始数据,并生成立体像对的左右影像原始数据对应的RPC参数文件;
步骤S200,基于所述垂轨环扫式卫星的下视遥感影像和前视遥感影像定义瞬时影像坐标系,在立体像对的左右影像内建立瞬时影像坐标系;
步骤S300,根据所述立体像对的左右影像原始数据构建像方和物方的关系;
根据所述立体像对的左右影像原始数据建立左像归一化坐标与物方归一化坐标之间的有理函数模型:
;
根据所述立体像对的左右影像原始数据建立右像归一化坐标与物方归一化坐标之间的有理函数模型:
;
其中,为物方点地理坐标的归一化坐标,/>和/>为左像点像方归一化坐标,和/>为右像点像方归一化坐标;
多项式、/>、/>和/>均为不超过3次的多项式;
步骤S400,根据像方和物方的关系,将所述立体像对的左右影像原始数据进行归一化后代入所述像方和物方的关系,求解获得每个扫描行对应的右像点归一化坐标;
在左像范围内选取一个像点,将像点坐标进行归一化处理后得到归一化左像点坐标;将归一化左像点坐标代入所述左像归一化坐标与物方归一化坐标之间的有理函数模型中;
遍历右像的每一个扫描行,将右像的每一个扫描行的y坐标进行归一化处理后代入所述右像归一化坐标与物方归一化坐标之间的有理函数模型中;
将所述归一化左像点坐标和归一化右像坐标后的有理函数模型和有理函数模型联立求解,获得每个扫描行对应的x坐标归一化值;
所述每个扫描行对应的x坐标归一化值为右像每个扫描行对应的右像点归一化坐标;
步骤S500,将每个扫描行对应的右像点归一化坐标去归一化并连线,获得左像点的核线,左像的同名点在该核线上;
所述同名点为同一物方在不同遥感影像上对应的像点;
所述核线为同名点轨迹,左像点在右像上的同名点就在核线上,在核线上搜索同名点,将二维搜索变为了一维搜索。
2.根据权利要求1所述的一种垂轨环扫式卫星影像核线生成方法,其特征在于,所述立体像对的左右影像原始数据,具体包括:
包括立体像对的左影像原始数据和立体像对的右影像原始数据;
所述立体像对的左影像原始数据和立体像对的右影像原始数据中包括原始控制点数据、传感器的姿态数据和遥感影像的内外方位元素。
3.根据权利要求1所述的一种垂轨环扫式卫星影像核线生成方法,其特征在于,所述RPC参数文件,获取过程包括:
根据所述立体像对的左右影像原始数据构建虚拟控制格网;
根据构建的虚拟控制格网进行RPC模型拟合,解算出RPC参数。
4.根据权利要求1所述的一种垂轨环扫式卫星影像核线生成方法,其特征在于,所述瞬时影像坐标系,具体包括:
瞬时影像坐标系原点规定为每个扫描行的中心,垂轨环扫式卫星环扫的方向为瞬时影像坐标系的y轴,垂轨环扫式卫星运行的方向为瞬时影像坐标系的y轴。
5.一种垂轨环扫式卫星影像核线生成系统,其特征在于,基于权利要求1-4任一项所述的垂轨环扫式卫星影像核线生成方法,所述系统包括原始数据获取模块、瞬时影像坐标系构建模块、像方物方关系构建模块、统计输出模块:
所述原始数据获取模块,用于获取包含待测地区的垂轨环扫式卫星的下视遥感影像和前视遥感影像,作为立体像对的左右影像原始数据,并生成立体像对的左右影像原始数据对应的RPC参数文件;
所述瞬时影像坐标系构建模块,用于基于所述垂轨环扫式卫星的下视遥感影像和前视遥感影像定义瞬时影像坐标系,在立体像对的左右影像内建立瞬时影像坐标系;
所述像方物方关系构建模块,用于根据所述立体像对的左右影像原始数据构建像方和物方的关系;
所述统计输出模块,用于根据像方和物方的关系,通过联立求解获得左像点的核线,左像的同名点即在该核线上。
6.根据权利要求5所述的一种垂轨环扫式卫星影像核线生成系统,其特征在于,所述统计输出模块,具体包括:
联立求解单元,用于根据像方和物方的关系,将所述立体像对的左右影像原始数据进行代入,求解获得每个扫描行对应的右像点归一化坐标;
核线构成单元,用于将每个扫描行对应的右像点归一化坐标去归一化并连线,获得左像点的核线,左像的同名点即在该核线上。
7.一种电子设备,其特征在于,包括:
至少一个处理器;以及
与至少一个所述处理器通信连接的存储器;其中,
所述存储器存储有可被所述处理器执行的指令,所述指令用于被所述处理器执行以实现权利要求1-4任一项所述的一种垂轨环扫式卫星影像核线生成方法。
8.一种计算机可读存储介质,其特征在于,所述计算机可读存储介质存储有计算机指令,所述计算机指令用于被所述计算机执行以实现权利要求1-4任一项所述的一种垂轨环扫式卫星影像核线生成方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202410132864.2A CN117664087B (zh) | 2024-01-31 | 2024-01-31 | 垂轨环扫式卫星影像核线生成方法、系统及设备 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202410132864.2A CN117664087B (zh) | 2024-01-31 | 2024-01-31 | 垂轨环扫式卫星影像核线生成方法、系统及设备 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN117664087A CN117664087A (zh) | 2024-03-08 |
CN117664087B true CN117664087B (zh) | 2024-04-02 |
Family
ID=90071673
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202410132864.2A Active CN117664087B (zh) | 2024-01-31 | 2024-01-31 | 垂轨环扫式卫星影像核线生成方法、系统及设备 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN117664087B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117928494B (zh) * | 2024-03-18 | 2024-05-24 | 中国人民解放军战略支援部队航天工程大学 | 一种光学卫星切片影像的几何定位测量方法、系统和设备 |
Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101799290A (zh) * | 2010-02-09 | 2010-08-11 | 武汉大学 | 基于高程基准面的卫星立体影像像对近似核线生成方法以及装置 |
CN102213762A (zh) * | 2011-04-12 | 2011-10-12 | 中交第二公路勘察设计研究院有限公司 | 基于rfm模型的多源星载sar影像自动匹配方法 |
CN103630120A (zh) * | 2013-07-16 | 2014-03-12 | 中国人民解放军信息工程大学 | 基于严密几何模型的火星表面线阵影像核线重采样方法 |
CN104236527A (zh) * | 2014-09-15 | 2014-12-24 | 中国人民解放军总参谋部测绘研究所 | 一种基于投影轨迹法的线阵卫星遥感影像核线特性定量分析方法 |
CN109829939A (zh) * | 2019-01-18 | 2019-05-31 | 南京泛在地理信息产业研究院有限公司 | 一种缩小多视影像匹配同名像点搜索范围的方法 |
CN110135455A (zh) * | 2019-04-08 | 2019-08-16 | 平安科技(深圳)有限公司 | 影像匹配方法、装置及计算机可读存储介质 |
CN110853140A (zh) * | 2019-10-11 | 2020-02-28 | 北京空间机电研究所 | 一种dem辅助的光学视频卫星影像稳像方法 |
CN111667533A (zh) * | 2019-03-08 | 2020-09-15 | 西安中科星图空间数据技术有限公司 | 一种基于有理函数模型获取卫星严格成像模型参数的方法 |
CN114332085A (zh) * | 2022-03-11 | 2022-04-12 | 西安中科西光航天科技有限公司 | 一种光学卫星遥感影像检测方法 |
CN115439571A (zh) * | 2022-09-19 | 2022-12-06 | 中科星图空间技术有限公司 | 一种适用于线阵推扫式卫星影像核线影像生成方法及装置 |
CN116188372A (zh) * | 2022-12-14 | 2023-05-30 | 自然资源部黑龙江基础地理信息中心(自然资源部黑龙江测绘资料档案馆) | 一种基于立体测绘的阻水构筑物识别方法 |
-
2024
- 2024-01-31 CN CN202410132864.2A patent/CN117664087B/zh active Active
Patent Citations (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101799290A (zh) * | 2010-02-09 | 2010-08-11 | 武汉大学 | 基于高程基准面的卫星立体影像像对近似核线生成方法以及装置 |
CN102213762A (zh) * | 2011-04-12 | 2011-10-12 | 中交第二公路勘察设计研究院有限公司 | 基于rfm模型的多源星载sar影像自动匹配方法 |
CN103630120A (zh) * | 2013-07-16 | 2014-03-12 | 中国人民解放军信息工程大学 | 基于严密几何模型的火星表面线阵影像核线重采样方法 |
CN104236527A (zh) * | 2014-09-15 | 2014-12-24 | 中国人民解放军总参谋部测绘研究所 | 一种基于投影轨迹法的线阵卫星遥感影像核线特性定量分析方法 |
CN109829939A (zh) * | 2019-01-18 | 2019-05-31 | 南京泛在地理信息产业研究院有限公司 | 一种缩小多视影像匹配同名像点搜索范围的方法 |
CN111667533A (zh) * | 2019-03-08 | 2020-09-15 | 西安中科星图空间数据技术有限公司 | 一种基于有理函数模型获取卫星严格成像模型参数的方法 |
CN110135455A (zh) * | 2019-04-08 | 2019-08-16 | 平安科技(深圳)有限公司 | 影像匹配方法、装置及计算机可读存储介质 |
WO2020206903A1 (zh) * | 2019-04-08 | 2020-10-15 | 平安科技(深圳)有限公司 | 影像匹配方法、装置及计算机可读存储介质 |
CN110853140A (zh) * | 2019-10-11 | 2020-02-28 | 北京空间机电研究所 | 一种dem辅助的光学视频卫星影像稳像方法 |
CN114332085A (zh) * | 2022-03-11 | 2022-04-12 | 西安中科西光航天科技有限公司 | 一种光学卫星遥感影像检测方法 |
CN115439571A (zh) * | 2022-09-19 | 2022-12-06 | 中科星图空间技术有限公司 | 一种适用于线阵推扫式卫星影像核线影像生成方法及装置 |
CN116188372A (zh) * | 2022-12-14 | 2023-05-30 | 自然资源部黑龙江基础地理信息中心(自然资源部黑龙江测绘资料档案馆) | 一种基于立体测绘的阻水构筑物识别方法 |
Non-Patent Citations (3)
Title |
---|
基于投影基准面的线阵推扫式卫星立体像对近似核线影像生成方法;胡芬 等;测绘学报;20091015;第38卷(第05期);第428-436页 * |
基于有理函数模型的线阵推扫式卫星近似核线影像的生成;聂倩 等;城市勘测;20131031(第05期);第60-62页 * |
基于有理函数模型的高分辨率卫星遥感影像匹配;袁修孝 等;武汉大学学报(信息科学版);20090605;第34卷(第06期);第671-674页 * |
Also Published As
Publication number | Publication date |
---|---|
CN117664087A (zh) | 2024-03-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN117664087B (zh) | 垂轨环扫式卫星影像核线生成方法、系统及设备 | |
US9219858B2 (en) | Generating a composite field of view using a plurality of oblique panoramic images of a geographic area | |
CN103310443B (zh) | 高分辨率遥感影像快速处理方法及系统 | |
CN114399588B (zh) | 三维车道线生成方法、装置、电子设备和计算机可读介质 | |
CN112561990B (zh) | 定位信息生成方法、装置、设备和计算机可读介质 | |
CN116182878B (zh) | 道路曲面信息生成方法、装置、设备和计算机可读介质 | |
CN114993328A (zh) | 车辆定位评估方法、装置、设备和计算机可读介质 | |
CN116245757A (zh) | 多模态数据的多场景通用性遥感影像云修复方法和系统 | |
CN117132652A (zh) | 基于三维网格的目标点云标记方法、装置、设备及介质 | |
CN117953063A (zh) | 数据处理方法、扫描方法、装置、设备、存储介质及系统 | |
CN113034679A (zh) | 一种用于输电线路的实时动态三维可视化展示方法 | |
CN113362090B (zh) | 一种用户行为数据处理方法和装置 | |
CN115357675B (zh) | 一种像控点标准化处理建设像控点数据库方法和系统 | |
Baltsavias | Integration of ortho-images in GIS | |
CN111738119A (zh) | 一种基于modis卫星监测的森林山火面积估算方法及装置 | |
CN114742707B (zh) | 多源遥感图像拼接方法、装置、电子设备和可读介质 | |
CN117664088B (zh) | 超宽幅垂轨环扫式卫星影像确定同名点方法、系统及设备 | |
Belmonte et al. | DEM generation from close-range photogrammetry using extended python photogrammetry toolbox | |
Liu et al. | A continuative variable resolution digital elevation model for ground-based photogrammetry | |
Geng et al. | A generic rigorous sensor model for photogrammetric processing of pushbroom planetary images | |
CN114972482B (zh) | 一种遥感影像应用面积统计方法 | |
Gong et al. | Quantitative Assessment of the Projection Trajectory-Based Epipolarity Model and Epipolar Image Resampling for Linear-Array Satellite Images | |
CN117197361B (zh) | 实景三维数据库构建方法、电子设备和计算机可读介质 | |
CN116819469B (zh) | 一种多雷达目标位置同步方法、装置、设备及存储介质 | |
Yang et al. | Three-dimensional lunar surface reconstruction based on remote sensing imagery from Chang E-2 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |