Mercator map projection from [X,Y] with known lat/lon coordinates?

Question

Peter Dorofy 2015-7-15

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/229984-mercator-map-projection-from-x-y-with-known-lat-lon-coordinates

评论： Peter Dorofy 2015-7-20

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 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

请先登录，再回答此问题。

Answer 1

Chad Greene 2015-7-15

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/229984-mercator-map-projection-from-x-y-with-known-lat-lon-coordinates#answer_186169

在 MATLAB Online 中打开

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.