def nearest(self, lat, lon, var):
output = {}
dims, idx = cropIndices(self.dims, lat, lon)
for v in var:
if v == 'height':
v = 'z'
subset = self.ncs[0].variables[v][idx['yn'], idx['xn']]
output[v] = ma.masked_all((lat.size, lon.size), dtype='f')
for yn_out, y in enumerate(lat):
yn_in = np.absolute(dims['lat']-y).argmin()
for xn_out, x in enumerate(lon):
xn_in = np.absolute(dims['lon']-x).argmin()
output[v][yn_out, xn_out] = subset[yn_in, xn_in]
return output
评论列表
文章目录