python
木易GIS
GIS遥感交流学习
展开
-
Arcpy给要素重新添加编码序号。ArcGIS的属性表中,由于编辑修改的原因,默认的FID或OID并不连续,经常需要给要素添加连读的序号,可使用Arc Python代码完成。
Arcpy给要素重新添加编码序号。ArcGIS的属性表中,由于编辑修改的原因,默认的FID或OID并不连续,经常需要给要素添加连读的序号,可使用Arc Python代码完成。原创 2023-09-12 11:10:55 · 416 阅读 · 0 评论 -
arppy gis 读取text 并批量添加字段 arcpy.AddField_management
arppy gis 读取text 并批量添加字段 arcpy.AddField_management原创 2023-09-11 09:21:07 · 295 阅读 · 0 评论 -
arcpy gis根据字段批量拆分(批量分析工具-筛选功能)
arcpy gis根据字段批量拆分(批量分析工具-筛选功能)原创 2023-09-11 09:12:50 · 856 阅读 · 0 评论 -
Python几何计算库读取几何,点线面要素游标搜索
Python几何计算库读取几何,点线面要素游标搜索原创 2023-08-17 17:14:10 · 307 阅读 · 0 评论 -
arcpy将xml文件转换为矢量shp格式文件
arcpy将xml文件转换为矢量shp格式文件原创 2023-06-07 10:29:01 · 1310 阅读 · 0 评论 -
arcpy批量合并对应gdb名称的要素,使用Append_management方法
arcpy批量合并对应gdb名称的要素,使用Append_management方法原创 2023-05-31 17:26:31 · 565 阅读 · 0 评论 -
arcpy判断featureclass要素类的几何类型,点线面判断
arcpy判断featureclass要素类的几何类型,点线面判断原创 2023-05-29 10:37:15 · 446 阅读 · 2 评论 -
arcpy提取符合gis条件的面要素的质心,并创建点状要素类,写入某些字段值:
arcpy提取符合gis条件的面要素的质心,并创建点状要素类,写入某些字段值原创 2023-05-29 10:31:28 · 541 阅读 · 0 评论 -
使用python实现ndvi归一化植被指数计算
使用python实现ndvi归一化植被指数计算原创 2023-02-07 08:45:12 · 1156 阅读 · 0 评论 -
arcpy批量修改shp文件的空间参考坐标系Spatial Reference
arcpy批量修改shp文件的空间参考坐标系Spatial Reference原创 2022-12-08 16:12:38 · 1072 阅读 · 0 评论 -
arcpy批量删除DeleteField/添加AddField shp文件字段:
arcpy批量删除DeleteField/添加AddField shp文件字段:原创 2022-12-06 15:41:44 · 1058 阅读 · 0 评论 -
arcPy通过经纬度坐标生成 Polygon的方法
arcPy通过经纬度坐标生成 Polygon的方法原创 2022-12-06 15:37:01 · 1060 阅读 · 1 评论 -
arcPy与添加 shp 文件(显示在 ArcMap),通过Arcpy打开mxd文档并添加数据的方式
arcPy与添加 shp 文件(显示在 ArcMap),通过Arcpy打开mxd文档并添加数据的方式原创 2022-12-06 15:32:47 · 1057 阅读 · 0 评论 -
arcpy gis栅格计算器RasterCalculator的实现
arcpy gis栅格计算器RasterCalculator的实现原创 2022-12-01 08:59:23 · 2518 阅读 · 2 评论 -
arcpy在gis中导入带经纬度坐标的 txt 文件
arcpy在gis中导入带经纬度坐标的 txt 文件原创 2022-12-01 08:48:02 · 1995 阅读 · 0 评论 -
arcpy实现gis同一图层的多个要素合并(Dissolve/Union)方法
arcpy实现gis同一图层的多个要素合并(Dissolve/Union)方法原创 2022-12-01 08:39:36 · 3755 阅读 · 0 评论 -
学习Arcpy,从零开始学gis
+” : 连接字符串(Eg: Part = “D:\Data\polygon” + “.shp”,输出Part:D:\Data\polygon.shp)4.3 列表(list,它是元素的有序集合,可以存放一种数据类型,也可以是多种数据类型(字符串,数字,其他列表,字典,对象等))4.5 元组(它是静态的,即创建后,不能添加也不能删除,也不能更改顺序)pop()删除列表中的元素(默认是最后一个)并返回该元素的值。4.6 字典(类似于列表,但字典是对象的无序集合,不能切片)学习Arcpy,从零开始学gis。原创 2022-12-01 08:30:21 · 267 阅读 · 0 评论 -
arcgis使用arcpy中按字段拆分要素到单独的shpfile文件 很好用,简单易懂
arcgis使用arcpy中拆分获取FeatureClass中各类型地物要素到单独的shpfile文件很好用,简单易懂。arcgis使用arcpy中按字段拆分要素到单独的shpfile文件:原创 2022-09-07 08:38:04 · 1681 阅读 · 2 评论 -
2022武汉大学资源与环境学院gis专业arcpy考试:arcpy代码答案及数据输入输出
2022武汉大学资源与环境学院gis专业arcpy考试:arcpy代码答案及数据输入输出原创 2022-07-21 11:06:52 · 536 阅读 · 0 评论 -
python批量转换:未知类型二进制净生产力(NEP)遥感定量反演产品数据的打开,转换为tif,可用于通用GIS、遥感软件打开。
python批量转换:未知类型二进制净生产力(NEP)遥感定量反演产品数据的打开,转换为tif,可用于通用GIS、遥感软件打开。原创 2022-05-20 08:23:45 · 1373 阅读 · 4 评论 -
python删除指定路径的目录中的空文件夹os.rmdir()
python删除指定路径的目录中的空文件夹os.rmdir():import osdef del_emp_dir(path): for (root, dirs, files) in os.walk(path): for item in dirs: dir = os.path.join(root, item) try: os.rmdir(dir) #os.rmdir() 方法用于删除指定路径的目录。仅原创 2021-11-03 11:48:54 · 420 阅读 · 0 评论 -
arcpy实现LUCC土地利用变化转移矩阵,2期影像或者矢量数据土地利用变化统计,生成统计表。使用“面积制表TabulateArea”轻松实现变化统计。
arcpy实现LUCC土地利用变化转移矩阵,2期影像或者矢量数据土地利用变化统计,生成统计表。使用“面积制表TabulateArea”轻松实现变化统计。coding = 'utf-8' #utf8编码import os #加载os库import arcpy #加载arcpy库from arcpy import env #导入环境变量from arcpy.sa import * #导入Spatial Analyst模块env.workspace = r"D:\Ex1" #设置工作空间为数据文件原创 2021-10-11 09:46:03 · 1297 阅读 · 0 评论 -
GIS遥感方面的数据处理、arcpy及批处理脚本的咨询、交流
应部分博友的需求,欢迎GIS遥感方面的数据处理、arcpy及批处理脚本的咨询、交流。博客中文章也有相关的应用案例,如果有特殊定制化需要的可以私信联系。(1)arcpy功能开发及批处理;(2)arctool box小工具制作;(3)栅格、矢量数据批处理;(4)影像数据下载;(5)python批处理实现...原创 2021-09-14 17:02:48 · 245 阅读 · 1 评论 -
python通过gdal读取、写入存储栅格数据。栅格数据tif转化为array数组、array数组转化为栅格tif的函数方法
python通过gdal读取、写入存储栅格数据。栅格数据tif转化为数组数组、array转化为栅格tif的函数方法:首先必须引用引用gdal的模块包。from osgeo import gdal, gdalconst(1)python读取tif转化为数组的函数 def read_img(self, img_path): """读取遥感数据信息""" dataset = gdal.Open(img_path, gdalconst.GA_ReadOnly)原创 2021-08-16 09:23:49 · 9120 阅读 · 7 评论 -
arcpy GIS提取线段要素的所有端点坐标(startPoint、endPoint):
arcpy提取线段要素的所有端点坐标(startPoint、endPoint):import arcpyfrom arcpy import env# Set environment settingsenv.workspace = "C:/data" inFeatures = "majorrds.shp"outFeatureClass = "c:/output/output.gdb/majorrds_midpt"# Execute FeatureVerticesToPoints# 参数为原创 2021-04-09 10:05:53 · 5621 阅读 · 0 评论 -
arcGIS python进行矢量数据裁剪Clip,要素类裁剪矢量数据。
arcGIS python进行矢量数据裁剪Clip,要素类裁剪矢量数据。#-*- coding:UTF-8 -*-import arcpyfrom arcpy import envenv.workspace = r"C:\Users\HP\Desktop\20210323\FC.gdb"out_clip = r"C:\Users\HP\Desktop\20210323\FC.gdb\result_clip"#矢量数据裁剪Clip FeatureClassFC_shp = r"C:\User原创 2021-03-24 17:19:40 · 2086 阅读 · 0 评论 -
arcGIS python进行栅格裁剪Clip,要素类裁剪栅格数据。
arcGIS python进行栅格裁剪Clip,要素类裁剪栅格数据。#-*- coding:UTF-8 -*-import arcpyfrom arcpy import envenv.workspace = r"C:\Users\HP\Desktop\20210323\tweets.gdb"result_KDE = r"C:\Users\HP\Desktop\20210323\tweets.gdb\result_KDE"#栅格裁剪Clip Raster原创 2021-03-23 16:03:52 · 887 阅读 · 0 评论 -
arcGIS python两个栅格数据,栅格计算器进行栅格代数计算。 信息熵(区位熵)计算
arcGIS python两个栅格数据,栅格计算器进行栅格代数计算。信息熵(区位熵)计算公式:#-*- coding:UTF-8 -*-import arcpyfrom arcpy import envfrom arcpy.sa import *env.workspace = r"C:\Users\HP\Desktop\20210323\tweets.gdb"#CheckOut liscen SpatialAnalysisarcpy.CheckOutExtension("ImageAn原创 2021-03-23 16:01:19 · 3676 阅读 · 0 评论 -
arcGIS python点要素类,求核密度KernelDensity
arcGIS python点要素类,求核密度#-*- coding:UTF-8 -*-import arcpyfrom arcpy import envfrom arcpy.sa import *env.workspace = r"C:\Users\HP\Desktop\20210323\tweets.gdb"#CheckOut liscen SpatialAnalysisarcpy.CheckOutExtension("ImageAnalyst")arcpy.CheckOutExten原创 2021-03-23 15:54:59 · 966 阅读 · 0 评论 -
arcgis中判断多个字段值的逻辑一致性(地理国情监测筛查要素FeatureClass数据):
arcgis中判断多个字段值的逻辑一致性(地理国情监测筛查要素FeatureClass数据):如GB码前四位为4201的,RN值首位必须是G,类型Type必须是国道。(1)新建字段“Hebing”文本型,进行字段计算。实例结果如’4201G国道’。输入语句:left([GB],4) & left([RN],4) &[Type](2)在select by attribute工具,进行选择不符合规范的要素:"Hebing" NOT IN ('4201G国道','4202S省道','原创 2021-03-17 11:32:45 · 1754 阅读 · 0 评论 -
arcgis中判断某一字段值是否在list数组集合里(地理国情监测):
arcgis中判断某一字段值是否在list数组集合里(地理国情监测):在select by attribute工具:"Hebing" NOT IN ('4201G国道','4202S省道','4203X县道','4204Y乡道')原创 2021-03-17 11:26:58 · 496 阅读 · 0 评论 -
arcpy中拆分获取FeatureClass中各类型地物要素到单独的shp中,类似于splitShp的功能(地理国情监测)
arcpy中分离获取FeatureClass中各类型地物要素到单独的shp中,类似于splitShp的功能:for str_Val in set(shp_JH_list): #循环该要素某一字段属性值list集(set去重,把重复的去掉了) shp_JH_list为某一字段的属性值集合(如土地利用一级类) where_Clause ='"CC" = \'%s\''%(str(str_Val)) outfc = "C:/Users/hp\Desktop/新建文件夹/" + str(str_V原创 2021-03-12 08:50:35 · 445 阅读 · 0 评论 -
Arcpy对要素类添加字段、计算统计要素面积Area、删除字段(地理国情监测)
Arcpy对要素类添加字段、计算统计要素面积Area、删除字段:import arcpyfrom arcpy import envenv.workspace = r"E:\001.yangdebin\DOM2020元数据\a.gdb"# 注意!!!!其中fc为FeatureClass的变量值# 添加double类型Area字段arcpy.AddField_management(fc, "Area", "DOUBLE", 20, 10)# 计算每个要素Area面积值,赋值到字段arcpy原创 2021-02-25 10:51:47 · 3456 阅读 · 0 评论 -
python获取List数组中重复元素的个数(arcpy中统计FeatureClass中各类型地物要素的图斑数)(地理国情监测)
python获取List数组中重复元素的个数(arcpy中统计FeatureClass中各类型地物要素的图斑数)for str_Val in set(shp_JH_list): #循环该要素某一字段属性值list集(set去重,把重复的去掉了) print("属性值为" + str_Val + "的字段值为" + shp_JH_list.count(str_Val)) # 打印要素类名称便于查看进度...原创 2021-02-19 15:05:03 · 689 阅读 · 0 评论 -
arcpy获取gdb/Dataset/featureClass中所有要素个数:(地理国情监测)
arcpy获取gdb/Dataset/featureClass中所有要素个数:coding = 'utf-8'import osimport arcpy#Author By AHCH_yangdebingdblist=[] #新建空白列表 file_dir=r"F:\test" #输入待合并GDB数据库所在的文件夹的路径dirlist=os.listdir(file_dir) #获取 D:\\Desktop 下一层级目录的路径for dir in dirlist: #循环dirlist列表原创 2021-02-19 14:57:54 · 3326 阅读 · 0 评论 -
arcpy使用SelectLayerByLocation进行空间位置选择、SelectLayerByAttribute进行属性选择,使用CopyFeatures_management进行选择要素导出。
arcpy使用SelectLayerByLocation进行空间位置选择、SelectLayerByAttribute进行属性选择,使用CopyFeatures_management进行选择要素导出。coding = 'utf-8'import os #加载os库import arcpy #加载arcpy库file_dir_Source ="E:\\001.yangdebin\\001.集水区\\国情监测水网数据\\test.gdb\\CHydDataset\\UV_HYDL" file_dir原创 2021-01-21 11:37:55 · 3162 阅读 · 0 评论 -
win7环境下安装Python时可能要到安装不成功的情况,解决问题步骤
安装Python时可能要到安装不成功的情况,解决问题步骤:(1)先看log日志;有的提示需要安装win7 SP1;(2)微软官方下载,先安装win7 SP1:windows6.1-KB976932-X64.exe;地址:https://www.microsoft.com/en-us/download/details.aspx?id=5842 ...原创 2020-05-28 11:58:00 · 1148 阅读 · 0 评论 -
python读取多层嵌套文件夹中的文件(zip文件嵌套在不同层级的文件夹中)
python读取多层嵌套文件夹中的文件(zip文件嵌套在不同层级的文件夹中):使用递归函数,使用python的os.path.isfile方法判断当前是不是.zip文件,如果不是再用os.listdir方法将子目录循环递归判断。import ospath = 'D:\关键\底图\vec_w'path_read = [] #path_read saves all executa...原创 2020-06-12 15:36:00 · 1218 阅读 · 0 评论 -
arcpy对要素类每个图斑批量出图。(地理国情监测)
每个图斑出一张图:这种情况有Data Driven Pages工具配合最好。打开ArcMap的Customize->Toolbars->Data Driven Pages,设置好图层、名称字段、排序字段、显示范围和比例尺,保存地图。# coding:utf-8import arcpymxdFile =r"F:\GeoData\ChinaArea\ChinaVect...原创 2020-12-25 17:16:00 · 760 阅读 · 0 评论 -
arcgis中shp文件的字符串string型字段转换为日期型Date格式。如“20190426”转为“2019/4/26”...
如果源字符串字段为“20190426”,则代码为:cdate(left( [DLG现势性] ,4) & "-" & mid( [DLG现势性],5,2) & "-" & right( [DLG现势性],2))在属性计算里输入以上代码即可。注意[DLG现势性]为源字符串字段。(1)如果源字符串字段为“2019-04-26”,则代码变为:cdate...原创 2020-12-24 08:47:00 · 3448 阅读 · 1 评论