• geotiffred函数

[A, R] = geotiffread(filename)

其中:

  • filename 是要读取的 GeoTIFF 文件的路径和名称。
  • A 是返回的地理影像数据,以 MATLAB 数组的形式表示。对于单波段影像,A 是一个二维数组;对于多波段影像,A 是一个三维数组,其中第三个维度表示波段数。
  • R 是返回的地理参考对象(Geographic Referencing Object),包含了地理坐标和像素坐标之间的映射关系,可以用于地理空间数据的处理和分析。

  • geotiffinfo函数

R_info = geotiffinfo ( 'boston.tif' )

其字段包含有关 GeoTIFF 文件的图像属性和制图信息,其中包含了最重要的影像的投影坐标参考系统。具体各个字段的详细描述可查看Information about GeoTIFF file - MATLAB geotiffinfo- MathWorks 中国 网页。

  • readgeoraster函数和georasterinfo函数

readgeoraster 函数是 MATLAB R2021a 中引入的新函数,它能够读取多种格式的栅格数据,包括 GeoTIFF、JPEG、PNG、BMP 等。可以说是geotiffread的升级版,用法与geotiffread函数一致。


% 读取 GeoTIFF 影像数据
[data, R] = readgeoraster('A.tif');
info = georasterinfo('A.tif'); %与geotiffinfo函数类似
  • ProjectedCRS

[A,R]=readgeoraster('A.tif');   
R_proinfo=R.ProjectedCRS; %影像使用的投影坐标系统使用ProjectedCRS查看

用于提取影像数据中的投影坐标参考系统。

  • map2pix函数

用于将投影坐标转换为像素坐标

[px, py] = map2pix(R, x_coordinate, y_coordinate) 
 %一般x_coordinate为经度lon,y_coordinate为维度lat

%worldToIntrinsic函数用法一致
[px, py] = worldToIntrinsic(R, x_coordinate, y_coordinate) 

其中:

  • R 是一个地理参考对象(Geographic Referencing Object),表示图像或地理数据的地理参考信息。R 可由geotiffred函数读取。
  • y_coordinate表示要转换为像素坐标的点的投影坐标系下的纬度lat。
  • x_coordinate 表示要转换为像素坐标的点的投影坐标系下的经度lon。
  • pxpy 是返回的像素坐标数组,与输入的经纬度数组大小相同,px为行号,py为列号

  • projfwd函数

对于有些坐标系为地理坐标系而非投影坐标系的坐标信息,可以使用projfwd函数将将地理坐标点(经度和纬度)投影到投影坐标系下

[x, y] = projfwd(proj, lat, lon)

其中:

  • proj 是投影结构体,包含了所使用的投影方法和相关参数。你可以使用 projcrs 函数创建投影结构体,或者使用geotiffinfo函数读取TIF影像的信息
  • lat 是一个标量或一个大小相同的经度数组,表示要进行投影的点的纬度。
  • lon 是一个标量或一个大小相同的经度数组,表示要进行投影的点的经度。
  • xy 是返回的平面坐标数组,与输入的经纬度数组大小相同。

  • projinv函数

% 假设你有一个投影坐标点的 x 和 y 值
x = 12345;
y = 67890;

% 创建一个地理坐标参考对象,指定投影信息
projcrs = projcrs('EPSG:4326'); % 使用 WGS84 地理坐标系

% 将投影坐标转换为经纬度坐标
[lat, lon] = projinv(projcrs, x, y);

这个函数用于执行逆投影,将投影坐标转换为地理坐标(经纬度)

Logo

魔乐社区(Modelers.cn) 是一个中立、公益的人工智能社区,提供人工智能工具、模型、数据的托管、展示与应用协同服务,为人工智能开发及爱好者搭建开放的学习交流平台。社区通过理事会方式运作,由全产业链共同建设、共同运营、共同享有,推动国产AI生态繁荣发展。

更多推荐