MATLAB Answers

Interpolation of Data on Different Grids

9 views (last 30 days)
David C
David C on 18 Nov 2014
Commented: Chad Greene on 20 Nov 2014
I have a 2-D data (ocean depth) on a latitude-longitude grid. (That is, the ocean depth is specified at each latitude and longitude. The latitude coordinates and longitude coordinates are separated by 1 deg each.) Question: How do I interpolate this data onto a fixed-distance grid (with an interval of, say, 100 m by 100 m)?
Thanks.

  0 Comments

Sign in to comment.

Accepted Answer

Chad Greene
Chad Greene on 20 Nov 2014
The best way is to project the elevation data using an appropriate projection. Then interpolation should be no problem with interp2. Choosing a projection will depend on where the data lie on the globe (near the equator?, near the poles?) and how large an area the data spans. If your bathymetry data spans only a few hundred kilometers, perhaps a UTM projection is appropriate. There's a pretty good explanation of projections here. If you have the mapping toolbox, you can project and unproject with projfwd and projinv.
Another option, for a small area of the globe, could be to use interpm2 in the lat dimension to get an array of lats, use it again for the lon dimension, and call that your new grid. But again, this would only be appropriate for a small area of the globe where spacing between meridians does not change very much over the north-south extents of the data.

  2 Comments

David C
David C on 20 Nov 2014
Thanks Chad - I knew about interp2, but did not know about the other functions. I'll look into them!
Chad Greene
Chad Greene on 20 Nov 2014
And if you don't have the mapping toolbox, m_map has several helpful functions.

Sign in to comment.

More Answers (0)


Translated by