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

worldGrid

栅格元素的世界坐标

全页折叠

语法

[X,Y] = worldGrid(R)
[X,Y] = worldGrid(R,gridOption)

说明

[X,Y] = worldGrid(R) 将栅格元素的世界坐标作为二维数组 X 和 Y 返回。栅格元素 (i,j) 的坐标为 (X(i,j),Y(i,j))。

示例

[X,Y] = worldGrid(R,gridOption)，其中 gridOption 即 "gridvectors"，返回 X 和 Y 作为行向量。栅格元素 (i,j) 的坐标为 (X(j),Y(i))。gridOption 的默认值为 "fullgrid"，其返回值为二维数组 X 和 Y。

示例

示例

全部折叠

打开实时脚本

导入高程数据，查找导入数据的坐标，然后将数据以曲面形式显示出来。

首先将华盛顿山周边区域的高程数据 [1] 作为数组和地图像元引用对象导入。通过将输出类型指定为 "double" 来准备绘图数据。

[Z,R] = readgeoraster("MtWashington-ft.grd",OutputType="double");

获取数组中每个元素的坐标。

[X,Y] = worldGrid(R);

使用适当的颜色图，将高程数据以表面形式显示出来。将 EdgeColor 属性设置为 none"，即可将网格从曲面上移除。

surf(X,Y,Z,EdgeColor="none")
demcmap(Z)
xlabel("x (meters)")
ylabel("y (meters)")
zlabel("Elevation (feet)")

Figure contains an axes object. The axes object with xlabel x (meters), ylabel y (meters) contains an object of type surface.

[1] The elevation data used in this example is courtesy of the U.S. Geological Survey.

打开实时脚本

为一个 3×4 的栅格创建地图像元引用对象，其中 x 值的范围为 [7000, 7400] 米，y 值的范围为 [2700, 3300] 米。获取栅格元素的坐标，并将其作为行向量返回。

R = maprefcells([7000 7400],[2700 3300],[3 4]);
[X,Y] = worldGrid(R,'gridvectors')
X = 1×4

        7050        7150        7250        7350


Y = 1×3

        2800        3000        3200

若未将第二个参数指定为 'gridvectors'，则 worldGrid 函数默认返回二维数组。

[xFull,yFull] = worldGrid(R)
xFull = 3×4

        7050        7150        7250        7350
        7050        7150        7250        7350
        7050        7150        7250        7350


yFull = 3×4

        2800        2800        2800        2800
        3000        3000        3000        3000
        3200        3200        3200        3200

输入参数

全部折叠

栅格参考，指定为 MapCellsReference 或 MapPostingsReference 对象。

若 R 是 MapCellsReference 对象，则 X 和 Y 均为单元中心。如果 R 是 MapPostingsReference 对象，则 X 和 Y 是发布点。

网格选项，指定为以下值之一：

  • "fullgrid" — 将 X 和 Y 作为二维数组返回，其中 X 的每一行内容相同，Y 的每一列内容相同。这是默认行为。

  • "gridvectors" - 将 X 和 Y 作为行向量返回。当您希望减少内存使用且无需二维数组时，请使用此选项，例如使用 surf 函数绘制大型数据集时。只有当 R 的 TransformationType 属性值为 "rectilinear" 时，才能将 gridOption 指定为 "gridvectors"。

此表显示了 "fullgrid" 与 "gridvectors" 之间的差异。

"fullgrid""gridvectors"
R = maprefcells([7000 7400],[2700 3300],[3 4]);
[X,Y] = worldGrid(R)
X =

        7050        7150        7250        7350
        7050        7150        7250        7350
        7050        7150        7250        7350


Y =

        2800        2800        2800        2800
        3000        3000        3000        3000
        3200        3200        3200        3200
 
R = maprefcells([7000 7400],[2700 3300],[3 4]);
[X,Y] = worldGrid(R,"gridvectors")
X =

        7050        7150        7250        7350


Y =

        2800        3000        3200

数据类型: char | string

输出参量

全部折叠

x-坐标，以二维数组或行向量的形式返回。默认情况下，X 是二维数组。要将 X 作为行向量返回，请将 gridOption 指定为 "gridvectors"。

默认情况下，当 gridOption 为 "fullgrid" 时，X 和 Y 的大小各等于 R 的 RasterSize 属性。当 gridOption 为 "gridvectors" 时，X 和 Y 的长度分别等于 R 的 RasterSize 属性的第二个和第一个元素。

y-坐标，以二维数组或行向量的形式返回。默认情况下，Y 是二维数组。要将 Y 作为行向量返回，请将 gridOption 指定为 "gridvectors"。

默认情况下，当 gridOption 为 "fullgrid" 时，X 和 Y 的大小各等于 R 的 RasterSize 属性。当 gridOption 为 "gridvectors" 时，X 和 Y 的长度分别等于 R 的 RasterSize 属性的第二个和第一个元素。

版本历史记录

在 R2021a 中推出

另请参阅

函数

对象