使用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中批处理了。