lat and lon from GEOTIFF when the matrix is not square

Question

Farz Tak 2016-5-20

0
链接

此问题的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/285204-lat-and-lon-from-geotiff-when-the-matrix-is-not-square

编辑： Farz Tak 2016-5-20

Hi, I have my geotiff file and want to extract lat and lon as: info = geotiffinfo(name); AX = info.Width; AY = info.Height; [y,x] = pix2map(info.RefMatrix, 1:AY, 1:AX); [lat,lon] = projinv(info, y,x);

But my AX and AY are not the same and it gives the error. Any help?

Thanks,

Farz

0 个评论
显示 -2更早的评论隐藏 -2更早的评论

请先登录，再进行评论。

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

Answer 1

Farz Tak 2016-5-20

0
链接

此回答的直接链接

https://ww2.mathworks.cn/matlabcentral/answers/285204-lat-and-lon-from-geotiff-when-the-matrix-is-not-square#answer_222928

编辑：Farz Tak 2016-5-20

在 MATLAB Online 中打开

 I found it after a day:
info = geotiffinfo(path11);
function [lat,lon] = geo_coordinate(info)
        AX = info.Width;
        AY = info.Height;
        [x1,y1] = pix2map(info.RefMatrix, 1:AY, ones(AY,1)');
        [lat1,lon1] = projinv(info, x1,y1);
        [x2,y2] = pix2map(info.RefMatrix, ones(AX,1)',1:AX);
        [lat2,lon2] = projinv(info, x2,y2);
        [lat,lon] = ndgrid(lat1,lon2);
        1;