项目场景:
根据经纬度数据和高程数据,进行数据采样,得到目标经纬度分辨率的高程数据
问题描述
用pandas和numpy都找不到实现该功能的函数
原因分析:
pandas和numpy没有实现该功能的函数
解决方案:
使用xarray
import numpy as np
import xarray as xr
lon_prov=np.array([105+i0.025 for i in range(0,30)])
lat_prov=np.array([37+i0.025 for i in range(0,30)])
realhei_prov = ‘dem.nc’ //分辨率0.00025
da = xr.open_dataset(realhei_prov)
dem_prov = da[‘dem’]
dem_prov = dem_prov.interp(lat=lat_prov, lon=lon_prov)