R语言krige克里金插值出现CRS对象注释丢失警告问题

该警告是空间包版本兼容问题导致：
新版sf/stars包采用WKT2格式存储CRS信息，内置了PROJ相关注释字段，而你使用的gstat包中krige函数内部仍沿用旧版proj4string接口读取坐标信息，该接口无法处理WKT2中的注释内容，每次转换CRS格式时就会抛出「CRS对象注释丢失」的提示。该警告仅为格式提示，不会影响克里金插值的计算结果。

• 方法1：临时屏蔽警告
  如果不需要处理格式问题，仅想消除警告提示，可直接用suppressWarnings包裹插值代码即可：

krige_result <- suppressWarnings(krige(formula = log(DO_S)~1, GGT, GGTgrid1, model = model_GGT))

• 方法2：统一CRS为PROJ4格式
  将所有空间对象的CRS赋值改为旧版PROJ4字符串，避免生成带注释的WKT2格式CRS：

# 替换原有st_crs赋值代码
proj_str <- "+proj=longlat +datum=WGS84 +no_defs"
st_crs(GGT) <- proj_str
st_crs(GGTgrid1) <- proj_str

• 方法3：升级依赖包
  新版gstat已经修复了该兼容问题，直接升级所有相关空间包即可解决：

install.packages(c("gstat", "sf", "stars", "sp"), dependencies = TRUE)

你当前使用的EPSG:4326为地理坐标系，距离单位是度，拟合变异函数时设置的range=200对应的实际距离约为200*111km，若你的研究区域范围较小，建议先将坐标系转换为以米为单位的投影坐标系（如对应区域的UTM坐标系），再执行插值，距离计算精度会更高。

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