数据范围转换回原CRS后未恢复的原因及等效处理方法咨询

嘿，我来帮你理清楚这个问题——这绝对不是bug，是投影转换的固有特性导致的，我来给你拆解明白，再给你正确的操作姿势：

为什么转回去坐标不一样了？ #

首先你得明白：不同坐标系之间的转换不是可逆的完美映射，尤其是你用的epsg:32198是兰伯特正形圆锥投影，epsg:4326是WGS84地理坐标系（经纬度），这俩的转换是基于椭球面的复杂计算，不是简单的线性变换。

更关键的是，你在转换整个SpatExtent的时候，其实是在处理矩形范围的投影变形：

1. 原投影坐标系里的矩形范围，转换到地理坐标系（经纬度）时，原本的直边会变成曲边，但SpatExtent只能用矩形来表示，所以terra会自动生成一个能包裹住所有曲边的最小外接矩形——这时候你的lulc_4326已经不是原范围的精确对应了，它是一个更大的近似矩形。
2. 再把这个近似矩形转回到原投影坐标系，自然就和原始范围完全不一样了，毕竟两次转换都有近似和计算精度的影响。

怎么实现你的需求（只改ymax，其他坐标和原范围一致）？ #

你的核心需求应该是：保持原范围的xmin、xmax、ymin不变，只把ymax改成对应经纬度下50.5°N的投影坐标值。根本不需要来回转换整个extent，直接针对单个坐标点转换就行，步骤如下：

library(terra)
# 原始范围
lulc.e <- ext(-830370, 787470, 117960, 2091060)

# 把经纬度50.5°N转换为epsg:32198下的y坐标
# 随便选一个x值就行，兰伯特投影的y坐标只和纬度相关
target_point <- vect(cbind(0, 50.5), crs = "epsg:4326")
projected_point <- project(target_point, to = "epsg:32198")
new_ymax <- y(projected_point)

# 直接修改原始范围的ymax，其他值保持不变
lulc_mer <- lulc.e
lulc_mer[4] <- new_ymax

这样得到的lulc_mer就完全符合你的要求：除了ymax是你想要的对应50.5°N的投影值，其他坐标和原始范围完全一致，不会出现转换带来的偏差。

内容来源于stack exchange