R读取Sentinel-1 TIFF文件时空间属性识别异常问题

问题原因 #

Sentinel-1 GRD格式的TIFF文件常使用**模型变换标签（ModelTransformationTag）**存储空间信息，而非标准GeoTIFF的地理变换标签。R中的raster/terra/stars包默认读取逻辑可能未解析该标签，导致投影、范围、分辨率识别错误，而ArcGIS会自动处理这类特殊元数据。

解决方案 #

方法1：强制GDAL解析模型变换标签（推荐） #

使用terra包的rast()函数，配合GDAL选项强制读取模型变换信息：

library(terra)
# 读取Sentinel-1 TIFF，强制GDAL使用模型变换
s1_rast <- rast("s1a-iw-grd-vh-20230208t052407-20230208t052432-047140-05a7dd-002.tiff", 
                options = c("GTIFF_FORCE_RGBA=NO", "PROJ_USE_GRIDS=NO"))

若仍无效，尝试指定GTiff驱动：

s1_rast <- rast("s1a-iw-grd-vh-20230208t052407-20230208t052432-047140-05a7dd-002.tiff", 
                driver = "GTiff")

方法2：手动提取并赋值空间属性 #

如果方法1失效，通过gdalUtils提取元数据后手动设置空间属性：

1. 安装并加载依赖包：

install.packages("gdalUtils")
library(gdalUtils)
library(stringr)

2. 编写批量读取函数：

read_sentinel1 <- function(file_path) {
  # 获取GDAL元数据
  gdal_meta <- gdalinfo(file_path, raw_output = TRUE)
  
  # 提取EPSG投影
  epsg_match <- str_extract(grep("EPSG", gdal_meta, value = TRUE), "\\d+")
  proj_crs <- paste0("EPSG:", epsg_match)
  
  # 提取地理范围
  ul_coords <- str_extract_all(grep("Upper Left", gdal_meta, value = TRUE), "-?\\d+\\.\\d+")[[1]]
  lr_coords <- str_extract_all(grep("Lower Right", gdal_meta, value = TRUE), "-?\\d+\\.\\d+")[[1]]
  raster_ext <- c(as.numeric(ul_coords[1]), as.numeric(lr_coords[1]),
                  as.numeric(lr_coords[2]), as.numeric(ul_coords[2]))
  
  # 读取栅格并设置空间属性
  s1_rast <- rast(file_path)
  crs(s1_rast) <- proj_crs
  ext(s1_rast) <- raster_ext
  res(s1_rast) <- c(10, 10) # Sentinel-1 GRD默认分辨率为10m
  
  return(s1_rast)
}

3. 批量处理文件：

# 获取指定目录下所有TIFF文件路径
tiff_files <- list.files(path = "./your_tiff_directory", pattern = "\\.tiff$", full.names = TRUE)
# 批量读取生成栅格列表
s1_raster_list <- lapply(tiff_files, read_sentinel1)

验证结果 #

读取完成后，用以下命令验证空间属性：

# 查看投影信息
crs(s1_rast)
# 查看地理范围
ext(s1_rast)
# 可视化确认位置
mapview::mapview(s1_rast)

内容的提问来源于stack exchange，提问作者jjunju