nmc_met_io.retrieve_micaps_server.get_fy_awx()读FY4A C012红外数据经纬度问题

[rossbyman]

你好，用get_fy_awx()读出来的数据，经纬度范围与MICAPS4显示的经纬度不一致，是什么问题？

• lat (lat) float64 55.0 54.95 54.9 54.85 54.8 ... -7.35 -7.4 -7.45 -7.5
• lon (lon) float64 70.0 70.05 70.1 70.15 ... 157.3 157.4 157.4 157.5

[heygy]

读取FY4A C009数据报错

directory='SATELLITE/FY4A/L1/CHINA/C009/'
filename='C009_20201117080000_FY4A.AWX'
IR = get_fy_awx(data_dir, filename=filename)

Exception has occurred: ValueError
Could not convert tuple of form (dims, data[, attrs, encoding]): ('channel', 2, {'long_name': 'channel', 'units': ''}) to Variable.
File "L:\py_develop\nmc_met_map\nmc_met_map\observation.py", line 41, in IR_Sounding_GeopotentialHeight
IR = get_fy_awx(data_dir[0], filename=filename_IR)
File "L:\py_develop\nmc_met_map\nmc_met_map\test_draw.py", line 24, in
draw_observation.IR_Sounding_GeopotentialHeight(Sounding_time='20111708',

[NMC-DAVE]

能够正常读取，并且经纬度范围也没有问题了。

from nmc_met_io.retrieve_micaps_server import get_fy_awx

directory='SATELLITE/FY4A/L1/CHINA/C009/'
filename='C009_20210114223418_FY4A.AWX'
IR = get_fy_awx(directory, filename=filename)
IR

<xarray.Dataset>
Dimensions: (channel: 1, lat: 1001, lon: 1751, time: 1)
Coordinates:

• time (time) datetime64[ns] 2021-01-14T14:34:00
• channel (channel) int16 2
• lat (lat) float64 15.0 15.04 15.08 15.12 ... 54.88 54.92 54.96 55.0
• lon (lon) float64 70.0 70.04 70.08 70.12 ... 139.9 139.9 140.0 140.0
  Data variables:
  image (time, channel, lat, lon) float64 243.1 243.1 243.1 ... 249.2 248.7
  Attributes:
  Conventions: CF-1.6
  Origin: MICAPS Cassandra Server