如何解决不同年份Sentinel-1同边界提取影像在Rasterio中的错位问题

问题原因 #

不同年份的Sentinel-1升降轨影像，因成像轨道方向、视角差异，即使覆盖同一区域，原始影像的**坐标系、栅格对齐参数（分辨率、左上角地理坐标）**也可能不一致。直接读取转NumPy数组会丢失地理信息，导致视觉上的错位；同年份数据因轨道参数更统一，所以无此问题。

解决步骤 #

1. 先确认影像的地理参数差异 #

用Rasterio检查两张影像的坐标系（CRS）和地理变换（Transform），定位错位根源：

import rasterio as rs

# 替换为你的两个文件路径
file_asc = "sentinel1_asc_2022.tif"
file_desc = "sentinel1_desc_2018.tif"

with rs.open(file_asc) as src_asc, rs.open(file_desc) as src_desc:
    print("升降轨坐标系:", src_asc.crs, src_desc.crs)
    print("升降轨地理变换:", src_asc.transform, src_desc.transform)
    print("升降轨分辨率:", src_asc.res, src_desc.res)

如果输出的CRS、Transform或分辨率不同，就是导致错位的直接原因。

2. 重投影配准到同一基准 #

选择其中一张影像作为参考（比如2022年升轨影像），将另一张错位影像重投影到参考影像的坐标系、分辨率和地理范围上：

from rasterio.warp import reproject, Resampling

# 参考影像（无错位的基准）
with rs.open(file_asc) as ref_src:
    ref_meta = ref_src.meta.copy()
    ref_crs = ref_src.crs
    ref_transform = ref_src.transform
    ref_height, ref_width = ref_src.shape

# 处理待配准的降轨影像
with rs.open(file_desc) as mis_src:
    mis_data = mis_src.read()
    # 执行重投影配准
    aligned_data, aligned_transform = reproject(
        source=mis_data,
        dst_crs=ref_crs,
        dst_transform=ref_transform,
        src_crs=mis_src.crs,
        src_transform=mis_src.transform,
        resampling=Resampling.bilinear  # 根据需求选择：最近邻、双线性等
    )
    # 更新元数据为参考影像的参数
    ref_meta.update({
        'height': ref_height,
        'width': ref_width,
        'transform': aligned_transform
    })
    # 保存配准后的影像
    aligned_file = "sentinel1_desc_2018_aligned.tif"
    with rs.open(aligned_file, 'w', **ref_meta) as dst:
        dst.write(aligned_data)

3. 优化GEE导出设置（从源头避免错位） #

在Google Earth Engine导出影像时，统一指定投影、分辨率和裁剪参数，避免自动生成不同的栅格对齐：

// GEE示例代码：导出时指定统一参数
var refImage = sentinel1Asc2022; // 参考影像
Export.image.toDrive({
  image: sentinel1Desc2018,
  description: "sentinel1_desc_2018_aligned",
  region: studyArea,
  crs: refImage.projection(), // 使用参考影像的投影
  crsTransform: refImage.projection().getInfo().transform, // 参考影像的栅格变换
  scale: 10, // 统一分辨率（Sentinel-1默认10米）
  maxPixels: 1e13
});

4. 处理影像时保留地理信息 #

原代码直接将影像转为数组后存为PNG，丢失了地理坐标信息。如果需要保存带地理信息的文件，应使用Rasterio写入TIFF，而非仅存PNG：

# 替换原代码中的保存步骤
with rs.open(file_path) as src:
    meta = src.meta.copy()
    # 处理影像数据...
    img_fixed = img_fixed.transpose(2,0,1)  # 转回Rasterio的通道优先格式
    meta.update(dtype='uint8', count=img_fixed.shape[0])
    with rs.open(png_file_path.replace('.png', '_geo.tif'), 'w', **meta) as dst:
        dst.write(img_fixed)

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