У меня есть 3 массива с одинаковыми размерами (3636,3636): 1 - массив со значениями высоты облаков (в метрах) 2 - массив со значениями долготы 3 - массив со значениями широты Я хочу отобразить значения из первого массива с картографическая проекция, континенты и границы стран и сетка карты с. Мне нужно использовать 2-й и 3-й массив. Я не знаю, как я могу сопоставить эти 3 массива. Я использовал много функций, например. image, map_image, contour, 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')
В этом коде я не использовал массивы с географическими координатами. Вот почему в этой проекции облака, которые должны быть над Центральной Европой, были над северной Скандинавией.
У вас есть идеи, как решить эту проблему?