Page 1 of 1

anti-meridian crossing with MATLAB

Posted: Fri Oct 29, 2021 10:08 am America/New_York
by dijudi89
Dear all,

I am wondering about a solution to work with a subset of data crossing anti-meridian. Currently, I am extracting data just like

J=find(lat > -30 & lat < 0);
K=find(lon > 170 & lon < 180);

and cannot cross < 180.

Help is very much appreciated!
Many thanks,
Judith

Re: anti-meridian crossing with MATLAB

Posted: Fri Oct 29, 2021 3:18 pm America/New_York
by ASDC - ingridgs
Hello Judith,

Thank you for contacting the Forum.
Can you please let us know what product are you working on so I can direct you to the right expert?

Thank you,

Re: anti-meridian crossing with MATLAB

Posted: Sat Oct 30, 2021 5:32 am America/New_York
by dijudi89
Many thanks for your quick reply. I am working on GPM IMERG Final Precipitation L3 Half Hourly 0.1 degree x 0.1 degree V06 (GPM_3IMERGHH) at GES DISC.

Re: anti-meridian crossing with MATLAB

Posted: Wed Nov 03, 2021 9:24 am America/New_York
by GES DISCx - zliu
Thank you for using IMERG from GES DISC.

Our data subsetter interface also has the similar restriction,

“East coordinate must be greater than west coordinate AND bounding box should not cross anti-meridian or geographical poles.”

However, Giovanni does allow this by using fill values such as,

https://giovanni.gsfc.nasa.gov/giovanni ... itationCal

Seems different software packages have different solutions to this question. In GrADS (http://cola.gmu.edu/grads/), one can use (0-360) for this situation, for example.

Have you contacted the Matlab Help Center? We have very limited support for Matlab. The HDF Group currently maintains a list of examples in different codes including Matlab for IMERG, https://hdfeos.org/zoo/index_openGESDISC_Examples.php

Hope these help.