![](https://img-blog.csdnimg.cn/e095013741684d93b78d58cdc2858918.png?x-oss-process=image/resize,m_fixed,h_224,w_224)
数据处理 代码
数据处理Arcgis Python
GeoLab 219
这个作者很懒,什么都没留下…
展开
-
CLM模型替换大气强迫数据 GSWP3---CMFD
CLM模型替换默认大气强迫数据 0.5°×0.5°GSWP3数据替换为0.1°×0.1°CMFD(中国区域地面气象要素驱动数据集)原创 2023-05-22 16:07:50 · 1339 阅读 · 11 评论 -
经验分享 | SEN+Mk趋势分析(matlab代码分享)
SEN+Mk趋势分析(matlab代码分享)原创 2022-06-10 17:23:46 · 14667 阅读 · 4 评论 -
代码分享 | Python对图像行列号、地理坐标和投影坐标进行转换
01代码描述:利用python对进行图像的行列号、地理坐标和投影坐标三者之间的互相转换工具:python,GDAL库操作:改变输入数据以及想要进行转换的坐标即可例子:下述代码是输入投影坐标为Lambert的栅格数据,读取数据投影进行坐标转换。小Tip:当代码出现下列错误时:PROJ:proj_create_from_database:Cannot find proj.db解决方法:在代码开头设置环境变量,添加你的proj.db所在的位置例如:os.environ['P原创 2022-05-14 15:07:53 · 2168 阅读 · 0 评论 -
代码分享 | Arcgis批量栅格值提取至点(自定义工具箱或代码)
任务举例:文件夹有20年某区域影像,需要提取该区域某点或多点值。数据准备:1.矢量点文件 2.需要处理的栅格影像(存在同一文件夹) 3.存储影像文件夹方法1工具准备:ARCGIS操作:自定义脚本工具1.我的工具箱添加脚本文件(后缀为.py文件)。2.注意存储相对路径记得打钩。3.导入脚本文件。4.脚本文件中的代码如下。# coding=gbkimport arcpyfrom arcpy import env #定义envfrom arcpy.sa impo..原创 2022-05-08 18:19:17 · 6734 阅读 · 5 评论 -
代码分享 | tiff,nc,hdf,mat,grib,raw遥感数据python读取
1. tiff数据read TIFF image#读取tiff数据需要GDAL库from osgeo import gdaldef Read_tiff(file_path): # img_file_path: 栅格数据路径 # return: 返回投影,几何信息,和转换后的数组dataset=gdal.Open(img_file_path)#读取栅格数据 # 判断是否读取到数据 if dataset is None: pri...原创 2022-04-22 20:04:10 · 1424 阅读 · 3 评论 -
基于第三方工具批量下载数据---以下载GLASS数据为例
GLASS数据处理(1)—— 数据下载1.根据下图确定研究区域2.下载网站GLASS数据下载中心(http://www.glass.umd.edu/Download.html)以GLASS FVC数据为例意思为分辨率为500米的GLASS数据–时间跨度00-18年–每一年是8天一景3.由于数据量太多,撰写代码进行下载...原创 2021-10-25 19:13:17 · 2568 阅读 · 5 评论 -
代码分享 Python | 批量nc文件转tif
01准备工作:查看nc文件属性等。工具:Panoply、Matlab等软件。操作:1.使用Panoply 软件。2.使用Matlab软件。即可查看nc文件内各种属性;如果想单独查看变量:(以经度为例)02任务举例:降水nc数据转tif,坐标系为WGS 84。工具准备:Python GDAL和netCDF4环境。操作:替换代码中的属性变量名和路径即可。任务举例:降水nc数据转tif,坐标系为WGS 84。工具准备:Python GDAL和netC.原创 2022-04-14 17:45:39 · 3600 阅读 · 8 评论 -
批量裁剪栅格文件(直接运行或自定义脚本工具)
批量裁剪栅格文件(直接运行或自定义脚本工具)利用数据:1.矢量文件 2.要裁剪呈矢量文件范围的栅格文件利用工具:arcgis中spatial analyst工具–提取分析–按掩膜提取方法一:arcpy代码# -*- coding:utf-8 -*-import arcpyfrom arcpy import env #定义envfrom arcpy.sa import *arcpy.env.workspace = "D:/GLASS/实验/" #定义工作空间及数据路径rasters =原创 2021-07-07 14:30:26 · 952 阅读 · 3 评论 -
代码分享 Python |批量删除栅格影像异常值
任务举例:下载的NPP影响最大值为63355,通过查看栅格数据直方图要删除异常值。数据准备:文件夹(工作空间)-含有需要批量处理的tif影像工具准备:Python arcpy 环境操作:替换代码中的路径即可# -*- coding: UTF-8 -*-import arcpyimport osimport globimport arcpyfrom arcpy.sa import *arcpy.CheckOutExtension("ImageAnalyst") # 检查许可a原创 2022-04-09 20:51:12 · 2060 阅读 · 0 评论 -
代码分享 Python | 批量计算栅格影像最大值,最小值,总和,均值等
01任务举例:文件夹内有12个月遥感影像求年最大值,累加值等(生成的是一幅栅格影像)数据准备:文件夹(工作空间)-含有需要批量处理的tif影像工具准备:Python arcpy 环境操作:替换代码中的路径即可# -*- coding: UTF-8 -*-import arcpyfrom arcpy import envfrom arcpy.sa import *# Set environment settingsenv.workspace = "E:/NPPresult.原创 2022-04-09 15:11:44 · 5943 阅读 · 3 评论 -
python安装包 whl文件rasterio包安装
whl 文件获取添加链接描述https://www.lfd.uci.edu/~gohlke/pythonlibs/出现错误,如下原因:GDAL包版本和rasterio包不匹配最后试出成功的匹配,我的python版本3.6![在这里插入图片描述](https://img-blog.csdnimg.cn/5637be6b83e14d0495deb680f42e1bd5.png...原创 2022-03-28 21:19:23 · 1144 阅读 · 0 评论 -
Arcgis/python批量栅格计算(乘法*缩放因子)
本文所用数据为一文件夹所有tif影像,进行乘0.01的缩放因子,处理后影像放入新文件夹中1.如果Python装好了arcgis运行环境,可以直接;python运行。# coding:utf-8import osimport arcpyfrom arcpy import envfrom arcpy.sa import *import stringarcpy.CheckOutExtension("spatial")arcpy.gp.overwriteOutput=1arcpy.env.wo原创 2022-03-03 09:40:05 · 3695 阅读 · 2 评论 -
arcpy/python循环批量裁剪tiff栅格文件
**arcgis利用一个矢量批量裁剪多个栅格,可参考**批量裁剪栅格文件2.arcgis利用多个矢量要裁剪多个栅格文件,裁剪后自动生成文件夹保存。例如:有全国每省的矢量和全国2000-2020年的植被覆盖图,现在要分省统计进行分析。不同省shp文件可以裁剪全国影像后,裁剪好的影像在路径中用省市名字保存。# -*- coding: UTF-8 -*-import arcpyimport osimport globimport arcpyfrom arcpy.sa import *a原创 2021-10-28 09:53:09 · 1711 阅读 · 1 评论