lat and lon from GEOTIFF when the matrix is not square

2 次查看（过去 30 天）
Farz Tak
Farz Tak 2016-5-20
编辑： 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

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

采纳的回答

Farz Tak
Farz Tak 2016-5-20
编辑：Farz Tak 2016-5-20
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;

更多回答（0 个）

类别

在 Help Center 和 File Exchange 中查找有关 Geographic Coordinate Reference Systems 的更多信息

标签

产品

Community Treasure Hunt

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

Start Hunting!

Translated by