缘起
做了好一阵子项目,每每要算面积总是怪自己记忆甚差,又忘了怎么做,作为一个GIS门外汉把解决过程记录一次,希望帮助和我一样初入门就被卡住的小夥伴。
座标
一开始总用偷懒的方式不想转换投影座标,直接在field用计算器算(记得选python)
!shape.geodesicArea@SQUAREKILOMETERS!
也可以直接用py执行
arcpy.CalculateField_management(shp, "AreaKM2", "!shape.geodesicArea@SQUAREKILOMETERS!")
但有几次用的我新作的shape file是得不到面积的,因为没有定义投影
定义投影
我一直没花时间弄懂投影这块知识点,每次都不知道要选哪种座标系,所以都选择跟其他图层一样。
首先,define projection,在转换 project找一个座标系即可。