CN112212833A - 机械拼接型tdi ccd推扫相机整体几何平差方法 - Google Patents
机械拼接型tdi ccd推扫相机整体几何平差方法 Download PDFInfo
- Publication number
- CN112212833A CN112212833A CN202010887236.7A CN202010887236A CN112212833A CN 112212833 A CN112212833 A CN 112212833A CN 202010887236 A CN202010887236 A CN 202010887236A CN 112212833 A CN112212833 A CN 112212833A
- Authority
- CN
- China
- Prior art keywords
- ccd
- formula
- slice
- orientation parameter
- inter
- 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
- 238000000034 method Methods 0.000 title claims abstract description 27
- 238000001444 catalytic combustion detection Methods 0.000 claims description 148
- 239000012634 fragment Substances 0.000 claims description 40
- 239000011159 matrix material Substances 0.000 claims description 39
- 238000004364 calculation method Methods 0.000 claims description 33
- 238000003384 imaging method Methods 0.000 claims description 31
- 239000000523 sample Substances 0.000 claims description 14
- 238000012935 Averaging Methods 0.000 claims description 3
- 238000004422 calculation algorithm Methods 0.000 claims description 3
- 238000013507 mapping Methods 0.000 abstract description 2
- 238000010276 construction Methods 0.000 abstract 1
- 230000003287 optical effect Effects 0.000 description 6
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Images
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
-
- 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/02—Picture taking arrangements specially adapted for photogrammetry or photographic surveying, e.g. controlling overlapping of pictures
Landscapes
- Engineering & Computer Science (AREA)
- Multimedia (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Image Processing (AREA)
Abstract
本发明适用于摄影测量与遥感测绘技术领域,具体提供了一种机械拼接型TDI CCD推扫相机整体几何平差方法,该方法提出了构建片间约束整体几何平差模块、片间约束外定向参数平差模型,片间约束内定向参数平差模型,通过计算外定向参数和内定向参数,可有效实现原始分片影像的高效拼接,形成连续一致的拼接影像,并辅以控制点进行整体几何平差,提高拼接影像的定位精度;同时,鉴于外定向参数与内定向参数之间具有很强的相关性,采用分开迭代的方式,确保机械拼接型TDI CCD推扫相机的高拼接精度。
Description
技术领域
本发明属于摄影测量与遥感测绘技术领域,具体涉及一种机械拼接型TDI CCD推扫相机整体几何平差方法。
背景技术
拼接型TDI CCD推扫相机是指按一定方式将多片TDI CCD拼接而形成的组合式并采用推扫成像方式的光学遥感相机。该类相机可兼顾地面分辨率、地面覆盖宽度、成像灵敏度以及信噪比等技术指标,因而在当前高分辨率光学遥感卫星上逐渐得到了广泛应用。与传统线阵推扫式相机相比,拼接型TDI CCD推扫式相机采用了CCD拼接技术。受制作工艺水平限制,单片TDI CCD器件的长度较短,无法满足星载光学相机的地面覆盖宽度要求,因此需采用一定技术手段将多片TDI CCD拼接在一起以使相机具备足够视场宽度。
根据拼接方案的不同,当前拼接型TDI CCD相机主要可分为三类:机械交错拼接相机、光学拼接相机以及特殊构型拼接相机。无论上述哪一类拼接型TDI CCD相机,其均可是空间分辨率、地面覆盖宽度等技术指标显著提升,但其特殊的物理构造和成像机理也使其所获取图像的几何处理呈现出一系列新特点。
机械拼接型TDI CCD推扫相机在某一成像瞬间,多片TDI CCD共享一套外方位元素,但由于焦平面上安置位置的不同,各片TDI CCD在地面上的成像不是一条连续的扫描行,而是多条不连续的“短扫描行”,如图1-1所示,为瞬间成像示意图。地面上“不连续”短扫描行在沿轨方向的错位大小取决于各分片CCD在焦平面上沿轨方向的错位程度。随着卫星平台的飞行,各TDI CCD分片分别推扫成像,形成多个连续的分片窄条带图像,由于相邻TDICCD分片之间的相对安置关系,将导致分片影像之间存在一定水平重叠和垂直错位,如图1-2所示为某各分片影像地面覆盖示意图。故有必要准确确定对机械拼接型TDI CCD推扫相机的内定向参数和外定向参数,确保相机的拼接精度。
作者为曹斌才的2014年硕士学位论文《多片TDI-CCD卫星影像拼接方法研究》提供了一种TDI-CCD卫星影像拼接方法,该方法采用了“先外后内”的在轨几何标定,以得到相机的内定向参数和外定向参数,对相机进行拼接。在轨几何标定的计算过程中,没有利用严格的片间约束几何条件,更没有推导得到片间约束外定向参数平差模型和片间约束内定向参数平差模型,只利用了地面坐标值相等的约束条件循环计算内定向参数和外定向参数,这样得到的内定向参数和外定向参数均不准确,进而无法保证相机的拼接和定位精度。
发明内容
本发明提供了一种机械拼接型TDI CCD推扫相机整体几何平差方法,用以解决现有技术得到的内定向参数和外定向参数均不准确的问题。
为解决上述技术问题,本发明的技术方案包括:
本发明提供了一种机械拼接型TDI CCD推扫相机整体几何平差方法,包括如下步骤:
S1:构建所述机械拼接型TDI CCD推扫相机的片间约束整体几何平差模型;
S2:根据所述机械拼接型TDI CCD推扫相机的影像获取地面控制点信息,并对不同分片CCD进行影像匹配确定相邻分片重叠区域内的同名点信息;
S3:基于所述片间约束整体几何平差模型,根据所述地面控制点信息与所述同名点信息构建片间约束外定向参数平差模型,解算外定向参数的外定向计算值;
S4:固定所述片间约束外定向参数平差模型中的所述外定向计算值,根据所述地面控制点信息与所述同名点信息构建片间约束内定向参数平差模型,解算内定向参数的内定向计算值;
S5:根据所述片间约束外定向参数平差模型与所述片间约束内定向参数平差模型,对所述内定向参数的内定向计算值进行二次平差;
S6:判断所述外定向参数的外定向计算值和二次平差后的内定向参数的内定向计算值是否满足相应的精度要求,若不满足,则执行步骤S3至步骤S5),直至满足精度要求;其中,满足相应的精度要求指的是:外定向参数的外定向计算值小于第一阈值且二次平差后的内定参数的内定向计算值小于第二阈值。
上述技术方案的有益效果为:本发明提出了构建片间约束整体几何平差模块、片间约束外定向参数平差模型,片间约束内定向参数平差模型,通过计算外定向参数和内定向参数,可有效实现原始分片影像的高效拼接,形成连续一致的拼接影像,并辅以控制点进行整体几何平差,提高拼接影像的定位精度;同时,鉴于外定向参数与内定向参数之间具有很强的相关性,采用分开迭代的方式,确保机械拼接型TDI CCD推扫相机的高拼接精度。
进一步的,步骤S1之前还具有以下步骤:
S01:构建所述机械拼接型TDI CCD推扫相机的几何成像模型;
S02:基于所述几何成像模型,构建相邻分片影像的片间约束几何模型。
进一步的,所述几何成像模型为公式1,公式1为:
其中,P(X,Y,Z)为一实际地面控制点的坐标;(XS,YS,ZS)为每一扫描行的实际外定向参数;(x,y)为实际地面控制点所对应的影像点的坐标;λ为投影缩放系数;R为由外定向参数(t,α,κv)确定的旋转矩阵,hij为旋转矩阵R的各矩阵元素,i=1,2,3,j=1,2,3;f为推扫相机参数;为影像点在推扫相机坐标系下探元指向角,且
进一步的,所述片间约束几何模型为公式2和公式3,公式2和公式3分别为:
其中,(XS1,YS1,ZS1)为第一CCD分片的第一外定向参数,λ1为第一CCD分片的第一投影缩放系数,R1为未知外定向参数(Δt,Δα,Δκv)确定的第一旋转矩阵,(x1,y1)为第一CCD分片上第一影像点坐标,为所述第一影像点对应的第一探元指向角;(XS2,YS2,ZS2)为第二CCD分片的第二外定向参数,λ2为第二CCD分片的第二投影缩放系数,R2为未知外定向参数(Δt,Δα,Δκv)确定的第二旋转矩阵,(x2,y2)为第二CCD分片上第二影像点坐标,为所述第二影像点对应的第二探元指向角。
进一步的,所述片间约束整体几何平差模型为公式4,公式4为:
其中,(Δt,Δα,Δκv)为待确定的外定向参数;A1为由未知外定向参数(Δt,Δα,Δκv)确定的第一系数矩阵;(Δx1,Δy1)为第一CCD分片的影像点的第一影像点坐标误差,由所述第一CCD分片的第一内定向参数确定;(Δx2,Δy2)为第二CCD分片的影像点的第二影像点坐标误差,由所述第二CCD分片的第二内定向参数确定;A2为未知外定向参数(Δt,Δα,Δκv)的第二系数矩阵,Δ(-f)为推扫相机参数的误差;(XS1,YS1,ZS1)为第一CCD分片的第一外定向参数,λ1为第一CCD分片的第一投影缩放系数,R1为未知外定向参数(Δt,Δα,Δκv)确定的第一旋转矩阵,(x1,y1)为第一CCD分片上第一影像点坐标,f为推扫相机参数,(XS2,YS2,ZS2)为第二CCD分片的第二外定向参数,λ2为第二CCD分片的第二投影缩放系数,R2为未知外定向参数(Δt,Δα,Δκv)确定的第二旋转矩阵,(x2,y2)为第二CCD分片上第二影像点坐标。
进一步的,所述片间约束外定向参数平差模型为公式5,公式5为:
进一步的,对所述片间约束外定向参数平差模型变形得到公式6,公式6为:
(U1,V1,W1)为第一CCD分片上的第一影像点的第一中间参数,(U2,V2,W2)为第一CCD分片上的第二影像点的第二中间参数,t1为第一CCD分片上扫描行的第一成像时间,t2为第二CCD分片上扫描行的第二成像时间;
将所述第一系数矩阵A1和第二系数矩阵A2代入至公式6,得到平差计算公式7,公式7为:
根据公式7,采用最小二乘平差算法进行计算,得到所述未知外定向参数(Δt,Δα,Δκ)的第一计算值。
进一步的,所述片间约束内定向参数平差模型为公式8,公式8为:
进一步的,采用的指向角模型为公式9,公式9为:
采用的指向角计算公式为公式10,公式10为:
则所述片间约束内定向参数平差模型变形为公式11和公式12,公式11和公式12分别为:
其中,
其中,公式12左边的两项分别为公式13和公式14,公式13和公式14分别为:
将公式13和公式14代入至公式12,进行最小二乘平差就可得到所述第一CCD分片在x方向指向角变化系数c0,c1,c2,c3和在y方向指向角变化系数r0,r1,r2,r3,得到各分片的指向角变化系数,即可计算得到每一CCD分片上每一影像点的影像点坐标。
附图说明
图1-1是现有技术的机械拼接型TDI CCD推扫相机的瞬间成像示意图;
图1-2是现有技术的机械拼接型TDI CCD推扫相机的各分片影像地面覆盖示意图;
图2是本发明的指向角表示视线向量示意图;
图3是本发明的机械拼接型TDI CCD推扫相机整体几何平差方法的流程图。
具体实施方式
本发明的一种机械拼接型TDI CCD推扫相机整体几何平差方法实施例,如图3所示,该机械拼接型TDI CCD推扫相机整体几何平差方法可包括以下步骤:
S1:构建所述机械拼接型TDI CCD推扫相机的片间约束整体几何平差模型;
S2:根据所述机械拼接型TDI CCD推扫相机的影像获取地面控制点信息,并对不同分片CCD进行影像匹配确定相邻分片重叠区域内的同名点信息;
S3:基于所述片间约束整体几何平差模型,根据所述地面控制点信息与所述同名点信息构建片间约束外定向参数平差模型,解算外定向参数的外定向计算值;
S4:固定所述片间约束外定向参数平差模型中的所述外定向计算值,根据所述地面控制点信息与所述同名点信息构建片间约束内定向参数平差模型,解算内定向参数的内定向计算值;
S5:根据所述外片间约束外定向参数平差模型与所述外片间约束内定向参数平差模型,对所述内定向参数所述内定向计算值进行二次平差;
S6:判断所述外定向参数的外定向计算值和二次平差后的内定向参数的内定向计算值是否满足相应的精度要求,若不满足,则执行步骤S3至步骤S5),直至满足精度要求;其中,满足相应的精度要求指的是:外定向参数的外定向计算值小于第一阈值且二次平差后的内定参数的内定向计算值小于第二阈值。
其中,第一阈值、第二阈值是可以根据设备精度、所实施环境、所要求精度进行设置的,并不是一固定单一数值。
在所述步骤S1之前还具有以下步骤:
S01:构建所述机械拼接型TDI CCD推扫相机的几何成像模型;
S02:基于所述几何成像模型,构建相邻分片影像的片间约束几何模型。
其中,对于机械拼接型TDI CCD推扫相机,建立几何成像模型需要的进行对片间约束外定向参数和片间约束内定向参数的确定。
确定外定向参数是确定成像瞬间相机投影中心位置和相机主光轴指向的过程。机械拼接型TDI CCD相机采用推扫方式成像,各扫描行具有不同的外定向参数,可采用轨道姿态模型,用模型参数对各扫描行的外定向参数进行统一描述。常采用一般多项式轨道姿态模型,每一扫描行的外定向参数在星历参数观测初值的基础上随时间线性变化。
其中,(XS,YS,ZS,t,α,κv)为每一扫描行的实际外定向参数,(XS0,YS0,ZS0,t0,α0,κv0)为该扫描行外定向参数的星历参数观测初值,(X′S,Y′S,Z′S,t′,α′,κ′v)为对应扫描行外定向参数的0阶变化量,(X″S,Y″S,Z″S,t″,α″,κ″v)为对应扫描行外定向参数的1阶变化量,t是该扫描行的成像时间,t0是起始成像时刻,t-t0是该扫描行相对于起始成像时刻的时间差。
当考虑相机畸变时,指向角计算公式描述为探元编号S的3阶多项式,如下式所示:
进一步的,作为本发明一优选实施方式,所述几何成像模型为:
其中,P(X,Y,Z)为一实际地面控制点的坐标;(XS,YS,ZS)为每一扫描行的实际外定向参数;(x,y)为实际地面控制点所对应的影像点的坐标;λ为投影缩放系数;R为由外定向参数(t,α,κv)确定的旋转矩阵;hij为旋转矩阵R的矩阵元素,i=1,2,3,j=1,2,3;f为推扫相机参数;为影像点在推扫相机坐标系下探元指向角,且
对于机械拼接型TDI CCD推扫相机,相邻CCD分片重叠区域对同一地面覆盖区进行成像,同一地面点在相邻CCD分片重叠区域形成两个对应像点。设地面实际控制点P的坐标为(X,Y,Z),在第一CCD分片上形成影像点p1,对应像平面坐标为(x1,y1),其对应的第一指向角为在第二CCD分片上形成影像点p2,对应像平面坐标为(x2,y2),其对应的第二指向角为分别对影像点p1和p2列立几何成像模型,如下式所示。
其中(XS1,YS1,ZS1,t1,α1,κv1)是第一CCD分片的第一外定向参数,(XS2,YS2,ZS2,t2,α2,κv2)是第二CCD分片的第二外定向参数。
进一步的,作为本发明一优选实施方式,所述几何约束模型为:
其中,(XS1,YS1,ZS1)为第一CCD分片的第一外定向参数,λ1为第一CCD分片的第一投影缩放系数,R1为未知外定向参数(Δt,Δα,Δκv)确定的第一旋转矩阵,(x1,y1)为第一CCD分片上第一影像点坐标,为第一影像点对应的第一探元指向角;(XS2,YS2,ZS2)为第二CCD分片的第二外定向参数,λ2为第二CCD分片的第二投影缩放系数,R2为未知外定向参数(Δt,Δα,Δκv)确定的第二旋转矩阵,(x2,y2)为第二CCD分片上第二影像点坐标,为第二影像点对应的第二探元指向角。
其中,对影像点p1和p2列立几何成像模型中所列出的两式进行联立,即可得到公式(2)、(3)。
进一步的,作为本发明一优选实施方式,所述片间约束整体几何平差模型为:
其中,(Δt,Δα,Δκv)为待确定的外定向参数;A1是由所述未知外定向参数(Δt,Δα,Δκv)确定的第一系数矩阵;(Δx1,Δy1)是所述第一CCD分片的影像点的第一影像点坐标误差;(Δx2,Δy2)是所述第二CCD分片的影像点的第二影像点坐标误差;A2是由所述未知外定向参数(Δt,Δα,Δκv)确定的第二系数矩阵;Δ(-f)是所述推扫相机参数误差。
其中,鉴于机械拼接型TDI CCD相机相邻CCD分片在焦平面沿轨方向的几何错位误差、推扫成像过程中的外定向参数误差、以及不同CCD分片自身的内定向参数误差都使拼接影像产生拼接错误,并降低定位精度。
为了提高定位精度,将公式(2)、(3)进行线性化,一阶泰勒级数展开后,得到:
其中,(ΔXS1,ΔYS1,ΔZS1,Δt1,Δα1,Δκv1)为第一CCD分片的第一外定向参数误差,(ΔXS2,ΔYS2,ΔZS2,Δt2,Δα2,Δκv2)为第二CCD分片的第二外定向参数误差,A1为第一CCD分片外定向参数(Δt1,Δα1,Δκv1)的第一系数矩阵,A2为第二CCD分片外定向参数(Δt2,Δα2,Δκv2)的第二系数矩阵。
机械拼接型TDI CCD相机各CCD分片共享一套光学系统,采用一组外定向参数模型描述外定向参数变化,因此(ΔXS1,ΔYS1,ΔZS1,Δt1,Δα1,Δκv1)=(ΔXS2,ΔYS2,ΔZS2,Δt2,Δα2,Δκv2)=(ΔXS,ΔYS,ΔZS,Δt,Δα,Δκv)。但是不同CCD分片的内定向参数互不相同,各分片的像点坐标误差是独立的,(Δx1,Δy1)为第一CCD分片上第一影像点p1的像点坐标误差,(Δx2,Δy2)为第二CCD分片上第二影像点p2的像点坐标误差。
对上述公式进行化简,即可得到公式(4):
进一步的,作为本发明一优选实施方式,所述片间约束外定向参数平差模型为:
其中,(Δt,Δα,Δκv)为所述待确定外定向参数。
进一步的,作为本发明一优选实施方式,对所述片间约束外定向参数平差模型变形得到:
(U1,V1,W1)为第一CCD分片上的第一影像点的第一中间参数,(U2,V2,W2)为第一CCD分片上的第二影像点的第二中间参数,t1为第一CCD分片上扫描行的第一成像时间,t2为第二CCD分片上扫描行的第二成像时间。
将所述第一系数矩阵A1和第二系数矩阵A2代入至公式(6),得到平差计算公式如下:
根据公式(7),采用最小二乘平差算法进行计算,得到所述未知外定向参数(Δt,Δα,Δκ)的第一计算值,即解算出外定向参数的外定向计算值。
进一步的,作为本发明一优选实施方式,所述片间约束内定向参数平差模型为:
其中,(Δx1,Δy1)是所述第一CCD分片的影像点的所述第一影像点坐标误差,由所述第一CCD分片的第一内定向参数确定;(Δx2,Δy2)是所述第二CCD分片的影像点的所述第二影像点坐标误差,由所述第二CCD分片的第二内定向参数确定。
其中,在计算得到外定向参数后,在外定向参数保持不变的情况下,构建片间约束内定向参数平差模型如公式(8),进行内定向参数的平差计算。如此,采用分开迭代方式进行内定向参数的计算。
进一步的,作为本发明一优选实施方式,采用指向角模型:
采用指向角计算公式:
则所述片间约束内定向参数平差模型变形为:
其中,
其中,公式(12)左边的两项分别为:
将公式(13)和公式(14)代入至公式(12),进行最小二乘平差就可得到所述第一CCD分片在x方向指向角变化系数c0,c1,c2,c3和在y方向指向角变化系数r0,r1,r2,r3,以及所述第二CCD分片在x方向指向角变化系数c0′,c1′,c2′,c3′和在y方向指向角变化系数r0′,r1′,r2′,r3′,得到各分片的指向角变化系数,即可计算得到每一CCD分片上每一影像点的影像点坐标。
其中,对公式(8)进行变形,移项得:
第一CCD分片与第二CCD分片的投影缩放系数一致,λ1=λ2=λ,则公式(15)可进一步变形为:
进一步,在指向角模型(9)、指向角计算公式(10)的基础上,对公式(16)作进一步变换为:
其中:
公式(12)左边的两项分别为:
将公式(13)和公式(14)代入公式(12)中,进行最小二乘平差就可得到第一CCD分片在x方向指向角变化系数c0,c1,c2,c3和在y方向指向角变化系数r0,r1,r2,r3,以及第二CCD分片在x方向指向角变化系数c0′,c1′,c2′,c3′和在y方向指向角变化系数r0′,r1′,r2′,r3′,得到各分片的指向角变化系数,即可计算得到每一CCD分片上每一影像点的影像点坐标。如此,即完成第一次机械拼接行TDI CCD推扫相机的内定向平差。
如此,本发明通过构建片间约束整体几何平差模块、片间约束外定向参数平差模型,片间约束内定向参数平差模型,通过计算外定向参数和内定向参数,可有效实现原始分片影像的高效拼接,形成连续一致的拼接影像,并辅以控制点进行整体几何平差,提高拼接影像的定位精度;同时,鉴于外定向参数与内定向参数之间具有很强的相关性,采用分开迭代的方式,确保机械拼接型TDI CCD推扫相机的高拼接精度。
以上所述实施例仅用以说明本发明的技术方案,而非对其限制;尽管参照前述实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围,均应包含在本发明的保护范围之内。
Claims (9)
1.一种机械拼接型TDI CCD推扫相机整体几何平差方法,其特征在于,包括如下步骤:
S1:构建所述机械拼接型TDI CCD推扫相机的片间约束整体几何平差模型;
S2:根据所述机械拼接型TDI CCD推扫相机的影像获取地面控制点信息,并对不同分片CCD进行影像匹配确定相邻分片重叠区域内的同名点信息;
S3:基于所述片间约束整体几何平差模型,根据所述地面控制点信息与所述同名点信息构建片间约束外定向参数平差模型,解算外定向参数的外定向计算值;
S4:固定所述片间约束外定向参数平差模型中的所述外定向计算值,根据所述地面控制点信息与所述同名点信息构建片间约束内定向参数平差模型,解算内定向参数的内定向计算值;
S5:根据所述片间约束外定向参数平差模型与所述片间约束内定向参数平差模型,对所述内定向参数的内定向计算值进行二次平差;
S6:判断所述外定向参数的外定向计算值和二次平差后的内定向参数的内定向计算值是否满足相应的精度要求,若不满足,则执行步骤S3至步骤S5),直至满足精度要求;其中,满足相应的精度要求指的是:外定向参数的外定向计算值小于第一阈值且二次平差后的内定参数的内定向计算值小于第二阈值。
2.根据权利要求1所述的机械拼接型TDI CCD推扫相机整体几何平差方法,其特征在于,步骤S1之前还具有以下步骤:
S01:构建所述机械拼接型TDI CCD推扫相机的几何成像模型;
S02:基于所述几何成像模型,构建相邻分片影像的片间约束几何模型。
4.根据权利要求2所述的机械拼接型TDI CCD推扫相机整体几何平差方法,其特征在于,所述片间约束几何模型为公式2和公式3,公式2和公式3分别为:
5.根据权利要求1所述的机械拼接型TDI CCD推扫相机整体几何平差方法,其特征在于,所述片间约束整体几何平差模型为公式4,公式4为:
其中,(Δt,Δα,Δκv)为待确定的外定向参数;A1为由未知外定向参数(Δt,Δα,Δκv)确定的第一系数矩阵;(Δx1,Δy1)为第一CCD分片的影像点的第一影像点坐标误差,由所述第一CCD分片的第一内定向参数确定;(Δx2,Δy2)为第二CCD分片的影像点的第二影像点坐标误差,由所述第二CCD分片的第二内定向参数确定;A2为未知外定向参数(Δt,Δα,Δκv)的第二系数矩阵,Δ(-f)为推扫相机参数的误差;(XS1,YS1,ZS1)为第一CCD分片的第一外定向参数,λ1为第一CCD分片的第一投影缩放系数,R1为未知外定向参数(Δt,Δα,Δκv)确定的第一旋转矩阵,(x1,y1)为第一CCD分片上第一影像点坐标,f为推扫相机参数,(XS2,YS2,ZS2)为第二CCD分片的第二外定向参数,λ2为第二CCD分片的第二投影缩放系数,R2为未知外定向参数(Δt,Δα,Δκv)确定的第二旋转矩阵,(x2,y2)为第二CCD分片上第二影像点坐标。
7.根据权利要求6所述的机械拼接型TDI CCD推扫相机整体几何平差方法,其特征在于,对所述片间约束外定向参数平差模型变形得到公式6,公式6为:
(U1,V1,W1)为第一CCD分片上的第一影像点的第一中间参数,(U2,V2,W2)为第一CCD分片上的第二影像点的第二中间参数,t1为第一CCD分片上扫描行的第一成像时间,t2为第二CCD分片上扫描行的第二成像时间;
将所述第一系数矩阵A1和第二系数矩阵A2代入至公式6,得到平差计算公式7,公式7为:
根据公式7,采用最小二乘平差算法进行计算,得到所述未知外定向参数(Δt,Δα,Δκ)的第一计算值。
9.根据权利要求8所述机械拼接型TDI CCD推扫相机整体几何平差方法,其特征在于,采用的指向角模型为公式9,公式9为:
采用的指向角计算公式为公式10,公式10为:
则所述片间约束内定向参数平差模型变形为公式11和公式12,公式11和公式12分别为:
其中,
其中,公式12左边的两项分别为公式13和公式14,公式13和公式14分别为:
将公式13和公式14代入至公式12,进行最小二乘平差就可得到所述第一CCD分片在x方向指向角变化系数c0,c1,c2,c3和在y方向指向角变化系数r0,r1,r2,r3,得到各分片的指向角变化系数,即可计算得到每一CCD分片上每一影像点的影像点坐标。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010887236.7A CN112212833B (zh) | 2020-08-28 | 2020-08-28 | 机械拼接型tdi ccd推扫相机整体几何平差方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010887236.7A CN112212833B (zh) | 2020-08-28 | 2020-08-28 | 机械拼接型tdi ccd推扫相机整体几何平差方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112212833A true CN112212833A (zh) | 2021-01-12 |
CN112212833B CN112212833B (zh) | 2021-07-09 |
Family
ID=74059125
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010887236.7A Expired - Fee Related CN112212833B (zh) | 2020-08-28 | 2020-08-28 | 机械拼接型tdi ccd推扫相机整体几何平差方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112212833B (zh) |
Citations (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6735348B2 (en) * | 2001-05-01 | 2004-05-11 | Space Imaging, Llc | Apparatuses and methods for mapping image coordinates to ground coordinates |
CN101750619A (zh) * | 2010-01-18 | 2010-06-23 | 武汉大学 | 自检校pos直接对地目标定位方法 |
CN101827223A (zh) * | 2010-04-20 | 2010-09-08 | 武汉大学 | 基于行频归一化的非共线tdi ccd成像数据内视场拼接方法 |
CN102778224A (zh) * | 2012-08-08 | 2012-11-14 | 北京大学 | 一种基于极坐标参数化的航空摄影测量光束法平差的方法 |
CN103400345A (zh) * | 2013-07-18 | 2013-11-20 | 西南交通大学 | 一种降低tdi-ccd相机图像模糊度的方法 |
CN105091906A (zh) * | 2015-06-30 | 2015-11-25 | 武汉大学 | 高分辨率光学推扫卫星稳态重成像传感器校正方法及系统 |
CN105551053A (zh) * | 2015-12-01 | 2016-05-04 | 中国科学院上海技术物理研究所 | 一种小面阵星载tdi ccd相机的快速几何精校正方法 |
CN106403902A (zh) * | 2016-08-31 | 2017-02-15 | 武汉大学 | 一种星地协同的光学卫星在轨实时几何定位方法及系统 |
DE102010021819B4 (de) * | 2010-05-28 | 2018-10-18 | Deutsches Zentrum für Luft- und Raumfahrt e.V. | Detektor-Modul und höchstauflösender Satellit |
CN109191532A (zh) * | 2018-08-09 | 2019-01-11 | 王涛 | 一种机载三线阵ccd相机标定方法 |
CN109696155A (zh) * | 2017-10-20 | 2019-04-30 | 武汉大学 | 光线共面约束的弱交会光学卫星影像联合平差方法及系统 |
CN111044037A (zh) * | 2019-12-26 | 2020-04-21 | 中国人民解放军战略支援部队信息工程大学 | 一种光学卫星影像的几何定位方法及装置 |
CN111473802A (zh) * | 2020-05-07 | 2020-07-31 | 武汉大学 | 一种基于线阵推扫的光学传感器内方位元素定标方法 |
CN111508029A (zh) * | 2020-04-09 | 2020-08-07 | 武汉大学 | 星载分片线阵ccd光学相机整体几何定标方法及系统 |
-
2020
- 2020-08-28 CN CN202010887236.7A patent/CN112212833B/zh not_active Expired - Fee Related
Patent Citations (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6735348B2 (en) * | 2001-05-01 | 2004-05-11 | Space Imaging, Llc | Apparatuses and methods for mapping image coordinates to ground coordinates |
CN101750619A (zh) * | 2010-01-18 | 2010-06-23 | 武汉大学 | 自检校pos直接对地目标定位方法 |
CN101827223A (zh) * | 2010-04-20 | 2010-09-08 | 武汉大学 | 基于行频归一化的非共线tdi ccd成像数据内视场拼接方法 |
DE102010021819B4 (de) * | 2010-05-28 | 2018-10-18 | Deutsches Zentrum für Luft- und Raumfahrt e.V. | Detektor-Modul und höchstauflösender Satellit |
CN102778224A (zh) * | 2012-08-08 | 2012-11-14 | 北京大学 | 一种基于极坐标参数化的航空摄影测量光束法平差的方法 |
CN103400345A (zh) * | 2013-07-18 | 2013-11-20 | 西南交通大学 | 一种降低tdi-ccd相机图像模糊度的方法 |
CN105091906A (zh) * | 2015-06-30 | 2015-11-25 | 武汉大学 | 高分辨率光学推扫卫星稳态重成像传感器校正方法及系统 |
CN105551053A (zh) * | 2015-12-01 | 2016-05-04 | 中国科学院上海技术物理研究所 | 一种小面阵星载tdi ccd相机的快速几何精校正方法 |
CN106403902A (zh) * | 2016-08-31 | 2017-02-15 | 武汉大学 | 一种星地协同的光学卫星在轨实时几何定位方法及系统 |
CN109696155A (zh) * | 2017-10-20 | 2019-04-30 | 武汉大学 | 光线共面约束的弱交会光学卫星影像联合平差方法及系统 |
CN109191532A (zh) * | 2018-08-09 | 2019-01-11 | 王涛 | 一种机载三线阵ccd相机标定方法 |
CN111044037A (zh) * | 2019-12-26 | 2020-04-21 | 中国人民解放军战略支援部队信息工程大学 | 一种光学卫星影像的几何定位方法及装置 |
CN111508029A (zh) * | 2020-04-09 | 2020-08-07 | 武汉大学 | 星载分片线阵ccd光学相机整体几何定标方法及系统 |
CN111473802A (zh) * | 2020-05-07 | 2020-07-31 | 武汉大学 | 一种基于线阵推扫的光学传感器内方位元素定标方法 |
Non-Patent Citations (11)
Title |
---|
LIANG-LIANG ZHENG ET AL.: "An improved adaptive preprocessing method for TDI CCD images", 《OPTOELECTRONICS LETTERS》 * |
LIHONG YANG ET AL.: "Precision Detection of CCD Splicing Based on Template matching Algorithm", 《INTERNATIONAL CONFERENCE ON COMPUTING CONTROL AND INDUSTRIAL ENGINEERING》 * |
MAOTENG ZHENG ET AL.: "Self-Calibration Adjustment of CBERS-02B Long-Strip Imagery", 《IEEE TRANSACTIONS ON GEOSCIENCE AND REMOTE SENSING》 * |
YONGJUN ZHANG ET AL.: "BUNDLE BLOCK ADJUSTMENT WITH SELF-CALIBRATION OF LONG ORBIT CBERS-02B IMAGERY", 《REMOTE SENSING AND SPATIAL INFORMATION SCIENCES》 * |
孟伟灿 等: "TDI CCD交错拼接推扫相机严格几何模型构建与优化", 《测绘学报》 * |
常琳 等: "航天TDI CCD相机成像拼接快速配准算法设计与分析", 《光学学报》 * |
张星祥 等: "TDI CCD焦平面的机械交错拼接", 《光学学报》 * |
潘红播 等: "资源三号测绘卫星传感器校正产品几何模型", 《测绘学报》 * |
王涛 等: "国产机载大视场三线阵CCD相机几何标定", 《遥感学报》 * |
王涛: "线阵CCD传感器实验场几何定标的理论与方法研究", 《中国博士学位论文全文数据库 基础科学辑》 * |
胡芬: "三片非共线TDI CCD成像数据内视场拼接理论与算法研究", 《中国博士学位论文全文数据库 信息科技辑》 * |
Also Published As
Publication number | Publication date |
---|---|
CN112212833B (zh) | 2021-07-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107014399B (zh) | 一种星载光学相机-激光测距仪组合系统联合检校方法 | |
US6735348B2 (en) | Apparatuses and methods for mapping image coordinates to ground coordinates | |
CN103557841B (zh) | 一种提高多相机合成影像摄影测量精度的方法 | |
Hu et al. | Understanding the rational function model: methods and applications | |
CN110006452B (zh) | 高分六号宽视场相机相对几何定标方法及系统 | |
CN107144293A (zh) | 一种视频卫星面阵相机的几何定标方法 | |
CN106871787A (zh) | 大空间线扫描成像三维测量方法 | |
CN110111274A (zh) | 一种星载推扫式光学传感器外方位元素定标方法 | |
CN105698766A (zh) | 考虑定向参数精度信息的卫星影像rfm模型区域网平差方法 | |
CN111524196B (zh) | 一种摆扫大幅宽光学卫星在轨几何标定方法 | |
CN113313769A (zh) | 一种光学卫星多面阵传感器片间无缝几何定标方法 | |
Daakir et al. | Improvement of photogrammetric accuracy by modeling and correcting the thermal effect on camera calibration | |
Wang et al. | Geometric calibration for the aerial line scanning camera GFXJ | |
CN114514412A (zh) | 误差估计装置、误差估计方法、误差估计程序 | |
CN114241064B (zh) | 一种遥感卫星内外方位元素实时几何定标方法 | |
CN115311366A (zh) | 基于rpc模型的星载分片线阵传感器几何定标方法及系统 | |
CN111156969A (zh) | 一种宽幅遥感影像立体测绘方法及系统 | |
Pi et al. | On-orbit geometric calibration using a cross-image pair for the linear sensor aboard the agile optical satellite | |
CN105628052B (zh) | 基于控制直线的光学卫星传感器在轨几何定标方法及系统 | |
CN113065277B (zh) | 协同多载荷数据的高分辨率遥感卫星颤振探测与建模方法 | |
CN112802118B (zh) | 一种光学卫星传感器在轨分时几何定标方法 | |
JP5619866B2 (ja) | 対称な露光写真を利用した地球観測システムに係るアライメント誤差のキャリブレーション方法 | |
CN112212833B (zh) | 机械拼接型tdi ccd推扫相机整体几何平差方法 | |
CN111275773A (zh) | 无场几何定标方法及系统 | |
CN109696155B (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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20210709 |
|
CF01 | Termination of patent right due to non-payment of annual fee |