本文整理汇总了Python中xarray.Dataset.attrs['lat']方法的典型用法代码示例。如果您正苦于以下问题:Python Dataset.attrs['lat']方法的具体用法?Python Dataset.attrs['lat']怎么用?Python Dataset.attrs['lat']使用的例子?那么恭喜您, 这里精选的方法代码示例或许可以为您提供帮助。您也可以进一步了解该方法所在类xarray.Dataset
的用法示例。
在下文中一共展示了Dataset.attrs['lat']方法的1个代码示例,这些例子默认根据受欢迎程度排序。您可以为喜欢或者感觉有用的代码点赞,您的评价将有助于系统推荐出更棒的Python代码示例。
示例1: radec2azel
# 需要导入模块: from xarray import Dataset [as 别名]
# 或者: from xarray.Dataset import attrs['lat'] [as 别名]
def radec2azel(scale: xarray.Dataset,
latlon: Tuple[float, float], time: datetime=None) -> xarray.Dataset:
if latlon is None or not isinstance(scale, xarray.Dataset):
return None
if time is None:
with fits.open(scale.filename, mode='readonly') as f:
try:
t = f[0].header['FRAME'] # TODO this only works from Solis?
except KeyError:
logging.error('no time given in file or manually, cannot compute az/el')
return None
time = parse(t)
logging.info('using FITS header for time')
elif isinstance(time, datetime):
pass
elif isinstance(time, (float, int)): # assume UT1_Unix
time = datetime.utcfromtimestamp(time)
else: # user override of frame time
time = parse(time)
print('image time:', time)
# %% knowing camera location, time, and sky coordinates observed, convert to az/el for each pixel
az, el = pymap3d.radec2azel(scale['ra'], scale['dec'], latlon[0], latlon[1], time)
# %% collect output
scale['az'] = (('y', 'x'), az)
scale['el'] = (('y', 'x'), el)
scale.attrs['lat'] = latlon[0]
scale.attrs['lon'] = latlon[1]
scale.attrs['time'] = time
return scale