本發(fā)明涉及熱紅外遙感監(jiān)測(cè)領(lǐng)域。尤其涉及一種基于靜止衛(wèi)星數(shù)據(jù)的林火和草原火遙感監(jiān)測(cè)方法。

背景技術(shù)：

林火和草原火災(zāi)會(huì)破壞林業(yè)或畜牧業(yè)生產(chǎn)，造成重大經(jīng)濟(jì)損失。同時(shí)，火災(zāi)釋放的大量溫室氣體和氣溶膠會(huì)改變大氣化學(xué)成分，對(duì)環(huán)境和氣候變化產(chǎn)生很大影響。因此，研究森林和草原火災(zāi)監(jiān)測(cè)和預(yù)警技術(shù)具有重要的現(xiàn)實(shí)意義和長(zhǎng)遠(yuǎn)的歷史意義。遙感技術(shù)手段客觀性好、整體性強(qiáng)、反應(yīng)速度快、耗費(fèi)省，可以在森林和草原火災(zāi)災(zāi)害監(jiān)測(cè)、預(yù)警及損失評(píng)估工作中發(fā)揮巨大的作用。相較于應(yīng)用極軌衛(wèi)星進(jìn)行林火和草原火遙感監(jiān)測(cè)技術(shù)，靜止衛(wèi)星始終與地球表面保持相對(duì)靜止同時(shí)，靜止衛(wèi)星覆蓋區(qū)域更廣，一顆衛(wèi)星可觀測(cè)地球中低緯度地區(qū)約1/3的范圍，具有很強(qiáng)的機(jī)動(dòng)靈活性，能夠?qū)μ囟▍^(qū)域進(jìn)行分鐘級(jí)高重復(fù)觀測(cè)，可快速監(jiān)測(cè)災(zāi)害目標(biāo)的動(dòng)態(tài)變化。利用靜止衛(wèi)星開展災(zāi)害監(jiān)測(cè)評(píng)估業(yè)務(wù)，特別適用于致災(zāi)因子危險(xiǎn)性分析、災(zāi)害范圍監(jiān)測(cè)和風(fēng)險(xiǎn)評(píng)估，在減災(zāi)救災(zāi)中具有廣闊的應(yīng)用前景。

目前，國(guó)內(nèi)外主要的在軌靜止衛(wèi)星有美國(guó)國(guó)家海洋和大氣管理局(NOAA)的地球靜止軌道業(yè)務(wù)環(huán)境衛(wèi)星(GOES)、日本多功能運(yùn)輸衛(wèi)星(MTSAT)、歐洲氣象衛(wèi)星組織(EUMETSAT)的第二代氣象衛(wèi)星(MSG)以及中國(guó)的風(fēng)云(FY)系列衛(wèi)星。日本新一代靜止氣象衛(wèi)星Himawari-8(葵花8號(hào))于2014年10月發(fā)射，于2015年7月正式開始廣播發(fā)送衛(wèi)星觀測(cè)數(shù)據(jù)。Himawari-8所搭載的輻射成像儀(AHI)包含16通道，由3個(gè)可見光、3個(gè)近紅外和10個(gè)紅外通道構(gòu)成。其中，可見光通道空間分辨率達(dá)到0.5km，近紅外和紅外通道空間分辨率達(dá)到1～2km，全盤觀測(cè)頻率達(dá)到每10min一次。最大程度地利用Himawari-8衛(wèi)星，將極大地提高對(duì)臺(tái)風(fēng)、暴雨、火災(zāi)等災(zāi)害的監(jiān)測(cè)和預(yù)報(bào)精度，并為相關(guān)的決策支持服務(wù)、公眾服務(wù)、科學(xué)研究工作提供支持。

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

本發(fā)明提出一種靜止衛(wèi)星火災(zāi)遙感監(jiān)測(cè)方法，解決基于靜止衛(wèi)星的森林和草原火點(diǎn)自動(dòng)提取、火災(zāi)蔓延態(tài)勢(shì)監(jiān)測(cè)和火災(zāi)損失評(píng)估的問題。

本申請(qǐng)實(shí)施例提供一種靜止衛(wèi)星火災(zāi)遙感監(jiān)測(cè)方法，包括以下步驟：

采集靜止衛(wèi)星數(shù)據(jù)，計(jì)算像元亮度溫度；

對(duì)陸地像元中非云、非耀斑的像元，用亮度溫度閾值的方法判別潛在火點(diǎn)像元；

對(duì)所述潛在火點(diǎn)像元，利用背景窗口像元波段信息統(tǒng)計(jì)判別真實(shí)火點(diǎn)像元；

基于固定火點(diǎn)庫和土地利用類型數(shù)據(jù)，剔除常年高溫火點(diǎn)，并提取林火和草原火位置、亮溫信息。

優(yōu)選地，所述靜止衛(wèi)星火災(zāi)遙感監(jiān)測(cè)方法還包括以下步驟：基于土地利用類型數(shù)據(jù)生成水陸掩碼，識(shí)別陸地像元。

優(yōu)選地，所述靜止衛(wèi)星火災(zāi)遙感監(jiān)測(cè)方法還包括以下步驟：基于靜止衛(wèi)星的云掩碼產(chǎn)品，識(shí)別云像元。

優(yōu)選地，所述靜止衛(wèi)星火災(zāi)遙感監(jiān)測(cè)方法還包括以下步驟：應(yīng)用中紅外4μm波段和熱紅外11μm波段的組合訂正水汽吸收。

優(yōu)選地，所述靜止衛(wèi)星火災(zāi)遙感監(jiān)測(cè)方法還包括以下步驟：同時(shí)滿足耀斑角小于10°、紅波段表觀反射率大于0.1和短波近紅外波段表觀反射率大于0.2三個(gè)條件的像元，確認(rèn)為太陽耀斑像元。

優(yōu)選地，所述亮度溫度閾值的方法是，滿足以下至少一項(xiàng)條件式的像元為潛在火點(diǎn)像元：

白天：T₄>[315+25cosθ]K 公式1

白天：T₄>[305+10cosθ]K&T₄-T₁₁>15K 公式2

夜晚：T₄>295K 公式3

夜晚：T₄>285K&T₄-T₁₁>10K 公式4

優(yōu)選地，所述利用背景窗口像元波段信息統(tǒng)計(jì)判別真實(shí)火點(diǎn)像元的步驟，進(jìn)一步包含：

以所述潛在火點(diǎn)像元為中心建立背景窗口，遍歷背景窗口中所有像元，將云像元、非陸地像元、耀斑像元和潛在火點(diǎn)像元作為無效像元，其他像元為有效像元，使背景窗口中無效像元數(shù)量小于像元總數(shù)的20％；

統(tǒng)計(jì)背景窗口內(nèi)波段信息，同時(shí)滿足以下兩個(gè)條件式的潛在火點(diǎn)像元為真實(shí)火點(diǎn)像元：

(T₄-T₁₁)>δT_4,11×C₁ 公式5

(T₄-Tb₄)>δT₄×C₂ 公式6

進(jìn)一步地，所述真實(shí)火點(diǎn)像元還需要滿足以下至少一個(gè)條件式：

(T₄-T₁₁)>δT_4,11&(T₁₁-Tb₁₁)>-20K 公式7

(R_4,11-Rb_4,11)>δR_4,11×C₃+0.5×C₄ 公式8

在公式1～8中，θ為太陽天頂角；T₄為中紅外4μm波段亮度溫度；T₁₁為熱紅外11μm波段亮度溫度；Tb₄和Tb₁₁分別為背景窗口中所有有效像元中紅外4μm波段和熱紅外11μm波段亮度溫度的平均值；δT_4,11表示背景窗口中所有有效像元中紅外4μm波段與熱紅外11μm波段亮度溫度差值的計(jì)算標(biāo)準(zhǔn)偏差，δT₄表示背景窗口中所有有效像元中紅外4μm波段亮度溫度的計(jì)算標(biāo)準(zhǔn)偏差；R_4,11＝(R₄-R₁₁)/R₄×10，R₄和R₁₁分別為中紅外4μm波段和熱紅外11μm波段輻亮度值；Rb_4,11為背景窗口中所有有效像元R_4,11的平均值；δR_4,11表示背景窗口中所有有效像元R_4,11的計(jì)算標(biāo)準(zhǔn)偏差，C₁、C₂、C₃、C₄均為經(jīng)驗(yàn)常數(shù)。

優(yōu)選地，本申請(qǐng)任意一項(xiàng)所述靜止衛(wèi)星火災(zāi)遙感監(jiān)測(cè)方法還包括以下步驟：根據(jù)靜止衛(wèi)星觀測(cè)頻次的林火和草原火遙感監(jiān)測(cè)結(jié)果，給出火災(zāi)蔓延態(tài)勢(shì)序列。

優(yōu)選地，本申請(qǐng)任意一項(xiàng)所述靜止衛(wèi)星火災(zāi)遙感監(jiān)測(cè)方法還包括以下步驟：基于火災(zāi)發(fā)生前期MODIS凈初級(jí)生產(chǎn)力產(chǎn)品，對(duì)林火或草原火災(zāi)所造成的損失進(jìn)行評(píng)估。

本申請(qǐng)實(shí)施例采用的上述至少一個(gè)技術(shù)方案能夠達(dá)到以下有益效果：本發(fā)明的方法憑借靜止衛(wèi)星觀測(cè)頻次高的優(yōu)勢(shì)，利用靜止衛(wèi)星數(shù)據(jù)和土地利用類型數(shù)據(jù)，針對(duì)我國(guó)人煙稀少的草原和森林地區(qū)進(jìn)行準(zhǔn)實(shí)時(shí)的火災(zāi)監(jiān)測(cè)，并對(duì)災(zāi)后損失進(jìn)行評(píng)估。由于靜止衛(wèi)星具有很強(qiáng)的機(jī)動(dòng)靈活性，能夠?qū)μ囟▍^(qū)域進(jìn)行分鐘級(jí)高重復(fù)觀測(cè)，可以快速監(jiān)測(cè)災(zāi)害目標(biāo)的動(dòng)態(tài)變化，有助于進(jìn)一步提高我國(guó)應(yīng)對(duì)重大災(zāi)害的預(yù)警和應(yīng)急處理率，為我國(guó)減災(zāi)救災(zāi)工作提供技術(shù)支持。

附圖說明

此處所說明的附圖用來提供對(duì)本申請(qǐng)的進(jìn)一步理解，構(gòu)成本申請(qǐng)的一部分，本申請(qǐng)的示意性實(shí)施例及其說明用于解釋本申請(qǐng)，并不構(gòu)成對(duì)本申請(qǐng)的不當(dāng)限定。在附圖中：

圖1為基于一種靜止衛(wèi)星數(shù)據(jù)的火災(zāi)遙感監(jiān)測(cè)方法實(shí)施例示意圖；

圖2為背景窗口波段信息統(tǒng)計(jì)示意圖；

圖3為背景窗口波段信息統(tǒng)計(jì)方法流程示意圖；

圖4為本發(fā)明包含火災(zāi)損失評(píng)估的實(shí)施例流程示意圖。

具體實(shí)施方式

為使本申請(qǐng)的目的、技術(shù)方案和優(yōu)點(diǎn)更加清楚，下面將結(jié)合本申請(qǐng)具體實(shí)施例及相應(yīng)的附圖對(duì)本申請(qǐng)技術(shù)方案進(jìn)行清楚、完整地描述。顯然，所描述的實(shí)施例僅是本申請(qǐng)一部分實(shí)施例，而不是全部的實(shí)施例?；诒旧暾?qǐng)中的實(shí)施例，本領(lǐng)域普通技術(shù)人員在沒有做出創(chuàng)造性勞動(dòng)前提下所獲得的所有其他實(shí)施例，都屬于本申請(qǐng)保護(hù)的范圍。

本發(fā)明的方法主要基于靜止衛(wèi)星數(shù)據(jù)，輔助以土地利用類型數(shù)據(jù)和歷史固定火點(diǎn)數(shù)據(jù)庫，自動(dòng)提取森林和草原火點(diǎn)，快速監(jiān)測(cè)火災(zāi)動(dòng)態(tài)蔓延趨勢(shì)，并在災(zāi)后對(duì)受災(zāi)損失作出評(píng)估。從而，本發(fā)明能夠針對(duì)林火和草原火的突發(fā)性特點(diǎn)，利用靜止衛(wèi)星實(shí)現(xiàn)準(zhǔn)實(shí)時(shí)的同步監(jiān)測(cè)，在第一時(shí)間獲取災(zāi)情發(fā)生地區(qū)空間分布范圍，并動(dòng)態(tài)把握災(zāi)害變化宏觀趨勢(shì)。同時(shí)，在應(yīng)急救災(zāi)階段可提供火燒跡地面積等用于評(píng)估災(zāi)情和指導(dǎo)救災(zāi)行動(dòng)，在災(zāi)情穩(wěn)定后，為災(zāi)害損失評(píng)估提供數(shù)據(jù)支持，以及為災(zāi)后重建和恢復(fù)工作提供參考依據(jù)。

本申請(qǐng)文件中，所述表觀反射率為大氣層頂?shù)姆瓷渎剩隽炼葴囟仁俏矬w的輻射能量用同輻射量的黑體溫度表示，所述輻亮度為單位投影面積、單位立體角上的輻射通量。

以下結(jié)合附圖，詳細(xì)說明本申請(qǐng)各實(shí)施例提供的技術(shù)方案。

圖1為基于一種靜止衛(wèi)星數(shù)據(jù)的火災(zāi)遙感監(jiān)測(cè)方法實(shí)施例示意圖；本申請(qǐng)實(shí)施例提供一種靜止衛(wèi)星火災(zāi)遙感監(jiān)測(cè)方法，包括以下步驟：

步驟11、采集靜止衛(wèi)星數(shù)據(jù)，計(jì)算像元亮度溫度；

步驟12、對(duì)陸地像元中非云、非耀斑的像元，用亮度溫度閾值的方法判別潛在火點(diǎn)像元；

步驟13、對(duì)所述潛在火點(diǎn)像元，利用背景窗口像元波段信息統(tǒng)計(jì)判別真實(shí)火點(diǎn)像元；

步驟14、基于固定火點(diǎn)庫和土地利用類型數(shù)據(jù)，剔除常年高溫火點(diǎn)，并提取林火和草原火位置、亮溫信息。

以下詳細(xì)說明以上各個(gè)步驟11～14的具體方案。其中，步驟11涉及公式9，步驟12涉及公式1～4，步驟13涉及公式5～8、圖2～3。

步驟11、采集靜止衛(wèi)星數(shù)據(jù)，計(jì)算像元亮度溫度。對(duì)靜止衛(wèi)星近紅外-紅外波段數(shù)據(jù)計(jì)算亮度溫度值，其求解公式如下：

T(L,B)＝C₂/L×ln(C₁/(L⁵×B(L,T)×10⁶)+1) 公式9

式中，B(L,T)為普朗克輻射值[Wm^-2sr^-1μm^-1]，C₁＝2hc²＝1.1910439×10^-16[W m²]，C₂＝(hc)/k，T(L,B)為亮度溫度[K]。其中，L為波長(zhǎng)[m]，h＝6.626×10^-34[J K]為普朗克常數(shù)，c＝3×10⁸[m/s]為光速，k＝1.38×10^-23[J/K]為玻爾茲曼常數(shù)，T為絕對(duì)溫度[K]。

優(yōu)選地，所述靜止衛(wèi)星為Himawari-8，所述數(shù)據(jù)為Himawari-8/AHI數(shù)據(jù)。用公式9對(duì)數(shù)據(jù)進(jìn)行預(yù)處理，計(jì)算Himawari-8/AHI的第7和14波段亮度溫度值。根據(jù)Himawari-8/AHI第7波段數(shù)據(jù)，計(jì)算中紅外4μm波段亮度溫度T₄；Himawari-8/AHI第14波段數(shù)據(jù)，計(jì)算熱紅外11μm波段亮度溫度T₁₁。

步驟12、對(duì)陸地像元中非云、非耀斑的像元，用亮度溫度閾值的方法判別潛在火點(diǎn)像元。

需要指出的是，要獲得陸地像元中的非云、非耀斑像元，需要對(duì)所述像元進(jìn)行預(yù)處理，首先識(shí)別陸地像元、云像元、耀斑像元。具體地，還包括以下步驟：

步驟12A、基于土地利用類型數(shù)據(jù)生成水陸掩碼，識(shí)別陸地像元。

步驟12B、基于靜止衛(wèi)星的云掩碼產(chǎn)品，識(shí)別云像元；優(yōu)選地，基于Himawari-8/AHI的云掩碼產(chǎn)品，識(shí)別云像元。

步驟12C、應(yīng)用中紅外4μm波段和熱紅外11μm波段的組合訂正水汽吸收；優(yōu)選地，應(yīng)用Himawari-8/AHI的第7波段(3.9μm)和第14波段(11.2μm)的組合來訂正水汽吸收。

應(yīng)用中紅外4μm波段和熱紅外11μm波段的組合訂正水汽吸收，主要是利用6S等國(guó)際通用輻射傳輸模型進(jìn)行大氣校正，訂正水汽吸收對(duì)上述兩個(gè)波段輻亮度的影響。

步驟12D、同時(shí)滿足耀斑角小于10°、紅波段表觀反射率大于0.1和短波近紅外波段表觀反射率大于0.2三個(gè)條件的像元，確認(rèn)為太陽耀斑像元。優(yōu)選地，所述紅波段采用Himawari-8/AHI的第3波段(0.64μm)，所述近紅外波段采用Himawari-8/AHI的第4波段(0.86)。

步驟12E、用亮度溫度閾值的方法判別潛在火點(diǎn)像元；

所述亮度溫度閾值的方法是，滿足以下至少一項(xiàng)條件式的像元為潛在火點(diǎn)像元：

白天：T₄>[315+25cosθ]K 公式1

白天：T₄>[305+10cosθ]K&T₄-T₁₁>15K 公式2

夜晚：T₄>295K 公式3

夜晚：T₄>285K&T₄-T₁₁>10K 公式4

在公式1～4中，θ為太陽天頂角；T₄為中紅外4μm波段亮度溫度；T₁₁為熱紅外11μm波段亮度溫度。

步驟13、對(duì)所述潛在火點(diǎn)像元，利用背景窗口像元波段信息統(tǒng)計(jì)判別真實(shí)火點(diǎn)像元。

圖2為背景窗口波段信息統(tǒng)計(jì)示意圖；圖3為背景窗口波段信息統(tǒng)計(jì)方法流程示意圖。優(yōu)選地，所述利用背景窗口像元波段信息統(tǒng)計(jì)判別真實(shí)火點(diǎn)像元的步驟，進(jìn)一步包含：

步驟13A、以所述潛在火點(diǎn)像元為中心建立背景窗口，遍歷背景窗口中所有像元，除潛在火點(diǎn)像元外，將云像元、非陸地像元和耀斑像元作為無效像元，其他像元作為有效像元，使背景窗口中無效像元數(shù)量小于像元總數(shù)的20％。

例如，圖2所示，以潛在火點(diǎn)像元1為中心，建立大小為w×w的背景窗口(w的初值設(shè)置為5)，背景窗口中所有像元，除潛在火點(diǎn)像元1外，將云像元2、水體像元3和耀斑像元4均判斷為無效像元，其他像元判斷為有效像元，若無效像元總數(shù)大于背景窗口中總像元數(shù)的20％，則窗口大小增加至(w+10)×(w+10)個(gè)像元，否則，窗口大小保持不變，其中，窗口大小調(diào)整的最大次數(shù)為20次。

步驟13B、統(tǒng)計(jì)背景窗口內(nèi)波段信息，選擇同時(shí)滿足以下兩個(gè)條件式的潛在火點(diǎn)像元為真實(shí)火點(diǎn)像元：

(T₄-T₁₁)>δT_4,11×C₁ 公式5

(T₄-Tb₄)>δT₄×C₂ 公式6

步驟13C、進(jìn)一步地，所述真實(shí)火點(diǎn)像元還需要滿足以下公式7、公式8中至少一個(gè)條件式：

(T₄-T₁₁)>δT_4,11&(T₁₁-Tb₁₁)>-20K 公式7

(R_4,11-Rb_4,11)>δR_4,11×C₃+0.5×C₄ 公式8

在公式4～8中，T₄為中紅外4μm波段亮度溫度；T₁₁為熱紅外11μm波段亮度溫度；Tb₄和Tb₁₁分別為背景窗口中所有有效像元中紅外4μm波段和熱紅外11μm波段亮度溫度的平均值；δT_4,11表示背景窗口中所有有效像元中紅外4μm波段與熱紅外11μm波段亮度溫度差值的計(jì)算標(biāo)準(zhǔn)偏差，δT₄表示背景窗口中所有有效像元中紅外4μm波段亮度溫度的計(jì)算標(biāo)準(zhǔn)偏差；R_4,11＝(R₄-R₁₁)/R₄×10，R₄和R₁₁分別為中紅外4μm波段和熱紅外11μm波段輻亮度值；δR_4,11表示背景窗口中所有有效像元R_4,11的計(jì)算標(biāo)準(zhǔn)偏差；Rb_4,11為背景窗口中所有有效像元R_4,11的平均值；C₁、C₂、C₃、C₄均為經(jīng)驗(yàn)常數(shù)。

針對(duì)Himawari-8衛(wèi)星的AHI傳感器，優(yōu)選地，取值為C₁＝2，C₂＝2.5，C₃＝2.5，C₄＝0.5。由于不同衛(wèi)星傳感器波段設(shè)置和儀器本身的差異，C₁、C₂、C₃、C₄會(huì)有所不同。

步驟14、基于固定火點(diǎn)庫和土地利用類型數(shù)據(jù)，剔除常年高溫火點(diǎn)，并提取林火和草原火位置、亮溫信息。具體包含以下步驟，

步驟14A、結(jié)合固定火點(diǎn)庫，剔除所述潛在火點(diǎn)像元和真實(shí)火點(diǎn)像元中的常年高溫火點(diǎn)；

步驟14B、將土地利用類型數(shù)據(jù)與所述火點(diǎn)像元進(jìn)行疊加，若所述火點(diǎn)像元對(duì)應(yīng)的土地利用類型為林地，則判斷為林火，若所述火點(diǎn)像元對(duì)應(yīng)的土地利用類型為草原，則判斷該火點(diǎn)為草原火；

步驟14C、輸出林火和草原火對(duì)應(yīng)的經(jīng)緯度、中紅外波段(4μm)和熱紅外波段(11μm)的亮度溫度值。

圖4為本發(fā)明包含火災(zāi)損失評(píng)估的實(shí)施例流程示意圖。

步驟11、采集靜止衛(wèi)星數(shù)據(jù)，計(jì)算像元亮度溫度；

步驟12、對(duì)陸地像元中非云、非耀斑的像元，用亮度溫度閾值的方法判別潛在火點(diǎn)像元；

步驟13、對(duì)所述潛在火點(diǎn)像元，利用背景窗口像元波段信息統(tǒng)計(jì)判別真實(shí)火點(diǎn)像元；

步驟14、基于固定火點(diǎn)庫和土地利用類型數(shù)據(jù)，剔除常年高溫火點(diǎn)，并提取林火和草原火位置、亮溫信息。

步驟15、重復(fù)步驟11～14，根據(jù)靜止衛(wèi)星觀測(cè)頻次的林火和草原火遙感監(jiān)測(cè)結(jié)果，給出火災(zāi)蔓延態(tài)勢(shì)序列。例如，針對(duì)Himawari-8/AHI的觀測(cè)頻次，每10分鐘輸出一次林火和草原或遙感監(jiān)測(cè)數(shù)據(jù)。

步驟16、基于火災(zāi)發(fā)生前期MODIS凈初級(jí)生產(chǎn)力(NPP)產(chǎn)品，對(duì)林火或草原火災(zāi)所造成的損失進(jìn)行評(píng)估。

具體地，包含以下步驟：

步驟16A、獲取火災(zāi)發(fā)生地區(qū)多個(gè)時(shí)段的靜止衛(wèi)星遙感影像，按時(shí)間序列構(gòu)建系列的光譜指數(shù)值，最后通過系列光譜指數(shù)值的合成，能夠得到更好的突出火燒跡地特征的光譜指數(shù)圖像，并進(jìn)而計(jì)算火燒跡地面積大小；

步驟16B、通過生物量估算模型，利用火災(zāi)發(fā)生前期的NPP和火燒跡地面積計(jì)算林地或草原的生物量損失；

步驟16C、根據(jù)林地或草原生物量損失，建立不同生態(tài)系統(tǒng)的經(jīng)濟(jì)模型，計(jì)算林火或草原火所造成的經(jīng)濟(jì)損失。

需要說明的是，步驟16中的凈初級(jí)生產(chǎn)力(NPP)，是用于表征植物把無機(jī)碳(CO₂)固定、轉(zhuǎn)化為有機(jī)碳這一過程的能力，即植被所固定的有機(jī)碳中扣除本身呼吸消耗的部分，這一部分用于植被的生長(zhǎng)和生殖。MODIS官方網(wǎng)站可提供覆蓋全球范圍的8天合成NPP產(chǎn)品。

還需要說明的是，在本說明書的全部公式中，邏輯運(yùn)算符“&”表示“并且”；本申請(qǐng)文件中，數(shù)學(xué)運(yùn)算符的優(yōu)先級(jí)高于邏輯運(yùn)算符。

本領(lǐng)域內(nèi)的技術(shù)人員應(yīng)明白，本發(fā)明的實(shí)施例可提供為方法、系統(tǒng)、或計(jì)算機(jī)程序產(chǎn)品。因此，本發(fā)明可采用完全硬件實(shí)施例、完全軟件實(shí)施例、或結(jié)合軟件和硬件方面的實(shí)施例的形式。而且，本發(fā)明可采用在一個(gè)或多個(gè)其中包含有計(jì)算機(jī)可用程序代碼的計(jì)算機(jī)可用存儲(chǔ)介質(zhì)(包括但不限于磁盤存儲(chǔ)器、CD-ROM、光學(xué)存儲(chǔ)器等)上實(shí)施的計(jì)算機(jī)程序產(chǎn)品的形式。

本發(fā)明是參照根據(jù)本發(fā)明實(shí)施例的方法、設(shè)備(系統(tǒng))、和計(jì)算機(jī)程序產(chǎn)品的流程圖和/或方框圖來描述的。應(yīng)理解可由計(jì)算機(jī)程序指令實(shí)現(xiàn)流程圖和/或方框圖中的每一流程和/或方框、以及流程圖和/或方框圖中的流程和/或方框的結(jié)合。可提供這些計(jì)算機(jī)程序指令到通用計(jì)算機(jī)、專用計(jì)算機(jī)、嵌入式處理機(jī)或其他可編程數(shù)據(jù)處理設(shè)備的處理器以產(chǎn)生一個(gè)機(jī)器，使得通過計(jì)算機(jī)或其他可編程數(shù)據(jù)處理設(shè)備的處理器執(zhí)行的指令產(chǎn)生用于實(shí)現(xiàn)在流程圖一個(gè)流程或多個(gè)流程和/或方框圖一個(gè)方框或多個(gè)方框中指定的功能的裝置。

還需要說明的是，術(shù)語“包括”、“包含”或者其任何其他變體意在涵蓋非排他性的包含，從而使得包括一系列要素的過程、方法、商品或者設(shè)備不僅包括那些要素，而且還包括沒有明確列出的其他要素，或者是還包括為這種過程、方法、商品或者設(shè)備所固有的要素。在沒有更多限制的情況下，由語句“包括一個(gè)……”限定的要素，并不排除在包括所述要素的過程、方法、商品或者設(shè)備中還存在另外的相同要素。

以上所述僅為本申請(qǐng)的實(shí)施例而已，并不用于限制本申請(qǐng)。對(duì)于本領(lǐng)域技術(shù)人員來說，本申請(qǐng)可以有各種更改和變化。凡在本申請(qǐng)的精神和原理之內(nèi)所作的任何修改、等同替換、改進(jìn)等，均應(yīng)包含在本申請(qǐng)的權(quán)利要求范圍之內(nèi)。