项目和显示栅格数据
项目或unproject均匀间隔的栅格地图与地理相关的数据或参考对象,您必须首先创建一个坐标网格光栅的大小相匹配。使用不同的网格创建函数根据投射的方式。当你项目经度坐标x- - - - - -y坐标,创建一个网格使用geographicGrid函数。当你unprojectx- - - - - -y坐标经度坐标,创建一个网格使用worldGrid函数。
将栅格数据转换后,可以显示在地图上使用可视化等功能mapshow和geoshow。使用mapshow为预测x- - - - - -y坐标和geoshow对于非预料的经度坐标。
项目栅格数据
项目与地理相关的数据栅格参考对象,首先创建一个网格光栅中的每个点的经度坐标。然后,项目的地理坐标x- - - - - -y地图坐标。
例如,导入高程栅格数据作为数组和一个地理单元的参考对象。得到每个点的经度坐标光栅使用geographicGrid函数。
(Z, R) = readgeoraster (“n39_w106_3arc_v2.dt1”);(纬度、经度)= geographicGrid (R);
现在,你有你的网格,选择地图投影投影坐标时使用。对于这个示例,创建一个projcrs对象UTM区13在北半球。然后,项目的经度坐标x- - - - - -y坐标。
p = projcrs (32613);(x, y) = projfwd (p,纬度,经度);
显示投影光栅表面通过调用mapshow和指定的x- - - - - -y坐标和高程数组。添加轴标签和应用colormap适合高程数据。
图mapshow (x, y, Z,“DisplayType”,“表面”)包含(“x(米))ylabel (“y”(米))demcmap (Z)
如果地理经度坐标不匹配的CRS地理CRS的投影CRS,那么投影坐标可能是不准确的。你可以找到的CRS地理projcrs对象或地理栅格引用对象通过查询它们GeographicCRS属性。
p.GeographicCRS.Name
ans = " WGS 84”
R.GeographicCRS.Name
ans = " WGS 84”
DTED文件中使用这个示例是由美国地质调查局提供。
Unproject栅格数据
unproject数据与地图有关光栅参考对象,首先创建一个网格x- - - - - -y每个点的坐标光栅。然后,unprojectx- - - - - -y坐标映射到地理坐标。
例如,波士顿的导入一个图像和地图细胞引用对象数组。地图投影的信息projcrs对象通过查询ProjectedCRS属性的引用对象。
(Z, R) = readgeoraster (“boston.tif”);p = R.ProjectedCRS;
得到了x- - - - - -y每个点的坐标光栅使用worldGrid函数。
(x, y) = worldGrid (R);
Unproject的x- - - - - -y坐标通过使用经度坐标projinv功能和指定projcrs对象和协调网格。
(纬度、经度)= projinv (p, x, y);
通过调用显示未计划的形象geoshow并指定经度坐标和图像数组。默认情况下,geoshow使用板Carree投影显示坐标。然后,添加轴标签。
图geoshow(纬度,经度,Z)包含(的经度(度))ylabel (的纬度(度))
