CN108593111A - 时空调制干涉成像光谱仪的运动成像仿真方法及装置 - Google Patents
时空调制干涉成像光谱仪的运动成像仿真方法及装置 Download PDFInfo
- Publication number
- CN108593111A CN108593111A CN201810550722.2A CN201810550722A CN108593111A CN 108593111 A CN108593111 A CN 108593111A CN 201810550722 A CN201810550722 A CN 201810550722A CN 108593111 A CN108593111 A CN 108593111A
- Authority
- CN
- China
- Prior art keywords
- imaging
- interference
- image
- space
- moment
- 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
- 238000003384 imaging method Methods 0.000 title claims abstract description 141
- 238000004088 simulation Methods 0.000 title claims abstract description 53
- 238000000034 method Methods 0.000 title claims abstract description 34
- 238000001228 spectrum Methods 0.000 claims abstract description 53
- 238000012952 Resampling Methods 0.000 claims abstract description 44
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 23
- 238000004611 spectroscopical analysis Methods 0.000 claims abstract description 8
- 238000004364 calculation method Methods 0.000 claims description 20
- 230000003595 spectral effect Effects 0.000 claims description 16
- VMXUWOKSQNHOCA-UKTHLTGXSA-N ranitidine Chemical compound [O-][N+](=O)\C=C(/NC)NCCSCC1=CC=C(CN(C)C)O1 VMXUWOKSQNHOCA-UKTHLTGXSA-N 0.000 claims description 8
- 238000006073 displacement reaction Methods 0.000 claims description 7
- 238000005315 distribution function Methods 0.000 claims description 7
- 238000010008 shearing Methods 0.000 claims description 7
- 238000006467 substitution reaction Methods 0.000 claims 1
- 230000003068 static effect Effects 0.000 description 10
- 230000008569 process Effects 0.000 description 6
- 238000005516 engineering process Methods 0.000 description 5
- 230000002441 reversible effect Effects 0.000 description 5
- 230000008901 benefit Effects 0.000 description 3
- 230000008859 change Effects 0.000 description 3
- 230000007547 defect Effects 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 230000003287 optical effect Effects 0.000 description 3
- HUTDUHSNJYTCAR-UHFFFAOYSA-N ancymidol Chemical compound C1=CC(OC)=CC=C1C(O)(C=1C=NC=NC=1)C1CC1 HUTDUHSNJYTCAR-UHFFFAOYSA-N 0.000 description 1
- 238000005314 correlation function Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000002427 irreversible effect Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000013178 mathematical model Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 238000010408 sweeping Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01J—MEASUREMENT OF INTENSITY, VELOCITY, SPECTRAL CONTENT, POLARISATION, PHASE OR PULSE CHARACTERISTICS OF INFRARED, VISIBLE OR ULTRAVIOLET LIGHT; COLORIMETRY; RADIATION PYROMETRY
- G01J3/00—Spectrometry; Spectrophotometry; Monochromators; Measuring colours
- G01J3/28—Investigating the spectrum
- G01J3/45—Interferometric spectrometry
Landscapes
- Physics & Mathematics (AREA)
- Spectroscopy & Molecular Physics (AREA)
- General Physics & Mathematics (AREA)
- Spectrometry And Color Measurement (AREA)
Abstract
本发明公开了一种时空调制干涉成像光谱仪的运动成像仿真方法及装置，其中，方法包括：将光谱数据立方体的位置作为物面光谱网格位置，并基于共线条件方程获取物面光谱网格位置在瞬时成像时刻所对应的成像面上的变形的像面光谱网格位置；基于傅里叶光谱学理论由光谱数据立方体的光谱数据获取瞬时像面位置对应的干涉强度信息；利用基于Sutherland‑Hodgman算法的面积配比算法对像面光谱网格位置和干涉强度信息进行重采样，以得到探测器标准网格位置处的干涉强度值，进而获取单帧推扫干涉图像；改变瞬时成像时刻的姿态参数和位置参数，以获取推扫干涉图像序列。该方法仿真度高、可广泛应用。
Description
技术领域
本发明涉及仪器数据仿真技术领域，特别涉及一种时空调制干涉成像光谱仪的运动成像仿真方法及装置。
背景技术
仪器仿真技术是通过建立数学模型模拟真实仪器的各种工作模式和工作状态，是人们在真实仪器使用前就能较准确地了解该仪器在一定的工作环境和工作条件下可能产生的结果和存在的问题，从而衡量真实仪器的性能，发现仪器潜在的设计缺陷，进一步指导仪器的设计，完善仪器的功能，该技术在航空航天遥感仪器的设计、研制及应用领域的作用尤为重大。
时空调制干涉型成像光谱仪作为一种新型成像光谱仪，克服了时间调制型干涉成像光谱仪稳定性差及空间调制型干涉成像光谱仪低通量的不足，具有多通道、更高的光通量、更好的光谱和空间分辨率等优点，已成为遥感领域的一大热点，正日益显示出重要的实用价值。
时空调制干涉成像光谱仪利用推扫方式扫描目标景物，在不同时刻目标成像在探测器的不同位置，遍历所有时刻得到目标点的完整点干涉图，再由点干涉图通过傅里叶变换复原产生该目标的光谱。因此，时空调制成像光谱仪遥感数据获取与搭载平台运动具有很强的相关性，平台的非平稳运动将导致成像位置的无规则变化。由于时空调制干涉成像光谱仪获得地物点的光谱曲线必须经过一次全视场的推扫过程，平台非平稳运动导致光轴指向位置与理想位置发生偏差，像面上的每个像素点的光谱信息是其周围理想位置的光谱的混合信息，因此，需要通过适当的方法仿真出合适的运动推扫图像，否则复原出的光谱就无法反映景物目标的真实属性。
发明内容
本申请是基于发明人对以下问题的认识和发现作出的：
相关技术中，(1)针对空间调制型干涉超光谱成像仪开展了仿真技术研究，该模型基于辐射能量传输，不涉及运动成像引入的位置偏差。(2)针对星载的干涉型成像光谱仪提出了一种基于归一化互相关系数的光谱数据立方体实现数据仿真的方法，由于星载平台运动较平稳，该方法按照理想运动成像仿真考虑，未考虑运动成像引入的成像位置变化。(3)数据仿真方法首先在光谱维度上即单谱段图像上进行图像重采样，对成像位置作光谱混合，然后基于傅里叶光谱学理论将数据转化到干涉维度上，得到单帧推扫图像。虽然考虑了位置偏差，但是在光谱反演时，需要在单帧推扫图像上进行图像重采样，两个过程不可逆，因此，这种仿真方法会引入不必要的误差项。
本发明旨在至少在一定程度上解决相关技术中的技术问题之一。
为此，本发明的一个目的在于提出一种时空调制干涉成像光谱仪的运动成像仿真方法，该方法不仅可以应用于时空调制干涉成像光谱仪的静态仿真和运动成像仿真，还可以应用于空间调制干涉成像光谱仪的静态仿真和运动成像仿真，应用广泛，仿真度高。
本发明的另一个目的在于提出一种时空调制干涉成像光谱仪的运动成像仿真装置。
为达到上述目的，本发明一方面实施例提出了一种时空调制干涉成像光谱仪的运动成像仿真方法，包括以下步骤：将光谱数据立方体的位置作为物面光谱网格位置，并基于共线条件方程获取所述物面光谱网格位置在瞬时成像时刻所对应的成像面上的变形的像面光谱网格位置；基于傅里叶光谱学理论由所述光谱数据立方体的光谱数据获取瞬时像面位置对应的干涉强度信息；利用基于Sutherland-Hodgman算法的面积配比算法对所述像面光谱网格位置和所述干涉强度信息进行重采样，以得到探测器标准网格位置处的干涉强度值，进而获取单帧推扫干涉图像；改变所述瞬时成像时刻的姿态参数和位置参数，以获取推扫干涉图像序列。
本发明实施例的时空调制干涉成像光谱仪的运动成像仿真方法，考虑了运动成像因素的影响，并将光谱维度重采样修改为在干涉强度维度重采样，与后续的数据反演过程可逆，避免了引入不可控误差，从而能更好衡量真实仪器的性能，发现仪器潜在的设计缺陷，仿真度高，从而不仅可以应用于时空调制干涉成像光谱仪的静态仿真和运动成像仿真，还可以应用于空间调制干涉成像光谱仪的静态仿真和运动成像仿真，应用广泛。
另外，根据本发明上述实施例的时空调制干涉成像光谱仪的运动成像仿真方法还可以具有以下附加的技术特征：
进一步地，在本发明的一个实施例中，所述干涉强度信息的计算公式为：
其中，B(v)为该像元对应地物的光谱分布函数，v为波数，v1、v2为波数范围，d为干涉仪的剪切量。
进一步地，在本发明的一个实施例中，所述获取单帧推扫干涉图像，进一步包括：根据中心位置按最近邻原则在变形网格中划出可能相交的区域；通过所述Sutherland-Hodgman算法获取所述可能相交的区域与变形网格的四边形的相交的图形的所有顶点位置，并根据所述所有顶点位置得到图形面积；根据所述图形面积和面积配比对干涉图像进行重采样，以得到每个探测器对应的重采样干涉强度值，获取所述单帧干涉推扫图像。
进一步地，在本发明的一个实施例中，其中，
所述面积的计算公式为：
其中，(x1,y1)(x2,y2)…(xn,yn)为相交面积的顶点坐标，；
重采样的计算公式为：
其中，Gi为标准网格正方形与变形网格四边形相交的图形面积所占标准网格的比例，Ii是每个变形网格对应的干涉强度值。
进一步地，在本发明的一个实施例中，所述获取推扫干涉图像序列，进一步包括：通过瞬时摄影测量中心与外方位角元素在不同成像时刻发生变化代入每一时刻的瞬时值，以得到所述每一时刻的干涉推扫图像；当代入预设连续时刻的位置和姿态数据时，获取所述推扫干涉图像序列。
为达到上述目的，本发明另一方面实施例提出了一种时空调制干涉成像光谱仪的运动成像仿装置，包括：第一获取模块，用于将光谱数据立方体的位置作为物面光谱网格位置，并基于共线条件方程获取所述物面光谱网格位置在瞬时成像时刻所对应的成像面上的变形的像面光谱网格位置；第二获取模块，用于基于傅里叶光谱学理论由所述光谱数据立方体的光谱数据获取瞬时像面位置对应的干涉强度信息；重采样模块，用于利用基于Sutherland-Hodgman算法的面积配比算法对所述像面光谱网格位置和所述干涉强度信息进行重采样，以得到探测器标准网格位置处的干涉强度值，进而获取单帧推扫干涉图像；推扫干涉模块，用于改变所述瞬时成像时刻的姿态参数和位置参数，以获取推扫干涉图像序列。
本发明实施例的时空调制干涉成像光谱仪的运动成像仿真装置，考虑了运动成像因素的影响，并将光谱维度重采样修改为在干涉强度维度重采样，与后续的数据反演过程可逆，避免了引入不可控误差，从而能更好衡量真实仪器的性能，发现仪器潜在的设计缺陷，仿真度高，从而不仅可以应用于时空调制干涉成像光谱仪的静态仿真和运动成像仿真，还可以应用于空间调制干涉成像光谱仪的静态仿真和运动成像仿真，应用广泛。
另外，根据本发明上述实施例的时空调制干涉成像光谱仪的运动成像仿装置还可以具有以下附加的技术特征：
进一步地，在本发明的一个实施例中，所述干涉强度信息的计算公式为：
其中，B(v)为该像元对应地物的光谱分布函数，v为波数，v1、v2为波数范围，d为干涉仪的剪切量。
进一步地，在本发明的一个实施例中，所述重采样模块进一步用于根据中心位置按最近邻原则在变形网格中划出可能相交的区域，通过所述Sutherland-Hodgman算法获取所述可能相交的区域与变形网格的四边形的相交的图形的所有顶点位置，并根据所述所有顶点位置得到图形面积，并根据所述图形面积和面积配比对干涉图像进行重采样，以得到每个探测器对应的重采样干涉强度值，获取所述单帧干涉推扫图像。
进一步地，在本发明的一个实施例中，其中，
所述面积的计算公式为：
其中，(x1,y1)(x2,y2)…(xn,yn)为相交面积的顶点坐标，；
重采样的计算公式为：
其中，Gi为标准网格正方形与变形网格四边形相交的图形面积所占标准网格的比例，Ii是每个变形网格对应的干涉强度值。
进一步地，在本发明的一个实施例中，所述推扫干涉模块进一步用于通过瞬时摄影测量中心与外方位角元素在不同成像时刻发生变化代入每一时刻的瞬时值，以得到所述每一时刻的干涉推扫图像，并当代入预设连续时刻的位置和姿态数据时，获取所述推扫干涉图像序列。
本发明附加的方面和优点将在下面的描述中部分给出，部分将从下面的描述中变得明显，或通过本发明的实践了解到。
附图说明
本发明上述的和/或附加的方面和优点从下面结合附图对实施例的描述中将变得明显和容易理解，其中：
图1为根据本发明一个实施例的时空调制干涉成像光谱仪的运动成像仿真方法的流程图；
图2为根据本发明一个具体实施例的时空调制干涉成像光谱仪的运动成像仿真方法的流程图；
图3为根据本发明一个实施例的共线条件方程示意图；
图4为根据本发明一个实施例的面积配比示意图；
图5为根据本发明一个实施例的时空调制干涉成像光谱仪的运动成像仿真装置的结构示意图。
具体实施方式
下面详细描述本发明的实施例，所述实施例的示例在附图中示出，其中自始至终相同或类似的标号表示相同或类似的元件或具有相同或类似功能的元件。下面通过参考附图描述的实施例是示例性的，旨在用于解释本发明，而不能理解为对本发明的限制。
下面参照附图描述根据本发明实施例提出的时空调制干涉成像光谱仪的运动成像仿真方法及装置，首先将参照附图描述根据本发明实施例提出的时空调制干涉成像光谱仪的运动成像仿真方法。
图1是本发明一个实施例的时空调制干涉成像光谱仪的运动成像仿真方法的流程图。
如图1所示，该时空调制干涉成像光谱仪的运动成像仿真方法包括以下步骤：
在步骤S101中，将光谱数据立方体的位置作为物面光谱网格位置，并基于共线条件方程获取物面光谱网格位置在瞬时成像时刻所对应的成像面上的变形的像面光谱网格位置。
可以理解的是，如图2所示，首先，本发明实施例计算物面光谱网格位置在瞬时成像时刻所对应的成像面上的变形的像面光谱网格位置。
具体而言，将已知的光谱数据立方体的位置作为物面光谱网格位置，基于共线条件方程，如图3所示，计算物面光谱网格位置在瞬时成像时刻所对应的成像面上的变形的像面光谱网格位置。
假设物面光谱网格中心位置为
基于共线条件方程计算物面网格中任一点A在瞬时像面上的位置：
其中，(XA,YA,ZA)为物面地物A在物空间坐标系中的坐标，(XS,YS,ZS)为瞬时摄影测量中心S在物空间坐标系中坐标位置，f为成像镜焦距，a1～a3、b1～b3、c1～c3为方向余弦，为瞬时成像时刻的外方位角元素，(x,y)为瞬时像面上对应点A1的位置坐标。
因此，成像瞬时时刻的像面上的位置为：
在步骤S102中，基于傅里叶光谱学理论由光谱数据立方体的光谱数据获取瞬时像面位置对应的干涉强度信息。
可以理解的是，如图2所示，计算对应的瞬时像面位置对应的干涉强度信息。
进一步地，在本发明的一个实施例中，干涉强度信息的计算公式为：
其中，B(v)为该像元对应地物的光谱分布函数，v为波数，v1、v2为波数范围，d为干涉仪的剪切量。
具体而言，基于傅里叶光谱学理论，由已知的光谱数据立方体的光谱数据计算对应的瞬时像面位置对应的干涉强度信息。利用上述瞬时成像位置数据，可以得到瞬时像面上每个像点与零光程差之间的距离。与零光程差位置相距x处的干涉强度可以通过下面公式计算得到：
其中，B(v)为该像元对应地物的光谱分布函数，v为波数，v1、v2为波数范围，d为干涉仪的剪切量。对瞬时像面上每个像点计算干涉强度，得到变形的像面干涉数据。
在步骤S103中，利用基于Sutherland-Hodgman算法的面积配比算法对像面光谱网格位置和干涉强度信息进行重采样，以得到探测器标准网格位置处的干涉强度值，进而获取单帧推扫干涉图像。
可以理解的是，对瞬时成像的像面变形网格位置和对应的干涉强度值进行重采样，得到单帧推扫干涉图像。利用基于Sutherland-Hodgman算法的面积配比算法对瞬时成像的像面变形网格位置和对应的干涉强度值进行重采样，得到探测器标准网格位置处的干涉强度值，即单帧推扫干涉图像。
进一步地，在本发明的一个实施例中，获取单帧推扫干涉图像，进一步包括：根据中心位置按最近邻原则在变形网格中划出可能相交的区域；通过Sutherland-Hodgman算法获取可能相交的区域与变形网格的四边形的相交的图形的所有顶点位置，并根据所有顶点位置得到图形面积；根据图形面积和面积配比对干涉图像进行重采样，以得到每个探测器对应的重采样干涉强度值，获取单帧干涉推扫图像。
在本发明的一个实施例中，其中，面积的计算公式为：
其中，(x1,y1)(x2,y2)…(xn,yn)为相交面积的顶点坐标，；
重采样的计算公式为：
其中，Gi为标准网格正方形与变形网格四边形相交的图形面积所占标准网格的比例，Ii是每个变形网格对应的干涉强度值。
具体而言，时空调制干涉成像光谱仪使用面阵CCD作为数据接收传感器，CCD探测器的像元排列在标准规则的网格位置，因此需要对上述得到变形的像面的干涉数据进行重采样。对任一标准网格正方形，网格中心位置坐标为(p,q)。首先根据其中心位置按最近邻原则在变形的网格中划出可能相交的区域G1～G9，如图4所示，然后利用Sutherland-Hodgman算法计算其与变形网格四边形G1的相交的图形的所有顶点位置，最后根据顶点位置计算面积，面积计算公式如下：
其中，(x1,y1)(x2,y2)…(xn,yn)为相交面积的顶点坐标，按逆时针顺序排列。同理计算该标准网格正方形与变形网格四边形G2～G9的相交的图形面积。按照面积配比对干涉图像进行重采样，重采样计算公式为：
对探测器面上的每一个标准网格正方形按照上述操作，计算出对应的重采样干涉强度值，得到单帧的干涉推扫图像。
在步骤S104中，改变瞬时成像时刻的姿态参数和位置参数，以获取推扫干涉图像序列。
可以理解的是，本发明实施例改变瞬时成像时刻的姿态参数和位置参数，重复步骤S101-S103，获取推扫干涉图像序列
进一步地，在本发明的一个实施例中，获取推扫干涉图像序列，进一步包括：通过瞬时摄影测量中心与外方位角元素在不同成像时刻发生变化代入每一时刻的瞬时值，以得到每一时刻的干涉推扫图像；当代入预设连续时刻的位置和姿态数据时，获取推扫干涉图像序列。
具体而言，瞬时摄影测量中心S(XS,YS,ZS)与外方位角元素在不同成像时刻会发生变化，代入每一时刻的瞬时值，重复步骤S101-S103，得到每一时刻的干涉推扫图像。当代入一段连续时刻的位置和姿态数据时，获取推扫干涉图像序列。
根据本发明实施例提出的时空调制干涉成像光谱仪的运动成像仿装置，考虑了运动成像因素的影响，并将光谱维度重采样修改为在干涉强度维度重采样，与后续的数据反演过程可逆，避免了引入不可控误差，从而能更好衡量真实仪器的性能，发现仪器潜在的设计缺陷，仿真度高，从而不仅可以应用于时空调制干涉成像光谱仪的静态仿真和运动成像仿真，还可以应用于空间调制干涉成像光谱仪的静态仿真和运动成像仿真，应用广泛。
其次参照附图描述根据本发明实施例提出的时空调制干涉成像光谱仪的运动成像仿装置。
图5是本发明一个实施例的时空调制干涉成像光谱仪的运动成像仿装置的结构示意图。
如图5所示，该时空调制干涉成像光谱仪的运动成像仿装置10包括：第一获取模块100、第二获取模块200、重采样模块300和推扫干涉模块400。
其中，第一获取模块100用于将光谱数据立方体的位置作为物面光谱网格位置，并基于共线条件方程获取物面光谱网格位置在瞬时成像时刻所对应的成像面上的变形的像面光谱网格位置。第二获取模块200用于基于傅里叶光谱学理论由光谱数据立方体的光谱数据获取瞬时像面位置对应的干涉强度信息。重采样模块300用于利用基于Sutherland-Hodgman算法的面积配比算法对像面光谱网格位置和干涉强度信息进行重采样，以得到探测器标准网格位置处的干涉强度值，进而获取单帧推扫干涉图像。推扫干涉模块400用于改变瞬时成像时刻的姿态参数和位置参数，以获取推扫干涉图像序列。本发明实施例的装置10充分考虑了时空调制干涉成像光谱仪的光谱反演过程，提出的运动成像仿真方法与其完全可逆，避免了在光谱维度上的光谱混合引入的误差项，仿真度高、可广泛应用。
进一步地，在本发明的一个实施例中，干涉强度信息的计算公式为：
其中，B(v)为该像元对应地物的光谱分布函数，v为波数，v1、v2为波数范围，d为干涉仪的剪切量。
进一步地，在本发明的一个实施例中，重采样模块进一步用于根据中心位置按最近邻原则在变形网格中划出可能相交的区域，通过Sutherland-Hodgman算法获取可能相交的区域与变形网格的四边形的相交的图形的所有顶点位置，并根据所有顶点位置得到图形面积，并根据图形面积和面积配比对干涉图像进行重采样，以得到每个探测器对应的重采样干涉强度值，获取单帧干涉推扫图像。
进一步地，在本发明的一个实施例中，其中，面积的计算公式为：
其中，(x1,y1)(x2,y2)…(xn,yn)为相交面积的顶点坐标，；
重采样的计算公式为：
其中，Gi为标准网格正方形与变形网格四边形相交的图形面积所占标准网格的比例，Ii是每个变形网格对应的干涉强度值。
进一步地，在本发明的一个实施例中，推扫干涉模块进一步用于通过瞬时摄影测量中心与外方位角元素在不同成像时刻发生变化代入每一时刻的瞬时值，以得到每一时刻的干涉推扫图像，并当代入预设连续时刻的位置和姿态数据时，获取推扫干涉图像序列。
需要说明的是，前述对时空调制干涉成像光谱仪的运动成像仿真方法实施例的解释说明也适用于该实施例的时空调制干涉成像光谱仪的运动成像仿真装置，此处不再赘述。
根据本发明实施例提出的时空调制干涉成像光谱仪的运动成像仿真装置，考虑了运动成像因素的影响，并将光谱维度重采样修改为在干涉强度维度重采样，与后续的数据反演过程可逆，避免了引入不可控误差，从而能更好衡量真实仪器的性能，发现仪器潜在的设计缺陷，仿真度高，从而不仅可以应用于时空调制干涉成像光谱仪的静态仿真和运动成像仿真，还可以应用于空间调制干涉成像光谱仪的静态仿真和运动成像仿真，应用广泛。
此外，术语“第一”、“第二”仅用于描述目的，而不能理解为指示或暗示相对重要性或者隐含指明所指示的技术特征的数量。由此，限定有“第一”、“第二”的特征可以明示或者隐含地包括至少一个该特征。在本发明的描述中，“多个”的含义是至少两个，例如两个，三个等，除非另有明确具体的限定。
在本说明书的描述中，参考术语“一个实施例”、“一些实施例”、“示例”、“具体示例”、或“一些示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本发明的至少一个实施例或示例中。在本说明书中，对上述术语的示意性表述不必须针对的是相同的实施例或示例。而且，描述的具体特征、结构、材料或者特点可以在任一个或多个实施例或示例中以合适的方式结合。此外，在不相互矛盾的情况下，本领域的技术人员可以将本说明书中描述的不同实施例或示例以及不同实施例或示例的特征进行结合和组合。
尽管上面已经示出和描述了本发明的实施例，可以理解的是，上述实施例是示例性的，不能理解为对本发明的限制，本领域的普通技术人员在本发明的范围内可以对上述实施例进行变化、修改、替换和变型。
Claims (10)
1.一种时空调制干涉成像光谱仪的运动成像仿真方法，其特征在于，包括以下步骤：
将光谱数据立方体的位置作为物面光谱网格位置，并基于共线条件方程获取所述物面光谱网格位置在瞬时成像时刻所对应的成像面上的变形的像面光谱网格位置；
基于傅里叶光谱学理论由所述光谱数据立方体的光谱数据获取瞬时像面位置对应的干涉强度信息；
利用基于Sutherland-Hodgman算法的面积配比算法对所述像面光谱网格位置和所述干涉强度信息进行重采样，以得到探测器标准网格位置处的干涉强度值，进而获取单帧推扫干涉图像；以及
改变所述瞬时成像时刻的姿态参数和位置参数，以获取推扫干涉图像序列。
2.根据权利要求1所述的时空调制干涉成像光谱仪的运动成像仿真方法，其特征在于，所述干涉强度信息的计算公式为：
其中，B(v)为该像元对应地物的光谱分布函数，v为波数，v1、v2为波数范围，d为干涉仪的剪切量。
3.根据权利要求1所述的时空调制干涉成像光谱仪的运动成像仿真方法，其特征在于，所述获取单帧推扫干涉图像，进一步包括：
根据中心位置按最近邻原则在变形网格中划出可能相交的区域；
通过所述Sutherland-Hodgman算法获取所述可能相交的区域与变形网格的四边形的相交的图形的所有顶点位置，并根据所述所有顶点位置得到图形面积；
根据所述图形面积和面积配比对干涉图像进行重采样，以得到每个探测器对应的重采样干涉强度值，获取所述单帧干涉推扫图像。
4.根据权利要求3所述的时空调制干涉成像光谱仪的运动成像仿真方法，其特征在于，其中，
所述面积的计算公式为：
其中，(x1,y1)(x2,y2)…(xn,yn)为相交面积的顶点坐标，；
重采样的计算公式为：
其中，Gi为标准网格正方形与变形网格四边形相交的图形面积所占标准网格的比例，Ii是每个变形网格对应的干涉强度值。
5.根据权利要求1-4任一项所述的时空调制干涉成像光谱仪的运动成像仿真方法，其特征在于，所述获取推扫干涉图像序列，进一步包括：
通过瞬时摄影测量中心与外方位角元素在不同成像时刻发生变化代入每一时刻的瞬时值，以得到所述每一时刻的干涉推扫图像；
当代入预设连续时刻的位置和姿态数据时，获取所述推扫干涉图像序列。
6.一种时空调制干涉成像光谱仪的运动成像仿装置，其特征在于，包括：
第一获取模块，用于将光谱数据立方体的位置作为物面光谱网格位置，并基于共线条件方程获取所述物面光谱网格位置在瞬时成像时刻所对应的成像面上的变形的像面光谱网格位置；
第二获取模块，用于基于傅里叶光谱学理论由所述光谱数据立方体的光谱数据获取瞬时像面位置对应的干涉强度信息；
重采样模块，用于利用基于Sutherland-Hodgman算法的面积配比算法对所述像面光谱网格位置和所述干涉强度信息进行重采样，以得到探测器标准网格位置处的干涉强度值，进而获取单帧推扫干涉图像；以及
推扫干涉模块，用于改变所述瞬时成像时刻的姿态参数和位置参数，以获取推扫干涉图像序列。
7.根据权利要求6所述的时空调制干涉成像光谱仪的运动成像仿真装置，其特征在于，所述干涉强度信息的计算公式为：
其中，B(v)为该像元对应地物的光谱分布函数，v为波数，v1、v2为波数范围，d为干涉仪的剪切量。
8.根据权利要求6所述的时空调制干涉成像光谱仪的运动成像仿真装置，其特征在于，所述重采样模块进一步用于根据中心位置按最近邻原则在变形网格中划出可能相交的区域，通过所述Sutherland-Hodgman算法获取所述可能相交的区域与变形网格的四边形的相交的图形的所有顶点位置，并根据所述所有顶点位置得到图形面积，并根据所述图形面积和面积配比对干涉图像进行重采样，以得到每个探测器对应的重采样干涉强度值，获取所述单帧干涉推扫图像。
9.根据权利要求8所述的时空调制干涉成像光谱仪的运动成像仿真装置，其特征在于，其中，
所述面积的计算公式为：
其中，(x1,y1)(x2,y2)…(xn,yn)为相交面积的顶点坐标，；
重采样的计算公式为：
其中，Gi为标准网格正方形与变形网格四边形相交的图形面积所占标准网格的比例，Ii是每个变形网格对应的干涉强度值。
10.根据权利要求6-9任一项所述的时空调制干涉成像光谱仪的运动成像仿真装置，其特征在于，所述推扫干涉模块进一步用于通过瞬时摄影测量中心与外方位角元素在不同成像时刻发生变化代入每一时刻的瞬时值，以得到所述每一时刻的干涉推扫图像，并当代入预设连续时刻的位置和姿态数据时，获取所述推扫干涉图像序列。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810550722.2A CN108593111B (zh) | 2018-05-31 | 2018-05-31 | 时空调制干涉成像光谱仪的运动成像仿真方法及装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810550722.2A CN108593111B (zh) | 2018-05-31 | 2018-05-31 | 时空调制干涉成像光谱仪的运动成像仿真方法及装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108593111A true CN108593111A (zh) | 2018-09-28 |
| CN108593111B CN108593111B (zh) | 2019-09-24 |
Family
ID=63630171
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810550722.2A Active CN108593111B (zh) | 2018-05-31 | 2018-05-31 | 时空调制干涉成像光谱仪的运动成像仿真方法及装置 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108593111B (zh) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109785339A (zh) * | 2018-12-29 | 2019-05-21 | 西安交通大学 | 一种同时分割和配准精密复杂曲面干涉条纹图像的方法 |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2012184962A (ja) * | 2011-03-03 | 2012-09-27 | Kagawa Univ | 分光特性測定装置及び分光特性測定方法 |
| CN102944308A (zh) * | 2012-10-17 | 2013-02-27 | 北京航空航天大学 | 一种时空联合调制干涉成像光谱仪姿态误差校正方法 |
| CN103728022B (zh) * | 2014-01-16 | 2015-11-04 | 北京航空航天大学 | 一种坏像元的校正方法 |
| CN206132219U (zh) * | 2016-09-29 | 2017-04-26 | 浙江大学 | 一种基于时空联合调制和mzi光波导阵列的成像光谱仪 |
| CN106908148A (zh) * | 2017-03-06 | 2017-06-30 | 中国科学院合肥物质科学研究院 | 一种红外干涉图重建方法 |
| CN107271040A (zh) * | 2017-05-18 | 2017-10-20 | 上海卫星工程研究所 | 时间调制型傅里叶变换光谱仪的双向干涉图协同使用方法 |
| CN107990982A (zh) * | 2017-09-25 | 2018-05-04 | 上海卫星工程研究所 | 傅里叶变换光谱仪光谱计算中的相位校正方法 |
-
2018
- 2018-05-31 CN CN201810550722.2A patent/CN108593111B/zh active Active
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2012184962A (ja) * | 2011-03-03 | 2012-09-27 | Kagawa Univ | 分光特性測定装置及び分光特性測定方法 |
| CN102944308A (zh) * | 2012-10-17 | 2013-02-27 | 北京航空航天大学 | 一种时空联合调制干涉成像光谱仪姿态误差校正方法 |
| CN103728022B (zh) * | 2014-01-16 | 2015-11-04 | 北京航空航天大学 | 一种坏像元的校正方法 |
| CN206132219U (zh) * | 2016-09-29 | 2017-04-26 | 浙江大学 | 一种基于时空联合调制和mzi光波导阵列的成像光谱仪 |
| CN106908148A (zh) * | 2017-03-06 | 2017-06-30 | 中国科学院合肥物质科学研究院 | 一种红外干涉图重建方法 |
| CN107271040A (zh) * | 2017-05-18 | 2017-10-20 | 上海卫星工程研究所 | 时间调制型傅里叶变换光谱仪的双向干涉图协同使用方法 |
| CN107990982A (zh) * | 2017-09-25 | 2018-05-04 | 上海卫星工程研究所 | 傅里叶变换光谱仪光谱计算中的相位校正方法 |
Non-Patent Citations (1)
| Title |
|---|
| 王倩: "卫星线阵影像处理关键技术研究", 《中国优秀硕士学位论文全文数据库信息科技辑》 * |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109785339A (zh) * | 2018-12-29 | 2019-05-21 | 西安交通大学 | 一种同时分割和配准精密复杂曲面干涉条纹图像的方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN108593111B (zh) | 2019-09-24 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Daneshmand et al. | 3d scanning: A comprehensive survey | |
| CN100489446C (zh) | 基于相位法的三维轮廓测量方法 | |
| CN104050716B (zh) | 一种海上多目标sar图像可视化建模方法 | |
| Haala et al. | Quality of 3D point clouds from highly overlapping UAV imagery | |
| CN104392435B (zh) | 鱼眼相机标定方法及标定装置 | |
| CN105510913B (zh) | 基于类光学像方改正的异源光学和sar遥感影像联合定位方法 | |
| Cavegn et al. | Benchmarking high density image matching for oblique airborne imagery | |
| CN107014312A (zh) | 一种振镜式线激光扫描三维测量系统的整体标定方法 | |
| CN107621636A (zh) | 一种基于psi的大型铁路桥梁健康监测方法 | |
| CN104658012A (zh) | 一种基于惯性与光学测量融合的运动捕捉方法 | |
| CN102927908A (zh) | 机器人手眼系统结构光平面参数标定装置及方法 | |
| Tournadre et al. | UAV linear photogrammetry | |
| CN101975947A (zh) | 二维镜像综合孔径辐射成像方法 | |
| CN110728671A (zh) | 基于视觉的无纹理场景的稠密重建方法 | |
| CN110191335A (zh) | 遥感相机成像质量评价与仿真系统 | |
| CN108955571A (zh) | 双频外差与相移编码相结合的三维测量方法 | |
| CN103616682B (zh) | 一种基于曲面投影的多基线InSAR处理方法 | |
| CN113390514B (zh) | 一种基于多传感器阵列的三维红外温度测量方法 | |
| CN104867179A (zh) | 一种全谱段光学成像仪遥感影像仿真方法 | |
| CN107102165A (zh) | 一种基于粒子图像测速的表面流场测量方法 | |
| CN108593111B (zh) | 时空调制干涉成像光谱仪的运动成像仿真方法及装置 | |
| Bybee et al. | Method for 3-D scene reconstruction using fused LiDAR and imagery from a texel camera | |
| CN105953921A (zh) | 气溶胶参数差异条件下对地观测辐射图像的快速仿真方法 | |
| CN107036555B (zh) | 一种交叉光轴光栅投影测量仿真系统及其实现方法 | |
| CN106650049A (zh) | 一种静轨面阵遥感相机分时动态成像模拟方法 |
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 |