将NOAA NetCDF文件转为GeoTIFF时坐标异常问题排查

你尝试将NOAA的NetCDF文件转成GeoTIFF时遇到两个问题：

• 用rioxarray的open_rasterio转换后，QGIS无法查看内容
• 用gdalwarp指定CRS后能查看，但数据显示在非洲而非正确的美国区域

问题原因 #

1. rioxarray读取方式错误：rioxarray.open_rasterio适用于符合栅格规范的文件，而该NetCDF的变量依赖独立的x/y坐标维度，需要用xarray读取后再赋予地理信息。
2. GDAL操作缺少地理变换信息：你仅设置了CRS，但未指定数据的实际地理范围（坐标偏移、分辨率），导致GDAL默认将数据放在投影原点附近（对应非洲区域）。

解决方案 #

方法一：用rioxarray正确处理

使用xarray读取NetCDF，手动赋予CRS和地理变换参数：

import rioxarray
import xarray as xr

# 读取NetCDF数据集
ds = xr.open_dataset('path/to/newhail.nc')
hail_data = ds['newhail']

# 从全局属性提取CRS并设置
crs_proj = ds.attrs['map_proj']
hail_data = hail_data.rio.set_crs(crs_proj)

# 计算并设置地理变换（基于x/y坐标的分辨率和起始值）
x_start = ds['x'].values[0]
x_res = ds['x'].values[1] - ds['x'].values[0]
y_start = ds['y'].values[0]
y_res = ds['y'].values[1] - ds['y'].values[0]
# 地理变换格式：(左上x, x分辨率, 旋转, 左上y, 旋转, y分辨率)，y分辨率为负因为纬度从北到南递减
geo_transform = (x_start, x_res, 0, y_start, 0, y_res)
hail_data = hail_data.rio.set_transform(geo_transform)

# 保存为GeoTIFF
hail_data.rio.to_raster('path/to/correct_newhail.tif')

方法二：用GDAL命令行直接转换

先通过gdal_translate读取NetCDF变量，同时指定CRS和地理范围：

# 先获取x/y的极值：可以用Python或gdalinfo查看
# 该数据集的x范围约为-3390000到2610000，y范围约为1920000到5880000

gdal_translate -of GTiff \
  -a_srs "+proj=lcc +lat_0=25 +lon_0=265 +lat_1=25 +lat_2=25 +x_0=0 +y_0=0 +R=6371200 +units=m +no_defs" \
  -a_ullr -3390000 5880000 2610000 1920000 \
  NETCDF:"path/to/newhail.nc":newhail \
  newhail_correct.tif

验证 #

转换完成后在QGIS中打开，数据应该正确覆盖美国区域，坐标匹配投影定义。

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