ArcGIS 10.2.1 for Desktop新特性
1 地理处理nArcGIS 10.2.1的地理处理新增和改进了一些工具，并且新增了ArcPy函数。n1.1 亮点n1) 重命名字段n新的Alter FieldProperties （修改字段属性）工具允许用户修改字段名和别名，输入数据必须是地理数据库属性表或要素类，shapefiles和coverages则不能修改字段属性。
Category: ArcGIS
读取excel 返回某一列 数据的list或者string
package com.ugo.task.util;nnimport java.io.File;nimport java.io.FileInputStream;nimport java.util.ArrayList;nnimport org.apache.poi.hssf.usermodel.HSSFWorkbook;nimport org.apache.poi.ss.usermodel.Cell…
在Arcgis中基于Python对地图分级别进行四色填充
四色填充是数学领域比较有名的定理，大概意思是说对于任意无飞地的多边形区域，总能选取四种颜色对每个多边形进行填充，保证相邻的多边形具有不同的颜色。在地图制图中，该定理被用于地图着色，保证只采用四种颜色而使得每个省/市/县与相邻区域具有不同的颜色。一、项目需求最近因项目需要，研究了地图四色填充算法。 n主要问题：在web制图中，地图被切成矢量瓦片用于渲染和显示，这就失去了地图的拓扑关系，如何来计算不同省
利用ajax请求arcgis发布的wfs地图要素服务
首先在写这篇博客之前，要感谢我的同事们，得到他们的指点，学到了不少的东西。可以说站在巨人的肩膀上，让你看得更高、更远。也荣幸遇到这样的巨人。nn最近开发的这个地理信息系统项目，我负责了三维平台部分。在三维平台中需要接入二维发布的数据，而我们二维地图使用的是esri的arcgis发布的wfs要素地图服务。之所以使用wfs地图服务，而不是使用存储在普通的数据库的原因，是因为如果用于存储在普通数据库中数…
SDE空间索引的内部运行机制2
2.2.1各种矢量数据存储类型的空间索引的实现方式nArcSDE从开始到现在主要支持以下几种存储类型：nnnnn存储类型nnn数据库nnnnnSDELOBnnnOracle,SQLServer,DB2nnnnnST_GEOMETRYnnnOracle,DB2,PostgresQL,InformixnnnnnSDO_GEOMETRY
EXCEL中的空间数据转换为shapefile和GDB中的feature class
在实际工作中，经常遇到需要将EXCEL格式的空间数据转换成shape格式的矢量图形数据文件，或者将其导入到GDB中的某个feature class中进行空间数据的入库工作。本文介绍如何通过使用ARCMAP软件实现上述过程。
ArcGIS ArcPy Python处理数据
1、使用搜索游标查看行中的字段值。import arcpynn# Set the workspacenarcpy.env.workspace = “c:/base/data.gdb”nn# Create the search cursorncursor = arcpy.SearchCursor(“roads”, ‘”TYPE” <> 4′)nn# Iterate through the…
点与多边形的叠加分析
一个点图层与一个多边形图层叠加，计算每个多边形中有多少个点，是一个很长见的操作。但是如果没有意识到点落在多边形边界上的特殊情况，很可能，你会得到一个错误的结果。rnrnrnrn比如看这个情景，多边形代表一些地区，点表示房产，每个点有个hhsize的属性，表示这个房产中有多少人。最后需要计算每个地区中的人数。rn上图，把hhsize符号化表示出来。rnrnrn具体方法很多，但比较“标准”的做法，是使
ArcSDE数据迁移
Arc SDE 数据迁移n题记n这篇文章主要是利用ArcGIS桌面、ArcSDE命令提供的多种ArcGIS数据导入导出数据（可以理解为数据迁移）的方式和数据库迁移的方式（以Oracle数据库为例），来介绍一下ArcGIS不同迁移方式、数据库不同迁移方式的优缺点，而且这些迁移方式也是根据不同的数据量、不同的运行环境、不同的业务需求来设计的。n注意n首先再次强调，迁移之前做好数据备份工作。
巧妙利用图廓进行图幅接边
原文发布时间：2012-08-20作者：怕冷的企鹅无论数据建库的结构有多么的千差万别，一旦涉及到地形图的数据入库，就会不可避免的涉及到数据接边问题。随着信息化、规范化的进程，数据接边也越来越趋向于程序的自动化。本篇文章就针对地形图的数据接边操作，提出一种新的思路和方法。利用内图廓的完成地形图中图幅间面要素基于图形、属性的接边工作的思路。注：本文中只涉及到地形图中的面要素。另外，在接边过程中不考虑容…
arcpy查找，修复丢失的数据源
#查找丢失数据源的图层nnimport arcpy.mapping as mappingnmxd = mapping.MapDocument(r”c:\ArcpyBook\Ch3\Crime_BrokenDataLinks.mxd”)nlistBrokenDS = mapping.ListBrokenDataSources(mxd)nfor layer in listBrokenDS:n p…
arcgis10.0及以上版本，使用arcpy加载在arctoolbox中实现对所有mxd中特定图层的定义查询
代码：#coding:gbkn#author:yangyongn#会报错，不过不影响正常使用n#cunMing表示村落名称，可用于定义查询，注意用%nimport arcpy,os,renimport arcpy.mapping as mappingnglobal cunMingnnimport sysnreload(sys)nsys.setdefaultencoding
arcpy 批量将线文件转换为面文件，解决中文乱码
# encoding: utf-8nimport osnimport arcpynpath=r”D:\data”nfor filename in os.listdir(path):n filename=filename.decode(‘GBK’).encode(‘utf-8’)n print filenamen if(filename!=”polygon”):n …
ArcGIS python tool没有许可，按掩膜提取
narcpy.CheckOutExtension(“spatial”)#打开软件许可narcpy.gp.overwriteOutput=1#覆盖之前文件noutExtractByMask = ExtractByMask(tifPath,shpPath)#按掩膜提取nn
【273】利用ArcPy建立处理数据的脚本
这个脚本可以直接运行处理程序，首先在 ArcPy 上面测试，成功后写入文件，下面的代码实现将指定文件夹内部的栅格数据进行 Calculate Statistics 操作，否则在进行专题图制作的时候会出现问题。
ArcPy将地图文档 (.mxd) 文件转换为服务定义草稿 (.sddraft) 文件
在上一篇中讲到了ArcPy的使用和主要参考，所以这里就不再重复，有不明白环境和使用的请参考上一篇博客 ArcPy创建ArcGIS Server连接文件。
arcmap导出数据为shapfile文件与另存为mxd文件
arcmap导出数据为shapfile文件与另存为mxd文件
Python叠加两幅栅格图像
目的 n 现有两幅栅格图像，一个是某地区道路栅格图，一个是某地区土地利用类型图，需要将道路叠加到土地利用类型图中，即叠加后，重合的像元值以道路图为准，其余的像元值仍是土地利用类型图原有的像元值。 n 图1 道路信息图 n n 图2 土地利用类型图 n n 图3 结果图 n n具体实现 from gdalconst import *nfrom osge
ArcGIS Python编程案例-电子资料链接
ArcGIS Python编程案例(1)-Python语言基础rnhttps://www.jianshu.com/p/dd90816d019brnArcGIS Python编程案例(2)-使用ArcPy编写脚本rnhttps://www.jianshu.com/p/932e83db7f59rnArcGIS Python编程案例(3)-管理地图文档和图层rnhttps
pyshp创建shp点文件
# coding:utf-8nimport shapefilennnw = shapefile.Writer()nw.autoBalance = 1nw = shapefile.Writer(shapefile.POINT)nw.field(‘x’, ‘C’)nw.field(‘y’, ‘C’, ’40’)nw.field(‘v’, ‘C’, ’40’)nnn# with open(‘data’)
python gbk编码问题
https://www.cnblogs.com/GUIDAO/p/6679574.htmlrnrnrnhttp://blog.csdn.net/a491057947/article/details/47292923rnrnrnrnhttps://www.jianshu.com/p/b0e386361570rnrnrnrnhttp://blog.csdn.net/xiaopihaierletian/
跟我学python之mysql操作详解 二、增删改操作
上一篇我们讲了python操作mysql的查询操作，链接为 跟我学python之mysql操作详解 一、查询操作 这一篇我们讲python操作mysql的增删改操作：n 1.获取连接对象n conn = pymysql.connect(host=’127.0.0.1′, port=3306, user=’root’, passwd=’mytestcon’, db=’prod’, char
python使用arcgis的nc栅格图层代码
使用创建 NetCDF 栅格图层nhttp://desktop.arcgis.com/zh-cn/arcmap/10.3/tools/multidimension-toolbox/make-netcdf-raster-layer.htmnn
客户端两种获取要素属性的方法
http://www.classicning.com/blog/post/454.htmlrn目标：在浏览器中点击一个地图要素（绝大多数情况下仅限于多边形）立刻得到该要素的属性信息。在Web端获得要素的信息通常情况下有两种方法可以使用，一是WMS的GetFeatureInfo，另一个是WFS的GetFeature。rnrn事实上WMS不仅仅可以用GetMap请求获得地图，还可以用GetFeatur…
pycharm中matplotlib.pyplot ginput()中获得鼠标点击的位置失败
最近需要使用matplot库来进行鼠标响应，在图片上点击鼠标左键获取点的位置。nn先在网上找了些代码运行，发现点击鼠标后没反映，然后找了解决方法。nn我用的是Windows :nnSettings / Tools / Python Scientific / Show Plots in Toolwindow，取消勾。nnMAC系统下：nnPreferences /Tools / Python Sc…
GlobalId与普通的字段转换
// Converts a GlobalID field to a GUID field.npublic static void ConvertGlobalIdToGuid(IWorkspace workspace, Stringn datasetName)nn // Open the table.n IFeatureWorkspace featureWorkspace = (IFeatu
【285】ArcPy 暗色窗体设置
预览图rnrn设置如下rnrnDefault：rn
基于FME的水库功能分区生态红线批量处理方案（从文本到gdb）
作者：凡江林 QQ：564438737 项目背景2017年7月，环境保护部办公厅、发展改革委办公厅共同印发《生态保护红线划定指南》（环办生态[2017]48号）。（1）问题描述部门同事接领导任务，我部门需要提交全省各大中型水库的水功能分区红线gdb。原始数据为xls表格文件，初步整理后的格式如图（1）所示。图（1）最终需要的成果如图（2）所示。图（2）（2）流程思路 原始收集的…
arcpython设置处理边界extent
转载自：https://blog.csdn.net/liyanzhong/article/details/52095172
合并两个shape文件并删除冗余信息
合并两个shape文件，并且重叠部分只保留一个rnrn首先 feature to polygon 把线图层转化成面图层rn之后 analysis Tools-》 Overlay-》 Spatial Joint rnrnrn打开生成文件的属性表 rnrnrn1 表示有重叠 0 表示没有重叠。 在其中一个线图层右键 join 通过OID把新生成的图层和原始线图层关联rnrnrn打开关联后的线图层， 选
计算图层面积（针对于面要素）ArcObject c++
ndouble getAreas() //得到图层面积nnIFeatureClassPtr ipFeatureClass = getFeatureClass();nIFeatureCursorPtr ipFeatureCursor;nipFeatureClass->Search(nullptr, VARIANT_TRUE, &ipFeatureCursor);nIFeaturePt…
【01】在ArcGIS中使用Python脚本小细节
以下将平时在ArcGIS中使用Python脚本遇到的一些小细节记录下来，同时也参考了一些大神的资料。nn1、关于“/”和”\”的使用nnnG:/exercise/projects/ProPythonAnalysis/Map.gdbnnnG:\\exercise\\projects\\ProPythonAnalysis\\Map.gdbnn在使用MapDocument.findAndReplaceW…
Unity Projector 改进：去除边缘拉伸
Unity Projector可以模拟投影机效果，但是在投影画面边缘透明度不是0的情况下经常出现拉伸现象，对投影机的shader做如下更改可以去除掉拉伸部分。nnUnity自带的投影shader：nnn// Upgrade NOTE: replaced ‘_Projector’ with ‘unity_Projector’n// Upgrade NOTE: replaced ‘_Projector…
ArcGis对面层批量进行孔洞、环岛处理，使用Intersect、Dissolve、Erase工具
ArcMap 对面层批量进行孔洞、岛屿处理的试验。
c#+AE将要素类导出为shapefile文件（两种方法）
方法一：使用IFeatureDataConverter.ConvertFeatureClass方法nn private int ExportToShapefile(IFeatureClass pFeatureClass, IWorkspace outWorkspace)nn nn trynn nn …
在 vscode 上 使用 pycharm
参考文档：https://zhuanlan.zhihu.com/p/44110740n我没有通过win上面的cmd命令行去创建虚拟环境，而是在最新版本的pycharm中创建项目，最新版本的pycharm默认是将项目放在一个虚拟环境中去执行的。n以此可以查看python解释器的位置：nn然后将这个位置配置到vscode的配置里面：nnn…
Arcengine获取图层某个字段唯一值的方法
private ArrayList GetLayerUniqueFieldValueByDataStatistics(IFeatureLayer rnrnpFeatureLayer, string fieldName)rn rn ArrayList arrValues = new ArrayList();rn IQueryFilter p
关于arcpy中的Cursor，for循环问题
MyDbf=ur’@.dbf’nNavi=ur’@.shp’nnMyCurs=arcpy.UpdateCursor(MyDbf)#第一个指针n#NaviCurs=arcpy.SearchCursor(Navi)#如果将第二个指针放在这里，则第一遍到头后，指针就到底了，因此，要将指针放在第一个循环内部print( time.strftime(‘%H:%M:%S’,time.localtime()))
arcpy读取地图文档中的图层
n# encoding: utf-8n#获取地图文档中的图层列表n#mxd可引用网络驱动器上的文档nimport arcpynfrom arcpy import mappingnmxdpath=r”C:\Users\h\Desktop\fp.mxd”nmxd=arcpy.mapping.MapDocument(mxdpath)nlayers=mapping.ListLayers
（七）arcpy开发&&利用arcpy解析json生成shp数据
最近项目有这样的需求，采集的是在android平台上开发的平板app的上获取的图形数据。我们采集了点(point)，线(polyline)，面(polygon)数据，然后利用arcgis for android的api接口将这些图形数据转为json以string类型数据存储到splite数据库中。外业采集完后，再转为相应的shapefile数据,那么现在遇到的问题是如何将json格式的数据转为sh…
shutil.rmtree()
描述如果存在以下树结构 | —- a.txt | —-picturc | — folder1 —-| —–folder_1 —–| —— entertainmentnG：—-test —–
Python 之 ImportError: No module named ***
如果想使用非当前模块中的代码，需要使用Import，这个大家都知道。 rnrnrn如果你要使用的模块（py文件）和当前模块在同一目录，只要import相应的文件名就好，比如在a.py中使用b.py： rnrnimport b rnrnrn但是如果要import一个不同目录的文件(例如b.py)该怎么做呢？ rn首先需要使用sys.path.append方法将b.py所在目录加入到搜素目录中。然后进
RasterConvolver示例：栅格坡度计算
RasterConvolver转换器在FME 2018.1+中可用。nn栅格坡度计算nnFME已经有从栅格要素直接计算坡度的转换器；但是，RasterConvolver也可以用于创建相同的结果，栅格DEM卷积是一个很好的例子。nn卷积是一个在栅格要素和数字矩阵（核）之间进行的数学运算。操作和核的定义通常叫做滤波。nn这个练习是用多种方法对栅格要素进行卷积，通过合并结果得到解决方案的很好的示例。nn…
基于ArcGIS 的GPS数据地图匹配
GPS数据由于信号不良、通信异常、定位误差等原因，不能很好地落在道路上，所以需要对GPS数据进行地图匹配，与道路进行关联。而空间数据由于其特殊的数据格式，最佳的处理途径仍旧是ArcGIS Desktop系列软件。一、路网数据检查根据实验的需要提取出路网数据，对路网进行拓扑检查，包括不能有重叠、不能有悬挂点等，对拓扑错误进行修改；同时，需要进行“在折点处打断”操作，形成正确可用的路网数据。 n二、GP
WFS服务转要素类矢量的方法
1.找到可以用的WFS服务，这里以Esri官网提供的WFS为例进行说明讲解nnnhttp://sampleserver6.arcgisonline.com/arcgis/services/SampleWorldCities/MapServer/WFSServer?request=GetCapabilities&service=WFSnn2.在工具箱中找到WFS转要素类工具nnnnWFS服务…
java调用python 第三方库
1.java调用python总出现python第三方库没有定义nnhttps://blog.csdn.net/Johline/article/details/65633466nn2.Jython：java调用python文件之第三方包路径问题nnhttps://blog.csdn.net/ztf312/article/details/51338060…
在ArcGIS中创建Python工具或者pyt工具箱
ArcGIS Desktop 自带有大量的工具，可以执行各种各样的地理处理任务。通过 nPython 我们可以灵活地调用这些工具，把工具组织成自己的工作流，甚至创建一些新的工具。今天就说说在 ArcGIS 中创建 npython工具这个话题。nn在 ArcGIS 中我们可以通过两种方式来创建python工具：nn一种是在标准工具箱下创建 脚本工具:nnnn一种是创建 脚本工具箱
arcgis 栅格计算器，img叠加运算，con函数，img转txt
# -*-coding:utf-8 -*-nnimport arcpynfrom arcpy import envnimport osnimport shutilnnclass GridCalculator(object):n “””n :栅格计算器 con函数的使用n “””n def __init__(self):n #img文件临时存储路径n …
关于ArcGIS10.0中的栅格计算中的函数
版本升级确实很重要，在ArcGIS10.1中计算成功的，在10.0中出了问题。rn问题rn在进行栅格计算时，计算公式很简单，包括两个Ln函数：rn”-22.73 + 11.1 * Ln(5) + 0.34 * Ln(“人口.tif”)”rnrn可计算结果报错：000539 : Python 错误消息。Error running expression: rcexec()
ArcGIS教程：坡度
判断栅格表面的各像元中的坡度（梯度或 z 值的最大变化率）。
最小生成树之——已知一条边求最小生成树
题目：hdu4463rn题意：已知一条边求最小生成树rn题解：用kruskal的话，先将该条边连接。再求最小生成树。rn注意：是双向边！rnrn#includen#includen#includen#includen#includenusing namespace std;nint n;nint p,q;nbool vis[55][55];nstruct pointnn double x,
干货—ARCGIS拓扑规则说明详细讲解
ARCGIS拓扑规则详细说明rnrnrn1.拓扑规则简介；在实际应用时，有时需要在要素之间保持某种特定的关；相交（Intersect）：线和线交叉，并且只有；接触（Touch）：某线段的端点和自身或其他线段；悬结点（DangleNode，Dangle）：线；伪结点（PseudoNode）：两个结点相互接触；拓扑规则的种类可以按点、线、面（多边形）来分；2.点拓扑规则举例；点拓扑规则一：Mustbe
Python计算任意多边形面积算法
多边形面积求解的方法很多，其中比较多见的就是增加一个点P，然后分别连接多边形各个定点与P点，然后计算每个三角形的符号面积（面积有正负之分），求和就可以计算出面积。rnrnrnrn 鉴于上面的方法，我们也可以用另外一种更简单的方法，将P点直接放在V0点，那么我们要求的面积就是V0V1V2,VoV2V3,…,V0Vn-2Vn-3的面积和。此处注意，如果我们假设逆时针的三角形（如上图V0V
使用ArcPy实现Geometry与WKT的相互转换
原文地址：http://blog.csdn.net/linghe301/article/details/29554257nn nn在Web GIS迅猛发展的今天，使用浏览器来进行交互以其方便性、快捷性被广大用户所推崇，那么在传输格式方面，都已比较简单的JSON或者WKT来解决网络带宽带来的数据压力。nn在ArcGIS10.2版本，已经提供了Geometry与JSON的相互转换，而且在现实业务中，G…
mapgis裁剪（矢量裁剪栅格）
1.图像处理-图像分析-打开影像-添加矢量（区文件）若提示不能加矢量更改设置2.影像制图-区文件裁剪影像-点击裁剪区-保存即可
python中for循环的小技巧
偶尔看到一个神奇的用法:rnrn[[row[i] for row in matrix] for i in range(4)]rn这是干嘛的？就是把一个矩阵转制。。转制。。转制。。我擦那么吊？rnrn然后就调用了一部分试试，发现乐是这么个意思：rnrn当for语句再外面有个[ ]的时候。可以在for前面添加一个简短的语句，然后通过for形成一个list。好厉害的感觉！rn如：rn [row[2]+
Geostatistical Analyst 地统计模块插值
用ArcGIS工具箱中的插值工具进行空间插值，可选的种类比较少，而且经常会报错。现学习地统计模块插值方法。nn加载模块： customize- extensions- geostatistical analystn 选择模块已加载的模块geostatistical wizard，选择一种插值方法，然后调整参数，看最后的一些指标。Finish,得到一个临时文件结果，但是发现该文件并没有包括我们要剪切…
【284】◀▶ arcpy.da & arcpy 数据访问模块
rn使用游标访问数据rn数据访问模块 (arcpy.da) rnrn 参考：rnrnArcGIS Python编程案例(9)-ArcPy数据访问模块rn读取几何rn写入几何rn使用 Python 指定查询rnrnrnrnrnrn01rn rnda.SearchCursorrn只读查看表或要素类数据。
[GIS算法] 2.1.1 空间关系的判定-知识储备
nnnn关系运算nn【关系运算】检验两个集合对象的特定的拓扑空间关系的逻辑方法nn【怎么比较两个对象的拓扑空间关系呢？】基本方法： n1. 比较两个几何对象的内部、边界和外部的交集 n2. 根据交集矩阵产生的实体就可以知道空间关系nnnn普通拓扑学nn【普通拓扑学】 n1. 定义了内部、边界和外部 n2. 但是这些概念适用于二维空间中的二维对象，要使用与二维空间中的一维和零维对象时，就需要组合拓扑…
arcgis var 遍历 Polygon 属性
[code=”java”]rnvar array4Ring:Array = (geometry as Polygon).rings;rnrn数据结构为 图片内格式rnrnrnrn[rn [rn “m”:null,”spatialReference”:null,”y”:45.76632218427849,”z”:null,”x”:85.09957964091112,”exte
