使用背景

上一階段完成了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进行按掩膜提取的批量操作的全部內容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網站內容還不錯，歡迎將生活随笔推薦給好友。