工作中需要做大量图层的拼接,为了在拼接完成后还能知道原始数据文件是什么,所以写了个Python脚本对每个图层的SOURCE字段进行赋值. 附上Python代码: # -*- coding: utf-8 -*- # nightroad import sys import arcpy path = "C:/Users/nightroad/Desktop/SHP" arcpy.env.workspace = path ff = arcpy.ListFeatureClasses() for f

PDO 查询mysql返回字段int变为String型解决方法使用PDO查询mysql数据库时,执行prepare,execute后,返回的字段数据全都变为字符型. 例如id在数据库中是Int的,查询后返回是String型. 对于PHP这种弱类型的语言,影响不大.在做API返回数据时,如果类型与数据库不一致,对于Java和Objective C这些强类型,影响就很大了.<pre><?php$pdo = new PDO($dsn, $user, $pass, $param); // 在创建

这两天遇到一个奇怪的问题,在整理项目中的建筑物数据时发现,有几个图层进行字段添加后出现问题,shp文件被损坏了.这问题很隐蔽,给shp添加字段后不报错,进行赋值,报错如下: 但是无论是选择“是”还是“否”都可以赋值成功,直到我对各个整理好的图层进行merger才发现,报如下错误: 即JN04_B图层无打开,不进行merger操作,只是在ArcMap中是看不出来的,只有重新打开这个shp才会发现.重新打开这个损坏的shp,会报下面的错误: 通过多次尝试,得出的结论是,JN04_B.shp在进行添加

利用代码可以进行批量处理,安装有10.5及以上版本ArcGIS可以使用工具Split by attributes完成上述任务 # -*- coding: utf-8 -*- # Import system modules import re import arcpy from arcpy import env # function for replace special characters def validateTitle(title): rstr = r"[\/\\\:\*\?\"

使用一个polygon矢量提取某个文件夹中所有的tif格式栅格数据 (要确保先安装好arcpy包) import arcpy arcpy.CheckOutExtension("spatial") arcpy.gp.overwriteOutput=1 arcpy.env.workspace = "E:\\data" #所有栅格影像所在文件夹 rasters = arcpy.ListRasters("*", "tif") mask

excel导入数据进行新增时,发现安装高度和可视距离在数据库创建都是double类型 程序跑完,执行成功后,数据库的数据是2,小数点后的数据没有了 打印sql并执行后发现sql并没有错误, 检查数据库结构,发现创建的height字段为double(4,0),4位长度,小数点后0位,所以没办法把0.5存入 修改数据库结构为double(4,2),保存后再次执行sql,成功存入2.5

注意,ArcGis中shp属性表至少要保留除OID(FID).Shape之外的1个字段. 例:批量删除“省级行政区域”除“FID”.“Shape”.“Area”之外的字段. import arcpy #如果要素类已加载至工作空间,可以直接使用其名字(此处是"省级行政区域"),否,则使用.shp文件的全路径. #code source: https://www.cnblogs.com/yzhyingcool/ QQ:975601416 fieldList=arcpy.ListFields

添加一个新的字段到shp文件中,并且从Excel里导入数据到该字段.原shp文件里的字段ID应该与Excel里的字段ID一一对应才能正确的导入.下图分别是shp的字段和Excel的字段 将class字段添加到shp中去: (1)从Excel中读取数据(为了读取方便,存为.csv或者txt文件) QStringList readFromCSV(QString mfilename) { QStringList readlist; if (mfilename !="") { QFileInf

目录: 1. 同一图层的多个要素合并(2种方法) 2. 导入带经纬度坐标的 txt 文件 3. 栅格计算器的实现 4. 添加 shp 文件(显示在 ArcMap) 5. 通过经纬度坐标生成 Polygon 6. 多边形重叠区域次数计算(赤潮发生频次) 7. 批量删除/添加 shp 文件字段 8. 批量修改文件的 Spatial Reference 9. 读取折线或面几何(点坐标) 10. 读取 Geometry(2种方法) 11. 批量修改图例显示内容 12. 批量引用图层样式(Symbolog

目录: 一.相关技巧 二.工具说明 一.相关技巧 技巧1:将工具从工具箱拖拽到 Python 窗体中自动会生成相应的函数,于是可以快速定位函数名称! 技巧2:将通过工具箱实现的操作结果拖拽到 Python 窗体,可以自动生成处理的代码,于是可以快速了解相关参数含义! 二.工具说明 01   CalculateStatistics 计算栅格数据集或镶嵌数据集的统计数据. 02   MakeXYEventLayer 根据源表中定义的X和Y坐标创建新的点要素图层. 03   Idw 使用反距离加权法

一.背景 沟壑密度是描述地面被水道切割破碎程度的一个指标.沟壑密度是气候.地形.岩性.植被等因素综合影响的反映.沟壑密度越大,地面越破碎,平均坡度增大,地表物质稳定性降低,且易形成地表径流,土壤侵蚀加剧.因此,沟壑密度的测定,对于了解区域地形发育特征,水土流失监测.水土保持规划有着重要的意义. 沟壑密度也称沟谷密度或沟道密度,指单位面积内沟壑的总长度.以公里/平方公里为单位.数学表达为 Ds为沟壑密度﹔>L为研究区域内的沟壑总长度(单位:公里);A为特定研究区域的面积(单位:平方公里). 二.目

关键词: Arcpy,ArcMap Toolbox, Coordinates of Polygon 需求分析: Arcmap中添加有三个图层:图层“GY乡级河道”(面图层)是全县13个乡镇全部的乡级河道,共计249条:图层“乡级河道polygon”(面图层)是其中某一个乡镇的乡级河道,有10条乡级河道,现在以此乡镇的10条乡级河道为例,求取这10个面要素的起终点坐标:图层“乡级河道polyline”是已经求好的这10条乡级河道的中心线,求取面图层的中心线的过程这里不再赘述. 因为面图层计算坐标时

# -*- coding: utf-8 -*- # --------------------------------------------------------------------------- # JSONToSHP.py # Created on: 2014-09-01 22:00:00.00000 # # Description: # --------------------------------------------------------------------------

上一篇介绍了shp文件的创建,接下来介绍dbf的创建. 推荐结合读取dbf的博客一起看! 推荐结合读取dbf的博客一起看! 推荐结合读取dbf的博客一起看! 1.Dbf头文件的创建 Dbf头文件的结构如下: 记录项数组说明: 字段类型说明: 关于每项的具体含义参照读取dbf文件的解释,这里重点解释几项: HeaderByteNum指dbf头文件的字节数,数值不用除于2,具体为:从version到Reserved2(共32) + n个字段 * 每一个字段长度 32 + terminator. Re

说明:首先要将sde要素类发布成对应的要素服务,通过对要素服务的操作,实现数据在sde要素类中的增删 //向服务器发出请求 public string getPostData(string postString, string url) { string result = ""; try { byte[] postData = Encoding.UTF8.GetBytes(postString); WebClient webClient = new WebClient(); strin

一.问题 现场有一张大数据量的分区表,数据量在10G以上.因某种原因需要删除其中的某些字段.如果直接用 alter table1 drop (column1,column2); 或者alter table1 drop column column1;和alter table1 drop column column2;  的话,需要执行很长时间,这期间该表被锁,会影响到其它应用. 二.解决方法 使用set unused,等系统空闲时再drop unused. 1. alter table table

最近整理空间数据库时,用sde比较多,发现在编辑sde中的数据时总是出现数据被锁或者是被其他应用程序占用.用了很多方法处理,但不是每个方法都实用.下面讲的是我在删除shp或者给shp增加字段时所遇到的三种处理成功的办法.我们的sde都在服务器上. 暂时已经没有办法将提示的错误重现了,下面只讲解决办法. (1)重启sde所在服务器.如果sde中的shp数据在server中发不成服务了,则需要在server停用相关服务.然后应该就可以操作sde中的shp了. (2)上面的方法不行,尝试在plsql(

本节将介绍如何利用python完成对shp的基本操作 1.读取shp四至 import shapefile sf = shapefile.Reader(r"E:\shp\1.shp") #读取shp四至 min_x, min_y, max_x, max_y = sf.bbox #读取每个图斑四至 shapes = sf.shapes() arr = [] for i in range(0, len(shapes)): arr.append(shapes[i].bbox) 利用GDAL

CSSDesk body { background-color: #2574b0; } /*! zybuluo */ article,aside,details,figcaption,figure,footer,header,hgroup,nav,section{display:block}audio,canvas,video{display:inline-block;*display:inline;*zoom:1}audio:not([controls]){display:none}html{

原理: 利用图斑最小外包矩形的左上角坐标(数学坐标)Y坐标将序.X坐标升序的方式获取自上而下,从左到右的要素记录排序,然后遍历编号. "!shape.extent.xmin!"计算要素最小外包矩形的左上角X坐标: "!shape.extent.ymax!" 计算要素最小外包矩形的左上角Y坐标: 字段计算器中可做Python脚本表达式用,Arcpy中其类型为Python_9.3,在“解析ArcGis字段计算器”相关博文中,已对上面代码信息进行过说明. 代码如下: #要

转换工具 1.Excel # Process: Excel 转表 arcpy.ExcelToTable_conversion("", 输出表, "") # Process: 表转 Excel arcpy.TableToExcel_conversion("", 输出_Excel_文件, "NAME", "CODE") 2.JSON # Process: JSON 转要素 arcpy.JSONToFeature