使用背景

上一阶段完成了nc转tiff,这一阶段就要开始准备裁剪了,我在ArcGIS上用模型构建器创建了操作模型,但对运行起来的速度不是很满意,又来寻求python的帮助啦
数据:大量的tif数据
参考文献:
《python|gdal实现按掩膜提取影像/裁剪影像》
《python gdal使用shp裁剪遥感影像》

1. 参考文献中的重点

观察了所有的参考文献,我发现按掩膜提取最重要的是这条语句:

gdal.Warp(outputImage, dataset, cutlineDSName=shp, cropToCutline=True) # 按掩膜提取

所以咱们的参数就跟着这个warp函数进行设置就好了

2. 修改

修改后的代码为:

#按掩膜提取import numpy as np
import netCDF4 as nc
from osgeo import gdal,osr,ogr
import os
import globshp = r"D:\DATA\TIFF\mask\poly.shp" # 圈选范围的路径
Input_folder=r"D:\DATA\nc\REA\2016\2016" #要裁剪的tif文件所在的文件夹
data_list = glob.glob(Input_folder + '/*.tif') #读取文件
for i in range(len(data_list)):   data = data_list[i]inputImage=data    # 遥感影像的路径dataset = gdal.Open(inputImage) # 打开遥感影像num=i+1 #为了方便写入文件的名字outputImage =r"D:\DATA\nc\REA\mask\2016-"+str(num)+".tif" # 按照圈选范围提取出的影像所存放的路径gdal.Warp(outputImage, dataset, cutlineDSName=shp, cropToCutline=True) # 按掩膜提取print(data + '-----转tif成功')

365张影像不到1分钟就处理完成了,速度真的是很快
我看到也有很多用ArcPy做的,也很不错

下一次,就是要用python进行重采样了

用python进行按掩膜提取的批量操作相关推荐

  1. MODIS数据下载,HEG批量拼接处理,arcmap批量掩膜提取,python栅格计算

    前言: 2022.12.30修改:补Python编程基础的过程中发现使用Print也可以向文件里写东西,查看写过的代码,然后发现HEG批量拼接处理代码中有的地方在复制粘贴至博客的时候出现了缩进错误,评 ...

  2. 【ArcGIS风暴】ArcGIS中国地表覆盖数据GlobeLand30预处理(批量投影、拼接、掩膜提取)附成品下载

    结果示意图: GlobeLand30是30米空间分辨率全球地表覆盖数据,目前可供下载的有3年的数据:2000-2010-2020.本文主要讲解在ArcGIS10.6平台下进行GlobeLand30的预 ...

  3. ArcGIS教程:按掩膜提取 (空间分析)

    摘要 提取所定义掩膜区域内的栅格像元. 插图 用法 · 按掩膜提取工具与设置掩膜环境大体相同,但二者的区别在于输入掩膜只适用于即时实例,而在环境中设置的掩膜可在发生更改或被禁用之前应用于所有工具. · ...

  4. python中需要用到的英文词汇-使用python从语料库中提取最常用的单词

    也许这是一个愚蠢的问题,但是我在使用Python从语料库中提取十个最常见的单词时遇到了问题.这就是到目前为止. (顺便说一句,我与NLTK一起阅读一个带有两个子类别的语料库,每个子类别有10个.txt ...

  5. python 从字符串中提取数字 re.findall()

    以前老用(.*?)提取数字,今天发现不对了,比如一行数字为: 0 0.248438 0.255556 0.128125 0.194444 用: re.findall('(.*?) (.*?) (.*? ...

  6. 使用OpenCV和Python从图像中提取形状

    Welcome to the first post in this series of blogs on extracting features from images using OpenCV an ...

  7. python实现视频关键帧提取(基于帧间差分)

    python实现视频关键帧提取(基于帧间差分) 在很多场景下,我们不想或者不能处理视频的每一帧图片,这时我们希望能够从视频中提取出一些重要的帧进行处理,这个过程我们称为视频关键帧提取. 关键帧提取算法 ...

  8. python从字符串中提取数字并转换为相应数据类型_python从PDF中提取数据的示例

    01 前言 数据是数据科学中任何分析的关键,大多数分析中最常用的数据集类型是存储在逗号分隔值(csv)表中的干净数据.然而,由于可移植文档格式(pdf)文件是最常用的文件格式之一,因此每个数据科学家都 ...

  9. python关键词提取源码,python实现textrank关键词提取

    python实现textrank关键词提取 来源:中文源码网    浏览: 次    日期:2018年9月2日 [下载文档:  python实现textrank关键词提取.txt ] (友情提示:右键 ...

最新文章

  1. 深入理解Linux中的文件权限
  2. HTML --- 盒模型
  3. stdthread(6)并发lockGuard
  4. 理解Go语言中的方法和接收者
  5. 配置 BeautifulSoup
  6. VGG Loss的Pytorch实现
  7. java swing如何设置jtextarea对齐方式_【爵士钢琴】一次搞懂爵士经典Swing节奏!
  8. 巴菲特如何滚雪球的?
  9. poj 3709 斜率优化
  10. js 骂人不带脏字 (!(~+[]) + {})[--[~+““][+[]] * [~+[]] + ~~!+[]] + ({} + [])[[~!+[]] * ~+[]] 图解
  11. 小记--解决方案: rsyslogd cpu占用率高
  12. 配置 已完成 请勿关闭计算机,win7系统关机提示“配置Windows Update已完成30%请勿关闭计算机”的解决方法...
  13. windows的由来与详细介绍
  14. PMP : PMP备考心得 (8)
  15. 数字图像频谱的中心化
  16. history of program atan2(y,x)和pow(x,y)
  17. Cocos Creator用cc.loader加载预制体资源和删除预制体资源
  18. PostgreSQL TID及tuple slot
  19. java用模板导出数据表格
  20. 史上最全 XMind 8 快捷键大全

热门文章

  1. 大一寒假周总结(一)
  2. 高速数据采集卡的数据参数表征
  3. JavaScript---常用的鼠标事件mouseover 和mouseenter的区别
  4. Cocos2dx 3.10音频无法正常播放问题的解决之路
  5. 组态王或SCADA与西门子200Smart通讯配置流程
  6. 100个IT人成为大牛必备的学习经验分享
  7. Android开发实验-- 电话拨号器实验报告
  8. C++学习(一五零)qt的公有类、私有类、Q_Q、Q_D、二进制兼容
  9. iOS巅峰之开发过程中的小知识点大全
  10. 嵌入式linux轻量级sshd服务工具Dropbear移植