基于遙感數(shù)據(jù)的區(qū)域大氣質(zhì)量空間分析模塊二次開發(fā)技術(shù)研究

1、    基于遙感數(shù)據(jù)的區(qū)域大氣質(zhì)量空間分析模塊二次開發(fā)技術(shù)研究    茅晶晶+姜晟【摘 要】顆粒物已經(jīng)成為我省主要城市的首要污染物，其作為氣溶膠物質(zhì)的一種長期以來一直是空氣環(huán)境質(zhì)量監(jiān)測的主要對象。利用衛(wèi)星遙感技術(shù)通過動態(tài)監(jiān)測氣溶膠光學(xué)厚度（aod）等關(guān)鍵指標(biāo)因子有利于進一步監(jiān)控灰霾、沙塵等空氣污染現(xiàn)象。本文利用gis二次開發(fā)手段將這種通過遙感技術(shù)與地面數(shù)據(jù)相結(jié)合的方式，將模型擬合后的環(huán)境空氣質(zhì)量分布與變化狀況以api指數(shù)的形式進行可視化表達，實時反映宏觀尺度全省各地區(qū)空氣質(zhì)量與污染分布的目的?！娟P(guān)鍵詞】ao；可視化；氣溶膠【abstract】partic

2、ulate matter has become the primary pollutant of the provinces major cities. as an aerosol substance， it has been a main target of ambient air quality monitoring. the technology of the satellite remote sensing can dynamic monitor the aerosol optical depth and other key index factors to further monit

3、or the haze， dust and other air pollution phenomena. complexed remote sensing technology and ground data， this paper gives the visualization of the process of air quality in the form of distribution and change of status api index， through the gis secondary development tools， reaching the goal of mac

4、roscale reflect the regional distribution of the provinces air quality and pollution.【key words】arc object； visualization； aerosol1 研究背景近年來，尤其是秋冬季節(jié)我省絕大多數(shù)城市的首要污染物是顆粒物，其作為氣溶膠物質(zhì)的一種長期以來一直是空氣環(huán)境質(zhì)量監(jiān)測的主要對象，同時也是空氣污染指數(shù)評價的重要因子。利用衛(wèi)星遙感技術(shù)通過動態(tài)監(jiān)測氣溶膠光學(xué)厚度（aerosol optical depth，簡稱aod）等關(guān)鍵指標(biāo)因子，能夠有效解決地面固定觀測站空間覆蓋度不夠等問題1，并

5、可進一步輔助監(jiān)控灰霾、沙塵等空氣污染現(xiàn)象。氣溶膠光學(xué)厚度（aod）被定義為介質(zhì)的消光系數(shù)在垂直方向上的積分，是描述氣溶膠對光的衰減作用的主要因子。它是目前可以得到的氣溶膠數(shù)據(jù)中覆蓋范圍最廣、較準(zhǔn)確的一種數(shù)據(jù)，也是推算氣溶膠含量、評估大氣污染程度、確定氣溶膠氣候效應(yīng)的關(guān)鍵因子2-3。gis可視化技術(shù)（gis visualization）是目前信息領(lǐng)域中廣泛應(yīng)用的一項技術(shù)，它通過強大、有效的地圖系統(tǒng)將復(fù)雜的空間數(shù)據(jù)和屬性數(shù)據(jù)以地圖學(xué)的形式進行描述，兼具人性化的界面風(fēng)格設(shè)計，實現(xiàn)了文本、圖形和圖像信息相結(jié)合的定位、查詢、檢索模式，信息表達形象直觀，使用操作簡單便捷4。通過遙感解譯數(shù)據(jù)與地面數(shù)據(jù)相結(jié)合

6、的方式，環(huán)境空氣質(zhì)量分布進行模型擬合利用arcgis平臺展示api指數(shù)擬合結(jié)果，在抽象的數(shù)值與具體的地理位置之間建立聯(lián)系，以特定的符號、顏色、圖表顯示在地圖上，實現(xiàn)衛(wèi)星遙感環(huán)境空氣質(zhì)量的可視化表達。同時基于地圖的可視化信息，查詢、統(tǒng)計功能也將變得更加簡單、直觀，通過可視化的查詢語言，實現(xiàn)對aod、api數(shù)據(jù)的圖表化查詢展現(xiàn)；通過建立專題地圖，幫助用戶加深理解，從而有助于發(fā)掘數(shù)據(jù)、信息間的內(nèi)在規(guī)律，實現(xiàn)輔助決策。2 需求分析通過提取美國對地觀測系統(tǒng)eos-modis衛(wèi)星遙感數(shù)據(jù)中的氣溶膠信息，結(jié)合地面數(shù)據(jù)及其他氣象資料，依據(jù)建立的線性回歸擬合模型演算出空氣污染指數(shù)api及對應(yīng)的環(huán)境空氣質(zhì)量等級，

7、并提供可視化表達界面，最終達到利用氣溶膠光學(xué)厚度aod等遙感技術(shù)監(jiān)測數(shù)據(jù)實時反映宏觀尺度全省各地區(qū)空氣質(zhì)量與污染分布狀況的目的。2.1 數(shù)據(jù)需求擬合模型所涉及的數(shù)據(jù)主要由三部分組成，其中輸入數(shù)據(jù)包括：衛(wèi)星遙感氣溶膠光學(xué)厚度數(shù)據(jù)，數(shù)據(jù)格式為柵格型，相對空氣濕度數(shù)據(jù)，可轉(zhuǎn)化為空間點圖層，數(shù)據(jù)格式為矢量型；輸出數(shù)據(jù)為空氣污染指數(shù)api，數(shù)據(jù)格式為柵格型。2.2 業(yè)務(wù)需求結(jié)合環(huán)境監(jiān)測具體業(yè)務(wù)需求，力求通過遙感技術(shù)與地面數(shù)據(jù)相結(jié)合的方式，將模型擬合后的環(huán)境空氣質(zhì)量分布與變化狀況以api指數(shù)的形式進行可視化表達，從而達到實時反映宏觀尺度全省各地區(qū)空氣質(zhì)量與污染分布的目的。3 設(shè)計思路3.1 模塊框架設(shè)計以

8、江蘇省作為主要研究區(qū)域，以氣溶膠光學(xué)厚度aod影像、相對空氣濕度rh數(shù)據(jù)為基礎(chǔ)，參照遙感aod數(shù)據(jù)處理流程，利用arcobject組件庫實現(xiàn)實時高效地遙感數(shù)據(jù)解析，并引入研究構(gòu)建的aod-api回歸統(tǒng)計模型，擬合計算得到全省各地區(qū)的api數(shù)值，同時利用arcgis軟件展現(xiàn)模擬結(jié)果，不僅大大提高了數(shù)據(jù)的處理效率，而且能夠直觀形象地表達全省范圍內(nèi)的空氣狀況模擬結(jié)果，并與自動站實測數(shù)據(jù)進行比對?？梢暬K的總體架構(gòu)分為界面層、業(yè)務(wù)邏輯層和數(shù)據(jù)訪問層三個層次。界面層用于采集用戶輸入的相關(guān)信息，展示地理信息和模擬結(jié)果。該層利用arcgis軟件界面，同時調(diào)用.net定制的winform界面。業(yè)務(wù)邏輯層在模

9、塊運行提供地理信息處理功能，根據(jù)具體業(yè)務(wù)需求，業(yè)務(wù)層把界面上用戶輸入的信息通過業(yè)務(wù)規(guī)則分解，需要進行地理信息處理的，把參數(shù)推送給ao組件的接口，并接受運行結(jié)果再傳遞給界面顯示；需要與屬性信息之間進行交互的，通過適配器讀取屬性信息進行運算，業(yè)務(wù)層起到承上啟下的作用。數(shù)據(jù)訪問層主要用于模塊運行時所要涉及的空間數(shù)據(jù)與屬性數(shù)據(jù)，提供接口供業(yè)務(wù)邏輯層訪問調(diào)用。 3.2 模塊開發(fā)環(huán)境氣溶膠數(shù)據(jù)反演模塊采用.net開發(fā)平臺，利用arcgis平臺提供的arcobject接口進行二次開發(fā)，實現(xiàn)aod影像數(shù)據(jù)解析、回歸模型擬合、結(jié)果統(tǒng)計渲染等功能。arcobject（簡稱ao）是esri公司arcgis家族中應(yīng)用

10、程序arcmap、arccatalog和arcscene的開發(fā)平臺，它是基于com（components object model 對象組件模型）技術(shù)所構(gòu)建的一系列com組件集，具有很強的gis功能和制圖功能，有1800多個組件、幾百個具有良好文檔說明的接口和數(shù)千個方法組成，作為arcgis可重用的通用二次開發(fā)組件集，以其強大的功能類庫，對地理信息系統(tǒng)（gis）的功能實現(xiàn)和擴展表現(xiàn)得十分出色5。3.3 模塊數(shù)據(jù)來源1）空間數(shù)據(jù)采用的是1：250，000江蘇省行政區(qū)劃圖、江蘇省省轄市行政區(qū)劃圖等。2）濕度數(shù)據(jù)模型構(gòu)建中涉及的濕度（rh）數(shù)據(jù)是由省氣象局提供，時間跨度由2008年7月1日至2009

11、年6月30日共計一年。在模塊計算過程中通過excel適配器讀取excel固定格式的濕度數(shù)據(jù)，根據(jù)時間、省轄市名稱獲得相應(yīng)值。4 開發(fā)過程根據(jù)課題的研究結(jié)果，設(shè)計參數(shù)的設(shè)置界面，在界面上選擇待處理的aod數(shù)據(jù)、研究區(qū)域、工作路徑等參數(shù)，根據(jù)參數(shù)首先對aod數(shù)據(jù)名稱進行解析，得出具體aod數(shù)據(jù)日期，同時進行格式轉(zhuǎn)換，剔除異常值并獲取研究區(qū)域內(nèi)的aod點圖層，通過idw插值補全全省aod數(shù)據(jù)。另一方面讀取excel統(tǒng)計表中各相關(guān)監(jiān)測點位的空氣相對濕度數(shù)據(jù)，根據(jù)樣本日期提取出當(dāng)日的相對濕度值，生成相對濕度點圖層，再通過idw插值方法生成全省濕度數(shù)據(jù)。之后，根據(jù)樣本日期所處的季節(jié)選擇相應(yīng)的aod-api

12、回歸統(tǒng)計模型，進行氣溶膠光學(xué)厚度aod與相對空氣濕度rh的公式擬合計算，得到全省范圍內(nèi)各地區(qū)的api估算值，并與實測數(shù)據(jù)進行比對，分別利用分層設(shè)色渲染，可視化展現(xiàn)模擬結(jié)果。5 關(guān)鍵技術(shù)應(yīng)用1）刪-矢數(shù)據(jù)轉(zhuǎn)換技術(shù)氣溶膠光學(xué)厚度（aod）數(shù)據(jù)是柵格格式，其中包含異常值，不利于進一步的空間分析，利用gis平臺提供的數(shù)據(jù)格式轉(zhuǎn)換conversion工具，將柵格數(shù)據(jù)轉(zhuǎn)換成矢量數(shù)據(jù)，在剔除異常值后，進行插值，獲得研究區(qū)域內(nèi)較為完整的數(shù)據(jù)分布。在多個圖層參與計算時，柵格數(shù)據(jù)較為簡便，此時將插值后的aod數(shù)據(jù)和濕度數(shù)據(jù)都轉(zhuǎn)換成柵格數(shù)據(jù)格式，利用raster calculator工具，對獨立兩個柵格圖層進行空間

13、數(shù)值計算，得到最終模型擬合結(jié)果。2）gis插值技術(shù)氣溶膠光學(xué)厚度（aod）數(shù)據(jù)解譯后，為了展示整個區(qū)域的空氣質(zhì)量狀況，需要將未知數(shù)據(jù)地區(qū)的污染物濃度值利用插值方法得出。arcgis提供了反距離權(quán)重插值法、樣條函數(shù)法、克里金法等一些特定用途的空間插值函數(shù)。經(jīng)過優(yōu)化模擬可知，反距離權(quán)重插值法的結(jié)果能更好地反映江蘇省地區(qū)的實際污染情況。反距離權(quán)重插值以插值點與實際觀測樣本點之間的距離為權(quán)重。反距離權(quán)重插值法要求對受影響的局部控制點有清楚的認(rèn)識，且其結(jié)論直接受到采樣點數(shù)值的影響，采樣點越多，對局部的真實反映越強，利用反距離權(quán)重插值法所確定的污染范圍就越準(zhǔn)確。3）分層設(shè)色渲染技術(shù)平臺的空間渲染模塊為管理

14、人員提供區(qū)域空氣質(zhì)量空間分析動態(tài)插值與渲染功能，直觀地展示空氣質(zhì)量實時空間分布特征，實現(xiàn)對監(jiān)測因子的空氣質(zhì)量空間分析動態(tài)插值渲染功能，直觀地展示了區(qū)域空氣質(zhì)量實時空間分布特征。平臺的空間渲染模塊需要根據(jù)各站點的小時值、日均值等實時生成污染物濃度空間渲染圖，利用gis技術(shù)實現(xiàn)這一目標(biāo)。6 可視化模塊展示1）加載預(yù)處理后的modis氣溶膠光學(xué)厚度aod柵格數(shù)據(jù)圖層。2）參數(shù)調(diào)整設(shè)置選擇氣溶膠模擬模塊，根據(jù)主要參數(shù)設(shè)置窗體的提示設(shè)置運算所需的參數(shù)。主要包括，待解析的遙感aod數(shù)據(jù)、參與運算的江蘇省行政區(qū)劃面狀圖層、根據(jù)選擇aod數(shù)據(jù)的名稱模塊會自動解析出aod數(shù)據(jù)所處的日期、確定運算過程中是否考慮濕

15、度因子、根據(jù)解析時間和是否考慮濕度因子自動顯示模塊運算的公式、參與模塊運算的外部數(shù)據(jù)（xls格式）的位置和模塊運行時產(chǎn)生中間和結(jié)果數(shù)據(jù)存儲的文件夾路徑。3）過程演算與專題制圖開始計算后，經(jīng)過后臺的數(shù)據(jù)轉(zhuǎn)換、idw插值、柵格運算、圖層渲染、結(jié)果輸出等步驟，最終生成一個系列圖層，主要包括aod參照圖層、api擬合圖層等。7 結(jié)論基于arcgis平臺二次開發(fā)了aod-api擬合與可視化表達模塊，實現(xiàn)遙感aod數(shù)據(jù)和rh數(shù)據(jù)自動推演出api數(shù)據(jù)的計算過程，并實現(xiàn)運行結(jié)果表明，rs/gis可視化技術(shù)支持下的全省空氣狀況及變化趨勢模擬與傳統(tǒng)表征手段相比，能夠更好的體現(xiàn)不同地區(qū)空氣環(huán)境質(zhì)量的區(qū)域分異性，更為準(zhǔn)確地顯示顆粒物污染的程度和變化情況?！緟⒖嘉墨I】1王家成，朱成杰，等.北京地區(qū)多氣溶膠遙感參量與pm2.5相關(guān)性研究j.中國環(huán)境科學(xué)，2015，35（7）：1947-1956.2單楠，楊曉暉，等.基于modis的中國陸地氣溶膠光學(xué)厚