CN114187345A - 一种基于星地相机成像系统的视轴指向解算方法 - Google Patents
一种基于星地相机成像系统的视轴指向解算方法 Download PDFInfo
- Publication number
- CN114187345A CN114187345A CN202010854464.4A CN202010854464A CN114187345A CN 114187345 A CN114187345 A CN 114187345A CN 202010854464 A CN202010854464 A CN 202010854464A CN 114187345 A CN114187345 A CN 114187345A
- Authority
- CN
- China
- Prior art keywords
- ground
- camera
- star
- imaging module
- 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.)
- Granted
Links
- 238000003384 imaging method Methods 0.000 title claims abstract description 113
- 230000000007 visual effect Effects 0.000 title claims abstract description 71
- 238000000034 method Methods 0.000 title claims abstract description 51
- 239000011159 matrix material Substances 0.000 claims abstract description 31
- 238000005070 sampling Methods 0.000 claims abstract description 18
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 12
- 238000005516 engineering process Methods 0.000 claims abstract description 9
- 238000013519 translation Methods 0.000 claims description 21
- 230000003287 optical effect Effects 0.000 claims description 10
- 238000006243 chemical reaction Methods 0.000 claims description 7
- HCUOEKSZWPGJIM-YBRHCDHNSA-N (e,2e)-2-hydroxyimino-6-methoxy-4-methyl-5-nitrohex-3-enamide Chemical compound COCC([N+]([O-])=O)\C(C)=C\C(=N/O)\C(N)=O HCUOEKSZWPGJIM-YBRHCDHNSA-N 0.000 claims description 5
- 101100162200 Aspergillus parasiticus (strain ATCC 56775 / NRRL 5862 / SRRC 143 / SU-1) aflD gene Proteins 0.000 claims description 5
- 238000004364 calculation method Methods 0.000 claims description 4
- 238000000605 extraction Methods 0.000 claims description 4
- 238000001514 detection method Methods 0.000 abstract description 9
- 238000013507 mapping Methods 0.000 abstract description 5
- 238000004088 simulation Methods 0.000 description 10
- 238000010586 diagram Methods 0.000 description 8
- 230000008859 change Effects 0.000 description 5
- 230000008569 process Effects 0.000 description 4
- 238000013461 design Methods 0.000 description 3
- 238000009434 installation Methods 0.000 description 2
- 230000009471 action Effects 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 230000005484 gravity Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000000737 periodic effect Effects 0.000 description 1
- 230000009897 systematic effect Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/60—Analysis of geometric attributes
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/30—Determination of transform parameters for the alignment of images, i.e. image registration
- G06T7/33—Determination of transform parameters for the alignment of images, i.e. image registration using feature-based methods
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10032—Satellite or aerial image; Remote sensing
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20048—Transform domain processing
- G06T2207/20056—Discrete and fast Fourier transform, [DFT, FFT]
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Computer Vision & Pattern Recognition (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Geometry (AREA)
- Image Processing (AREA)
Abstract
本发明提供一种基于星地相机成像系统的视轴指向解算方法,首先结合共光路一体化设计星地相机对天成像模块和对地成像模块,且对天成像模块与对地成像模块之间的夹角经过地面标定,然后通过互相关算法实现模板图像和检测图像之间的粗配准,再通过上采样矩阵、卷积核以及傅里叶变换等一系列操作实现模板图像和检测图像之间的精配准,得到主相机视轴相对于星地相机对地成像模块视轴的偏移量,最后基于偏移量得到主相机的视轴指向;由此可见,本发明通过互相关算法的亚像元细分技术与上采样矩阵、傅里叶变换等高精度图像配准技术,可以计算出主相机的高精度的视轴指向,进而提高主相机的测绘精度,具有较强的易用性和实用性。
Description
技术领域
本发明属于遥感卫星成像技术领域,尤其涉及一种基于星地相机成像系统的视轴指向解算方法。
背景技术
光学遥感相机受到重力变化和微小温度场变化,都会引起视轴漂移;为保证成像质量,相机在轨实时调焦,也会引起视轴指向变化,故需实时确定相机的视轴指向。星敏感器指向与相机视轴之间的夹角存在系统周期性变化,这种变化如不加以预先消除或控制,会降低遥感图像的地面定位精度,制约遥感相机的应用性能。现有方法涉及星敏感器、星场与相机的安装关系,监测视轴夹角的变化,但受到星敏的精度限制,及长时间安装矩阵的系统误差导致于其计算精度不高。
发明内容
为解决上述问题,本发明提供一种基于星地相机成像系统的视轴指向解算方法,可以计算出主相机的高精度的视轴指向,进而提高主相机的测绘精度,具有较强的易用性和实用性。
一种基于星地相机成像系统的视轴指向解算方法,所述成像系统涉及星敏感器、主相机以及星地相机,所述星地相机包括共光路一体化设计、且视轴夹角固定的星地相机对天成像模块和星地相机对地成像模块,所述方法具体包括以下步骤:
S1:在轨获取星敏感器的指向,然后根据星敏感器的指向调节星地相机对天成像模块的指向,使得两者的指向差值小于设定值;
S2:星地相机对天成像模块获取星空图像Img1;
S3:根据星空图像Img1提取星点质心及星点质心的赤经、赤纬,然后根据星点质心及星点质心的赤经、赤纬获取星地相机对天成像模块的视轴矢量v天;
S4:根据地面标定的星地相机对天成像模块与星地相机对地成像模块的视轴转换矩阵C1及视轴矢量v天,获取星地相机对地成像模块的指向角α地和β地;
S5:星地相机对地成像模块获取地面某一特征景物图像f(x,y),同时,主相机在同一时刻也对包含同一特征景物的场景进行拍摄,得到图像g(x,y);
S6:将图像f(x,y)和图像g(x,y)进行配准,得到主相机相对于星地相机对地成像模块的指向角偏移量α1和β1;
S7:根据星地相机对地成像模块的指向角与指向角偏移量获取主相机的视轴指向:
α主=α地+α1
β主=β地+β1
其中,α主为弧矢方向上的视轴指向,β主为子午方向上的视轴指向。
进一步地,步骤S6中所述指向角偏移量α1和β1的获取方法包括以下步骤:
S61:采用互相关算法获取图像f(x,y)与图像g(x,y)之间的整数级平移量(row_shift,col_shift),其中,row_shift为行方向上的平移量,col_shift为列方向上的平移量;
S62:计算上采样矩阵in,其中,上采样矩阵的行数和列数分别为nr和nc:
in=buf2ft.*conj(buf1ft)
其中,buf2ft为图像f(x,y)的傅里叶变换,buf1ft为图像g(x,y)的傅里叶变换,conj为matlab中的复共轭函数,.*表示两个矩阵按位置相乘;
S63:计算行方向卷积核kernr和列方向卷积核kernc:
kernr=exp((-i×2π/(nr×usfac))×([0:nor-1].'-roff)×(ifftshift([0:nr-1])-floor(nr/2)))
kernc=exp((-i×2π/(nc×usfac))×(tmpnc-floor(nc/2))×([0:noc-1]-coff))
其中,i为虚部单位,roff=dftshift-row_shift×usfac,coff=dftshift-col_shift×usfac,dftshift为中间变量,usfac为设定的上采样因子,且dftshift=fix(ceil(usfac×1.5)/2),fix和ceil为matlab中的取整函数,floor为matlab中的向下取整函数,iffshift为matlab中的傅里叶反变换函数,nor=noc=usfac×1.5,exp为matlab中的底数为e的指数函数,tmpnc为ifftshift([0:nc-1])的转置,.'表示转置;
S64:计算上采样矩阵in的快速傅里叶变换out:
out=kernr×in×kernc
S65:计算上采样相位相关结果CC:
CC=conj(out)/(md2×nd2×usfac2)
其中,md2=m/2,nd2=n/2,m和n分别为模板图像f(x,y)的行数和列数,conj为matlab中的复共轭函数;
S66:获取图像f(x,y)与图像g(x,y)之间的最终平移量:
row_shift^=row_shift+rloc/usfac
col_shift^=col_shift+cloc/usfac
其中,row_shift^为行方向上的最终平移量,col_shift^为列方向上的最终平移量,rloc与cloc分别为上采样相位相关结果CC中最大值所在的行号和列号;
S67:根据最终平移量计算主相机相对于星地相机对地成像模块的指向角偏移量α1和β1:
α1=θ×row_shift^
β1=θ×col_shift^
其中,θ为主相机的单位像元的指向角。
进一步地,步骤S67中单位像元的指向角θ的计算方法为:
θ=arctan(a/f)
其中,f为主相机的焦距,a为主相机的单位像元大小。
进一步地,步骤S3中所述的视轴矢量v天的获取方法包括以下步骤:
S31:基于内插细分技术,采用加权平方的质心提取方法从星空图像Img1提取N个亚像素级的星点质心位置;
S32:获取N个星点质心位置在导航星表中的赤经、赤纬坐标(αi,δi),i=1,2,...,N,并将赤经、赤纬坐标转换为单位向量坐标(cosαicosδi,sinαicosδi,sinδi)
S33:采用QUEST算法,基于N个星点质心位置的单位向量坐标解算出星地相机相对于天体坐标系J2000的姿态矩阵C;
S34:根据姿态矩阵获取星地相机对天成像模块在天体坐标系J2000下的视轴矢量v天=C-1×v1,其中,C-1为姿态矩阵C的逆矩阵,v1为星地相机对天成像模块视轴的单位矢量。
进一步地,步骤S4中所述星地相机对地成像模块的指向角α地和β地的获取方法为:
S41:根据地面标定的星地相机对天成像模块与星地相机对地成像模块的视轴转换矩阵C1,得到星地相机对地成像模块视轴在天体坐标系J2000下的单位矢量v地=C1×v天=(x地,y地,z地);
S42:将单位矢量v地=C1×v天=(x地,y地,z地)转换为星地相机对地成像模块在球体坐标系下的指向角:
α地=arctan(y地/x地)
有益效果:
本发明提供一种基于星地相机成像系统的视轴指向解算方法,首先基于星敏感器引导星地相机对天成像、计算星地相机对天成像视轴指向、结合共光路一体化设计,对天成像模块与对地成像模块之间的夹角经过地面标定且保持不变,计算星地相机对地成像视轴指向;通过星地相机对地成像,然后将对地成像模块得到的图像和主相机得到的图像进行配准,得到主相机视轴相对于星地相机对地成像模块视轴的偏移量,最后基于偏移量得到主相机的视轴指向;由此可见,本发明通过星敏感器引导星地相机对天成像、星地相机对天成像视轴解算、星地相机对地成像视轴解算、星地相机图像与主相机图像配准,解算出主相机的高精度的视轴指向,进而提高主相机的测绘精度,具有较强的易用性和实用性。
附图说明
图1为本发明提供的一种基于星地相机成像系统的视轴指向解算方法的流程图;
图2为本发明提供的一种基于星地相机成像系统的视轴指向解算中高精度配准方法的原理框图;
图3为本发明提供的主相机与星地相机配准示意图;
图4为本发明提供的星点成像模拟中加权质心法-误差随着噪声均值影响的曲线示意图;
图5为本发明提供的主相机采集的图像;
图6为本发明提供的星地相机采集的图像;
图7为本发明提供的仿真结果示意图;
图8为本发明提供的另一个仿真结果示意图。
具体实施方式
为了使本技术领域的人员更好地理解本申请方案,下面将结合本申请实施例中的附图,对本申请实施例中的技术方案进行清楚、完整地描述。
本发明为针对于星地相机成像光学设计系统开发的视轴解算方法,星地相机采用离轴TMA光学系统,其特点是相机的一部分视场对天,一部分视场对地,而且两部分共焦面,进而可以用于高精度的标定主相机的视轴指向。在光学遥感相机上设置星地成像相机,星地相机具有同时对星空和地面成像能力,且对地成像视场与主相机视场重合,同时,星地相机两个视场通过统一光路系统和成像焦面,具有固定的角度保持能力,在此基础上,采用高精度图像配准方法能够标定主相机对地指向精度。
如图1所示,为一种基于星地相机成像系统的视轴指向解算方法的流程图,星地相机包括共光路一体化设计、且视轴夹角固定的星地相机对天成像模块和星地相机对地成像模块。如图2所示,为一种基于星地相机成像系统的视轴指向解算方法的原理框图,该方法流程主要包括:星敏感器引导星地相机对天成像、星地相机对天成像视轴解算、星地相机对地成像视轴解算、星地相机图像与主相机图像配准、主相机视轴指向解算,具体包括以下步骤:
S1:在轨获取星敏感器的指向,然后根据星敏感器的指向调节星地相机对天成像模块的指向,使得两者的指向差值小于设定值;
S2:星地相机对天成像模块获取星空图像Img1,其中,星空图像Img1包括N个星点;
S3:根据星空图像Img1获取星地相机对天成像模块的视轴矢量v天,具体包括以下步骤:
S31:基于内插细分技术,采用加权平方的质心提取方法从星空图像Img1提取N个亚像素级的星点质心位置(xi,yi);
S32:获取N个星点质心位置在导航星表中的赤经、赤纬坐标(αi,δi),i=1,2,...,N,并将赤经、赤纬坐标转换为单位向量坐标(cosαicosδi,sinαicosδi,sinδi)
S33:采用QUEST算法,基于N个星点质心位置的单位向量坐标解算出星地相机相对于天体坐标系J2000的姿态矩阵C;
S34:根据姿态矩阵获取星地相机对天成像模块在天体坐标系J2000下的视轴矢量v天=C-1×v1,其中,C-1为姿态矩阵C的逆矩阵,v1为星地相机对天成像模块视轴的单位矢量;
S4:根据地面标定的星地相机对天成像模块与星地相机对地成像模块的视轴转换矩阵C1及视轴矢量v天,获取星地相机对地成像模块的指向角α地和β地,具体包括以下步骤:
S41:根据地面标定的星地相机对天成像模块与星地相机对地成像模块的视轴转换矩阵C1,得到星地相机对地成像模块视轴在天体坐标系J2000下的单位矢量v地=C1×v天=(x地,y地,z地),v1=(0,0,1)';
S42:将单位矢量v地=C1×v天=(x地,y地,z地)转换为星地相机对地成像模块在球体坐标系下的指向角:
α地=arctan(y地/x地)
需要说明的是,在星地相机对天成像模块对星成像前,在轨获取星敏感器的指向,据此调整星地相机的指向,然后再对星成像,是为了进一步提高成像精度。
S5:星地相机对地成像获取地面某一特征景物图像作为图像f(x,y),同时,主相机在同一时刻也对同一特征景物进行拍摄,得到检测图像g(x,y);
如图3所示,实线区域代表不同时刻主相机探测器线阵成像示意图,虚线区域代表不同成像时刻星地相机面阵探测器成像示意图。当然,主相机探测器不仅限于推扫线阵成像,星地相机也不局限于面阵探测器,此处只是为了举例说明。。
由此可见,星地相机成像系统对地成像模块与主相机存在视场重叠部分,也就是说,同一个时刻可以获取到对一特征景物的图像,则通过后续步骤中的高精度配准方法,可以完成星地相机图像与主相机图像的配准。
S6:将图像f(x,y)和图像g(x,y)进行配准,得到主相机相对于星地相机对地成像模块的指向角偏移量α1和β1,具体包括以下步骤:
S61:采用互相关算法获取模板图像f(x,y)与检测图像g(x,y)之间的整数级平移量(row_shift,col_shift),其中,row_shift为行方向上的平移量,col_shift为列方向上的平移量;
需要说明的是,步骤S61为初始定位,得到的是“粗定位”结果,鉴于互相关算法为图像配准领域中的常规方法,本发明对此不赘述;下面详细阐述后续的“细定位”过程。
S62:计算上采样矩阵in,其中,上采样矩阵的行数和列数分别为nr和nc:
in=buf2ft.*conj(buf1ft)
其中,buf2ft为模板图像f(x,y)的傅里叶变换,buf1ft为检测图像g(x,y)的傅里叶变换,conj为matlab中的复共轭函数,.*表示两个矩阵按位置相乘;
需要说明的是,上采样网格区域的大小的计算方式如下:
row=round(row_shift×usfac)/usfac
col=round(col_shift×usfac)/usfac
其中,row为上采样网格区域的行大小,col为上采样网格区域的列大小,usfac为设定的上采样因子,可取usfac=100,round为matlab中的四舍五入函数;
S63:计算行方向卷积核kernr和列方向卷积核kernc:
kernr=exp((-i×2π/(nr×usfac))×([0:nor-1].'-roff)×(ifftshift([0:nr-1])-floor(nr/2)))
kernc=exp((-i×2π/(nc×usfac))×(tmpnc-floor(nc/2))×([0:noc-1]-coff))
其中,i为虚部单位,roff=dftshift-row_shift×usfac,coff=dftshift-col_shift×usfac,dftshift为中间变量,且dftshift=fix(ceil(usfac×1.5)/2),fix和ceil为matlab中的取整函数,floor为matlab中的向下取整函数,iffshift为matlab中的傅里叶反变换函数,nor=noc=usfac×1.5,exp为matlab中的底数为e的指数函数,tmpnc为ifftshift([0:nc-1])的转置,.'表示转置;[0:nor-1]表示0,1,2,...nor-1,[0:nr-1]表示0,1,2,...nr-1,[0:noc-1]表示0,1,2,...noc-1;
S64:计算上采样矩阵in的快速傅里叶变换out:
out=kernr×in×kernc
S65:计算上采样相位相关结果CC:
CC=conj(out)/(md2×nd2×usfac2)
其中,md2=m/2,nd2=n/2,m和n分别为模板图像f(x,y)的行数和列数,conj为matlab中的复共轭函数;其中,检测图像g(x,y)与模板图像f(x,y)的行数和列数相同;
S66:获取模板图像f(x,y)与检测图像g(x,y)之间的最终平移量:
row_shift^=row_shift+rloc/usfac
col_shift^=col_shift+cloc/usfac
其中,row_shift^为行方向上的最终平移量,col_shift^为列方向上的最终平移量,rloc与cloc分别为上采样相位相关结果CC中最大值所在的行号和列号;
S67:根据最终平移量计算主相机相对于星地相机对地成像模块的指向角偏移量α1和β1:
α1=θ×row_shift^
β1=θ×col_shift^
其中,θ为主相机的单位像元的指向角,且θ=arctan(a/f),f为主相机的焦距,a为主相机的单位像元大小;
S7:根据星地相机对地成像模块的指向角与指向角偏移量获取主相机的视轴指向:
α主=α地+α1
β主=β地+β1
其中,α主为弧矢方向上的视轴指向,β主为子午方向上的视轴指向。
由此可见,本发明结合共光路一体化设计星地相机对天和对地成像,通过亚像元细分技术、高精度图像配准技术,计算出主相机的高精度的视轴指向进而可以提高主相机的测绘精度,具有较强的易用性和实用性。
需要说明的是,虽然本发明提供了一种高精度图像配准方法,星地相机对天指向求解方法采用的QUEST方法及质心提取方法,但是对于相关应用,只要保证精度,配法可以根据需求选择。同样地,对于各环节采用的方法,对于本领域普通技术人员来讲,很显然可以采用其他方法代替。但是,本发明提出的主相机视轴指向解算的主体思路和流程是不变的。
下面结合具体实施例对本发明提供的解算方法进行实验仿真:
首先进行星点成像模拟,星的能量分布可以近似为二维高斯分布,在实验仿真过程,使用以下星点灰度模型:
在CCD面阵上模拟像素的中心坐标值在(3,3)的位置加上[0,1]偏移量,图像的系统误差主要表现在灰度值上受噪声的均值影响。因此在图像中加入均方差为0.02的高斯噪声,均值μ范围设置为0~20,每种噪声进行1000次随机实验,求加权平方亚像元星点定位的均方根值RMS,结果如下图4所示。加权权质心算法的精度,在灰度系统误差均值在17以内(仿真过程中图像采用8bit),平方加权灰度质心法的精度保持在0.1像素以内,且在灰度系统误差在10时,精度小于0.05像素,工程上保守估计0.1个像素,采用设计的星地相机光学系统,角分辨率为0.66”,则星地相机指向角精度为0.066”。
对于图像配准,采用多组进行仿真,图5和图6为其中一组,在仿真过程,将一幅图像加入不同程度的高斯白噪声模糊,以模拟两个传感器由于自身设计或者成像工况对配准误差的影响。加入均值从0到30,方差值0.5的噪声。附图7和图8为对应的仿真结果,可见图像配准误差保守估计为0.1个像元,采用设计的星地相机光学系统,角分辨率为0.66”,则配准精度为0.066”。
综上所述,本发明提供一种基于星地相机成像系统的视轴指向解算方法,首先结合共光路一体化设计星地相机对天成像模块和对地成像模块,且对天成像模块与对地成像模块之间的夹角经过地面标定,然后通过互相关算法实现模板图像和检测图像之间的粗配准,再通过上采样矩阵、卷积核以及傅里叶变换等一系列操作实现模板图像和检测图像之间的精配准,得到主相机视轴相对于星地相机对地成像模块视轴的偏移量,最后基于偏移量得到主相机的视轴指向;由此可见,本发明通过互相关算法的亚像元细分技术与上采样矩阵、傅里叶变换等高精度图像配准技术,可以计算出主相机的高精度的视轴指向,进而提高主相机的测绘精度,具有较强的易用性和实用性。
当然,本发明还可有其他多种实施例,在不背离本发明精神及其实质的情况下,熟悉本领域的技术人员当然可根据本发明作出各种相应的改变和变形,但这些相应的改变和变形都应属于本发明所附的权利要求的保护范围。
Claims (5)
1.一种基于星地相机成像系统的视轴指向解算方法,所述成像系统涉及星敏感器、主相机以及星地相机,其特征在于,所述星地相机包括共光路一体化设计、且视轴夹角固定的星地相机对天成像模块和星地相机对地成像模块,所述方法具体包括以下步骤:
S1:在轨获取星敏感器的指向,然后根据星敏感器的指向调节星地相机对天成像模块的指向,使得两者的指向差值小于设定值;
S2:星地相机对天成像模块获取星空图像Img1;
S3:根据星空图像Img1提取星点质心及星点质心的赤经、赤纬,然后根据星点质心及星点质心的赤经、赤纬获取星地相机对天成像模块的视轴矢量v天;
S4:根据地面标定的星地相机对天成像模块与星地相机对地成像模块的视轴转换矩阵C1及视轴矢量v天,获取星地相机对地成像模块的指向角α地和β地;
S5:星地相机对地成像模块获取地面某一特征景物图像f(x,y),同时,主相机在同一时刻也对包含同一特征景物的场景进行拍摄,得到图像g(x,y);
S6:将图像f(x,y)和图像g(x,y)进行配准,得到主相机相对于星地相机对地成像模块的指向角偏移量α1和β1;
S7:根据星地相机对地成像模块的指向角与指向角偏移量获取主相机的视轴指向:
α主=α地+α1
β主=β地+β1
其中,α主为弧矢方向上的视轴指向,β主为子午方向上的视轴指向。
2.如权利要求1所述的一种基于星地相机成像系统的视轴指向解算方法,其特征在于,步骤S6中所述指向角偏移量α1和β1的获取方法包括以下步骤:
S61:采用互相关算法获取图像f(x,y)与图像g(x,y)之间的整数级平移量(row_shift,col_shift),其中,row_shift为行方向上的平移量,col_shift为列方向上的平移量;
S62:计算上采样矩阵in,其中,上采样矩阵的行数和列数分别为nr和nc:
in=buf2ft.*conj(buf1ft)
其中,buf2ft为图像f(x,y)的傅里叶变换,buf1ft为图像g(x,y)的傅里叶变换,conj为matlab中的复共轭函数,.*表示两个矩阵按位置相乘;
S63:计算行方向卷积核kernr和列方向卷积核kernc:
kernr=exp((-i×2π/(nr×usfac))×([0:nor-1].'-roff)×(ifftshift([0:nr-1])-floor(nr/2)))
kernc=exp((-i×2π/(nc×usfac))×(tmpnc-floor(nc/2))×([0:noc-1]-coff))
其中,i为虚部单位,roff=dftshift-row_shift×usfac,coff=dftshift-col_shift×usfac,dftshift为中间变量,usfac为设定的上采样因子,且dftshift=fix(ceil(usfac×1.5)/2),fix和ceil为matlab中的取整函数,floor为matlab中的向下取整函数,iffshift为matlab中的傅里叶反变换函数,nor=noc=usfac×1.5,exp为matlab中的底数为e的指数函数,tmpnc为ifftshift([0:nc-1])的转置,.'表示转置;
S64:计算上采样矩阵in的快速傅里叶变换out:
out=kernr×in×kernc
S65:计算上采样相位相关结果CC:
CC=conj(out)/(md2×nd2×usfac2)
其中,md2=m/2,nd2=n/2,m和n分别为模板图像f(x,y)的行数和列数,conj为matlab中的复共轭函数;
S66:获取图像f(x,y)与图像g(x,y)之间的最终平移量:
row_shift^=row_shift+rloc/usfac
col_shift^=col_shift+cloc/usfac
其中,row_shift^为行方向上的最终平移量,col_shift^为列方向上的最终平移量,rloc与cloc分别为上采样相位相关结果CC中最大值所在的行号和列号;
S67:根据最终平移量计算主相机相对于星地相机对地成像模块的指向角偏移量α1和β1:
α1=θ×row_shift^
β1=θ×col_shift^
其中,θ为主相机的单位像元的指向角。
3.如权利要求2所述的一种基于星地相机成像系统的视轴指向解算方法,其特征在于,步骤S67中单位像元的指向角θ的计算方法为:
θ=arctan(a/f)
其中,f为主相机的焦距,a为主相机的单位像元大小。
4.如权利要求1所述的一种基于星地相机成像系统的视轴指向解算方法,其特征在于,步骤S3中所述的视轴矢量v天的获取方法包括以下步骤:
S31:基于内插细分技术,采用加权平方的质心提取方法从星空图像Img1提取N个亚像素级的星点质心位置;
S32:获取N个星点质心位置在导航星表中的赤经、赤纬坐标(αi,δi),i=1,2,...,N,并将赤经、赤纬坐标转换为单位向量坐标(cosαicosδi,sinαicosδi,sinδi)
S33:采用QUEST算法,基于N个星点质心位置的单位向量坐标解算出星地相机相对于天体坐标系J2000的姿态矩阵C;
S34:根据姿态矩阵获取星地相机对天成像模块在天体坐标系J2000下的视轴矢量v天=C-1×v1,其中,C-1为姿态矩阵C的逆矩阵,v1为星地相机对天成像模块视轴的单位矢量。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010854464.4A CN114187345B (zh) | 2020-08-24 | 2020-08-24 | 一种基于星地相机成像系统的视轴指向解算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010854464.4A CN114187345B (zh) | 2020-08-24 | 2020-08-24 | 一种基于星地相机成像系统的视轴指向解算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114187345A true CN114187345A (zh) | 2022-03-15 |
CN114187345B CN114187345B (zh) | 2024-08-20 |
Family
ID=80600640
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010854464.4A Active CN114187345B (zh) | 2020-08-24 | 2020-08-24 | 一种基于星地相机成像系统的视轴指向解算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114187345B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114088060A (zh) * | 2020-08-24 | 2022-02-25 | 中国科学院长春光学精密机械与物理研究所 | 一种用于光学遥感卫星指向测量的星地相机成像系统 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105371870A (zh) * | 2015-11-27 | 2016-03-02 | 北京空间飞行器总体设计部 | 一种基于星图数据的星敏感器在轨精度测量方法 |
CN105444780A (zh) * | 2015-11-30 | 2016-03-30 | 上海卫星工程研究所 | 星载摆扫式光学相机图像定位验证系统及处理方法 |
CN105547258A (zh) * | 2016-01-26 | 2016-05-04 | 北京空间飞行器总体设计部 | 一种遥感卫星tdiccd相机积分时间星上优化计算方法 |
CN105698764A (zh) * | 2016-01-30 | 2016-06-22 | 武汉大学 | 一种光学遥感卫星影像时变系统误差建模补偿方法及系统 |
CN106767714A (zh) * | 2016-11-29 | 2017-05-31 | 上海卫星工程研究所 | 提高卫星图像定位精度的等效失配模型多级标校方法 |
-
2020
- 2020-08-24 CN CN202010854464.4A patent/CN114187345B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105371870A (zh) * | 2015-11-27 | 2016-03-02 | 北京空间飞行器总体设计部 | 一种基于星图数据的星敏感器在轨精度测量方法 |
CN105444780A (zh) * | 2015-11-30 | 2016-03-30 | 上海卫星工程研究所 | 星载摆扫式光学相机图像定位验证系统及处理方法 |
CN105547258A (zh) * | 2016-01-26 | 2016-05-04 | 北京空间飞行器总体设计部 | 一种遥感卫星tdiccd相机积分时间星上优化计算方法 |
CN105698764A (zh) * | 2016-01-30 | 2016-06-22 | 武汉大学 | 一种光学遥感卫星影像时变系统误差建模补偿方法及系统 |
CN106767714A (zh) * | 2016-11-29 | 2017-05-31 | 上海卫星工程研究所 | 提高卫星图像定位精度的等效失配模型多级标校方法 |
Non-Patent Citations (1)
Title |
---|
张廷华;高昆;倪国强;卢岩;: "多传感器的天文标定与空间配准", 光学精密工程, no. 06, 15 June 2017 (2017-06-15) * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114088060A (zh) * | 2020-08-24 | 2022-02-25 | 中国科学院长春光学精密机械与物理研究所 | 一种用于光学遥感卫星指向测量的星地相机成像系统 |
CN114088060B (zh) * | 2020-08-24 | 2024-04-30 | 中国科学院长春光学精密机械与物理研究所 | 一种用于光学遥感卫星指向测量的星地相机成像系统 |
Also Published As
Publication number | Publication date |
---|---|
CN114187345B (zh) | 2024-08-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Jiang et al. | Geometric calibration and accuracy assessment of ZiYuan-3 multispectral images | |
McCormac et al. | DONUTS: A science frame autoguiding algorithm with sub-pixel precision, capable of guiding on defocused stars | |
Piatek et al. | Proper motions of dwarf spheroidal galaxies from Hubble Space Telescope imaging. I. Method and a preliminary measurement for fornax | |
CN107144293A (zh) | 一种视频卫星面阵相机的几何定标方法 | |
Jiang et al. | Detection and correction of relative attitude errors for ZY1-02C | |
CN101750067B (zh) | 一种成像式地球敏感器地球扁率修正方法 | |
Xiong et al. | A generic method for RPC refinement using ground control information | |
CN110211054A (zh) | 一种星载推扫式光学传感器无畸变影像制作方法 | |
Chen et al. | Calibration and validation of ZY-3 optical sensors | |
CN107018402B (zh) | 裸眼三维显示器的图像误差修正方法及装置 | |
CN111238540A (zh) | 一种基于恒星拍摄的珞珈一号相机-星敏安装标定方法 | |
US9794483B1 (en) | Video geolocation | |
CN110751601A (zh) | 一种基于rc光学系统的畸变校正方法 | |
CN114187345B (zh) | 一种基于星地相机成像系统的视轴指向解算方法 | |
CN105222777B (zh) | 用假设经纬度法观测天体定位的方法 | |
Klaus et al. | Camera calibration from a single night sky image | |
CN105023281B (zh) | 基于点扩散函数波前修正的星点像质心计算方法 | |
KR100870894B1 (ko) | 선형 스캐닝 영상의 자동 영상기하 보정 방법 | |
Meyer et al. | Astrometry with the MCAO instrument MAD-An analysis of single-epoch data obtained in the layer-oriented mode | |
CN111366149B (zh) | 一种基于月相的太阳矢量方向的提取方法 | |
CN111879299B (zh) | 一种地基望远镜全自动指星方法 | |
Zacharias et al. | Accurate optical positions of extragalactic radio reference frame sources | |
CN111895968B (zh) | 一种高精度复合水平测角系统 | |
Ren et al. | Advanced Multiple Aperture Seeing Profiler | |
CN109029379A (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 |