马赛克空间参考光栅瓷砖
地理空间栅格数据提供者通常包数据作为相邻的瓷砖。例如,SRTM高程数据打包成1-degree-by-1-degree瓷砖来填补。分析数据分布在多个瓷砖,如计算视线可见性时,你必须马赛克瓷砖到单个光栅。
结合细胞和位图的位图的流程发布点是不同的。为区分细胞组成的位图信息和发布点,看到的空间参考导入位图。
mosaicking瓷砖之前,您必须确保瓷砖相邻,瓷砖边界对齐。否则,张拼接光栅的元素可以引用空间错误的位置。
马赛克的细胞位图
这个例子展示了如何导入两个相邻栅格数据文件由细胞,马赛克数据到一个光栅,并显示地图上张拼接光栅。
在这个示例中,使用的文件concord_ortho_e.tif和concord_ortho_w.tif,包含东西方相邻图像concord_ortho_e.tif的东部concord_ortho_w.tif。细胞区段和世界y限制的图片是相同的。
导入图片使用readgeoraster函数。输出eastA和westA包含图像数据的数组,输出westR和eastR是MapCellsReference对象包含空间参考信息。
[eastA, eastR] = readgeoraster (“concord_ortho_e.tif”);[westA, westR] = readgeoraster (“concord_ortho_w.tif”);
创建一个张拼接光栅结合数组。图像相邻,是由细胞组成的,所以东部边界westA与西方的边界eastA。
mosaicA = [westA eastA];
空间参考张拼接光栅通过创建一个光栅的参考对象。您可以创建一个引用对象的投影光栅细胞利用maprefcells函数。
指定的世界x- - -y光栅的限制。的x -张拼接光栅的最低限制x西方的光栅和最大的极限x限制的东部光栅。世界y张拼接光栅的限制是一样的y限制进口的位图。创建参考对象。
xlimits = [westR.XWorldLimits (1) eastR.XWorldLimits (2)];ylimits = westR.YWorldLimits;mosaicR = maprefcells (xlimits ylimits、大小(mosaicA));
数组的列使用进口readgeoraster函数从北方。因此,设置ColumnsStartFrom引用对象的属性“北”。
mosaicR。ColumnsStartFrom =“北”;
扫描线的图像显示在地图上使用mapshow函数。
mapshow (mosaicA mosaicR,“DisplayType”,“图像”)
马赛克位图的帖子
这个例子展示了如何导入两个相邻栅格数据文件的发布点,马赛克数据到一个光栅,并显示地图上张拼接光栅。
在这个示例中,使用的文件n39_w106_3arc_v2.dt1和n40_w106_3arc_v2.dt1,包含南北相邻海拔网格n39_w106_3arc_v2.dt1向南的n40_w106_3arc_v2.dt1。网格的北纬限制n39_w106_3arc_v2.dt1是一样的网格的南纬限制n40_w106_3arc_v2.dt1。发布点的间距和网格的经度范围是相同的。
导入网格使用readgeoraster函数。输出southZ和northZ是数组类型双包含高程数据,输出southR和northR是GeographicPostingsReference对象包含空间参考信息。
[southZ, southR] = readgeoraster (“n39_w106_3arc_v2.dt1”,“OutputType”,“双”);[northZ, northR] = readgeoraster (“n40_w106_3arc_v2.dt1”,“OutputType”,“双”);
创建一个张拼接光栅结合数组。发布的边界点的位图是由最外层的发布点,和进口使用数组的列readgeoraster函数从北方。因此,最南端的行northZ和北行southZ都是一样的。以避免重复行张拼接光栅,除去最南端的行northZ在结合数组。
northZ (,) = [];mosaicZ = [northZ;southZ];
空间参考张拼接光栅通过创建一个光栅的参考对象。您可以创建一个引用对象的地理栅格点使用georefpostings函数。
指定的纬度和经度范围光栅。纬度限制最低纬度的光栅南部和北部的最大纬度光栅。张拼接光栅的经度范围是一样的经度导入位图的局限性。创建参考对象。
latlim = [southR.LatitudeLimits (1) northR.LatitudeLimits (2)];lonlim = southR.LongitudeLimits;mosaicR = georefpostings (latlim lonlim、大小(mosaicZ));
设置ColumnsStartFrom和GeographicCRS引用对象的属性,这样他们进口引用对象的属性相匹配。
mosaicR。ColumnsStartFrom = southR.ColumnsStartFrom;mosaicR。GeographicCRS = southR.GeographicCRS;
显示地图上张拼接光栅。创建一个地图轴通过指定的纬度和经度范围的数据。然后,通过使用表面显示数据geoshow函数。应用colormap适合高程数据使用demcmap函数。
usamap (mosaicR.LatitudeLimits mosaicR.LongitudeLimits) geoshow (mosaicZ mosaicR,“DisplayType”,“表面”)demcmap (mosaicZ)
这个示例中使用的海拔数据来自美国地质调查局(U.S. Geological Survey)。
