.nc文件根据经纬度提取点上数值python

1、数据集

VOD Climate Archive （VODCA） 数据由维也纳工业大学提供（https://doi.org/
10.5281/zenodo.2575599）。VODCA是一种空间分辨率为0.25°的全球每日VOD产品，其
基于多个传感器 （SSM/I、TMI、AMSR-E、WINDSAT 和 AMSR2） 的 VOD 检索结果的
融合，使用AMSR-E作为缩放参考，通过累积分布函数匹配进行共同校准，创建了C波
段 （~6.9 GHz、2002.6—2018.12）、X 波段 （10.7 GHz、1997.12—2018.12） 和 Ku 波段
（~19 GHz、1987.7—2017.7） 的长期 VOD 数据集。
VOD数据集下载

可以看出来，这个数据是三维数据(time维、lat纬度、lon经度)。

2、代码部分

首先，我们加载一下刚刚下载好的数据

# 提取VOD数据
import pandas as pd
import