利用shp点文件对影像按指定尺寸进行裁剪
nnnn ArcGIS在安装时已自动安装Python2.7以及arcpy包。在用编辑器编辑的时候需要将编译环境设置为python2.7。本文使用Pycharm进行编辑。设置为python2.7的操作如下…
nnnn ArcGIS在安装时已自动安装Python2.7以及arcpy包。在用编辑器编辑的时候需要将编译环境设置为python2.7。本文使用Pycharm进行编辑。设置为python2.7的操作如下…
// gdal_sample.cpp : Defines the entry point for the console application.//#include “stdafx.h”#include “gdal_version.h”#include “gdal_priv.h”#include “gdal.h”#include #include “ogrsf_f…
利用python脚本和ArcMap的字段计算机对属性表中的数据进行分类赋值。
(2017/2/20:这是我的第二篇博客,继续探寻上一篇未解决的利用arcpy.BuildPyramids_management()函数批量建立金字塔并加载遥感影像的问题, )
import osrnrn… target_folder = ‘D:\T20161202’rnrn… if not os.path.isdir(target_folder):rnrn… print(“The path is not a valid folder path: 0”.format(target_folder))rnrn… if not os.path.exi
n# encoding: utf-8n#获取地图文档中的图层列表n#mxd可引用网络驱动器上的文档nimport arcpynfrom arcpy import mappingnmxdpath=r”C:\Users\h\Desktop\fp.mxd”nmxd=arcpy.mapping.MapDocument(mxdpath)nmxd.activeDataFrame.zoomToSelectedF
一、给字段修改和添加值用到“arcpy.UpdateCursor”方法,加载shp图层shppath = r’G:\Desktop\Point.shp’ncursor = arcpy.UpdateCursor(shppath)二、将准备修改/添加的值放到list中listc = [1, 5, 8, 9, 13]三、循环游标,逐行更新新的字段值、此处以shp属性表的“test”字段为例,用到“upd…
策略:先创建一个GDB,将遍历到的第一个GDB里面所有的要素图层
复制带有要素数据集的数据库,在复制要素数据集的时候要更改工作空间,xiande
import arcpynfrom arcpy import envnfrom arcpy.sa import *nnarcpy.CheckOutExtension(“spatial”)nenv.overwriteOutput=1nnworkingDir=r”D:\\lyj\\code\\”nenv.workspace=workingDir+”globe_temp_raster”nnin_rast
import sys nreload(sys) nsys.setdefaultencoding(‘utf-8’) nnnimport arcpynfrom arcpy import envn n# Allow overwriting of output n# nenv.overwriteOutput = 1n n# Set current workspacen#ninPath = arcpy
看书的时候,遇见了如下数据:
今天试验了一下编辑的python脚本在
未来Python3会成为主流,但当前仍有部分项目停留在Python2阶段,如ArcGIS 10.2的ArcPy模块。本文以Windows下的Python3安装为例,介绍Python开发环境的搭建及PyCharm的安装。 这里首先安装Anaconda,这里选择安装Anaconda可以方便安装Python第三方库。…
https://blog.csdn.net/zglybl/article/details/752141831、在打开ArcMap工程文档的前提下,使用内嵌在ArcMap里的命令行窗口,将数据添加到当前的mxd文档中,参考代码如下:mxd = arcpy.mapping.MapDocument(“CURRENT”)df = arcpy.mapping.ListDataFrames(mxd, “Lay…
最近为了做点小事情,需要对12个Excel文件进行处理,每个Excel文件的数据量都非常庞大,所处理的数据是全国196个站点2000年至2011年每天的降水数据,需要根据这些数据
一、题目:二、代码(内有注释,易懂)# coding:utf-8nimport arcpynfrom arcpy.sa import *nimport mathnn# 烟囱高度(m)nH = 50n# 排放源强(mg/s)nQ = 150n# 风向,风速(m/s),正西风不用进行坐标轴变换nu = 2n# 点源坐标x0,y0nx0, y0 = 0.0, 0.0n# 存放预测点坐标的listnpoi…
通过sde文件链接,sde可在catalog中添加数据库连接,连接成功后,右键属性可以看到sde文件被创建了。在引用此文件即可nnimport arcpynarcpy.env.workspace=r”C:\Users\Administrator\AppData\Roaming\ESRI\Desktop10.1\ArcCatalog\连接到 localhost.sde”nfdlist=arcpy.L…
在使用ArcGIS日常的数据处理过程中,经常会用到一些批处理的功能,虽然ArcGIS工具箱本身提供了批处理命令,但是每一次都需要单独输入参数,这是比较烦人的。笔者以前经常使用AE编程的方式,调用GeoProcessor来进行输出处理,但是这种方式,对ArcGIS的版本有要求。所以综合来看,使用Python脚本的方式,编写一些批处理工具,更加的方便快捷。以下是在ArcGIS创建、使用自定义工具的基本…
准备工作nnPython安装n pip安装n 将pip所在的文件夹添加到环境变量path路径中n 下载相应的安装包,.whl格式。链接:https://www.lfd.uci.edu/~gohlke/pythonlibs/#wordcloudn pip3.4 install D:\Python34\Scripts\numpy-1.11.2-cp34-none-win32.whl#md5=532f1…
之前根据具体的应用情景写过一些关于Python模块的文章 nPython中的模块与包 nPython中第三方模块的安装 nPython中如何解决ImportError:No Module Named arcpy问题今天需要讨论一下Python模块内容如何编写的问题。模块Python将代码拆分为模块最初的目的是功能的独立,所以说模块提供的是某个独立的功能。
由于工作原因,需要在服务器中用python脚本相隔固定时间出密度图供网页调用,一直在摸索方法,现在将相应的python脚本粘贴如下:n# -*- coding: utf-8 -*-nimport arcpynmxd = arcpy.mapping.MapDocument(“D:\ServicePro\MD_Make.mxd”)narcpy.Extent(8000244.717300,327
参考书籍《空间分析》- 邓敏 刘启亮 吴静 编著 – 测绘出版社nn
利用GDAL/OGR提供的Python 开发包,可以快速、自主的进行矢量数据处理,环境搭建的方法,详见博客python地理数据处理环境搭建。下面对利用GDAL/OGR进行地理数据处理的知识点进行总结。1、OGR类结构图如下图所示,OGR的类结构图,与地理空间数据的组织方式基本一致,由数据源、图层、要素、集合对象和属性等组成,此外还包括驱动、属性定义及一些常量等。2、数据源与驱动数据源,顾名思义,即…
length函数rn(length ‘(“M20” “M24” “M30”))rn返回数组长度rnrn来源: http://www.afralisp.net/autolisp/tutorials/list-manipulation.php
环境:win10+geany+python2.7rnimport arcpy # ——根据ArcGis版本引用ArcGis10.2的python版本为2.7rnrnout_folder_path = ‘C:\Project’ #——路径任意设置(避免python识别的 \t \n等字符)rnrnout_name = ‘server.ags’ #——
今天遇到了一个问题,在使用a.Shape = b as IGeometry 时,经常报the geometry has no z values错误。。。查了一下,解决方法如下:rnrnrn首先,要判断a.Shape是不是含Z值的,如果把一个不含Z值的geometry赋给含Z值的a.Shape,是会报错的,反过来,也报一样的错。。。。。rn int index = pFeature.Fields.F
arcpy.mapping-认识arcpy.mapping by 李远祥 ArcMap提供了强大的地图制图功能,在实际的工作中,只要涉及到专题地图和地图册的制作,难免会遇到非常多的出图工作和地图图面元素的调整,这部分的工作虽然简单,但基本上是繁琐的重复操作,这些操作不单是耗时长,而且容易出错,往往需要打开ArcMap界面进行多次的操作。而对于一些数据量比较大的地图来说,通过ArcMap界面来逐
java/linux执行Python脚本工作原因,需要通过Python脚本进行接口访问获取接口数据。n分析网上java执行Python脚本无外乎三种方法: n 1、直接执行Python脚本代码 n 引用 org.python包PythonInterpreter interpreter = new PythonInterpreter(); n interpreter.
nn二、Pycharm导入arcpy站点包nn于是这里尝试着用python写一个工具。由于…
网络分析(Network Analysis)rn【背景】现实世界中,若干线状要素相互连接成网状结构,资源沿着这个线性网流动,这样就构成了一个网络rn【网络】就是指现实世界中,由链和结点组成的、带有环路,并伴随着一系列支配网络中流动之约束条件的线网图形,它的基础数据是点与线组成的网络数据rn【网络分析】rnrn网络分析的基础是:图论和运筹学rn通过研究网络的状态以及模拟和分析资源在网络上的流动和分配情况,对网络…
在arcgis的arcpy二次开发中,往往我们需要设置工具的数据类型。如下图的Add Geometry Attributes工具中Input Features选项。在这里我们需要设置的数据类型为矢量属性数据。nnnn又如下面的工具Points to line。
记录Python启动问题解决过程和方法。n解决问题:启动Python时总是报“ImportError: No module named site”这个错误。
nGISutil.createJPEG(mxdPath, outputPath, outTIFPath, issue, cycle)nn###################下面代码新建文件命名为 GISutil ###############n# -*- coding: utf-8 -*-nimport arcpy nimport os,sysnimport
按单个值进行重分类rn “重分类”工具可更改栅格中的各个值。按单个值进行重分类工具将以一对一的方式将一个值更改为另一个值。例如,执行鹿栖息地分析时,土地利用栅格中的各个值表示土地利用的不同类型,需要将其更改为一个优先级范围(如,1 到10),从而使各种土地利用类型对鹿而言均具有一定意义。将鹿相对偏爱的土地类型重分类为较高值,而那些不受偏爱的类型则重分类为较低值。例如,森林土地利用类型将重分类为
众所周知,使用arcgis toolbox的tools》》proximity》》near可以用来计算点与点之间的最近距离,但是通常是分为两个点图层。比如,一个是居民位置点图层,一个是超市点图层,则可通过上述方法超找到距离每个居民点最近的超市,并计算这个最短距离。rn但是如何查找同一图层内,距离某点最近的点及计算最短距离呢?比如,在居民位置点图层内,查找距离任意一个居民点最近的居民点(有点绕口)。其
按字段属性值分类导出数据,比如,左侧表中lev值共有 a b c d 四种,要求按照lev值分类输出 a.shp,b.shp,c.shp,d.shp。也就是把 lev =”a” 的所有的导出为 a.shp,其余类推。
我说重覆劳动事情就要arcpy帮我们做把。她说那就交给你。。。。。。。。。。。。rn 一下懵懂了,于是就拼命的找资料终于找到一份源码,像我这种懒人直接复制——粘贴,在执行时不报错就是不出结果。又懵懂了!于是只能一遍又一遍看源码,不懂查看
要求:出某县各个村的功能划分图,数据源有地籍子区,地籍区,功能分布图,影像。nn流程:用模型构建器搭建工作流导出一个村的一个要素类的数据库。建立一个mxd文档,加载源数据,设置好标题,图例,比例尺等。再用idle编写脚本。nn代码很简单,但是花费了很多时间调试。遇到了很多问题,比如python2 中文乱码问题,layer.replaceDataSource的第三个参数只能为名字,不能带文件名的后缀…
遥感,植被指数,NDVI
在arcgis的ArcTool工具箱中,提供创建栅格图层来提取单波段影像数据。rnrn选择原始多波段影像数据,和需要提取的波段,完成波段提取。
nimport requestsnimport jsonnnagisurl = “http://sampleserver1.arcgisonline.com/ArcGIS/rest/services/Demographics/ESRI_Census_USA/MapServer/1″npayload = ‘where’: ‘STATE_FIPS = \’48\’ and CNTY_FIPS = …
ArcGIS提供获取栅格数据最大最小值的接口,使用IRasterStatistics接口,其包含Maximum、Minimum属性代表最大最小值。nn由于最大最小值获取与波段相关,每个波段都有其最大最小值,因此最大最小值获取要使用各个波段接口获得。nn方法步骤:nn获取栅格图层IRasterLayern 从栅格图层获得栅格数据IRastern 将栅格数据转为栅格波段集IRasterBandColl…
参考:http://pro.arcgis.com/zh-cn/pro-app/arcpy/classes/spatialreference.htmnnhttp://desktop.arcgis.com/zh-cn/arcmap/10.3/tools/data-management-toolbox/define-projection.htmnn创建坐标系nnnn
1、前言 n最近项目需要通过软件出宗地图,但是软件输出的mxd未能达到成果要求,但是图的数量很大,人工修改会增加工作量,故想到通过arcpy批处理进行宗地图的整饰。 n2、整饰要素 n主要通过三个字段计算到需标注字段内,并进行符号编辑和标注。
n 这篇文章将主要讲述在进行脚本编程时用到的一部分ArcGIS模块和工具nnnnn目录nnnnn目录n1. arcpy和arcgisscripting站点包n2. ArcGIS中地理处理工具n2.1使用地理处理工具语法:n2.2 地理处理工具参数顺序nnn3.ArcGIS中地理处理工具箱n4. arcpy中的函数与类n4.1 arcpy中函数调用方法n4.2 arcpy中类的使用方法nnn5….
arcpy.mapping实战-专题图制图自动化 by 李远祥 最初对arcpy.mapping感兴趣是因为一次大规模的专题地图调整的需要,由于某某单位利用ArcEngine编写的专题图出图系统,出现了一些小问题,导致地图整饰元素的位置出现了误差。由于最终的生产任务已经到了尾声才发现一些问题,所以需要在出图系统之外使用ArcMap进行修改,这样针对上千幅成果文档进行修改,是个不小的工作量。后来
需求描述:现有一个矢量图层,其中ACRES字段记录了每一个面要素的面积,如何分段对该字段求和,如该字段数值在以下范围时0-50,50-300,300-1000,>4000,其面积分别是多少?nnnn解决方案:nn(1) n对于数据量少的,可以采用此方法;nnn可以通过属性查询(select bynattributes),通过表达式:ACRES >0 AND ACRES 将分段范围的要素类
这个例子是学习使用arcpy来统计面要素和线要素的个数,以及顶点(端点)个数的。利用到pointCount以及令牌方面的知识。首先还是对要素进行了遍历,然后统计相应的顶点数。下面我们来看一下实现代码。
本次我们来学习一下关于使用arcpy将hdf转tif数据。hdf格式数据还是常见的。那么什么事hdf格式的数据呢,来看一下相关资料。
今天我们要学习这个功能是利用arcpy在arcgis实现对某一个要素属性数据shapefile的某一字段进行数据的筛选。实现的过程使用了SelectLayerByAttribute_management函数将该字段中的某一个字段值分出来,然后将数据另存为shapefile文件。
n更改线坐标的教程nhttps://www.cnblogs.com/beautifulplanet/p/4317724.htmln#reusult为结果图层路径ncursor2 = arcpy.UpdateCursor(result)nn#inRow5[1][0] inRow5[1][1]为X,y坐标nn nnnvertex = arcpy.CreateObject(“Point”)nvertex…
表格有几行: var trCnt = table.rows.length; (table为Idrn)rnrn每行有几列:for (var i=0; i
欢迎使用Markdown编辑器写博客本Markdown编辑器使用StackEdit修改而来,用它写博客,将会带来全新的体验哦:nMarkdown和扩展Markdown简洁的语法n代码块高亮n图片链接和图片上传nLaTex数学公式nUML序列图和流程图n离线写博客n导入导出Markdown文件n丰富的快捷键n快捷键n加粗 Ctrl + B n斜体 Ctrl + I n引用 Ctrl
import arcpy.da as dann… import osnn… print(“os walk”)nn… for dirpath, dirnames, filenames in os.walk(“C:\Users\Administrator\Documents\ArcGIS\Default.gdb”):nn… for filename in filenames:n…
在学习arcpy脚本时准备使用pydev进行编写py脚本出现了错误首先是找不到 arcpy包 原因在于 默认检索目录中不存在arcpy的目录 可以查看sys.path进行查看。rn之后我使用了 program files(X86)目录中的arcpy 但是除了问题报错DLL load failed: %1 is not a valid Win32 application 问题就在于默认的python
nimport arcpynfrom arcpy import envnnsourceGDB=r”d:\polygon.gdb”nenv.workspace=sourceGDBntarget=”merge”nout_feature_class=”p2″narcpy.MakeFeatureLayer_management(out_feature_class, ‘
用python在图层里建立一个新的要素选择:对于字段类型是数字:arcpy.SelectLayerByAttribute_management(“pp”,”NEW_SELECTION”,'”ID” = 1′)对于字段类型是字符串:arcpy.SelectLayerByAttribute_management(“国道”,”NEW_SELECTION”,'”ROADCODE” = \’G107
将某一文件夹下的shp文件全部加载:rnimport arcpyrnmxd = arcpy.mapping.MapDocument(r”D:\tes\Operation.mxd”)rndf = arcpy.mapping.ListDataFrames(mxd, “Layers”)[0]rnrn# set workspace to directory of interestrnarcpy.env.workspa…
本次我们学习的是如利用arcgis提供的python接口arcpy在arcgis中创建gdb文件。在代码中,多处利用检查机制,首先检查路径传入的路径是否存在,如果传入的需要创建的gdb文件已经存在,那么代码会说明已经存在了。当路径已经创建好后,接下来要创建的gdb文件检查是否创建好了。接下来我们看一下实现的具体代码。