【arcgis】使用ArcGIS完成空间插值——趋势面分析法

一、目标

制作一副年平均降雨量的地图。

二、数据

某地的175个气象站数据的shp文件station.shp，以及这个地方轮廓的栅格数据idoutlgd。

链接：https://pan.baidu.com/s/1SES7nu_wg8uGBhidTVLSTA 密码：khik

三、制作方法

1.首先加载数据。

2.在菜单栏/customize/toolbars/中找到geostatistical analyst工具 

然后，点击下拉菜单，点击geostatistical wizard，进入界面

然后在右侧的source dataset选择station文件，datefield选择ann_prec，同时右侧选择global polynomial interpolation，效果如下图：

然后点击next，此时选择趋势面模型的幂。通过改变幂的值，找到趋势面均方根的最小值。

通过调节可得，当幂为5时，效果较好。

最后点击finsh，生成Global Polynomial Interpolation，这是一张地统计分析生成的地图，地图范围与station一致。

然后右击这个图层，选择属性（properites）。symbology的按钮下的show栏有四个选项。选中filled contours，清除其余的，然后点击分类（classify），将method改为manual，将class改为7。也就是说我们的降雨分级是10,15,20,25,30,35，然后点击ok，完成标记。.

然后要对这张生成的图进行切割，使其分为与我们之前给出的栅格数据的范围要一致才行。

首先要把这张生成的图进行栅格化。

选中图层，右键data—export to raster。在出现的对话框中，输入2000m作为像元的大小，然后命名即可。

然后就可以开始进行切割。

双击spatital analyst tools/Extraction工具下的extract by mask工具。

选择刚才生成的栅格图像作为输入，idotlgd作为掩膜数据，指定输出名称。

完成切割。切割效果如下。

然后用户可以修改一下配色方案。

直接点击图册修改颜色即可。

我大概改了一下，效果如下：

然后我们还可以生成等值线。

点击spatital analyst下拉菜单中的surface analyst，选择contour。

出现的对话框中选择好输入图层，在contour interval中输入5，base contour输入10，设置好输入图层。

右键生成的图层，选择属性，在labels按钮下，点击label features。从label字段中选择contour，至此，等值线添加完成。