arcgis中python计算面积的表达式_解析ArcGis的字段计算器（一）——数值型数据计算，从“面积计算”开始...

先来点儿背景知识铺垫：

ArcMap的字段计算器提供了两种脚本语言的支持用以计算，两种脚本语言是VBScript与Python。

多数人选择使用前者，因为它的基本函数和Excel的函数貌似一样。注意我这里用了一个“貌似”，虽然Excel函数与VB函数有着千丝万缕的关系，但它毕竟不是VB函数（ArcMap里用VBScript），把Excel函数照搬进ArcMap的计算器，许多是不可以运行的。使用VBScript函数正确的打开方式应该是查询VBScript帮助。至于VB、VBA、VBScript啥关系，不讲。

在ArcMap的字段计算器里，VBScript语句中字段（Field）应该用“[field_Name]”的样式输入，Python语句中字段应该用“!field_Name!”的样式输入。（提醒“新来的”，①不要把这玩意儿“”抄过去，②所有符号应为英文的）。

Tip：数值型与字符串型拼接的常见错误

不论在VBScript还是Python，"面积："+1.2 这样的写法都是错误的。

+适用于字符串类型的连接，在VB中，字符串与数值连接可以使用 & ；而在Python中， "面积："+str(1.2) 这样一个转换是必要的。

获取图斑（面要素）面积：

面积值是double型的，所以要放到Double型字段计算。

方法一：常用的，是在字段标题上右键——计算几何（Calculate Geometry）

方法二：使用Python语句的方式，在字段计算器中“解析语言”勾选Python，在语句块输入框填 !shape.area!

方法二中使用Python语句，在ArcMap中要素的一些几何（Geometry）属性，可以通过“!shape.？!”的方式进行计算，这个我们以后继续分解。

有“新来的”可能会遇到面积不能计算的情况，这一般是因为数据的空间参考没有设置好，至于“空间参考”是神马玩意，不知道的请百度或者发红包给作者，包教不会。

面积的单位转换（平方米转亩）： 直接[MJ]*0.0015

面积的舍入与取整：

1、舍入

一般我们会采用四舍五入的方式进行小数舍入，下面我们来贴图

VBScript，使用round( [MJ],3)对面积字段的值取小数后3位，结果如图：

Python（ArcGis10.1使用Python2.7），使用round( !MJ! ,3)对面积字段的值取小数后3位，结果就不贴图了，和上边一样。

有点意思吧？和你想的结果一样吗？一般使用中这样不会有问题，下面部分是闲扯扩展，有兴趣可读。

在编程语言中小数舍入一般采用“四舍六入五凑双”的方式，说的通俗点就是看舍入位数之后的数

距离哪个近就进位位给谁：

保留1位的话2.449是2.4 2.451是2.5

如果一样远（5），那就凑双数：

如果取3位，第4位是末位值为5，进位后第3位不是双数就不进，所以2.2325是2.232,而2.2315也是2.232。

那么，怎么实现四舍五入呢？

①VBScript用 FormatNumber( [MJ],3)试试；

②Python用下面的函数试试：

def getRound(val,n):

"""

:function:替换内置round函数,实现保留n位小数的精确四舍五入

:param val: 类型string或数值型

:param n: 需要保留小数的位数，要求大于0

:return: 四舍五入后的数据，类型string

"""

n = int(n)

f = '{:.%sf}'%n

return f.format(round(float(val) * pow(10,n)) / pow(10.0,n))

不知道怎么使用？给你个见下图

③确定没问题的方法——Excel函数处理后挂接回来。

2、 取整

取整一般有两种，一种是直接取整数位，一种是舍入取整。

第一种在VBScript或者Python里都可以使用int函数实现

第二种使用round函数，这时候就有比较大的差距了。

题外话，再来点随机凑数

vb的随机数函数是rnd，它返回一个小于1，大于等于0的数

如果生成2-10之间的数，可以用 rnd()*2+8 不过有一个问题是它永远取不到10。