How to recenter spatial data?

6 次查看（过去 30 天）
JMSE
JMSE 2021-11-2
评论： JMSE 2021-11-11
Ran in:
I want to recenter spatial data from a grid of -180 to 180 lon to a grid of 0 - 360 lon. Using the command wrapTo360 leads to the issue below (see 2nd figure - you can see that the data > 180 shows an error). How can I reshape the rain/precipitation data accordingly so that the data is shifted correctly?
Help is very much appreciated.
dlat = load('latitude.mat');
lat = dlat.lat;
dlon = load('longitude.mat');
lon = dlon.lon;
dprec = load('rain.mat');
prec = dprec.prec;
lonW=wrapTo360(lon);
pcolor(lon,lat,prec);shading flat
pcolor(lonW,lat,prec);shading flat

请先登录，再回答此问题。

采纳的回答

Konrad
Konrad 2021-11-2
hi,
maybe I didn't get it, but why not just add 180?
lonW = lon + 180;
  3 个评论
显示 1更早的评论
Konrad
Konrad 2021-11-2
编辑：Konrad 2021-11-2
I'm still not sure what you aim to achieve. To you want to preserve the order of your x-data (longitudes)? Then
pcolor(lon+180,lat,prec);shading flat
will give you longitudes from 0 to 360 on the x axis; everything else looks exactly like in your first plot.
In contrast, wrapTo360() changes e.g. -1 to 359, -180 to 180. So basically, the left half of your plot is moved to the right. pcolor() seems to expect sorted values, so you would have to do:
[lonW,idx] = sort(wrapTo360(lon));
precW = prec(:,idx);
figure;pcolor(lonW,lat,precW);shading flat;
Best, Konrad
JMSE
JMSE 2021-11-11
Many thanks, it helped a lot.

请先登录，再进行评论。

更多回答（0 个）

类别

在 Help Center 和 File Exchange 中查找有关 Lengths and Angles 的更多信息

标签

产品


版本

R2021a

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by