嘿，你已经搞定了RasterBrick到NetCDF的转换，就差把米制坐标转成十进制度了对吧？这事儿不难，核心就是先给栅格补对原始投影信息，再转成地理坐标系（WGS84）就行，下面给你两种常用的实现方式：

方法一：用terra包（推荐，替代旧的raster包） #

terra是现在R空间数据处理的主流工具，速度更快，功能更全。

# 先加载包
library(terra)

# 假设你的原始RasterBrick叫raster_brick，先转成SpatRaster对象
spat_rast <- rast(raster_brick)

# 关键步骤：给栅格设置原始的米制投影（必须替换成你数据实际的投影！比如UTM 50N是EPSG:32650）
crs(spat_rast) <- "EPSG:32650"

# 转换到WGS84地理坐标系（十进制度）
spat_rast_latlon <- project(spat_rast, "EPSG:4326")

# 导出成带正确经纬度的NetCDF
writeCDF(spat_rast_latlon, "output_latlon.nc", overwrite = TRUE)

方法二：用传统的raster包 #

如果你还在使用旧版的raster包，也可以这么操作：

library(raster)

# 给原始RasterBrick设置投影（同样要替换成实际的EPSG码）
proj4string(raster_brick) <- CRS("+init=EPSG:32650")

# 转换到WGS84
raster_brick_latlon <- projectRaster(raster_brick, crs = CRS("+init=EPSG:4326"))

# 导出NetCDF
writeRaster(raster_brick_latlon, "output_latlon.nc", format = "CDF", overwrite = TRUE)

重要提醒 #

• 必须确定原始栅格的正确投影：上面代码里的EPSG:32650只是示例，你得根据数据来源找到对应的投影参数（比如UTM带号、其他投影的PROJ字符串），不然转换出来的经纬度肯定不对。如果不知道原始投影，得去查原始栅格的文档或者数据提供者。
• 如果已经有了生成好的NetCDF文件，不想重新从RasterBrick处理，也可以用ncdf4包直接修改坐标，但这种方法只适合规则网格，不如重新处理栅格再导出可靠。

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