本發(fā)明涉及聲像技術(shù)領(lǐng)域，特別是涉及一種遙感圖像相對輻射校正方法及系統(tǒng)。

背景技術(shù)：

隨著遙感技術(shù)飛速發(fā)展，利用多傳感器、多時相遙感數(shù)據(jù)進(jìn)行土地利用和土地覆蓋變化監(jiān)測、全球資源環(huán)境分析、氣候變化監(jiān)測等應(yīng)用越來越廣泛，拓寬了遙感數(shù)據(jù)在時間維度和空間維度的利用，這種大區(qū)域大維度的數(shù)據(jù)應(yīng)用，往往對數(shù)據(jù)集一致性要求非常高，甚至高于對單幅圖像的精度要求。

對于多景數(shù)據(jù)處理，由于每景數(shù)據(jù)的處理誤差、獲取時間的差異以及傳感器差異等問題，造成拼接后的在數(shù)據(jù)拼接線附近出現(xiàn)“邊界”現(xiàn)象，這種“邊界”現(xiàn)象增大了后續(xù)工程應(yīng)用的難度。因此，現(xiàn)有多數(shù)采用相對輻射校正技術(shù)進(jìn)行輻射處理，能最大限度保證圖像間的輻射一致性，不僅能夠糾正大氣狀況變化帶來的差異，而且能夠削弱傳感器等其它原因產(chǎn)生的噪音。

雖然現(xiàn)有的相對輻射校正方法能較好的保證圖像之間輻射一致性，但對圖像時相有嚴(yán)格的約束條件：要求參考圖像與目標(biāo)圖像的成像時間相近，即要求可以是不同年份的，但月份要求在一個月以內(nèi)，成像時刻在1小時以內(nèi)，這一約束大大限制了全國、甚至全球范圍內(nèi)遙感數(shù)據(jù)進(jìn)行土地利用、資源分析等應(yīng)用，特別是針對時相差別較大、傳感器成像幾何差異等原因造成的圖像本身輻射差異較大，很難保證圖像之間輻射一致性，因此大大降低了遙信數(shù)據(jù)的利用率。

技術(shù)實現(xiàn)要素：

本發(fā)明的目的是提供一種遙感圖像相對輻射校正方法及系統(tǒng)，以保證圖像間的輻射一致性，以提高遙信數(shù)據(jù)的利用率。

為實現(xiàn)上述目的，本發(fā)明提供了一種遙感圖像相對輻射校正方法，該包括以下步驟：

獲取目標(biāo)圖像；

根據(jù)所述目標(biāo)圖像確定與所述目標(biāo)圖像相對應(yīng)的參考圖像；

分別對所述參考圖像、所述目標(biāo)圖像進(jìn)行預(yù)處理后，獲得對應(yīng)所述參考圖像的參考圖像像素點集和對應(yīng)所述目標(biāo)圖像的目標(biāo)圖像像素點集；

基于核典型相關(guān)分析方法根據(jù)所述參考圖像像素點集和所述目標(biāo)圖像像素點集確定非線性輻射控制點；

根據(jù)所述非線性輻射控制點對所述目標(biāo)圖像各階段像素點進(jìn)行校正，獲得地表反射率值。

可選的，所述根據(jù)所述目標(biāo)圖像確定與所述目標(biāo)圖像相對應(yīng)的參考圖像的具體步驟包括：根據(jù)與所述目標(biāo)圖像具有相同重疊區(qū)域、近似成像幾何、坐標(biāo)范圍的圖像確定為與所述目標(biāo)圖像相對應(yīng)的參考圖像。

可選的，所述基于核典型相關(guān)分析方法根據(jù)所述參考圖像像素點集和所述目標(biāo)圖像像素點集確定非線性輻射控制點的具體步驟包括：

分別將所述參考圖像像素點集、所述目標(biāo)圖像像素點集映射到核空間，利用核典型相關(guān)分析方法確定參考圖像核典型變量、目標(biāo)圖像核典型變量；

基于mad中的權(quán)重公式對所述參考圖像核典型變量、所述目標(biāo)圖像核典型變量進(jìn)行篩選，確定非線性輻射控制點。

可選的，所述權(quán)重公式為：

其中：u為映射到高維特征空間的目標(biāo)圖像核典型變量集合，ui為u的第i個元素，v為映射到高維特征空間的參考圖像核典型變量集合，vi是v的第i個元素，var()為方差，z為標(biāo)準(zhǔn)化后的方差，ρλi為第i個特征根，為卡方檢驗，pr(nochange)為不變像素點的概率值，n為重疊區(qū)域的點的個數(shù)。

可選的，所述根據(jù)所述非線性輻射控制點對所述目標(biāo)圖像各階段像素點進(jìn)行校正，獲得地表反射率值的具體步驟包括：

對所述非線性輻射控制點進(jìn)行非線性關(guān)系擬合獲得非線性擬合公式；

根據(jù)所述非線性擬合公式對所述目標(biāo)圖像各波段像素點進(jìn)行校正，獲得地表反射率。

可選的，所述對所述非線性輻射控制點進(jìn)行非線性關(guān)系擬合獲得非線性擬合公式具體步驟包括：

根據(jù)偏差平方和最小公式對所述非線性輻射控制點進(jìn)行非線性關(guān)系擬合獲得非線性擬合公式；

所述偏差平方和最小公式為：

其中，δi為第i個非線性輻射控制點的偏差，ri為第i個非線性輻射控制點的原始目標(biāo)圖像像素值，f(ri)為將第i個非線性輻射控制點的原始目標(biāo)圖像像素值代入擬合的非線性關(guān)系表達(dá)式后，計算得到的像素值，si為第i個非線性輻射控制點對應(yīng)的參考圖像像素值，r為非線性輻射控制點的總個數(shù)。

本發(fā)明還提供了一種遙感圖像輻射校正系統(tǒng)，該系統(tǒng)包括：

獲取模塊，用于獲取目標(biāo)圖像；

匹配模塊，用于根據(jù)所述目標(biāo)圖像確定與所述目標(biāo)圖像相對應(yīng)的參考圖像；

處理模塊，用于分別對所述參考圖像、所述目標(biāo)圖像進(jìn)行預(yù)處理后，獲得對應(yīng)所述參考圖像的參考圖像像素點集和對應(yīng)所述目標(biāo)圖像的目標(biāo)圖像像素點集；

非線性輻射控制點確定模塊，用于基于核典型相關(guān)分析方法根據(jù)所述參考圖像像素點集和所述目標(biāo)圖像像素點集確定非線性輻射控制點；

校正模塊，用于根據(jù)所述非線性輻射控制點對所述目標(biāo)圖像各階段像素點進(jìn)行校正，獲得地表反射率值。

可選的，所述非線性輻射控制點確定模塊的具體包括：

第一確定單元，用于分別將所述參考圖像像素點集、所述目標(biāo)圖像像素點集映射到核空間，利用核典型相關(guān)分析方法確定參考圖像核典型變量、目標(biāo)圖像核典型變量；

第二確定單元，用于基于mad中的權(quán)重公式對所述參考圖像核典型變量、所述目標(biāo)圖像核典型變量進(jìn)行篩選，確定非線性輻射控制點。

可選的，所述權(quán)重公式為：

其中：u為映射到高維特征空間的目標(biāo)圖像核典型變量集合，ui為u的第i個元素，v為映射到高維特征空間的參考圖像核典型變量集合，vi是v的第i個元素，var()為方差，z為標(biāo)準(zhǔn)化后的方差，ρλi為第i個特征根，為卡方檢驗，pr(nochange)為不變像素點的概率值，n為重疊區(qū)域的點的個數(shù)。

可選的，所述校正模塊具體包括：

擬合單元，用于對所述非線性輻射控制點進(jìn)行非線性關(guān)系擬合獲得非線性擬合公式；

校正單元，用于根據(jù)所述非線性擬合公式對所述目標(biāo)圖像各波段像素點進(jìn)行校正，獲得地表反射率。

根據(jù)本發(fā)明提供的具體實施例，本發(fā)明公開了以下技術(shù)效果：

1)、本發(fā)明基于核典型相關(guān)分析方法自動確定非線性輻射控制點，保證目標(biāo)圖像和參考圖像間光譜一致，消除圖像之間的輻射差異，提高圖像間的輻射一致性；特別是針對時相差別較大、傳感器成像幾何差異等原因造成的圖像本身輻射差異較大，也能保證圖像之間輻射一致性，大大提高了遙信數(shù)據(jù)的利用率。

2)、本發(fā)明采用核典型相關(guān)分析的辦法自動確定非線性輻射控制點，避免了選擇參考圖像時對時相的約束和受傳感器的限制、通用性強(qiáng)、無需人工交互、處理流程單、運算速度快、算法自動性高、適用性強(qiáng)等特點。

附圖說明

為了更清楚地說明本發(fā)明實施例或現(xiàn)有技術(shù)中的技術(shù)方案，下面將對實施例中所需要使用的附圖作簡單地介紹，顯而易見地，下面描述中的附圖僅僅是本發(fā)明的一些實施例，對于本領(lǐng)域普通技術(shù)人員來講，在不付出創(chuàng)造性勞動性的前提下，還可以根據(jù)這些附圖獲得其他的附圖。

圖1為本發(fā)明實施例遙感圖像相對輻射校正方法的具體流程圖；

圖2為本發(fā)明實施例遙感圖像相對輻射校正系統(tǒng)的結(jié)構(gòu)示意圖。

具體實施方式

下面將結(jié)合本發(fā)明實施例中的附圖，對本發(fā)明實施例中的技術(shù)方案進(jìn)行清楚、完整地描述，顯然，所描述的實施例僅僅是本發(fā)明一部分實施例，而不是全部的實施例?；诒景l(fā)明中的實施例，本領(lǐng)域普通技術(shù)人員在沒有做出創(chuàng)造性勞動前提下所獲得的所有其他實施例，都屬于本發(fā)明保護(hù)的范圍。

本發(fā)明的目的是提供一種遙感圖像相對輻射校正方法及系統(tǒng)，以保證圖像間的輻射一致性，以提高遙信數(shù)據(jù)的利用率。

為使本發(fā)明的上述目的、特征和優(yōu)點能夠更加明顯易懂，下面結(jié)合附圖和具體實施方式對本發(fā)明作進(jìn)一步詳細(xì)的說明。

圖1為本發(fā)明實施例遙感圖像相對輻射校正方法的具體流程圖；詳見圖1。其具體的步驟包括：

步驟100：根據(jù)所述目標(biāo)圖像確定與所述目標(biāo)圖像相對應(yīng)的參考圖像。

根據(jù)與所述目標(biāo)圖像具有相同重疊區(qū)域、近似成像幾何、坐標(biāo)范圍的圖像確定為與所述目標(biāo)圖像相對應(yīng)的參考圖像，所述參考圖像與所述目標(biāo)圖像無需時相相近，大大提高輻射處理圖像間的輻射一致性，提高了遙信數(shù)據(jù)的利用率。

步驟200：分別對所述參考圖像、所述目標(biāo)圖像進(jìn)行預(yù)處理后，分別獲得參考圖像像素點集和目標(biāo)圖像像素點集。

所述預(yù)處理包括坐標(biāo)重投影、幾何配準(zhǔn)、計算地理重疊區(qū)域、各波段像素點對應(yīng)、并使用掩膜技術(shù)排除極值點、過飽和點、云、云下陰影及水體的像素點。對所述參考圖像、所述目標(biāo)圖像進(jìn)一步處理，為了獲得較為精確的參考圖像像素點集和目標(biāo)圖像像素點集。

坐標(biāo)重投影與幾何配準(zhǔn)屬通常的遙感圖像處理步驟。

計算地理重疊區(qū)域：對于單個目標(biāo)圖像，有可能找到多景參考圖像，此時就需要進(jìn)行重疊區(qū)域內(nèi)數(shù)據(jù)的裁切與拼接整合。

比各波段像素點對應(yīng)：對于目標(biāo)圖像與參考圖像為同種傳感器數(shù)據(jù)，分辨率相同，當(dāng)坐標(biāo)重投影、幾何配準(zhǔn)精度足夠高時，像素點本身就是對應(yīng)的。對于不同傳感器的遙感圖像，分辨率差異超過3倍以上時，就需要插值。

掩膜過濾異常點：該步先通過統(tǒng)計目標(biāo)圖像與參考圖像重疊區(qū)內(nèi)的直方圖，過濾掉過亮和過暗的像素點，然后檢測參考圖像中的云、云下陰影、水體并將其掩膜掉。具體的檢測算法依據(jù)不同的傳感器數(shù)據(jù)而不同。

對掩膜過濾后得到的數(shù)據(jù)點集進(jìn)行典型相關(guān)分析。

對掩膜過濾后得到的數(shù)據(jù)點集進(jìn)行典型相關(guān)分析，進(jìn)一步從這些點集中選出典型相關(guān)點集。所選用的典型相關(guān)分析算法應(yīng)當(dāng)具有不受總體的大氣狀況和傳感器標(biāo)定導(dǎo)致的線性變換的影響的特性。

典型相關(guān)分析算法為了遮蔽兩個不同時相圖像中的變化像素，首先形成兩幅圖像n個波段內(nèi)像素值的線性組合。

步驟300：基于核典型相關(guān)分析方法根據(jù)所述參考圖像像素點集和所述目標(biāo)圖像像素點集確定非線性輻射控制點。

步驟301：分別將所述參考圖像像素點集、所述目標(biāo)圖像像素點集映射到核空間，利用核典型相關(guān)分析方法確定參考圖像核典型變量、目標(biāo)圖像核典型變量。

核典型相關(guān)分析方法是一種常用的非線性典型相關(guān)方法，將核函數(shù)的思想引入典型相關(guān)分析當(dāng)中，把低維數(shù)據(jù)映射到高維的特征空間中，為在高維特征中解決復(fù)雜關(guān)聯(lián)分析奠定基礎(chǔ)。

常用的核函數(shù)種類非常多，本次實施中選取的核函數(shù)為anova核函數(shù)該核函數(shù)屬于徑向基的一種，非常適合多維回歸的問題。其具體規(guī)則為：設(shè)x,z∈x，x屬于rⁿ空間，非線性函數(shù)實現(xiàn)輸入空間x到特征空間f的映射，其中f屬于r^m，n<<m。則核函數(shù)公式為：

其中<,>為內(nèi)積，k(x,z)即為核函數(shù)。

x為目標(biāo)圖像像素點集，用目標(biāo)圖像像素點矩陣表示，此處作為樣本矩陣輸入，n為重疊區(qū)域的點的個數(shù)，n為每個像素點的數(shù)據(jù)維度，即為各波段的灰度值，則x可為n個n維向量，xn×n＝(x1,x2...,xn)，xi為第i個n維向量，y為參考圖像像素點集，用參考圖像像素點矩陣表示，yi為第i個n維向量，此處作為樣本矩陣輸入，n為重疊區(qū)域的點的個數(shù)，n為每個像素點的數(shù)據(jù)維度，即為各波段的灰度值，y可為n個n維向量，yn×n＝(y1,y2...,yn)。

依據(jù)核典型相關(guān)分析方法，將經(jīng)過預(yù)處理的目標(biāo)圖像和參考圖像的n個波段內(nèi)像素值分別映射為高維空間的變量，獲取其在高維特征空間的線性組合，具體公式為：

其中，c和d分別為高維空間(m維)的常向量矩陣，c^t為c的轉(zhuǎn)置，d^t為d的轉(zhuǎn)置，將目標(biāo)圖像的向量x由低維空間(rⁿ)映射到高維空間(r^m)的非線性變換函數(shù)(n<<m)，為將參考圖像的向量y由低維空間(rⁿ)映射到高維空間(r^m)的非線性變換函數(shù)，變換后的和均為m×n維矩陣，u為映射到高維特征空間的目標(biāo)圖像核典型變量集合，v為映射到高維特征空間的參考圖像核典型變量集合。

常向量矩陣c和d可以用n維空間的常向量矩陣α和β與和為：

其中，α＝(α1,α2,,...,αn)，αi是α的第i個元素，每一個αi是一個n維向量；同理βi是β的第i個元素。此時將求解高維空間常向量c和d的問題轉(zhuǎn)化為求解低維空間常向量α和β的問題。將上述公式(5)代入公式(4)，則可得到綜合變量集合u和v為：

步驟302：基于mad中的權(quán)重公式對所述參考圖像核典型變量、所述目標(biāo)圖像核典型變量進(jìn)行篩選，確定非線性輻射控制點，其中、mad為多元變化檢測算法(multivariatealterationdetection)。

根據(jù)核函數(shù)公式(3)確定在x和y上定義核矩陣kx和ky為：

其中，kx(i,j)表示核矩陣kx的第i行j列元素，ky(i,j)表示核矩陣ky的第i行j列元素。

通過綜合變量u和v變換公式(6)和核矩陣公式(7)，我們可以計算出u和v的方差和協(xié)方差為公式為：

其中，var()為方差，cov()為協(xié)方差。

結(jié)合上述公式(8)，可得到相關(guān)系數(shù)的表達(dá)公式(9)：

其中，ρ為pearson相關(guān)系數(shù)。

令相關(guān)系數(shù)ρ取最大值時，求解時固定分母來求分子，即該問題簡化轉(zhuǎn)換為一個優(yōu)化問題，即優(yōu)化公式為：

這個問題為有等式約束條件的優(yōu)化問題，可以采用拉格朗日乘數(shù)法(lagrangemultiplier)求解。因此根據(jù)公式(10)構(gòu)建的拉格朗日公式為：

其中，ρα和ρβ都是待求特征根。

分別對拉格朗日公式中的α、β求偏導(dǎo)，且令偏導(dǎo)分別等于0，具體詳見公式(12)：

令ρλ＝ρα＝ρβ，同時防止矩陣不可逆對其正則化，上述問題可轉(zhuǎn)換為求解公式(12)的特征根問題：

其中，ε1和ε2是正則化參數(shù)，取值范圍[0，1]。

根據(jù)公式(13)求解不同的ρλ以及對應(yīng)的典型向量α和β。

本發(fā)明依據(jù)nielsen和canty在低維空間中提出的概念，在var(u)＝var(v)＝1條件下，使var(u-v)最大的情況下，ui-vi中最后幾個變量集中包含著兩組數(shù)據(jù)之間的線性不變關(guān)系，將上述規(guī)則映射到高位空間進(jìn)行求解，以選取圖像間的非線性不變點，采用迭代的辦法選取不變信息，因為mad變量是測量變量的線性組合，由中心極限定理可得它們呈現(xiàn)的高斯分布。每一個獨立的mad變量互相正交，標(biāo)準(zhǔn)化后的像素方差符合自由度為n的分布，其具體權(quán)重公式為：

其中，z為標(biāo)準(zhǔn)化后的方差，ui為u中的第i個元素，vi是v中的第i個元素，ρλi為第i個特征根，為卡方檢驗，pr(nochange)為不變像素點的概率值，用來確定哪些像素可以最終被選為不變特征的，可以設(shè)定一個固定的域值t(比如t＝95％)，當(dāng)pr(nochange)＞t時，判斷為特征不變像素點，以此獲取非線性輻射控制點。

步驟400：所述根據(jù)所述非線性輻射控制點對所述目標(biāo)圖像各階段像素點進(jìn)行校正，獲得地表反射率值。

步驟401：采用最小二乘法對所述非線性輻射控制點進(jìn)行非線性關(guān)系擬合獲得非線性擬合公式。

擬合公式為：

f(r)＝a0+a1r+...+apr^p(14)；

其中，f(r)為待擬合的表達(dá)式，將輻射控制點的原始像素值輸入該式可獲取校正后的像素值，a0、a1...ap為待求解的常系數(shù)，r為非線性輻射控制點的原始像素值。

偏差平方最小公式為：

其中，δi為第i個非線性輻射控制點的偏差，ri為第i個非線性輻射控制點的原始目標(biāo)圖像像素值，f(ri)為將第i個非線性輻射控制點的原始目標(biāo)圖像像素值代入擬合的非線性關(guān)系表達(dá)式后，計算得到的像素值，si為第i個非線性輻射控制點的參考圖像像素值，r為非線性輻射控制點的總個數(shù)。

令p為一個常數(shù)，將第i個非線性輻射控制點的原始目標(biāo)圖像像素值ri和第i個非線性輻射控制點對應(yīng)的參考圖像目標(biāo)像素值si帶入擬合公式(14)便可確定待求解的常系數(shù)a0、a1...ap，進(jìn)而求解出具體擬合公式。

在實驗中取不同p值，可以得到不同的非線性關(guān)系表達(dá)式，依據(jù)校正結(jié)果進(jìn)行調(diào)整，選取校正結(jié)果最好時的p值。

將第i個非線性輻射控制點的原始目標(biāo)圖像像素值ri帶入多個不同的擬合公式便可得出多個不同的f(ri)，再將多個不同的f(ri)和多個第i個非線性輻射控制點對應(yīng)的參考圖像目標(biāo)像素值si帶入偏差平方最小公式2便可求取所有非線性輻射控制點總體最小偏差。

根據(jù)最小偏差便可確定該擬合公式為非線性擬合公式。

步驟402：根據(jù)所述非線性擬合公式對所述目標(biāo)圖像各波段像素點進(jìn)行校正，獲得地表反射率，后續(xù)可作為地表反射率產(chǎn)品的應(yīng)用。

根據(jù)上述方法獲取的非線性關(guān)系表達(dá)式，對所述目標(biāo)圖像的各波段像素點帶入非線性擬合公式便可求出校正后的像素點，所述校正后的像素點即為地表反射率。

本發(fā)明還提供了一種遙感圖像相對輻射校正系統(tǒng)，以保證圖像間的輻射一致性，以提高遙信數(shù)據(jù)的利用率。

圖2為本發(fā)明實施例遙感圖像相對輻射校正系統(tǒng)的結(jié)構(gòu)示意圖，如圖2所示，所述系統(tǒng)包括：獲取模塊1、匹配模塊2、處理模塊3、非線性輻射控制點確定模塊4、校正模塊5。

獲取模塊1，用于獲取目標(biāo)圖像。

匹配模塊2，用于根據(jù)所述目標(biāo)圖像確定與所述目標(biāo)圖像相對應(yīng)的參考圖像。

根據(jù)與所述目標(biāo)圖像具有相同重疊區(qū)域、近似成像幾何、坐標(biāo)范圍的圖像確定為與所述目標(biāo)圖像相對應(yīng)的參考圖像，所述參考圖像與所述目標(biāo)圖像無需時相相近，大大提高輻射處理圖像間的輻射一致性，提高了遙信數(shù)據(jù)的利用率。

處理模塊3，用于分別對所述參考圖像、所述目標(biāo)圖像進(jìn)行預(yù)處理后，獲得對應(yīng)所述參考圖像的參考圖像像素點集和對應(yīng)所述目標(biāo)圖像的目標(biāo)圖像像素點集。

所述預(yù)處理包括坐標(biāo)重投影、幾何配準(zhǔn)、計算地理重疊區(qū)域、各波段像素點對應(yīng)、并使用掩膜技術(shù)排除極值點、過飽和點、云、云下陰影及水體的像素點。對所述參考圖像、所述目標(biāo)圖像進(jìn)一步處理，為了獲得較為精確的參考圖像像素點集和目標(biāo)圖像像素點集。

非線性輻射控制點確定模塊4，用于基于核典型相關(guān)分析方法根據(jù)所述參考圖像像素點集和所述目標(biāo)圖像像素點集確定非線性輻射控制點。

所述非線性輻射控制點確定模塊4，具體包括：第一確定單元、第二確定單元。

第一確定單元分別將所述參考圖像像素點集、所述目標(biāo)圖像像素點集映射到核空間，利用核典型相關(guān)分析方法確定參考圖像核典型變量、目標(biāo)圖像核典型變量；

第二確定單元基于mad中的權(quán)重公式對所述參考圖像核典型變量、所述目標(biāo)圖像核典型變量進(jìn)行篩選，確定非線性輻射控制點。

所述權(quán)重公式為：

其中：u為映射到高維特征空間的目標(biāo)圖像核典型變量集合，ui為u的第i個元素，v為映射到高維特征空間的參考圖像核典型變量集合，vi是v的第i個元素，var()為方差，z為標(biāo)準(zhǔn)化后的方差，ρλi為第i個特征根，為卡方檢驗，pr(nochange)為不變像素點的概率值，n為重疊區(qū)域的點的個數(shù)。

校正模塊5，用于根據(jù)所述非線性輻射控制點對所述目標(biāo)圖像各階段像素點進(jìn)行校正，獲得地表反射率值；所述校正模塊具體包括：擬合單元、校正單元。

擬合單元根據(jù)偏差平方和最小公式對所述非線性輻射控制點進(jìn)行非線性關(guān)系擬合獲得非線性擬合公式；

所述偏差平方和最小公式為：

其中，δi為第i個非線性輻射控制點的偏差，ri為第i個非線性輻射控制點的原始目標(biāo)圖像像素值，f(ri)為將第i個非線性輻射控制點的原始目標(biāo)圖像像素值代入擬合的非線性關(guān)系表達(dá)式后，計算得到的像素值，si為第i個非線性輻射控制點對應(yīng)的參考圖像像素值，r為非線性輻射控制點的總個數(shù)。

校正單元根據(jù)所述非線性擬合公式對所述目標(biāo)圖像各波段像素點進(jìn)行校正，獲得地表反射率。

本說明書中各個實施例采用遞進(jìn)的方式描述，每個實施例重點說明的都是與其他實施例的不同之處，各個實施例之間相同相似部分互相參見即可。對于實施例公開的系統(tǒng)而言，由于其與實施例公開的方法相對應(yīng)，所以描述的比較簡單，相關(guān)之處參見方法部分說明即可。

本文中應(yīng)用了具體個例對本發(fā)明的規(guī)則及實施方式進(jìn)行了闡述，以上實施例的說明只是用于幫助理解本發(fā)明的方法及其核心思想；同時，對于本領(lǐng)域的一般技術(shù)人員，依據(jù)本發(fā)明的思想，在具體實施方式及應(yīng)用范圍上均會有改變之處。綜上所述，本說明書內(nèi)容不應(yīng)理解為對本發(fā)明的限制。