geotiff或geojp2 GEOIMREAD读取子地区的形象。


%的语法

一个= geoimread(文件名)
一个= geoimread(文件名,xlim ylim)
一个= geoimread(文件名,latlim lonlim)
缓冲= geoimread (……)
[A, x, y,我]= geoimread (…)
geoimread (…)


%的描述

A = geoimread(文件名)返回完整的图像由一个文件名。这
语法是相当于= geotiffread(文件名)。

= geoimread(文件名,xlim ylim)限制了geotiff区域文件
的限制由xlim和ylim地图单位(通常是米)相对的
对数据的投影。例如,如果geotiff预计在德克萨斯州为中心的
映射系统/兰伯特正形坐标,xlim和ylim单位
米相对于原点(100 W, 18 N)。xlim和ylim multimensional,
在这种情况下,地图的限制将被视为的极限的极限
所有点的分布xlim ylim。

= geoimread(文件名,latlim lonlim)如果没有值在xlim, ylim超过
正常的价值观纬度和经度,geoimread假设你输入
限制度经度和纬度的地理坐标。第一个
输入纬度,第二个输入经度。

= geoimread(…,缓冲)添加一个缓冲在地图单位(通常是米和英尺)
感兴趣的区域的局限性。这可能是有用的,如果你想加载一个图像
周围分散纬度/经度数据。如果你想要一个额外的图像周围2公里
你的数据,输入2000作为缓冲。如果缓冲区是一个双元素向量,第一
元素应用于图像的左和右区段,第二个元素
应用于图像的顶部和底部区段。

[A, x, y,我]= geoimread(…)也返回像素中心坐标(x, y)
输出图像和geotiff信息结构。我是一个有用的输入projfwd和projinv。

geoimread(…)没有任何输出显示了输出图像没有装载
任何数据到工作区中。


%的例子:

%显示整个geotiff:
geoimread (“boston.tif”);

%比较结果从上面geotiff的一个子集:
mapx = (765884 766035 766963);%单位是英尺
mapy = (2959218 2957723 2959218);
geoimread (boston.tif, mapx mapy)

%或如果你有在纬度/经度坐标,你想要一个500英尺的缓冲区:
lat = (42.3675288 42.3634246 42.3668397);
朗= (-71.0940009 -71.0934685 -71.0900125);

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