Matlab快速将影像的二维坐标转换为经纬度坐标

A-Chin

已于 2022-10-25 16:57:00 修改

阅读量3.2k

点赞数 5

分类专栏： matlab 数据处理文章标签： matlab 坐标转换二维坐标经纬度坐标

于 2022-06-27 00:43:23 首次发布

本文链接：https://blog.csdn.net/l_j_kin/article/details/125476617

版权

matlab 同时被 2 个专栏收录

35 篇文章 20 订阅

订阅专栏

数据处理

11 篇文章 4 订阅

订阅专栏

这篇博客介绍了如何通过避免双重循环，利用Matlab的矩阵运算来大幅提升pix2latlon函数处理大量影像坐标转换的效率。作者提供了一种创建X和Y矩阵的方法，然后一次性运行pix2latlon函数，从而显著提高了代码运行速度。

摘要由CSDN通过智能技术生成

一般使用 [lat, lon] = pix2latlon(R,row,col) 逐个将影像的二维坐标转换为经纬度坐标，但这样用到双重循环，运行较慢，代码如下：

[A, RA] = readgeoraster('D:\test.tif');
info = geotiffinfo('D:\test.tif');

A_Lat = zeros(size(A));
A_Lon = zeros(size(A));
for i = 1 : size(A, 1)
    for j = 1 : size(A, 2)
        [A_Lat(i, j), A_Lon(i, j)] = pix2latlon(info.RefMatrix, i, j);
    end
end

为了快速得到影像的经纬度坐标，要避免使用双重循环，利用好Matlab矩阵运算的优势；而且 pix2latlon 函数可以对矩阵进行处理。

所以，在这里一种较快的方法为，先建立二维坐标的X矩阵和Y矩阵（代码中的A_I和A_J），然后仅需运行一次 pix2latlon 函数即可，运行效率能够大幅提高，代码如下：

[A, RA] =  geotiffread('D:\test.tif');
info = geotiffinfo('D:\test.tif');

A_I = zeros(size(A));
A_J = zeros(size(A));

for i = 1 : size(A_I, 1)
    A_I(i, :) = i;
end
for j = 1 : size(A_J, 2)
    A_J(:, j) = j;
end

[A_Lat, A_Lon] = pix2latlon(info.RefMatrix, A_I, A_J);