Pandas DataFrame导出NetCDF文件时出现AttributeError问题求助

我来帮你分析这个问题并给出解决方案：

首先，你遇到的错误原因很明确：Pandas的DataFrame对象并没有to_netcdf这个方法，这个导出NetCDF的方法是xarray库中Dataset和DataArray的专属方法，所以当你调用Newdf.to_netcdf()时，Python会告诉你找不到这个属性。

另外，你的原代码还有两个小问题需要优化：

1. 把to_netcdf放在循环内部，每次循环都重新保存一次文件，这会重复写入、浪费性能，应该把保存操作移到循环外面，只执行一次。
2. 文件路径里混合了/和\，容易出现转义问题，建议用原始字符串r"路径"或者统一用/。

下面给你两种解决思路，推荐第一种（更符合气象数据处理的最佳实践）：

方案1：全程用xarray处理（推荐） #

xarray是专门为NetCDF这类多维气候数据设计的库，直接用它处理和导出会更高效，生成的NetCDF文件结构也更规范（带有站点维度，方便后续分析）：

import xarray as xr 

# 读取原始NetCDF文件
nc_file = r"C:\Users\DELL 3090\Desktop\Projection\RF_58\pr_Amon_AWI-CM-1-1-MR_ssp585_mm_month.nc" 
NC = xr.open_dataset(nc_file)

# 定义站点信息
lat_list = [13.3,7.5,8.5,7.9,9,7.8,11.3,14.3,14.2,7]
lon_list = [39.8,34.3,39.8,38.7,38.8,39.9,37.5,39.5,38.9,39.9]
station_names = ['Abala','Abobo','Abomsa','Adamitulu','AddisAbabaBole','Adele','Adet','Adigrat','Adwa','Agarfa']

# 批量获取每个站点的最近邻数据，并添加站点维度
site_datasets = []
for lat, lon, name in zip(lat_list, lon_list, station_names):
    # 选择最近邻格点
    site_data = NC.sel(lat=lat, lon=lon, method='nearest')
    # 给当前站点数据添加站点名称的坐标
    site_data = site_data.assign_coords(station=name)
    site_datasets.append(site_data)

# 把所有站点数据合并成一个Dataset（新增station维度）
combined_ds = xr.concat(site_datasets, dim='station')

# 导出为NetCDF文件
combined_ds.to_netcdf(
    r'C:\Users\DELL 3090\Desktop\Projection\RF_58\pr_Amon_AWI-CM-1-1-MR_ssp585_mm_month111.nc',
    mode='w'  # 覆盖写入
)

方案2：基于你原代码修改（保留DataFrame处理流程） #

如果你一定要用DataFrame来中转数据，可以先把最终的DataFrame转换成xarray的Dataset，再调用to_netcdf方法：

import pandas as pd
import xarray as xr 

nc_file = r"C:\Users\DELL 3090\Desktop\Projection\RF_58\pr_Amon_AWI-CM-1-1-MR_ssp585_mm_month.nc" 
NC = xr.open_dataset(nc_file)

# 定义站点信息
lat=[13.3,7.5,8.5,7.9,9,7.8,11.3,14.3,14.2,7]
lon=[39.8,34.3,39.8,38.7,38.8,39.9,37.5,39.5,38.9,39.9]
name=['Abala','Abobo','Abomsa','Adamitulu','AddisAbabaBole','Adele','Adet','Adigrat','Adwa','Agarfa']

Newdf = pd.DataFrame([])

for i,j,id in zip(lat,lon,name): 
    dsloc = NC.sel(lat=i,lon=j,method='nearest') 
    DT = dsloc.to_dataframe()
    # 添加站点名称列，方便后续区分不同站点的数据
    DT['station'] = id
    Newdf = Newdf._append(DT, sort=True)

# 关键：把DataFrame转成xarray Dataset再导出NetCDF
xr.Dataset.from_dataframe(Newdf).to_netcdf(
    r'C:\Users\DELL 3090\Desktop\Projection\RF_58\pr_Amon_AWI-CM-1-1-MR_ssp585_mm_month111.nc',
    mode='w'
)

备注：内容来源于stack exchange，提问作者mohammed abera