arcpy
A873054267
这个作者很懒,什么都没留下…
展开
-
python实现空间数据聚类-DBSCAN
主要思想就是读取坐标值进行聚类,坐标来源包括shp,txt,数据库# encoding: utf-8import arcpyimport numpy as npimport pandas as pdfrom sklearn.cluster import DBSCANoutgdb=r"D:\84.gdb"lsls=[]with open(r"C:\Users\h\Desktop\...原创 2019-05-22 14:55:41 · 3029 阅读 · 6 评论 -
arcpy将csv转换为shp
#和txt的处理方式大同小异#实现思路为列出文件夹中所有的文件,然后读出每个csv中的记录,再添加shp# encoding: utf-8import csvimport arcpyimport oswenjianjia=r"C:\Users\h\Downloads\csv"outGDB=r"d:\csv.gdb"parents = os.listdir(wenjianjia)...原创 2019-02-15 22:57:16 · 2812 阅读 · 1 评论 -
arcpy批量转出为kml文件
# edfncoding: utf-8from __future__ import divisionfrom sklearn.cluster import KMeansimport numpy as npimport arcpyimport pandas as pdfrom arcpy import envfrom sklearn.cluster import DBSCANimp...原创 2019-01-17 09:04:28 · 1842 阅读 · 0 评论 -
arcpy获得图层要素数量
方法1#fc为要素的路径arcpy.GetCount_management(fc)返回值类型是result,要获取具体的数字可以使用arcpy.GetCount_management(fc).getOutput(0)方法2 #cursor遍历的方式。方法3 转为numpy的方式,以后用得着,栅格矢量都可以my_array = arcpy.RasterToNumPyArra...原创 2019-01-22 12:29:37 · 4020 阅读 · 0 评论 -
arcmap实现大批量数据的空间校正,saptial adjust
arcmap对应大批量数据的渲染体验很差,在空间校正中想要从几十万个面块中去定位到单个面块是特别困难的。上一篇博文写到如何使用程序去 进行空间校正,但是对于单纯的使用者来说门槛较高,因此本文采用arcmap自带的方式来进行。首先从几十万个图斑中抽取要校正的单个面块出来,然后与待校正的面块利用工具进行校正。这样可以得到校正文件。然后利用此校正文件来校正原来的几十万个面块。具体的操作过程,...原创 2019-01-15 15:54:23 · 1563 阅读 · 7 评论 -
arcpy spatial adjust更改构成面.线的坐标点
之前写过更改点的坐标点和读取构成线面的坐标# encoding: utf-8import arcpydata=r"d:\JZ\jz\jianzhuclip.shp"featureList=[]point = arcpy.Point()array = arcpy.Array()feacount=0with arcpy.da.UpdateCursor(data, ["SHAPE...原创 2019-01-15 15:45:08 · 884 阅读 · 1 评论 -
arcpy获取组成点,线,面要素的坐标
点要素简单就不说了,线面同理# encoding: utf-8import arcpydata=r"C:\Users\h\Desktop\nb\Export_Output_2.shp"with arcpy.da.SearchCursor(data, ["SHAPE@","FID"]) as cursor: for row in cursor: for pt...原创 2019-01-15 10:59:11 · 3442 阅读 · 0 评论 -
arcmap发布服务卡死,预览界面无反应
原因是由于前两天删除python路径的时候,将server使用的64位python给删除了,只需要将server中带的64位python2.7.3安装在注册表目录下即可,此路径要与注册表下面的路径对应。 ...原创 2019-01-18 19:50:39 · 1877 阅读 · 0 评论 -
arcpy,python读取xml文件
from xml.dom import minidomxmldoc = minidom.parse("WitchFireResidenceDestroyed.xml")childNodes = xmldoc.childNodeseList = childNodes[0].getElementsByTagName("fire")for e in eList: if e.hasA...原创 2019-01-09 21:39:13 · 332 阅读 · 0 评论 -
arcpy。python创建zip文件
import osimport zipfile#create the zip filezfile = zipfile.ZipFile("shapefiles2.zip", "w", zipfile.ZIP_STORED)files = os.listdir("c:/ArcpyBook/data")for f in files: if f.endswith("shp") or ...原创 2019-01-09 21:39:22 · 275 阅读 · 0 评论 -
python arcpy检索ftp中的文件
import ftplibimport osimport socketHOST = 'ftp.nifc.gov'DIRN = '/Incident_Specific_Data/2012 HISTORIC/ROCKY_MTN/Arapaho/GIS/20120629'FILE = '20120629_0600_Arapaho_PIO_0629_8x11_land.pdf'try: ...原创 2019-01-07 15:47:46 · 245 阅读 · 0 评论 -
arcpy valuetable对象的使用,多值参数指定方式
#多值参数指定方式1 python list类型2 字符串类型,以逗号隔开3 valuetable类型import arcpytry: arcpy.env.workspace = r'c:\ArcpyBook\data' vTab = arcpy.ValueTable() vTab.setRow (0, "5") vTab.setRow (...原创 2019-01-07 15:42:46 · 684 阅读 · 0 评论 -
arcpy fieldmap,合并字段不匹配的图层
import arcpytry: # Local variables arcpy.env.workspace = r"c:\ArcpyBook\data" outFeatureClass = r"c:\ArcpyBook\data\AllTracts.shp" # Create a fieldmappings adding the three new fi...原创 2019-01-07 15:40:43 · 835 阅读 · 0 评论 -
arcmap实现最小成本
1 首先根据起点 和成本栅格,使用path cost工具计算出距离栅格和成本回溯链栅格然后再用成本路径选择对应的文件即可计算出最小成本路径原创 2019-02-21 20:05:38 · 708 阅读 · 0 评论 -
arcpy实现 kml批量转出为shp
分两步走:首先转出到layerimport arcpy, osarcpy.env.workspace = r"D:\sq\ddx"outLocation = "D:\sq\ddx"MasterGDB = 'sq.gdb'MasterGDBLocation = os.path.join(outLocation, MasterGDB)arcpy.CreateFileGDB_manage...原创 2019-02-23 13:00:51 · 3132 阅读 · 1 评论 -
arcpy创建gdb,批量删除gdb中的所有文件
# encoding: utf-8import arcpyimport pandas as pdfrom arcpy import envimport numpy as npimport osdef deleteGDBFile(gdbpath): env.workspace=gdbpath fcs=arcpy.ListFeatureClasses() f...原创 2019-02-18 15:24:17 · 2609 阅读 · 0 评论 -
arcpy批量镶嵌
# encoding: utf-8import arcpyimport arcpyfrom arcpy import envenv.workspace = r"C:\Users\h\Desktop\cs\cs"rasters=arcpy.ListRasters()#存储所有的栅格文件名data=[]for raster in rasters: print raster ...原创 2019-05-08 16:27:18 · 1396 阅读 · 0 评论 -
arcpy批量发布影像服务-10.2以后,不用再创建mxd
importarcpyimportosimportarcpy.mappingasmappingdefcheck_PubMXD(imgpath,imgName,wrkspc,service,con):sddraft=wrkspc+service+'.sddraft'arcpy.CreateImageSDDraft(imgpath,sddraft...原创 2019-04-01 22:01:15 · 522 阅读 · 0 评论 -
arcpy将shp属性表导入到数据库,sqlserver
由于表记录太长,导出到txt会出现乱码,无法用工具导入,且记录会莫名其妙的减少,shp导入postgis也会出问题因此写代码来导入吧...# encoding: utf-8import pymssqlimport arcpyimport osfrom arcpy import envimport numpy as npimport pandas as pdserver='l...原创 2019-03-11 10:53:25 · 1425 阅读 · 0 评论 -
arcpy批量近邻分析,不同的道路
# encoding: utf-8import osimport arcpyfrom arcpy import envpath=r"D:\data\polygon"dl=r"D:\data"for fc in os.listdir(path): if(fc.endswith(".shp")): fcname=fc.decode('GBK').encode('u...原创 2019-03-03 14:41:27 · 778 阅读 · 1 评论 -
arcpy 批量将线文件转换为面文件,解决中文乱码
# encoding: utf-8import osimport arcpypath=r"D:\data"for filename in os.listdir(path): filename=filename.decode('GBK').encode('utf-8') print filename if(filename!="polygon"): ...原创 2019-03-03 14:14:35 · 799 阅读 · 0 评论 -
arcpy创建带高程的点图层
def createShpFromTxt(txtPath,outPath): #打开txt文件 allLine = [] with open(txtPath,'r') as file: #按行读取 for line in file: #存储到列表中 allLine.append(line)...原创 2019-03-05 22:02:11 · 745 阅读 · 0 评论 -
arcpy批量复制栅格数据
# encoding: utf-8import arcpyimport osfrom arcpy.sa import *arcpy.CheckOutExtension("spatial")arcpy.gp.overwriteOutput=1raster=r"F:\raster2.gdb"raster2=r"F:\raster2.gdb"arcpy.env.workspace=ra...原创 2019-02-20 20:46:32 · 1134 阅读 · 0 评论 -
arcpy 最短路径
# encoding: utf-8import arcpyfrom arcpy import envtry: # Check out the Network Analyst extension license arcpy.CheckOutExtension("Network") # Set environment settings env.workspace...原创 2019-02-20 15:12:52 · 1490 阅读 · 4 评论 -
arcpy删除GDB文件
此前做的工作是将GDB里面的文件遍历删除再重新创建def deleteGDBFile(gdbpath): env.workspace=gdbpath fcs=arcpy.ListFeatureClasses() for fc in fcs: arcpy.Delete_management(fc) fcs = arcpy.ListTables()...原创 2019-02-25 08:33:34 · 2146 阅读 · 2 评论 -
arcpy批量合并shp
# encoding: utf-8import csvimport arcpyimport osshuju=r"D:\short2.gdb"arcpy.env.workspace=shujufcs=arcpy.ListFeatureClasses()lspt = []for fc in fcs: lspt.append(fc)arcpy.Merge_manageme...原创 2019-02-19 10:59:13 · 2790 阅读 · 0 评论 -
arcpy批量更新shp某一字段的值
#提取文件名更新到dlbh字段# encoding: utf-8import csvimport arcpyimport osshuju=r"D:\short2.gdb"wenjianjia=r"D:\Desktop\csv.gdb"fish=r"D:\Desktop\yw\fishnet.shp"result=r"D:\Desktop\kongjian.gdb"arcpy.e...原创 2019-02-19 09:57:22 · 2408 阅读 · 0 评论 -
arcpy进行地理编码
import requestsimport jsonagisurl = "http://geocode.arcgis.com/arcgis/rest/services/World/GeocodeServer/find"payload = { 'text': '1202 Sand Wedge, San Antonio, TX, 78258', 'f': 'pjson'}r = requ...原创 2019-01-07 15:38:23 · 345 阅读 · 0 评论 -
arcpy查询地图服务
import requestsimport jsonagisurl = "http://sampleserver1.arcgisonline.com/ArcGIS/rest/services/Demographics/ESRI_Census_USA/MapServer/1/query"payload = { 'where': 'STATE_FIPS = \'48\' and CNTY_FI...原创 2019-01-07 15:36:56 · 387 阅读 · 0 评论 -
arcpy缩放至要素的另一种方法
# encoding: utf-8#获取地图文档中的图层列表#mxd可引用网络驱动器上的文档import arcpyfrom arcpy import mappingmxdpath=r"C:\Users\h\Desktop\fp.mxd"mxd=arcpy.mapping.MapDocument(mxdpath)df=mapping.ListDataFrames(mxd,"Cri...原创 2019-01-05 19:58:11 · 894 阅读 · 0 评论 -
arcpy缩放至所选要素
# encoding: utf-8#获取地图文档中的图层列表#mxd可引用网络驱动器上的文档import arcpyfrom arcpy import mappingmxdpath=r"C:\Users\h\Desktop\fp.mxd"mxd=arcpy.mapping.MapDocument(mxdpath)mxd.activeDataFrame.zoomToSelectedF...原创 2019-01-05 19:48:43 · 2173 阅读 · 0 评论 -
arcpy读取地图文档中的图层
# encoding: utf-8#获取地图文档中的图层列表#mxd可引用网络驱动器上的文档import arcpyfrom arcpy import mappingmxdpath=r"C:\Users\h\Desktop\fp.mxd"mxd=arcpy.mapping.MapDocument(mxdpath)layers=mapping.ListLayers(mxd)for ...原创 2019-01-05 19:41:45 · 1784 阅读 · 1 评论 -
arcpy另存地图文档
其中mxdpath不得有中文# encoding: utf-8import arcpymxdpath=r"C:\Users\h\Desktop\fp.mxd"mxd=arcpy.mapping.MapDocument(mxdpath)mxd.title="xuexi"mxd.saveACopy(r"C:\Users\h\Desktop\cs.mxd")原创 2019-01-05 19:37:46 · 698 阅读 · 0 评论 -
arcpy更改点坐标
更改线坐标的教程https://www.cnblogs.com/beautifulplanet/p/4317724.html#reusult为结果图层路径cursor2 = arcpy.UpdateCursor(result)#inRow5[1][0] inRow5[1][1]为X,y坐标 vertex = arcpy.CreateObject("Point")vertex...原创 2018-12-24 15:40:25 · 1410 阅读 · 0 评论 -
arcpy SelectByLocation_management 按位置选择
import arcpyfrom arcpy import envsourceGDB=r"d:\polygon.gdb"env.workspace=sourceGDBtarget="merge"out_feature_class="p2"arcpy.MakeFeatureLayer_management(out_feature_class, 'newlinks_lyr')ar...原创 2018-12-11 13:26:22 · 2603 阅读 · 1 评论 -
arcpy强制升级pip,numpy
可先卸载pip,然后至官网pip下载,用python setup.py installnumpy 可以使用 pip install --ignore-installed numpy命令来强行升级。升级matplotlib,由于此模块3.0后不再支持python2,因此需要指定低于3.0的版本pip install --ignore-installed matplotlib==2.0即可...原创 2018-11-26 15:07:30 · 1021 阅读 · 2 评论 -
python2.7 arcpy安装pandas
网上说的镜像源全部没有成功。下面的步骤不一定要全部尝试,楼主的成功方式卸载pip,更换为pip18版本下载pandas 0.12版本,不要太高,因为依赖于numpy的版本,升级卸载此包均失败,所以只有妥协。使用python setup.py instasll方式安装,exe安装的import就会报错然后安装pandas会提示各种包缺失,只有一个一个在pypi上下载下来再安装...原创 2018-11-18 16:12:12 · 1948 阅读 · 0 评论 -
基于arcpy实现空间数据聚类,kmeans
并不能直接进行空间数据的聚类,原理是读取要素的x,y坐标来进行聚类,然后将聚类中心保存为空间数据以达到效果# encoding: utf-8from sklearn.cluster import KMeansimport numpy as npimport arcpyimport pandas as pdfrom arcpy import envenv.workspace=r"D...原创 2018-11-20 16:43:49 · 956 阅读 · 0 评论 -
arcpy根据txt创建点图层
sr = arcpy.SpatialReference(4326)fc=arcpy.CreateFeatureclass_management( r"D:\cs","test.shp", "POINT", "", "","", sr)arcpy.AddField_management(r"D:\cs\test.shp", "leibie", "TEXT")cursor=arcpy.Inse...原创 2018-11-20 15:00:13 · 1633 阅读 · 4 评论 -
arcpy 连接postresql,Oracle,sqlserver,并获取数据库中的所有数据
通过sde文件链接,sde可在catalog中添加数据库连接,连接成功后,右键属性可以看到sde文件被创建了。在引用此文件即可import arcpyarcpy.env.workspace=r"C:\Users\Administrator\AppData\Roaming\ESRI\Desktop10.1\ArcCatalog\连接到 localhost.sde"fdlist=arcpy.L...原创 2018-11-01 14:33:19 · 1227 阅读 · 0 评论