熊文兵， 莫偉華， 丁美花， 王君華， 莫建飛

(1．廣西來賓市氣象局，來賓 546100； 2.廣西氣象科技服務(wù)中心，南寧 530022；3.廣西氣象減災(zāi)研究所，南寧 530022)

0 引言

地理信息系統(tǒng)(geographic information system，GIS)及衛(wèi)星遙感技術(shù)在廣西氣象行業(yè)的應(yīng)用研究已經(jīng)開展多年，成果豐厚。近年來，廣西利用遙感資料開展了廣西生態(tài)質(zhì)量氣象評價方面的研究，取得了很好的效果。經(jīng)過逐步的試驗研究，在生態(tài)質(zhì)量氣象評價方面的成果已經(jīng)投入業(yè)務(wù)應(yīng)用[1]。但是，各種評價業(yè)務(wù)都需要完成多種指數(shù)模型的計算，部分指數(shù)需要應(yīng)用GIS計算，圖形產(chǎn)品也需要在ENVI和ArcGIS等平臺上來處理、設(shè)計和完成?？傊?，一份產(chǎn)品的制作過程要經(jīng)過繁瑣的手工操作才能完成，產(chǎn)品風格由工作人員臨時設(shè)計，業(yè)務(wù)流程沒有標準化。針對以上問題和需求，本研究利用GIS技術(shù)，設(shè)計開發(fā)了廣西生態(tài)質(zhì)量氣象評價業(yè)務(wù)系統(tǒng)，讓該項業(yè)務(wù)工作人員能夠方便地利用該系統(tǒng)對氣象觀測數(shù)據(jù)及MODIS遙感數(shù)據(jù)進行處理，對廣西生態(tài)質(zhì)量氣象評價指標中的濕潤指數(shù)、植被覆蓋指數(shù)、水體密度指數(shù)、土地退化指數(shù)及災(zāi)害指數(shù)5個指標信息進行提取與計算，并利用GIS生成業(yè)務(wù)產(chǎn)品，實現(xiàn)業(yè)務(wù)工作的流程化與標準化。

GIS是一種采集、存儲、管理、分析、顯示與應(yīng)用地理信息的計算機技術(shù)系統(tǒng)，是分析和處理海量地理數(shù)據(jù)的通用技術(shù)。所有具有空間分布的數(shù)據(jù)都可以使用GIS進行處理，在國土、林業(yè)、農(nóng)業(yè)、水利、測繪、交通、生態(tài)、規(guī)劃和物流等很多領(lǐng)域有著廣泛的應(yīng)用[2-7]。近年來GIS技術(shù)在我國氣象上的應(yīng)用迅速發(fā)展，方法日趨成熟，成果喜人。如武漢市利用GIS技術(shù)建立了長江中游氣象水文預(yù)報系統(tǒng)[8]； 湖南省利用GIS技術(shù)建立了山洪災(zāi)害預(yù)警系統(tǒng)[9]； 沙莎等[10]研發(fā)了“氣象災(zāi)害評估系統(tǒng)(干旱模塊)”； 吳盛洪[11]研發(fā)了“氣象服務(wù)監(jiān)控管理平臺”； 董學士等[12]設(shè)計開發(fā)了“氣象查詢服務(wù)系統(tǒng)”； 李軒等[13]開發(fā)了“農(nóng)作物病蟲害氣象等級預(yù)報系統(tǒng)”； 熊文兵等[14]與施益強等[15]分別在基于多普勒雷達數(shù)據(jù)可視化的基礎(chǔ)上研發(fā)了預(yù)警系統(tǒng)； 高大偉等[16]則研究了風云氣象衛(wèi)星數(shù)據(jù)的應(yīng)用。

常見GIS應(yīng)用開發(fā)使用的地圖組件有很多種，功能和效果各有不同。目前免費的開源GIS平臺越來越受到開發(fā)人員的青睞。開源GIS組件具有小巧方便、效率高和靈活性強等特點，特別適合輕量級業(yè)務(wù)應(yīng)用方面的開發(fā)。本研究選擇的開源GIS——MapWinGIS是一個高效的地圖引擎，放大、縮小、漫游速度很快，開發(fā)人員對其評價較高。MapWinGIS是一套二次開發(fā)組件庫，其核心庫是一個名字為MapWinGIS的ActiveX控件，開發(fā)人員可以利用這個組件在已有系統(tǒng)中添加GIS相關(guān)功能，其二次開發(fā)環(huán)境支持VB和C#等語言。

1 系統(tǒng)功能設(shè)計

系統(tǒng)主體框架包含3部分： GIS功能模塊、氣候條件分析模塊和生態(tài)質(zhì)量指數(shù)計算模塊。GIS功能模塊主要包括地圖管理、圖層管理、地圖縮放和漫游等； 氣候條件分析模塊包括月份與季度氣候資料的累加計算、距平計算和數(shù)據(jù)分析功能，主要任務(wù)是生成平均氣溫分布圖、平均氣溫距平分布圖、降水量分布圖、降水量距平百分率分布圖、日照時數(shù)分布圖和日照時數(shù)距平分布圖等圖件產(chǎn)品； 生態(tài)質(zhì)量指數(shù)計算模塊實現(xiàn)濕潤指數(shù)、植被覆蓋指數(shù)、水體密度指數(shù)、土地退化指數(shù)和災(zāi)害指數(shù)等的計算及評價分析。

系統(tǒng)功能布局設(shè)計為常規(guī)的帶菜單欄的模式窗體(圖1)。布局設(shè)計包括功能菜單欄+工具條、地圖框+操作欄選項卡。地圖圖層操作設(shè)計為常用的樹目錄結(jié)構(gòu)。主功能菜單包括文件、氣候條件分析、評價指標、地理信息、窗口和幫助等(圖2)。

圖1 系統(tǒng)運行界面

圖2 系統(tǒng)功能菜單結(jié)構(gòu)

2 主要功能實現(xiàn)

系統(tǒng)的功能是依據(jù)業(yè)務(wù)流程進行設(shè)計的，功能的編排也按照業(yè)務(wù)操作順序進行。按照功能屬性的不同可以分為數(shù)據(jù)分析處理功能、GIS操作功能和產(chǎn)品制作功能3個類別。

2.1 數(shù)據(jù)分析處理功能

數(shù)據(jù)分析處理主要是為了得到分析結(jié)果的表格或圖形等業(yè)務(wù)產(chǎn)品而進行的初步工作，包括氣候條件數(shù)據(jù)即氣溫、降雨和日照等的網(wǎng)格化分析，柵格距平圖的計算及MODIS衛(wèi)星數(shù)據(jù)的計算處理。

2.1.1 氣候條件數(shù)據(jù)網(wǎng)格化

廣西各氣象臺站的常規(guī)觀測資料(氣溫、降雨和日照等)僅代表該站點所在地的氣象要素值，不能全面真實反映廣西全區(qū)范圍內(nèi)氣候資源的整體分布特征，為比較客觀真實地反映廣西不同區(qū)域的氣候資源狀況，需建立氣候要素隨地理因子變化的空間分析模型，即

Y=f(φ,λ,h,β,θ)+ε，

(1)

式中：Y為氣候要素(如氣溫、降雨和日照時數(shù)等)；φ,λ,h,β,θ分別為緯度、經(jīng)度、海拔高度、坡向、坡度等地理因子；ε為余差項，稱為綜合地理殘差，可視為φ,λ,h,β,θ所擬合的氣候?qū)W方程的殘差部分。若用Y*表示臺站實測值，則

ε=Y*-f(φ,λ,h,β,θ) 。

(2)

將廣西各臺站(92個常規(guī)站及1 500多個自動站)所測的氣溫、降雨及日照時數(shù)等氣候要素與各站點所在地的緯度、經(jīng)度、海拔高度、坡向和坡度等地理因子采用回歸分析方法建立各氣候要素的空間分析模型f(φ,λ,h,β,θ)。

根據(jù)各氣候要素的回歸分析模型，利用事先準備好的1 km×1 km網(wǎng)格上的φ,λ,h,β,θ等地理數(shù)據(jù)，推算出各氣候要素在1 km×1 km網(wǎng)格上的分布值，再進行殘差訂正。

殘差訂正方法是： 首先，將所有氣象站點的φ,λ,h,β,θ等地理數(shù)據(jù)分別代入各回歸方程，得到各氣候要素在站點的推算值； 再根據(jù)式(2)，用各站點的實測值與推算值得到殘差項ε； 運用反距離權(quán)重法，將ε插值到1 km×1 km網(wǎng)格點上，生成氣候要素的殘差柵格圖； 最后將該柵格圖與前面計算好的氣候資源分布圖進行疊加運算，即得到各氣候要素的實際分布圖。

本系統(tǒng)為了方便業(yè)務(wù)使用，設(shè)計了分析向?qū)В瑢⒎治鲞^程流程化。向?qū)Х?步：

1)選擇氣候資料，建立擬合方程。在使用對話框?qū)我坏臍夂驍?shù)據(jù)文件導入后，系統(tǒng)自動提供要素選擇項供用戶選擇作業(yè)的要素(如氣溫)，在執(zhí)行線性回歸方法后，很快得出擬合方程。系統(tǒng)專門編寫了線性回歸分析的類LineRegress(圖3)。

圖3 系統(tǒng)開發(fā)用到的組件(箭頭連接線表示類間的被調(diào)用邏輯關(guān)系)

2)數(shù)據(jù)網(wǎng)格化。根據(jù)擬合方程，可以得到網(wǎng)格化數(shù)據(jù)，并按照反距離權(quán)重法進行殘差訂正，系統(tǒng)提供選擇項，可以選擇訂正與否，并提供反距離權(quán)重的參數(shù)設(shè)置等。系統(tǒng)專門編寫了反距離權(quán)重插值的類IDW(圖3)。

3)出圖。以上計算結(jié)果以TIF格式輸出成柵格圖，此步驟提供選擇輸出結(jié)果導入系統(tǒng)地圖顯示界面，將TIF文件作為新圖層添加到地圖。此步驟可以獨立使用。

2.1.2 距平計算

為了得到氣候條件分析的距平圖，系統(tǒng)需要提供距平計算功能。為便于使用，專門開發(fā)了可以單獨使用的距平計算子模塊——柵格圖相減，即將歷史數(shù)據(jù)分布預(yù)先計算，也保存為柵格圖。計算距平時利用MapWinGIS的Grid接口同時打開2個柵格文件，執(zhí)行逐個格點的柵格值相減(若被減值為空值，則不執(zhí)行)，并將結(jié)果另存為新的柵格文件。

2.1.3 MODIS遙感數(shù)據(jù)處理

系統(tǒng)的綜合評價指數(shù)由5個評價指標計算得出，其中植被覆蓋指數(shù)需要從MODIS單通道數(shù)據(jù)(*.Vix)提取并計算得到。為了使用的規(guī)范化，專門編寫了此類遙感數(shù)據(jù)處理的類MODIS。Vix文件為二進制格式，包含了250個字節(jié)的文件頭，文件頭包含衛(wèi)星標識、日期時間、通道、投影方式、分辨率和經(jīng)緯度范圍等基本信息，以及火災(zāi)、水體和海水等識別閾值等數(shù)據(jù)信息。本系統(tǒng)用3個結(jié)構(gòu)體來存放文件頭。處理好文件頭后，用一個二維矩陣來存放指數(shù)數(shù)據(jù)。數(shù)據(jù)讀取按照文件頭獲取的行列數(shù)進行。

為了獲得更加可靠的數(shù)據(jù)，還需對得到的矩陣數(shù)據(jù)進行降噪處理。系統(tǒng)以直方圖統(tǒng)計的方法執(zhí)行降噪，編制了接口，并且提供動態(tài)的步長和閾值等參數(shù)設(shè)置。

2.2 GIS操作功能

系統(tǒng)的GIS操作實現(xiàn)了GIS應(yīng)用系統(tǒng)的常規(guī)功能，即地圖的漫游、縮放和移動等操作，并通過編程實現(xiàn)了圖層管理和屬性設(shè)置功能，如圖層的添加和刪除，矢量要素圖層的顏色、大小、形狀以及標注等的設(shè)定。

2.2.1 圖層管理

系統(tǒng)利用MapWinGIS提供的Map和LegendControl控件，設(shè)計了帶地圖的窗體，并以樹狀目錄的形式設(shè)計了圖層管理界面(圖1)。設(shè)計快捷工具條，提供圖層管理按鈕和地圖放大、縮小和漫游等基本操作按鈕，并提供圖層屬性查看、圖像輸出和工作向?qū)У瓤旖莅粹o。在菜單欄也可以添加各種預(yù)設(shè)的邊界地圖，以便使業(yè)務(wù)人員操作盡量便捷。LegendControl控件負責管理地圖中的圖層，實現(xiàn)圖層顯隱和屬性設(shè)置。

2.2.2 圖層屬性設(shè)置

MapWinGIS是開源的GIS引擎，提供了豐富的GIS功能接口，能實現(xiàn)大部分地圖和圖層功能，但其并沒有像ArcGIS那樣給出預(yù)開發(fā)的組件(或控件)，需要開發(fā)人員完成圖層管理、屬性和風格設(shè)置等功能的交互界面開發(fā)。系統(tǒng)針對上述功能，專門開發(fā)了圖層屬性設(shè)置窗體，將所有的圖層屬性管理集成到一個界面(圖4)。

圖4 地圖圖層屬性設(shè)置對話框

除常規(guī)的圖層名稱、要素形狀、大小、顏色和標注等，還完成了圖層的分色渲染功能，并實現(xiàn)了常用的多種字段值多種方式的分色，屬性設(shè)置比較完善。MapWinGIS的二次開發(fā)是開放性的，要實現(xiàn)的交互功能雖然比較麻煩，但是對開發(fā)人員來說非常自由。系統(tǒng)運用到的MapWinGIS提供的主要接口或類見圖3。

2.3 產(chǎn)品制作功能

按照業(yè)務(wù)流程，系統(tǒng)需要根據(jù)獲取的各地市濕潤指數(shù)、植被覆蓋指數(shù)、水體密度指數(shù)、土地退化指數(shù)和災(zāi)害指數(shù)等5種生態(tài)質(zhì)量指數(shù)的原始統(tǒng)計數(shù)據(jù)，分別計算各生態(tài)質(zhì)量指數(shù)，最后將各指數(shù)加權(quán)得出綜合評價指數(shù)。

2.3.1各生態(tài)質(zhì)量指數(shù)的計算

為了方便使用，我們將各生態(tài)質(zhì)量指數(shù)的計算界面設(shè)計為風格統(tǒng)一的窗體，根據(jù)計算的具體對象，對界面做相應(yīng)的細節(jié)調(diào)整。

各生態(tài)質(zhì)量指數(shù)的具體計算步驟為： 打開相應(yīng)的資料文件，設(shè)置好相應(yīng)的參數(shù)，給定結(jié)果保存的路徑，即可執(zhí)行計算。計算結(jié)果保存在文件中，同時顯示在列表控件中，并可以即時加入地圖中瀏覽，也可以單獨打開以前計算好的指數(shù)執(zhí)行瀏覽操作。每個生態(tài)質(zhì)量指數(shù)分別計算，其中植被覆蓋指數(shù)需要打開遙感數(shù)據(jù)并執(zhí)行降噪操作，其余指數(shù)的數(shù)據(jù)文件均為文本文件； 各生態(tài)質(zhì)量指數(shù)的數(shù)據(jù)包含輕、中、重3種程度的統(tǒng)計值，計算時需要按照設(shè)定的權(quán)重加權(quán)，權(quán)重在界面中可調(diào)。

生態(tài)質(zhì)量指數(shù)計算窗口如圖5所示。

圖5 生態(tài)質(zhì)量指數(shù)計算對話框

2.3.2 綜合評價指數(shù)

綜合評價指數(shù)的計算相對簡單。將之前計算得到的各個生態(tài)質(zhì)量指數(shù)文件打開，進行加權(quán)平均(各評價權(quán)重數(shù)值的設(shè)置不屬本研究內(nèi)容)，即可得到綜合評價指數(shù)。界面提供對各個權(quán)重的設(shè)置，也提供即時加入地圖中瀏覽和風格設(shè)置(圖6)。

圖6 綜合評價指數(shù)計算對話框

2.3.3產(chǎn)品制作

系統(tǒng)提供2種方式的產(chǎn)品制作。一種是在上述計算界面得到的指數(shù)數(shù)值列表中提供右鍵菜單的表格復制功能，表格復制后即可粘貼到服務(wù)產(chǎn)品的文檔中； 另一種為圖片產(chǎn)品，即在將指數(shù)加入地圖后，可以點擊“圖片輸出”按鈕，調(diào)用圖片輸出窗體，窗體即可將地圖中設(shè)置好的當前圖層制作為圖片，同時提供標題設(shè)置； 還可以將圖片復制到剪貼板，或者將圖片保存為文件，供服務(wù)產(chǎn)品的編輯使用。圖片輸出對話框如圖7所示。

圖7 圖片輸出對話框

3 結(jié)論

開源GIS在當前的地圖應(yīng)用開發(fā)中日益受到歡迎，其組件小巧方便，具有二次開發(fā)方便快捷、效率高和靈活性強等特點，特別適合輕量級業(yè)務(wù)應(yīng)用方面的開發(fā)。廣西生態(tài)質(zhì)量氣象評價業(yè)務(wù)系統(tǒng)的研發(fā)基本掌握了MapWinGIS的接口和有關(guān)類庫的使用方法，為進一步使用該開源組件打下很好的基礎(chǔ)。將開源GIS組件應(yīng)用到本業(yè)務(wù)平臺的研發(fā)中，是研究業(yè)務(wù)系統(tǒng)自主開發(fā)的一次有益嘗試。

本研究自主開發(fā)了多元回歸、反距離插值和遙感數(shù)據(jù)處理等C#類程序，實現(xiàn)了氣候數(shù)據(jù)的分析處理、MODIS數(shù)據(jù)讀取及處理和氣象評價指標分析等功能。利用MapWinGIS的接口和類庫實現(xiàn)了業(yè)務(wù)平臺中GIS功能的開發(fā)，包括地圖的縮放及漫游、圖層管理、數(shù)據(jù)導入、圖層屬性設(shè)置和評價結(jié)果顯示等，并利用GIS生成一系列業(yè)務(wù)產(chǎn)品。該平臺從業(yè)務(wù)流程需要出發(fā)，實現(xiàn)了氣象評價指標中濕潤指數(shù)、植被覆蓋指數(shù)、水體密度指數(shù)、土地退化指數(shù)及災(zāi)害指數(shù)5個指標數(shù)據(jù)的處理及計算，最后得到綜合評價指標。該評價系統(tǒng)為業(yè)務(wù)工作提供了標準化的工作流程。

參考文獻(References)：

[1] 王君華,莫偉華,陳燕麗,等.基于遙感的廣西生態(tài)質(zhì)量氣象評價方法研究[J].氣象科學,2009,29(6):833-837.

Wang J H,Mo W H,Chen Y L,et al.Research on meteorological evaluation of ecological quality in Guangxi based on satellite data[J].Scientia Meteorologica Sinica,2009,29(6):833-837.

[2] 肖建民.基于Web與GIS技術(shù)的黑龍江省防汛指揮信息服務(wù)子系統(tǒng)的設(shè)計與實現(xiàn)[D].大連:大連理工大學,2002.

Xiao J M.Design and Realization of Flood-Control Information Service Subsystem with Web and GIS in Heilongjiang Province[D].Dalian: Dalian University of Technology,2002.

[3] 戴長雷,遲寶明,林 嵐，等.基于GIS的地下水監(jiān)測管理信息系統(tǒng)(GMSMIS)分析與設(shè)計[J].遙感技術(shù)與應(yīng)用,2005,20(6):625-629.

Dai C L,Chi B M,Lin L,et al.Analysis and design of management information systems of groundwater monitoring based on GIS[J].Remote Sensing Technology and Application,2005,20(6):625-629.

[4] 王學登.廣東省水資源管理地理信息系統(tǒng)[J].中國農(nóng)村水利水電,2005(11):35-36.

Wang X D.Guangdong Province water resources management geographic information system[J].China Rural Water and Hydropower,2005(11):35-36.

[5] 閭國年,李 峰,邵進達.江蘇省公路地理信息系統(tǒng)建設(shè)規(guī)劃(上),(下)[J].交通科技,2000(1):18-22,2000(2):31-33.

Lyu G N,Li F,Shao J D.GIS construction plan of highway in Jiangsu Province,GIS construction plan of highway in Jiangsu Province (continued)[J].Transportation Science and Technology,2000(1):18-22,2000(2):31-33.

[6] 石紅兵,魏慶朝,曾學貴.地理信息系統(tǒng)在勘測設(shè)計一體化軟件開發(fā)過程中的應(yīng)用[J].鐵路航測,2000(2):20-22.

Shi H B,Wei Q C,Zeng X G.GIS application in the developing process of survey and design integration software[J].Railway Air Survey,2000(2):20-22.

[7] 程 玉,張 群,熊 英.基于GIS現(xiàn)代物流信息系統(tǒng)的設(shè)計[J].湖北工業(yè)大學學報,2006,21(4):95-97.

Cheng Y,Zhang Q,Xiong Y.A design of logistic information system based on GIS technology[J].Journal of Hubei University of Technology,2006,21(4):95-97.

[8] 黃治勇,袁正騰,萬 軍,等.基于GIS長江中游氣象水文預(yù)報系統(tǒng)的設(shè)計[J].暴雨災(zāi)害,2007,26(4):316-318.

Huang Z Y,Yuan Z T,Wan J,et al.Design of meteorological and hydrological forecasting system based on GIS in the middle reaches of Yangtze River[J].Torrential Rain and Disasters,2007,26(4):316-318.

[9] 李 良,黃生葉.基于GIS的山洪災(zāi)害預(yù)警系統(tǒng)的設(shè)計與實現(xiàn)[J].科學技術(shù)與工程,2006,6(17):2712-2715.

Li L,Huang S Y.Design and implementation of forewarning decision system of mountain torrent base on GIS[J].Science Technology and Engineering,2006,6(17):2712-2715.

[10] 沙 莎,曹 蕓,朱曉晨.基于GIS的氣象災(zāi)害評估系統(tǒng)——干旱評估模塊的研究[J].科技信息,2009(23):5-6,13.

Sha S,Cao Y,Zhu X C.Meteorological disaster evaluation system based on GIS:Drought assessment module research[J].Science and Technology Information,2009(23):5-6,13.

[11] 吳盛洪.基于GIS的氣象服務(wù)監(jiān)控管理平臺的設(shè)計與實現(xiàn)[D].成都:電子科技大學,2010.

Wu S H.The Design and Implementation of Meteorological Services Monitoring Management Platform based on GIS[D].Chengdu:University of Electronic Science and Technology of China,2010.

[12] 董學士,畢碩本,郭文政.基于GIS氣象查詢服務(wù)系統(tǒng)的設(shè)計和實現(xiàn)[J].微計算機信息,2010,26(7):146-147,137.

Dong X S,Bi S B,Guo W Z.The design and implementation of meteorology inquiry service system based on GIS[J].Control and Automation,2010,26(7):146-147,137.

[13] 李 軒,郭安紅,莊立偉.基于GIS的主要農(nóng)作物病蟲害氣象等級預(yù)報系統(tǒng)研究[J].國土資源遙感,2012,24(1):104-109.doi:10.6046/gtzyyg.2012.01.19.

Li X,Guo A H,Zhuang L W.The main crop diseases and pests meteorology level forecast system based on GIS[J].Remote Sensing for Land and Resources,2012,24(1):104-109.doi:10.6046/gtzyyg.2012.01.19.

[14] 熊文兵,孫 瑩,王艷蘭,等.基于GIS技術(shù)和天氣雷達的短臨預(yù)警平臺研究[J].遙感技術(shù)與應(yīng)用,2012,27(3):472-478.

Xiong W B,Sun Y,Wang Y L,et al.The study of early-warning platform based GIS technology and weather Radar[J].Remote Sensing Technology and Application,2012,27(3):472-478.

[15] 施益強,吳麗娜,吳陳鋒.基于GIS的雷達數(shù)據(jù)三維可視化與預(yù)警系統(tǒng)設(shè)計與實現(xiàn)[J].國土資源遙感,2013,25(1):171-175.doi:10.6046/gtzyyg.2013.01.30.

Shi Y Q,Wu L N,Wu C F.Design and implementation of a 3D visualization and early warning system for Radar data based on GIS[J].Remote Sensing for Land and Resources,2013,25(1):171-175.doi:10.6046/gtzyyg.2013.01.30.

[16] 高大偉,樊高峰,胡永亮,等.FY-3氣象衛(wèi)星微波水汽三維可視化及其在強天氣監(jiān)測中的應(yīng)用[J].國土資源遙感,2014,26(1):139-143.doi:10.6046/gtzyyg.2014.01.24.

Gao D W,Fan G F,Hu Y L,et al.3D visualization of humidity distribution from FY-3 MWHS and its application to severe weather monitoring[J].Remote Sensing for Land and Resources,2014,26(1):139-143.doi:10.6046/gtzyyg.2014.01.24.