Mercator map projection from [X,Y] with known lat/lon coordinates?
8 views (last 30 days)
Show older comments
Peter Dorofy
on 15 Jul 2015
Commented: Peter Dorofy
on 20 Jul 2015
I am working with a satellite image in simple binary, pixels are stored as 2-byte integers. I've used the code below to build a matrix. I am able to display the satellite's "native" view. What I need to do is map the X,Y to a mercator projection. I have the upperleft and lowerright lat/lon coordinates, so I am thinking this could probably be done, but I am not sure how or what functions I should be looking at. Any suggestions?
Thanks.
img = multibandread(file, [3712 3712 1], 'int16', 0, 'bil', 'ieee-le');
m = img(ROWa:ROWb,COLa:COLb);
[X,Y] = ind2sub(size(m), 1:numel(m));
Z = m(:);
0 Comments
Accepted Answer
Chad Greene
on 15 Jul 2015
If you have the Mapping Toolbox, you could initialize a mercator map, then use projfwd to get the x/y locations of the corner coordinates. If you need gridded x/y coordinates, you could then use
[X,Y] = meshgrid(linspace(minx,maxx,cols),linspace(miny,maxy,rows))
where minx, maxx are the corner x coordinates and cols is the number of columns.
More Answers (0)
See Also
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!