Arcpy与ArcGIS
Arcpy与ArcGIS教程
Chaoying.
这个作者很懒,什么都没留下…
展开
-
arcpy表连接与字段名更改
arcgis基于字段的表连接以及连接后字段名更改原创 2022-06-03 09:16:19 · 950 阅读 · 1 评论 -
Arcpy自动批量化制图:修改图层数据源 修改布局元素 导出布局地图
先上代码# coding:utf-8_import arcpyimport osif __name__ == "__main__": tifdir = r'I:\DInSAR\PJQZ\result' shpdir = r'I:\DInSAR\reference' tifs = [] with open("./filter_list.txt", mode='r') as f: tifs = f.readlines() tifs = [tif.原创 2022-04-14 12:09:48 · 1515 阅读 · 0 评论 -
Arcpy计算字段
# coding=utf-8# 字段分类import arcpyimport os# 获取当前shp的所有字段def getFieldNames(fc): fieldList = [] for f in arcpy.ListFields(fc): fieldList.append(str(f.name)) return fieldListRecalss_area = """def Reclass(area1): if (area1 <原创 2022-04-10 10:16:32 · 1225 阅读 · 0 评论 -
Arcpy读取SHP
#coding:utf-8_import arcpyimport osdef sim(shpFile1, shpFile2): ds1 = arcpy.CopyFeatures_management(shpFile1, arcpy.Geometry()) ds2 = arcpy.CopyFeatures_management(shpFile2, arcpy.Geometry()) e1 = ds1[0] e2 = ds2[0] # 求两个椭圆的相交部分原创 2022-04-08 14:04:23 · 1348 阅读 · 0 评论 -
Arcgis批量渲染图层脚本
脚本文件 render.pyimport arcpydef renderbatch(rlayer): # rLayer arcpy.mapping.Layer rtype = None if rlayer.isFeatureLayer: rtype = 0 elif rlayer.isRasterLayer: rtype = 1 else: arcpy.AddMessage("Reference Layer ca原创 2022-02-03 18:14:45 · 1534 阅读 · 0 评论 -
ArcGIS二次开发法——图层相关操作(添加、修改)
import arcpyimport oswaveHs = os.listdir('./tifs1')waveHs = filter(lambda waveH: waveH.endswith('.tif'), waveHs)waveCs = os.listdir('./contours1')waveCs = filter(lambda waveC: waveC.endswith('.shp'), waveCs)waves = zip(waveHs, waveCs)tifpath = r原创 2022-01-09 15:17:31 · 1424 阅读 · 0 评论 -
leaflet加载ArcGIS for Server地图图层
leaflet 加载ArcGIS for Server地图服务的两种方式基于ESRI leaflet插件var map = L.map('map').setView([37.7614, -122.3911], 12);L.esri.tiledMapLayer({ url: 'https://services.arcgisonline.com/ArcGIS/rest/services/USA_Topo_Maps/MapServer', maxZoom: 15}).addTo(map);原创 2022-01-08 21:19:13 · 1871 阅读 · 1 评论 -
ArcGIS二次开发——Arcpy批量自动化发布地图服务
参考博客相关函数说明见ArcGIS帮助文档import arcpyimport osdef publishServer(mxd): # define local variables wrkspc = 'C:/Users/user/Desktop/SeaGIS' mapDoc = arcpy.mapping.MapDocument(wrkspc + '/mxds' + '/{}'.format(mxd)) con = wrkspc + '/Server' + '/c原创 2022-01-08 22:28:44 · 826 阅读 · 0 评论 -
ArcGIS批量导出图层脚本
该脚本会将ArcMap中的每一个图层单独导出为图片import arcpyimport osdef exportbatch(outdir, export_width, export_height, resolution): # outdir the directory to store the exported pngs mxd = arcpy.mapping.MapDocument('CURRENT') df = arcpy.mapping.ListDataFrames(原创 2022-02-03 18:24:39 · 3561 阅读 · 0 评论 -
ArcGIS二次开发——Arcpy删除图层
import arcpyimport osmxds = os.listdir('./mxds')for mxd in mxds: Mapname = os.path.join(r'C:\Users\user\Desktop\SeaGIS\mxds', mxd) Map = arcpy.mapping.MapDocument(Mapname) df = arcpy.mapping.ListDataFrames(Map)[0] layers = arcpy.mappin原创 2022-01-08 21:05:10 · 1618 阅读 · 0 评论 -
Arcpy创建shp文件
import netCDF4import numpy as npimport arcpyarcpy.env.workspace = './shps/'arcpy.env.overwriteOutput = Trueds = netCDF4.Dataset('./Wave_2.nc')nodes_x = ds.variables['nodes_x'][:]nodes_y = ds.variables['nodes_y'][:]time_str = ds.variables['time_str原创 2021-10-24 21:55:41 · 881 阅读 · 0 评论 -
Arcpy自动化批量制图
import arcpymxd = arcpy.mapping.MapDocument('./demo_bak.mxd')lyrs = arcpy.mapping.ListLayers(mxd)time = arcpy.mapping.ListLayoutElements(mxd, 'TEXT_ELEMENT', 'time')[0]for i in range(121): shp = './shps/{}.shp'.format(i) cursor = arcpy.SearchC原创 2021-10-25 17:01:18 · 1448 阅读 · 0 评论