nc文件的xarray Dataset转pandas DataFrame报NotImplementedError如何处理

报错原因 #

该报错由坐标维度的非法值触发，除了你发现的-243069-06-09 17:00:00异常时间戳外，经度维度末尾的35.0、纬度维度末尾的9.9498e-40也不属于正常范围，非法值导致xarray转换为pandas MultiIndex时触发了类型处理的底层bug。

解决方案 #

步骤1：清理异常坐标值

import xarray as xr
import pandas as pd
import numpy as np

# 1. 清理异常时间：仅保留可正常解析为datetime的时间值
valid_time = DS1.time.astype('datetime64[ns]', errors='coerce').notnull()
DS_clean = DS1.sel(time=valid_time)

# 2. 清理异常经纬度：根据数据集正常范围（纬度30~35、经度15~20）过滤
valid_lat = (DS_clean.latitude >= 30) & (DS_clean.latitude <= 35)
valid_lon = (DS_clean.longitude >= 15) & (DS_clean.longitude <= 20)
DS_clean = DS_clean.sel(latitude=valid_lat, longitude=valid_lon)

步骤2：转换为DataFrame

清理完成后直接调用内置方法即可：

df = DS_clean.to_dataframe()

如果仍有报错，可以手动展开三维数据构造DataFrame，绕过xarray内置转换逻辑：

# 生成三维网格坐标
time_grid, lat_grid, lon_grid = np.meshgrid(
    DS_clean.time.values, 
    DS_clean.latitude.values, 
    DS_clean.longitude.values, 
    indexing='ij'
)

# 构造DataFrame
df = pd.DataFrame({
    'time': time_grid.ravel(),
    'latitude': lat_grid.ravel(),
    'longitude': lon_grid.ravel(),
    'u10': DS_clean['u10'].values.ravel(),
    'v10': DS_clean['v10'].values.ravel(),
    't2m': DS_clean['t2m'].values.ravel(),
    'tp': DS_clean['tp'].values.ravel()
})

# 可选：去除全为空值的无效行
df = df.dropna(subset=['u10', 'v10', 't2m', 'tp'], how='all')

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