arcgis python添加几何属性
import arcpy
import numpy
import math def AddGeometryAttributes(fc, geomProperties, lUnit, aUnit, cs):
"""-------------------------------------------------------------------------
Tool: Add Geometry Attributes (Data Management Tools)
Source Name: AddGeometryAttributes.py
Version: ArcGIS 10.2.1
Author: Esri, Inc.
Usage: arcpy.AddGeometryAttributes_management(
Input_Features,
Geometry_Properties,
{Length_Unit},
{Area_Unit},
{Coordinate_System})
Required Arguments: Input Features
Geometry Properties
Optional Arguments: Length Unit
Area Unit
Coordinate System
Description: Adds attribute fields to the input features containing
measurements and coordinate properties of the feature
geometries (for example, length or area).
Updated: Not yet.
------------------------------------------------------------------------""" desc = arcpy.Describe(fc)
hasZ, hasM = desc.hasZ, desc.hasM
if cs:
sr = arcpy.SpatialReference()
sr.loadFromString(cs)
try:
srMetersPerUnit = sr.metersPerUnit
except:
srMetersPerUnit = 1
else:
try:
srMetersPerUnit = desc.spatialReference.metersPerUnit
except:
srMetersPerUnit = 1
shapeDict = {"POINT":1,
"MULTIPOINT":1.5,
"POLYLINE":2,
"POLYGON":3,}
shapeDim = shapeDict[str(desc.shapeType).upper()]
del desc fields = CreateOutputFields(fc, geomProperties, hasZ, hasM) arcpy.SetProgressor("STEP", arcpy.GetIDMessage(86174), 0,
int(arcpy.GetCount_management(fc).getOutput(0)), 1) hasNulls = False
# Calculate geometry properties into new fields
with arcpy.da.UpdateCursor(fc,fields + ["SHAPE@"],"",cs) as ucur:
ucurFields = ucur.fields
global ucurFields
for row in ucur:
global row
geom = row[ucurFields.index("SHAPE@")]
if geom:
if shapeDim == 1:
if "POINT_X_Y_Z_M" in geomProperties:
Update("POINT_X", geom.firstPoint.X)
Update("POINT_Y", geom.firstPoint.Y)
if hasZ:
Update("POINT_Z", geom.firstPoint.Z)
if hasM:
Update("POINT_M", geom.firstPoint.M)
if shapeDim>1:
if "PART_COUNT" in geomProperties:
Update("PART_COUNT", geom.partCount)
if "CENTROID" in geomProperties:
Update("CENTROID_X", geom.trueCentroid.X)
Update("CENTROID_Y", geom.trueCentroid.Y)
if hasZ:
Update("CENTROID_Z", geom.trueCentroid.Z)
if hasM:
Update("CENTROID_M", geom.trueCentroid.M)
if "EXTENT" in geomProperties:
Update("EXT_MIN_X", geom.extent.XMin)
Update("EXT_MIN_Y", geom.extent.YMin)
Update("EXT_MAX_X", geom.extent.XMax)
Update("EXT_MAX_Y", geom.extent.YMax)
if shapeDim>=2:
if "POINT_COUNT" in geomProperties:
Update("PNT_COUNT", geom.pointCount)
if "LINE_START_MID_END" in geomProperties:
Update("START_X", geom.firstPoint.X)
Update("START_Y", geom.firstPoint.Y)
if shapeDim ==2:
midPoint = geom.positionAlongLine(0.5,
True).firstPoint
else:
line = arcpy.Polyline(geom.getPart(0), "#",
hasZ, hasM)
if line.length > 0:
midPoint = line.positionAlongLine(0.5,
True).firstPoint
else:
hasNulls = True
del line
Update("MID_X", midPoint.X)
Update("MID_Y", midPoint.Y)
Update("END_X", geom.lastPoint.X)
Update("END_Y", geom.lastPoint.Y)
if hasZ:
Update("START_Z", geom.firstPoint.Z)
Update("MID_Z", midPoint.Z)
Update("END_Z", geom.lastPoint.Z)
if hasM:
Update("START_M", geom.firstPoint.M)
Update("MID_M", midPoint.M)
Update("END_M", geom.lastPoint.M)
del midPoint
if "CENTROID_INSIDE" in geomProperties:
Update("INSIDE_X", geom.centroid.X)
Update("INSIDE_Y", geom.centroid.Y)
if hasZ:
Update("INSIDE_Z", geom.centroid.Z)
if hasM:
Update("INSIDE_M", geom.centroid.M)
if shapeDim==2:
if "LINE_BEARING" in geomProperties:
lat1 = geom.firstPoint.Y
lon1 = geom.firstPoint.X
lat2 = geom.lastPoint.Y
lon2 = geom.lastPoint.X
Update("BEARING", (90 - math.degrees(math.atan2(lat2-lat1, lon2-lon1))) % 360)
del lat1, lon1, lat2, lon2
if "LENGTH" in geomProperties:
Update("LENGTH", ConvertFromMeters("LINEAR",
geom.length,
lUnit,
srMetersPerUnit))
if "LENGTH_3D" in geomProperties:
Update("LENGTH_3D", ConvertFromMeters("LINEAR",
geom.length3D,
lUnit,
srMetersPerUnit))
if "LENGTH_GEODESIC" in geomProperties:
Update("LENGTH_GEO", ConvertFromMeters("LINEAR",
geom.getLength("PRESERVE_SHAPE"),
lUnit,
srMetersPerUnit))
if shapeDim==3:
if "PERIMETER_LENGTH" in geomProperties:
Update("PERIMETER", ConvertFromMeters("LINEAR",
geom.length,
lUnit,
srMetersPerUnit))
if "AREA" in geomProperties:
Update("POLY_AREA", ConvertFromMeters("AREA",
geom.area,
aUnit,
srMetersPerUnit))
if "AREA_GEODESIC" in geomProperties:
Update("AREA_GEO", ConvertFromMeters("AREA",
geom.getArea("PRESERVE_SHAPE"),
aUnit,
srMetersPerUnit))
if "PERIMETER_LENGTH_GEODESIC" in geomProperties:
Update("PERIM_GEO", ConvertFromMeters("LINEAR",
geom.getLength("PRESERVE_SHAPE"),
lUnit,
srMetersPerUnit))
ucur.updateRow(row)
else:
hasNulls = True
arcpy.SetProgressorPosition()
if hasNulls:
arcpy.AddIDMessage("WARNING", 957) def CreateOutputFields(fc, geomProperties, hasZ, hasM):
propDict = {"POINT_X_Y_Z_M": ["POINT_X",
"POINT_Y",
"POINT_Z",
"POINT_M"],
"PART_COUNT": ["PART_COUNT"],
"CENTROID": ["CENTROID_X",
"CENTROID_Y",
"CENTROID_Z",
"CENTROID_M"],
"EXTENT": ["EXT_MIN_X",
"EXT_MIN_Y",
"EXT_MAX_X",
"EXT_MAX_Y"],
"POINT_COUNT": ["PNT_COUNT"],
"LINE_START_MID_END": ["START_X",
"START_Y",
"START_Z",
"START_M",
"MID_X",
"MID_Y",
"MID_Z",
"MID_M",
"END_X",
"END_Y",
"END_Z",
"END_M"],
"LINE_BEARING": ["BEARING"],
"CENTROID_INSIDE": ["INSIDE_X",
"INSIDE_Y",
"INSIDE_Z",
"INSIDE_M"],
"LENGTH": ["LENGTH"],
"PERIMETER_LENGTH": ["PERIMETER"],
"AREA": ["POLY_AREA"],
"LENGTH_GEODESIC": ["LENGTH_GEO"],
"AREA_GEODESIC": ["AREA_GEO"],
"LENGTH_3D": ["LENGTH_3D"],
"PERIMETER_LENGTH_GEODESIC":["PERIM_GEO"],
}
if not hasZ:
propDict["POINT_X_Y_Z_M"].remove("POINT_Z")
propDict["CENTROID"].remove("CENTROID_Z")
propDict["CENTROID_INSIDE"].remove("INSIDE_Z")
propDict["LINE_START_MID_END"].remove("START_Z")
propDict["LINE_START_MID_END"].remove("MID_Z")
propDict["LINE_START_MID_END"].remove("END_Z")
if not hasM:
propDict["POINT_X_Y_Z_M"].remove("POINT_M")
propDict["CENTROID"].remove("CENTROID_M")
propDict["CENTROID_INSIDE"].remove("INSIDE_M")
propDict["LINE_START_MID_END"].remove("START_M")
propDict["LINE_START_MID_END"].remove("MID_M")
propDict["LINE_START_MID_END"].remove("END_M") addList = []
geomPropertiesList = []
currentFields = [field.name for field in arcpy.ListFields(fc)]
for prop in geomProperties:
for field in propDict[prop.upper()]:
geomPropertiesList.append(field)
if field not in currentFields:
addList.append(field)
else:
arcpy.AddIDMessage("WARNING", 1097, field) if addList:
narray = numpy.array([], numpy.dtype([("_ID", numpy.int)] +
[(n, numpy.float) for n in addList]))
arcpy.da.ExtendTable(fc, "OID@", narray, "_ID") return geomPropertiesList def ConvertFromMeters(type, value, unit, metersPerUnit):
if not unit:
return value
else:
distanceUnitInfo = {"METERS": 1.0,
"FEET_US": 0.304800609601219,
"NAUTICAL_MILES": 1852.0,
"MILES_US": 1609.34721869444,
"KILOMETERS": 1000.0,
"YARDS": 0.914401828803658,} areaUnitInfo = {"ACRES": 4046.86,
"HECTARES": 10000.0,
"SQUARE_METERS": 1.0,
"SQUARE_FEET_US": 0.09290341161327473,
"SQUARE_NAUTICAL_MILES": 3429904.0,
"SQUARE_MILES_US": 2589998.4703195295,
"SQUARE_KILOMETERS": 1000000.0,
"SQUARE_YARDS": 0.8361307045194741,}
if type == "LINEAR":
return (value*metersPerUnit)/distanceUnitInfo[unit]
if type == "AREA":
return (value*math.pow(metersPerUnit,2))/areaUnitInfo[unit] def Update(field, value):
if value:
if not math.isnan(value):
row[ucurFields.index(field)] = value #run the script
if __name__ == '__main__':
# Get Parameters
fc = arcpy.GetParameterAsText(0)
if arcpy.GetParameterAsText(1).find(";") > -1:
geomProperties = arcpy.GetParameterAsText(1).upper().split(";")
else:
geomProperties = [arcpy.GetParameterAsText(1).upper()]
lUnit = arcpy.GetParameterAsText(2)
aUnit = arcpy.GetParameterAsText(3)
cs = arcpy.GetParameterAsText(4)
if not cs:
cs = arcpy.env.outputCoordinateSystem # Run the main script
AddGeometryAttributes(fc, geomProperties, lUnit, aUnit, cs)
arcgis python添加几何属性的更多相关文章
- python 添加类属性
类属性必须赋值. 创建类属性 类是模板,而实例则是根据类创建的对象. 绑定在一个实例上的属性不会影响其他实例,但是,类本身也是一个对象,如果在类上绑定一个属性,则所有实例都可以访问类的属性,并且,所有 ...
- ArcGIS + Python 批量裁剪、添加X/Y坐标脚本
前言 前一段时间,同事拿来的数据范围太大,用不了那么多(只需要一个乡镇的,结果拿来区县的),太多了加载也是问题.所以就让我给处理下. 由于文件较多,手动裁剪的话,我一个一个用ArcGIS工具箱中的工具 ...
- arcgis python arcpy add data script添加数据脚本
arcgis python arcpy add data script添加数据脚本mxd = arcpy.mapping.MapDocument("CURRENT")... df ...
- python装饰器、继承、元类、mixin,四种給类动态添加类属性和方法的方式(一)
介绍装饰器.继承.元类.mixin,四种給类动态添加类属性和方法的方式 有时候需要給类添加额外的东西,有些东西很频繁,每个类都需要,如果不想反复的复制粘贴到每个类,可以动态添加. # coding=u ...
- Python小白学习之如何添加类属性和类方法,修改类私有属性
如何添加类属性和类方法,修改类私有属性 2018-10-26 11:42:24 类属性.定义类方法.类实例化.属性初始化.self参数.类的私有变量的个人学习笔记 直接上实例: class play ...
- Arcgis Engine 添加一个Symbol符号样式步骤
public static void DrawPictureMarkerSymbol(IGlobe globe, String layerName) { //添加一个图层 ESRI.ArcGIS.Ca ...
- ArcGIS Python人门到精通目录基于ArcGIS10.2,100以上案例15章42个视频806分钟,51GIS网站上线
ArcGIS Python人门到精通目录 闫老师 QQ:276529800 微信13108507190 1. ArcGIS Python基础 1.1 ArcGIS为什么学习Python 1.2 A ...
- python限定类属性的类属性:__slots__
__slots__ 由于Python是动态语言,任何实例在运行期都可以动态地添加属性. 如果要限制添加的属性,例如,Student类只允许添加 name.gender和score 这3个属性,就可以利 ...
- Python中的属性管理
Python管 理属性的方法一般有三种:操作符重载(即,__getattr__.__setattr__.__delattr__和 __getattribute__,有点类似于C++中的重载操作符).p ...
随机推荐
- LeetCode828. Unique Letter String
https://leetcode.com/problems/unique-letter-string/description/ A character is unique in string S if ...
- Effective STL 学习笔记 Item 21:Comparison Function 相关
Effective STL 学习笔记 Item 21:Comparison Function 相关 */--> div.org-src-container { font-size: 85%; f ...
- mongo oplog 整理
首先需要介绍一下mongodb local库的作用 local库是MongoDB的系统库,记录着时间戳和索引和复制集等信息 test:PRIMARY> use local switched to ...
- Decorator 装饰
意图 动态地给一个对象添加一些额外的职责.就增加功能来说,Decorator模式相比生成子类更为灵活. 结构 Component:定义一个对象接口,可以给这些对象动态地添加职责:(纯虚函数) Conc ...
- 使用VSCode配置简单的vue项目
由于最近要使用的项目框架为前后端分离的,采用的是vue.js+webAPI的形式进行开发的.因为之前我没有接触过vue.js,也只是通过视频文档做了一些简单的练习.今天技术主管说让大家熟悉下VSCod ...
- 老的API实现WordCount
使用Hadoop版本0.x实现单词统计 package old; import java.io.IOException; import java.net.URI; import java.util.I ...
- 四、django rest_framework源码之频率控制剖析
1 绪言 权限判定之后的下一个环节是访问频率控制,本篇我们分析访问频率控制部分源码. 2 源码分析 访问频率控制在dispatch方法中的initial方法调用check_throttles方法开始. ...
- ActiveMQ (三):项目实践
1. 简单项目demo Com.hoo.mq路径下(除了com.hoo.mq.spring)是普通java中使用activeMQ. Com.hoo.mq.spring路径下是非web环境spring集 ...
- iOS 11开发教程(九)iOS11数据线连接真机测试
iOS 11开发教程(九)iOS11数据线连接真机测试 在Xcode 7.0之后,苹果公司在开发许可权限上做了很多的改变,在测试App方面取消了一些限制.在Xcode7.0之前的版本,苹果公司只向注册 ...
- luoguP4555 [国家集训队]最长双回文串 manacher算法
不算很难的一道题吧.... 很容易想到枚举断点,之后需要处理出以$i$为开头的最长回文串的长度和以$i$为结尾的最长回文串的长度 分别记为$L[i]$和$R[i]$ 由于求$R[i]$相当于把$L[i ...