主要内容

本页采用了机器翻译。点击此处可查看英文原文。

mapcrop

裁剪平面地图栅格

说明

[B,RB] = mapcrop(A,RA,xlimits,ylimits) 对使用栅格引用 RA 的平面地图栅格 A 进行裁剪,并返回裁剪后的栅格 B 和栅格引用 RB。返回的栅格被裁剪为世界坐标系中的边界,该边界与 xlimitsylimits 指定的边界接近。

若要裁剪由点、线或多边形形状表示的、位于平面坐标系中的向量数据,请改用 mapclip 函数。

示例

示例

全部折叠

裁剪投影栅格图像,并显示裁剪后的图像。

首先,导入一张波士顿的投影栅格图像和一个地图像元参考对象。然后,将图像裁剪至 xlimitsylimits 指定的边界。

[A,RA] = readgeoraster("boston.tif");
xlimits = [771660  773290];
ylimits = [2953410 2955240];
[B,RB] = mapcrop(A,RA,xlimits,ylimits);

显示裁剪后的图像。

mapshow(B,RB)

输入参数

全部折叠

平面地图栅格数据,指定为 M × NM × N × P 格式的数值数组或逻辑数组。

A 的栅格参考,指定为 MapCellsReference 对象或 MapPostingsReference 对象。RATransformationType 属性必须为 'rectilinear'

x 的最小值和最大值限制,指定形式为 [xmin xmax] 的两个元素的数值向量,其中 xmax 大于 xmin

y 的最小值和最大值限制,指定形式为 [ymin ymax] 的两个元素的数值向量,其中 ymax 大于 ymin

输出参量

全部折叠

裁剪后的平面地图栅格数据,以数值数组或逻辑数组的形式返回。B 的数据类型和大小与 A 的数据类型和大小一致。

如果由 xlimitsylimits 指定的边界与由 ARA 指定的栅格没有交集,则 B 为空。

B 的栅格参考,以 MapCellsReference 对象或 MapPostingsReference 对象的形式返回。RB 的对象类型与 RA 的对象类型一致。

RB 的确切 xy 边界与 xlimitsylimits 指定的边界不一致,除非它们与单元格边界或发布位置重合。否则,RB 的极限值会比 xlimitsylimits 稍大一些。

如果由 xlimitsylimits 指定的边界与由 ARA 指定的栅格没有交集,则 RB 为空。

替代功能

  • 要在地理坐标系中裁剪栅格,请使用 geocrop 函数。

  • 要创建一个边界与 xlimitsylimits 匹配的栅格,请使用 mapresample 函数对该栅格进行重采样。

版本历史记录

在 R2020a 中推出