前言:(个人建议) 1.进行数据处理工作,由于通常数据量比较大.所以在个人电脑中,要将基础数据.不断增长的过程数据.结果数据等分门别类,使得简单易找. 2.将重要数据备份,因为在数据处理过程中,由于出错使重要数据损失,勿让重要数据成绝版.所以生成或保存压缩文件是很好的选择. 3.当数据量大,需要进行批量处理时.模板文档或代码要杜绝出现细节错误.容易一锅端,费时间且低效率. 处理工具:ArcMap 处理数据:MOD13Q1的NDVI成品(.tif)文件 研究区矢量文件(以河南省省界为例) 处理流程

有个处理数据的需求是把一个图层中的要素单独显示在底图上,设置固定的比例尺,并打印出图片. 考虑到后续会有重复的大量的数据要处理,决定使用arcpy处理. 首先新建一个mxd底图文档,把需要打印的地图都添加进来.在layout视图配置好出图样式,添加标题比例尺等挂件.标题中设置标识字段,在程序运行过程中会替换掉. 在layout视图下打开Python工具. 运行以下代码: # 第一步,布局打印格式,必须添加标题(插入文本,内容“title”),可以设置比例尺,标注等.必须切换到布局窗口,否则比例尺

工作中需要做大量图层的拼接,为了在拼接完成后还能知道原始数据文件是什么,所以写了个Python脚本对每个图层的SOURCE字段进行赋值. 附上Python代码: # -*- coding: utf-8 -*- # nightroad import sys import arcpy path = "C:/Users/nightroad/Desktop/SHP" arcpy.env.workspace = path ff = arcpy.ListFeatureClasses() for f

很久没用arcpy了,碰了好几次壁,把这次做的贴上来,以备下次可以跳过这些简单的问题 import arcpy arcpy.env.workspace = 'C:\Users\Qian\Documents\ArcGIS\Default.gdb' a=['sichuan1990','sichuan2000','sichuan2010'] b=['藏族','彝族','羌族','回族'] for x in a: for y in b: arcpy.MeanCenter_stats(x, x+y+"1.

在GIS数据处理中,选择数据是十分频繁的操作,常用的是"按属性选择"和"按位置选择",这两个功能虽然比较强大,但有时也不能满足实际需求.比如可能时常会遇到这样一种情景:将指定OID(假设3和6)的要素选择出来. 1.按属性SQL选择 最容易想到的是使用按属性选择构造WHERE子句("OBJECTID=3 OR OBJECTID=6")即可通过SQL选择出来. 2.属性连接 那么问题来了,如果给定的ID有100个,而且它们没有什么规律,构造SQL语

先看如何增加一个字段 函数:arcpy.AddField_management 语法:AddFields_management (in_table, field_description) 参数 说明 数据类型 in_table 要添加指定字段的输入表.该字段将被添加到现有输入表,并且不会创建新的输出表. 可将字段添加到 ArcSDE 的要素类.文件或个人地理数据库的要素类.coverage.shapefile.栅格目录.独立表.带属性表的栅格和/或图层. Mosaic Layer; Raster

目录: 1. 同一图层的多个要素合并(2种方法) 2. 导入带经纬度坐标的 txt 文件 3. 栅格计算器的实现 4. 添加 shp 文件(显示在 ArcMap) 5. 通过经纬度坐标生成 Polygon 6. 多边形重叠区域次数计算(赤潮发生频次) 7. 批量删除/添加 shp 文件字段 8. 批量修改文件的 Spatial Reference 9. 读取折线或面几何(点坐标) 10. 读取 Geometry(2种方法) 11. 批量修改图例显示内容 12. 批量引用图层样式(Symbolog

  出图是项目里常见的任务,有的项目甚至会要上百张图片,所以批量出土工具很有必要.arcpy.mapping就是ArcGIS里的出图模块,能快速完成一个出图工具. arcpy.mapping模块里常用的类有MapDocument.DataFrame.Layer.DataDrivenPages和TextElement. MapDocument类是地图文档(.mxd文件)对应的类.初始化参数是一个字符串,一般是.mxd文件的路径: mxd=arcpy.mapping.MapDocument(r"F:

ArcPy地理处理工具案例教程-批量添加栅格数据 商务合作,科技咨询,版权转让:向日葵,135-4855__4328,xiexiaokui#qq.com 关键字: Arcpy,python,地理处理工具 目的:批量添加任意工作空间(例如文件夹)里的栅格数据集,迭代处理子空间. 特点: 1. 支持任意深层的子空间,包括当前工作空间和任意深的子工作空间 2. 支持任意类型工作空间,包括文件夹,个人数据库(mdb),文件数据库(gdb) 3. 支持任意栅格类型,包括tif, jpg,grid, png

使用一个polygon矢量提取某个文件夹中所有的tif格式栅格数据 (要确保先安装好arcpy包) import arcpy arcpy.CheckOutExtension("spatial") arcpy.gp.overwriteOutput=1 arcpy.env.workspace = "E:\\data" #所有栅格影像所在文件夹 rasters = arcpy.ListRasters("*", "tif") mask

最大合成法(MVC)能够在Envi中的Band Math中进行,式子是B1>B2,可是无法批量化.本文实如今ArcGIS中利用Python代码批量进行,例如以下: 用到的Modis NDVI数据是在MRT中进行拼接与转投影后的月数据,一个月有两期,数据格式是.tif,文件名称的格式如:20040101.1_km_16_days_NDVI.tif.20040102.1_km_16_days_NDVI.tif代表2004年1月份的两期数据.本次处理的为2004-2013年的10年的数据. Pytho

曾经写过<使用Python脚本批量裁切栅格>,但今天又遇到这个情况则发现了问题.我们遇到的实际问题往往是有一个需要裁剪的影像(大块的),另外有一个矢量面,现在需要按矢量面每一个要素进行裁剪,无奈arcgis里的工具无法方便地做到.只能自己写工具,这次使用了clip而不是ExtractByMask,因为ExtractByMask有很多限制! 下面是工具的操作示例:按每一个要素进行裁剪栅格,输出栅格以选择的字段命名,前提是字段的每个值是唯一的. 其中,输出类型这个combox设置方法是: 下面是消

arcpy.mapping常用四大件-MapDocument by 李远祥 点开arcpy.mapping的帮助,可以看到其有限的几个类,看起来东西不是很多,但却是非常的使用.由于arcpy定位就是粗粒度的调用,因此,它的类和接口结对不可能像AO那么多,但可别小看了arcpy.mapping的这聊聊无几的类,它几乎封装了ArcMap的所有地图交互操作.下图是arcpy.mapping的类 其中包含了最为常用的几个大类,笔者认为最常用的主要入口的四类分别为MapDocument .Mapsurro

arcpy.mapping-认识arcpy.mapping by 李远祥 ArcMap提供了强大的地图制图功能,在实际的工作中,只要涉及到专题地图和地图册的制作,难免会遇到非常多的出图工作和地图图面元素的调整,这部分的工作虽然简单,但基本上是繁琐的重复操作,这些操作不单是耗时长,而且容易出错,往往需要打开ArcMap界面进行多次的操作.而对于一些数据量比较大的地图来说,通过ArcMap界面来逐个调整,很容易会因为调整过程中地图的刷新带来卡机或软件等待响应. 虽然可视化交互操作是大众基础,但面对巨

注意,ArcGis中shp属性表至少要保留除OID(FID).Shape之外的1个字段. 例:批量删除“省级行政区域”除“FID”.“Shape”.“Area”之外的字段. import arcpy #如果要素类已加载至工作空间,可以直接使用其名字(此处是"省级行政区域"),否,则使用.shp文件的全路径. #code source: https://www.cnblogs.com/yzhyingcool/ QQ:975601416 fieldList=arcpy.ListFields

这一段是批量定义要素类(FeatureClasses)投影的ArcPy代码: 把要处理的要素类塞进一个文件夹(工作空间,workspace),然后将代码开头的路径换成这个“文件夹”的路径,处理完后再做其他处理(比如批量塞回原来文件夹,vba,python啥的代码.小工具处理即可). 或者,把要处理的要素类都加载到ArcMap工作空间,然后把下面代码set environment workspace(设置工作空间)的这一行代码干掉,代码执行的默认工作空间就是当前ArcMap的工作空间了. defi

在实际操作中,经常对数据库文件进行合并.裁切等.如果遇到gdb比较多,要素层比较多,而且还存在数据集.虽然ArcGIS中的批量处理的功能,但填写参数过程也比较麻烦,如果一次性处理过多,程序容易停止工作. 所以,选用Python来处理这样的问题,优势就非常明显了.把路径设置好,读取数据库及要素类,然后进行合并操作.下面是Python示例代码: # ---------------------------------------------------------------------------

对栅格的裁切,我们通常使用裁切(数据管理-栅格-栅格处理)或按掩膜提取(空间分析-提取分析)来裁切,裁切的矢量要素通常是一个要素图层或Shape文件.如果要进行批量处理,可以使用ToolBox中的批量处理工具.但是,有时我们选择的裁切矢量要素不是一个要素图层,而是要素图层中的一个要素,如按标准图幅裁切栅格影像,以前的做法是将要素一个一个的导出来,再进行批量处理.快速导出方法:ArcGIS按字段属性分割文件. 现在,可以使用Python脚本来批量裁切,在网上找了位大神的原始代码学习一下: #-*-

 Howto: 在ArcGIS10中将地图文档(mxd文档)批量保存到之前版本 文章编号 : 38783 软件: ArcGIS - ArcEditor 10 ArcGIS - ArcInfo 10 ArcGIS - ArcView 10 操作系统: N/A 摘要: 以下说明展示了如何在ArcGIS10中通过简单的脚本将地图文档(mxd文档)保存成之前版本. 内容: 在ArcGIS 10中,地图文档被保存为10的版本.为了与使用早期版本ArcGIS的用户共享地图文档,可以使用File->Save

友好阅读链接:(http://swj.me/2015/08/26/batchPublishtools/) update0918: 修复了创建链接文件时而出错的bug 修复了在24011的警告已处理的情况下,依然出现提示的异常. 将多个脚本文件,压缩成一个脚本文件.压缩后的脚本文件名称为Publishservice.py,使用如下命令执行工具 python Publishservice.py 重新创建了一个repo,更改了脚本git 存放的地址,新地址如下: BathPublishMapServi