Re-gridding matrix from one projection to another (an equal area projection to lat/long projection)

2 views (last 30 days)
I am currently working on some remote sensing datasets and I am attempting an intercomparison between them, thus I hope to regrid all of them onto the same projection. I am very much a beginner to image analysis and Matlab, so please pardon any weird-sounding or non-technical terms used here!
Datasets A and B are on latitudinal/longitudinal projection, where one pixel in a matrix represents 0.05 degrees x 0.05 degrees (A) and 0.04 x 0.04 degrees (B) of lat/long area respectively. The final datasets, Dataset C (a 586 x 1383 size matrix), is on an "EASE-GRID" projection that is an equal area projection, ie. one pixel = 25km by 25km of area on Earth. I wish to change Dataset A into a projection which is a latitude/longtitude projection, ie. one pixel = a certain lat/long spatial resolution (e.g. degrees x degrees) of the Earth's surface). that preserves as much as possible the resolution of the current dataset.
For dataset C, I have two separate matrices that describes the corresponding latitutdes and longitudes of the midpoint of each 25km by 25km pixel. Due to the equal area projection, the distance between each latitudinal figure (vertical axis) varies. On the vertical axis: towards the poles (top and bottom of the matrix), there is a difference of up to nearly 2 degrees between each pixel, but this shrinks towards the midpoint of the axis with only 0.1952 degrees between each pixel. On the horizontal axis, the distance between each pixel in terms of longtitudinal degrees is fixed at 0.4 degrees.
How can I regrid dataset C to be on a lat/long projection in Matlab, such that it also possesses a spatial resolution where one pixel is equal to x degrees by x degrees? In addition, how can this be done in a way that minimises the loss of resolution in my data, perhaps by ensuring the finest spatial resolution of 0.1952 degrees is preserved?
Please let me know if I can clarify this any further, thank you!

Answers (1)

Mateusz
Mateusz on 2 Jan 2014
During my searching over the Internet, I found many solutions which allow conversion to EASE-Grid (for example: 1, 2)
Is there any ready-solution for transforming data from original EASE-Grid (in my case EPSG:3410, saved in NetCDF file) or should I try to prepare script myself (do you have any hints for that process)?
Unfortunately, I have no Mapping Toolbox installed, so I can't check if it provides any solutions for my issue.

Community Treasure Hunt

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

Start Hunting!