nmc_met_io icon indicating copy to clipboard operation
nmc_met_io copied to clipboard

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

Open rossbyman opened this issue 5 years ago • 2 comments

你好,用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

rossbyman avatar Oct 22 '20 04:10 rossbyman

读取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',

heygy avatar Nov 17 '20 08:11 heygy

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

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

NMC-DAVE avatar Jan 14 '21 15:01 NMC-DAVE