项目光栅数据

要投影与地理栅格引用对象相关联的数据，首先为栅格中的每个点创建一个纬度-经度坐标网格。然后，将地理坐标投影到x-y地图坐标。

例如，导入高程栅格数据作为数组和地理单元格引用对象。方法获取栅格中每个点的经纬度坐标geographicGrid函数。

[Z,R] = readgeoraster(“n39_w106_3arc_v2.dt1”）;[lat,lon] = geoticgrid (R);

有了网格之后，选择投影坐标时要使用的映射投影。对于本例，创建一个projcrs为北半球UTM第13区观测。然后，将经纬度坐标投影到x-y坐标。

P = projects (32613);[x,y] = projfwd(p,lat,lon);

通过调用将投影光栅显示为曲面mapshow并指定x-y坐标和仰角数组。添加轴标签并应用适合高程数据的色图。

图mapshow (x, y, Z,“DisplayType”，“表面”)包含(“x(米)) ylabel (“y”(米)) demcmap (Z)

如果纬度-经度坐标的地理CRS与投影CRS的地理CRS不匹配，则投影坐标可能不准确。您可以找到一个的地理CRSprojcrs对象或地理光栅引用对象GeographicCRS属性。

p.GeographicCRS.Name

ans = "WGS 84"

R.GeographicCRS.Name

ans = "WGS 84"

本例中使用的DTED文件由美国地质调查局提供。

取消光栅数据投影

要取消与映射光栅引用对象关联的数据的投影，首先要创建的网格x-y栅格中每个点的坐标。然后，取消投影x-y将坐标映射为地理坐标。

例如，导入波士顿的图像作为数组和映射单元格引用对象。获取有关映射投影的信息projcrs对象。ProjectedCRS引用对象的属性。

[Z,R] = readgeoraster(“boston.tif”）;p = R.ProjectedCRS;

得到了x-y控件为光栅中每个点的坐标worldGrid函数。

[x,y] = worldGrid(R);

Unproject的x-y方法将坐标转换为经纬度坐标projinv函数，并指定projcrs对象和坐标网格。

[lat,lon] = projinv(p,x,y);

通过调用显示未投影的图像geoshow并指定经纬度坐标和图像数组。默认情况下,geoshow使用Plate Carrée投影显示坐标。然后，添加轴标签。

图geoshow(lat,lon,Z) xlabel(的经度(度)) ylabel (的纬度(度)）