ArcGIS中的空间插值和面积计算

说明:本文阐述了空间插值和污染面积估算的方法，供群内交流学习用，若要用于商业用途或转载，请与原作者联系。本文若有不正确之处，敬请指出！

一、空间插值

插值方法种类很多，每种插值方法里参数也很多，至于哪种最好，没有定论，只能根据需求以及制图的效果来选定。建议：插值效果图与网格图进行对比，哪种效果最接近网格图（能体现局部）而且又能反映整体趋势就取哪种。

1.1、

1.2、以“反距离权重法，1次方”为例：

请问：此处有可选smooth ，可以做进行平滑处理吗？

可以，但精度会受到影响，看平滑后的效果来决定是否进行平滑处理。建议不做

3、扩展研究区域

4、至此可以制作分层设色图filled contours/等值线图contours

为减少误差，还可以对分级进行设置

请问：此处分级该如何设置？有无相应依据？

含量图主要根据百分含量，如果作图效果不好，适当调整

评价图根据污染等级

5、这是采用“反距离权重法，1次方”来插值的。

可选用“局部多项式”或“普通克里格插值”方法来试试，看哪种和网格分级图更接近些。但无论哪种方法聚类误差可能都较大，一部分高值可能被掩盖。

二、下面转成栅格图层再进行分层设色图制作，这样精度较高，且图层可用来进行面积估算

2.1、导出成栅格图层

2.2、设置格网大小，一般在50到100左右(本次都设为100)

（2.3和2.4均非必要步骤，只是为了另外的处理或制图的美观性。如果是为了制图的美观性有可能这两个步骤会弄巧成拙，是否须要请根据具体需要和效果来定）

2.3、并可对栅格图层重分类，生成新的栅格图层如（ah_cd）

请问：此处重分类又该如何设置？有无相应依据？同上

注：此处生成的文件名请以p开头的各元素名，如p_AS，后面面积估算时会用上一定要放在E:\SOIL\BACKUP目录下

2.4 转成矢量图层

2.4 做平滑处理（可选）

2.5 至此可进行分层设色，再成图

2.6 再掩膜（或切割），最后成图

掩膜可以提高作图效率

切割

三、面积计算

准备两部分数据：

遥感解译的土地利用数据经过大类合并后转成的grid 数据，文件名useland；以及插值生成的各元素污染程度的格网数据，文件名以p开头,如p_AS。

两部分数据坐标系必须是等积投影而且必须一致

上面两个文件都放在E:\SOIL\BACKUP目录下,

另外，把tjgdldwrmj.aml和SmrzPollueteArea.exe两个文件也考到此目录下

然后程序—附件—命令提示符，在DOS窗口运行

>ARC

再输入ARC:&W O:\SOIL\BACKUP

再输入ARC:&run tjgdldwrmj.aml

再运行SmrzPollueteArea.exe（需要.NET2.0支持，如无，请安装）

最后得到污染面积汇总表

将用地类型数据先转成GeoDatabase格式，然后采用下面命令整合成六大类，并进行投影转换

update ahld set tmpcl=left(LD2006_ID,1)

——————————————

（由于开始是直接在DOS下运行，不能拷贝，我自己根据做的在你的文档上修改的，可能大小写等有点不一样，只是表达了意思，这些也不需给作图人员，请你看看我的做法可对）

3.1 Cover坐标转换（如果转换不了，先不做跳过这个步骤，导出成栅格图层后在toolbox

下进行转换）

Arc: project cover ld2006 ld2006_new

**************************************************

* The INPUT projection has been defined. *

**************************************************

Use OUTPUT to define the output projection and END

to finish.

Project: output

Project: PROJECTION ALBERS

Project: Spheroid Krasovsky

Project: UNITS METERS

Project: PARAMETERS

1st standard parallel [ 0 0 0.000 ]: 30.5 0 0 2nd standard parallel [ 0 0 0.000 ]: 33.5 0 0 central meridian [ 0 0 0.000 ]: 117 0 0 latitude of projection's origin [ 0 0 0.000 ]: 0 0 0

false easting (meters) [ 0.00000 ]: 500000 false northing (meters) [ 0.00000 ]: 0 Project: end

Arc: build ld2006_new poly

3.1在workstation下增加一个字段tmpcl,用来存储六大类信息

Arc: w d\soil

Arc: ae

//Arcedit: &wo e:\ah

Arcedit: edit ld2006_new

The edit coverage is now E:\AH\HLD2006_new

Arcedit: ef poly

75061 element(s) for edit feature POLY

Arcedit: additem tmpcl 8 8 i

Adding TMPCL to PAT...

Arcedit: items

COLUMN ITEM NAME WIDTH OUTPUT TYPE N.DEC ALTERNATE NAME INDEXED?

1 AREA 8 18 F 5 -

9 PERIMETER 8 18 F 5 -

17 LD2006_new# 4 5 B - -

21 LD2006_NEW-ID 4 5 B - -

25 TMPCL 8 8 I - - Arcedit: save

Saving changes for E:\AH\LD2006_new

Saving polygons...

For PAL of E:\AH\LD2006_new:

75061 polygon records(s) written.

from the original 75061, 75061 added and 75061 deleted.

Saving arcs...

141234 arc(s) written to E:\AH\LD2006_new

from the original 141234, 0 added and 0 deleted

Reopening arcs...

Saving labels...

75061 polygon attribute record(s) written to E:\AH\LD2006_new

75060 label(s) written to E:\AH\LD2006_new

from the original 75060, 0 added and 0 deleted

Reopening labels...

Reopening polygons...

BND replaced into E:\AH\LD2006_new

Saving set tolerances to TOL file...

Re-establishing edit feature POLYGON

3.2整合成六大类

Arc: ae

Arcedit: edit ld2006_new

The edit coverage is now E:\AH\HLD2006

Arcedit: ef poly

75061 element(s) for edit feature POLY