專利名稱：：一種基于查找表的遙感影像逐像元大氣校正方法
技術(shù)領(lǐng)域：
：本發(fā)明涉及一種遙感影像的大氣校正方法，具體涉及一種基于查找表的遙感影像逐像元大氣校正方法。
背景技術(shù)：
：隨著遙感理論和技術(shù)的不斷發(fā)展及其在應(yīng)用中的不斷完善，遙感已經(jīng)從開始的定性解釋發(fā)展到了現(xiàn)在的定量分析。大多數(shù)地表參數(shù)的定量遙感模型如反照率、葉面積指數(shù)等等都是建立在地表反射率的基礎(chǔ)上，通過(guò)大氣校正獲取精確的地表反射率信息是遙感定量分析的關(guān)鍵問(wèn)題，也是主要難點(diǎn)之一。隨著定量遙感技術(shù)迅速發(fā)展，特別是利用多傳感器、多時(shí)相遙感數(shù)據(jù)進(jìn)行土地利用和土地覆蓋變化監(jiān)測(cè)、全球資源環(huán)境分析、氣候變化監(jiān)測(cè)等的需要，使得遙感圖像大氣校正方法的研究越來(lái)越受到重視，對(duì)遙感圖像大氣校正的研究具有重大意義。太陽(yáng)輻射從穿過(guò)大氣層到達(dá)地表然后被反射再次穿過(guò)大氣層到達(dá)遙感傳感器的過(guò)程中，受到大氣分子、氣溶膠和云粒子等大氣成份吸收與散射的影響，使其獲取的遙感信息中帶有一定的非目標(biāo)地物的成像信息，因此需要進(jìn)行大氣校正去除遙感圖像中的大氣影響。大氣對(duì)太陽(yáng)輻射的影響主要來(lái)源于大氣的吸收、大氣分子的瑞利散射和氣溶膠的米散射作用，大氣校正的目標(biāo)就在于從遙感傳感器觀測(cè)到的輻射信號(hào)中去除這些作用的影響，得到地面目標(biāo)物的真實(shí)反射率。目前大氣校正的方法有暗目標(biāo)像元法、地表實(shí)測(cè)線性回歸法、大氣輻射傳輸模型法等等。其中輻射傳輸模型法由于普適性好、精度較高在近些年來(lái)得到廣泛應(yīng)用。輻射傳輸模型法的基本原理是通過(guò)大氣狀況、衛(wèi)星和太陽(yáng)位置信息等輸入?yún)?shù)，利用合適的大氣輻射傳輸模型計(jì)算出大氣校正參數(shù)的數(shù)值解(如大氣透過(guò)率、程輻射等)，以此為依據(jù)計(jì)算地表真實(shí)反射率。目前已有多位學(xué)者利用此方法對(duì)遙感影像進(jìn)行了大氣校正處理。由于前人的研究主要針對(duì)范圍較小的TM、ASTER、CBERS等遙感影像，考慮的變化因素相對(duì)較少，或者把整幅影像的輸入?yún)?shù)看成均一，或者僅僅考慮氣溶膠光學(xué)厚度的變化。也有學(xué)者利用大氣輻射傳輸模型對(duì)MODIS影像進(jìn)行了大氣校正，但是所選研究區(qū)范圍都不大，并非針對(duì)全景影像，也都是采用了統(tǒng)一的輸入?yún)?shù)進(jìn)行大氣校正。但是對(duì)于大范圍的遙感影像而言，各個(gè)像元的差異比較大，用統(tǒng)一的輸入?yún)?shù)進(jìn)行校正無(wú)法考慮各個(gè)像元大氣狀況、觀測(cè)角度等條件的差異，但是如果每個(gè)像元都在線運(yùn)行輻射傳輸模型計(jì)算將耗費(fèi)大量的計(jì)算機(jī)資源和時(shí)間，如何在保證精度的條件下快速實(shí)現(xiàn)遙感影像的逐像元大氣校正成為迫切需要解決的問(wèn)題。
發(fā)明內(nèi)容本發(fā)明目的是針對(duì)現(xiàn)有技術(shù)存在的缺陷，在分析大氣校正結(jié)果對(duì)于大氣氣溶膠光學(xué)厚度、衛(wèi)星和太陽(yáng)位置、海拔高度輸入?yún)?shù)變化敏感性的基礎(chǔ)上，通過(guò)6S大氣輻射傳輸模型離線計(jì)算建立了不同條件下大氣校正系數(shù)的查找表，基于查找表對(duì)MODIS影像進(jìn)行逐像元的大氣校正。本發(fā)明的校正結(jié)果與利用統(tǒng)一輸入?yún)?shù)校正結(jié)果進(jìn)行比較分析表明，本發(fā)明的方法具有更高的精度和實(shí)用性。本發(fā)明為實(shí)現(xiàn)上述目的，采用如下技術(shù)方案本發(fā)明一種基于查找表的遙感影像逐像元大氣校正方法包括如下步驟(一)離線模擬(1)采用6S模型進(jìn)行大氣校正，模擬太陽(yáng)輻射在大氣中的傳輸狀況得到地表真實(shí)反射率Ps’并給出大氣校正系數(shù)xa、xb和XC,=_^__pr+a(^sA)ps~tgescos{es)t{es)t{ev)t{0s)t{0v)<{，其中，Tg為大氣中03、H2O,N2,CO2,02、CH4氣體對(duì)太陽(yáng)輻射的吸收透過(guò)率，PE+A為大氣層輻射反射率即氣溶膠程輻射與瑞利散射程輻射之和，θ,、θ別為太陽(yáng)天頂角和傳感器天頂角，τ(θs)、τ(θν)分別為太陽(yáng)至地面、地面至傳感器的大氣路徑透過(guò)率，S為大氣底層向下的半球反射率，L為表觀輻亮度，Es為對(duì)應(yīng)波長(zhǎng)的大氣層頂太陽(yáng)輻照度；(2)敏感性分析遙感影像，只考慮氣溶膠光學(xué)厚度以及傳感器天頂角的變化，太陽(yáng)天頂角及地表海拔設(shè)為常量；(3)利用6S模型根據(jù)步驟⑴所述的大氣校正系數(shù)xa、xb、xc公式離線計(jì)算不同氣溶膠光學(xué)厚度、太陽(yáng)天頂角、傳感器天頂角以及地表海拔情況下的大氣校正系數(shù)xa、xb、xc，并建立針對(duì)氣溶膠光學(xué)厚度、太陽(yáng)天頂角、傳感器天頂角以及地表海拔相對(duì)于大氣校正系數(shù)xa、xb、xc的查找表；(二)在線校正(4)對(duì)遙感影像進(jìn)行大氣校正時(shí)，從遙感影像數(shù)據(jù)中導(dǎo)出表觀輻亮度、太陽(yáng)天頂角、傳感器天頂角數(shù)據(jù)，從遙感影像的氣溶膠光學(xué)厚度中讀出光學(xué)厚度數(shù)據(jù)，結(jié)合Ikm分辨率的SRTMDEM數(shù)據(jù)，逐像元根據(jù)步驟(3)所述的查找表進(jìn)行線性插值，計(jì)算出當(dāng)前像元對(duì)應(yīng)波段的大氣校正系數(shù)xa、xb、xc。本發(fā)明具有如下優(yōu)點(diǎn)(1)本發(fā)明依據(jù)MODIS圖像自身的信息，不需要地表實(shí)測(cè)數(shù)據(jù)，對(duì)于輸入?yún)?shù)的依賴性較小。(2)本發(fā)明的校正方法可以有效改善由于大氣條件、傳感器位置等空間分布差異對(duì)MODIS影像大氣校正的影響。(3)本發(fā)明相對(duì)于6S輻射傳輸模型在線校正法而言，精度相近，速度要快捷很多；相對(duì)于統(tǒng)一輸入?yún)?shù)法而言，精度有明顯提高，能夠體現(xiàn)不同像元大氣性質(zhì)、觀測(cè)條件和地表海拔的空間分布差異。圖1輸入?yún)?shù)變化對(duì)地表反射率校正結(jié)果的影響；圖2基于查找表的MODIS逐像元大氣校正流程圖；圖3掃描線上各像元的氣溶膠光學(xué)厚度、太陽(yáng)天頂角、傳感器天頂角、地表海拔；圖4=MODIS第14波段三種不同方法校正結(jié)果的比較；具體實(shí)施例方式下面結(jié)合附圖對(duì)本發(fā)明進(jìn)一步詳述本發(fā)明一種基于查找表的遙感影像逐像元大氣校正方法，包括如下步驟1)對(duì)輻射傳輸方程進(jìn)行分析推導(dǎo)6S(SecondSimulationoftheSatelliteSignalintheSolarSpectrumradiativecode)模型是目前世界上發(fā)展得比較完善的大氣校正模型之一，由Vermote等人在5S(theSimulationoftheSatelliteSignalintheSolarSpectrumradiativecode)模型的基礎(chǔ)上改進(jìn)而來(lái)，適用于0.254μπι波長(zhǎng)范圍內(nèi)電磁波的大氣輻射傳輸模擬。利用6S模型進(jìn)行大氣校正的工作流程是將大氣參數(shù)、幾何參數(shù)、觀測(cè)波段和海拔高度等條件輸入6S模型之后，通過(guò)計(jì)算模擬太陽(yáng)輻射在大氣中的傳輸狀況，并給出大氣校正系數(shù)xa、xb和xc。根據(jù)下面的公式就可以計(jì)算得到經(jīng)過(guò)大氣校正的地表真實(shí)反射率。y=xa^(measuredradiance)-xb,acr=y/(1.+xc^y)(1)公式1雖然形式簡(jiǎn)單，但是并不是一個(gè)單純的經(jīng)驗(yàn)公式。通過(guò)對(duì)輻射傳輸方程進(jìn)行推導(dǎo)和變換可知，系數(shù)xa、xb、xc和方程是有著具體物理意義的。在6S模型中，表觀反射率可以描述為(2)式中，ρΤ0Α為表觀反射率，Tg為大氣中03、H2O等氣體對(duì)太陽(yáng)輻射的吸收透過(guò)率，PΕ+Α為大氣程輻射反射率(氣溶膠程輻射與瑞利散射程輻射之和)，θs、θν分別為太陽(yáng)天頂角和衛(wèi)星天頂角，Τ(θs)、Τ(θν)分別為太陽(yáng)-地面、地面-傳感器的大氣路徑透過(guò)率，S為大氣底層向下的半球反射率，Ps為地表真實(shí)反射率。上式可改寫為其中表觀反射率PTQA=πL/Escos(θs)，式中，L為表觀輻亮度，Es為對(duì)應(yīng)波長(zhǎng)的Ptoa、H、_D(00\大氣層頂太陽(yáng)輻照度。另外，再設(shè)則上式可改寫為⑷將公式1和4對(duì)照，發(fā)現(xiàn)憂)7^v);xb=r^)m);xc=s。xa、Xb,XC三個(gè)系數(shù)與大氣狀況、太陽(yáng)與傳感器的位置、地面的海拔、觀測(cè)波段等有關(guān)，而與地表反射特性無(wú)關(guān)。在其它條件相同的情況下改變目標(biāo)的表觀反射率或者輻亮度，xa、xb、xc三個(gè)校正系數(shù)不發(fā)生改變。2)敏感性分析通常在利用6S模型進(jìn)行大氣校正時(shí)，往往假設(shè)整個(gè)研究區(qū)大氣條件均一、太陽(yáng)和傳感器的觀測(cè)角度一致、地表具有相同的海拔高度，然后利用6S模型計(jì)算出大氣校正系數(shù)，并應(yīng)用于整個(gè)圖像，得到整個(gè)研究區(qū)的地表真實(shí)反射率。對(duì)于小范圍的遙感影像(如TM或者裁切為小塊的MODIS數(shù)據(jù))而言，這種方法是可行的，但是對(duì)于大范圍的遙感影像(如整幅MODIS或者AVHRR數(shù)據(jù))而言，這樣的簡(jiǎn)單假設(shè)就顯得不合理了，勢(shì)必會(huì)影響校正得到的地表反射率精度。以MODIS數(shù)據(jù)為例，一景影像范圍超過(guò)400萬(wàn)平方公里，各個(gè)像元的大氣條件、太陽(yáng)位置、傳感器位置、海拔高度等均不相同，甚至有較大差異。本發(fā)明首先分析了6S模型地表反射率校正結(jié)果對(duì)于這些參數(shù)的敏感性，從而更深層次的了解各種參數(shù)對(duì)于大氣校正結(jié)果的影響程度，為后面在大氣校正過(guò)程中的參數(shù)選取提供參考。本發(fā)明選擇了大氣氣溶膠光學(xué)厚度、太陽(yáng)天頂角、傳感器天頂角以及地表海拔這4個(gè)參數(shù)來(lái)分析這些參數(shù)的變化對(duì)6S模型大氣校正結(jié)果的影響。6S模型的初始輸入?yún)?shù)為太陽(yáng)天頂角設(shè)為30度，方位角為0度，衛(wèi)星天頂角設(shè)為30度，方位角為90度，成像日期7月1日，大氣模式為中緯度夏季，氣溶膠類型選擇大陸型氣溶膠，550nm光學(xué)厚度為0.5，地表海拔設(shè)為1km，波段選用MODIS的第1、2、3、4波段，表觀反射率設(shè)為0.3。然后分別計(jì)算氣溶膠光學(xué)厚度在01之間變化、太陽(yáng)天頂角在535度之間變化、傳感器天頂角在065度之間變化、海拔高度在04000m之間變化對(duì)于校正得到的地表反射率的影響，結(jié)果如圖1所示。從圖1中可以看出，氣溶膠光學(xué)厚度對(duì)于大氣校正的影響明顯要高于其它因素，傳感器天頂角對(duì)于大氣校正的影響也比較大，太陽(yáng)天頂角以及海拔高度的變化對(duì)于6S模型大氣校正的影響相對(duì)較小。如果對(duì)于大氣校正精度要求不高的話，可以只考慮氣溶膠光學(xué)厚度以及傳感器天頂角的變化，其它要素設(shè)為常量，這樣查找表比較簡(jiǎn)單而且運(yùn)算速度大大提高。3)建立查找表利用6S模型離線計(jì)算不同氣溶膠光學(xué)厚度、太陽(yáng)天頂角、傳感器天頂角以及地表海拔情況下的校正系數(shù)xa、xb、xc，建立針對(duì)這4個(gè)因素的查找表。6S模型的初始參數(shù)設(shè)置基本同前面敏感性分析時(shí)的輸入?yún)?shù)，不過(guò)氣溶膠光學(xué)厚度、太陽(yáng)天頂角、傳感器天頂角和海拔的取值范圍有所變化。在敏感性分析中4個(gè)參數(shù)取的都是比較常見(jiàn)的值，而在建立查找表時(shí)還要考慮一些不常見(jiàn)的情況，所以此處參數(shù)取值范圍要大于前面敏感性分析時(shí)的取值。光學(xué)厚度取值02，在0到1之間以0.1步長(zhǎng)變化，在12之間以0.2步長(zhǎng)變化共16個(gè)值；太陽(yáng)天頂角從0度至80度以10度的步長(zhǎng)變化，共有9個(gè)值；感器天頂角從0度至70度以10度的步長(zhǎng)變化，共有8個(gè)值；地表海拔選取0m、200m、500m、1000m、2000m和4000m共6個(gè)值。將上述參數(shù)組合代入6S模型循環(huán)計(jì)算，得到MODIS前7個(gè)波段在不同氣溶膠光6學(xué)厚度、太陽(yáng)方位角、傳感器方位角、地表海拔下的大氣校正系數(shù)xa、xb、XC查找表。4)逐像元校正方法對(duì)MODIS影像進(jìn)行大氣校正時(shí)，從MODISLlB數(shù)據(jù)中導(dǎo)出表觀輻亮度、太陽(yáng)天頂角、傳感器天頂角數(shù)據(jù)，從MODIS的氣溶膠光學(xué)厚度產(chǎn)品中讀出光學(xué)厚度數(shù)據(jù)，結(jié)合Ikm分辨率的SRTMDEM數(shù)據(jù)，逐像元根據(jù)查找表進(jìn)行線性插值，計(jì)算出當(dāng)前像元對(duì)應(yīng)波段的大氣校正系數(shù)xa、xb、xc,并應(yīng)用公式1計(jì)算該像元的地表反射率。整個(gè)校正流程見(jiàn)圖2。5)實(shí)例驗(yàn)證MODIS影像范圍很大，很難像TM大氣校正那樣在衛(wèi)星過(guò)境同時(shí)進(jìn)行地表同步觀測(cè)來(lái)對(duì)校正結(jié)果進(jìn)行驗(yàn)證。前人的研究成果已經(jīng)表明，在大氣參數(shù)、幾何位置信息等輸入?yún)?shù)較為準(zhǔn)確的情況下，小范圍的遙感影像經(jīng)過(guò)6S模型的大氣校正之后得到的地表反射率信息比較準(zhǔn)確。本發(fā)明主要驗(yàn)證的是在比較大的范圍內(nèi)，利用查找表逐像元進(jìn)行大氣校正與整景影像利用統(tǒng)一輸入?yún)?shù)進(jìn)行大氣校正的差別。為了對(duì)本方法的大氣校正結(jié)果進(jìn)行驗(yàn)證，本發(fā)明對(duì)比了利用6S模型逐像元在線計(jì)算、基于查找表逐像元校正以及整景影像利用統(tǒng)一輸入?yún)?shù)的6S模型校正方法對(duì)于M0DISL1B影像的校正結(jié)果。為了陳述方便，三種方法分別簡(jiǎn)稱為Online法(6Sonlinecorrection)、LUT法(Look-uptablemethod)禾口PU法(parameter-uniformmethod)0因?yàn)?S模型在線計(jì)算耗時(shí)太長(zhǎng)，本文取了一景2005年10月17日的AQUA/M0DISLlB數(shù)據(jù)的一個(gè)掃描行(2137個(gè)像元)來(lái)進(jìn)行驗(yàn)證。圖3給出了掃描線上自西向東各個(gè)像元的氣溶膠光學(xué)厚度、太陽(yáng)天頂角、傳感器天頂角和地表海拔。從圖上可以看出，氣溶膠變化比較劇烈，有些像元低于0.2，有些像元高于0.6；太陽(yáng)天頂角自西向東逐漸增大；而傳感器天頂角則從星下點(diǎn)開始向兩邊逐漸增大；海拔高度總體上來(lái)說(shuō)西高東低，西部總體上在1500m至3000m之間，東部海拔一般低于50m，為平原區(qū)域。6S在線逐像元校正過(guò)程中，各個(gè)像元的氣溶膠光學(xué)厚度數(shù)據(jù)、太陽(yáng)天頂角、傳感器天頂角和地表海拔分別來(lái)自各像元對(duì)應(yīng)的數(shù)據(jù)。整景影像利用統(tǒng)一輸入?yún)?shù)的6S模型校正過(guò)程中，氣溶膠光學(xué)厚度取中心像元值0.15，太陽(yáng)天頂角取值49，傳感器天頂角取值0，地表海拔取值30m。三種不同校正方法得到的MODIS第1、2、3、4波段地表反射率見(jiàn)圖4，從圖中可以看出，第1、3和4波段的LUT法和Online法校正結(jié)果很接近，而P-U法到的地表反射率與前兩者的差別要大一些，第2波段(近紅外波段)受到氣溶膠散射影響相對(duì)較小，LUT法、Online法和P-U法這三種方法計(jì)算得到的地表反射率比較接近。以O(shè)nline法的校正計(jì)算結(jié)果為基準(zhǔn)計(jì)算了LUT法和P-U法計(jì)算得到的第14波段反射率的均方根誤差RMSE(表1)來(lái)定量評(píng)價(jià)這兩種方法。從表中也看出，除了第2波段兩種方法的誤差相差較小之外，其余波段查找表法的校正精度要顯著高于統(tǒng)一參數(shù)法，表明了本發(fā)明大氣校正結(jié)果相對(duì)于常用的統(tǒng)一參數(shù)法校正結(jié)果的優(yōu)越性。表ILUT法和P-U法校正結(jié)果的均方根誤差權(quán)利要求一種基于查找表的遙感影像逐像元大氣校正方法，其特征在于包括如下步驟(一)離線模擬(1)采用6S模型進(jìn)行大氣校正，模擬太陽(yáng)輻射在大氣中的傳輸狀況得到地表真實(shí)反射率ρs，并給出大氣校正系數(shù)xa、xb和xc<mrow><mfencedopen='{'close=''><mtable><mtr><mtd><msub><msup><mi>ρ</mi><mo>′</mo></msup><mi>s</mi></msub><mo>=</mo><mfrac><mi>πL</mi><mrow><msub><mi>T</mi><mi>g</mi></msub><msub><mi>E</mi><mi>s</mi></msub><mi>cos</mi><mrow><mo>(</mo><msub><mi>θ</mi><mi>s</mi></msub><mo>)</mo></mrow><mi>T</mi><mrow><mo>(</mo><msub><mi>θ</mi><mi>s</mi></msub><mo>)</mo></mrow><mi>T</mi><mrow><mo>(</mo><msub><mi>θ</mi><mi>v</mi></msub><mo>)</mo></mrow></mrow></mfrac><mo>-</mo><mfrac><mrow><msub><mi>ρ</mi><mrow><mi>R</mi><mo>+</mo><mi>A</mi></mrow></msub><mrow><mo>(</mo><msub><mi>θ</mi><mi>s</mi></msub><mo>,</mo><msub><mi>θ</mi><mi>v</mi></msub><mo>)</mo></mrow></mrow><mrow><mi>T</mi><mrow><mo>(</mo><msub><mi>θ</mi><mi>s</mi></msub><mo>)</mo></mrow><mi>T</mi><mrow><mo>(</mo><msub><mi>θ</mi><mi>v</mi></msub><mo>)</mo></mrow></mrow></mfrac></mtd></mtr><mtr><mtd><msub><mi>ρ</mi><mi>s</mi></msub><mo>=</mo><mfrac><msub><msup><mi>ρ</mi><mo>′</mo></msup><mi>s</mi></msub><mrow><mn>1</mn><mo>+</mo><mi>S</mi><msub><msup><mi>ρ</mi><mo>′</mo></msup><mi>s</mi></msub></mrow></mfrac></mtd></mtr></mtable></mfenced><mo>,</mo></mrow><mrow><mi>xa</mi><mo>=</mo><mfrac><mi>π</mi><mrow><msub><mi>T</mi><mi>g</mi></msub><msub><mi>E</mi><mi>s</mi></msub><mi>cos</mi><mrow><mo>(</mo><msub><mi>θ</mi><mi>s</mi></msub><mo>)</mo></mrow><mi>T</mi><mrow><mo>(</mo><msub><mi>θ</mi><mi>s</mi></msub><mo>)</mo></mrow><mi>T</mi><mrow><mo>(</mo><msub><mi>θ</mi><mi>v</mi></msub><mo>)</mo></mrow></mrow></mfrac><mo>,</mo></mrow><mrow><mi>xb</mi><mo>=</mo><mfrac><mrow><msub><mi>ρ</mi><mrow><mi>R</mi><mo>+</mo><mi>A</mi></mrow></msub><mrow><mo>(</mo><msub><mi>θ</mi><mi>s</mi></msub><mo>,</mo><msub><mi>θ</mi><mi>v</mi></msub><mo>)</mo></mrow></mrow><mrow><mi>T</mi><mrow><mo>(</mo><msub><mi>θ</mi><mi>s</mi></msub><mo>)</mo></mrow><mi>T</mi><mrow><mo>(</mo><msub><mi>θ</mi><mi>v</mi></msub><mo>)</mo></mrow></mrow></mfrac><mo>,</mo></mrow>xc＝S；其中，Tg為大氣中O3、H2O、N2、CO2、O2、CH4氣體對(duì)太陽(yáng)輻射的吸收透過(guò)率，ρR+A為大氣層輻射反射率即氣溶膠程輻射與瑞利散射程輻射之和，θs、θv分別為太陽(yáng)天頂角和傳感器天頂角，T(θs)、T(θv)分別為太陽(yáng)至地面、地面至傳感器的大氣路徑透過(guò)率，S為大氣底層向下的半球反射率，L為表觀輻亮度，Es為對(duì)應(yīng)波長(zhǎng)的大氣層頂太陽(yáng)輻照度；(2)敏感性分析遙感影像，只考慮氣溶膠光學(xué)厚度以及傳感器天頂角的變化，太陽(yáng)天頂角及地表海拔設(shè)為常量；(3)利用6S模型根據(jù)步驟(1)所述的大氣校正系數(shù)xa、xb、xc公式離線計(jì)算不同氣溶膠光學(xué)厚度、太陽(yáng)天頂角、傳感器天頂角以及地表海拔情況下的大氣校正系數(shù)xa、xb、xc，并建立針對(duì)氣溶膠光學(xué)厚度、太陽(yáng)天頂角、傳感器天頂角以及地表海拔相對(duì)于大氣校正系數(shù)xa、xb、xc的查找表；(二)在線校正(4)對(duì)遙感影像進(jìn)行大氣校正時(shí)，從遙感影像數(shù)據(jù)中導(dǎo)出表觀輻亮度、太陽(yáng)天頂角、傳感器天頂角數(shù)據(jù)，從遙感影像的氣溶膠光學(xué)厚度中讀出光學(xué)厚度數(shù)據(jù)，結(jié)合1km分辨率的SRTMDEM數(shù)據(jù)，逐像元根據(jù)步驟(3)所述的查找表進(jìn)行線性插值，計(jì)算出當(dāng)前像元對(duì)應(yīng)波段的大氣校正系數(shù)xa、xb、xc。全文摘要本發(fā)明公布了一種基于查找表的遙感影像逐像元大氣校正方法，本發(fā)明方法利用6S大氣輻射傳輸模型離線計(jì)算建立了不同氣溶膠光學(xué)厚度、太陽(yáng)天頂角、傳感器天頂角以及地表海拔條件下大氣校正系數(shù)的查找表，基于該查找表對(duì)MODIS遙感影像進(jìn)行逐像元大氣校正。本發(fā)明可以有效改善由于大氣條件、傳感器位置等空間分布差異對(duì)MODIS圖像大氣校正的影響。文檔編號(hào)G01S7/487GK101915914SQ20101024050公開日2010年12月15日申請(qǐng)日期2010年7月30日優(yōu)先權(quán)日2010年7月30日發(fā)明者徐永明,趙巧華,陳愛(ài)軍申請(qǐng)人:南京信息工程大學(xué)