利用python将所有时序NDVI影像转换为DataFrame数据并进行归一化
代码如下：nnnn# coding:utf-8nimport gdalnimport pandas as pdnfrom pandas import DataFramen#首先读取tif文件npath = r”D:\NDVI.tif”ndataset = gdal.Open(path)nim_width = dataset.RasterXSize;#获取宽度nim_height = dataset…
Category: ArcPy
python 使用matplotlib 实现从文件中读取x，y坐标的可视化
1. test.txt文件，数据以逗号分割，第一个数据为x坐标，第二个为y坐标，数据如下：1.1,2n2.1,2n3.1,3n4.1,5n40,38n42,41n43,422. python部分代码#!/usr/bin/pythonn# coding: utf-8nimport matplotlib.pyplot as pltnimport numpy as npnimport matplotli…
ArcGIS教程：图斑整理之字段计算器使用技巧
1.字段计算器简介rn 在数据整理过程中经常要用到对属性表的处理，即为字段进行赋值或运算。字段计算器(Field Calculator)是一个强大的处理字段值的工具，不仅可以实现快速批量赋值，还支持Python和VBScript，可以通过代码进行复杂条件的赋值工作，并且字段计算器还可以在Model Builder中调用，构建空间模型。rn rnrnrn 在某个属性字段的右键菜单中即可调出字段
ArcPy脚本——优化版删除GIS数据库空层（支持GDB,mdb,shp文件夹）
之前有人写了一个批量删除gdb，mdb和shp当中数据为空的图层的脚本，原文参考：https://blog.csdn.net/wl05031/article/details/50476251n 使用过程中发现，某些mdb中的空图层无法删除，该mdb并非是由arcgis创建的。分析mdb发现，使用arcgi…
Anacoda安装—python3.6
1.anaconda官网下载win64版本 nn傻瓜式安装 点击下一步nn记住勾选环境配置nn nn安装完后nn打开Anaconda prompt n输入命令 conda list显示已经安装库nnumpy pandas等重要库nn如果Anaconda官网下载比较慢时，建议采用清华镜像下载n点击Jupyter Notebook 提供浏览器编译环境n新建点击new按钮n复杂调试还需要IDEn还想安…
【基础功能】ArcPy批量栅格转矢量
批量处理100个栅格转矢量，用ACRGIS得弄到手残，所以~不如跳舞~~~哈哈~~~当让是不如码代码rn这是用python语言，Arcpy的强大助力批量让100个名称带有数字，从3-300，间隔为3的tif转为shp矢量文件rnimport arcpynfrom arcpy import envnnenv.workspace = “E:/data/temp/data”nfield
Arcpy中用pathon写一个DEM生成等高线再用面裁剪的Toolbox工具
工具是将DEM生成等高线，然后用clip进行裁剪rnarcpy的代码如下rn# Import system modulesnimport arcpynfrom arcpy import envnfrom arcpy.sa import *nn# 设置环境nenv.workspace = “C:/data”nn# 设置参数值ncontourInterval =
ArcGIS如何获取影像的像素值
首先将栅格数据转为点raster to point , 将每个像元转为点， 之后用sample 或 extract values to points n输入栅格和转化的点数据， 即可生成包含每个像元属性值得表
【arcgis】文件地理数据库，个人地理数据库与ArcSDE的局别
地理数据库的类型rnrnrnrnrn地理数据库是用于保存数据集集合的“容器”。有以下三种类型：rnrnrn文件地理数据库 – 在文件系统中以文件夹形式存储。每个数据集都以文件形式保存，该文件大小最多可扩展至 1 TB。建议使用文件地理数据库而不是个人地理数据库。rn个人地理数据库 – 所有的数据集都存储于 Microsoft Access 数据文件内，该数据文件的大小最大为 2 GB。rnArcS
SELECTION-SCREEN选择界面字段值的实时获取-函数DYNP_VALUES_READ
项目中ABAP程序会遇到一种需求：实时获取屏幕字段的值。如下图所示：nnnn根据销售组织/分销渠道/产品组来限制凭证类型的选择，类似与系统VA01操作。因为存在大量的必输字段，导致一个问题：只有把所有的必输字段维护完成之后，程序才会经过AT SELECTION-SCREEN的节点，但是实际的业务需求：未维护完整其他必输的情况下，也可以约束到凭证类型的选择范围。nn为此，我们首先采用了MEMORY …
脚本加密与解密
此代码可以进行Encode加解密，比如下面这样的asp代码：rnrnrnnn98062–脚本天地nnnn
当ArcGIS10.2遇到SQLite
在刚刚结束的2013Esri开发者大会中，有好多用户都会咨询ArcGIS10.2什么时候推出，可以见得大家对ArcGIS新版本的期待，今天就给大家介绍一下一个新特性，使用SQLite来存储GIS数据。大家都知道ArcGIS支持传统的关系型数据库（oracle、SQL Server、IBM DB2、Informix、PostgreSQL），还支持现在比较火的云数据库（SQL Azure），而且还有一
ArcGIS10.2+ 解决shp文件属性表中文乱码问题
ArcGIS 10.2以后，shp文件属性表或者导出的表格中文会有乱码的情况，如下图所示nnnn解决方案：nn在CMD命令行中，输入一下命令：nnreg add HKEY_CURRENT_USER\Software\ESRI\Desktop10.4\Common\CodePage /v dbfDefault /t REG_SZ /d 936 /fnn版本按照自己安装的写对应的即可nnnn问题分析：…
ImportError: No module named ‘xlrd’
在用pandas处理Excel文件的时候，nnnimport pandas as pdndf = pd.read_excel(‘excelname.xlsx’) # 读取数据nn如上所示读入数据后可能会出现这样的错误提示：nnnImportError: No module named ‘xlrd’nn其原因是，我们需要单独安装一个xlrd模块才能支持对Excel的处理。nn在anaconda中安…
吧啦下 Background Geoprocessing
今天和大家说一下后台地理处理（GP），说到后台，自然会想到前台。rn 其实在 ArcGIS Desktop 9.X 以及更早的年代，ArcToolbox 中的工具运行的模式只有前台模式，也就是大家熟悉的情景：运行某个工具，然后工具运行窗口挡在 ArcMap的前面，窗口中的Log不停的滚动，直至工具运行完毕，ArcMap 才能使用。并且，由于rn ArcMap 和 ArcCatalog
栅格那点儿事（一）—Raster是个啥子东西
Raster是个啥子东西rn rn现如今，不仅是在遥感应用中能看到花花绿绿的影像了，在GIS应用中也能随处看到她们的身影。在各种在线地图中，卫星影像作为底图与矢量的道路层叠加；高程DEM作为高程来源实现地形的山影效果。相信大家对影像这个词，应该一点儿也不感到陌生。rn rnrn那么栅格呢？相对与我们熟悉的影像，这俩只是名词上的区别么？rn栅格（Raster）和影像（Images/Imagery）在
删除GIS数据库空层(GDB,MDB,Shape)
批量删除GIS数据库空层。rn# -*- coding: utf-8 -*-n# —————————————————————————n# Merge.pyn# Created on: 2013-01-21 10:25:22.00000n# (generated by WangLin_TJCH)n
windows 环境下基于Python 的GDAL 安装
最近由于需要利用pytho处理地理空间数据，但是python本身并没有访问和处理地理空间数据的包，只能借助于GDAL（Geospatial Data Abstraction Library）来进行访问和处理。GDAL并非纯净python脚本的包，所以需要通过其他途径进行安装。具体安装步骤如下：rn1.检查windows下python 安装版本，确定以后下载相应的GDAL安装文件。我的python
关于arcgis中按掩膜提取的一个问题
在Arcgis中进行按掩膜提取操作后，有时候会得到比较奇怪的结果，这种奇怪的结果可能会导致后续的操作出错。以下数据是珠江三角洲的土地利用数据。然后使用珠江三角洲的行政边界的shp文件对其进行掩膜提取，会得到如下结果，即右下角会多一块绿色的区域，如果使用识别工具对其进行识别，并不是noData，而是分类的值。造成这种错误结果的原因，是因为上图林地地类和背景值的值都为0，arcgis掩膜提取时，可能在…
1、Python(Anaconda)集成开发环境搭建
Anaconda里面集成了很多关于python科学计算的第三方库，安装方便。由于接下来要使用开源的Python科学计算库NumPy。所以选择了Anaconda（一个开源的Python发行版本）。nn nnAnaconda安装步骤如下：nn（1）、到官网下载Anaconda安装包。包含Windows/Mac/Linux三个版本，可以根据自己需要进行选择。nn https://www.anaconda…
python绘制excel数据热力图
1.首先安装PIP n2.安装插件 xlrd numpy matplotlib（不多做描述）大体思路就是 通过xlrd读取excel的数据，默认存放在list当中，再将list数组转换为python矩阵。通过二维图形显示出来from numpy import *nimport matplotlib.pyplot as pltnimport matplotlib.cm as cmnimport xlr
计算机几何基础，计算多边形的面积, HDU-2036
还记得我们是怎么计算多边形面积的么？nnnn让你先思考10秒钟，Any good idea?nn好吧相信你心中已经有了想法。nn还记得我们学过的叉乘么，两个向量同一起点的叉乘积再乘1/2就是三角形的面积（哦哦哦，我知道，把多边形分成三角形就可以了）nnnnnn然后很自然的我们就能联想到多边形的面积，就是把多边形分成N-2个三角形就能算多边形的面积了nnnn那凹多边形的面积呢？nnnnnn我们计算时…
ArcGIS将带有X、Y坐标的Excel数据转成shp点数据
1、首先，确保Excel中有两列是X、Y坐标，如下图所示：2、打开ArcMap，在工具箱中找到Conversion Tools –>Excel –>Excel To Table工具。3、在工具中选择要导入的Excel文件，选择表导出的位置，和数据所在的sheet页，点击OK。4、数据转换完成后，表会自动添加到TOC窗口中。右击表，选择Display XY Data…选项。5、在…
Arcgis中设置.mxd相对路径
复制拷贝地图文档后再次打开会遇到图层前一个！的标识符号，是因为加载图层的相对路径变化了，需要重新指定图层位置。对此将.mxd文档设置为相对路径起到作用，保证文件到任何地方都可以打开,具体做法如下：rn（1）讲文档和文档含有的数据源都放在一个大文件夹下面；rn（1）File–>Document Properties–>打开对话框中，点击Data Source Options…按钮；rn（2）
python-UnicodeEncodeError: ‘ascii’ codec can’t encode characters in position
python2和python3都是用的办法, 非常推荐的方法:import sysnimport codecsnsys.stdout=codecs.getwriter(‘utf8’)(sys.stdout.detach())适用于python2的还有一种办法:import sysnreload(sys)nsys.setdefaultencoding(‘utf8’)…
如何导出SHP文件中的点坐标(x,y)
arccatalog建一个personal geodatabase,加入要素，其中就有x y坐标了rnrnarctoolbox里有工具rnfeatures->Add XY coordinatesrnrn1.利用系统自带的方法求出各节点的坐标就是所求坐标：rn1).右击点图层rn2).选择Open Attribute Tablern3).在打开的表格右下角有个Options按钮，点击rn4).选择
ArcGIS中把img图层文件转成tif图层文件的方法
不清楚在ArcGIS中如何把img图层文件转成tif图层文件？
在C#中调用Python方法
1. 安装IronPythonnn到http://ironpython.codeplex.com/下载IronPython。nn安装下载下来的安装包（要先装VS啊）。nn nn2. 创建项目nn创建一个C#的控制台应用程序。nn添加引用： 浏览到IronPython的安装目录中，添加对IronPython.dll，Microsoft.Scripting.dll 两个dll的
获取arcgis克里格插值后的每个点的插值结果
最近在arcgis克里格插值后需要得到插值后的矩阵。找了很久才找到工具箱中的calculateZvalue，但是一次只能计算一个点，因此需要用脚本去提取。nnimport arcpynarcpy.CheckOutExtension(“GeoStats”)nnntotal=0npoints=[]nfor width in range(37578564,37636464,50):
在各种WorkSpace中判断要打开的文件(如:栅格,要素)是否存在
用IWorkspace2.NameExists可以判断工作空间是否存在指定名称和类型的数据集..rn rn引用MSDN的例子,获得工作空间的文件名字,历遍:rn//IWorkspace DatasetNames Examplern public void IWorkspace_DatasetNames_Example(IWorkspace workspace) rn r…
ARCGIS中按照某字段属性批量导出shp ，即将一个shp图层分割成单个的shp
如何通过某字段属性批量导出各个字段对应的shp文件？比如我现有一副国家的的行政区划图，想要提取各个省份的shp文件。nn通过Analysis Tools——Extract——Split工具实现。nn步骤如下：nn在“Input Features”和“Split Features”中都选择需要分割的行政区划图，相当于用自己切割自己。在“Split Field”中选择具有唯一标识的字段，例如省级行政区…
ArcGIS10.0操作之11—同时矢量化同一类要素的多个不同图层
ArcGIS10.0中可以同时矢量化同一类要素的多个不同图层，刚发现这个功能，在这里展示一下：rn1、多个图层但都是Polygonrnrn rn2、组织编辑要素模板，如下，弹出组织要素模板对话框rnrn rnrn单击New Template,弹出如下对话框rnrn单击Select All，全选，单击Finish,单击Close,既可以同时编辑这三个去，下图是创建要素对话框，编辑相应图层是，选中相
Python SyntaxError: EOL while scanning string literal
异常描述：SyntaxError: EOL while scanning string literal解决方法：将单引号’改成改成三引号”’（或”””）,示例：
The Geometry has no Z values 解决办法
我们在ArcEngine开发过程中进行创建要素时，简单的IFeatureClass.CreateFeature之后，然后把IGeometry对象赋给Feature.shape，会引发“The Geometry has no Z values”的错误提示。类似的错误在做地图标注过程中也会出现。n参考网上的解决方案如下，nn基本思路就是在CreateFeature之后，设置Z值（顺带将M值也一并
如何将python程序封装成exe可执行文件
如何将python程序封装成exe可执行文件将python程序直接转成exe程序的方法很多，主流的工具有PyInstaller、cx_freeze、py2exe等，论及优缺点，可谓各有千秋。然而，再好的工具在使用的时候也会遇到各种各样的问题。这里无法比较哪个好那个坏，不同工具对不同类型程序的处理能力不同，也许你用这个方法转化出来的程序出问题了，换个方法就行了。本文主要介绍py2exe的安装和使用。p
python使用xlrd读取excel数据时，整数和日期变小数的解决办法
想让excel数据返回一个字符串字典list，造了一个excel。数据如下： 用代码读取后发现整型和日期变成了小数(123456变成了123456.0，2017/5/20 5:20变成了42875.22222222222)： 有个比较简单的解决办法就是在数字和日期的单元格内容前加上一个英文的引号即可。如果数据比较多，也可以批量加英文引号的前缀（网上都有方法）。 …
TypeError: ‘zip’ object is not callable
>>> list=zip(‘1234′,’abc’)n>>> list(list)nTraceback (most recent call last):n File “”, line 1, in n list(list)nTypeError: ‘zip’ object is not callablen>>> l=zip(‘ab’,’123′)n>>> lnn>>> lis
ArcGIS 10.1 for Desktop SP1新特性
1 SP1综述rn 总体来说，ArcGIS 10.1 SP1使得整个ArcGIS产品线在稳定性和功能上都有所改进和增强，以下是软件功能改变的部分内容，详细内容可到微博下载，地址：http://vdisk.weibo.com/s/lwke5rnrn2 Geoprocessing and Arcpyrn 简言之，SP1中地理处理和Arcpy改进的重点内容为
ArcGIS如何获取地理要素的几何属性
ArcGIS如何获取地理要素的几何属性rnrnrn在ArcGIS中，将向输入要素添加新的属性字段以存储各种属性，如长度、面积、四至、中心点或 x、y、z 和 m 坐标。rn首先打开ArcMap，加载需要获取几何属性的要素图层。rn rn打开ArcToolBox工具rnrn在ArcToolBox工具中，找到“数据管理工具”里面的“要素”下面的“添加几何属性”，双击弹出“添加结合属性”窗口rn rn可
批量投影栅格
本篇博客使用modelbuilder来对大批量的栅格文件进行批量投影。nnnnnn链接：https://pan.baidu.com/s/1THc_54tItzew0pt-TXLjsg n提取码：gdjo nn（示例中的模型是基于ArcGIS10.6.1制作的）nn注意事项：nn如果输出位置是要将投影后的栅格数据存储到gdb中，请注意gdb中的文件不能以数字开头。n 如果是将投影后的栅格数据存储到文…
解决PyCharm Please select a valid Python interpreter问题
之前PC上安装的是Python2，后来工作需要转成Python3了。然后在用pycharm运行Python2的程序时发现源程序运行报错（出去语法错误）error:please select a valid Python interpreter问题原因：在pycharm导入源项目的时候没有选择运行Python的程序解决办法：1、更改Pycharm的设置 打开settings（CTRL + ALT +…
解决问题：启动Python时报“ImportError: No module named site”错误
记录Python启动问题解决过程和方法。n解决问题：启动Python时总是报“ImportError: No module named site”这个错误。
ArcGIS 10.2中栅格缺失值处理–Nodata in raster
问题描述：rn需要对同一区域的多图层的栅格数据进行叠加操作，而相同位置，如果有一个图层的值为Nodata，则计算结果也是Nodata,即：Nodata+value = Nodatarnrn需要将Nodata转换成0，以便进行后续操作。rn工具：Spatial Analyst Tools / Map Algebra / Raster Calculationrn首先利用Con() 函数进行条件的限制；
同一点图层内最近点的查询与距离计算
众所周知，使用arcgis toolbox的tools》》proximity》》near可以用来计算点与点之间的最近距离，但是通常是分为两个点图层。比如，一个是居民位置点图层，一个是超市点图层，则可通过上述方法超找到距离每个居民点最近的超市，并计算这个最短距离。rn但是如何查找同一图层内，距离某点最近的点及计算最短距离呢？比如，在居民位置点图层内，查找距离任意一个居民点最近的居民点（有点绕口）。其
点数据转成Arcgis线、面文件
[转载]点数据转成Arcgis线、面文件rn rn(2012-08-02 16:26:58)rn转载▼rnrnrnrnrnrnrn标签： rnrn转载rnrnrn分类： 学习资料rnrnrnrnrnrnrn原文地址：点数据转成Arcgis线、面文件作者：哈哈猫rn一、创建execl点文件rnrn要求格式：第一列为点号，接着两列为X、Y坐标，表格第一行最好为ID、X、Yrn二、导入ARCGISrn点
调用gp工具把shp、mdb、gdb中的table导入到sde（oracle）中
添加这两个引用nnusing ESRI.ArcGIS.ConversionTools;//nnusing ESRI.ArcGIS.Geoprocessor;nn nn调用gp 实现shp入库nn//调用gp工具nnn public void shpFeatureClassToGeodatabase(string inputPath, string sdePath)n n …
ArcPy水文分析（河网分级、流域、集水区）
n好在有model builder和ArcPy，最近刚好用到，写好脚本测试跑了一把，效果不错，做个记录。其中使用的DEM数据是网上下载的aster那个，后面flowacc的自定义参数是90。# -*- coding: utf-8 -*-n# ————————————
ArcGIS中求多边形的“四至点”
ArcGIS10.0中求多边形的四至点可以直接根据字段计算器，利用Python语句直接求得rnminX = !shape.extent.xmin!rnmaxX= !shape.extent.xmax!rnminY= !shape.extent.ymin!rnmaxY = !shape.extent.ymax!
ARCGIS中按某字段属性值分类导出
按字段属性值分类导出数据，比如，左侧表中lev值共有 a b c d 四种，要求按照lev值分类输出 a.shp，b.shp，c.shp，d.shp。也就是把 lev =”a” 的所有的导出为 a.shp，其余类推。
arcgis脚本Python批量裁切与统计实现
首先通过批量裁切为arcgis的extractbymask工具，代码如下：n n其次，利用arcgis的Python脚本的GetRasterProperties_management方法，代码如下：
（拾肆）规划师方向技能：使用ArcGIS批量创建公共服务设施点和最近道路距离之间的连线
应用场景：n进行公共和市政设施布局规划时，如何合理地确定设施的数量、位置、规模和服务范围，以实现设施布局的优化时，利用ArcGIS创建构建基础网络要素集，如何批量将服务
ArcGIS教程：栅格重分类的类型
按单个值进行重分类rn “重分类”工具可更改栅格中的各个值。按单个值进行重分类工具将以一对一的方式将一个值更改为另一个值。例如，执行鹿栖息地分析时，土地利用栅格中的各个值表示土地利用的不同类型，需要将其更改为一个优先级范围(如，1 到10)，从而使各种土地利用类型对鹿而言均具有一定意义。将鹿相对偏爱的土地类型重分类为较高值，而那些不受偏爱的类型则重分类为较低值。例如，森林土地利用类型将重分类为
PyCharm Process finished with exit code -1073741819 (0xC0000005)
错误如下n Process finished with exit code -1073741819 (0xC0000005)nn1. h5py这个坑n之前使用TensorFlow好好的，为了装cv2不小心装了一个h5py包，之后死活报上面的错误，解决方案是卸载该包：nC:\Users\XXX>pip uninstall h5pynUninstalling h5py-2.7.0:n Woul…
（三十）arcpy开发&pycharm导入arcpy
nn在编写arcpy脚本之前我们需要将arcpy环境设置好。下面我们来设置一些arcpy需要的环境。nnn 从创建工程开始路径设置n n工程导入arcpy站点包，首先我们使用pycharm新创建一个工程,如下图所示。nnnn在上图中我们可以看到有一个闪电号的标记【No Python interpreter selected】，…
如何去除DEM数据中的缺失值（NoData）
最近在工作中遇到这样一个问题：DEM数据中有Nodata(我们这里所说的情况并非数据边缘部分的nodata,而是数据中间类似数据缺失的情况),见附图1.n 那么，在这种情况下，如何消除nodata,n 使这个数据保持完整性呢？nnnnnnnn其实，利用Arcgis软件能很容易的解决这个问题。这里，我们主要用到con和nibble两个命令。可能大家对Con命令已经非常熟悉了，那么
（4）ArcGIS 10.2 去除要素的Z值和M值
nn1.引言n2. ArcGIS 去除Z值和M值nnnnnnnn1.引言nn ArcGIS要素类是可以包含Z值M值的，但是我们在一些环境中是不需要该值的，该值的存在可能会产生一些bug现象，因此在本篇博客中主要介绍，如何去除要素类中的Z值和M值。nnnn2. ArcGIS 去除Z值和M值nnn打开toolbox的数据管理工具->要素->复制要素n…
python 读取线 的shp文件
import shapefilernsf = shapefile.Reader(“E:\\1.2\\cs\\DX_CSL.shp”)rnshapes = sf.shapes()rnprint shapes[1].partsrnprint len(shapes) #79条记录rn#print len(list(sf.iterShapes())) #79条记录rn#for name in dir(sh
ArcGIS教程：图斑整理之字段计算器使用技巧
1.字段计算器简介rn 在数据整理过程中经常要用到对属性表的处理，即为字段进行赋值或运算。字段计算器(Field Calculator)是一个强大的处理字段值的工具，不仅可以实现快速批量赋值，还支持Python和VBScript，可以通过代码进行复杂条件的赋值工作，并且字段计算器还可以在Model Builder中调用，构建空间模型。rn rnrnrn 在某个属性字段的右键菜单中即可调出字段
ArcGIS属性表中对某一字段进行分段…
需求描述：现有一个矢量图层，其中ACRES字段记录了每一个面要素的面积，如何分段对该字段求和，如该字段数值在以下范围时0-50,50-300,300-1000,>4000，其面积分别是多少？nnnn解决方案：nn（1） n对于数据量少的，可以采用此方法；nnn可以通过属性查询（select bynattributes），通过表达式：ACRES >0 AND ACRES 将分段范围的要素类
【arcgis】文件地理数据库，个人地理数据库与ArcSDE的局别
地理数据库的类型rnrnrnrnrn地理数据库是用于保存数据集集合的“容器”。有以下三种类型：rnrnrn文件地理数据库 – 在文件系统中以文件夹形式存储。每个数据集都以文件形式保存，该文件大小最多可扩展至 1 TB。建议使用文件地理数据库而不是个人地理数据库。rn个人地理数据库 – 所有的数据集都存储于 Microsoft Access 数据文件内，该数据文件的大小最大为 2 GB。rnArcS
