![](https://img-blog.csdnimg.cn/20210829092513582.jpg?x-oss-process=image/resize,m_fixed,h_224,w_224)
Python(arcpy) 学习日记
记录学习python的心得
天外来客Jonty
超越过去,一切自今日始
展开
-
Python(arcpy) 批量shp转raster
# -*- coding: utf-8 -*-'''批量实现shp转raster2021/12/24'''import arcpyimport osimport re#检查Arcgis的权限arcpy.CheckOutExtension("spatial")arcpy.gp.overwriteOutput = 1#设置了输出范围及输出坐标系,这三行也可注释掉arcpy.env.overwriteOutput = Truearcpy.env.extent = arcpy.Ext原创 2021-12-24 14:31:35 · 1222 阅读 · 2 评论 -
Python(arcpy) 批量计算流域面积
#-*- coding: UTF-8 -*-'''计算子流域面积'''import arcpyimport globimport osimport xlrdimport rearcpy.env.workspace = "F:WATA_polygon_shp\\"outpath = "F:WATA_shp_with_area_field\\"shps = arcpy.ListFiles("*.shp")for shp in shps: oldname = [substr.sta原创 2021-11-29 22:20:56 · 1412 阅读 · 0 评论 -
Python(arcpy) 对多个流域、多年图像上进行批量掩膜裁减
进行批量掩模裁剪的准备工作:1.将边界文件(线元素)转为掩膜文件(面元素)2.将遥感图像和掩膜文件存储在同一文件夹下(路径一定要一致!!)3.将需要进行操作的流域名汇总在EXEL中第一列4.检查Arcgis的Arcgis administrator 中各个工具箱的授权情况,如果过期要重新授权好,接下来是第一步:#-*- coding: UTF-8 -*- #识别中文'''这篇...原创 2019-03-27 10:35:54 · 246 阅读 · 1 评论 -
Python(arcpy) 根据站点经纬度(坐标)批量提取对应格点值
任务描述:#-*- coding: UTF-8 -*-'''这篇代码是用于根据站点的shp文件批量提取降雨图中该点所在像元的值并存为新的点元素文件'''import arcpyfrom arcpy import envfrom arcpy.sa import *import osarcpy.CheckOutExtension("Spatial") #仅在使用浮动板许可情况下需要 CheckOutExtension获取许可#先遍历点元素文件,再遍历栅格文件,做提取cFi原创 2021-04-07 16:03:07 · 6499 阅读 · 19 评论 -
Python(arcpy) 批量将线元素转化为面元素
我所参考的原文链接如下:https://blog.csdn.net/lucky51222/article/details/49699847 《Python实现PointToRaster》#-*- coding: UTF-8 -*- #识别中文'''这篇代码是在尝试将线文件通过调用FeatureToPolygon工具转化成面文件以便于后期的掩模裁剪时使用'''impor...原创 2019-03-26 22:11:31 · 1515 阅读 · 0 评论 -
Python(arcpy) 从EXEL中读取文件名,并通过模糊搜索实现批量转移指定文件到指定文件夹
背景介绍:想将G盘All_catchment文件夹中的部分文件转移到F盘中,筛选条件是Desktop / name.xlsx,即从EXEL中读取文件名作为筛选条件软件:Sublime Text#导入模块#!/usr/bin/env python # -*- coding: utf-8 -*- #识别中文import osimport shutilimport x...原创 2019-03-26 17:18:17 · 2388 阅读 · 25 评论 -
Python(arcpy) 批量获取并输出栅格均值
首先声明,这篇代码是改自他处,链接是:https://blog.csdn.net/liyanzhong/article/details/44935885 《python获取栅格点和面值》直接贴修改后的代码:import arcpyOutputFile = open('statisticResult.txt', 'w') arcpy.env.workspace = ('G:\Glob...原创 2019-03-20 08:34:20 · 6273 阅读 · 8 评论 -
Python(arcpy) 重采样
任务描述:有高精度的数据,想要重采样降低精度以便跟其它数据匹配import arcpyfrom arcpy import env #导入 env 类(env 类包含所有地理处理环境)#定义工作环境arcpy.env.workspace = "F:\\" #重采样时的捕捉栅格(定义环境),这一句可要可不要,取决于你有没有需要对齐的目标栅格 arcpy.env.snapRaster = "F:\\sample_std.tif" #输入栅格原创 2021-04-05 23:49:55 · 2484 阅读 · 0 评论 -
Python(arcpy) NameError: name ‘RasterCalculator‘ is not defined
想用栅格计算器进行批量处理数据,结果提示未定义该函数报错原因:栅格计算器I hope this has given you some insight原创 2021-04-05 23:08:27 · 728 阅读 · 0 评论 -
Python(arcpy) 批处理多波段HDF文件为TIF文件
研究需要批量处理多年雪盖数据,NASA给出的数据源是包含两个波段(0和1)的Hdf文件,需要转换为Tif文件才能方便在ArcGIS中提取更多信息,为方便读者查看,现将数据源及其路径展示如下:#-*- coding: UTF-8 -*- #识别中文'''这篇代码是...原创 2019-12-09 17:27:33 · 6939 阅读 · 5 评论 -
Python(arcpy) 转换投影坐标系
任务描述:转换投影坐标系为WGC_1984# -*- coding:utf-8 -*-'''用于全球分布图重投影'''import arcpy#工作区arcpy.env.workspace = r"F:\\MODIS_NDVI_Monthly_1km_v006\\tif"#输入栅格in_raster = 'F:\\MODIS_NDVI_Monthly_1km_v006\\tif\\201001.tif'#输出栅格out_raster = 'F:\\MODIS_NDVI_Mont原创 2021-04-05 20:50:42 · 2883 阅读 · 3 评论 -
MODIS数据下载及图像处理教程
任务描述:如题,以2010年月尺度1km的MODIS的植被覆盖度(NDVI)数据为例第一步 :获得MODIS数据下载链接Earthdata Search(下载地址)可以选择自己想要的时间空间范围,NASA官网会自动生成下载链接,点击保存得到txt文件第二步:下载数据(使用MATLAB)那么如何根据Txt文件中的url链接下载数据呢?%% 代码用于下载MODIS NDVI月尺度数据,全部存在默认文件夹中,需要将系统默认浏览器设置为Chrome,存储路径在浏览器中修改clcclearclos原创 2021-04-05 14:44:49 · 12422 阅读 · 7 评论 -
Python(arcpy) 批量裁减TIF图像
任务描述:已有流域边界的面文件,从长时间序列的全球尺度栅格图像裁减出流域的部分效果如下:import arcpyimport globimport osimport xlrd#检查Arcgis的权限arcpy.CheckOutExtension("spatial")arcpy.gp.overwriteOutput = 1##################################裁减径流图##################################读取掩膜文件(矢原创 2021-03-27 11:57:07 · 4121 阅读 · 3 评论 -
Python 对ASCII文件(及TIF文件)进行批量掩模裁剪
首先声明,这篇文章是关于网上诸多方法的测试和修正,所以是基于大佬的文章之上的,参考文章的链接如下: https://blog.csdn.net/gisboygogogo/article/details/75195760利用ArcGIS Python批量拼接裁剪遥感影像(arcpy batch processing) https://blog.csdn.net/dou3516/article...原创 2019-03-13 23:50:56 · 4558 阅读 · 5 评论 -
Python(arcpy) 批量读取点元素文件(shp)属性表的字段
任务描述:我有3年内某站点的逐日降雨数据,shp文件的属性表中包含经纬度和降雨量数据(就是下图右边的Lat、LonRASTERBALU),想把他们全部提取到一个TXT中在网上找了好久能读取点元素文件属性的代码,眼花缭乱,最后终于找到个能用的:python/arcpy提取shp文件属性表中的字段值我的代码是在他的基础上完成的,主要加了批量操作,感谢这位博主!#-*- coding:utf-8 -*-'''这篇代码是用于批量读取点元素的属性值'''#导入模块,使用该模块里定义的类、方法或者原创 2020-10-08 11:20:23 · 5955 阅读 · 3 评论 -
Python(arcpy) 提取包含目标流域的矩形区域
任务描述:有一张全国的DEM图以及一个流域的边界文件,想要从全国DEM图中裁出一个包含该流域的矩形区域,以便将流域以外的邻近区域也考虑进来# -*- coding: utf-8 -*-"""这篇代码是用于提取比目标流域大一点(包含目标流域)的矩形区域"""import sysimport osimport shapefileimport arcpyfrom arcpy import env #导入 env 类(env 类包含所有地理处理环境)from arcpy.sa impo原创 2020-11-11 10:36:49 · 1169 阅读 · 2 评论 -
Sublime 如何手动添加python库(以shapefile.py为例)
Sublime挺好用的,但很多时候需要用到的库不是自带的,试过很多办法,最后发现还是手动安装来的快,以下以shapefile.py为例讲解:第一步:找到目标库的网上资源我在Github上找到了shapefile.py第二步:确定你的Sublime的系统变量import sysprint(sys.path)运行以上代码,会出如下结果:可以看到大部分包都是存在F:\Python27\ArcGIS10.2\lib路径下,所以把下好的shapefile.py复制到该路径即可...原创 2020-11-11 10:47:23 · 2214 阅读 · 0 评论