[发明专利]一种线阵CCD的行进畸变校正方法

权利要求书

1.一种线阵CCD的行进畸变校正方法，其特点在于，按照以下步骤具体实施：

步骤1、去除标准图像与待校正图像中差异较大的图像列

为了描述方便起见，将去除了较大差异列的标准图像和待校正图像仍旧以和表示，其Sobel边缘图仍旧以和表示；

步骤2、分别提取标准图像与待校正图像特征线元

2a）对步骤1差异处理后的标准图像在水平方向sobel边缘和待校正图像在水平方向sobel边缘分别进行行方向投影，得到和如下式：

hk(x)=Σy=1nfx_Edgek(x,y),k=O,R,x=1,2,...,mk,mk=mO,mR,---(3)]]>

2b）对步骤2a）得到的投影曲线和求其均值μ_O和μ_R：

μk=1mkΣx=1mkhk(x),k=O,R,mk=mO,mR,---(4)]]>

2c）按照所求出的均值μ_O和μ_R，对投影曲线和进行分段，曲线的段边界点为均值与投影曲线的交点x_l^O，l＝1,2,...,n_O，曲线的段边界点为均值与投影曲线的交点x_l^R，l＝1,2,...,n^R，即：

x_l^k＝x|h(x)＝μ_k，k＝O,R，l＝1,2,...,n_k，n_k＝n_O,n_R，    （5）

其中的n_O是投影曲线与均值μ_O的相交点的个数，n_R是投影曲线与均值μ_R的相交点的个数，

分段后的曲线段Cu^O(i)，i＝1,2,...,n_O+1，Cu^R(i)，i＝1,2,...,n_R+1为：

Cuk(i)={x|x∈[xi-1k,xik]},x0k=1,xnk+1k=mk,k=O,R,]]>

i＝1,2,...,n_k+1，n_k＝n_O,n_R，m_k＝m_O,m_R；    （6）

2d）分别求出每个曲线段Cu^O(i)和Cu^R(i)的极大值点x_max^O(i)，i＝1,2,...,n_O+1，x_max^R(i)，i＝1,2,...,n_R+1，如下式：

x_max^k(i)＝argmax{Cu^k(i)}，k＝O,R，i＝1,2,...,n_k+1，n_k＝n_O,n_R，    （7）

这些极大值点分别对应了标准图像以及标准图像的边缘待校正图像以及待校正图像的边缘的特征线元所在行的位置；

步骤3、计算标准图像和待校正图像的特征线元的相似性，获得两者相匹配的特征线元对(x_max^O(i),x_max^R(i))，i＝1,2,...,N，其中N为配对成功的特征线元对的个数，具体步骤是：

3a）对步骤1a）得到的标准图像的边缘为用步骤2得到的极大值点x_max^O(i)，i＝1,2,...,n_O+1，作为各个子图的边界点，分成子图i＝1,2,...,n_O+1；

同样地，对步骤1a）得到的待校正图像的边缘为用步骤2得到的极大值点x_max^R(i)，i＝1,2,...,n_R+1，作为各个子图的边界点，分成子图子图i＝1,2,...,n_R+1，

各个子图的分割计算公式如下：

i＝1,2,...,n_k+1，n_k＝n_O,n_R，    （8）

3b）计算步骤3a）得到的各个子图的投影曲线i＝1,2,...,n_O+1和i＝1,2,...,n_R+1，公式如下：

hki(y)=Σx=1xmaxk(i)fx_Edgeki(x,y),ki=Oi,Ri,k=O,R,]]>

i＝1,2,...,n_k+1，n_k＝n_O,n_R（9）

3c）对子图的投影曲线i＝1,2,...,n_O+1和i＝1,2,...,n_R+1计算其相关系数矩阵公式如下：

cor(i,j)=Σy=1n[hOi(y)-h‾Oi][hRj(y)-h‾Rj]Σy=1n[hOi(y)-h‾Oi]2·Σy=1n[hRj(y)-h‾Rj]2,]]>

i＝1,2,...,n_O+1，j＝1,2,...,n_R+1，     （10）

其中的是i＝1,2,...,n_O+1的均值，是j＝1,2,...,n_R+1的均值；

3d）找出中的最大值点(i_max,j_max)，公式如下：

(imax,jmax)=argmaxi=1,...,nO+1maxj=1,...,nR+1{cor(i,j)},---(11)]]>

由此得到标准图像与待校正图像的配对特征线元位置坐标为(i_max⁽¹⁾,j_max⁽¹⁾)=(i_max,j_max)，即标准图像的第i_max⁽¹⁾行与待校正图像的第j_max⁽¹⁾行为相同的特征线元位置坐标；

3e）根据(i_max⁽¹⁾,j_max⁽¹⁾)将矩阵分为两个小矩阵和即：

cor₁(i,j)＝cor(i,j)，i＝1,2,...,i_max⁽¹⁾，j＝1,2,...,j_max⁽¹⁾，    （12）

cor₂(i,j)＝cor(i_max+i,j_max+j)，i＝1,2,...,n_O+1-i_max，j＝1,2,...,n_R+1-j_max，    （13）

3f）将这两个小矩阵分别代入步骤3d）和步骤3e），进行同样的处理后，得到标准图像与待校正图像的配对特征线元位置坐标为(i_max⁽²⁾,j_max⁽²⁾)和(i_max⁽³⁾,j_max⁽³⁾)，并且又分别得到两个小矩阵

[cor11(i,j)]imax(2)×jmax(2),[cor12(i,j)](imax(1)+1-imax(2))×(imax(1)+1-jmax(2)),]]>

[cor21(i,j)]imax(3)×jmax(3),[cor12(i,j)](nO+1-imax(1)+1-imax(3))×(nR+1-imax(1)+1-jmax(3)),]]>

之后，再将得到的四个小矩阵分别代入步骤3d）和步骤3e），直到所有的小矩阵只有两行两列为止，再设定得到了N对标准图像与待校正图像的配对特征线元位置序列为(i_max^(k),j_max^(k))，k＝1,2,,...,N；

步骤4、采用匹配的特征线元对，对待校正图像进行缩放

由特征线元位置坐标对(i_max^(k),j_max^(k))，k＝1,2,...,N，

将待校正图像分割为子图像块[fRk(x,y)](jmax(k)-jmax(k-1))×n=]]>fR(x,y)|x∈[jmax(k-1),jmax(k)],]]>y＝1,2,...,n，其大小为(j_max^(k)-j_max^(k-1))×n，j_max⁽⁰⁾＝1，j_max^(N+1)＝m_R，k＝1,2,...,N+1；

将标准图像分割为子图像块[fOk(x,y)](imax(k)-imax(k-1))×n=]]>fO(x,y)|x∈[imax(k-1),imax(k)],]]>y＝1,2,...,n，其大小为(i_max^(k)-i_max^(k-1))×n，i_max⁽⁰⁾＝1，i_max^(N+1)＝m_O，k＝1,2,...,N+1；

按照特征线元位置坐标对(i_max^(k),j_max^(k))，k＝1,2,...,N的对应关系，将带校正图的子图像块k＝1,2,...,N从大小为(j_max^(k)-j_max^(k-1))×n，缩放至与标准图像对应子块k＝1,2,...,N+1相同大小，即(i_max^(k)-i_max^(k-1))×n，得到计算公式如下：

f^Rk(x,y)=fRk(ck,x,y),]]>

x＝i_max^(k-1)，i_max^(k-1)+1，...，i_max^(k)，y＝1，2，...，n，k＝1，2，...，N+1，(14)

其中，ck=imax(k)-imax(k-1)jmax(k)-jmax(k-1),---(15)]]>

由此，完成对待校正图像的初匹配，处理后的待校正图像与标准图像大小一致，即为：

f^R(x,y)=f^Rk(x,y),x=imax(k-1),imax(k-1)+1,...,imax(k),k=1,2,...,N+1;---(16)]]>

步骤5、对得到的初匹配图像采用交叉寻找的方法进行匹配，得到细匹配的匹配线元

5a）对经过步骤4得到的图像以及标准图像分别求其Sobel水平边缘图，并分别以和表示；

5b）逐行计算带校正图与标准图像的行累积差e(x)，x＝1,2,...,m_O，公式如下：

e(x)=1n·Σy=1n|f^y_EdgeR(x,y)-fy_EdgeO(x,y)|,x=1,2,...,mO;---(17)]]>

5c）将计算得到的行累积差e(x)，x＝1,2,...,m_O，设定阈值是的均值，α是调整系数；

步骤6、在得到的匹配线元的基础上，对标准图像和待校正图像，按照步骤4的方法进行缩放，便可得到最终的畸变校正结果，即可进行对应的比对，得出车辆底盘异物的检测结论。