多尺度光谱TDICMOS点线配准的电荷累积方向测试系统

多尺度光谱tdicmos点线配准的电荷累积方向测试系统
技术领域
1.本发明涉及空间光学遥感技术领域，具体涉及一种相机焦面tdicmos多尺度图像配准与电荷累积量计算的方法。


背景技术：

2.多光谱遥感技术广泛应用于农业监测、矿物探测、军事侦察等领域，是一种高精度地物探测手段。近年来，天基遥感探测不断朝着大视场、高宽幅的方向拓展，对cmos的成像幅宽有了更高的要求。为了满足该要求，高分辨率相机一般采用多个片元错位拼接的方式来实现大视场成像。但错位拼接会引入偏流角、相移等一系列误差。
3.由于各个谱段感光尺寸相同且间距相等，卫星采用正向推扫的方式成像，会导致各谱段时域和空域并不一一对应，各谱段对同一目标点的成像的时间并不一致。因此，需一种高精度配准方式去实现多尺度、多光谱之间的时空域的配准。现有的图像配准算法大多基于特征点匹配，容易忽略图像中直线、几何形状等特征，造成图像的直线发生曲折，画面有撕裂感。多光谱的图像配准容易出现特征点失调，经常会出现匹配失败的情况，因此，需要一种全新的单谱拼接和多谱间配准方法。
4.像移是造成遥感成像质量下降的重要因素，依据配准结果反解出像移之间的关系，为消除像移影响提供重要的数据支撑。电荷累积方向是像移直接表现在cmos上的形式，分析其大小对成像质量的影响意义重大。


技术实现要素：

5.本发明为解决现有图像配准方法忽略图像中直线、几何形状等特征，造成图像的直线发生曲折，画面有撕裂感，导致图像配准易出现特征点失调，匹配失败等问题，提供一种多尺度光谱tdicmos点线配准的电荷累积方向测试系统。
6.多尺度光谱tdicmos点线配准的电荷累积方向测试系统，该系统包括分光式多光谱tdicmos多级积分成像单元，多尺度片内点线结合配准单元，光谱维度的多尺度轮廓特征配准单元以及电荷累积方向计算单元；
7.所述分光式多光谱tdicmos多级积分成像单元由三段cmos探测器的片元拼接组成，采用单摆运动模式推扫采集每个cmos探测器的五个谱段的高分辨率图像；
8.所述多尺度片内点线结合配准单元接收所述五个谱段的高分辨率图像，采用特征点和特征线相结合配准的方法分别实现相同谱段图像精确配准拼接，获得五个独立谱段的拼接图像；
9.所述光谱维度的多尺度轮廓特征配准单元接收五个独立谱段的拼接图像，采用多尺度图像轮廓特征匹配的方式，实现各个谱段之间的配准；
10.所述电荷累积方向计算单元根据所述光谱维度的多尺度轮廓特征配准单元输出的配准图像与原始图像的特征点位置偏移量，计算每个谱段每个位置曝光时间偏移量，并作为自变量，导出像移与时间的函数关系，进一步获得电荷累积方向与时间的关系。
11.本发明的有益效果：本发明所述的测试系统，在模拟真实的推扫模式下，采用反射分光排布式tdicmos采集原始的多光谱图像，相比传统错位拼接方式，采取反射拼接的方式会有效减少三段图像的偏流角误差，提高电荷矢量方向计算精度；依据片元图像空间特征的多样性，采用特征点线结合的单谱配准方法对原始图像进行配准拼接，解决只基于特征点匹配方式所导致拼接图像几何特征扭曲的问题；依据多光谱图像分辨率差异较大以及辐射强度多样性较强的特征，采用基于css检测改进型多尺度多光谱配准方法实现对于频谱的特征点精确对准，解决多光谱图像配准容易失调的问题；最终根据各光谱图像匹配点的坐标精确计算出不同时差所产生的像移大小，分析像移实时变化，根据实时像移矢量的大小计算出电荷累积方向。依据先后顺序评估每一个环节的精度，实现多光谱tdicmos图像空域和频域的联合配准与电荷累积方向快速测量。
附图说明
12.图1为本发明所述的多尺度光谱tdicmos点线配准的电荷累积方向测试系统的原理图；
13.图2为cmos探测器的结构示意图；
14.图3为卫星成像模式原理图；
15.图4为单谱图像拼接原理图；
16.图5为基于点线结合的图像配准原理图；
17.图6为基于轮廓特征多光谱图像配准原理图。
具体实施方式
18.结合图1至图6说明本实施方式，如图1所示，多尺度光谱tdicmos点线配准的电荷累积方向测试系统，包括分光式多光谱tdicmos多级积分成像单元，多尺度片内点线结合配准单元，光谱维度的多尺度轮廓特征配准单元以及电荷累积方向计算单元；
19.所述分光式多光谱tdicmos多级积分成像单元的单个tdicmos的感光区域如图2所示。每一个cmos的面元的五个谱段呈纵向排布，分别是p、b1、b2、b3和b4谱段，其中全色谱段p对应的最小像元尺寸为7um，分辨率为8224*1；其它四个为多光谱段，最小像元尺寸为14um，分辨率为2056*1。遥感卫星依照图3所示的摆扫模式进行时间积分成像。所述分光式多光谱tdicmos多级积分成像单元结构包含一个全色谱段p，四个多光谱段b。其中p谱与b谱分辨率相差4倍。光学系统像面采用三段cmos片元拼接而成，相机按单摆运动模式推扫采集五个谱段的三节高分辨率图像。
20.所述多尺度片内点线结合配准单元接收成像单元所成的高分辨率图像，采用特征点和特征线相结合配准的算法，分别对五个谱段的三节图像进行重叠区域的匹配点计算，实现每个谱三段子图像精确配准拼接，合成五幅相互独立的不同谱段的拼接图像。
21.本实施方式中，每个谱的图像均由三节待拼接子图组成。若将三节图像直接按几何关系拼接时，由装调或者机械振动引起的误差会导致图像配准精度下降。因此，需要在图像处理阶段采用高精度配准算法对图像进行配准。
22.遥感成像的地面目标中包含大量的规则的直线或者矩形特征，为保证两张子图拼接后直线不会发生曲折。采用特征点和特征线相结合配准的算法，特征点用于计算全局仿
射变换矩阵，特征线用于精细分割区域，调整每个区域内的变换矩阵系数，实现各个部位精确配准。对于两张待拼接的图像，采用surf快速检测出特征点对集合p1,p2,p3,p4、p1',p2',p3',p4'，其中p
i
(i＝1，2，3，4)与p
i
'一一对应，在实际的对应的空间位置上它们属于同一点。求解出全局仿射变换矩阵。采用最小二乘法拟合特征点集合，找出全局最大的直线特征矢量l
i
、l
j
'，依据点线关系做进一步配准。如图5所示，若右图特征点序列p
i
(i＝1，2，3，4)构成的矢量l
i
在右图的投影矢量l
i
'必须旋转一定的角度才能与左图全局最大特征矢量l
j
'匹配，则两张图像间出现了配准直线失准的情况。针对全局仿射变换直线失准的问题，必须建立以直线约束为依据的网格能量函数，对特征直线l
i
'左右两侧分别求取子区域的坐标变换矩阵，进而估计图像之间的最优配准变换。
23.所述光谱维度的多尺度轮廓特征配准单元接收点线配准单元所成的像，针对不同光谱图像的存在像移问题，采用多尺度图像轮廓特征匹配的方式实现各个谱段之间的配准。多尺度图像轮廓特征匹配为对分辨率较低b谱图像超分辨重构，然后与p谱图像进行基于最大梯度轮廓特征点插值算法的精细匹配。
24.本实施方式中，面元p谱比b谱的分辨率高四倍，即使两谱段在空间位置上基本配准，但配准精度很低，会出现重影现象。采用超分辨技术对b谱图像进行重构，将问题转换为同尺度多光谱图像间的配准；地面目标的光谱特性复杂多样，不同谱段图像的特征点容易失调，但空间信息保留较好。采用改进型的css轮廓及角点检验算法，找出位于各谱图像相似轮廓上的特征点。虽然该特征点数目会大幅下降，但错误匹配的概率也会缩小。最后，为了进一步提高精度，将每相邻四个特征匹配点对角相连，令其交点作为最终的拟合特征点。
25.具体采取的方法为：
26.首先，对b谱图像进行超分辨率重构。采用pixelshuffle上采样方法，将b谱低分辨特征进行卷积重组，得到高分辨率的特征；采用改进型的css检测轮廓及角点，使用适用于多光谱图像特征主方向计算方法计算角点的主方向；提取角点及其描述子，按照序列提取匹配点并去除显著错误匹配；采用特征点插值的方式拟合最终的匹配点，提高匹配精度。
27.所述电荷累积方向计算单元接收光谱图像配准单元的图像，根据配准图像与原始图像的特征点位置偏移量，算出每个谱段每个位置曝光时间偏移量，作为自变量，导出像移与时间的函数关系，进而导出电荷累积方向与时间的关系。
28.如图6所示，由于像移的存在，同一特征点在不同谱段的位置会发生变换。依据推扫成像的模式，面元沿垂直于轨道方向运动，对于地面同一个点目标，b1谱要比b2谱先接收到该点的光强信息，其他谱段依次类推。时间差会导致各谱段的对目标点所成的像存在像移量。建立像插值匹配点追踪模型。记录各谱段匹配点原始坐标a0(λ)＝(x0,y0)，经过图像配准后，配准点的坐标为a(λ)＝(x,y)。则像移量为各谱段图像时间偏移量的计算如下，其中v是cmos推扫的速度。
[0029][0030]
式中，h＝0,h,2h,3h,4h,h为p谱与b1谱感光区域的垂直间隔，λ为各个感光区域所对应中心波段波长。
[0031]
最终，可得到δa(λ)和t(λ)的对应关系。最终依照电荷累积方向与像移量的线性关系得到电荷累积方向矢量的大小。