專利名稱：一種dmc彩色合成影像的輻射再處理方法
技術(shù)領(lǐng)域：
本發(fā)明屬于圖像處理與分析領(lǐng)域，涉及一種DMC影像的輻射再處理方法。
背景技術(shù)：
航空數(shù)碼相機的應(yīng)用，使得從影像存儲到測圖完全采用數(shù)字的形式進行成為可能，與膠片相機相比，極大的降低了成本，縮短了信息處理的周期，提高了對地觀測的實效性，使其面臨著前所未有的需求和廣闊的應(yīng)用前景。然而在有限的飛行控制條件下，為了達到航空攝影重疊度以及獲取影像高分辨率的要求，出于技術(shù)水平的限制以及成本的考慮，目前直接生產(chǎn)大幅面的面陣CCD相機還有困難，因此出現(xiàn)了一些基于面陣組合技術(shù)的航空數(shù)碼相機，如Z/I Imaging的DMC(Digital Mapping Camera)，它通過將多臺小面陣CCD組合構(gòu)成虛擬的大面陣CCD來達到在具有高幾何分辨率的同時具有較大的地面覆蓋范圍。
DMC航空數(shù)碼相機由八個同步工作的CCD相機組合而成，4個多光譜相機(紅、綠、藍和近紅外)和4個全色相機(其獲得的4幅全色影像之間有一定的重疊，便于鑲嵌處理)，每個CCD相機都有自己的鏡頭。DMC系統(tǒng)在首先對各鏡頭獲取的影像進行幾何和輻射校正，將四幅全色影像鑲嵌成一個大幅面的虛擬影像；然后將虛擬影像和紅、綠、藍三個波段的影像進行融合處理獲得DMC彩色合成(Color composite)影像。
DMC彩色合成影像由DMC系統(tǒng)提供的后處理軟件包獲得，在大多數(shù)情況下影像質(zhì)量較好，但在有些情況下由于不同相機的CCD在成像時響應(yīng)的差異性等各種因素的影響，生成的DMC彩色合成影像在鑲嵌處理的拼接線附近的色彩過渡有時仍然不夠平滑，存在殘余輻射差異，這給影像的進一步應(yīng)用造成了不同程度的影響。由于存在這種現(xiàn)象的影像沒有明顯的接縫，常用的鑲嵌接縫消除方法并不適用，因此需要研究新的處理方法以改善DMC彩色合成影像的質(zhì)量。

發(fā)明內(nèi)容
本發(fā)明所要解決的問題是提供一種DMC彩色合成影像的輻射再處理方法，該方法可以有效的解決DMC彩色合成影像存在的色彩過渡不夠平滑的問題，消除殘余輻射差異，改善影像質(zhì)量。
本發(fā)明提供的技術(shù)方案是一種DMC彩色合成影像的輻射再處理方法，對DMC彩色合成影像中的每一個色彩過渡區(qū)域按如下步驟分別進行處理一、色彩過渡區(qū)域和拼接線的自動定位a、建立影像金字塔在DMC彩色合成影像中，對每兩個相鄰的CCD影像之間的矩形區(qū)域，采用像元平均的方法建立3-5層金字塔影像層；每兩個相鄰的CCD影像之間的矩形區(qū)域，當(dāng)兩個相鄰的CCD影像為左右相鄰時，矩形區(qū)域的高度為CCD影像的高度，寬度取相鄰CCD影像中間線兩側(cè)各30-300像素，當(dāng)兩個CCD影像為上下相鄰時，矩形區(qū)域?qū)挾葹镃CD影像的寬度，高度為相鄰CCD影像中間線兩側(cè)各30-300像素；
b、影像分塊將每兩個相鄰的CCD影像之間由上一步確定的矩形區(qū)域沿拼接線方向，將其分成大小相同、互不重疊的1-50個矩形影像塊；c、利用金字塔影像層確定過渡區(qū)域和拼接線的具體位置a)、對每個金字塔影像層，在每個影像塊中，取平行于拼接線的各影像行或者影像列的像素的亮度平均值形成一條一維曲線；b)、對每條一維曲線進行二階差分計算；c)、根據(jù)最上層的金字塔影像層計算色彩過渡區(qū)域和拼接線位置對于最上層的金字塔影像層，將相鄰CCD影像的中間線作為本次計算的拼接線的初始位置，中間線兩側(cè)各1-10個像素寬度范圍為色彩過渡區(qū)域的初始范圍；在每個影像塊中，分別計算當(dāng)前金字塔影像層中色彩過渡區(qū)域兩側(cè)區(qū)域的像素值的平均值；若左側(cè)平均值大于右側(cè)平均值，在該影像塊當(dāng)前金字塔影像層對應(yīng)的一維曲線的二階差分的計算結(jié)果中，在拼接線的初始位置左側(cè)搜索最小值，右側(cè)搜索最大值，若左側(cè)平均值小于右側(cè)平均值，則在該影像塊當(dāng)前金字塔影像層對應(yīng)的一維曲線的二階差分的計算結(jié)果中，在拼接線的初始位置左側(cè)搜索最大值，右側(cè)搜索最小值；d)、利用金字塔影像層，將上一層計算得到的色彩過渡區(qū)域和拼接線位置作為本次計算的色彩過渡區(qū)域和拼接線的初始位置；在每個影像塊中，分別計算當(dāng)前金字塔影像層中色彩過渡區(qū)域兩側(cè)區(qū)域的像素值的平均值；設(shè)色彩過渡區(qū)域的初始范圍在一維曲線中對應(yīng)的區(qū)間為[xA，xB]，xA、xB為對應(yīng)的與拼接線平行的影像列或影像行的編號，若左側(cè)平均值大于右側(cè)平均值，則在該影像塊當(dāng)前金字塔影像層對應(yīng)的一維曲線的二階差分的計算結(jié)果中，在xA左右各1-10個像素范圍內(nèi)搜索最小值，同時保證搜索范圍在初始拼接線的左側(cè)，在xB左右各1-10個像素范圍內(nèi)搜索最大值，同時保證搜索范圍在初始拼接線的右側(cè)，若左側(cè)平均值小于右側(cè)平均值，則在該影像塊當(dāng)前金字塔影像層對應(yīng)的一維曲線的二階差分的計算結(jié)果中，在xA左右各1-10個像素范圍內(nèi)搜索最大值，同時保證搜索范圍在初始拼接線的左側(cè)，在xB左右各1-10個像素范圍內(nèi)搜索最小值，同時保證搜索范圍在初始拼接線的右側(cè)；e)根據(jù)各影像塊計算得到的色彩過渡區(qū)域和拼接線的位置決定相鄰CCD影像間色彩過渡區(qū)域和拼接線的位置對每個影像塊色彩過渡區(qū)域位置進行投票，取出現(xiàn)頻數(shù)最高的色彩過渡區(qū)域和拼接線的位置作為最后的相鄰CCD影像間的色彩過渡區(qū)域[XA，XB]和拼接線的位置(XA+XB)/2，XA，XB為對應(yīng)的出現(xiàn)頻數(shù)最高的與拼接線平行的影像列或影像行的編號；二、“打斷”處理“打斷”處理按下面步驟在紅、綠、藍三個波段上分別進行a、計算相鄰CCD影像之間的灰度差異值將色彩過渡區(qū)域沿拼接線方向上分成大小相等、互不重疊的矩形區(qū)域，這樣拼接線也同時被分成了大小相等、互不重疊的小段，使每小段拼接線長度為10-100像素，然后分別計算每個小段拼接線所在的色彩過渡區(qū)域兩側(cè)3-100像素范圍內(nèi)灰度值的均值CiL、CiR，其中i為該小段拼接線的編號，則第i段拼接線兩側(cè)CCD影像之間存在的灰度差異值為ΔCi＝CiR-CiLb、平滑各段拼接線兩側(cè)CCD影像之間存在的灰度差異值采用一維均值濾波器對求得的灰度差異值進行平滑，濾波器窗口大小取3至51；c、在拼接線兩側(cè)分別改正相鄰CCD影像之間的灰度差異值以垂直于拼接線方向的影像行或影像列為單位，在拼接線兩側(cè)各40-3000像素范圍內(nèi)對每個像素的灰度值按下述公式進行改正 其中Iin表示處理前像素灰度值；Iout表示處理后像素灰度值；d為處理像素到拼接線的距離；w為拼接線一側(cè)的改正寬度；ΔC′為該像素所在的垂直于拼接線方向的影像行或影像列需要改正的灰度差異值，該值通過與該影像行或影像列最鄰近的兩個小段拼接線兩側(cè)CCD影像之間的灰度差異值ΔCi及ΔCi+1線性內(nèi)插得到；三、“接縫”消除處理“接縫”消除處理按下面步驟在紅、綠、藍三個波段上分別進行a、將色彩過渡區(qū)域沿拼接線方向上分成大小相等、互不重疊的矩形區(qū)域，這樣拼接線也同時被分成了大小相等、互不重疊的小段，使每小段拼接線長度為1-20像素，然后分別計算每一小段拼接線兩側(cè)1-10像素范圍內(nèi)灰度值的均值FjL、FjR，其中j為該小段拼接線的編號，則第j段拼接線兩側(cè)的灰度差異值為ΔFj＝FjR-FjL。
b、在每小段拼接線所在的色彩過渡區(qū)域內(nèi)，分別計算平行于拼接線的各影像行或影像列的像素灰度值的平均值Mjn，其中n為平行于拼接線的影像行或影像列的編號；c、由下述公式計算各小段拼接線所在的色彩過渡區(qū)域的灰度改正值Δgj=MjXA+ΔFj(n-XA)(XB-XA)-Mjn]]>色彩過渡區(qū)域內(nèi)的每個像素的灰度改正值由與該像素最鄰近的兩個小段拼接線所在的色彩過渡區(qū)域的灰度改正值Δgj和Δgj+1線性內(nèi)插得到；將色彩過渡區(qū)域內(nèi)的每一個像素的灰度值加上該像素的灰度改正值就得到“接縫”消除處理的結(jié)果影像。
本發(fā)明可以有效的解決由于不同相機的CCD在成像時響應(yīng)的差異性等各種因素的影響導(dǎo)致的DMC彩色合成影像存在的色彩過渡不夠平滑的問題，消除殘余輻射差異，明顯改善影像質(zhì)量，使DMC彩色合成影像在其鑲嵌拼接線附近同種地物的色調(diào)區(qū)域一致，為無縫影像數(shù)據(jù)庫的建立和參考判讀、變化檢測分析、分類、目標識別、特征提取等進一步應(yīng)用提供高質(zhì)量的影像數(shù)據(jù)。


圖1為本發(fā)明DMC彩色合成影像的輻射再處理方法的處理流程圖；圖2為DMC彩色合成影像；圖3為采用本發(fā)明對圖2進行處理的效果圖；圖4為圖2的局部放大圖；圖5為圖3對應(yīng)于圖4相應(yīng)區(qū)域的效果圖。
具體實施例方式
參見圖1，本發(fā)明提供一種DMC彩色合成影像的輻射再處理方法，該方法可有效的解決DMC彩色合成影像存在的色彩過渡不夠平滑的問題，消除殘余輻射差異，明顯改善影像質(zhì)量。本發(fā)明對DMC彩色合成影像中的每一個色彩過渡區(qū)域按如下步驟分別進行處理一、色彩過渡區(qū)域和拼接線的自動定位a、建立影像金字塔在DMC彩色合成影像中，對每兩個相鄰的CCD影像之間的矩形區(qū)域，采用像元平均的方法建立3-5層金字塔影像層，如采用9像元平均的方法建立4層金字塔(包括原始影像層)；每兩個相鄰的CCD影像之間的矩形區(qū)域，當(dāng)兩個相鄰的CCD影像為左右相鄰時，矩形區(qū)域的高度為CCD影像的高度，寬度取相鄰CCD影像中間線兩側(cè)各30-300像素，如取兩側(cè)各121像素，當(dāng)兩個CCD影像為上下相鄰時，矩形區(qū)域?qū)挾葹镃CD影像的寬度，高度為相鄰CCD影像中間線兩側(cè)各30-300像素，如取兩側(cè)各121像素；b、影像分塊將每兩個相鄰的CCD影像之間由上一步確定的矩形區(qū)域沿拼接線方向?qū)⑵浞殖纱笮∠嗤⒒ゲ恢丿B的1-50個矩形影像塊(即拼接線也被平均分成了1-50等份)，如平均分成10塊；c、利用金字塔影像層確定過渡區(qū)域和拼接線的具體位置a)、對每個金字塔影像層，在每個影像塊中，取平行于拼接線的各影像行或者影像列的像素的亮度平均值形成一條一維曲線；b)、對每條一維曲線，進行二階差分計算；設(shè)經(jīng)平滑處理的一維曲線為y＝f(x)，則其二階差分為y″＝f(x-1)+f(x+1)-2f(x)，其中x為與拼接線平行的影像列或影像行的編號；c)、根據(jù)最上層的金字塔影像層計算色彩過渡區(qū)域和拼接線位置對于最上層的金字塔影像層，將相鄰CCD影像的中間線作為本次計算的拼接線的初始位置，中間線兩側(cè)各1-10個像素寬度范圍為色彩過渡區(qū)域的初始范圍，如取2個像素；在每個影像塊中，分別計算當(dāng)前金字塔影像層中色彩過渡區(qū)域兩側(cè)區(qū)域的像素值的平均值；若左側(cè)平均值大于右側(cè)平均值，在該影像塊當(dāng)前金字塔影像層對應(yīng)的一維曲線的二階差分的計算結(jié)果中，在拼接線的初始位置左側(cè)搜索最小值，右側(cè)搜索最大值，若左側(cè)平均值小于右側(cè)平均值，則在該影像塊當(dāng)前金字塔影像層對應(yīng)的一維曲線的二階差分的計算結(jié)果中，在拼接線的初始位置左側(cè)搜索最大值，右側(cè)搜索最小值；設(shè)其分別在xA、xB處取得最大值和最小值，xA、xB為對應(yīng)的與平接線平行的影像列或影像行的編號，則xA、xB之間區(qū)域即為色彩過渡區(qū)域，拼接線所在的影像列或影像行的編號為(xA+xB)/2；d)、利用金字塔影像層，將上一層計算得到的色彩過渡區(qū)域和拼接線位置作為本次計算的色彩過渡區(qū)域和拼接線的初始位置；在每個影像塊中，分別計算當(dāng)前金字塔影像層中色彩過渡區(qū)域兩側(cè)區(qū)域的像素值的平均值；設(shè)色彩過渡區(qū)域的初始范圍在一維曲線中對應(yīng)的區(qū)間為[xA，xB]，xA、xB為對應(yīng)的與拼接線平行的影像列或影像行的編號，若左側(cè)平均值大于右側(cè)平均值，則在該影像塊當(dāng)前金字塔影像層對應(yīng)的一維曲線的二階差分的計算結(jié)果中，在xA左右各1-10個像素范圍內(nèi)搜索最小值，如取4個像素，同時保證搜索范圍在初始拼接線的左側(cè)，在xB左右各1-10個像素范圍內(nèi)搜索最大值，如取4個像素，同時保證搜索范圍在初始拼接線的右側(cè)，若左側(cè)平均值小于右側(cè)平均值，則在該影像塊當(dāng)前金字塔影像層對應(yīng)的一維曲線的二階差分的計算結(jié)果中，在xA左右各1-10個像素范圍內(nèi)搜索最大值，如取4個像素，同時保證搜索范圍在初始拼接線的左側(cè)，在xB左右各1-10個像素范圍內(nèi)搜索最小值，如取4個像素，同時保證搜索范圍在初始拼接線的右側(cè)；e)根據(jù)各影像塊計算得到的色彩過渡區(qū)域和拼接線的位置決定相鄰CCD影像間色彩過渡區(qū)域和拼接線的位置對每個影像塊色彩過渡區(qū)域位置進行投票，取出現(xiàn)頻數(shù)最高的色彩過渡區(qū)域和拼接線的位置作為最后的相鄰CCD影像間的色彩過渡區(qū)域[XA，XB]和拼接線的位置(XA+XB)/2，XA，XB為對應(yīng)的出現(xiàn)頻數(shù)最高的與拼接線平行的影像列或影像行的編號；二、“打斷”處理a、計算相鄰CCD影像之間的灰度差異值將色彩過渡區(qū)域沿拼接線方向上分成大小相等、互不重疊的矩形區(qū)域，這樣拼接線也同時被分成了大小相等、互不重疊的小段，使每小段拼接線長度為10-100像素，如取25像素，然后分別計算每個小段拼接線所在的色彩過渡區(qū)域兩側(cè)3-100像素范圍內(nèi)灰度值的均值CiL、CiR，如取30像素，其中i為該小段拼接線的編號，則第i段拼接線兩側(cè)CCD影像之間存在的灰度差異值為ΔCi＝CiR-CiLb、平滑各段拼接線兩側(cè)CCD影像之間存在的灰度差異值采用一維均值濾波器對求得的灰度差異值進行平滑，濾波器窗口大小取3至51，如取33；c、在拼接線兩側(cè)分別改正相鄰CCD影像之間的灰度差異值以垂直于拼接線方向的影像行或影像列為單位，在拼接線兩側(cè)各40-3000像素范圍內(nèi)對每個像素的灰度值按下述公式進行改正，如取1000像素 其中Iin表示處理前像素灰度值；Iout表示處理后像素灰度值；d為處理像素到拼接線的距離；w為拼接線一側(cè)的改正寬度；ΔC′為該像素所在的垂直于拼接線方向的影像行或影像列需要改正的灰度差異值，該值通過與該影像行或影像列最鄰近的兩個小段拼接線兩側(cè)CCD影像之間的灰度差異值ΔCi及ΔCi+1線性內(nèi)插得到；三、“接縫”消除處理a、計算色彩過渡區(qū)域兩側(cè)的灰度差異色彩過渡區(qū)域兩側(cè)的灰度差異通過色彩過渡區(qū)域之外較窄的區(qū)域計算；將色彩過渡區(qū)域沿拼接線方向上分成大小相等、互不重疊的矩形區(qū)域，這樣拼接線也同時被分成了大小相等、互不重疊的小段，使每小段拼接線長度為1-20像素，如取5像素，然后分別計算每一小段拼接線兩側(cè)1-10像素范圍內(nèi)灰度值的均值FjL、FjR，如取3像素，其中j為該小段拼接線的編號，則第j段拼接線兩側(cè)的灰度差異值為ΔFj＝FjR-FjL。
b、在每小段拼接線所在的色彩過渡區(qū)域內(nèi)，分別計算平行于拼接線的各影像行或影像列的像素灰度值的平均值Mjn，其中n為平行于拼接線的影像行或影像列的編號；c、由下述公式計算各小段拼接線所在的色彩過渡區(qū)域的灰度改正值Δgj=MjXA+ΔFj(n-XA)(XB-XA)-Mjn]]>色彩過渡區(qū)域內(nèi)的每個像素的灰度改正值由與該像素最鄰近的兩個小段拼接線所在的色彩過渡區(qū)域的灰度改正值Δgj和Δgj+1線性內(nèi)插得到；將色彩過渡區(qū)域內(nèi)的每一個像素的灰度值加上該像素的灰度改正值就得到“接縫”消除處理的結(jié)果影像。
采用本發(fā)明的上述方法，對圖2的DMC彩色合成影像進行處理后解決了色彩過渡區(qū)域存在的色彩過渡不夠平滑的問題，消除了其存在的殘余輻射差異，在其鑲嵌拼接線附近，使同種地物的色調(diào)區(qū)域一致，明顯改善了影像質(zhì)量(如圖3所示)；從圖4和圖5可更清楚的看出本發(fā)明的處理效果。
權(quán)利要求
1.一種DMC彩色合成影像的輻射再處理方法，其特征在于，對DMC彩色合成影像中的每一個色彩過渡區(qū)域按如下步驟分別進行處理一、色彩過渡區(qū)域和拼接線的自動定位a、建立影像金字塔在DMC彩色合成影像中，對每兩個相鄰的CCD影像之間的矩形區(qū)域，采用像元平均的方法建立3-5層金字塔影像層；每兩個相鄰的CCD影像之間的矩形區(qū)域，當(dāng)兩個相鄰的CCD影像為左右相鄰時，矩形區(qū)域的高度為CCD影像的高度，寬度取相鄰CCD影像中間線兩側(cè)各30-300像素，當(dāng)兩個CCD影像為上下相鄰時，矩形區(qū)域?qū)挾葹镃CD影像的寬度，高度為相鄰CCD影像中間線兩側(cè)各30-300像素；b、影像分塊將每兩個相鄰的CCD影像之間由上一步確定的矩形區(qū)域沿拼接線方向，將其分成大小相同、互不重疊的1-50個矩形影像塊；c、利用金字塔影像層確定過渡區(qū)域和拼接線的具體位置a)、對每個金字塔影像層，在每個影像塊中，取平行于拼接線的各影像行或者影像列的像素的亮度平均值形成一條一維曲線；b)、對每條一維曲線進行二階差分計算；c)、根據(jù)最上層的金字塔影像層計算色彩過渡區(qū)域和拼接線位置對于最上層的金字塔影像層，將相鄰CCD影像的中間線作為本次計算的拼接線的初始位置，中間線兩側(cè)各1-10個像素寬度范圍為色彩過渡區(qū)域的初始范圍；在每個影像塊中，分別計算當(dāng)前金字塔影像層中色彩過渡區(qū)域兩側(cè)區(qū)域的像素值的平均值；若左側(cè)平均值大于右側(cè)平均值，在該影像塊當(dāng)前金字塔影像層對應(yīng)的一維曲線的二階差分的計算結(jié)果中，在拼接線的初始位置左側(cè)搜索最小值，右側(cè)搜索最大值，若左側(cè)平均值小于右側(cè)平均值，則在該影像塊當(dāng)前金字塔影像層對應(yīng)的一維曲線的二階差分的計算結(jié)果中，在拼接線的初始位置左側(cè)搜索最大值，右側(cè)搜索最小值；d)、利用金字塔影像層，將上一層計算得到的色彩過渡區(qū)域和拼接線位置作為本次計算的色彩過渡區(qū)域和拼接線的初始位置；在每個影像塊中，分別計算當(dāng)前金字塔影像層中色彩過渡區(qū)域兩側(cè)區(qū)域的像素值的平均值；設(shè)色彩過渡區(qū)域的初始范圍在一維曲線中對應(yīng)的區(qū)間為[xA，xB]，xA、xB為對應(yīng)的與拼接線平行的影像列或影像行的編號，若左側(cè)平均值大于右側(cè)平均值，則在該影像塊當(dāng)前金字塔影像層對應(yīng)的一維曲線的二階差分的計算結(jié)果中，在xA左右各1-10個像素范圍內(nèi)搜索最小值，同時保證搜索范圍在初始拼接線的左側(cè)，在xB左右各1-10個像素范圍內(nèi)搜索最大值，同時保證搜索范圍在初始拼接線的右側(cè)，若左側(cè)平均值小于右側(cè)平均值，則在該影像塊當(dāng)前金字塔影像層對應(yīng)的一維曲線的二階差分的計算結(jié)果中，在xA左右各1-10個像素范圍內(nèi)搜索最大值，同時保證搜索范圍在初始拼接線的左側(cè)，在xB左右各1-10個像素范圍內(nèi)搜索最小值，同時保證搜索范圍在初始拼接線的右側(cè)；e)根據(jù)各影像塊計算得到的色彩過渡區(qū)域和拼接線的位置決定相鄰CCD影像間色彩過渡區(qū)域和拼接線的位置對每個影像塊色彩過渡區(qū)域位置進行投票，取出現(xiàn)頻數(shù)最高的色彩過渡區(qū)域和拼接線的位置作為最后的相鄰CCD影像間的色彩過渡區(qū)域[XA，XB]和拼接線的位置(XA+XB)/2，XA，XB為對應(yīng)的出現(xiàn)頻數(shù)最高的與拼接線平行的影像列或影像行的編號；二、“打斷”處理“打斷”處理按下面步驟在紅、綠、藍三個波段上分別進行a、計算相鄰CCD影像之間的灰度差異值將色彩過渡區(qū)域沿拼接線方向上分成大小相等、互不重疊的矩形區(qū)域，這樣拼接線也同時被分成了大小相等、互不重疊的小段，使每小段拼接線長度為10-100像素，然后分別計算每個小段拼接線所在的色彩過渡區(qū)域兩側(cè)3-100像素范圍內(nèi)灰度值的均值CiL、CiR，其中i為該小段拼接線的編號，則第i段拼接線兩側(cè)CCD影像之間存在的灰度差異值為ΔCi＝CiR-CiLb、平滑各段拼接線兩側(cè)CCD影像之間存在的灰度差異值采用一維均值濾波器對求得的灰度差異值進行平滑，濾波器窗口大小取3至51；c、在拼接線兩側(cè)分別改正相鄰CCD影像之間的灰度差異值以垂直于拼接線方向的影像行或影像列為單位，在拼接線兩側(cè)各40-3000像素范圍內(nèi)對每個像素的灰度值按下述公式進行改正 其中Iin表示處理前像素灰度值；Iout表示處理后像素灰度值；d為處理像素到拼接線的距離；w為拼接線一側(cè)的改正寬度；ΔC′為該像素所在的垂直于拼接線方向的影像行或影像列需要改正的灰度差異值，該值通過與該影像行或影像列最鄰近的兩個小段拼接線兩側(cè)CCD影像之間的灰度差異值ΔCi及ΔCi+1線性內(nèi)插得到；三、“接縫”消除處理“接縫”消除處理按下面步驟在紅、綠、藍三個波段上分別進行a、將色彩過渡區(qū)域沿拼接線方向上分成大小相等、互不重疊的矩形區(qū)域，這樣拼接線也同時被分成了大小相等、互不重疊的小段，使每小段拼接線長度為1-20像素，然后分別計算每一小段拼接線兩側(cè)1-10像素范圍內(nèi)灰度值的均值FjL、FjR，其中j為該小段拼接線的編號，則第j段拼接線兩側(cè)的灰度差異值為ΔFj＝FjR-FjL。b、在每小段拼接線所在的色彩過渡區(qū)域內(nèi)，分別計算平行于拼接線的各影像行或影像列的像素灰度值的平均值Mjn，其中n為平行于拼接線的影像行或影像列的編號；c、由下述公式計算各小段拼接線所在的色彩過渡區(qū)域的灰度改正值Δgj=MjXA+ΔFj(n-XA)(XB-XA)-Mjn]]>色彩過渡區(qū)域內(nèi)的每個像素的灰度改正值由與該像素最鄰近的兩個小段拼接線所在的色彩過渡區(qū)域的灰度改正值Δgj和Δgj+1線性內(nèi)插得到；將色彩過渡區(qū)域內(nèi)的每一個像素的灰度值加上該像素的灰度改正值就得到“接縫”消除處理的結(jié)果影像。
全文摘要
本發(fā)明屬于圖像處理與分析領(lǐng)域，針對DMC彩色合成影像由于不同相機的CCD在成像時響應(yīng)的差異性以及后處理等各種因素的影響存在的殘余輻射差異問題，公開了一種DMC彩色合成影像的輻射再處理方法。該方法采用金子塔策略以及差分影像的邊緣檢測原理來自動定位拼接線以及色彩過渡區(qū)域的位置，并采用了分塊方法以及投票機制來提高方法的穩(wěn)健性，確保能可靠的檢測出拼接線以及色彩過渡區(qū)域的準確位置。在此基礎(chǔ)上，該方法采用了“先破后立”的思想，首先打斷原影像灰度值在色彩過渡區(qū)域的連續(xù)性，產(chǎn)生新的接縫，然后再通過消除這種接縫來重構(gòu)圖像，以達到消除DMC彩色合成影像存在的殘余輻射差異問題。
文檔編號H04N1/56GK1924927SQ20061012454
公開日2007年3月7日 申請日期2006年9月18日 優(yōu)先權(quán)日2006年9月18日
發(fā)明者潘俊, 王密 申請人:武漢大學(xué)