[GIS]计算面积

缘起

做了好一阵子项目，每每要算面积总是怪自己记忆甚差，又忘了怎么做，作为一个GIS门外汉把解决过程记录一次，希望帮助和我一样初入门就被卡住的小夥伴。

座标

一开始总用偷懒的方式不想转换投影座标，直接在field用计算器算（记得选python）

!shape.geodesicArea@SQUAREKILOMETERS!

也可以直接用py执行

arcpy.CalculateField_management(shp, "AreaKM2",  "!shape.geodesicArea@SQUAREKILOMETERS!")

但有几次用的我新作的shape file是得不到面积的，因为没有定义投影

定义投影

我一直没花时间弄懂投影这块知识点，每次都不知道要选哪种座标系，所以都选择跟其他图层一样。
首先，define projection，在转换 project找一个座标系即可。