arcpyIDW插值生成等值线等值面
读取grid浓度、温度、风场等数据
IDW插值
插值数据根据插值范围生成等值线
等值线平滑处理
等值线根据插值范围clip生成等值面
等值面与等值线spatial_join赋等值面范围
输出json结果
Tagged: ArcPy
栅格那点儿事（一）—Raster是个啥子东西
Raster是个啥子东西rn rn现如今，不仅是在遥感应用中能看到花花绿绿的影像了，在GIS应用中也能随处看到她们的身影。在各种在线地图中，卫星影像作为底图与矢量的道路层叠加；高程DEM作为高程来源实现地形的山影效果。相信大家对影像这个词，应该一点儿也不感到陌生。rn rnrn那么栅格呢？相对与我们熟悉的影像，这俩只是名词上的区别么？rn栅格（Raster）和影像（Images/Imagery）在
将Visual Studio Code配置成好用的Python IDE
VSCode是一款非常好用的编辑器（或者IDE），具有很好的可扩展性，功能比较强大，占用的系统资源也适中，启动速度较快，而且支持全平台，比较适合作为Python开发用的IDE。本文针对Linux（主要是Ubuntu，其他发行版类似），整合一些Python开发相关的配置，仅供刚入坑Linuxer参考。一、VSCode与其他编辑器（或IDE）的比较（1）VSCode与Atom的比较：Atom是一款由g…
热力图原理及ArcGIS Pro 制作热力图
热力图（Heat Map）是通过密度函数进行可视化用于表示地图中点的密度的热图。它使人们能够独立于缩放因子感知点的密度。rn (Heat Map):一般是基于离散点、线或面的分析与表达，或者基于连续表面的密度分析得到的热力图或热度图。强调空间位置和基本的空间分布特征。rnrnrnrnMH370 location probability heat maprnrn 这里主要以A
安装PostgresSQl 并创建SDE数据库
本博客编写的软件环境如下 rn操作系统： win7 rnArcGIS版本 10.1rn数据库： postgresql-9.0.6-1-windows-x64 以及连接postgresql数据库的32位客户端的5个dll 下载地址（libeay32.dll，libiconv-2.dll，libintl-8.dll，libpq.dll，ssleay32.dll）。rn rn 之前写在线编
ArcGIS地图发布与数据切片
ArcGIS地图发布与数据切片n 下面介绍如何发布ArcGIS Server服务并对其进行切片n地图发布在发布服务之前首先要对地图服务按照标准进行数据整备，当一个MXD文档中的地图数据整备完毕之后就可以将，MXD文档发布为ArcGIS地图服务了。步骤如下：1、首先点击 文件（File） (Share As) 服务Service 进入发布对话框图6.12、这里有三个选项
使用ArcGIS for Desktop的Select by attribute进行多条件查询的问题解决
测试环境：ArcGIS10.2.2 for Desktop 、ArcSDE10.2.2、Oracle11.2.0.3（个人认为该问题跟环境无关）ce
五分钟学GIS | GPU在GIS技术上的应用
在GIS中图形处理能力尤为重要，特别是在三维GIS技术中，为了让三维场景能够更逼真、更流畅地显示，往往需要为计算机配置一个独立显卡，利用其GPU技术来满足GIS的图形运算需求。
使用python制作ArcGIS插件（6）案例分析
利用ArcPy制作航空制图插件 nBy 李远祥 n这是笔者两年多前写的一个面向航路图做的一个插件，基本上将航路图的制作进行流程化，制作成为可交互的插件，只要有航路和机场的信息，就可以直接生成一个航路图，每个图的时间成本大概为3分钟左右。n整个插件都是使用Python进行制作，其设计思路是做好制图模版，剩下的工作就交给插件结合ArcPy调用ArcGIS工具进行数据的生成和处理，直到最后出图。
将Visual Studio Code配置成好用的Python IDE
VSCode是一款非常好用的编辑器（或者IDE），具有很好的可扩展性，功能比较强大，占用的系统资源也适中，启动速度较快，而且支持全平台，比较适合作为Python开发用的IDE。本文针对Linux（主要是Ubuntu，其他发行版类似），整合一些Python开发相关的配置，仅供刚入坑Linuxer参考。一、VSCode与其他编辑器（或IDE）的比较（1）VSCode与Atom的比较：Atom是一款由g…
密度分析（七） Python实现
密度分析（七） Python实现rnrnR语言的核密度渲染弄完了，今天来看看Python的核密度。本来没准备些Python，但是今天在玩seaborn包的时候，发现了强大的kedplot这样一样工具，让虾神我惊为天人，不写感觉都对不起这个方法。rnrn在讲python的核密度分析工具工具之前，前置条件是先说说seaborn包，这个包是matplotlib的扩展包，基本是在matplotlib上继承
[Python] python基础 – 安装拓展包(pip | 手动) | 虚拟环境virtualenv
python安装pip，修改pip下载源，使用Pip下载拓展包。
sql 中CURSOR 的使用
CURSOR是游标，常用于脚本处理。这里主要介绍自己常用的方法，同时也会把网上的一般格式进行解释。一、游标一般格式：DECLARE 游标名称 CURSOR FOR SELECT 字段1,字段2,字段3,… FROM 表名 WHERE …OPEN 游标名称FETCH NEXT FROM 游标名称 INTO 变量名1,变量名2,变量名3,…WHILE @@FETCH_STATUS=0 …
通过源代码分析Servlet的service方法和doXXX方法
转载自http://blog.csdn.net/xuniailian/article/details/8239751（虚拟恋爱）rnrn首先，对于我们创建的任意一个Servlet都必须要实现Servlet接口，GenericServlet是一个通用的、不特定与任何协议的Servlet，它实现了Servlet接口，而HttpServlet继承于GenericServlet，因此HttpServlet
R-ArcGIS探秘（1）安装以及Sample运行
在今年的全球用户大会上，Esri官方公布了R-ArcGIS的官方示例，在ArcMap和ArcGISrn pro中，直接通过Toolbox可以调用R的分析工具包，将R的分析能力直接作用在ArcGIS上面。打通了两大工具之间的壁垒，不能不说，是分析界的一大盛事啊。rn rn视频截图如下：rnrnrn rn总所周知，R在数据分析领域的地位，犹如JAVA在编程语言领域的地位，简直就占据了开源软件的头把交椅
导出mxd专题图到图片
直接上代码：
python2.7 arcpy安装pandas
网上说的镜像源全部没有成功。nn下面的步骤不一定要全部尝试，楼主的成功方式nn卸载pip，更换为pip18版本nn下载pandas 0.12版本，不要太高，因为依赖于numpy的版本，升级卸载此包均失败，所以只有妥协。nn使用python setup.py instasll方式安装，exe安装的import就会报错nn然后安装pandas会提示各种包缺失，只有一个一个在pypi上下载下来再安装nn…
ArcPy遍历Dataset中的FeatureClass添加GLOBALID字段
某项目的要素类按不同的数据集存放，现需为每个要素类添加一个唯一标识符字段（若已有则跳过）。代码如下：env.workspace = r’D:\RK20180201.gdb’nprint env.workspacendss = arcpy.ListDatasets()nfor ds in dss:n fcs = arcpy.ListFeatureClasses(feature_dataset …
（十一）arcpy开发&利用arcpy实现在arcgis中要素类的遍历与记录数获取
在arcpy中可以使用GetCount_management，该函数可以返回表记录数，但是需要注意的是该返回结果返回类型为Result类型，具体使用type函数打印的结果如下所示。nntype a：
arcpy更换符号系统
#原理是基于一个图层的符号系统去替换目标图层中的符号系统nnimport arcpy.mapping as mappingnmxd = mapping.MapDocument(“CURRENT”)ndf = mapping.ListDataFrames(mxd, “Crime”)[0]nupdateLayer = mapping.ListLayers(mxd,”Crime Density by S…
arcpy提取线段/道路起点，终点
# encoding: utf-8nimport arcpyndata=r”D:\删减后的道路网\delbjroad.shp”nls=[]nstrout=””nwith arcpy.da.SearchCursor(data, [“SHAPE@”,”ID”]) as cursor:n for row in cursor:n for row in cursor:n …
利用Arcpy发布地图服务，制作切片
1. 制作地图文档（*.mxd）;n2. 发布地图文档;n3. 制作服务器缓存;n4. 生成切片；n5. 打包成zip
下载文件时空格变加号问题
URLEncoder.encode(fileName,”UTF-8″).replaceAll(“//+”, “%20”)
由坐标数据生成点SHP文件，并由多组点生成线SHP文件（上）
由坐标数据生成点SHP文件，并由多组点生成线SHP文件（上）在最近的公司项目中，老板要求将外业采集的管线连接点坐标表格数据生成线状的shp文件，上网查了下都是生成点状shp文件的例子，作为一个小白根据例子苦逼摸索了半个星期终于搞定了，放上来和大家分享一下。实现思路n连接Excel表格，将表格中需要的数据生成DataTablen获取DataTable中的坐标数据，根据要求获取坐标对n创建空的点shp文
如何制作nc文件（小插曲）
借处理CMORPH降水数据，写一篇关于python制作nc文件的文章。制作nc的模块Matlab， R ，Python 都有，任意选择一个使用即可。之前尝试看过CDO软件，是一个很强大的软件，命令简洁高效，有兴趣可自行google。https://wenku.baidu.com/view/867650d7647d27284a7351ff.htmlnn原始矩阵是一个1440*720的二维矩阵，以这个…
MATLAB读入并处理全球DAT格式气温数据
该脚本读取从 https://crudata.uea.ac.uk/cru/data/hrg/下载的全球.dat格式的气温数据。网站对源数据进行了说明：360*720是全球一个月的数据，然后第二个月接在后面。本次读的是1901-1910年全球气温数据，所以dat文件为43200行*720列。
ArcPy数据处理应用之批量字段赋值
ArcPy处理代码如下： … n = 0 … arcpy.env.workspace=r’D:/data.g…
python 报错’tuple’ object does not support item assignment
ntuple是只可读，而不支持写的，当试图去修改元组的数据时就会报此错误n
Sqlite数据的制作
移动端数据本地化解决方案之一
ArcGIS小知识（三）——矢量面转点、点转面（arcgis），面矢量转点矢量
有时需要把面矢量打散成散点图，有时又需要合并有序点生成一个完整的面矢量，就是这篇的主要内容了。nn我用这个功能是先将面转点，而后在转回面。因为点转面时必须导入某一格式的txt才能识别并将相同图斑的点聚集在一起。因此，在面转点时需要分两步进行：首先将面矢量打散得到折点，并导出其相应的属性和序号信息；然后将点矢量导出为特殊格式的txt，然后利用Excel进行分块。nn一、面转点nn1.1、打散得到折点…
arcpy spatial adjust更改构成面.线的坐标点
之前写过更改点的坐标点和读取构成线面的坐标nnn# encoding: utf-8nimport arcpyndata=r”d:\JZ\jz\jianzhuclip.shp”nfeatureList=[]npoint = arcpy.Point()narray = arcpy.Array()nfeacount=0nnwith arcpy.da.UpdateCursor(data, [“SHAPE…
ArcGIS Python arcpy 批量创建SDE空间索引
1、 SDE连接文件建立使用Desktop系统工具箱建立SDE连接文件，如果catalog已经直连到SDE可以省略此步骤。2、Python脚本导入Toolbox设置脚本中使用变量为SDE表空间，在导入脚本过程中只需要指定这一个变量即可（1）新建工具箱（2）添加脚本（3）变量设置输入名称，在数据类型下拉列表中选择“工作空…
GSMaP 降水数据下载指南
1打开官网：http://sharaku.eorc.jaxa.jp/GSMaP/index.htm 点击注册，填好邮箱等。nnnn2：收到确认邮件，会给一个ftp地址，进入即可下载。nn3：可提供数据一览nnnn注意;nn00Z-23Z代表从0点到23点的降水。根据需要，自行选择。nn …
ArcGIS批量添加字段 AddField
前些日子做了一个数据库的修改工作，要为数据添加字段，数据为Geodatabase类型，里面有90多个图层，每个图层添加6个字段，一共要做2个小时左右，而且效率比较低，容易出错，因此想到了通过程序来解决。rnrn 实现如下，首先先建立程序主窗体，包括TocControl，AxMapControl和ToolBarControl，设置好他们之间的Buddy属性，在ToolBarControl
在Anaconda2中使用arcpy站点包
由于arcpy非开源且只能在ArcGIS环境下使用，那么自然无法利用Anaconda下载使用了，但当有ArcGIS环境时利用arcpy确实非常方便，因此便可以尝试在Anaconda中导入arcpy。
用Arcpy批量对raster文件进行区统计
import arcpynfrom arcpy import envnfrom arcpy.sa import *nimport arcpy.da as dannnarcpy.CheckOutExtension(“spatial”)nenv.overwriteOutput=1nnworkingDir=”D:\\lyj\\code\\”nenv.workspace
ArcPy批量掩膜裁剪栅格/图像
我在测绘院实习的时候，根据某些实际的工作需求，从批量裁剪矢量图层想到也想搞一个栅格/影像批量裁剪，在ArcGIS中掩膜裁剪栅格/图像是非常好用的，掩膜裁剪工具很好用，既能矢量裁剪，也能栅格/图像裁剪。实际工作中，可能会面临大批量数据处理，ArcGIS现有的工具箱不能满足要求的情况下，不得不想想构建模型迭代器，或者Python脚本用以处理批量处理的问题了。nn假如是栅格/影像标…
SELECTION-SCREEN选择界面字段值的实时获取-函数DYNP_VALUES_READ
项目中ABAP程序会遇到一种需求：实时获取屏幕字段的值。如下图所示：nnnn根据销售组织/分销渠道/产品组来限制凭证类型的选择，类似与系统VA01操作。因为存在大量的必输字段，导致一个问题：只有把所有的必输字段维护完成之后，程序才会经过AT SELECTION-SCREEN的节点，但是实际的业务需求：未维护完整其他必输的情况下，也可以约束到凭证类型的选择范围。nn为此，我们首先采用了MEMORY …
Arcgis RasterCalculator
要命nnnprivate void SetGPEnvironment()n n object obj = null;n if (m_GP == null) return;nn m_GP.ClearMessages();n m_GP.OverwriteOutput = base.Overwrite…
ArcGIS中把img图层文件转成tif图层文件的方法
不清楚在ArcGIS中如何把img图层文件转成tif图层文件？
如何导出SHP文件中的点坐标(x,y)
arccatalog建一个personal geodatabase,加入要素，其中就有x y坐标了rnrnarctoolbox里有工具rnfeatures->Add XY coordinatesrnrn1.利用系统自带的方法求出各节点的坐标就是所求坐标：rn1).右击点图层rn2).选择Open Attribute Tablern3).在打开的表格右下角有个Options按钮，点击rn4).选择
获取arcgis克里格插值后的每个点的插值结果
最近在arcgis克里格插值后需要得到插值后的矩阵。找了很久才找到工具箱中的calculateZvalue，但是一次只能计算一个点，因此需要用脚本去提取。nnimport arcpynarcpy.CheckOutExtension(“GeoStats”)nnntotal=0npoints=[]nfor width in range(37578564,37636464,50):
利用ArcPy获取每个Toolbox下面工具的个数
有时，我们想准确的知道各个工具箱工具的数量，怎么办尼，可以通过ArcPy来做，利用arcpy的list函数和一点python代码就可以完成。rnrnrnrnrnrn通过这两个函数就可以遍历所有的工具箱，然后检索某个工具箱的工具了，主要一个注意的地方就是ListTools函数的参数需要时通配符，需要匹配…
arcpy限制返回地图元素
import arcpy.mapping as mappingnmxd = mapping.MapDocument(“CURRENT”)nn#仅仅返回图例nfor el in mapping.ListLayoutElements(mxd,”LEGEND_ELEMENT”,”*Crime*”):n print el.namen#更改布局元素的属性nnimport arcpy.mapping a…
ArcGIS API for Python(2)
在这一篇文章里面，准备简单的使用此API。用例就是绘制朝阳区出租房价格热力图。
栅格那点儿事（一）—Raster是个啥子东西
Raster是个啥子东西rn rn现如今，不仅是在遥感应用中能看到花花绿绿的影像了，在GIS应用中也能随处看到她们的身影。在各种在线地图中，卫星影像作为底图与矢量的道路层叠加；高程DEM作为高程来源实现地形的山影效果。相信大家对影像这个词，应该一点儿也不感到陌生。rn rnrn那么栅格呢？相对与我们熟悉的影像，这俩只是名词上的区别么？rn栅格（Raster）和影像（Images/Imagery）在
ArcGIS中各种合并要素（Unio…
原文地址：ArcGIS中各种合并要素（Union、Merge、Append、Dissolve）的异同点分析作者：傻子舵手nnArcGIS中将两个要素类合并成一个要素有Union、Dissolve、Append、Merge等，在Arctoolbox中均有相应工具，但功能上有所不同：nnUnion：合并输入要素类到新的要素类中。n图形：union只能合并polygon类型的要素类。两个要素类合并时会处
arcpy.mapping 类列表
arcpy.mapping 类可使用地图文档 (.mxd) 或图层文件 (.lyr) 中的不同对象类型的各种方法和属性。此文档可专门用作快速参考。有关详细信息，请使用链接跳转至各帮助页面。nnarcpy.mapping 类 nDataDrivenPagesnn方法
arcgis给每个要素添加坐标字段并自动计算出每个色块的质心xy
n n n nnnnn已经添加的xy字段效果n做法：1 添加x y字段。nnnnnn2 在x字段上右击计算几何nnnnnnnnnnnn选择质心x坐标，点击确定；即可自动生成每一个色块面的x坐标。3 同理生成y坐标即可。n n nn…
使用arcpy.mapping 更新和修复数据源
在许多情况下，您都可能需要修复数据源或重定向数据源至其他位置。然而，如果是在每个相关的地图文档中手动进行更改，则会显得异常麻烦。arcpy.mapping 脚本环境提供了多种方法使得您无需打开地图文档即可自动进行更改。您可针对各个图层逐一更新数据源，也可一次更新同一工作空间中的所
Python实现遥感生态指数计算
最近在做一些遥感相关的图像处理项目，涉及到遥感生态指数的计算。由于项目要求Python实现，搜索互联网关于Python实现的遥感生态指数计算程序资料很少，于是就自己实现了一个并分享在这里，供需要的朋友参考。n首先需要了解遥感生态指数是什么，不是很清楚的朋友可以参考下面的几篇文章：n城市遥感生态指数的创建及其应用n区域生态环境变化的遥感评价指数n基于遥感生态指数的南京市生态变化分析n知道了遥感生态指…
使用ArcGis工具合并图层
1.开始编辑rn 2.选中要合并的图层：根据图层的编号，点击Selected后，再点击Merge，合并rn
Arcpy 帮助文档
https://resources.arcgis.com/zh-cn/help/main/10.2/index.html ArcGIS 帮助（10.2、10.2.1 和 10.2.2）帮助文档
arcpy.mapping常用四大件-Layer
arcpy.mapping常用四大件-Layer by 李远祥 图层是地图里面非常重要的组成，几乎所有的制图显示都与图层有关。首先它是连接数据与符号渲染的重要桥梁，其次在出版地图中，图层又与图例关联在一起，当然，图层最重要的就是组织数据，数据的分类、上下叠加、标注的显示等等，全部都跟图层扯上关系。在arcpy.mapping 中，Layer是一个非常重要的类型，如果说MapDo
ArcPy修改字段值
修改字段值（注册版本时需要开启编辑）rnrnrnrnrnrnrnrnrnrnfile_names=(‘BOUAN’,’BOULK’,’BOUNT’,’BOUPT’,’CTRLK’,’CTRPT’,’HYDAN’,’HYDLK’,’HYDNT’,’HYDPT’,’PIPAN’,’PIPLK’,’PIPNT’,’PIPPT’,’RESAN’,’RESLK’,’RESNT’,’RESPT’,’ROAAN
Arcpy 复制数据库
分别按照要素类、表格、要素数据集、栅格数据集分别将数据存储到新的数据库
Arcpy基础方法之Tools&Toolboxes
rnrnrnc:\python >>>import arcpy rn>>>arcpy.ImportToolbox(“c:/PCovertL.tbx”) rnrnrnc:\pyth
利用shp点文件对影像按指定尺寸进行裁剪
nnnn ArcGIS在安装时已自动安装Python2.7以及arcpy包。在用编辑器编辑的时候需要将编译环境设置为python2.7。本文使用Pycharm进行编辑。设置为python2.7的操作如下…
读取shp的几何坐标
// 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
基于Python脚本的ArcMap字段计算器分类赋值
利用python脚本和ArcMap的字段计算机对属性表中的数据进行分类赋值。