arcpy专题

GIS之arcgis系列09：arcpy实现克里金差值

矢量点数据经过克里金差值后可以转换成栅格数据，那么就需要了解一下什么是克里金差值。什么是克里金法? 　　IDW(反距离加权法)和样条函数法插值工具被称为确定性插值方法，因为这些方法直接基于周围的测量值或确定生成表面的平滑度的指定数学公式。第二类插值方法由地统计方法(如克里金法)组成，该方法基于包含自相关(即，测量点之间的统计关系)的统计模型。因此，地统计方法不仅具有产生预测表面的功能，而且能

GIS之arcgis系列10：arcpy实现批量掩膜提取

按掩膜提取 (Spatial Analyst) 提取掩膜所定义区域内的相应栅格像元。 OutRas = ExtractByMask(InRas1, InMsk1, "INSIDE") 使用情况输入栅格中的其他属性（若有的话）将按照原样添加到输出栅格属性表。根据所记录的属性，某些属性值可能需要重新计算。将多波段栅格指定为输入栅格（Python 中的 in_raster）值时

GIS之arcgis系列08：arcpy实现批量excel转矢量点

文件夹内有很多excel表，每个excel表内有很多数据，爬取每条数据的日期、精度、纬度，根据经纬度将数据转换成矢量点，并指定坐标系。完整版代码如下：代码均已经过测试，可直接copy至arcgis工具箱中使用。（copy后注意改成适用于自己数据的代码，微调相关信息） # -*- coding: utf-8 -*-import pandas as pdimport arcpy

GIS之arcgis系列07：conda环境下安装arcpy环境

首先将python27环境下的“Desktop10.8.pth”拷贝到anaconda环境下。路径如下（仅参考）： C:\Python27\ArcGIS10.8\Lib\site-packages\Desktop10.8.pth D:\Anaconda\Lib\site-packages 在anaconda prompt中穿创建一个新环境 conda create -n arc

网络分析(ArcPy)

一.前言 GIS中的网络分析最重要的便是纠正拓扑关系，建立矫正好的网络数据集，再进行网络分析，一般大家都是鼠标在arcgis上点点点，今天说一下Arcpy来解决的方案，对python的要求并不高,具体api参数查询arcgis帮助文档即可。二.数据资源在我的资源发布里，下载即可三.步骤新建数据库和数据集，并将数据导入数据集中建立拓扑—导入要素进拓扑—对拓扑添加规

arcpy批量导出图且图名为shp属性值

1.打开arcmap加载需要导出的图。需求是逐村显示“村界内图斑”并导出为图，在导出每个村时不显示周围的村和“村界内图斑” 2.arcmap上方空白处右键打开“数据驱动页面” 3.在“数据驱动页面”工具条点击第一个图标，打开“设置数据驱动页面” 4.在“设置数据驱动页面” 勾选“启用数据驱动页面”，在“索引图层”中选择需要驱动的图层（按照哪个图层的来批量出图），名称字段（按照

Arcpy安装和环境配置

一、前言 ArcPy 是一个以成功的arcgisscripting 模块为基础并继承了arcgisscripting 功能进而构建而成的站点包。目的是为以实用高效的方式通过 Python 执行地理数据分析、数据转换、数据管理和地图自动化创建基础。该包提供了丰富纯正的 Python 体验，具有代码自动完成功能（输入关键字和点即可获得该关键字所支持的属性和方法的弹出列表；从中选择一个属性或方法即可将

在Pycharm中使用ARCGIS Pro提供的python.exe或者ArcPy

在安装了ArcGIS Pro后，如果需要在Pycharm中使用Pro提供的python.exe或者ArcPy，还需要进行相关设置，主要是3步：1、克隆python环境；2、Pycharm选择python；3、导入ArcPy等。详细内容可见：公众号GISideas 1、克隆python环境 ArcGIS Pro自带的Python版本是无法被修改的，因此需要先使用Pro自带的包管理器进行

ArcPy高级开发教程—要素操作

今天，介绍使用ArcPy对矢量数据的要素feature进行操作方法。目录 1. 空间分析相关建立缓冲区计算要素面积和周长按属性选择 2. 要素遍历 3. 导出每个要素 4. 生成shpfile文件 1. 空间分析相关建立缓冲区 # _*_ coding: UTF-8 _*_import arcpy, osdef buildBufffer(shpPath,

ArcPy高级开发教程—属性表字段操作

今天，将介绍使用ArcPy对矢量数据的属性表的字段进行增删改查操作。 1. 属性表字段查询 import arcpyfields = arcpy.ListFields(shpPath) shpPath是shpfile的路径，fields得到的是包含Field对象的列表，Field对象具有字段名称、别名、长度、精度等属性信息，例如，你只想得到属性表中包含的字段名称，那么可以直接获取Field对

arcpy多线程热力图

起因是这样一段对话，领导：你会用脚本生成热力图图片吗？我：可以研究下。领导：那这个需求就给你了。我：… 经过一番研究，研究出大概的思路，先将有经纬度的表中的数据筛选出表并生成xy事件，接着利用核密度工具生成栅格，最后呢裁剪栅格通过mapping包出图。出图过程比较慢，所以考虑使用多线程(一张图五分钟，十九张图多线程九分钟你说吼不吼)。多线程第一次用，一开始打算用数据库存放

arcpy_根据excel批量修改字段别名

主要实现从Excel表中读取gdb表需要修改的字段原始名称，以及需要修改为的字段别名，下面是Excel表的结构。下面是代码 #coding:UTF-8#python版本：2.7#作者：K#使用前需要安装pandas库；#pandas版本要与python27配套，否则无法使用；#使用时务必保证被修改的图层处于被打开的状态，建议关闭arcgis;#参照ModifyFieldNa

用arcpy创建图层、定义坐标系、添加字段

针对单一图层来创建的，前提是已经存在gdb文件。对于批量创建，或者根据excel、txt创建，肯定是有方法，只是还不会。学习中。 # -*- coding: utf-8 -*-# 创建图层，定义坐标系，并添加字段# Python version:2.7# Arcmap version:10.2.2# Author: ****# 参照arcgis原生帮助文件：https://resou

【Arcpy学习实践教程】wgs84坐标系和火星坐标系的转换中demo的对与错

度娘和谷哥已经变成了我们学习工作生活中必不可少的工具。更有甚者，甚至已经不用输入法来搜索，而直接通过语音识别来搜索。但是我们搜索的结果真的可靠？我们在找到我们想要的资源之后是否有认真检验一下，我们找到的代码，找到的资料是否正确。当没有思考力和鉴别力的搬运工进入大众视野时，我们就需要谨慎起来了。最近，因为工作的原因需要对高德坐标（即火星坐标）和wgs84坐标系实现互转。一直记得有位大神曾经在

arcpy实现kml批量转出为shp 包括shp合并

# -*- coding: utf-8 -*-import syssys.path.append("D:\\Python27\\ArcGIS10.4\\python")sys.path.append("D:\\Program Files (x86)\\ArcGIS\\Desktop10.4\\arcpy\\arcpy")import arcpy, os# Set workspace (wh

【ArcPy】游标访问几何数据

访问质心坐标相关数据结果展示代码 import arcpyshppath =r"C:\Users\admin\Desktop\excelfile\a2.shp"with arcpy.da.SearchCursor(shppath, ["SHAPE@","SHAPE@XY","SHAPE@TRUECENTROID","SHAPE@X","SHAPE@Y","SHAPE@Z","SH

【ArcPy】游标访问数据

游标类型类型方法说明搜索游标arcpy.da.SearchCursor检索行更新游标arcpy.da.UpdateCursor更新和删除行插入游标arcpy.da.InsertCursor插入行使用搜索游标遍历所有结果展示代码 import arcpyshppath =r"C:\Users\admin\Desktop\excelfile\1.shp"wit

【ArcPy】简化ArcGISPro默认Python环境体量

参考文献安装 ArcPy—ArcGIS Pro | 文档

Python中ArcPy实现对不同时相的栅格遥感影像依据其成像时间分别批量拼接

本文介绍基于Python中Arcpy模块，对大量不同时相的栅格遥感影像按照其成像时间依次执行批量拼接的方法。在前期的文章Python中arcpy栅格创建与遥感影像多景数据批量拼接Mosaic（https://blog.csdn.net/zhebushibiaoshifu/article/details/118901707）中，我们已经介绍了利用Python实现栅格遥感影像批量拼接的方法

[ArcPy] 批量波段合成 Landsat8数据为例

# -*- encoding:utf-8 -*-import osimport sysimport arcpyfrom arcpy import envworkspace = r'F:\landsat8' #文件outpath = r'D:\layer_stacking' #结果存放# function:按顺序拼接字符串def getInputParam(files):ret = ""

[ArcPy] 遥感影像去黑边-第六届全国大学生GIS技能大赛试题

完整试题：http://www.higis.cn/753 文章目录问题分析ArcPy暴力剔除边缘的像元值代码将Py脚本制定成自定义工具结果找出边界后掩膜问题去除遥感影像黑边数据：分析【方法一】将边界进行矢量化，然后用掩膜提取缺点：遥感影像边缘是锯齿形，矢量化边界精度不高【方法二】边界像元值为0，将0值重分类为NoData 缺点：将内部像元值为0

[ArcPy] 将道路(shp)平均分配给志愿者生成6*6格网 | 第六届全国大学生GIS技能大赛试题

第六届全国大学生GIS应用技能大赛下午题目完整试题：http://www.higis.cn/753 文章目录前言题目数据说明分析结果 ArcPy实现代码制作ArcGIS自定义工具生成的6*6格网连接志愿者表所有文件下载： http://download.csdn.net/download/summer_dew/10164828 前言题目创建一个6*6

基于arcpy api的RPC校正

本文是来填坑的，之前写的RPC的基本原理，从底层写RPC校正还是很有难度的，目前为了快速的工程化，先写了一个基于arcpy中RPC校正的脚本。核心代码如下： import os import arcpydef OrthoCorrected(in_ras, out_ras_path):'''Incorporates elevation data and image metadata to ac

arcpy高德爬取路况信息数据json转shp

最近工作上遇到爬取的高德路况信息数据需要在地图上展示出来，由于json数据不具备直接可视化的能力，又联想到前两个月学习了一点点arcpy的知识，就花了一些时间去写了个代码，毕竟手动处理要了老命了。 1、json文件解读 json文件显示如下：通过api和json组织结构解读，我们只需要道路的polyline坐标点空间信息，道路名称，道路方向、道路限速和道路拥堵程度等属性信息。所以我们需

Arcpy基础入门-6、arcpy要素属性操作

空间数据具有三大特征：空间特征、属性特征和时间特征，同样组成一个几何要素的对象也包括空间几何、属性以及时间表达. 这一节我们以一个实例来说一下arcpy对几何对象属性的操作，属性的操作和数据库的操作基本相同，包括查询、添加、更改和删除，添加和删除主要是针对整个属性添加或者删除一个字段，查询和更改是针对一个要素的属性进行操作，我们来先看一个实例： import time,osimport

Arcpy基础入门-5、读取自定义格式xml

arcpy作为一种脚本开发语言，其功能的目的是为用户提供强有力的可扩展性。这一节通过示例，介绍比较下使用python和arcpy读取自定义格式的xml并构建要素的方法。先简单描述下自定义的xml格式： j 自定义了一种北京54坐标系的航迹数据，每一个航迹包括了属性和对应的几何点，由于通用的航迹格式GPS只支持经纬度坐标的读取，所以需要自己写脚本，对要素进行处理将xml格