Блог

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

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

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

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

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

Вот график on ray в моем массиве долгот

 plot(lon[100,:])
  

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

Вы могли бы попробовать использовать 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,:] .