当ArcGIS10.2遇到SQLite
在刚刚结束的2013Esri开发者大会中，有好多用户都会咨询ArcGIS10.2什么时候推出，可以见得大家对ArcGIS新版本的期待，今天就给大家介绍一下一个新特性，使用SQLite来存储GIS数据。大家都知道ArcGIS支持传统的关系型数据库（oracle、SQL Server、IBM DB2、Informix、PostgreSQL），还支持现在比较火的云数据库（SQL Azure），而且还有一
Category: ArcPy
利用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…
脚本加密与解密
此代码可以进行Encode加解密，比如下面这样的asp代码：rnrnrnnn98062–脚本天地nnnn
（十六）arcpy开发&利用arcpy在arcgis中创建gdb文件
本次我们学习的是如利用arcgis提供的python接口arcpy在arcgis中创建gdb文件。在代码中，多处利用检查机制，首先检查路径传入的路径是否存在，如果传入的需要创建的gdb文件已经存在，那么代码会说明已经存在了。当路径已经创建好后，接下来要创建的gdb文件检查是否创建好了。接下来我们看一下实现的具体代码。
（二十五）arcpy开发&实现对面和线要素端点的统计
这个例子是学习使用arcpy来统计面要素和线要素的个数，以及顶点（端点）个数的。利用到pointCount以及令牌方面的知识。首先还是对要素进行了遍历，然后统计相应的顶点数。下面我们来看一下实现代码。
ArcGIS如何获取影像的像素值
首先将栅格数据转为点raster to point , 将每个像元转为点， 之后用sample 或 extract values to points n输入栅格和转化的点数据， 即可生成包含每个像元属性值得表
欢迎使用CSDN-markdown编辑器
欢迎使用Markdown编辑器写博客本Markdown编辑器使用StackEdit修改而来，用它写博客，将会带来全新的体验哦：nMarkdown和扩展Markdown简洁的语法n代码块高亮n图片链接和图片上传nLaTex数学公式nUML序列图和流程图n离线写博客n导入导出Markdown文件n丰富的快捷键n快捷键n加粗 Ctrl + B n斜体 Ctrl + I n引用 Ctrl
windows 环境下基于Python 的GDAL 安装
最近由于需要利用pytho处理地理空间数据，但是python本身并没有访问和处理地理空间数据的包，只能借助于GDAL（Geospatial Data Abstraction Library）来进行访问和处理。GDAL并非纯净python脚本的包，所以需要通过其他途径进行安装。具体安装步骤如下：rn1.检查windows下python 安装版本，确定以后下载相应的GDAL安装文件。我的python
Arcgis中设置.mxd相对路径
复制拷贝地图文档后再次打开会遇到图层前一个！的标识符号，是因为加载图层的相对路径变化了，需要重新指定图层位置。对此将.mxd文档设置为相对路径起到作用，保证文件到任何地方都可以打开,具体做法如下：rn（1）讲文档和文档含有的数据源都放在一个大文件夹下面；rn（1）File–>Document Properties–>打开对话框中，点击Data Source Options…按钮；rn（2）
1、Python(Anaconda)集成开发环境搭建
Anaconda里面集成了很多关于python科学计算的第三方库，安装方便。由于接下来要使用开源的Python科学计算库NumPy。所以选择了Anaconda（一个开源的Python发行版本）。nn nnAnaconda安装步骤如下：nn（1）、到官网下载Anaconda安装包。包含Windows/Mac/Linux三个版本，可以根据自己需要进行选择。nn https://www.anaconda…
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问题分析：…
计算机几何基础，计算多边形的面积, HDU-2036
还记得我们是怎么计算多边形面积的么？nnnn让你先思考10秒钟，Any good idea?nn好吧相信你心中已经有了想法。nn还记得我们学过的叉乘么，两个向量同一起点的叉乘积再乘1/2就是三角形的面积（哦哦哦，我知道，把多边形分成三角形就可以了）nnnnnn然后很自然的我们就能联想到多边形的面积，就是把多边形分成N-2个三角形就能算多边形的面积了nnnn那凹多边形的面积呢？nnnnnn我们计算时…
删除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
关于arcgis中按掩膜提取的一个问题
在Arcgis中进行按掩膜提取操作后，有时候会得到比较奇怪的结果，这种奇怪的结果可能会导致后续的操作出错。以下数据是珠江三角洲的土地利用数据。然后使用珠江三角洲的行政边界的shp文件对其进行掩膜提取，会得到如下结果，即右下角会多一块绿色的区域，如果使用识别工具对其进行识别，并不是noData，而是分类的值。造成这种错误结果的原因，是因为上图林地地类和背景值的值都为0，arcgis掩膜提取时，可能在…
arcgis for python 使用os.walk 和da.walk读取gdb(mdb)
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…
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改进的重点内容为
（六）arcpy开发&利用arcpy在arcgis中批量将hdf数据转tif数据(python文件的变量，数据转换)
本次我们来学习一下关于使用arcpy将hdf转tif数据。hdf格式数据还是常见的。那么什么事hdf格式的数据呢，来看一下相关资料。
ArcPy脚本——优化版删除GIS数据库空层（支持GDB,mdb,shp文件夹）
之前有人写了一个批量删除gdb，mdb和shp当中数据为空的图层的脚本，原文参考：https://blog.csdn.net/wl05031/article/details/50476251n 使用过程中发现，某些mdb中的空图层无法删除，该mdb并非是由arcgis创建的。分析mdb发现，使用arcgi…
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
吧啦下 Background Geoprocessing
今天和大家说一下后台地理处理（GP），说到后台，自然会想到前台。rn 其实在 ArcGIS Desktop 9.X 以及更早的年代，ArcToolbox 中的工具运行的模式只有前台模式，也就是大家熟悉的情景：运行某个工具，然后工具运行窗口挡在 ArcMap的前面，窗口中的Log不停的滚动，直至工具运行完毕，ArcMap 才能使用。并且，由于rn ArcMap 和 ArcCatalog
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’)…
Python之Arcpy按位置选择批处理(SelectLayerByLocation)
工作中，往常使用FME批量压盖选取不同城市的Link、Node数据（存储于GDB地理数据库），能否通过Python引用Arcpy顺利实现上述功能？
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、在…
ArcGIS10.0操作之11—同时矢量化同一类要素的多个不同图层
ArcGIS10.0中可以同时矢量化同一类要素的多个不同图层，刚发现这个功能，在这里展示一下：rn1、多个图层但都是Polygonrnrn rn2、组织编辑要素模板，如下，弹出组织要素模板对话框rnrn rnrn单击New Template,弹出如下对话框rnrn单击Select All，全选，单击Finish,单击Close,既可以同时编辑这三个去，下图是创建要素对话框，编辑相应图层是，选中相
table总结insertRow、deleteRow
表格有几行: var trCnt = table.rows.length; (table为Idrn)rnrn每行有几列：for (var i=0; i
ARCGIS中按照某字段属性批量导出shp ，即将一个shp图层分割成单个的shp
如何通过某字段属性批量导出各个字段对应的shp文件？比如我现有一副国家的的行政区划图，想要提取各个省份的shp文件。nn通过Analysis Tools——Extract——Split工具实现。nn步骤如下：nn在“Input Features”和“Split Features”中都选择需要分割的行政区划图，相当于用自己切割自己。在“Split Field”中选择具有唯一标识的字段，例如省级行政区…
如何将python程序封装成exe可执行文件
如何将python程序封装成exe可执行文件将python程序直接转成exe程序的方法很多，主流的工具有PyInstaller、cx_freeze、py2exe等，论及优缺点，可谓各有千秋。然而，再好的工具在使用的时候也会遇到各种各样的问题。这里无法比较哪个好那个坏，不同工具对不同类型程序的处理能力不同，也许你用这个方法转化出来的程序出问题了，换个方法就行了。本文主要介绍py2exe的安装和使用。p
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…
Python SyntaxError: EOL while scanning string literal
异常描述：SyntaxError: EOL while scanning string literal解决方法：将单引号’改成改成三引号”’（或”””）,示例：
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 =
【基础功能】ArcPy批量栅格转矢量
批量处理100个栅格转矢量，用ACRGIS得弄到手残，所以~不如跳舞~~~哈哈~~~当让是不如码代码rn这是用python语言，Arcpy的强大助力批量让100个名称带有数字，从3-300，间隔为3的tif转为shp矢量文件rnimport arcpynfrom arcpy import envnnenv.workspace = “E:/data/temp/data”nfield
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
The Geometry has no Z values 解决办法
我们在ArcEngine开发过程中进行创建要素时，简单的IFeatureClass.CreateFeature之后，然后把IGeometry对象赋给Feature.shape，会引发“The Geometry has no Z values”的错误提示。类似的错误在做地图标注过程中也会出现。n参考网上的解决方案如下，nn基本思路就是在CreateFeature之后，设置Z值（顺带将M值也一并
在各种WorkSpace中判断要打开的文件(如:栅格,要素)是否存在
用IWorkspace2.NameExists可以判断工作空间是否存在指定名称和类型的数据集..rn rn引用MSDN的例子,获得工作空间的文件名字,历遍:rn//IWorkspace DatasetNames Examplern public void IWorkspace_DatasetNames_Example(IWorkspace workspace) rn r…
在C#中调用Python方法
1. 安装IronPythonnn到http://ironpython.codeplex.com/下载IronPython。nn安装下载下来的安装包（要先装VS啊）。nn nn2. 创建项目nn创建一个C#的控制台应用程序。nn添加引用： 浏览到IronPython的安装目录中，添加对IronPython.dll，Microsoft.Scripting.dll 两个dll的
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中安…
ArcGIS如何获取地理要素的几何属性
ArcGIS如何获取地理要素的几何属性rnrnrn在ArcGIS中，将向输入要素添加新的属性字段以存储各种属性，如长度、面积、四至、中心点或 x、y、z 和 m 坐标。rn首先打开ArcMap，加载需要获取几何属性的要素图层。rn rn打开ArcToolBox工具rnrn在ArcToolBox工具中，找到“数据管理工具”里面的“要素”下面的“添加几何属性”，双击弹出“添加结合属性”窗口rn 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点
（拾肆）规划师方向技能：使用ArcGIS批量创建公共服务设施点和最近道路距离之间的连线
应用场景：n进行公共和市政设施布局规划时，如何合理地确定设施的数量、位置、规模和服务范围，以实现设施布局的优化时，利用ArcGIS创建构建基础网络要素集，如何批量将服务
（三十）arcpy开发&pycharm导入arcpy
nn在编写arcpy脚本之前我们需要将arcpy环境设置好。下面我们来设置一些arcpy需要的环境。nnn 从创建工程开始路径设置n n工程导入arcpy站点包，首先我们使用pycharm新创建一个工程,如下图所示。nnnn在上图中我们可以看到有一个闪电号的标记【No Python interpreter selected】，…
解决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 +…
ArcGIS 10.2中栅格缺失值处理–Nodata in raster
问题描述：rn需要对同一区域的多图层的栅格数据进行叠加操作，而相同位置，如果有一个图层的值为Nodata，则计算结果也是Nodata,即：Nodata+value = Nodatarnrn需要将Nodata转换成0，以便进行后续操作。rn工具：Spatial Analyst Tools / Map Algebra / Raster Calculationrn首先利用Con() 函数进行条件的限制；
Python2和Python3之间的str处理方式导致乱码的讲解
Python字符串问题nn在arcpy中版本为 python2.xnn在QGIS中版本为 python2.x 或者 python3.xnnpython2 和python3 之间的str处理方式经常会导致乱码，故出此文nnpython3版本nnnnpython2版本nn- 在python2版本中使用unicode方式nnnn备注nn在python中不管任何版本，
python gdal投影坐标系转换问题
要将xian80地理坐标系转换成投影坐标系：nnxian1980 = “””nGEOGCS[“GCS_Xian_1980”,n DATUM[“Xian_1980”,n SPHEROID[“Xian_1980”,6378140.0,298.257]],n PRIMEM[“Greenwich”,0.0],n UNIT[“Degree”,0.017453292519943
Python中如何创建文件夹和嵌套文件夹的方法
1、mkdir()—》用于创建单个文件夹。rn f=os.mkdir(“e:\\wuxiaobing”)rn2、makedirs–》创建嵌套的文件夹rnrn f=os.makedirs(“e:\\wuxiaobing1\\1\\2\\”)
错误AttributeError: ‘dict’ object has no attribute ‘iteritems’
Python3.5中：iteritems变为items找到iteritems所在的位置，改成items
arcgis脚本Python批量裁切与统计实现
首先通过批量裁切为arcgis的extractbymask工具，代码如下：n n其次，利用arcgis的Python脚本的GetRasterProperties_management方法，代码如下：
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水文分析（河网分级、流域、集水区）
n好在有model builder和ArcPy，最近刚好用到，写好脚本测试跑了一把，效果不错，做个记录。其中使用的DEM数据是网上下载的aster那个，后面flowacc的自定义参数是90。# -*- coding: utf-8 -*-n# ————————————
调用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 …
SyntaxError: EOL while scanning string literal问题及解决方法
path1 = ‘F:\wamp1\www\stats\’如上面代码所示，在路径的末尾有反斜杠‘\’的时候，会提示你一条错误：【SyntaxError: EOL while scanning string literal】这应该是因为转义字符的问题。所以尝试了一下几种修改方式。正确且是我们想要的： n输入：path1 = ‘F:\\wamp1\\www\\stats\\’nprint(path1)输
ArcGIS中求多边形的“四至点”
ArcGIS10.0中求多边形的四至点可以直接根据字段计算器，利用Python语句直接求得rnminX = !shape.extent.xmin!rnmaxX= !shape.extent.xmax!rnminY= !shape.extent.ymin!rnmaxY = !shape.extent.ymax!
DEM提取一定高程值范围像元
最近在做水文分析，需要提取一定高程一下的像元，但是碰到个问题。rn如果直接利用栅格计算器来提取，系统自动赋值。比如“DEM”rnrnrn那么该怎么做呢？rn做了很久发现了这么一个问题，原始DEM是没有属性表的。所以想着可以先在Raster工具下给DEM附上属性值。但是在操作过程中发现了报错，说是只接受单波段图层。但是我这DEM的确是单波段的，所以非常疑惑。rn接下来，发现此DEM数据是32位的，所
解决Python调试OSError: [WinError 193] %1 不是有效的 Win32 应用程序
Python调试出现OSError: [WinError 193] %1 不是有效的 Win32 应用程序，一般来讲是因为python位数和dll位数不一样导致的，dll是32位的，而python是64位的。nn我的系统是64位Win7安装了64位的python，当调用32位dll（比如32位的ODBC驱动）的时候，就报错。nn如何解决：nn1、安装64位的dll（一般很难找到）nn2、安装32…
python连接SQLServer数据库，执行给定的查询SQL语句，并返回查询结果数据
今天在使用SQL语句做查询任务的时候发现，在数据库中可以直接执行的语句，放在程序中结果却是None，于是乎，就实践了一下，下面是具体的实践。nnn#!usr/bin/env pythonn# encoding:utf-8nnn”’n__Author__:沂水寒城n功能：执行给定的查询SQL语句，并返回查询结果数据n”’nnnimport sysnimport pymssqlnreloa…
将Excel数据导入arcgis中，并进行插值分析（根据经纬度坐标匹配）
最近在准备论文开题，用到了腾讯定位大数据，数据是腾讯定位平台上爬取下来的，整理成Excel格式，导入arcgis中进行分析nn数据刚爬取下来是这样的，存储方式是txt格式，而且腾讯定位数据是5分钟更新一次，所以爬取下来以后还要对数据进行整理，nn比如想获取8：:00~9:00的数据，就要将每个点8：:00~9:00每隔5分钟的定位次数相加，然后整理到Excel表格中，这部分我还没nn想好要怎么处理…
（二）ArcGIS Server之发布切片地图服务
什么是切片地图服务n1切片地图服务优点n2切片地图服务缺点nnn利用ArcMap发布切片地图服务n1用ArcMap打开我们想要发布为服务的数据n2为了发布出来的数据好看一点我们可以适量的将我们的地图数据进行符号化可选此处省略n3利用ArcMap连接ArcGIS Servern4修改地图服务属性修改为切片地图服务n5发布服务
sql按字段值进行统计
用group byrn如有个student表里有性别sex来统计rnselect sex,count(*) from student group by sex;
两条线段求交点编程
今天需要验证两条线段是否相交，突然发现两点求一条直线不会了，然后我就自己在草稿上求了一番。rn已知两个线段分别是（x1,y1）(x2,y2)以及（x3,y3）(x4,y4)，假设这两条线段不平行。rn设线段所在的直线为ax-by+c=0rna1=y2-y1,b1=x2-x1,c1=b1*y1-a1*x1rna2=y4-y3,b2=x4-x3,c2=b2*y3-a1*x3rnd=a1*b2-a2*b
