确定R中的两个格网是否完全匹配

library(raster)

# 假设r1和r2是两个raster对象
if (dim(r1) == dim(r2)) {
  print("维度匹配")
} else {
  print("维度不匹配")
}

if (res(r1) == res(r2)) {
  print("分辨率匹配")
} else {
  print("分辨率不匹配")
}

if (crs(r1) == crs(r2)) {
  print("坐标参考系统匹配")
} else {
  print("坐标参考系不匹配")
  # 尝试将r2的CRS转换为r1的CRS
  r2 <- projectRaster(r2, crs = crs(r1))
}

library(raster)

# 假设r1和r2是两个raster对象
check_grid_match <- function(r1, r2) {
  if (dim(r1) == dim(r2) && res(r1) == res(r2) && crs(r1) == crs(r2)) {
    print("两个格网完全匹配")
  } else {
    print("两个格网不完全匹配")
    if (dim(r1) != dim(r2)) {
      print("维度不匹配")
    }
    if (res(r1) != res(r2)) {
      print("分辨率不匹配")
    }
    if (crs(r1) != crs(r2)) {
      print("坐标参考系统不匹配")
      # 尝试将r2的CRS转换为r1的CRS
      r2 <- projectRaster(r2, crs = crs(r1))
    }
  }
}

# 示例调用
check_grid_match(r1, r2)