一种窄重叠双视场高光谱成像仪的图像配准的方法

摘要

一种窄重叠双视场高光谱成像仪的图像配准的方法，首先利用几何校正后图像的地面坐标粗略的计算左右视场对应关系，以此对应关系为基准在一定的相对位置范围内计算互信息，并去除掉重叠像元数较少时计算的互信息，在剩下的所有互信息中找到最大互信息，则最大互信息所对应的左右视场图像的对应关系即为图像最佳匹配的对应关系。本方法利用了图像几何校正后的地面坐标和互信息来计算图像的匹配关系，可以解决几何校正后地面坐标不精确带来的拼接误差问题和基于互信息配准时误匹配的问题，快速准确的得到图像的配准关系。

主权项

一种窄重叠双视场高光谱成像仪左右视场图像配准的方法，通过计算机进行图像数据读入、分析及计算处理，该方法步骤如下：1)读取所需参数：从图像的头文件中分别读取左、右视场两幅图像左上角的地面坐标(X₁,Y₁)和(X₂,Y₂)、像元数Samples_L和Samples_R、行数Lines_L和Lines_R及地面分辨率r；2)计算两幅图像左上角点在X方向和Y方向上相差的像元数(SX_d,SY_d)：$\left\{\begin{array}{c}{\mathrm{SX}}_d=\frac{X_1-X_2}{r}\\ {\mathrm{SY}}_d=\frac{Y_1-Y_2}{r}\end{array}\right.$3)设置两幅图像的配准范围：设X方向的循环范围为(minoffset_x,maxoffset_x)，XX为X方向的循环范围的循环变量，Y方向的循环范围为(minoffset_y,maxoffset_y)，YY为Y方向的循环范围的循环变量；4)计算两幅图像最佳匹配位置：a)选取待计算图像数据：在左、右视场两幅图像中选取对应的N个波段，该N个波段应为有效且图像清晰的波段，将这N个波段的图像数据读入内存；b)设波段数的循环变量bandx的初始值为1；c)读取N个波段数据中第bandx波段数据；d)设XX的初始值为minoffset_x；e)设YY的初始值为minoffset_y；f)设左右视场两幅图像中的重合像元数count的初始值为0，左右视场两幅图像中的像元的最小最大值分别为L_min＝0、L_max＝99999999和R_min＝0、R_max＝99999999；g)获取计算互信息的图像范围：若SX_d＜0则得到X方向范围取值的最小值(即计算互信息图像范围的X方向范围的最小值)x_min＝‑SX_d+1+XX，最大值为x_max＝Samples_L，否则X方向范围取值的最小值x_min＝1，最大值为x_max＝Samples_R‑SX_d+XX；若SY_d＞0则得到Y方向范围的最小值(即计算互信息图像范围的Y方向范围的最小值)y_min＝SY_d+1+YY，最大值为y_max＝Lines_L，否则Y方向范围的最小值y_min＝1，最大值为y_max＝Lines_R+SY_d+YY；h)设Y方向范围的循环变量linex的初始值为y_min；i)设X方向范围的循环变量samplex的初始值为x_min；j)若左视场图像第linex行第samplex个像元的值不为0且右视场图像的第linex‑SY_d‑YY行第samplex+SX_d‑XX个像元的值不为0，则转入k)，否则则转入l)；k)重合像元数count加1，若左视场图像第linex行第samplex个像元的值大于左视场像元最大值L_max，则此像元值赋给L_max，若左视场图像第linex行第samplex个像元的值小于左视场像元最小值，则此像元值赋给L_min，若右视场图像的第linex‑SY_d‑YY行第samplex+SX_d‑XX个像元的值大于右视场像元最大值，则此像元值赋给R_max，若右视场图像的第linex‑SY_d‑YY行第samplex+SX_d‑XX个像元的值小于右视场像元最小值，则此像元值赋给R_min；l)samplex加1，若samplex大于x_max则转入m)，否则则转入j)；m)linex加1，若linex大于y_max则转入n)，否则则转入j)；n)判断是否符合计算互信息的条件：若重合像元数count小于Lines_L*Lines_R*CHSampls的90％，其中CHSamples是已知参数由窄重叠双视场高光谱成像仪硬件所决定的左右视场一行像元的重合像元数，则转入(r)；否则则转入o)；o)设左视场重合像元直方图级数为Llevel，右视场重合像元直方图级数为Rlevel；计算得到左视场重合像元直方图步长为Lstep＝(L_max–L_min)/(Llevel‑1)，右视场重合像元直方图步长为Rstep＝(R_max–R_min)/(Rlevel‑1)，则可以把左视场重合像元的像元值从小到大分为Llevel个区间，把右视场重合像元的像元值从小到大分为Rlevel个区间；p)计算重合像元落在左视场像元值的Llevel个区间中每一个区间的像元数，并将每个区间的像元数除以count，即得到像元落在不同区间的概率，记为p_i,(1≤i≤Llevel),同理计算重合像元落在右视场不同区间的概率，记为p_j,(1≤j≤Rlevel)；计算重合像元值落在左视场像元值第i个区间同时落在右视场第j个区间的像元数，并处以count，记为p_ij,(1≤i≤Llevel,1≤j≤Rlevel)；q)根据以下公式计算互信息MI_XX,YY,(minoffset_x≤XX≤maxoffset_x，minoffset_y≤YY≤maxoffset_y)：${\mathrm{MI}}_{XX,YY}=(-\underset{i=0}{\overset{Llevel-1}{\Sigma }}{p}_i{\mathrm{logp}}_i)+(-\underset{j=0}{\overset{Rlevel-1}{\Sigma }}{p}_j{\mathrm{logp}}_j)-(-\underset{i=0}{\overset{Llevel-1}{\Sigma }}\underset{j=0}{\overset{Rlevel-1}{\Sigma }}{p}_{ij}{\mathrm{logp}}_{ij})$r)YY加1，若YY大于maxoffset_y则转入s)，否则则转入f)；s)XX加1，若XX大于minoffset_x则转入t)，否则则转入f)；t)在记录下的所有互信息MI_XX,YY，(minoffset_x≤XX≤maxoffset_x，minoffset_y≤YY≤maxoffset_y)中，找到最大的互信息，记录下最大互信息的下标XX和YY的值即为匹配位置，用(X_mi_bandx,Y_mi_bandx)表示，其中(1≤bandx≤N)；u)bandx加1，若bandx大于N则转入v)，否则则转入c)；v)取所有记录下的匹配位置的平均值，即为所求的最佳匹配位置，我们记为X_last和Y_last；5)计算左右视场对应像元关系：①设X方向偏移像元数为offset_x，则offset_x＝SX_d‑X_last，设Y方向偏移像元数为offset_y，则offset_y＝SY_d+Y_last；②设左视场和右视场的第m列、第n行点的坐标分别为(m_L,n_L)和(m_R,n_R)，则得到左右视场像元对应关系为：$\left\{\begin{array}{c}{m}_L=m_R-offset\_x+1\\ n_L=n_R-offset\_y+1\end{array}\right.$上述得到的左右视场对应关系即为本发明方法所求的配准关系。