projfwd

使用Proj.4地图投影库

描述

例子

[[X,,,,y] = projfwd(Proj,,,,拉特,,,,返回Xy从正向投影转换中绘制坐标。Proj是定义地图投影的结构。Proj可以是MS结构或Geotiff信息结构体。拉特是纬度和经度坐标的阵列。

有关Geotiff的完整列表信息以及您可以与之一起使用的图投影结构Projinv,请参阅参考页Projlist

例子

全部收缩

在波士顿的正骨上显示马萨诸塞州的状态边界。

首先,读取马萨诸塞州状态边界的地理矢量数据。

s = shaperead('usastatehi',,,,'usegeocoords',真的,...“选择器”,{@(name)strcmpi(名称,“马萨诸塞州”),'姓名'});

获取有关正常光的投影信息Geotiffinfo。然后,投影状态边界向量。

proj = geotiffinfo('boston.tif');[x,y] = projfwd(proj,s.lat,s.lon);

阅读并显示波士顿的正射图像。

[a,r] = readgeoraster('boston.tif');MapShow(A,R)Xlabel(“马地国家飞机升级,调查脚”)ylabel(“马地国家飞机北方,调查脚”

显示状态边界。通过设置地图限制来显示更多详细信息。

mapshow(x,y,'颜色',,,,'K',,,,'行宽',2)Xlim([645000 895000]);Ylim([2865000 3040000]);

输入参数

全部收缩

地图投影,指定为标量图投影结构(MSTRUCT)或Geotiff信息结构。

数据类型:结构

一个或多个点的大地纬度(以标量值,矢量,矩阵或n-d数组为单位)为单位。大小必须与大小相匹配输入。

数据类型:单身的|双倍的

一个或多个点的纵向,以标量值,向量,矩阵或n-d数组为单位。大小必须与大小相匹配拉特输入。

数据类型:单身的|双倍的

输出参数

全部收缩

投影坐标系中一个或多个点的X坐标,作为标量值,向量,矩阵或N-D数组返回。

投影坐标系中一个或多个点的y坐标坐标值,作为标量值,向量,矩阵或n-D数组返回。

在R2006a之前引入