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

geographicToDiscrete

将地理坐标转换为离散下标

全页折叠

语法

[I,J] = geographicToDiscrete(R,lat,lon)

说明

[I,J] = geographicToDiscrete(R,lat,lon) 根据栅格参考对象 R 定义的关系，根据地理坐标 lat 和 lon 计算离散行下标 I 和 J。

对 I 和 J 的解释取决于 R 的 RasterInterpretation 属性。

  • 'cells' — I 和 J 分别是栅格像元或图像像素的行下标和列下标。

  • 'postings' - I 和 J 分别是最近的发布点样本的行下标和列下标。

示例

示例

全部折叠

打开实时脚本

将科罗拉多州某地区的 elevation 数据 [1] 作为数组和 GeographicPostingsReference 对象导入。找出距离南博尔德峰纬度和经度最近的发布点的索引，该峰的纬度和经度分别为 39.9539 度和-105.2992 度。

[Z,R] = readgeoraster('n39_w106_3arc_v2.dt1');
[I,J] = geographicToDiscrete(R,39.9539,-105.2992)
I = 
56

J = 
842

该结果表明，最接近该地理坐标的发布点位于数组的第 56 行第 842 列。

[1] 本示例中使用的高程数据来自美国地质调查局。

输入参数

全部折叠

栅格参考，指定为 GeographicCellsReference 或 GeographicPostingsReference 对象。

纬度坐标，指定为数值数组。lat 的大小必须与 lon 的大小相匹配。

数据类型: single | double

经度坐标，指定为一个数值数组。lon 的大小必须与 lat 的大小相匹配。

数据类型: single | double

输出参量

全部折叠

行索引，以整数数组的形式返回。I 与 lat 大小相同。

对于一个 m-by-n 的栅格，且地理坐标位于栅格边界内时，I 的值位于 [1, m] 范围内。

如果第 k 个输入坐标 (lat(k),lon(k)) 超出了地理坐标系中的栅格边界，则该函数将相应的下标 I(k) 和 J(k) 设置为 NaN 的值。

数据类型: double

列索引，以整数数组的形式返回。J 与 lon 的大小相同。

对于一个 m-by-n 的栅格，且地理坐标位于栅格边界内时，J 的值范围为 [1, n]。

如果第 k 个输入坐标 (lat(k),lon(k)) 超出了地理坐标系中的栅格边界，则该函数将相应的下标 I(k) 和 J(k) 设置为 NaN 的值。

数据类型: double

版本历史记录

在 R2013b 中推出

另请参阅

函数

对象