python存储CSV二维数据——利用arcpy进行要素差分量统计
# coding=utf-8n#要素差分量统计nimport arcpyn#from arcpy import envnnf = open(‘count.csv’,’w’)ntitle = [‘CITY’,’GREENURBAN_SymDiff’,’WATER_SymDiff’]nf.write(‘,’.join(title)+’\n’) nn#全部城市60个
Category: ArcGIS
no model named shapefile
使用arcpy批处理时，经常使用shapefile.reader等包，但直接import shapefile 模块却报错。原因是shapefile在模块pyshp中，因此导入包pyshp后，import pyshp即可。
Arcpy基础入门-5、读取自定义格式xml
arcpy作为一种脚本开发语言，其功能的目的是为用户提供强有力的可扩展性。n这一节通过示例，介绍比较下使用python和arcpy读取自定义格式的xml并构建要素的方法。n先简单描述下自定义的xml格式：njnnn自定义了一种北京54坐标系的航迹数据，每一个航迹包括了属性和对应的集合点。
如何用Arcpy批量提取输入要素和邻近要素的最短距离？
我先用Arcgis点距离分析生成了下面的图表，如何用Arcpy批量提取INPUT_FID和NEAR_FID之间的最短距离？（INPUT_FID有817条数据，每条对应1504个NEAR_FID数据。）求大神指点！ …
arcengine获取gdb中所有的数据，数据名
FileGDBWorkspaceFactoryClass fac = new FileGDBWorkspaceFactoryClass();nn IWorkspace workspace = fac.OpenFromFile(pathname, 0);nnIEnumDatasetName enumDatasetName = workspace.get_DatasetNames(esriDat…
arcgis的arcmap中导入pandas失败，ImportError: No module named pandas
转载请注明作者（独孤尚良dugushangliang）出处：https://blog.csdn.net/dugushangliang/article/details/83015364nn nn打开arcmap的python窗口，试着导入pandas，如下图：nnnn这表明pandas这个模块没有被找到。nn我曾多次尝试把anaconda和arcgis联系在一起。找了好多blog，做了好多尝试，遇到…
ERROR 000824: The tool is not licensed.
接上一篇arcgis 10.4，在python 中可以导入arcpy成功后还不能保证你代码中调用arcpy函数能够使用，这里涉及到你的arcgis server 是否安装正确，上一篇博客的安装中arcgis server 可能说的不够清楚，这里补上。
ArcGIS中使用python实现:将面要素类的所有节点信息存储在属性表的字段中
总结两种方法实现：将面要素类的所有节点信息，按节点顺序存储在属性表的字段中一种是使用字段计算器的python代码块，在网上找到了实现代码，但在运行中报如下错误，最后确定原因在于所添加的字段长度不够导致的:首先需要通过add field添加一个string类型的字段，注意：字段长度必须足够大（如果面的节点很多），修改如下截图中Length值。然后右键字段打开字段计算器，勾选显示代码块，如下截图所示：…
SDE open rastercatalog，遍历rasterCatalog下的所有RasterDataset
IRasterWorkspaceEx pRasterWex=Utilities.OpenSDERasterWSEx();//访问SDE 中的 RasterCatalog//连接SDE 中的RasterCatalog public static IRasterWorkspaceEx OpenSDERasterWSEx()n n ESRILicense();//L
table总结insertRow、deleteRow
rnrn表格有几行: var trCnt = table.rows.length; rn(table为Idrnrn)rnrn每行有几列：for (var i=0; i
镶嵌数据集工具小结（四）镶嵌数据集的属性 Ⅰ
修改镶嵌数据集类
python3.6 pycharm anaconda3基础名片管理
n# 名片管理功能nprint(‘#’*50)nprint(‘名片管理系统V8.1’)nprint(‘#’*50)nprint(‘1 添加’)nprint(‘2 删除’)nprint(‘3 修改’)nprint(‘4 查询’)nprint(‘5 显示所有的名片’)nprint(‘6 退出’)n# 用来存储名片ncard =[]nwhile True:n num = int(input(‘请输…
欢迎使用CSDN-markdown编辑器
arcgis的python环境搭建若干问题-arcpy是用哪种python编写的？ n Cpython nhttps://geonet.esri.com/thread/50292
arcpy批量转出为kml文件
n# edfncoding: utf-8nfrom __future__ import divisionnfrom sklearn.cluster import KMeansnimport numpy as npnimport arcpynimport pandas as pdnfrom arcpy import envnfrom sklearn.cluster
python 常见问题汇总（待续）
1. No module named skimagepip install scikit-image –upgrade2. No module named dicomsudo pip install pydicom3. python name ‘os’ is not definedimport osThis will import the python’s module os, which app
ArcGIS栅格数据SetNull IsNull Con函数
当两个栅格进行叠加，有时会有一部分没有数据，即用identify点击该区域，Value为NoData，而不是像其他非空区域一样有值。如果不进行处理的话，nodata数据跟其他图层做运算是会变为nodata。所以nodata区域要赋予0值，因为nodata任何数=nodata，因此当需要对栅格空值数据赋值时，要采用条件查询函数将NoData的地方进行赋值。以下以arcgis10.0为例：（1）方法是…
Python在FME中的使用（一）——启动和关闭Python脚本
主要包括： FME中Python使用，包括使用简介，连接Oracle数据库、引入arcgis中python处理地理数据的模块（arcgisscripting 或arcpy）等的简单用法。首先说一下FME中Python的简单用法。nn一、关于帮助文档nn我们可以…
python遍历查找文件夹中符合要求的文件
python遍历查找文件夹中符合要求的文件n样例1：n#! /usr/bin/env pythonn# -*- coding: utf-8 -*-nn# Name: searchfile.pyn# Description: Traverse to find the required files in the foldern# Requirements: Spatial Analyst Extens…
105.忙碌：arcpy
n n n arcpy原来这么好用，arcmap的功能很是强大，那么多工具包，当然对shp的更改可以保存到文件，不用担心格式支持的问题。n虽然arcpy只支持10.2和python2的功能，但是语法并没有差太多。用起来上手很是有意思，但是自己的进度有点落后，所以明天就一定要抓紧进度。n当然什么事是最重要的，什么事要先满足，但是满足的时间不一定要那么长。
连接SDE空间数据库——SQL Server篇
一、软件配置rn要从 ArcMap 连接至 SQL Server 中的数据库或地理数据库，请在 ArcMap 计算机上安装 SQL Server native client 或 Microsoft ODBC driver 11 for SQL Serverrnrnrn二、手动链接rn打开Catalog，选择数据库连接下的添加数据库连接。rnrn在弹出的设置窗口中，设置如下参数：rn（1）数据库平台…
安装Anaconda3 后，怎样使用 Python 2.7？
conda create -n py2 python=2.7activate py2activate py2
Python地学分析 — GDAL批量计算MODIS NDVI
欢迎关注博主的微信公众号：“智能遥感”。nnnn该公众号将为您奉上Python地学分析、爬虫、数据分析、Web开发、机器学习、深度学习等热门源代码。nn本人的GitHub代码资料主页（持续更新中，多给Star，多Fork）：nnhttps://github.com/xbr2017nnCSDN也在同步更新：nnhttps://blog.csdn.net/XBR_2014nn“ NDVI是遥感领域最常…
‘ascii’ codec can’t encode characters in position 0-2: ordinal not in range~~
使用百度的api写了一个翻译shp属性表里韩语的脚本，python2.7需要进行编码的转换，跟着网上的操作做了之后发现整个程序好像不太正常了，看了这篇博客之后才发现“import sys ”之后还有别的注意的地方nn原文——>>>>reload(sys)后print失效问题解决nnpython版本: python2.7.6nn#查看python默认编码格式nnn>&g…
安装 PythonWin
Python IDE
使用python进行合并数据集
1、对问题进行描述，手里有两张表，第一张表是手机号码归属地的数据库，里面包含手机号，手机号码前七位，手机号的归属地省份，手机号归属地城市以及运营商等；另一张表中只有手机号和手机号前七位，现在要匹配第二张表中的手机号的归属地 n2、用python进行实现#####导入归属地数据库import pandas as pd nimport numpy as np nfr=open(r’C:\Users\so
arcpy。python创建zip文件
import osnimport zipfilenn#create the zip filenzfile = zipfile.ZipFile(“shapefiles2.zip”, “w”, zipfile.ZIP_STORED)nfiles = os.listdir(“c:/ArcpyBook/data”)nnfor f in files:n if f.endswith(“shp”) or …
地形图的分幅与编号(一)
地形图的分幅与编号n地形图的分幅方法有两种：一种是经纬网梯形分幅法或国际分幅法；另一种是坐标格网正方形或矩形分幅法。前者用于国家基本比例尺地形图，后者用于工程建设大比例尺地形图。n一、经纬网国际分幅法：n1) 1∶100万比例尺地形图的分幅和编号n1∶100万地形图分幅和编号是采用国际标准分幅的经差6°、纬差4°为一幅图。从赤道起向北或向南至纬度88°止，按纬差每4°划作22个横列，依次
ArcGIS对未知坐标系要素建立拓扑查错修错
ArcGIS的拓扑工具很强大，城市规划使用AutoCAD绘制的重叠、间隙等错误百出dwg文件在导入ArcGIS后就可以建立拓扑来自动查错，批量修错，非常方便。n文件地理数据库下面的要素数据集对数据的规范程度似乎要求很高，我一直没搞明白的是，即便新建同样的未知坐标系的要素数据集，要素类还是偶尔能导出进来（貌似只有第一个可以导出过来），很多时候导入报错。n最近突然想到不如直接在要素数据集里新建空白
迭代要素类进行批量添加字段和面积
转载自：https://blog.csdn.net/sinat_25006091/article/details/83040111
最详细的ArcGIS工具集制作流程（普通工具集和python工具集）
1概述n本文旨在介绍如何在ArcGIS pro中如何构建toolbox或者python toolbox,下面将分别从这两个工具创建进行介绍，每一步都有详细截图，方便查看学习！n2 普通toolbox创建n在Catalog面板下面的toolboxes，选中右键，选择 New Toolboxnn选择存放的路径nn选中刚刚创建的工具集，右键Newnn从上图我们可以看出、普通工具集可以通过您的模型构建器、…
缓冲区分析
n n n 一.距离制图－创建栅格缓冲区n软件准备：ArcGIS（10,0以上的版本）数据准备：图层文件point.shp，lline.shp，polygon.shp(http://pan.baidu.com/s/1o7CczYa)1.1点要素图层的缓冲区分析1)在ArcMap中新建地图文档，加载图层：StudyArea，point2)打开工具[Spatial Ana…
python实现栅格数据最大值合成法
哈哈，先用arcpy实现吧，有余力了再自己写（估计猴年马月了）。n根据多个栅格数据计算每个像元的统计数据。可用的统计数据有：众数、最大值、均值、中位数、最小值、少数、范围、标准差、总和及变异度。n在本示例中，将针对多个输入 Grid 栅格中的每个像元计算标准差，并将结果输出为 Grid 栅格。nCe…
创建要FeatureDataset素数据集
创建矢量要素数据集
py2exe打包的简单使用
一. 安装py2exenn1. 安装Python2.7.9, 配置好环境变量n2. 安装py2exe(可以自动检测到Python的安装路径)n下载地址nhttp://sourceforge.net/projects/py2exe/files/py2exe/nnn二. 打包步骤nn2.1 建一个目录, 例如 py2exeTestn2.2 新建如下两个文件n#-*- c
栅格数据创建与保存
栅格数据创建与保存nn作者：阿振nn邮箱：tanzhenyugis@163.comnn博客：https://blog.csdn.net/theonegis/article/details/80089375nn修改时间：2018-05-24nn声明：本文为博主原创文章，转载请注明原文出处nnnnnn思路与方法nn使用Python进行栅格数据处理，很多时候，我们会将GDAL的Dataset对象转化为N…
python脚本交换要素XY坐标
在实际工作中你可能会遇到这样的问题，将空间图形的XY坐标进行反转，如下图所示：rn三种类型的要素都可以完成XY坐标点的互换操作。rn整个处理的思路为：rn1.获取输入要素，然后逐个进行遍历rn2.同时新建空间类型相同的要素类，然后创建要素rn3.将遍历要素中所有节点的XY反转后赋予新要素（注意点、线、面的操作略有不同）rn4.保存新建要素，并输出结果。rn参考代码如下（在arcgis9.3环境下测
Python实现批量执行ArcGIS版本的协调和提交 – ArcGIS技术研究 – 博客频道 – CSDN.NET…
http://t.cn/zls43hc
python三维模型批量导入
import osrnimport os.pathrnrootdir=”D:\test”rnoutputFC=”C:\Users\yao\Documents\ArcGIS\Default.gdb\tt”rninputfeature=[1]*100rni=0rnfor parent,dirnames,filenames in os.walk(rootdir):rn for filename
shp加载并合并
private void simpleButton4_Click(object sender, EventArgs e)rn rn //创建一个文件打开对话框来打开源文件rn OpenFileDialog pDlg = new OpenFileDialog();rn pDlg.Filter = “Shape(*.sh
python_进行栅格计算_加法_批量处理
nimport arcpy nfrom arcpy.sa import * narcpy.CheckOutExtension(“spatial”) narcpy.gp.overwriteOutput=1 n#定义输出路径 noutpath=”e:/ccc/” noutfilename=”mean”n#更改工作空间将，将”e:/bbb”更改为处理cdf文件时的输出路径
python用updatecursor删除行
arcpy.env.workspace = “C:/ArcpyBook/Ch8/WildfireData/WildlandFires.mdb”nn… try:nn… with arcpy.da.UpdateCursor(“FireIncidents”,(“CONFID_RATING”),'[CONFID_RATING] = \’1\”) as cursor:nn… c…
利用Python实现ARCGIS中栅格计算器的简单功能
1.加法操作import arcpynfrom arcpy import envenv.workspace = “E:\\work\\cppservice\\RasterCalc” #工作空间可以设置也可不设置ninRaster1 = “E:\\work\\cppservice\\finalResult\\propertyGower.tif”ninRaster2
数据转换：Excel转shp时常见的问题（2）
在Excel转shp过程中，发现打开【从XY表创建要素类】对话框后，无法获取数据的X字段和Y字段，如下图所示。但自己在Excel中新建一个字段却可以读取，这是什么原因呢？rn rn在Excel中打开表格，发现在x和y列的数值左上角有个绿色三角形，如下图所示，rn rn这个标志表示当前的值为文本格式，单击数据，单击叹号下拉框。选择【转换为数字】如下图所示。依次将文本数据转换为数字即可。
CSV文件转Shp文件
准备工作n环境：Python3.7，需要引入shapefile、csv两个库n开发工具：Pycharmn示例文件：链接：https://pan.baidu.com/s/16xbjqLMaTiVy2hB8FGUh4wn提取码：ch2dn原始教程网站(Pyshp文档)：https://pypi.org/project/pyshp/nn代码实现nimport shapefilenimport csvnn…
Python图像批处理（图像重采样）
本文为JPG格式图像重采样批处理，可根据ListRasters选择其他格式进行处理，可选择重采样方式，输入文件路径为”E:/aaa”，输出文件路径为”E:/bbb”
一个SDE无法启动的问题以及解决方案
环境： SDE9.3.1+SP2(64bit) OS: Redhat 5.5 X64 DBMS: Oracle10.2.0.4(64bit) 周末，同事将Oracle的plsql_compiler_flags的值改成NATIVE后，再启动SDE的时候会报一些对包没有执行权限的错误，今天来到公司，查看了一下资料，发现从oracle10g开始，plsql_compiler_flags参数已经被oracle废掉了，替代它的参数为plsql_code_type ，我将plsql_code_type 参数改
Python3.6之Windows下Anaconda3和PyCharm 开发环境安装
继上一篇https://blog.csdn.net/qq_30993595/article/details/80344312博客后，初次见识Python开发后，那我们这篇就来总结下Python的开发环境，毕竟工欲善其事，必先利其器嘛。nn我这里所用到的Python管理系统是Anaconda，集成开发工具用的是Pycham，这两个是什么先不说，先上张图看看nnnn nnnn nn这两个用起来简直嗨的…
（八）arcpy开发&利用arcpy在arcgis中创建属性表
今天我们来学习一下，关于在arcgis中，使用arcpy来创建属性表。其中使用到的函数为nnCreateFeatureclass_management，在使用该函数创建属性表需要注意的是需要定义好投影，我们来看一下实现的代码。
（十四）arcpy开发&利用arcpy遍历shp文件并获取要素记录数
这是一个小实例程序，首先是遍历某路径下所有.shp文件，然后创建临时图层。然后遍历临时图层获取相应的属性记录数。最后的代码如下所示。nnnimport osnimport arcpyndirectory = r”C:\Users\qrb_PC\Desktop\fast\shp”nnfor root, dirs, files in os.walk(directory):n ct = 0n …
Python/ArcPy遍历指定目录中的MDB文件
#遍历指定目录中的MDB文件，构造FeatureClass名rnrn>>> target_folder = ‘D:\T20161202’rnrn… file_names=(‘BOUAN’,’BOULK’,’BOUNT’,’BOUPT’,’CTRLK’,’CTRPT’,’HYDAN’,’HYDLK’,’HYDNT’,’HYDPT’,’PIPAN’,’PIPLK’,’PIPNT’,’PIPPT’,
ArcPy创建ArcGIS Server连接文件
关于更多的ArcPy的参考，可以查阅ArcGIS的帮助文档，里面都会有很详细的代码示例，只需要修改对应的参数即可。
使用arcpy进行简单的要素编辑
类似上篇博客内容，这次直接对要素属性内容考试编辑rnimport arcpyrnfrom arcpy import envrnenv.workspace=r”C:UsersyaoDocumentsArcGISDefault.gdb”rnprint “修改之前：”rnCur=arcpy.SearchCursor(“地级市”,”回族rnfor row in Cur
Follow Me: ArcPy(1)
下面，让我来对ESRI的新产品ArcPy，进行一番深度吐槽！
csv文件转换为shp(Python实现) python——csv文件转换为shp
python——csv文件转换为shpn 转载自https://blog.csdn.net/qq_23926575/article/details/82316965 n学地信的经常会有这样的需求，即将csv格式的文件转换为shapefile格式加载到arcMap中进行显示与分析，arcMap本身提供了这样的功能，即nn 文件->添加数据->添加XY数据…
arcpy批处理，使用栅格的范围裁剪矢量，并将裁剪后的矢量栅格化
使用arcpy实现，代码如下：import arcpynarcpy.env.workspace = ‘F:/DL_DATA/mass_buildings/test/sat’nrasters = arcpy.ListRasters(“*”, “tiff”)nOutput_Workspace = “F:/DL_DATA/ma…
批量为shp文件加字段并赋值文件名
import sys, string, os, arcgisscriptingrnrnrn# Create the Geoprocessor objectrngp = arcgisscripting.create()rnrnrn# Load required toolboxes…rngp.AddToolbox(“C:/Program Files (x86)/ArcGIS/ArcToolbox/
ArcEngine 清除影像黑边
通常许多的影像有黑边，也就是value=0。通过以下代码可以清除Valu=0的值。 IRasterLayer pRasterLayer = new RasterLayerClass(); n //strFileFullName 输入文件n pRasterLayer.CreateFromFilePath(strF
python读取和生成excel文件
今天来看一下如何使用python处理excel文件，处理excel文件是在工作中经常用到的，python为我们考虑到了这一点，python中本身就自带csv模块…rnrnrn 今天来看一下如何使用python处理excel文件，处理excel文件是在工作中经常用到的，python为我们考虑到了这一点，python中本身就自带csv模块。rnrnrn 1.用python读取csv文件：rnrnr
C#调用python
首先下载最新的IrorPython（点击打开链接），目前最新的是2.7.8。安装。引用安装路径下的DLL（C:\Program Files\IronPython 2.7）新建一个Winfrom程序调用写好的Python代码ScriptEngine pyEngine = Python.CreateEngine();//创建Python解释器对象n ScriptSource so…
ArcGIS10.1正式版全套下载、安装与破解（附ArcGIS10.1全套下载地址）
一、准备文件1.ArcGIS10.1安装包：ArcGIS_Desktop_10.1_129026(en)n2.认证服务：Pre-release_license_managern3.注册机：arcgis10.1KeyGen（据说可以用10.0的代替n nn 以上软件在本文最后一部分均附有下载地址，不能保证长期有效，欲下请从速！n 二、安装步骤n 1.安装Pre-release_l
