1. 程式人生 > >使用ArcGIS批處理txt格點資料為柵格資料並插值

使用ArcGIS批處理txt格點資料為柵格資料並插值

   資料:現在有國家氣候中心釋出的中國地區氣候變化預估資料集,為RCP4.5、RCP8.5情景下未來氣候變化預估資料,例如RCP45情景下2006年1月份的降雨量資料為

這樣的格式,有全國的格點資料,格點解析度為0.5*0.5度,時間解析度為月。

目的:用於模型,則還需要將其按1KM*1KM插值成柵格資料。

方法:(以降雨資料為例)

一、ToRaster

1.使用ArcToolBox,開啟conversion tools-toraster-ascii to raster,右鍵點選Batch(批處理)。


2.先新增5個做個實驗,然後全選。如圖所示:


3.使用批處理,我的方法是用Excel完成,在桌面新建一個Excel表格,裡面按arcGis上面所示格式輸入,如圖:


4.然後貼上到ArcGis,點選OK即可。這樣就完成資料轉換為柵格的了。結果命名為pr_year+month.

二、投影

1.在ArcToolbox裡面找到Data Management tools-Projections and transformations-define projection工具,使用它Batch處理。(具體操作同上就不貼圖了)

2.使用Data Management tools-Projections and transformations-raster-project raster工具,batch批處理用來投影轉換。投影之後可能並不會看到變化的結果,不要擔心,這是arcgis的Data frame問題,你可以重啟arcgis然後開啟投影之後的檔案就會看到變化了,不過我覺得這樣沒必要,因為只是顯示問題並不影響後面的處理。

3.由上面柵格投影生成的柵格的命名方式為:prp_year+month.

三、裁剪

1.使用data managementtool-raster-raster processing-clip工具,也是用Excel貼上到batch的框框中批處理。

2.通過clip工具剪下的目標區域柵格的命名方式為:prd_year+month.效果如下:


四、生成point點檔案,及插值

1.使用conversion Tools-Fromraster-raster to point工具,生成點要素檔案,由以上投影柵格生成點要素檔案的命名方式為prs_year+month.shp。

2.使用spaticalAnalysis-interpolation-idw工具來插值,你也可以用其他方法插值(像克里金法,樣條曲線法,最小距離法等等),本實驗要求使用IDW而已。生成的IDW插值柵格的命名方式為:prw_year+month。插值效果如圖:


3.矯正,經過以上步驟發現最後的結果prw200601,檢視其Layer properties的columns and rows 與我所需要的行列數不一致,因此需要矯正。此處用到的方法是將標準檔案r2001_1加到arcgis。

4.然後使用Spatial analyst工具,先設定其option的analyst extent為r2001_1,然後點選raster  calculator,輸入r2006_1=[prw_200601] + [r_2000_1] * 0,得到的r2006_1即為最終結果。最後一步也可以在Excel上做好,然後複製到新建txt檔案中,替換Excel自動生成的空格即可把一堆公式加到raster calculator中批處理了。