第一次用字段计算器算面积的时候我在地图上看到结果完全懵了——一个看起来接近两万平方米的地块算出来只有0.0004。后来才反应过来这块数据本身是经纬度坐标系字段计算器返回的是平方度不是平方米。这个坑很典型而且和今天要讲的几何属性提取密不可分。在ArcGIS里处理点、线、面数据很多时候最终需求都落在几个几何量上面积、周长、长度、坐标、质心、范围。这些属性其实一直存在要素里只是默认不写到属性表中。而字段计算器就是把它们批量落到表里的最直接工具学会之后不用再一遍遍打开计算几何对话框也不用手动量算一行表达式就能跑完整张表。这篇文章会把几何属性提取这件事讲透彻为什么要用字段计算器而不是右键计算几何、shape对象里到底藏着哪些属性、坐标系和字段类型会在什么地方坑你以及几个我一直觉得很有用的进阶写法。不管你是刚接触ArcGIS的新手还是已经用了几年但习惯停留在鼠标操作的老用户这份内容应该都能让人有点收获。1. 为什么要把几何属性写进属性表三个真实场景1.1 场景一属性表里没有面积字段统计和筛选寸步难行做宗地图、地块统计、林地小班汇总这类工作最常见的一个尴尬是图形画得清清楚楚但属性表里就是没有一个现成的面积字段。没有面积字段就意味着你没法在属性表里直接排序看哪个地块最大没法用SQL语句选出面积大于某个阈值的所有要素也没法在符号化的时候按面积区间直接做分级渲染。很多人遇到这种情况第一反应是右键打开「计算几何」把面积算出来临时在表里看一眼。这当然可以但计算几何的结果如果不刻意保留本质上只是一个临时查看动作而且它只适合单要素逐条地看。如果你想按面积做筛选比如把面积超过十公顷的图斑全部挑出来你必须在属性表里有一个真实存在的数值字段SQL查询才能生效。字段计算器恰恰是干这个的先把面积写进字段后面筛选、排序、制图、汇总全部水到渠成。1.2 场景二坐标、质心是空间分析的基础数据另一种高频需求是把坐标和质心变成普通属性。比如外业采集回来的点要素你想在Excel报表里输出每个点的经纬度或者做国土调查时要给每个图斑生成一个质心点再和附近设施做空间连接。这些操作的本质都是把几何属性比如点位坐标、面要素质心坐标提取出来从空间域转化成属性域。坐标和质心这类数据你在地图上可以看到但拿不到表格里。字段计算器里写!shape.X!、!shape.centroid.X!就可以把几十万条要素的坐标和质心一次性写入字段分分钟变成可统计、可导出的普通数值。后续做转换、做分析、做报表都比在图形界面里一个个查要高效得多。1.3 计算几何对话框与字段计算器到底该用哪个能力点计算几何字段计算器结果写入属性表可以可以支持只计算选中要素支持支持且会弹出确认提示支持自定义逻辑固定选项基本无扩展可写Python/VBScript代码块自由度极高批量处理能力一个图层逐次操作可配合脚本、模型构建器批量跑坐标系选择可选要素类或数据框坐标系按要素类自身坐标系不可动态切换结果持久化写入所选字段写入所选字段我的习惯是只是偶尔看一眼面积用计算几何没问题。但如果这个字段之后要参与筛选、制图、报表或下游分析就老老实实建字段用字段计算器。一次建好后续所有流程都不用来回折腾。2. 字段计算器背后的几何对象shape到底能暴露什么2.1 先看解析程序Python和VBScript的语法差异字段计算器里新建完字段、右键打开以后最先要注意的就是解析程序Parser选项。ArcGIS 10.x里面常见的两个解析程序是VBScript和PythonArcGIS Pro里主要用Python 3。两个语法在引用字段和几何对象上的写法差异很大选错了表达式直接报错。Python引用字段是!字段名!引用几何对象是!shape!访问属性用点号。VBScript引用字段是[字段名]引用几何对象是[shape]。这里最容易犯的错就是把Python的语法写进VBScript或者反过来。比如想在Python里求面要素的面积应该是!shape.area!如果解析程序是VBScript就得写[shape].area。一旦搞混计算器会告诉你表达式无效。我自己的建议是从今天起统一用Python。VBScript在ArcGIS Pro里已经被边缘化Python生态里可以做条件判断、字符串处理、甚至写多个函数VBScript能做的事情它全部能做而且写出来的东西更接近现代脚本语言的感觉以后过渡到ArcPy也不会陌生。2.2 面、线、点要素的几何属性一览几何对象在字段计算器里并不仅仅是一个图形框它能暴露出很多有用的属性。下面这个表是我实际工作中最常用的几种要素类型想获取的属性Python表达式VBScript表达式面面积!shape.area![shape].area面周长!shape.length![shape].length线长度!shape.length![shape].length点X坐标!shape.X![shape].X点Y坐标!shape.Y![shape].Y面/点质心X!shape.centroid.X![shape].centroid.x面/点质心Y!shape.centroid.Y![shape].centroid.y面/线范围最小X!shape.extent.XMin![shape].extent.XMin面/线范围最小Y!shape.extent.YMin![shape].extent.YMin这里面有几个初学者容易绕晕的地方。首先面要素也有length属性它指的就是多边形的周长别因为名字叫length就觉得只适用于线。其次线要素没有area所以别试图给线图层算面积先想清楚自己拿到的数据类型是什么。另外点要素没有length、area能用的核心就是X、Y坐标。2.3 一个容易忽略的细节!shape!的嵌套属性访问字段计算器里一个很好用的特性是几何属性可以层层往下访问。!shape!是一个完整的地理对象你可以沿着它去取子对象的子属性比如!shape.centroid.X!就是先取质心这个点再取这个点的X坐标。同样!shape.extent.XMin!是先取要素的包围盒范围再取范围的最小X值。这种嵌套逻辑如果你能理解很多需求就不必借助别的工具。比如要判断一条道路是南北走向还是东西走向可以比较!shape.extent.XMin!和!shape.extent.XMax!的差值大小也可以比较起点和终点的坐标差。我甚至见过有人用嵌套属性计算线要素的走向方位角虽然表达式写起来长一些但效果相当稳定。嵌套访问的本质是面向对象逻辑每个几何对象里有子对象每个子对象还有属性。把这个思路理清字段计算器对你来说就从“一个算面积的对话框”变成“一个花园口感的表达式引擎”。3. 实操新建双精度字段用一行表达式批量算完3.1 新建字段时类型选错后面全白干字段计算器能不能正确拿到几何属性的结果很大程度取决于你新建字段时的数据类型。我去很多项目组培训的时候发现一个问题反复出现大家觉得面积是个数值就随手选了个长整型结果几百个图斑的面积全被截成整数还有人选了文本型算出来的面积变成科学计数法的字符串后续一排序全乱。正确的做法是选双精度Double。面积和长度这类几何量经常有大数和小数双精度可以存15位有效数字浮点型只有7位左右双精度在绝大多数情况下的精度都够了。如果是存坐标我也建议双精度不然小数点后几位容易被丢掉。总之和几何计算挂钩的新建字段一律双精度起步不要给自己埋坑。3.2 字段计算器界面与表达式录入的完整步骤直接走一遍完整流程右键图层打开属性表。点击表视图右上角的菜单按钮选择“添加字段”名称填Area_Ha类型选“双精度”。在字段名Area_Ha上右键选择“字段计算器”。解析程序选择 Python勾选“显示代码块”可以后续放复杂函数。表达式框里输入!shape.area! / 10000这样就把平方米转成了公顷。点击确定全表瞬间填满。如果当前有选中要素打开字段计算器的时候ArcGIS会弹出一个提示问你是想让这个计算作用于全部记录还是只作用于选中记录。很多人没仔细看就点了确定结果发现只算了一部分然后丈二和尚摸不着头脑。这个提示一定看清楚你的本意是什么就选什么。表达式录入的时候建议先在记事本或代码编辑器里写好再粘贴进去。字段计算器的输入框有一些历史遗留的怪癖比如某些版本的输入框对中文引号特别敏感你在Word里打了一对中文双引号粘贴进去十有八九会报错。先在纯文本环境里写好能省掉一堆莫名其妙的麻烦。3.3 两个可以直接抄的作业面积周长、经纬度坐标给面要素同时计算面积和周长最简单直接的方式就是建两个字段然后分别设置表达式。比如面积字段Area_SQM表达式!shape.area!周长字段Perim_M表达式!shape.length!如果想要自动换算成公顷就在表达式里写!shape.area! / 10000。给点要素提取经纬度坐标建字段Lon和Lat表达式分别写!shape.X!和!shape.Y!。需要注意的是如果数据本身是投影坐标系算出来的坐标是投影坐标比如米不是经纬度。想输出经纬度先把要素类投影到WGS84或CGCS2000地理坐标系再执行字段计算器。这两个作业覆盖了日常最频繁的一类需求。只要把流程走通一次后面无非是换字段名、换表达式整套操作逻辑是完全一致的。4. 算出来的面积对不上号坐标系、单位与动态投影的坑4.1 为什么有人算出来是0.000324而不是3245平方米回到开头那个例子。一个地块地图上看过去挺大字段计算器算出来却是0.000324这个结果本身没错错的是坐标系。当要素类所处的是地理坐标系如CGCS2000、WGS84它的单位是度字段计算器返回的shape.area单位就是平方度。而平方度和平方米没有一个固定的换算系数因为1度代表的实际地面长度在不同纬度上是不一样的所以千万不要拿平方度结果去乘一个常数往平方米上凑。正确的做法是先把数据投影到适合当地范围的投影坐标系比如全国范围可以选择Albers等积投影中高纬度地区可以用UTM或高斯-克吕格投影。投影之后字段计算器算出来的shape.area才是平方米或该投影的单位平方。这一步很多教程都会提但强调得不够投影这个动作一定要在字段计算之前完成而不是在地图上选了一个显示坐标系就完事。还有一种情况是要素类压根没定义坐标系。这时用字段计算器算几何属性轻则返回空值重则直接报错。先打开图层属性查看源选项卡确认坐标系那一栏不是未知。如果是未知先用“定义投影”工具Define Projection给它指定正确的坐标系然后再投影到目标坐标系。4.2 动态投影不影响字段计算器这是最常见误解ArcMap和ArcGIS Pro的地图窗口中有一个东西叫“数据框坐标系”也就是你用视图 → 数据框属性里设置的坐标系。它只决定这个地图窗口用什么投影来显示数据并不修改数据本身。很多人在数据框里把坐标系切换成Web墨卡托之后发现地图上多边形的形状和面积好像变了于是想当然地认为字段计算器也算的是墨卡托面积。实际上字段计算器读的是要素类自身的坐标系统。你对地图窗口做的动态投影只是换了一个“观看视角”底层数值一个字节都没动。所以如果你在地图上看到的面积极大或者极小但字段计算器的结果又是另一个数别迷信眼睛先去看数据的真实坐标系。这是我在项目里见过最普遍的现象数据本身是WGS84地图框改成Web墨卡托用户拿地图上估算的面积去比对字段计算器结果然后怀疑自己表达式写错了。如果你确实需要用地图框的坐标系来计算面积字段计算器做不到但“计算几何”对话框里可以勾选“使用数据框的坐标系”这种情况下建议优先用计算几何。如果非要在字段计算器里实现就只能先把要素类投影到目标坐标系再执行计算。4.3 字段类型对结果的精度影响坐标系决定了算出来的数值是什么单位字段类型则决定这个数值在表里最终存成什么样。前面提过双精度这里再展开一点。ArcGIS的数值型字段分短整型、长整型、浮点型、双精度几种字段类型存储范围适用场景短整型-32768 到 32767小整数、分类编码长整型约 -21亿 到 21亿整数ID、人口数浮点型约7位有效数字只需要粗略数值时双精度约15位有效数字面积、周长、坐标等精确几何量面积字段如果选了浮点型图斑大一点就可能在个位甚至十位就开始偏差这在后续汇总统计时会造成看起来不大但实际不可忽略的误差。我在做地籍属性整理时通常所有几何计算字段全部用双精度宁可存储空间多占一点也不给精度留隐患。单精度浮点型的设计初衷是节省空间但在今天的硬件条件下双精度的存储消耗完全可以忽略这个便宜不占也罢。5. 进阶玩法质心坐标、Python逻辑与批量处理5.1 把面要素质心坐标写进属性表给地块加质心坐标是我做专题图时最常干的活儿。多边形要素的质心在字段计算器里对应!shape.centroid.X!和!shape.centroid.Y!当图层有几十万个要素时这个表达式仍然能一批跑完速度非常可观。质心对不规则多边形尤其有意义。比如做农区调查时一个图斑可能是狭长的河滩地或凹多边形质心不等于形心它是按面积加权计算出来的几何中心比简单地取顶点平均更符合“这个图斑的中心在哪里”的直觉。写进字段之后可以在符号化时直接按质心放置标注也可以生成点图层做空间连接。生成点图层那一步用“要素转点”工具也能做但如果你只是想拿到坐标数值而不是真的生成一份点数据字段计算器更轻量也不产生额外文件。5.2 用Python代码块做条件运算面积分级自动打标签字段计算器并不仅限于写一行表达式勾选“显示代码块”之后可以在输入框里写完整的Python函数再在下方表达式里调用。比如要根据面积给图斑打“大中小”标签def classify(area): if area 10000: return 大 elif area 1000: return 中 else: return 小表达式框里写classify(!shape.area!)这样每个要素都会根据面积字段动态判断自动写入大、中、小。整个过程不需要先算面积再做属性分类两个动作一次完成。同理你可以在代码块里定义多个函数或者用round函数对数值做四舍五入。有网友问过“字段计算器怎么只保留后三位数字”这里说一下如果是对计算出来的数值保留三位小数直接在表达式里写round(!shape.area!, 3)即可如果你是从一个文本字段里提取最后三位字符串可以用Python的切片语法!编码字段![-3:]。两者应对的是不同需求先分清你的字段是数值型还是文本型再选对应写法。5.3 大文件与多图层场景下的批量字段计算一个图层手动操作还算轻松如果工程里有几十个图层都需要批量补面积和周长逐一手动点就很浪费时间。这时可以用ArcPy脚本来做核心函数是CalculateField_managementimport arcpy fc rC:\data\地块.shp field_name Area_Ha expression !shape.area! / 10000 code_block arcpy.CalculateField_management(fc, field_name, expression, PYTHON, code_block)放在循环里就可以一次性对多个要素类执行同样的字段计算。这个脚本也可以挂到模型构建器里做成一个通用工具输入任意要素类、指定字段名和表达式剩下的交给模型跑。批量计算场景下我更推荐在脚本开头检查要素类的坐标系一旦发现是地理坐标系脚本直接跳出提示防止全部算出色调离谱的平方度结果。宁可多写两行判断也好过处理完几十个文件之后才发现结果全部作废。5.4 一个值得记住的扩展测地线面积最后聊一个容易被大多数人忽略的细节。前面说投影后再算面积这在绝大多数场景下都对但如果数据跨越的经纬度范围特别大比如全球尺度的格网或多省大范围空间分析任何投影平面都存在一定变形。ArcGIS Pro的字段计算器里可以试试!shape.geodesicArea!它直接按椭球面计算测地线面积不依赖局部投影的精度通常用于跨大区域或对精度要求较高的分析。ArcMap老版本里这个属性不一定能用我更常用的折中方案是先把数据投影到Albers等积投影再字段计算。等积投影专门保证了面积关系不变结果相对可靠。说到底每一个几何值都有它的适用边界清楚自己在什么坐标系、什么分析尺度下工作比死记表达式更重要。我在实际项目里养成的一个习惯是任何要素类入库的第一步先建好面积和周长字段并计算完毕再谈后续分析。这样后续所有筛选、制图、共享数据的人都能直接在属性表里拿到关键几何量不需要每个人重复踩一遍坐标系的坑。希望这篇围绕字段计算器提取几何属性的分享能让你的ArcGIS操作路径短一点、结果稳一点。
