Использование python Basemap.pcolormesh с немонотонными прыжками долготы

Question

У меня есть данные спутниковой развертки, которые я пытаюсь построить на базовой карте с помощью pcolormesh.

Данные организованы как 2d-матрица (bTemp) с двумя соответствующими 2D-массивами lat и lon, которые дают соответствующую широту и долготу в каждой точке. Таким образом, мой код построения выглядит как

m = Basemap() 
m.drawcoastlines() 
m.pcolormesh(lon, lat, bTemp) 

Однако это не совсем дает мне результат, который я ищу. Большая карта проходит по карте.

Я думаю, что причиной этого является то, что моя долгота увеличивается немонотонно вдоль данного луча в точке обертывания.

Вот сюжет на луче в моей долготе массиве

plot(lon[100,:]) 

Что бы лучший способ исправить это будет, так что pcolormesh сюжет просто переходит к другой стороне карту без соединения двух точек с заполненной областью?

Answer 1

Чтобы развернуть данные, вы можете использовать numpy.unwrap. Он был разработан именно для этого: развернуть угловые прыжки (в радианах), которые являются артефактами из-за того, как мы используем однозначные функции.

Вам необходимо конвертировать ваши долгот назад и вперед между степенями и радианах для этой работы, но суть заключается в следующем:

import numpy as np 
lon_unwrapped = np.rad2deg(np.unwrap(np.deg2rad(lon))) 

Найдет эти скачки в ваших данных, размер которых превышает 180 градусов, и принять их 360-градусный дополнение. Пример:

np.rad2deg(np.unwrap(np.deg2rad([150,-150]))) 
# array([ 150., 210.]) 

Кроме того, unwrap принимает ndarray данных, и вы можете сказать, что, по которой оси для выполнения распутывания. По умолчанию это последнее измерение, которое, кажется, подходит для вашего случая, где вы хотите развернуть lon[i,:].