TypeError: Cannot handle this data type: (1, 1, 72), |u1还是这样报错

这个错误提示是由于Pillow库无法处理多波段的数据类型。在这种情况下，你可以将多波段的数据拆分成单波段的数据，分别转换为Pillow库可以处理的数据类型，然后再将它们合并成多波段数据。可以按照以下步骤修改代码：

1. 将多波段的数据拆分成单波段的数据，使用numpy库的split方法，将数据按照波段维度拆分：

import numpy as np

# 按照波段维度拆分数据
bands = np.split(da_resampled.values, da_resampled.shape[-1], axis=-1)

2. 将单波段的数据转换为Pillow库可以处理的数据类型，使用numpy库的astype方法将数据类型转换为uint8：

# 将每个单波段数据转换为Pillow库可以处理的数据类型
bands = [Image.fromarray(band.squeeze().astype('uint8')) for band in bands]

3. 将转换后的单波段数据合并成多波段数据，使用Pillow库的merge方法将单波段数据合并：

# 将单波段数据合并成多波段数据
img = Image.merge('RGB', bands)

最后，使用rio库的write方法将数据和缩略图一起写入到文件中：

with rasterio.open(out_root / f'{nc_path.stem}.tif', 'w', driver='GTiff', crs=ref_da.crs, transform=ref_da.transform) as dst:
    dst.write(da_resampled.values, 1)
    dst.update_tags(ns='rio', res=str(ref_da.res[0])) # 将分辨率写入文件属性
    dst.write(img.tobytes(), 1, photometric='RGB', thumbnail=(256, 256)) # 将缩略图写入文件属性

这样就可以解决多波段数据类型无法处理的问题了。