Проблемы с координатами и перевернутая карта
Я использую данные о солености ORAS5, а также переменные и информацию предварительного просмотра:
salinity =xr.open_dataset('....salinity_ORAS5.nc')
sal = salinity['sosaline']
Итак, большая проблема заключается в том, что каждый «nav_lon» и «nav_lat» имеют 2D по x и y. Где в обычном случае, например, данные широты имеют 1D.
sal.nav_lon
В предварительном просмотре карта перевернута:
plt.imshow(sal)
Я попытался преобразовать координаты x/y в широту/долготу, но мои данные сохраняются в 2D, меняются только значения.
Большой вопрос заключается в том, как работать с моими координатами, с долготой (x, y) и широтой (x, y)? Как преобразовать lon в 0-360 и lat в 0-180? И визуализируйте мою карту в нормальной проекции, а не в перевернутой.
Любая помощь приветствуется!
Большое спасибо.
РЕДАКТИРОВАТЬ
Мне нужно решение для работы с сеткой модели Projection: Tripolar. Или работать с трехполярной сеткой модели или интерполировать данные. Каков наилучший способ?
Может ли кто-нибудь помочь мне?
@Giacomo Мне нужно преобразовать координаты в «действительные числа» (от 0 до 180 и от 0 до 360) и правильно визуализировать мою карту, а не перевернутую.
Да, но нужно проверять каждый шаг. Как вы можете себе представить, есть ошибка (и такие ошибки часто делаются из-за того, что я писал ранее). Вы должны просто проверять каждый шаг. Поскольку мы не знаем, что вы делаете на каждом этапе, и какой источник вы используете, мы не можем помочь.
Я не знаю, поможет ли, потому что я просто открываю и рисую свои данные. Но я включу код в пост.
Ответы 1
Мне нужно преобразовать координаты в «нормальную» широту/долготу, используя CDO.
cdo remapbil
Надеюсь, это поможет :)
Трудно понять, что вы делаете, и в чем ошибка: вы не предоставили никакого кода. В любом случае: широта/долгота – это вертикаль/горизонт. x/y является горизонтальным/вертикальным, а geoJSON использует длину/широту (поэтому горизонтальный/вертикальный). Просто внимательно проверьте происхождение ваших данных, а также порядок или координаты.