基于被動微波遙感的青藏高原雪水當(dāng)量估算方法及系統(tǒng)的制作方法
【專利摘要】本發(fā)明提供一種基于被動微波遙感的青藏高原雪水當(dāng)量估算方法及系統(tǒng)，方法包括：S1、從被動微波遙感數(shù)據(jù)中剔除包含各種干擾因子的像元，保留可能存在積雪像元；S2、通過閾值法將可能存在積雪像元分為中厚雪區(qū)域、淺雪區(qū)域及無雪區(qū)域；S3.確定各地物下墊面的雪深經(jīng)驗(yàn)反演算法；S4、對于中厚雪區(qū)域，檢索為某一地物時調(diào)用該地物下墊面的雪深經(jīng)驗(yàn)反演算法計(jì)算雪水當(dāng)量，得到中厚雪區(qū)域的雪水當(dāng)量；S5、設(shè)定淺雪區(qū)域的雪水當(dāng)量；S6、根據(jù)得到的中厚雪區(qū)域的雪水當(dāng)量和設(shè)定的淺雪區(qū)域的雪水當(dāng)量，得到青藏高原雪水當(dāng)量數(shù)據(jù)?？筛纳魄嗖馗咴┧a(chǎn)品精度、實(shí)現(xiàn)對長時間系列遙感數(shù)據(jù)自動批量處理，提供青藏高原時間序列的相關(guān)積雪產(chǎn)品數(shù)據(jù)集。
【專利說明】
基于被動微波遙感的青藏高原雪水當(dāng)量估算方法及系統(tǒng)
技術(shù)領(lǐng)域
[0001] 本發(fā)明涉及遙感數(shù)據(jù)地學(xué)參數(shù)反演技術(shù)領(lǐng)域，具體涉及一種基于被動微波遙感的 青藏高原雪水當(dāng)量估算方法及系統(tǒng)。
【背景技術(shù)】
[0002] 目前國內(nèi)外有許多與積雪有關(guān)的系統(tǒng)，大體可以分為兩類，一類是專門的積雪面 積、覆蓋度和積雪雪水當(dāng)量產(chǎn)品生產(chǎn)系統(tǒng)，往往積雪面積和覆蓋度的生產(chǎn)較多（可靠）。其 中，國外的產(chǎn)品生產(chǎn)系統(tǒng)有:美國國家航空和宇宙航行局(NASA)的基于Linux的中分辨率成 像光譜儀的數(shù)據(jù)處理系統(tǒng)(M0DAPS)和先進(jìn)微波掃描輻射計(jì)(AMSR-E)產(chǎn)品生產(chǎn)系統(tǒng);美國國 家海洋和大氣局的交互式多傳感器冰雪制圖系統(tǒng)（MS)、微波綜合反演系統(tǒng)(MIRS)和自動 雪圖生產(chǎn)系統(tǒng);歐空局的基于ENVI/IDL的積雪覆蓋度產(chǎn)品生產(chǎn)系統(tǒng)和采用MATLAB語言利用 微波和地面觀測點(diǎn)生產(chǎn)雪圖系統(tǒng)等。國內(nèi)的主要是風(fēng)云三號衛(wèi)星的產(chǎn)品生產(chǎn)系統(tǒng)(PGS)和 北半球積雪監(jiān)測診斷業(yè)務(wù)系統(tǒng)。另一類是基于積雪產(chǎn)品的預(yù)警監(jiān)測系統(tǒng)，這類應(yīng)用在國內(nèi) 外比較多見，涉及氣候、洪水、干旱、雪災(zāi)的預(yù)警、作物估產(chǎn)等方面。國內(nèi)部分有:鄧曉東等開 發(fā)的基于AVHRR(The Advanced Very High Resolution Radiometer，先進(jìn)甚高分辨率福射 儀)資料的內(nèi)蒙古積雪監(jiān)測業(yè)務(wù)系統(tǒng)，傅華等人針對北疆的積雪監(jiān)測系統(tǒng)，邊多等人的基于 MODIS(Moderate_Resolution Imaging Spectroradiometer，中分辨率成像光譜儀）資料的 西藏遙感積雪監(jiān)測業(yè)務(wù)化系統(tǒng)，候慧珠等的基于M0DIS影像的內(nèi)蒙古草原監(jiān)測系統(tǒng)等。國外 有:加拿大的氣象中心業(yè)務(wù)化系統(tǒng)、利用AVHRR數(shù)據(jù)的阿爾卑斯山監(jiān)測系統(tǒng)、瑞典洪水預(yù)警 系統(tǒng)等。這些系統(tǒng)的特點(diǎn):考慮降低系統(tǒng)運(yùn)營成本和提高系統(tǒng)的運(yùn)行速度，多采用Linux系 統(tǒng)平臺。系統(tǒng)穩(wěn)健、可靠、易于維護(hù)和更新、支持全天24小時不間斷工作。編譯語言主要有 F0RTRAN77、VisualC++6.0、MicrosoftVisualC++.NE。系統(tǒng)多采用Windows 或 Linux 操作系統(tǒng) 平臺。以上產(chǎn)品生產(chǎn)系統(tǒng)，除了北半球積雪監(jiān)測診斷業(yè)務(wù)系統(tǒng)的中國區(qū)雪深產(chǎn)品和積雪日 數(shù)產(chǎn)品采用常規(guī)觀測數(shù)據(jù)生產(chǎn)外，其他產(chǎn)品生產(chǎn)系統(tǒng)的大部分產(chǎn)品是通過遙感數(shù)據(jù)反演獲 得。
[0003] 陸地積雪是全球變化的敏感因子，是冰凍圈中季節(jié)性變化最為快速的因子之一。 青藏高原積雪變化不僅對高原本身，還對下游乃至全球的氣候、環(huán)境、水資源等產(chǎn)生重大影 響。雪災(zāi)是青藏高原牧區(qū)的主要自然災(zāi)害，青藏高原地區(qū)現(xiàn)有積雪監(jiān)測能力薄弱，致使牧區(qū) 人民生命財(cái)產(chǎn)遭受重大損失。因此，準(zhǔn)確獲取不同尺度、不同范圍的高精度積雪雪水當(dāng)量信 息不僅能夠提高現(xiàn)有青藏高原地區(qū)的雪災(zāi)預(yù)警和監(jiān)測能力，降低災(zāi)害帶來的損失，具有巨 大的社會經(jīng)濟(jì)效益，同時對于研究青藏高原地區(qū)的大氣科學(xué)研究和全球水體和能量循環(huán)具 有重大的科學(xué)價值。
[0004] 然而，現(xiàn)有已發(fā)布積雪數(shù)據(jù)集在青藏高原地區(qū)精度有限，特別青藏高原積雪積雪 分布的時空規(guī)律性較北半球高煒度地區(qū)復(fù)雜且薄雪分布較多，不能滿足研究應(yīng)用要求。當(dāng) 前在國內(nèi)外有多套積雪數(shù)據(jù)集，如北半球頂S日冰雪蓋數(shù)據(jù)(美國冰學(xué)數(shù)據(jù)中心，1997年至 今）、全球SSM/1(Special Sensor Microwave/Image，特殊傳感器微波/圖片）冰密集度和雪 蓋范圍數(shù)據(jù)(嫩34，1995年至今）、全球41〇?4積雪數(shù)據(jù)(嫩3六，2002年至今）、中國區(qū)域3]\?? (Scanning Multichannel Microwave Radiometer，多通道掃描微波福射計(jì)）和SSM/I 日雪 深/雪水當(dāng)量數(shù)據(jù)（中科院寒旱所，1978~2005年)等等。但是，這些數(shù)據(jù)集是針對全球或中 國區(qū)域的，并未開發(fā)專門針對青藏高原地區(qū)的數(shù)據(jù)產(chǎn)品。同時，由于微波積雪產(chǎn)品空間分辨 率低，又受地面觀測資料不足的影響，數(shù)據(jù)集在青藏高原地區(qū)缺乏系統(tǒng)性的驗(yàn)證，在青藏高 原地區(qū)的精度有限，業(yè)務(wù)化程度不高等特點(diǎn)。
[0005] 國際上的積雪遙感算法在全球尺度上取得了一定的效果，但這些算法在青藏高原 地區(qū)則出現(xiàn)很大的估算誤差。國外的多數(shù)算法適用于高煒度地區(qū)，且算法中各參數(shù)具有地 域局限性，因此這些算法用到青藏高原地區(qū)就會出現(xiàn)SWE(Snow Water Equivalent，雪水當(dāng) 量)被高估的現(xiàn)象。且國際算法多考慮的是中厚雪的反演(主要是利用37GHz和18/19GHZ的 極化差來估算相對較厚的積雪），但是考慮到青藏高原地區(qū)積雪的物理性質(zhì)，其并不滿足國 際算法的前提，且雪深基本小于30cm，當(dāng)前算法大多是在除青藏高原外的地區(qū)有所成效，具 體應(yīng)用到青藏高原也會出現(xiàn)高估現(xiàn)象。加之在青藏高原的地面觀測點(diǎn)稀疏，較難驗(yàn)證算法 的可靠性和反演精度。故國內(nèi)外適用于青藏高原地區(qū)的積雪微波遙感算法需加以改進(jìn)。
[0006] 目前，國內(nèi)大量與青藏相關(guān)的專題系統(tǒng)所生產(chǎn)的產(chǎn)品多依據(jù)系統(tǒng)的設(shè)計(jì)需求定 位，即產(chǎn)品進(jìn)行了特化處理，不適于作為研究該地區(qū)的獨(dú)立數(shù)據(jù)，因此當(dāng)前缺乏專門生產(chǎn)青 藏高原地區(qū)積雪研究所需具有普適性的基礎(chǔ)數(shù)據(jù)生產(chǎn)系統(tǒng)。當(dāng)今大范圍積雪監(jiān)測手段的主 要手段是通過衛(wèi)星遙感監(jiān)測方法實(shí)現(xiàn)，然而隨著衛(wèi)星遙感監(jiān)測方法的廣泛應(yīng)用，又產(chǎn)生了 的新的問題，即海量的長時間序列的遙感數(shù)據(jù)處理問題。同時，研究青藏高原地區(qū)積雪的時 空分布特征又需要長時間序列的積雪產(chǎn)品提供數(shù)據(jù)支持，結(jié)合積雪災(zāi)害預(yù)警預(yù)測相關(guān)研究 工作需要具有較高時效性的積雪產(chǎn)品。因此迫切需要結(jié)合高原區(qū)的積雪識別和反演過程算 法，開發(fā)一個可實(shí)現(xiàn)對海量遙感數(shù)據(jù)快速處理系統(tǒng)，提供長時間序列產(chǎn)品的業(yè)務(wù)化系統(tǒng)。
[0007] 綜上所述，現(xiàn)有技術(shù)中缺乏針對青藏高原地區(qū)的業(yè)務(wù)化積雪產(chǎn)品生產(chǎn)系統(tǒng)，因此， 亟待提供一種青藏高原雪水當(dāng)量估算方法及系統(tǒng)，在保證一定精度的情況，也可以實(shí)現(xiàn)對 長時間序列數(shù)據(jù)開展快速處理，為積雪的業(yè)務(wù)化和歷史數(shù)據(jù)的處理提供算法和工具。

【發(fā)明內(nèi)容】

[0008] 本發(fā)明的目的在于提供一種基于被動微波遙感的青藏高原雪水當(dāng)量估算方法及 系統(tǒng)，用于改善青藏高原雪水產(chǎn)品精度、實(shí)現(xiàn)對海量遙感數(shù)據(jù)的自動批量處理，供產(chǎn)時間序 列的相關(guān)積雪產(chǎn)品數(shù)據(jù)集，從而作為研究青藏高原地區(qū)積雪研究所需基礎(chǔ)數(shù)據(jù)，進(jìn)一步提 高現(xiàn)有青藏高原地區(qū)積雪監(jiān)測業(yè)務(wù)化水平。
[0009] 第一方面，本發(fā)明提供一種基于被動微波遙感的青藏高原雪水當(dāng)量估算方法，包 括步驟：
[0010] S1、從被動微波遙感數(shù)據(jù)中剔除包含各種干擾因子的像元，保留可能存在積雪像 元；
[0011] S2、通過閾值法將所述可能存在積雪像元分為中厚雪區(qū)域、淺雪區(qū)域以及無雪區(qū) 域；
[0012] S3、確定各地物下墊面的雪深經(jīng)驗(yàn)反演算法；
[0013] S4、對于所述中厚雪區(qū)域，檢索為某一地物時調(diào)用該地物下墊面的雪深經(jīng)驗(yàn)反演 算法計(jì)算雪水當(dāng)量，得到所述中厚雪區(qū)域的雪水當(dāng)量；
[0014] S5、設(shè)定所述淺雪區(qū)域的雪水當(dāng)量；
[0015] S6、根據(jù)得到的所述中厚雪區(qū)域的雪水當(dāng)量和設(shè)定的所述淺雪區(qū)域的雪水當(dāng)量， 得到青藏高原雪水當(dāng)量數(shù)據(jù)。
[0016] 優(yōu)選的，對于任意像元，所述步驟S1進(jìn)一步包括：
[0017] S11、判識該像元是否為水體;若否，則繼續(xù)步驟S12;
[0018] S12、判識該像元是否從氣候?qū)W角度經(jīng)驗(yàn)性地認(rèn)為不可能為積雪;若否，則繼續(xù)步 驟S13;
[0019] S13、判識該像元是否為冰川;若否，則繼續(xù)步驟S14;
[0020] S14、判識該像元是否為降雨區(qū);若否，則繼續(xù)步驟S15;
[0021] S15、判識該像元是否可能為積雪:若是，則判識為可能為積雪像元。
[0022]優(yōu)選的，所述步驟S15之后還包括：
[0023] S16.對于任意所述可能為積雪像元:通過高分辨率的輔助數(shù)據(jù)檢索到該像元中水 體大于百分之五十或者其他覆蓋大于百分之五十時，則剔除該像元，得到最終的可能存在 積雪像元。
[0024]優(yōu)選的，所述步驟S11-S15在第一分辨率下進(jìn)行，所述步驟S16在第二分辨率下進(jìn) 行，所述第二分辨率高于所述第一分辨率。
[0025] 優(yōu)選的，所述步驟S2采用標(biāo)準(zhǔn)的引入微波高頻的薄厚雪的判別方法、并通過閾值 法將所述可能存在積雪像元分為中厚雪區(qū)域、淺雪區(qū)域以及無雪區(qū)域，具體包括：
[0026] 依據(jù)增加高頻數(shù)據(jù)和中低頻率數(shù)據(jù)形成閾值的方式開展判別，使得積雪像元分為 中厚雪區(qū)域、淺雪區(qū)域以及無雪區(qū)域。
[0027]優(yōu)選的，所述步驟S3進(jìn)一步包括：
[0028]依據(jù)青藏高原觀測臺站實(shí)測數(shù)據(jù)結(jié)合不同地物下不同頻率亮溫組合進(jìn)行對比，提 取各地物相關(guān)性最好的亮溫組合，采用線性回歸方法，確定各地物下墊面的雪深經(jīng)驗(yàn)反演 算法。
[0029]優(yōu)選的，所述地物包括森林、灌木、草地以及裸土，所述各地物下墊面的雪深經(jīng)驗(yàn) 反演算法的系數(shù)如下：
[0030]森林下墊面的雪深經(jīng)驗(yàn)反演算法：
[0031 ] SD = 0 ? 023*tl036v*pol36*pol36+l ? 5;
[0032] 灌木下墊面的雪深經(jīng)驗(yàn)反演算法：
[0033] SD = -0.32519*tl036v*invlogpol36+9.16511；
[0034] 草地下墊面的雪深經(jīng)驗(yàn)反演算法：
[0035] SD = 0 ? 161*tl836v+0 ? 0516*pol89+0 ? 478;
[0036]裸土下墊面的雪深經(jīng)驗(yàn)反演算法：
[0037] SD = 0 ? 22*tl836v*invlogpol36*invlogpol36+0 ? 425;
[0038] 其中：tl036v表示亮溫梯度，是傳感器的10.7GHz垂直極化亮溫?cái)?shù)據(jù)與36.5GHz垂 直極化亮溫?cái)?shù)據(jù)的差值；
[0039] tl836v表示亮溫梯度，是傳感器的18.7GHz垂直極化亮溫?cái)?shù)據(jù)與36.5GHz垂直極化 亮溫?cái)?shù)據(jù)的差值；
[0040] P0136表示極化差，是傳感器的36.5GHz的垂直極化亮溫?cái)?shù)據(jù)與36GHz水平極化亮 溫?cái)?shù)據(jù)的差值；
[00411 pol89表示極化差，是傳感器的89GHz的垂直極化亮溫?cái)?shù)據(jù)與89GHz水平極化亮溫 數(shù)據(jù)的差值；
[0042] invlogpol36表示傳感器36.56他垂直極化亮溫(￥)和水平極化亮溫(田差值對數(shù) 的倒數(shù)；
[0043] SD代表計(jì)算的雪深，算法中的其他數(shù)據(jù)為相應(yīng)的系數(shù)。
[0044]優(yōu)選的，所述步驟S4可以進(jìn)一步包括：
[0045] S41.對于所述中厚雪像元區(qū)域:檢索為森林時，調(diào)用森林下墊面的雪深經(jīng)驗(yàn)反演 算法計(jì)算雪水當(dāng)量;檢索為灌木時，調(diào)用灌木下墊面的雪深經(jīng)驗(yàn)反演算法計(jì)算雪水當(dāng)量;檢 索為草原時，調(diào)用草原下雪深反演算法計(jì)算雪水當(dāng)量;檢索為裸土?xí)r，調(diào)用裸土下墊面的雪 深經(jīng)驗(yàn)反演算法計(jì)算雪水當(dāng)量；
[0046] S42.根據(jù)各地物下雪水當(dāng)量以及各地物在所述中厚雪區(qū)域中所占百分比計(jì)算所 述中厚雪區(qū)域的雪水當(dāng)量。
[0047] 第二方面，本發(fā)明還提供了實(shí)現(xiàn)上述任意一種基于被動微波遙感的青藏高原雪水 當(dāng)量估算方法的青藏高原雪水當(dāng)量估算系統(tǒng)，包括：
[0048] 像元獲取模塊，用于從被動微波遙感數(shù)據(jù)中剔除包含各種干擾因子的像元，保留 可能存在積雪像元；
[0049] 雪型劃分模塊，用于通過閾值法將所述可能存在積雪像元分為中厚雪區(qū)域、淺雪 區(qū)域以及無雪區(qū)域；
[0050] 算法生成模塊，用于確定各地物下墊面的雪深經(jīng)驗(yàn)反演算法；
[0051] 中厚雪區(qū)域雪水當(dāng)量計(jì)算模塊，用于對于所述中厚雪區(qū)域，檢索為某一地物時調(diào) 用該地物下墊面的雪深經(jīng)驗(yàn)反演算法計(jì)算雪水當(dāng)量，得到所述中厚雪區(qū)域的雪水當(dāng)量； [0052]淺雪區(qū)域雪水當(dāng)量設(shè)定模塊，用于設(shè)定所述淺雪區(qū)域的雪水當(dāng)量；
[0053]結(jié)果輸出模塊，用于根據(jù)得到的所述中厚雪區(qū)域的雪水當(dāng)量和設(shè)定的所述淺雪區(qū) 域的雪水當(dāng)量，得到青藏高原雪水當(dāng)量數(shù)據(jù)。
[0054] 優(yōu)選地，所述青藏高原雪水當(dāng)量估算系統(tǒng)采用基于可執(zhí)行算法包的集成和系統(tǒng)擴(kuò) 展策略。
[0055] 由上述技術(shù)方案可知，在本發(fā)明的基于被動微波遙感的青藏高原雪水當(dāng)量估算方 法及系統(tǒng)中，所采用的產(chǎn)品算法系數(shù)是基于青藏高原地區(qū)的實(shí)測數(shù)據(jù)開發(fā)獲得，產(chǎn)品質(zhì)量 控制決策樹依據(jù)青藏高原地區(qū)實(shí)際地物分布特征創(chuàng)建，形成符合青藏高原積雪的算法流 程，并完成數(shù)據(jù)產(chǎn)品的標(biāo)準(zhǔn)化輸出，并且上述估算系統(tǒng)采用模塊化設(shè)計(jì)，基于可執(zhí)行算法包 的集成和拓展方式，具有算法和產(chǎn)品的快速升級和更新的特點(diǎn)。因此，通過本發(fā)明所提供的 估算方法及系統(tǒng)所生產(chǎn)的產(chǎn)品精度較已有產(chǎn)品精度有大幅度提高;并且本發(fā)明實(shí)施例中所 提供的估算方法及系統(tǒng)具有產(chǎn)品生產(chǎn)速度快，具備近實(shí)時業(yè)務(wù)化處理能力，可完成長時間 序列積雪產(chǎn)品的生產(chǎn)，系統(tǒng)易于算法易于集成和拓展，系統(tǒng)生產(chǎn)的產(chǎn)品具有較高的精度，滿 足青藏高原地區(qū)積雪監(jiān)測業(yè)務(wù)化的要求，達(dá)到應(yīng)用要求。
【附圖說明】
[0056] 圖1是本發(fā)明實(shí)施例提供的基于被動微波遙感的青藏高原雪水當(dāng)量估算方法的流 程不意圖；
[0057] 圖2是本發(fā)明實(shí)施例提供的基于被動微波遙感的青藏高原雪水當(dāng)量估算方法中步 驟S1的流程示意圖；
[0058]圖3是本發(fā)明實(shí)施例提供的基于被動微波遙感的青藏高原雪水當(dāng)量估算方法中步 驟S2-S6的流程不意圖；
[0059] 圖4是本發(fā)明實(shí)施例中青藏高原雪水當(dāng)量估算系統(tǒng)的目錄結(jié)構(gòu)示意圖。
【具體實(shí)施方式】
[0060] 下面結(jié)合附圖和實(shí)施例，對本發(fā)明的【具體實(shí)施方式】做進(jìn)一步描述。以下實(shí)施例僅 用于說明本發(fā)明，但不用來限制本發(fā)明的范圍。
[0061] 圖1示出了本發(fā)明實(shí)施例提供的基于被動微波遙感的青藏高原雪水當(dāng)量估算方法 的流程示意圖，如圖1所示，本發(fā)明的基于被動微波遙感的青藏高原雪水當(dāng)量估算方法主要 包括以下步驟：
[0062] S1、從被動微波遙感數(shù)據(jù)中剔除包含各種干擾因子的像元，保留可能存在積雪像 元；
[0063] S2、通過閾值法將所述可能存在積雪像元分為中厚雪區(qū)域、淺雪區(qū)域以及無雪區(qū) 域；
[0064] S3、確定各地物下墊面的雪深經(jīng)驗(yàn)反演算法；
[0065] S4、對于所述中厚雪區(qū)域，檢索為某一地物時調(diào)用該地物下墊面的雪深經(jīng)驗(yàn)反演 算法計(jì)算雪水當(dāng)量，得到所述中厚雪區(qū)域的雪水當(dāng)量；
[0066] S5、設(shè)定所述淺雪區(qū)域的雪水當(dāng)量；
[0067] S6、根據(jù)得到的所述中厚雪區(qū)域的雪水當(dāng)量和設(shè)定的所述淺雪區(qū)域的雪水當(dāng)量， 得到青藏高原雪水當(dāng)量數(shù)據(jù)。
[0068] 下面對本發(fā)明的各個步驟加以詳細(xì)的說明。
[0069] S1.從被動微波遙感數(shù)據(jù)中剔除包含各種干擾因子的像元，保留可能存在積雪像 元，即基于青藏高原積雪分布特點(diǎn)設(shè)計(jì)產(chǎn)品質(zhì)量控制決策樹。
[0070] 水體、冰川和植被覆蓋以及降水的判別等因素都是影響反演雪深和雪水當(dāng)量估算 的干擾因子，在建立積雪算法前都需要進(jìn)行判別剔除，以免出現(xiàn)誤判而高估了雪深和雪水 當(dāng)量?，F(xiàn)有技術(shù)中，通常是利用積雪分類樹和采用輔助數(shù)據(jù)來剔除上述干擾因子。例如，K. G議ta利用土地覆蓋分類圖分析土地覆蓋對積雪累積和分布的影響；車濤利用1:4000000 中國湖泊分布圖判別受水體影響的像元，并利用高程數(shù)據(jù)將研究區(qū)分為四類分別進(jìn)行結(jié)果 驗(yàn)證；孫之文利用IGBP(International Geosphere Biosphere Program，國際地圈生物圈 計(jì)劃)分類圖消除大型水體的影響，獲取MOD 12Q1V004植被覆蓋數(shù)據(jù)和MOD IS NDVI產(chǎn)品調(diào)查 植被覆蓋情況等。
[0071] 本實(shí)施例中，使用青藏高原地物分類數(shù)據(jù)剔除這些干擾因子，實(shí)際處理中影響積 雪反演的一個重要因素是混合像元問題。由于通常采用的AMSR-E亮溫?cái)?shù)據(jù)的空間分辨率只 有25km，這樣一個像元就對應(yīng)一個25km*25km的區(qū)域，該區(qū)域可能包含多種地物，采用適合 單一地物的算法，反演得到的結(jié)果誤差很大，本實(shí)施例中的解決方是對這個區(qū)域采用青藏 高原地區(qū)高精度的地物分類數(shù)據(jù)進(jìn)行再次分類，針對分類的結(jié)果，不同的地物采用不同的 反演算法，從而達(dá)到提高該區(qū)域反演的精度。
[0072]例如:為了提高反演速度以及反演精度，本實(shí)施例中在制定質(zhì)量控制決策樹時，采 用了兩步分類策略：即采用兩種分辨率的地物分類數(shù)據(jù)，先采用粗分辨率的地物分類數(shù)據(jù) (第一分辨率，可優(yōu)選為25km)對像元進(jìn)行初步分類，對判斷是可能存在積雪像元采用青藏 高原地區(qū)高精度的地物分類數(shù)據(jù)(第二分辨率，可優(yōu)選為〇.5km)進(jìn)行再次細(xì)分。通過實(shí)際計(jì) 算產(chǎn)品以及測試結(jié)果表明，兩步分類的方法具有執(zhí)行速度快，分類準(zhǔn)確的特點(diǎn)。
[0073]如圖2中所示，本實(shí)施例中的步驟S1具體步驟如下：
[0074] 首先，利用積雪判識文件(第一分辨率），對產(chǎn)品像元進(jìn)行逐一判別，剔除含干擾因 子的像元，僅對可能有積雪存在的像元進(jìn)行積雪存在判斷，若有散射體存在，則進(jìn)行積雪計(jì) 算操作。在積雪判識中，利用判識文件首先判識該像元是否為水體，若是水體，則標(biāo)識為"水 體"，并給一個標(biāo)識值"254"，若不是水體，則繼續(xù)判識;然后判識該像元是否從氣候?qū)W角度 經(jīng)驗(yàn)性地認(rèn)為不可能為積雪，若是，則標(biāo)識為"不可能為積雪"，并賦標(biāo)識值"252"，若不是， 則繼續(xù)判識;接著判識該像元是否為冰川，若是，則標(biāo)識為"冰川"，并賦標(biāo)識值"253"，若不 是，則繼續(xù)判識;接著判識該像元是否為降雨區(qū)，若是，則標(biāo)識為"降雨區(qū)"，若不是，則繼續(xù) 判識;最后判識該像元是否可能為積雪:若判識該像元為可能為積雪，則標(biāo)識為"可能為積 雪"；
[0075] 其次，利用高分辨率（第二分辨率）的青藏高原地區(qū)地物分類數(shù)據(jù)，對"可能為積 雪"的像元進(jìn)行再次分類。例如，當(dāng)檢索到該像元水體>50%時，算法標(biāo)記為"水體"，賦值為 "254"，不進(jìn)行后續(xù)的積雪反演;若檢索為其他覆蓋>50%時，則也不進(jìn)行后續(xù)的積雪反演， 并標(biāo)記為"其他"，賦值為"249"。上述提到的高分辨率（分辨率0.5km)的地分類數(shù)據(jù)采用 M0DIS 土地覆蓋數(shù)據(jù)，土地覆蓋類型按照國際陸界生物圈方案（IGBP)劃分，永久濕地、雪和 冰(指的是山區(qū)永久性積雪和冰川）、未分類數(shù)據(jù)在本發(fā)明中定義為"其他"；
[0076] 最后，被保留的像元即最終的可能存在積雪像元。
[0077]需要說明的是，本實(shí)施例中的第二分辨率高于第一分辨率，本實(shí)施例中的粗分辨 率(即第一分辨率)和高分辨率(即第二分辨率)的取值僅僅是根據(jù)原始數(shù)據(jù)而言;在實(shí)際應(yīng) 用中，當(dāng)然也可以是其他分辨率。
[0078] 步驟S2-步驟S6如圖3中所示：
[0079] S2、通過閾值法將所述可能存在積雪像元分為中厚雪區(qū)域、淺雪區(qū)域以及無雪區(qū) 域。
[0080] 在具體應(yīng)用中，所述步驟S2采用標(biāo)準(zhǔn)的引入微波高頻的薄厚雪的判別方法、并通 過閾值法將所述可能存在積雪像元分為中厚雪區(qū)域、淺雪區(qū)域以及無雪區(qū)域，具體包括： [0081 ]依據(jù)增加高頻數(shù)據(jù)和中低頻率數(shù)據(jù)形成閾值的方式開展判別，使得積雪像元分為 中厚雪區(qū)域、淺雪區(qū)域以及無雪區(qū)域。
[0082] S3、確定各地物下墊面的雪深經(jīng)驗(yàn)反演算法；即依據(jù)青藏高原地區(qū)實(shí)測數(shù)據(jù)開發(fā) 的算法系數(shù)。例如，本實(shí)施例中，所述地物包括森林、灌木、草地以及裸土;依據(jù)青藏高原觀 測臺站實(shí)測數(shù)據(jù)結(jié)合不同地物下不同頻率亮溫組合進(jìn)行對比，提取各地物相關(guān)性最好的亮 溫組合，采用線性回歸方法，確定各地物下墊面的雪深經(jīng)驗(yàn)反演算法。最終確定的各地物下 墊面的雪深經(jīng)驗(yàn)反演算法系統(tǒng)如下：
[0083] 森林下墊面的雪深經(jīng)驗(yàn)反演算法：
[0084] SD = 0 ? 023*tl036v*pol36*pol36+l ? 5;
[0085]灌木下墊面的雪深經(jīng)驗(yàn)反演算法：
[0086] SD = -0.32519*tl036v*invlogpol36+9.16511；
[0087] 草地下墊面的雪深經(jīng)驗(yàn)反演算法：
[0088] SD = 0 ? 161*tl836v+0 ? 0516*pol89+0 ? 478;
[0089]裸土下墊面的雪深經(jīng)驗(yàn)反演算法：
[0090] SD = 0 ? 22*tl836v*invlogpol36*invlogpol36+0 ? 425;
[0091] 其中：
[0092] 1 )tl036v表示亮溫梯度，是傳感器的10.7GHz垂直極化亮溫?cái)?shù)據(jù)與36.5GHz垂直極 化亮溫?cái)?shù)據(jù)的差值。
[0093] 2) tl836v表示亮溫梯度，是傳感器的18.7GHz垂直極化亮溫?cái)?shù)據(jù)與36.5GHz垂直極 化亮溫?cái)?shù)據(jù)的差值。
[0094] 3)p〇136表示極化差，是傳感器的36.5GHz的垂直極化亮溫?cái)?shù)據(jù)與36GHz水平極化 亮溫?cái)?shù)據(jù)的差值。
[0095] 4)p〇189表示極化差，是傳感器的89GHz的垂直極化亮溫?cái)?shù)據(jù)與89GHz水平極化亮 溫?cái)?shù)據(jù)的差值。
[0096] 5)invlogpol36表示傳感器36.56他垂直極化亮溫(￥)和水平極化亮溫(田差值對 數(shù)的倒數(shù)，等于 1 ./(l〇gl〇(P〇136))。
[0097] 6) SD代表計(jì)算的雪深，算法中的其他數(shù)據(jù)為相應(yīng)的系數(shù)。
[0098] S4.對于所述中厚雪區(qū)域，檢索為某一地物時調(diào)用該地物下墊面的雪深經(jīng)驗(yàn)反演 算法計(jì)算雪水當(dāng)量，得到所述中厚雪區(qū)域的雪水當(dāng)量;例如，對于所述中厚雪像元區(qū)域:檢 索為森林時，調(diào)用森林下墊面的雪深經(jīng)驗(yàn)反演算法計(jì)算雪水當(dāng)量;檢索為灌木時，調(diào)用灌木 下墊面的雪深經(jīng)驗(yàn)反演算法計(jì)算雪水當(dāng)量;檢索為草原時，調(diào)用草原下雪深反演算法計(jì)算 雪水當(dāng)量;檢索為裸土?xí)r，調(diào)用裸土下墊面的雪深經(jīng)驗(yàn)反演算法計(jì)算雪水當(dāng)量；
[0099] 根據(jù)各地物下雪水當(dāng)量以及各地物在所述中厚雪區(qū)域中所占百分比計(jì)算所述中 厚雪區(qū)域的雪水當(dāng)量。
[0100] 具體的，采用上述四種地物下墊面的雪深經(jīng)驗(yàn)反演算法求取中厚雪區(qū)域雪水當(dāng)量 數(shù)據(jù)說明具體如表1中所示。
[0101] 衷1青藏高原中厚雪岡域各下墊面的雪水當(dāng)量筧法
[0103] 上表中，SWE'為中厚雪區(qū)域總的雪水當(dāng)量，單位為mm;Ff、Fs、Fg、Fb分別代表森林、 灌木、草原、裸地覆蓋類型在中厚雪區(qū)域內(nèi)所占的百分比；SD f、SDs、SDg、SDb分別表示森林、 灌木、草原、裸地覆蓋類型的雪深反演算法。
[0104] ps是積雪雪密度，從該算法所需的輔助數(shù)據(jù)中得到。
[0105] S5、設(shè)定所述淺雪區(qū)域的雪水當(dāng)量SWE^。
[0106] 可理解的是，由于所述淺雪區(qū)域的被動微波普遍難于監(jiān)測，可將所述淺雪區(qū)域積 雪的雪深估算設(shè)定為5cm，SWE淺=5*密度。
[0107] S6、根據(jù)得到的所述中厚雪區(qū)域的雪水當(dāng)量和設(shè)定的所述淺雪區(qū)域的雪水當(dāng)量， 得到青藏高原雪水當(dāng)量數(shù)據(jù)，BP:
[0108] SWE=Ff X SWEf+Fs X SWEs+Fg X SffEg+Fb X SWEb 〇
[0109] 本實(shí)施例的基于被動微波遙感的青藏高原雪水當(dāng)量估算方法，采用的產(chǎn)品算法系 數(shù)是基于青藏高原地區(qū)的實(shí)測數(shù)據(jù)開發(fā)獲得，產(chǎn)品質(zhì)量控制決策樹依據(jù)青藏高原地區(qū)實(shí)際 地物分布特征創(chuàng)建，形成符合青藏高原積雪的算法流程，并完成數(shù)據(jù)產(chǎn)品的標(biāo)準(zhǔn)化輸出，通 過本實(shí)施例所提供的估算方法所生產(chǎn)的產(chǎn)品精度較已有產(chǎn)品精度有大幅度提高;并且本發(fā) 明實(shí)施例中所提供的估算方法具有產(chǎn)品生產(chǎn)速度快，具備近實(shí)時業(yè)務(wù)化處理能力，可完成 長時間序列積雪產(chǎn)品的生產(chǎn)，算法易于集成和拓展，生產(chǎn)的產(chǎn)品具有較高的精度，滿足青藏 高原地區(qū)積雪監(jiān)測業(yè)務(wù)化的要求，達(dá)到應(yīng)用要求。
[0110] 本發(fā)明還提供了實(shí)現(xiàn)圖1所示實(shí)施例所述的基于被動微波遙感的青藏高原雪水當(dāng) 量估算方法的青藏高原雪水當(dāng)量估算系統(tǒng)，包括：
[0111] 像元獲取模塊，用于從被動微波遙感數(shù)據(jù)中剔除包含各種干擾因子的像元，保留 可能存在積雪像元；
[0112] 雪型劃分模塊，用于通過閾值法將所述可能存在積雪像元分為中厚雪區(qū)域、淺雪 區(qū)域以及無雪區(qū)域；
[0113] 算法生成模塊，用于確定各地物下墊面的雪深經(jīng)驗(yàn)反演算法；
[0114] 中厚雪區(qū)域雪水當(dāng)量計(jì)算模塊，用于對于所述中厚雪區(qū)域，檢索為某一地物時調(diào) 用該地物下墊面的雪深經(jīng)驗(yàn)反演算法計(jì)算雪水當(dāng)量，得到所述中厚雪區(qū)域的雪水當(dāng)量；
[0115] 淺雪區(qū)域雪水當(dāng)量設(shè)定模塊，用于設(shè)定所述淺雪區(qū)域的雪水當(dāng)量；
[0116]結(jié)果輸出模塊，用于根據(jù)得到的所述中厚雪區(qū)域的雪水當(dāng)量和設(shè)定的所述淺雪區(qū) 域的雪水當(dāng)量，得到青藏高原雪水當(dāng)量數(shù)據(jù)。
[0117] 進(jìn)一步的，為了建立靈活，簡易的系統(tǒng)算法更新機(jī)制和方便系統(tǒng)拓展，所述青藏高 原雪水當(dāng)量估算系統(tǒng)可采用基于可執(zhí)行算法包的集成和系統(tǒng)擴(kuò)展策略。例如：
[0118] 建立如圖4中所示系統(tǒng)目錄結(jié)構(gòu)，作為青藏高原雪水當(dāng)量估算系統(tǒng)的目錄結(jié)構(gòu)。系 統(tǒng)目錄設(shè)計(jì)方案為：系統(tǒng)主文件是名為"tibet"文件夾，下面是二級目錄，"prouducts"（產(chǎn) 品）存放各產(chǎn)品生產(chǎn)子系統(tǒng)生產(chǎn)的產(chǎn)品；"waring"存放批處理生成的警告日志，"pr 〇dUCt_ system"（產(chǎn)品生產(chǎn)系統(tǒng)）用于存放產(chǎn)品生產(chǎn)子系統(tǒng)。"TibetSnow"可執(zhí)行腳本文件，功能是 調(diào)度整個生產(chǎn)系統(tǒng)；"bin"目錄下的"AMSR-E""算法腳本實(shí)現(xiàn)對不同算法的調(diào)用。三四級目 錄作用來對產(chǎn)品生產(chǎn)子系統(tǒng)進(jìn)行分類，形如"AMSR-E/SWE"分類原則是生產(chǎn)用原始數(shù)據(jù)和產(chǎn) 品類型分類。
[0119] 目錄說明見表2。
[0120]表2系統(tǒng)目錄結(jié)構(gòu)說明表
[0122] 系統(tǒng)添加的新算法子系統(tǒng)策略：
[0123] 首先建立一個與發(fā)布系統(tǒng)相同的目錄結(jié)構(gòu)，在"SWE"目錄下建立一個名為"src"的 目錄，編寫新算法的代碼，新算法所需的輔助數(shù)據(jù)保存在"anc"文件夾，生產(chǎn)日志保存在 "log"文件夾下，系統(tǒng)生產(chǎn)產(chǎn)品時中間計(jì)算數(shù)據(jù)臨時數(shù)據(jù)保存在"tmp"文件夾，代碼編譯生 成的可執(zhí)行程序保存在"bin"文件夾。新算法編譯完成，測試可以發(fā)布之后，將"bin"文件夾 中的可執(zhí)行程序復(fù)制到發(fā)布系統(tǒng)對應(yīng)"bin"文件夾下，新建文件"SWE##"內(nèi)。系統(tǒng)默認(rèn)設(shè)計(jì) 了通用的接口，這些接口包括:產(chǎn)品生產(chǎn)模式接口，產(chǎn)品投影方式選擇接口，產(chǎn)品算法選擇 接口，輸入數(shù)據(jù)路徑，輸出數(shù)據(jù)路徑。當(dāng)新添加算法采用默認(rèn)接口時，此時系統(tǒng)即完成了新 算法的添加工作。當(dāng)新算法有新的接口需求時，可在"AMSR-E"算法腳本和"TibetSnow"系統(tǒng) 總調(diào)度腳本中增加部分語句，靈活的定制新算法系統(tǒng)所需接口，實(shí)現(xiàn)系統(tǒng)算法的添加。由于 和發(fā)布系統(tǒng)采用相同的目錄結(jié)構(gòu)，所以只需在發(fā)布系統(tǒng)中增加相應(yīng)的調(diào)用語句即可。新算 法的編譯獨(dú)立與發(fā)布系統(tǒng)，各個算法之間相互獨(dú)立，互不影響。
[0124] 系統(tǒng)拓展策略：
[0125] 首先復(fù)制一個發(fā)布系統(tǒng)，根據(jù)系統(tǒng)拓展需求建立相應(yīng)的子目錄，如圖4所示帶填充 顯示文件夾部分，拓展部分算法添加更新部分參照系統(tǒng)添加的新算法子系統(tǒng)策略，在復(fù)制 發(fā)布系統(tǒng)的"TibetSnow"生產(chǎn)系統(tǒng)腳本中增加相應(yīng)調(diào)用語句，編譯完成以后，在發(fā)布系統(tǒng)中 增加相關(guān)目錄，替換相應(yīng)的腳本文件，即可實(shí)現(xiàn)系統(tǒng)的拓展。
[0126] 上述青藏高原雪水當(dāng)量估算系統(tǒng)制定了基于可執(zhí)行程序的集成模式，系統(tǒng)算法的 添加和拓展方便靈活，易于實(shí)現(xiàn);并且，所采用的產(chǎn)品算法系數(shù)是基于青藏高原地區(qū)的實(shí)測 數(shù)據(jù)開發(fā)獲得，系統(tǒng)產(chǎn)品質(zhì)量控制決策樹依據(jù)青藏高原地區(qū)實(shí)際地物分布特征創(chuàng)建，系統(tǒng) 所生廣廣品精度較已有廣品精度有大幅度提尚。
[0127] 本發(fā)明所提供的青藏高原雪水當(dāng)量估算系統(tǒng)具有產(chǎn)品生產(chǎn)速度快，具備近實(shí)時業(yè) 務(wù)化處理能力，可完成長時間序列積雪產(chǎn)品的生產(chǎn)，系統(tǒng)易于算法易于集成和拓展，系統(tǒng)生 產(chǎn)的產(chǎn)品具有較高的精度，滿足青藏高原地區(qū)積雪監(jiān)測業(yè)務(wù)化的要求，達(dá)到應(yīng)用要求。
[0128] 以上實(shí)施方式僅用于說明本發(fā)明，而并非對本發(fā)明的限制，有關(guān)技術(shù)領(lǐng)域的普通 技術(shù)人員，在不脫離本發(fā)明的精神和范圍的情況下，還可以做出各種變化和變型，因此所有 等同的技術(shù)方案也屬于本發(fā)明的保護(hù)范疇。
【主權(quán)項(xiàng)】
1. 一種基于被動微波遙感的青藏高原雪水當(dāng)量估算方法，其特征在于，包括步驟： 51、 從被動微波遙感數(shù)據(jù)中剔除包含各種干擾因子的像元，保留可能存在積雪像元； 52、 通過閾值法將所述可能存在積雪像元分為中厚雪區(qū)域、淺雪區(qū)域以及無雪區(qū)域； 53、 確定各地物下墊面的雪深經(jīng)驗(yàn)反演算法； 54、 對于所述中厚雪區(qū)域，檢索為某一地物時調(diào)用該地物下墊面的雪深經(jīng)驗(yàn)反演算法 計(jì)算雪水當(dāng)量，得到所述中厚雪區(qū)域的雪水當(dāng)量； 55、 設(shè)定所述淺雪區(qū)域的雪水當(dāng)量； 56、 根據(jù)得到的所述中厚雪區(qū)域的雪水當(dāng)量和設(shè)定的所述淺雪區(qū)域的雪水當(dāng)量，得到 青藏高原雪水當(dāng)量數(shù)據(jù)。2. 根據(jù)權(quán)利要求1所述的方法，其特征在于，對于任意像元，所述步驟Sl進(jìn)一步包括： SI 1、判識該像元是否為水體;若否，則繼續(xù)步驟S12; 512、 判識該像元是否從氣候?qū)W角度經(jīng)驗(yàn)性地認(rèn)為不可能為積雪；若否，則繼續(xù)步驟 S13； 513、 判識該像元是否為冰川;若否，則繼續(xù)步驟S14; 514、 判識該像元是否為降雨區(qū);若否，則繼續(xù)步驟S15; 515、 判識該像元是否可能為積雪:若是，則判識為可能為積雪像元。3. 根據(jù)權(quán)利要求2所述的方法，其特征在于，所述步驟S15之后還包括： 516、 對于任意所述可能為積雪像元:通過高分辨率的輔助數(shù)據(jù)檢索到該像元中水體大 于百分之五十或者其他覆蓋大于百分之五十時，則剔除該像元，得到最終的可能存在積雪 像元。4. 根據(jù)權(quán)利要求3所述的方法，其特征在于，所述步驟S11-S15在第一分辨率下進(jìn)行，所 述步驟S16在第二分辨率下進(jìn)行，所述第二分辨率高于所述第一分辨率。5. 根據(jù)權(quán)利要求1所述的方法，其特征在于，所述步驟S2采用標(biāo)準(zhǔn)的引入微波高頻的薄 厚雪的判別方法、并通過閾值法將所述可能存在積雪像元分為中厚雪區(qū)域、淺雪區(qū)域以及 無雪區(qū)域，具體包括： 依據(jù)增加高頻數(shù)據(jù)和中低頻率數(shù)據(jù)形成閾值的方式開展判別，使得積雪像元分為中厚 雪區(qū)域、淺雪區(qū)域以及無雪區(qū)域。6. 根據(jù)權(quán)利要求1任一項(xiàng)所述的方法，其特征在于，所述步驟S3進(jìn)一步包括： 依據(jù)青藏高原觀測臺站實(shí)測數(shù)據(jù)結(jié)合不同地物下不同頻率亮溫組合進(jìn)行對比，提取各 地物相關(guān)性最好的亮溫組合，采用線性回歸方法，確定各地物下墊面的雪深經(jīng)驗(yàn)反演算法。7. 根據(jù)權(quán)利要求6所述的方法，其特征在于，所述地物包括森林、灌木、草地以及裸土， 所述各地物下墊面的雪深經(jīng)驗(yàn)反演算法的系數(shù)如下； 森林下墊面的雪深經(jīng)驗(yàn)反演算法： SD = O .023*tl036v*pol 36*po 136+1 · 5; 灌木下墊面的雪深經(jīng)驗(yàn)反演算法： SD = -0.32519*tl036v*invlogpol36+9.16511; 草地下墊面的雪深經(jīng)驗(yàn)反演算法： SD = O. 161*tl836v+0.0516*pol89+0.478; 裸土下墊面的雪深經(jīng)驗(yàn)反演算法： SD = O. 22*tl836v*invlogpol36*invlogpol36+0.425; 其中：tl036v表示亮溫梯度，是傳感器的10.7GHz垂直極化亮溫?cái)?shù)據(jù)與36.5GHz垂直極 化亮溫?cái)?shù)據(jù)的差值； tl836v表示亮溫梯度，是傳感器的18.7GHz垂直極化亮溫?cái)?shù)據(jù)與36.5GHz垂直極化亮溫 數(shù)據(jù)的差值； p〇136表示極化差，是傳感器的36.5GHz的垂直極化亮溫?cái)?shù)據(jù)與36GHz水平極化亮溫?cái)?shù) 據(jù)的差值； p〇189表示極化差，是傳感器的89GHz的垂直極化亮溫?cái)?shù)據(jù)與89GHz水平極化亮溫?cái)?shù)據(jù) 的差值； invl〇gp〇136表示傳感器36.5GHz垂直極化亮溫(V)和水平極化亮溫(H)差值對數(shù)的倒 數(shù)； SD代表計(jì)算的雪深，算法中的其他數(shù)據(jù)為相應(yīng)的系數(shù)。8. 根據(jù)權(quán)利要求1所述的方法，其特征在于，所述步驟S4包括：541. 對于所述中厚雪像元區(qū)域:檢索為森林時，調(diào)用森林下墊面的雪深經(jīng)驗(yàn)反演算法 計(jì)算雪水當(dāng)量;檢索為灌木時，調(diào)用灌木下墊面的雪深經(jīng)驗(yàn)反演算法計(jì)算雪水當(dāng)量;檢索為 草原時，調(diào)用草原下雪深反演算法計(jì)算雪水當(dāng)量;檢索為裸土?xí)r，調(diào)用裸土下墊面的雪深經(jīng) 驗(yàn)反演算法計(jì)算雪水當(dāng)量；542. 根據(jù)各地物下雪水當(dāng)量以及各地物在所述中厚雪區(qū)域中所占百分比計(jì)算所述中 厚雪區(qū)域的雪水當(dāng)量。9. 一種實(shí)現(xiàn)權(quán)利要求1-8中任一項(xiàng)所述的基于被動微波遙感的青藏高原雪水當(dāng)量估算 方法的青藏高原雪水當(dāng)量估算系統(tǒng)，其特征在于，包括： 像元獲取模塊，用于從被動微波遙感數(shù)據(jù)中剔除包含各種干擾因子的像元，保留可能 存在積雪像元； 雪型劃分模塊，用于通過閾值法將所述可能存在積雪像元分為中厚雪區(qū)域、淺雪區(qū)域 以及無雪區(qū)域； 算法生成模塊，用于確定各地物下墊面的雪深經(jīng)驗(yàn)反演算法； 中厚雪區(qū)域雪水當(dāng)量計(jì)算模塊，用于對于所述中厚雪區(qū)域，檢索為某一地物時調(diào)用該 地物下墊面的雪深經(jīng)驗(yàn)反演算法計(jì)算雪水當(dāng)量，得到所述中厚雪區(qū)域的雪水當(dāng)量； 淺雪區(qū)域雪水當(dāng)量設(shè)定模塊，用于設(shè)定所述淺雪區(qū)域的雪水當(dāng)量； 結(jié)果輸出模塊，用于根據(jù)得到的所述中厚雪區(qū)域的雪水當(dāng)量和設(shè)定的所述淺雪區(qū)域的 雪水當(dāng)量，得到青藏高原雪水當(dāng)量數(shù)據(jù)。10. 根據(jù)權(quán)利要求9所述的系統(tǒng)，其特征在于，所述青藏高原雪水當(dāng)量估算系統(tǒng)采用基 于可執(zhí)行算法包的集成和系統(tǒng)擴(kuò)展策略。
【文檔編號】G06F19/00GK105893744SQ201610188859
【公開日】2016年8月24日
【申請日】2016年3月29日
【發(fā)明人】邱玉寶
【申請人】中國科學(xué)院遙感與數(shù)字地球研究所