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

projfwd

将经纬度坐标投影至 x-y 地图坐标系

全页折叠

语法

[x,y] = projfwd(proj,lat,lon)

说明

[x,y] = projfwd(proj,lat,lon) 将由 lat 和 lon 指定的经纬度坐标转换为 x 和 y 地图坐标，该坐标系采用由 proj 指定的投影坐标参考系。使用 projcrs 对象 (自 R2020b 起)、地图投影结构或 GeoTIFF 信息结构来指定 proj。

示例

示例

全部折叠

打开实时脚本

通过指定地图投影，将经纬度坐标投影至 x-y 坐标系。然后，将投影后的坐标显示在地图上。

为此，首先需指定波士顿地标的经纬度坐标。请使用 NAD83 地理坐标参考系指定坐标。

lat = [42.3604 42.3691 42.3469 42.3480 42.3612];
lon = [-71.0580 -71.0710 -71.0623 -71.0968 -71.0941];

然后，将波士顿的 GeoTIFF 图像导入为数组和地图单元参考对象。通过查询参考对象的 ProjectedCRS 属性获取地图投影信息。验证投影 CRS 所基于的地理坐标系是否为 NAD83。

[A,R] = readgeoraster("boston.tif");
proj = R.ProjectedCRS;
proj.GeographicCRS.Name
ans = 
"NAD83"

将经纬度坐标投影到 x-y 坐标系，使用与 GeoTIFF 图像相同的投影坐标参考系 (CRS)。

[x,y] = projfwd(proj,lat,lon);

在同一张地图上同时显示 GeoTIFF 图像和投影坐标。更改坐标的标记符号和颜色，使其更醒目。然后，添加坐标轴标签。

mapshow(A,R)
mapshow(x,y,DisplayType="point",Marker="o", ...
    MarkerFaceColor="y",MarkerEdgeColor="none")
xlabel("x (survey feet)")
ylabel("y (survey feet)")

输入参数

全部折叠

地图投影，指定为 projcrs 对象 (自 R2020b 起)、标量地图投影结构 (mstruct) 或 GeoTIFF 信息结构。有关地图投影结构的更多信息，请参阅 defaultm。有关 GeoTIFF 信息结构的更多信息，请参阅 geotiffinfo。

数据类型: struct

大地纬度，指定为标量值、向量、矩阵或 N 维数组，单位为度。lat 和 lon 的大小必须匹配。

数据类型: single | double

大地经度，指定为标量值、向量、矩阵或 N 维数组，单位为度。lat 和 lon 的大小必须匹配。

数据类型: single | double

输出参量

全部折叠

投影后的 x 坐标，以标量值、向量、矩阵或 N 维数组形式返回。

投影后的 y 坐标，以标量值、向量、矩阵或 N 维数组形式返回。

提示

如果 lat 和 lon 的地理坐标参考系与 proj 的地理坐标参考系不匹配，则 x 和 y 的数值可能不准确。当 proj 是 projcrs 对象时，可通过查询其 GeographicCRS 属性获取其地理坐标系。例如，此代码演示了如何根据 EPSG 代码 32610 创建一个 projcrs 对象，并查找其关联的地理坐标参考系。

proj = projcrs(32610);
proj.GeographicCRS.Name
ans = 

    "WGS 84"

版本历史记录

在 R2006a 之前推出

全部展开

另请参阅

函数

对象

主题