IDL:显示带有地图投影的卫星数据

时间:2020-03-16 20:55:36

标签: arrays dictionary projection idl satellite

我有3个尺寸相同的数组(3636,3636): 1-具有云顶高度值的数组(以米为单位) 2-具有经度值的数组 3-具有纬度值的数组 我想显示带有地图投影,大洲和国家/地区边界以及地图网格的第一个数组的值。我需要使用第二和第三数组。我不知道如何匹配这3个数组。 我使用了很多功能,例如。图片,map_image,轮廓,map_proj_image,iimage,但它们都无法正常工作。
在这里,这是我使用的代码的一部分:

lat = MAKE_ARRAY(3636,3636)
lon = MAKE_ARRAY(3636,3636)
cth = MAKE_ARRAY(3636,3636)

wsp = NCDF_OPEN(plik, /nowrite)
wys = NCDF_OPEN(ctx, /nowrite)
NCDF_VARGET, wsp, 'lat', lat
NCDF_VARGET, wsp, 'lon', lon
NCDF_VARGET, wys, 'cth', cth
NCDF_CLOSE, wsp
NCDF_CLOSE, wys

image2 = IMAGE(re, RGB_TABLE=40, MARGIN=[0.1,0.02,0.08,0.02], MAPPROJECTION='GOES-R' $        GRID_UNITS='degrees', DIMENSIONS=[1000,600], IMAGE_DIMENSIONS=[162,162], IMAGE_LOCATION=[-82,-82])

map = MAP('orthographic', CENTER_LONGITUDE=0, CENTER_LATITUDE=0, /OVERPLOT, COLOR='white')
mc = MAPCONTINENTS(/COUNTRIES, COLOR='white')

在这段代码中,我没有使用具有地理坐标的数组。因此,在该投影中,必须在中欧上空的云层才在斯堪的纳维亚半岛上空。

您是否有解决此问题的想法?

0 个答案:

没有答案