基于遙感數據的區域大氣質量空間分析模塊二次開發技術研究

茅晶晶+姜晟

【摘 要】顆粒物已經成為我省主要城市的首要污染物，其作為氣溶膠物質的一種長期以來一直是空氣環境質量監測的主要對象。利用衛星遙感技術通過動態監測氣溶膠光學厚度（AOD）等關鍵指標因子有利于進一步監控灰霾、沙塵等空氣污染現象。本文利用GIS二次開發手段將這種通過遙感技術與地面數據相結合的方式，將模型擬合后的環境空氣質量分布與變化狀況以API指數的形式進行可視化表達，實時反映宏觀尺度全省各地區空氣質量與污染分布的目的。

【關鍵詞】AO；可視化；氣溶膠

【Abstract】Particulate 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 monitor 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 macroscale reflect the regional distribution of the provinces air quality and pollution.

【Key words】Arc Object； Visualization； Aerosol

1 研究背景

近年來，尤其是秋冬季節我省絕大多數城市的首要污染物是顆粒物，其作為氣溶膠物質的一種長期以來一直是空氣環境質量監測的主要對象，同時也是空氣污染指數評價的重要因子。利用衛星遙感技術通過動態監測氣溶膠光學厚度（Aerosol Optical Depth，簡稱AOD）等關鍵指標因子，能夠有效解決地面固定觀測站空間覆蓋度不夠等問題[1]，并可進一步輔助監控灰霾、沙塵等空氣污染現象。

氣溶膠光學厚度（AOD）被定義為介質的消光系數在垂直方向上的積分，是描述氣溶膠對光的衰減作用的主要因子。它是目前可以得到的氣溶膠數據中覆蓋范圍最廣、較準確的一種數據，也是推算氣溶膠含量、評估大氣污染程度、確定氣溶膠氣候效應的關鍵因子[2-3]。

GIS可視化技術（GIS Visualization）是目前信息領域中廣泛應用的一項技術，它通過強大、有效的地圖系統將復雜的空間數據和屬性數據以地圖學的形式進行描述，兼具人性化的界面風格設計，實現了文本、圖形和圖像信息相結合的定位、查詢、檢索模式，信息表達形象直觀，使用操作簡單便捷[4]。

通過遙感解譯數據與地面數據相結合的方式，環境空氣質量分布進行模型擬合利用ArcGIS平臺展示API指數擬合結果，在抽象的數值與具體的地理位置之間建立聯系，以特定的符號、顏色、圖表顯示在地圖上，實現衛星遙感環境空氣質量的可視化表達。同時基于地圖的可視化信息，查詢、統計功能也將變得更加簡單、直觀，通過可視化的查詢語言，實現對AOD、API數據的圖表化查詢展現；通過建立專題地圖，幫助用戶加深理解，從而有助于發掘數據、信息間的內在規律，實現輔助決策。

2 需求分析

通過提取美國對地觀測系統EOS-MODIS衛星遙感數據中的氣溶膠信息，結合地面數據及其他氣象資料，依據建立的線性回歸擬合模型演算出空氣污染指數API及對應的環境空氣質量等級，并提供可視化表達界面，最終達到利用氣溶膠光學厚度AOD等遙感技術監測數據實時反映宏觀尺度全省各地區空氣質量與污染分布狀況的目的。

2.1 數據需求

擬合模型所涉及的數據主要由三部分組成，其中輸入數據包括：①衛星遙感氣溶膠光學厚度數據，數據格式為柵格型，②相對空氣濕度數據，可轉化為空間點圖層，數據格式為矢量型；輸出數據為空氣污染指數API，數據格式為柵格型。

2.2 業務需求

結合環境監測具體業務需求，力求通過遙感技術與地面數據相結合的方式，將模型擬合后的環境空氣質量分布與變化狀況以API指數的形式進行可視化表達，從而達到實時反映宏觀尺度全省各地區空氣質量與污染分布的目的。

3 設計思路

3.1 模塊框架設計

以江蘇省作為主要研究區域，以氣溶膠光學厚度AOD影像、相對空氣濕度RH數據為基礎，參照遙感AOD數據處理流程，利用ArcObject組件庫實現實時高效地遙感數據解析，并引入研究構建的AOD-API回歸統計模型，擬合計算得到全省各地區的API數值，同時利用ArcGIS軟件展現模擬結果，不僅大大提高了數據的處理效率，而且能夠直觀形象地表達全省范圍內的空氣狀況模擬結果，并與自動站實測數據進行比對。

可視化模塊的總體架構分為界面層、業務邏輯層和數據訪問層三個層次。界面層用于采集用戶輸入的相關信息，展示地理信息和模擬結果。該層利用ArcGIS軟件界面，同時調用.Net定制的WinForm界面。業務邏輯層在模塊運行提供地理信息處理功能，根據具體業務需求，業務層把界面上用戶輸入的信息通過業務規則分解，需要進行地理信息處理的，把參數推送給AO組件的接口，并接受運行結果再傳遞給界面顯示；需要與屬性信息之間進行交互的，通過適配器讀取屬性信息進行運算，業務層起到承上啟下的作用。數據訪問層主要用于模塊運行時所要涉及的空間數據與屬性數據，提供接口供業務邏輯層訪問調用。

3.2 模塊開發環境

氣溶膠數據反演模塊采用.Net開發平臺，利用ArcGIS平臺提供的ArcObject接口進行二次開發，實現AOD影像數據解析、回歸模型擬合、結果統計渲染等功能。

ArcObject（簡稱AO）是ESRI公司ArcGIS家族中應用程序ArcMap、ArcCatalog和ArcScene的開發平臺，它是基于COM（Components object Model 對象組件模型）技術所構建的一系列COM組件集，具有很強的GIS功能和制圖功能，有1800多個組件、幾百個具有良好文檔說明的接口和數千個方法組成，作為ArcGIS可重用的通用二次開發組件集，以其強大的功能類庫，對地理信息系統（GIS）的功能實現和擴展表現得十分出色[5]。

3.3 模塊數據來源

1）空間數據

采用的是1：250，000江蘇省行政區劃圖、江蘇省省轄市行政區劃圖等。

2）濕度數據

模型構建中涉及的濕度（RH）數據是由省氣象局提供，時間跨度由2008年7月1日至2009年6月30日共計一年。在模塊計算過程中通過Excel適配器讀取Excel固定格式的濕度數據，根據時間、省轄市名稱獲得相應值。

4 開發過程

根據課題的研究結果，設計參數的設置界面，在界面上選擇待處理的AOD數據、研究區域、工作路徑等參數，根據參數首先對AOD數據名稱進行解析，得出具體AOD數據日期，同時進行格式轉換，剔除異常值并獲取研究區域內的AOD點圖層，通過IDW插值補全全省AOD數據。另一方面讀取Excel統計表中各相關監測點位的空氣相對濕度數據，根據樣本日期提取出當日的相對濕度值，生成相對濕度點圖層，再通過IDW插值方法生成全省濕度數據。之后，根據樣本日期所處的季節選擇相應的AOD-API回歸統計模型，進行氣溶膠光學厚度AOD與相對空氣濕度RH的公式擬合計算，得到全省范圍內各地區的API估算值，并與實測數據進行比對，分別利用分層設色渲染，可視化展現模擬結果。

5 關鍵技術應用

1）刪-矢數據轉換技術

氣溶膠光學厚度（AOD）數據是柵格格式，其中包含異常值，不利于進一步的空間分析，利用GIS平臺提供的數據格式轉換Conversion工具，將柵格數據轉換成矢量數據，在剔除異常值后，進行插值，獲得研究區域內較為完整的數據分布。

在多個圖層參與計算時，柵格數據較為簡便，此時將插值后的AOD數據和濕度數據都轉換成柵格數據格式，利用Raster Calculator工具，對獨立兩個柵格圖層進行空間數值計算，得到最終模型擬合結果。

2）GIS插值技術

氣溶膠光學厚度（AOD）數據解譯后，為了展示整個區域的空氣質量狀況，需要將未知數據地區的污染物濃度值利用插值方法得出。ArcGIS提供了反距離權重插值法、樣條函數法、克里金法等一些特定用途的空間插值函數。經過優化模擬可知，反距離權重插值法的結果能更好地反映江蘇省地區的實際污染情況。反距離權重插值以插值點與實際觀測樣本點之間的距離為權重。

反距離權重插值法要求對受影響的局部控制點有清楚的認識，且其結論直接受到采樣點數值的影響，采樣點越多，對局部的真實反映越強，利用反距離權重插值法所確定的污染范圍就越準確。

3）分層設色渲染技術

平臺的空間渲染模塊為管理人員提供區域空氣質量空間分析動態插值與渲染功能，直觀地展示空氣質量實時空間分布特征，實現對監測因子的空氣質量空間分析動態插值渲染功能，直觀地展示了區域空氣質量實時空間分布特征。平臺的空間渲染模塊需要根據各站點的小時值、日均值等實時生成污染物濃度空間渲染圖，利用GIS技術實現這一目標。

6 可視化模塊展示

1）加載預處理后的MODIS氣溶膠光學厚度AOD柵格數據圖層。

2）參數調整設置

選擇氣溶膠模擬模塊，根據主要參數設置窗體的提示設置運算所需的參數。主要包括，待解析的遙感AOD數據、參與運算的江蘇省行政區劃面狀圖層、根據選擇AOD數據的名稱模塊會自動解析出AOD數據所處的日期、確定運算過程中是否考慮濕度因子、根據解析時間和是否考慮濕度因子自動顯示模塊運算的公式、參與模塊運算的外部數據（xls格式）的位置和模塊運行時產生中間和結果數據存儲的文件夾路徑。

3）過程演算與專題制圖

開始計算后，經過后臺的數據轉換、IDW插值、柵格運算、圖層渲染、結果輸出等步驟，最終生成一個系列圖層，主要包括AOD參照圖層、API擬合圖層等。

7 結論

基于ArcGIS平臺二次開發了AOD-API擬合與可視化表達模塊，實現遙感AOD數據和RH數據自動推演出API數據的計算過程，并實現運行結果表明，RS/GIS可視化技術支持下的全省空氣狀況及變化趨勢模擬與傳統表征手段相比，能夠更好的體現不同地區空氣環境質量的區域分異性，更為準確地顯示顆粒物污染的程度和變化情況。

【參考文獻】

[1]王家成，朱成杰，等.北京地區多氣溶膠遙感參量與PM2.5相關性研究[J].中國環境科學，2015，35（7）：1947-1956.

[2]單楠，楊曉暉，等.基于MODIS的中國陸地氣溶膠光學厚度時空分布特征[J].中國水土保持科學，2010，10（5）：24-30.

[3]李成才，毛節泰，等.利用MODIS研究中國東部地區氣溶膠光學厚度的分布和季節變化[J].科學通報，2003，48（19）：2094-2100.

[4]孫玉鋒，郝天曙，章彧，陸惠勇.基于GIS的統計數據可視化技術研究與實踐[J].電子技術與軟件工程，2015，16：109-111.

[5]賴劍菲，江舟.基于ArcObjects的組件式GIS的開發與研究[J].國土資源信息化，2005：29-32.

[責任編輯：湯靜]