yGIS

在文章开始之前，大家可以参考学习视频。https://edu.csdn.net/course/detail/25535首先，大家需要安装好pycharm,使用社区版的pycharm就足够了。在编写arcpy脚本之前我们需要将arcpy环境设置好。下面我们来设置一些arcpy需要的环境。 从创建工程开始路径设置 工程导入arcpy站点包，首先我们使用pycharm新创建一个工程...

在安装arcgis10.5之前，可以微信扫一下图的二维码，或是搜索公众号“测绘科技”，回复：arcgis10.5，获取下载链接 arcgis license部分 安装license 选择arcgis license的安装包setu...

现在我们来学习一下关于arcpy.da.Walk这个api，首先呢，该函数需要提供三个参数，分别是工作空间、数据类型、类型。现在我们测试一下下面这段代码，注意使用的过滤条件。设置数据类型、类型都是any。#coding=utf-8import arcpyimport osws = r"D:\Data\中国国界和省界的SHP格式数据\省界\demo.gdb"walk = arcpy....

今天在学习关于工作空间相关知识，用于了下面的函数来实现判断当前输入的要素参数在哪个工作空间。其中主要运用到了arcpy.Describe函数。def get_workspace(featureClass): catalogPath = os.path.dirname(featureClass.catalogPath) if arcpy.Describe(catalogPath...

在arcgis中字段修改工具，其中的位置如下图所示。我们可以来查看函数对应的参数。arcpy.AlterField_management(in_table="test",field="YY",new_field_name="zz",new_field_alias="#")而我们现在需要使用arcpy的方式来修改我们的创建好的字段，其中我们可以使用提供的api中的alter函数，现...

这个例子是学习使用arcpy来统计面要素和线要素的个数，以及顶点（端点）个数的。利用到pointCount以及令牌方面的知识。首先还是对要素进行了遍历，然后统计相应的顶点数。下面我们来看一下实现代码。#coding=utf-8import arcpydef getVertCount(featureclass): desc = arcpy.Describe(featurecla...

首先在测试这段代码之前，我们这里新建一个字段用于测试。然后将FID的值相等过来。最后我们的属性表如下图所示。那么我们来看一下代码的实现过程，首先我们读取了一个要素类，然后设置我们需要修改长度的对应字段，然后设置对应的长度。在没有将字段值删除之前，我们需要将原有的数据给保存起来。然后再将我们的shapefile数据需要删除的字段给删除点，之后就是重新创建字段了。创建完字段后，...

今天我们来学习一下关于获取栅格数据的一些信息，首先我在这里呢是遍历某个工作空间下面所有数据，然后读取这些栅格信息，最后将这些信息写入到txt文件中。其中有关arcpy获取栅格信息的标记可以参考esri的相关文档。http://resources.arcgis.com/zh-cn/help/main/10.1/index.html#//0017000000m7000000，下面我们来看一下具体的实现...

今天我们来学习一个关于使用shapefile文件导出csv文件的功能。首先，这里需要遍历shapefile文件。然后将字段值的写入到csv文件中。具体的实现代码如下所示。#coding=utf-8import arcpyimport codecsimport sysdef getCount( fc): r = arcpy.GetCount_management(fc)...

首先小伙伴们可以参考arcpy开发学习视频。https://edu.csdn.net/course/detail/25535本次学习，我们在arcgis中利用arcpy来创建渔网工具，其中利用到CreateFishnet_management函数，该函数有多个参数。可以参考一下arcgis的帮助文档。http://resources.arcgis.com/zh-cn/help/main...

在arcpy开发中，我们会使用到关于信息提示，现在为大家提供一个关于这么方面的信息类。首先我们来看一下相关的函数。获取时间，其中格式为%d%b_%H%M，来看一下函数定义。 def getTS(self): # Gets a timestamp in a default format (day month time) startTime = date...

本次我们学习的是关于如何利用arcpy实现，给定一个字段名，和一些给定关于该字段名一些具体值。然后统计关于该字段名中不是在给定具体值中其他列值。具体来说就是，我们现在有一个字段名GBCODE，然后里面有很多值，你比如说61030,61031,61032、61033，61034等等，而现在呢，我们需要统计除了61030，61031，61032其他值的数目。那么如和实现呢，首先我们需要验证一些输入的字...

今天我们要学习这个功能是利用arcpy在arcgis实现对某一个要素属性数据shapefile的某一字段进行数据的筛选。实现的过程使用了SelectLayerByAttribute_management函数将该字段中的某一个字段值分出来，然后将数据另存为shapefile文件。我们来看一下这里的测试功能，我们对如下图的数据操作，从中可以看出属性记录一共有800多条，现在对FID字段进行分类选择，那...

这几天，有小伙伴有这么样的需求。首先为了便于说明，将数据展示出来，能够直观了解要描述的问题。来看一下下面的数据。他有很多条记录数据，其中字段INPUT_FID按照一大类进行标记，比如0为一类，1为一类，依次类推。然后现在要求出，每一类中对应字段DISTANCE的最小值的记录，不知道我的描述够清楚了没。就好像，我们现在有很多个班级，在一次考试中，我们要找出每一班中，语文成绩最低的那位同学。那...

arcgis开发python视频教程https://edu.csdn.net/course/detail/25535最近在群里看到了小伙伴们在做三调，有需要锐角检测这样的工具。于是试着写了一个，这是利用arcpy写的一个锐角检测工具，最后做成界面如下图所示。那么代码是如何实现呢。首先是对每个面进行遍历，分别找出面中的顶点。对顶点进行编号，然后就是利用平面向量乘积等数学相关内容，来计...

UpdateCursor游标用于对记录的更新，当然也可以用于数据的删除。我们只要正确的读取好数据后，获取游标，然后在设置我们需要取出的字段，当然也可以设置过滤条件。最后使用updateRow函数来更新即可。现在我们来做一个测试案例。如下图所示，是我们准备好的一份测试数据。现在我们设置要获取字段TDYT，也可同时获取其他的数据。我们在这里做这样的操作，如果TDYT这个字段里面的数据值为空，那...

今天我们来学习一下关于字段映射方面的知识，来看一下esri关于字段映射的介绍。FieldMappings对象是一组FieldMap对象，它用作执行字段映射的工具的参数值，如Merge。要处理这些对象，最简单的方法就是先创建FieldMappings对象，然后通过添加要组合的输入要素类或表对FieldMap对象进行初始化。提供了所有输入后，FieldMappings对象将为所有...

今天看到群里有小伙伴，有这么一份数据，他需要对数据逐一，即按每一条数据导出来。这个属性表中，有多少条记录，就能导出多少条shapefile数据。其实这样的功能在之前的博客中就有提到过了，那么其实在ArcGIS Desktop中是有这么的一个工具来实现的，在你的分析工具中有分割【split】可用来实现这样的功能。在开始之前我想着要使用代码来实现，是这样的思路：先遍历属性数据，然后创建属性数据，再导出...

大家是否记得之前，写过一篇关于使用python脚本对数据进行分类的例子，如果不清楚的话，可以查阅下面的链接。https://blog.csdn.net/u010608964/article/details/89929059在上面的例子中，需要对数据进行读取，然后利用更新游标对属性表遍历，按照我们的逻辑对数据进行分类，然后再更新分类字段。而今天我们使用的是ArcGIS Desktop的字段计...

今天我们要学习的是，使用我们的数据对某个字段进行分类，有点类似分级渲染的感觉。如下图所示的数据。现在对字段BOU2_4W_ID依次以等于1，为A类，在1至2000为B类，2000至4000为C类，4000至6000为D类，大于6000为E类。在这个过程中，我们会用到字段的创建，以及遍历，和更新。好了，下面是我们的实现代码。#encoding: utf-8import arcpy...

今天我们来学习一下，使用arcpy实现多幅影像批量拼接。首先在这里，大家需要知道使用python的for循环，与数组的遍历。在python有这样的遍历的方式。方式一:whu_gis = ["大地测量","摄影测量","地理信息系统"]for itemGis in whu_gis: print itemGis#****运行结果*****#大地测量#摄影测量#地理信息系统...

今天我们来学习一下，使用arcpy来寻找polyline要素类的起始端点和结束端点。首先，程序开始会对要素类polyline进行遍历，然后利用半正矢公式对两个端点进行计算。这里主要利用到半正矢公式，以及在arcpy中创建要素类字段，并添加相应的值。我们来看一下具体的实现代码。import arcpyimport sysimport mathdef haversine(point1, ...

今天我们来学习wicket的事件，最后做出的页面如下。我们来看homepage.javapublic class HomePage extends WebPage { public HomePage() { final WebMarkupContainer containerInTheMiddle = new WebMarkupContainer("conta...

这次我们来学习一下，关于为polygon面计算几何中心,我们指定这个中心为点对象，然后将polygon中所具有的属性一起赋值给这些点对象，最后创建相应的shapefile数据。在这个例子中，将会用到字段的创建，要素创建，游标的创建。首先我们做一份数据，如下图所示。由几个polygon组成，相应的我们为每个polygon面要素添加属性字段，以及一些对应的值。来看一下编码，首先是读取pol...

这次我们又来学习一下arcgis的python脚本开发。今天我们要学习的这个例子是利用arcpy读取某一工作空间下面的所有shapefile数据，然后将读到的矢量shapefile数据一些信息，写入到txt文件中。在这里使用了，arcpy.ListFeatureClasses函数来获取目录下所有的要素类。arcpy.Describe获取文件的描述信息，这里为什么不指定shapefile，是因为...

今天我们来学习这个例子是创建一个gdb、然后在该gdb目录下创建一个shapefile文件，再为这个矢量文件定义投影，相应的添加到字段。这过程中使用的函数有arcpy.CreateFileGDB_management、arcpy.CreateFeatureclass_management、arcpy.Project_management、arcpy.AddField_management。最后创建...

本次我们学习的是如利用arcgis提供的python接口arcpy在arcgis中创建gdb文件。在代码中，多处利用检查机制，首先检查路径传入的路径是否存在，如果传入的需要创建的gdb文件已经存在，那么代码会说明已经存在了。当路径已经创建好后，接下来要创建的gdb文件检查是否创建好了。接下来我们看一下实现的具体代码。#coding=utf-8import arcpyimport osi...

本次我们来学习一下在arcpy中如何将.csv格式的数据转为shapefile。其中csv为点数据，这些点数据都是如下图所示。现在利用arcpy首先要定义要投影，然后创建好字段。在遍历.csv文件。将数据填入到shapefile文件中。具体的实现代码如下图所示。import arcpy, systry: path = "C:\\Users\\qin\\Desktop\\...

在发布地图服务的时候，我们需要针对我们的数据制作一个切片缓存（方案）。在arcgis中切片缓存(方案)格式为.xml。下面我们来看一下在arcgis中是如何制作切片缓存(方案)。首先需要将我们的数据另存为一个.mxd，然后在arcgis的工具箱如下图所示，依次选择弹出如下图的设置框。我们这里一一对其进行设置。首先选择mxd文件的路径，然后设置Number of Scales，我们这里设...

首先这是一个利用arcpy求交、以及数据融合、文件复制操作例子。大概需求是这样的，这里有一个图幅结合表的shapefile面数据，该面属性表数据对应有图幅号，现在利用某一区域（暂且一个县）的面shapefile数据与图幅结合表相交，计算出相交结果的图幅结合表，然后遍历该结果图幅结合表的图幅字段，最后利用利用这些字段信息在指定文件将多个分幅好的DEM数据拼接为该区域（一个县）一大块DEM数据，整个工...

首先小伙伴们，可以看一下我们的学习视频https://edu.csdn.net/course/detail/25535明天就国庆了，刚好最近在使用python在arcgis方面批处理开发工作，写了一些总结，以备后用。这次我们先做了一个测试工具，其主要目的是利用toolbox工具箱的Mosaic To New Raster将多个栅格数据合成一个栅格数据。如果说栅格数量少的话，手动选择并没...

今天遇到了pyt后缀的pyt脚本声明加载到Pycharm是灰色的情况，具体如下图所示。那么如何解决呢。首先在pycharm中找到setting。如下图所示，以次选择【File Types】、在Recongnize File Type选项中选择Python,然后在Registered Patterns中追加.pyt后缀文件。操作界面如下图所示。重新设置后点击【appl...

值中秋之际，和小伙伴出去后，借点时间看一下关于多实体实时拖拽效果（jig）,在AutoCAD中，通过鼠标的移动，动态的展示效果会给绘图者直观的感觉。在AutoCAD中关于拖拽有两个实现类，第一是EntityJig,另外一个是DrawJig类。其中前者只能针对一个实体，而另外一个是针对多个实体，可以实现拖拽效果。而对于前者者继承类中需要实现Sampler和Updata两个函数，Sampler用于数据...

这次来学习一下关于点是否在多段线内部，首先要说明一下的是，该多段线为闭合的多段线。下面使用两种方法：一、使用多重多边形。使用函数AppendLoopFromBoundary(pline, true, tolerance)传入多段线范围线和容差创建多重多边形，然后根据函数IsPointInsideMPolygon((point, tolerance). Count返回值是否为1来判断是否在多段线...

实体是在idata软件中展现图形的载体，通过设置实体的相关属性，让各种实体表现出各自的差异。可以说实体在idata中算是比较复杂对象。下面来学习一下idata的实体相关函数。. 函数名 说明 addToScene() 将实体加入到容器中 boundingRect() 获得实体外包...

类似于AutoCAD的命令行中选择提示，让用户选择，或者输入一定信息的用户交互，idata也有用户交互。查阅了idata相关的交互函数，和AutoCAD相比，要稍微简单得多，函数的数量也少了很多。下面我们来看一些函数，以及了解函数的具体说明。如下表格所示。 函数 说明 iDataGetInt 获取用户的整数数据输入 ...

本篇文章是利用teigha这个库来执行cad的dwg文件相关操作。利用该库读取了dwg数据库，然后对数据库复制块，设置块的单位为无。开发这样功能工具的缘由，来源于dwg文件是由arcgis转换过来的，其中带有投影信息，在进行坐标转换的时候出现由于投影信息的缘故，使文件数据库dwg文件的单位为英寸，导致了坐标转换出现了很大的偏差。而我们有相当数量这样的文件，为此急需开发这样的小工具。而如果是使用...

南方数码的idata数据工厂目前有针对第三次全国土地调查的版本。下面我们来学习一下有关图层的一个小案例。这里案例是通过指定图层名称，然后选择图层。如果数据库中存在指定名称的图层，然后遍历图层的所有实体，让实体高亮显示。最后计算该图层所有的实体个数。写这个代码，大家需要在Visual Studio中添加CSInterface.dll的引用。整个代码如下。using Southgis.iData...

这是一个与学习创建和获取扩展字典对象与扩展记录数据相关的内容。我们知道在南方cass中有一部分数据是存在扩展数据xdata属性中，其中宗地的宗地基本属性就是印证了这一点。而我们接下了的学习是字典Dictionary扩展数据、扩展记录数据Xrecord。这对了解AutoCAD扩展数据具有一定帮助。工程的代码是摘自相关网站，源代码使用的是..net6.0，由于本次使用的是AutoCAD2010，因...

作为和地理信息相关的工作者，总是免不了和制图相关，这次是和arcgis制图相关，使用arcgis制作属性矢量面，即polygon。下面是具体的步骤。首先如下图所示，在arcgis catalog连接的任意文件夹右击，依次选择【New】、【Shapefile】。进入如下图的界面，该界面让用户选择属性数据类型、以及设置该属性数据的投影信息。如下图所示、在1部分是地理坐标系、而2是投...