arcpy
YuanYWRS
展开
-
Arcpy 操作数据框中启用时间的图层
【介绍】:【DataFrameTime】对象可以执行时间管理,用来管理数据框中启用时间的图层。它是引用DataFrame.time属性返回的结果,可以检索当前时间【currentTime】、结束时间【endTime】、开始时间【startTime】、时间步长间隔【timeStepInterval】,数据框中的图层必须启用了时间属性才可以实现以上功能。【操作步骤】:在【1】处的【layer time】里选择【each feature has a single time field】;在【2】处的【time原创 2022-07-07 21:43:37 · 332 阅读 · 0 评论 -
Arcpy 利用updatelayer函数改变图层的符号系统
【介绍】:arcpy.mapping 里的updatelayer函数具有更新图层符号系统的功能,也可以用来更新各种图层属性,它作为一个多功能的函数,需要掌握其各个参数的使用情况。【处理前】:【代码示例】:【代码分析】:[第一行]:导入arcpy.mapping模块;[第二行]:引用当前的地图文档,并赋值给变量;[第三行]:获取对数据框的引用;[第四行]:定义将要更新的图层;如果更新图层的符号系统,必须保证更新图层和源图层具有相同的几何图形(点、线、面),还要检查属性定义是否相同,分级颜色和分原创 2022-07-03 22:24:35 · 831 阅读 · 0 评论 -
Arcpy——InsertLayer()函数的使用:掺入图层到地图文档里
【声明】:本书代码和数据来源于【基于arcgis的python编程秘笈】。【思路】:insertlayer()函数可以很准确的将数据添加插入到图层的某个位置中,它需要使用一个参考图层来指定位置,新图层将会添加在指定参考图层的上方或者下方。【注意】:InsertLayer()函数需要使用参考图层,所有不能对空的数据框使用该函数。【操作】:【结果展示】:...原创 2022-06-24 23:21:07 · 541 阅读 · 0 评论 -
Arcpy 添加图层到地图文档
【原理】:使用【AddLayer()函数】把图层添加到地图文档里。【说明】:本实验的代码思路以及数据均来自于牟乃夏老师的书【基于Arcgis的python编程秘笈】。【思路】:1:导入【arcpy.mapping模块】,并获取对当前活动地图文档的引用。2:创建一个新的变量【df】存储对crime数据框的引用,该引用是通过【ListDataFrames()函数】返回数据框列表使用【0】访问列表中的第1项元素来获取。3:创建一个实例【layer】,layer变量引用存储在本地的为school_Distric原创 2022-06-20 23:07:56 · 1114 阅读 · 0 评论 -
ArcGIS基础:基于python代码简单计算坡度、坡向和洪水淹没
1、坡度分析:读取栅格影像:rast_file=r'F:\arcpy_data\imagery\钦州湾.tif'... rast_obj=arcpy.Raster(rast_file)... rast_obj计算坡度:slp=arcpy.sa.Slope(rast_obj)输出结果:2、计算坡向:代码:aspect=arcpy.sa.Aspect(rast_obj)结果输出:3、洪水淹没分析:代码:flood=rast_obj<0.5结果输出:4、保存原创 2022-05-24 17:01:41 · 1521 阅读 · 0 评论 -
Arcpy-Pycharm 创建批量裁剪自定义工具
本实验讲解了三个裁剪方式:手动裁剪、代码裁剪和创建批量裁剪工具进行裁剪。下面是手动裁剪:这是手动裁剪的常见方式,功能位于分析工具下。上图是手动裁剪结果图。下面是代码裁剪:import arcpy>>> arcpy.env.overwriteOutput=True>>> out_features=r'F:\arcpy_data\钦州市\railway_qinnanqu.shp'>>> arcpy.Clip_analysis("铁路","原创 2022-05-20 17:27:38 · 423 阅读 · 0 评论 -
arcpy 实现新增字段合并字段
新增字段area_ha(公顷),并计算出来:arcpy.AddField_management('区界',"area_ha","DOUBLE")>>> arcpy.CalculateField_management("区界","area_ha","!AREA!/10000","PYTHON_9.3")运算界面:输出结果:新增字段area_MU(亩),并计算出来:arcpy.AddField_management('区界',"area_MU","DOUBLE")... a原创 2022-04-29 18:00:40 · 791 阅读 · 0 评论 -
arcpy 使用ListFields函数获取字段列表的名称、类型和长度
>>> import arcpy>>> arcpy.env.workspace=r'F:\arcpy_data\钦州市'>>> flds=arcpy.ListFields(r'兴趣点.shp')#获取要素类的字段列表>>> fld=flds[0]#获取第一个字段>>> print fld.name,fld.type,fld.length#打印字段名称、类型和长度FID OID 4#输出结果>>&原创 2022-04-28 17:45:22 · 2744 阅读 · 0 评论 -
arcpy 实现列举目录下的要素类与描述矢量数据要素类
>>> import arcpy#导入arcpy包... arcpy.env.workspace=ur'F:\arcpy_data\钦州市'#设置好工作空间... fcs=arcpy.ListFeatureClasses()#使用列举要素类的函数ListFeatureClass()... for fc in fcs:... print fc,#使用for循环打印出来。... print (fc)... BeibuGulfUniversity.shp POI_tmp.s原创 2022-04-28 13:06:38 · 398 阅读 · 0 评论 -
arcpy 基于python制作六十进制与十进制互换的工具箱
#将六十进制的度分秒转换为十进制的度:import mathdms=[1,2,3]deg=dms[0]+dms[1]/60.0+dms[2]/3600.0print (deg)#上述算法写为函数形式:def dmsdeg(dms): deg=dms[0]+dms[1]/60.0+dms[2]/3600.0 return degdmsdeg(dms)print (deg)开发为脚本:import arcpyimport mathdef dms2deg(dms):原创 2022-04-27 17:22:52 · 1764 阅读 · 0 评论 -
arcpy 创建弧度制与十进制互换的个人工具箱
import arcpy... rad=math.pi/2... deg=math.degrees(rad)... rad=math.radians(deg)... print (deg,rad)... (90.0, 1.5707963267948966)可以在arcgis里实现,。1:python窗口2:输入代码。3:输出结果。也可以使用pycharm实现。下面是编写弧度制转为十进制的arcpy脚本。import mathimport arcpyarcpyprint.原创 2022-04-27 15:14:17 · 2377 阅读 · 0 评论 -
arcpy 验证中心点是否位于图层之内
for attr in dir(mxd):... if not attr.startswith("_")and hasattr(mxd,attr):... print(attr,getattr(mxd,attr))... 查看地图文档的所有属性。(‘activeDataFrame’, <DataFrame object at 0x2581d690[0x25832b60]>)(‘activeView’, u’\u77e2\u91cf\u5730原创 2022-04-26 18:04:24 · 535 阅读 · 0 评论 -
arcpy 查看arcgis的产品类型、产品许可与安装信息
1:打开地理处理下的python窗口。2:输入import arcpy ;arcpy,输出结果跟我的一样,就代表已经导入了arcpy模块。输入上述代码可以知道安装的产品类型。desktop包括:arcmap arccatalog arctoolbox输入上述代码可以得到desktop的许可类型arcinfo(最高行政级别的许可)。infos=arcpy.GetInstallInfo()>>> infos=arcpy.GetInstallInfo()... for key,.原创 2022-04-26 14:08:55 · 1380 阅读 · 0 评论