根据要素创建要素类
//根据模型构建要素类n private IFeatureClass PRV_FeatureClassBuild(IFeature IN_Feature,string IN_FeatureCLassName,IFeatureDataset IN_FeatureDataset)n n //作为搭建模型的要素n IFeature
Category: ArcPy
Script & Tool之Appenddata
# -*- coding:gb2312 -*-n#Author:tony.gu email: tony.gu@gisinfo.comn#function: append data from one workspace to anothern#目前只支持shp到数据库中的某个数据集nimport arcpy,osnimport fnmatchnfrom arcpy
ArcGIS Python脚本之Arcpy数据访问模块
目录nn一、使用SearchCursor检索要素类中的要素nn二、使用几何令牌改进游标性能nn三、使用InsertCursor插入行nn四、使用UpdateCursor更新行nn五、使用UpdateCursor删除行nn六、在编辑会话中插入和更新行nn七、读取要素类的几何信息nn八、使用Walk()遍历目录nn 游标是包含表格或要素类中的一行或多行数据的内存对象。每一行都包含数据源中每个字段…
ArcGIS Python 获取GDB中的要素类信息
我们经常会有这样的需求，想看下数据库中有多少个数据集多少个要素类，如果同数的方式会发现比较麻烦，那么我们可以通过一小段python脚本完成这样的任务，可以在桌面集成的窗口中完成一可以在单独python环境中完成只是没有太大的区别，下面是我实验的一段代码仅供参考。rn#coding=utf-8rnimport arcpyrnimport osrnarcpy.env.workspace = ‘E:\d
【Arcpy】Arcpy核心
易智瑞（中国）技术支持中心 穆天龙nn
Arcpy 学习笔记1 – 加载栅格数据并显示在数据框中
本文使用的是ArcGIS pro 1.4中的Arcpy，基于python3.4，支持python3语法。 0. 导入Arcpy包 >import arcpy as ap 1. 加载栅格数据，并显示在数据框中 从本地文件中加载栅格数据，这一步将栅格数据加载到缓存中，并不会显示在数据框中。
通过Python脚本自动发布ArcGIS服务
ArcGIS的强大之处就在于其为我们提供了多种选择，我们可以通过Python编写脚本，调用ArcPy模块中的功能来实现自动扫描文件夹并将其中的地图文档发布为地图服务。
Python与GIS之arcpy基础代码详解（一）
最近在学习arcpy，参考书籍为《面向ArcGIS的Python脚本编程》，把书中的重点Python脚本代码整理了一下，仅供参考。nn第一部分为访问空间数据，脚本代码包括检查某个数据集是否存在，使用Describe函数描述数据集的属性，使用列表函数处理不同类型的元素，包括工作空间、字段、数据集、要素类、文件、栅格、属性表等。
发布要素地图服务（Feature Layer）
软件环境：rn操作系统：win7 64位rnGIS平台：ArcGIS 10.1rnGIS数据服务：ArcSDE for SQLSERVERrn数据库：SQLSERVER 2008rnrn如何实现要素的编辑并且保存到服务中，最简单的一种方式就是发布要素服务（Feature Layer），要素服务添加到地图文档的图层和表格通过服务提供。以下将介绍要素服务数据的要求：rnrn所有数据必须来自单个 Arc
ArcGIS教程：模糊隶属
摘要rn 根据指定的模糊化算法，将输入栅格转换为 0 到 1 数值范围以指示其对某一集合的隶属度。rn 值 1 表示完全隶属于模糊集，而当值降为 0 时，则表示不是模糊集的成员。rn rnrnrn 用法rn · 此工具无法对分类数据进行度量。要将分类数据用于模糊叠加分析，需要执行预处理操作。您可以创建一个模型或运行下列地理处理工具。首先，使用重分类工具获得一个新的值域范围(例如，1 到
Python工具分享
rn ArcPy Mapping脚本工具rn 版本：rn ArcGIS 10rn 介绍：rn 这些工具举例说明了使用arcpy.mapping能够执行的各种任务，这个脚本工具分为3大类：制图工具、输出和打印工具、.mxd和.lyr管理工具。几乎所有这些工具都可以在ArcCatal
【短篇】Python下如何把带有千位分隔符的浮点数字符串形式转化为数字
ok，标题写的有点晦涩，其实意思很简单就是如何把rn”、“12,345.678”rn
python arcpy检索ftp中的文件
import ftplibnimport osnimport socketnnHOST = ‘ftp.nifc.gov’nDIRN = ‘/Incident_Specific_Data/2012 HISTORIC/ROCKY_MTN/Arapaho/GIS/20120629’nFILE = ‘20120629_0600_Arapaho_PIO_0629_8x11_land.pdf’nntry:n …
arcpy站点包中常用方法与函数
http://blog.csdn.net/alexbnlee/article/details/6965276
【python】创建目录/文件夹
python可以使用mkdir()和makedirs()两种方法来创建目录，mkdir()创建一层目录，makedirs()创建多层目录。nn创建一层目录：nnn#引入模块nimport osn#创建目录’D:\xxoo’nos.mkdir(‘d:\\xxoo’)n#创建目录’D:\test’nos.makedirs(‘d:\\test’)nn这里面使用mkdir()和makedirs()两种方法…
ArcGIS Python编程案例(2)-使用ArcPy编写脚本
n n n 第一章 Python语言基础 ||| 第三章 管理地图文档和图层nn本章将涉及以下案例：nn使用ArcGIS Python窗体n访问ArcPyn从脚本中执行工具n使用ArcGIS桌面软件帮助系统n使用变量存储数据n访问ArcPy中的模块nnn引言n地理处理任务往往是耗时且重复性的工作，通常又需要定期执行。
（二十八）arcpy开发&关于arcgis中的arcpy.da.Walk
现在我们来学习一下关于arcpy.da.Walk这个api，首先呢，该函数需要提供三个参数，分别是工作空间、数据类型、类型。现在我们测试一下下面这段代码，注意使用的过滤条件。设置数据类型、类型都是any。nnn#coding=utf-8nimport arcpynimport osnnws = r”D:\Data\中国国界和省界的SHP格式数据\省界\demo.gdb”nwalk = arcpy….
ArcGIS for Desktop 图层符号化
添加至ArcMap中的数据是以图层的形式来展现，样式单一且颜色随机，往往不能满足实际工作中出图的要求。因此，我们通常需要对图层的显示样式进行设置，包括颜色、线型等等，这称为图层符号化。 图层符号化设置是图层的属性，因此在图层的属性对话框中进行设置。右键单击图层，在弹出的菜单中选择“属性Properties”，就打开了图层的属性对话框。然后选择“符号系统Symbology”选项卡。 图层符号化…
一个使用C#和ArcPy实现的版本压缩工具（一）
目录ArcGIS版本压缩功能设计与开发传统做法预期功能界面设计nArcGIS版本压缩功能设计与开发nArcGIS版本化数据管理机制中，增量数据最初是保存在 增量表 即A表和D表中的。当我们压缩版本时，增量表中的数据才会进入到基表（对归档后的数据而言，实质上是视图，只不过给用户的感觉还是一个数据表）。n传统做法n我们固然可以使用ArcGIS Engine去实现版本压缩功能，ArcGIS Engine…
python 读取各种 文件
使用mapPartitions()来重用解析器读取 非结构化的pythonimport jsonndata = input.map(lambda x : json.loads(x))n读取CSV 文件 用pythonimport csvnimport StringIOn…ndef loadRecord(line):n “””解析一行csv记录”””n input = StringIO.Stri…
ArcPy添加字段
#添加字段（如存在，则跳过。注意区分大小写）nnfile_names=(‘BOUAN’,’BOULK’,’BOUNT’,’BOUPT’,’CTRLK’,’CTRPT’,’HYDAN’,’HYDLK’,’HYDNT’,’HYDPT’,’PIPAN’,’PIPLK’,’PIPNT’,’PIPPT’,’RESAN’,’RESLK’,’RESNT’,’RESPT’,’ROAAN’,’ROALK’,’R
[Python] 坐标转换|火星坐标与WGS84
# -*- coding: utf-8 -*-nimport jsonnimport mathnnx_pi = 3.14159265358979324 * 3000.0 / 180.0npi = 3.1415926535897932384626 # πna = 6378245.0 # 长半轴nee = 0.00669342162296594323 # 扁率nndef wgs84togcj02…
要素外接矩形的四个角点坐标、长度、宽度、面积如何计算到要素属性表中？
使用工具 Calculate Field。nnnnnnn9.3以上版本可按下图的方法分别求出envelope的x、y坐标的最大、最小值，由此可得出你要的那些值，需注意的的是求长度、面积等要在投影坐标系下进行。nnnX的最大值：!shape.extent.XMax!nX的最小值：!shape.extent.XMin!nY的最大值：!shape.extent.YMax!nY的
在MXD空间中从文件名中提取行政村名称（用于定义查询）
代码：#encoding=gbkn#在python中正确，在arcgis中出错nimport re,os,stringnmz = ‘街道’nfor root,dirs,files in os.walk(r’O:\arcgis_directory\MXD\RHYXT’):n for afile in files:n if afile[-4:].lower() == ‘.mxd’:n new_…
Python自动发布Image service的实现
使用Python自动发布地图服务已经在上一篇博客中讲到，使用Python创建.sd服务定义文件，实现脚本自动发布ArcGIS服务，下面是利用Python自动发布Imagern service的实现。rn———–华丽的分割线———–rn# -*- coding: utf-8 -*-n# ——————————————-
arcpy将shp属性表导入到数据库，sqlserver
由于表记录太长，导出到txt会出现乱码，无法用工具导入，且记录会莫名其妙的减少，shp导入postgis也会出问题nn因此写代码来导入吧…nnn# encoding: utf-8nimport pymssqlnimport arcpynimport osnfrom arcpy import envnimport numpy as npnimport pandas
Anacoda安装—python3.6
1.anaconda官网下载win64版本 nn傻瓜式安装 点击下一步nn记住勾选环境配置nn nn安装完后nn打开Anaconda prompt n输入命令 conda list显示已经安装库nnumpy pandas等重要库nn如果Anaconda官网下载比较慢时，建议采用清华镜像下载n点击Jupyter Notebook 提供浏览器编译环境n新建点击new按钮n复杂调试还需要IDEn还想安…
道克拉斯-普克轨迹抽稀算法
最近项目有这样的一个需求，国土局工作人员需要拿着我们开发地图应用app，围绕某个地块走一圈，然后将轨迹记录下来，形成一个面，填写相应的属性。我们知道在使用GPS定位的时候，获取定位点是密密麻麻的，现在我们需要将定位点数据抽稀。然后将点绘制在arcgis for android开发的地图应用上，于是网络上寻找了相关的资料，使用的是道克拉斯-普克抽稀算法，我们来看一下网络相关该算法的资料。nn主要是依…
Python 实现ARCGIS中的ExtractByMask
1.单个文件的extractimport arcpynfrom arcpy import envnfrom arcpy.sa import *arcpy.env.workspace = “E:\\Data_Ex\\FJJXtif”narcpy.env.extent = “E:\\Data_Ex\\FJJXtif\\20006.tif”narcpy.env.snapRaster
Python求多幅图像栅格值的平均值
本程序所采用的方法并不是最优方法，ARCGIS已经提供了相关的函数供调用。本程序仅供参考。n程序说明： n 文件夹E://work//EVI_Data_tif中存放的是某地区2000-2010年的EVI图像，其中每个年份共13幅。目的是将每年的13幅图像的每个栅格相加求均值，生成相应年份的tif。例如，将2000年的13幅图像相加求均值生成2000.tif，里面的每个栅格的值就是13幅
[ArcPy] Excel转shp文件 | ArcPy与要素类
Python使用xlrd解析Excel，提取信息，保存要素类
ArcPy脚本——根据指定关键字批量导出
nnimport arcpy, osnn#输入文件夹nworkspace = arcpy.GetParameterAsText(0)n#查询条件，例如[系统级别]=’一级’nwhereclause=arcpy.GetParameterAsText(1)n#输出数据的后缀…
ArcMap中使用Python遇到中文字段名错误的解决方法
前一段时间遇到一个客户需要在ArcMap中批量的通过等高线图生成DEM，但是等高线图中的高程字段名称全部是中文，大家都知道Python遇到中文就会出现各种诡异状况，于是我被用户找去解决问题。rn先看看Python的版本，是2.6.5，比用户的2.5.1要高，不知道这个问题能不能在自己电脑上重现。rnrn好了，创建一个ChineseTest.py，内容如下：rnrn测试一下，结果：rnrn好吧，悲剧
GIS批量合并文件夹里的shp文件
最近需要大批量地把一些shp文件合并起来，GIS自带的模型构建里的迭代器功能并不能对文件夹里的所有shp文件进行迭代，并且放入合并工具需要的数据集中，所以想到用Python来实现。由于是Python初学者，所以在网上找了一个怎么批量合并的博客看，发现源代码打进去并没有用。。然后又查了一些别的博客，再自己试验了一下，总结出这个比较简洁的脚本文件，计算的速度比较快。利用了GIS里自带的参…
ArcGIS中python对图层固定边界字段进行插值、渲染并输出图片
最近项目中有需要对shp中某一字段进行自动进行出图，在平时制作相应的模型然后输出python或发布成地图服务即可，但是因为涉及到出题出图要具体边界及渲染，功能实现比较简单，代码尽管各处都有可以作参考，但是还是都是零散的，还得花点时间去查找和重新组织，这里楼主对python代码进行组织和整理，希望能对有需要的人提供帮助。其他解释不多说，直接上代码：rn# -*- coding: utf-8 -*-r
ArcGIS 二次开发系列之Python版示例代码——矢量栅格数据批量投影转化
说明：本示例代码使用最简单的python语法，简要介绍如何使用Python做ArcGIS二次开发。适合初学者阅读nn# -*- coding: utf-8 -*-nnimport arcpynInput_folder = “H:\\Original\\qhnu\\zx20170318\\MAP”nfileExtent = [“.shp”,”.tif”,”.tiff”,”.img”]nnimpo
导入要素类到要素数据集当中（C++）（ArcObject）史上最快
n //导入要素类到要素数据集n IDatasetContainerPtr ipDatasetContainer = ipFeatureDataset;n IDatasetPtr ipInDataset = ipFeatureClass; //ipFeatureClass为同一工作空间下的要素类;n ipDatasetContainer->AddDataset(ipInDa…
在ArcGIS中创建Python工具（二）
上一篇中我们了解到有两种方式在 ArcGIS 中创建 Python工具，这一篇就来看看如何在标准工具箱中创建脚本工具。ArcGIS Help 中脚本工具的帮助过于枯燥，在这里，我以一个具体的实例来总结构建脚本工具的过程，我要实现的需求是做个快速实现羽化边界效果的小工具，预期得到如下的效果：上面效果在ArcMap中完全可以手工运行几个工具实现，但是过程稍微繁琐，那么需求来了，如何做个自定的一键生成羽化
arcgis python get raster extent and set envs.extent (环境设置中的范围)
在Arcgis python函数操作前想要设置范围可以通过arcpy.envs.extent实现，对于栅格来说，获取范围的方法是定义栅格变量，获取其范围。nnnimport arcpynnelevRaster = arcpy.sa.Raster(‘C:/data/elevation’)nmyExtent = elevRaster.extentnnprint myExtentnnarcpy.envs…
ArcEngine调用cmd执行python
介绍经常使用ArcGIS的小伙伴估计经常需要用到ArcToolBox，这个工具箱提供了丰富的工具为我们进行数据处理等操作，但是我们在程序中如何使用这些工具呢？Esri提供了GP工具为我们执行这些操作，但是使用过GP工具的同学都知道这个接口并不是特别的好用，而Esri也主推大家使用Python去处理我们的数据，下面介绍如何调用cmd执行我们的python.如果你不知道python如何写，Esri的官网
Arcpy神器之NumPyArrayToFeatureClass
Arcpy神器之NumPyArrayToFeatureClassrn之前在“ofo数据获取&坐标转换&生成shp数据”一文中涉及到用Arcpy创建点要素数据，最近再次使用时发现怎么那么麻烦rnrn于是有另辟新路，翻阅Arcpy接口文档，终于找到了一个神器NumPyArrayToFeatureClass。
（一）arcpy开发&利用arcpy在arcgis中批量裁剪影像
最近有这么一个需求：我有一个面文件数据，该面文件shapefile有多条记录，具体如下图所示。现在需要遍历出图层中的每个面记录数据，然后作为范围面和一个大的tif影像裁剪。nnnn我们知道，在arcgis中可以使用裁剪工具对影像裁剪。具体步骤，依次选择【Data Management Tools】、【Raster】、【Raster】、【Raster Processing】、【Clip】，最后弹出如…
安装的python导入arcpy
/usr/bin/env pythonn#coding=utf-8nimport archooknarchook.get_arcpy()nimport arcpynfrom arcpy import envnfrom arcpy.sa import *nn# Check out the ArcGIS Spatial Analyst extension licensenarcpy.CheckOu
Arcpy基础入门-6、arcpy要素属性操作
空间数据具有三大特征：空间特征、属性特征和时间特征，同样组成一个几何要素的对象也包括空间几何、属性以及时间表达.n这一节我们以一个实例来说一下arcpy对几何对象属性的操作，属性的操作和数据库的操作基本相同，包括查询、添加、更改和删除，添加和删除主要是针对整个属性添加或者删除一个字段，查询和更改是针对一个要素的属性进行操作，我们来先看一个实例：nnimport
Python中如何解决ImportError：No Module Named arcpy问题
Python中如何解决ImportError：No Module Named arcpy问题rn2015-05-11 09:23本站整理浏览(44)rnPython是ArcGIS的脚本语言，ArcGIS安装时会安装ArcPy，即以成功的arcgisscripting模块为基础并继承了arcgisscripting功能进而构建的站点包，其目的是以高效实用的方式通过
arcgis的arcpy写入几何怎么创建一个空心面要素并读取
nn做这个肯定要参考官方文档，下面附上两个链接，是arcpy的读取几何和写入几何。nnhttp://pro.arcgis.com/zh-c…
（十八）arcpy开发&利用arcpy实现在arcgis中对字段的删除
本节，我们将带领大家一起来学习一下关于字段删除方面的函数。我们将会删除字段中除了FID shape外的其他字段。其中利用到函数字段列举函数。arcpy.ListField,下面定义的这个函数将返回所有的字段。nnndef getFieldNames( fc):nn fieldList = []n for f in arcpy.ListFields(fc):n …
（二十九）arcpy开发&利用arcpy在arcgis中实现批量计算图斑四至坐标（polygon顶点计算，范围获取、pyhon类创建）
一、问题与解决方法nn这个方法的灵感来自昨天昨天，下面来说一下实现功能。nn一直以来都有对图斑(polygon)多边形四至坐标提取的需求，其实是这样的，之前在群里看到了林业部门的一个小伙伴有这样的需求，就是提取图斑(polygon)多边形最四至坐标，即一个多边形的最北点、最东点、最西点、最南点。一直以来，我也想着解决这个问题，都没有想到解决的办法，因为在多边形(polygon)中，arcgis的a…
[GIS原理] 8 GIS基本空间分析-叠置分析|缓冲区分析|窗口分析
文章目录叠置分析(overlay processing)问题背景传统的叠置分析–透图法叠置分析矢量数据的叠加分析点+多边形叠置线+多边形叠置多边形+多边形叠置栅格数据的叠加分析空间逻辑运算数学运算复合法缓冲区分析(buffer analysis)窗口分析三个要素中心点分析窗口大小与类型运算方式例子rn【空间分析】rnrn是从空间数据中获取有关地理对象的空间位置、分布、形态、形成和演变等信息的分析技术rn…
解决GP服务产生的结果无法自动发布为地图服务的问题
在ArcGIS for Javascript API或REST API调用GP服务时，常常会遇到这样一种情况：GP服务运行以后，执行成功，也能够生成结果，然而结果并没有直接产生动态的地图服务供API调用（即使勾选了”view results with a map service”，如下图）。nnnn那么产生这一现象的原因可能为：n1、 GP服务输出的结果路径存在人为干预，也就是生成的结果
[GIS算法] DEM – 种子搜索法替换面属性 – C
【题目】种子搜索法替换面属性rn#include
arcpy更改图层属性
import arcpy.mapping as mappingnmxd = mapping.MapDocument(“CURRENT”)ndf = mapping.ListDataFrames(mxd, “Crime”)[0]nupdateLayer = mapping.ListLayers(mxd,”Crimes2009″,df)[0]nsourceLayer = mapping.Layer(r…
利用arcpy和GP构建批量提取字段名的工具
这时就可以利用arcpy和GP构建工具，已达到批量提取、并且输出到excel表格中，这样可以直观的查阅。
ArcPy复制数据
> import arcpyrnrn… import osrnrn… import timernrn… def export_target_fields(path_input,folder_output,file_output):rnrn… if os.path.isfile(path_input):rnrn… file_path_input=os.path
ArcEngine GP 资料（持续补充中……)
ArcGIS Engine中调用GP的两种方法：http://blog.csdn.net/lysc_forever/article/details/7674332怎样执行GP工具(官方教材）：http://resources.esri.com/help/9.3/ArcGISDesktop/dotnet/c4ff8b68-0410-435f-b8e5-682d5cea47cf.htmGeoproce…
ArcGIS mdb数据库写入失败问题
多次GP之后，即使删除了数据图层，其空间占用也不会释放，需要使用catalog的compact工具进行压缩。如果你的数据库选择为mdb，建议使用之后进行compact
Python图像批处理（图像非监督分类）
本文为JPG格式图像非监督分类，可根据ListRasters选择其他格式进行处理，输入文件路径为”E:/aaa”，输出文件路径为”E:/bbb”
windows系统下给python安装gdal、ogr的操作
转载请注明作者（独孤尚良dugushangliang）出处：https://blog.csdn.net/dugushangliang/article/details/83055390nn nnogr附带在gdal，所以我们只要安装了gdal就可以用ogr了。nn资料找了好久。https://blog.csdn.net/u014177758/article/details/73250888nn打开c…
如何用anaconda安装不在anaconda环境中的python包
近日想要使用python的分词包，但是使用conda install jieba命令却提示没有找到这个包。但是如果用传统的pip install，会与anaconda中的pip命令混淆。解决方案：1、在 D:\ProgramData\Anaconda3\Scripts 中 把 pip.exe 和 pip-script.py 改成 condapip.exe 和 condapip-script.py …
ArcGIS & python学习：（一）批处理修改文件名（NDVI数据）
下载的NDVI数据，一共26年，每年有24个文件夹，
