数据:现在有国家气候中心发布的中国地区气候变化预估数据集,为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中批处理了。