使用 Shapefile 提取多个 .tiff 文件

发布于 2025-01-16 21:06:42 字数 737 浏览 3 评论 0原文

我有一个名为 Mar22_073685 的文件夹，其中包含 1372 tiff 图像。我需要使用另一个名为“Export_Output_2_Project.shp”的文件夹中的 shapefile 裁剪这些 tiff 文件。我在循环的帮助下裁剪这些 tiff 文件，但是当我打开输出文件时，这些文件中没有数据。它是空白的。此外，文件不会以与输入相同的名称保存。

我试过：

require(rgeos)
require(sp)
require(sf)
require(maptools)
require(raster)
require(rgdal)
Shape <- readOGR("F:/NSSS_2019/Export_Output_2_Project.shp")
plot(Shape)
flname1 <- list.files(path = 'F:/FOG_ISLANDS/DEC_2015/Mar22_073685/', pattern = '.tif',full.names=TRUE)
nfile1 <- length(flname1)
for (i in 1:nfile1) {
r1 <- raster(flname1[i])
cropimage <- mask(r1, Shape)
writeRaster(cropimage, paste0('F:/FOG_ISLANDS/DEC_2015/test/', flname1[i], '.tiff', sep = ''), format = 'GTiff', overwrite = T) 
}

原文

I have a folder named Mar22_073685 with 1372 tiff image . I need to crop these tiff files using a shapefile in another folder named 'Export_Output_2_Project.shp'. I am cropping these tiff files with the help of loop but when I am opening the output files there is no data in those files. It is blank. Also, the files are not saved with the same name as the input.

I tried:

require(rgeos)
require(sp)
require(sf)
require(maptools)
require(raster)
require(rgdal)
Shape <- readOGR("F:/NSSS_2019/Export_Output_2_Project.shp")
plot(Shape)
flname1 <- list.files(path = 'F:/FOG_ISLANDS/DEC_2015/Mar22_073685/', pattern = '.tif',full.names=TRUE)
nfile1 <- length(flname1)
for (i in 1:nfile1) {
r1 <- raster(flname1[i])
cropimage <- mask(r1, Shape)
writeRaster(cropimage, paste0('F:/FOG_ISLANDS/DEC_2015/test/', flname1[i], '.tiff', sep = ''), format = 'GTiff', overwrite = T) 
}

分享到QQ

分享到微博