使用背景

上一階段完成了nc轉(zhuǎn)tiff，這一階段就要開始準(zhǔn)備裁剪了，我在ArcGIS上用模型構(gòu)建器創(chuàng)建了操作模型，但對(duì)運(yùn)行起來(lái)的速度不是很滿意，又來(lái)尋求python的幫助啦
數(shù)據(jù)：大量的tif數(shù)據(jù)
參考文獻(xiàn)：
《python|gdal實(shí)現(xiàn)按掩膜提取影像/裁剪影像》
《python gdal使用shp裁剪遙感影像》

1. 參考文獻(xiàn)中的重點(diǎn)

觀察了所有的參考文獻(xiàn)，我發(fā)現(xiàn)按掩膜提取最重要的是這條語(yǔ)句：

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

所以咱們的參數(shù)就跟著這個(gè)warp函數(shù)進(jìn)行設(shè)置就好了

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 + '-----轉(zhuǎn)tif成功')

365張影像不到1分鐘就處理完成了，速度真的是很快
我看到也有很多用ArcPy做的，也很不錯(cuò)

下一次，就是要用python進(jìn)行重采樣了

總結(jié)

以上是生活随笔為你收集整理的用python进行按掩膜提取的批量操作的全部?jī)?nèi)容，希望文章能夠幫你解決所遇到的問題。

如果覺得生活随笔網(wǎng)站內(nèi)容還不錯(cuò)，歡迎將生活随笔推薦給好友。