GEOIMREAD读取geotiff或geojp2图像的子区域。


%的语法

A = geoimread(文件名)
A = geoimread(filename,xlim,ylim)
A = geoimread(filename,latlim,lonlim)
A = geoimread(…，buffer)
[A,x,y,I] = geoimread(…)
geoimread(…)


%的描述

A = geoimread(filename)返回文件名给出的完整图像。这
语法相当于A = geotiffread(filename)。

A = geoimread(filename,xlim,ylim)将geotiff文件的区域限制为
由xlim和ylim给出的极限，它们是相对的地图单位(通常是米)
到数据投影。例如，如果geotiff投影在德州中心
映射系统/兰伯特共形坐标，xlim和ylim将有单位
米相对原点(100 W, 18 N)。xlim和ylim可以是多维的，
在这种情况下，地图的极限将被视为极限的极限
xlim, ylim中所有点的分布。

A = geoimread(filename,latlim,lonlim)如果xlim, ylim中没有值
纬度和经度的正常值，geoimread假定您已经输入
纬度和经度在地理坐标上的界限。第一个
输入是纬度，第二个输入是经度。

将一个以地图单位(通常是米或英尺)为单位的缓冲区添加到
感兴趣区域的界限。如果您想加载一个图像，这可能是有用的
周围分散的晚/长数据。如果你想要额外的2公里的图像
您的数据，输入2000作为缓冲区。如果buffer是一个双元素向量，则第一个
元素应用于图像的左和右范围，第二个元素
应用于图像的顶部和底部区段。

[A,x,y,I] = geoimread(…)还返回
I是projfwd和projinv的有用输入。

没有输出的geoimread(…)显示没有加载的输出图像A
工作区中的任何数据。


%的例子:

显示一个完整的geotiff:
geoimread(“boston.tif”);

将上述结果与一个子集geotiff进行比较:
Mapx = [765884 766035 766963];%单位为英尺
map = [2959218 2957723 2958972];
geoimread (boston.tif, mapx mapy)

%或者如果你有纬度/纬度坐标，你想要一个500英尺的缓冲区:
Lat = [42.3675288 42.3634246 42.3668397];
Lon = [-71.0940009 -71.0934685 -71.0900125];

geoimread (boston.tif,纬度,经度,500);