Некогерентная обработка координат в базовой карте addcyclic - PullRequest
Новая
       28

Некогерентная обработка координат в базовой карте addcyclic

0 голосов
/ 18 января 2019

Я успешно использовал все функции mpl_toolkits.basemap в нескольких глобальных наборах данных GCM netcdf. Пока я не встретил эту сетку, с долготами, начинающимися в 0.9375 (вместо 0, как я всегда видел) и заканчивающимися в 359.062.

Чтобы подготовить сюжет, мне нужно:

  1. сделать сюжет непрерывным с: 

    # input_var is a 2D numpy array
var_cyclicDUMMY, lons_cyclicDUMMY = addcyclic(input_var, lons)

    Таким образом, я получаю двумерный массив var_cyclicDUMMYwith дополнительный столбец (одна дополнительная долгота) и одномерный массив lons_cyclicDUMMY с одним дополнительным элементом в конце, т.е. одной дополнительной долготой, но с 0,9375 вместо 360, которые нужно. Действительно на следующем шаге, где я

  2. смещение сетки, поэтому долгота изменяется от -180 до 180 вместо 0 до 360, с: 

    var_cyclic, lons_cyclic = shiftgrid(180., var_cyclicDUMMY, 
                                    lons_cyclicDUMMY, start=False)

    Я получаю ValueError: lon0 outside of range of lonsin

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

...