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!
- Subject Matter Expert
- Posts: 138
- Joined: Fri Apr 23, 2021 9:14 am America/New_York
- Has thanked: 14 times
- Been thanked: 7 times
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?
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.