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

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

2、e 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 monitor the haze， dust and othe

3、r 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 macroscale reflect the region

4、al distribution of the provinces air quality and pollution.【Key words】Arc Object； Visualization； Aerosol1 研究背景近年來，尤其是秋冬季節(jié)我省絕大多數(shù)城市的首要污染物是顆粒物，其作為氣溶膠物質(zhì)的一種長期以來一直是空氣環(huán)境質(zhì)量監(jiān)測的主要對象，同時(shí)也是空氣污染指數(shù)評價(jià)的重要因子。利用衛(wèi)星遙感技術(shù)通過動(dòng)態(tài)監(jiān)測氣溶膠光學(xué)厚度（Aerosol Optical Depth，簡稱AOD）等關(guān)鍵指標(biāo)因子，能夠有效解決地面固定觀測站空間覆蓋度不夠等問題1，并可進(jìn)一步輔助監(jiān)控灰霾、沙塵等空氣污染現(xiàn)象。氣溶膠光學(xué)

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

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

7、D等遙感技術(shù)監(jiān)測數(shù)據(jù)實(shí)時(shí)反映宏觀尺度全省各地區(qū)空氣質(zhì)量與污染分布狀況的目的。2.1 數(shù)據(jù)需求擬合模型所涉及的數(shù)據(jù)主要由三部分組成，其中輸入數(shù)據(jù)包括：衛(wèi)星遙感氣溶膠光學(xué)厚度數(shù)據(jù)，數(shù)據(jù)格式為柵格型，相對空氣濕度數(shù)據(jù)，可轉(zhuǎn)化為空間點(diǎ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ù)的形式進(jìn)行可視化表達(dá)，從而達(dá)到實(shí)時(shí)反映宏觀尺度全省各地區(qū)空氣質(zhì)量與污染分布的目的。3 設(shè)計(jì)思路3.1 模塊框架設(shè)計(jì)以江蘇省作為主要研究區(qū)域，以氣溶膠光學(xué)厚度AOD影像、

8、相對空氣濕度RH數(shù)據(jù)為基礎(chǔ)，參照遙感AOD數(shù)據(jù)處理流程，利用ArcObject組件庫實(shí)現(xiàn)實(shí)時(shí)高效地遙感數(shù)據(jù)解析，并引入研究構(gòu)建的AOD-API回歸統(tǒng)計(jì)模型，擬合計(jì)算得到全省各地區(qū)的API數(shù)值，同時(shí)利用ArcGIS軟件展現(xiàn)模擬結(jié)果，不僅大大提高了數(shù)據(jù)的處理效率，而且能夠直觀形象地表達(dá)全省范圍內(nèi)的空氣狀況模擬結(jié)果，并與自動(dòng)站實(shí)測數(shù)據(jù)進(jìn)行比對?？梢暬K的總體架構(gòu)分為界面層、業(yè)務(wù)邏輯層和數(shù)據(jù)訪問層三個(gè)層次。界面層用于采集用戶輸入的相關(guān)信息，展示地理信息和模擬結(jié)果。該層利用ArcGIS軟件界面，同時(shí)調(diào)用.Net定制的WinForm界面。業(yè)務(wù)邏輯層在模塊運(yùn)行提供地理信息處理功能，根據(jù)具體業(yè)務(wù)需求，業(yè)務(wù)層

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

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

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

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

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

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

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