投射CRS坐标变换到另一个

打开生活的脚本

如果你直接比较数据集与不同的投影坐标参考系统(crs),那么结果是不准确的。因此,在比较数据集,首先验证crs是相同的。如果不同的投影CRS有相同的底层地理CRS,然后你可以转换从一个投影坐标CRS到另一个。一旦数据集引用相同的投影CRS,你可以进行比较。

改变计划x- - - - - -y坐标不同投影CRS,首先unprojectx- - - - - -y坐标通过使用经度坐标projinv函数。然后,项目的经度坐标x- - - - - -y坐标投影CRS通过使用不同projfwd函数。

例如,包含导入shapefilex- - -y坐标在波士顿的道路。也导入shapefile信息结构。找到坐标的投影CRS访问CoordinateReferenceSystem场的结构。

s = shaperead (“boston_roads.shp”);x1 = [s.X];日元= [s.Y];信息= shapeinfo (“boston_roads.shp”);p1 = info.CoordinateReferenceSystem;

Unproject的x- - - - - -y坐标和返回经度坐标。

(纬度、经度)= projinv (p1 (x1, y1);

选择一个新的投影CRS目标投影。对于这个示例,创建一个projcrs对象UTM区19 n。确认两个投影CRS地理CRS相同。如果地理crs是不同的,那么投影坐标可能是不准确的。你可以找到地理CRS通过查询GeographicCRS财产的projcrs对象。

p2 = projcrs (26919);p2.GeographicCRS.Name

ans = " NAD83 "

p1.GeographicCRS.Name

ans = " NAD83 "

项目的经度坐标x- - - - - -y通过指定的坐标projcrs对象创建。

(x2, y2) = projfwd (p2,纬度,经度);

比较原始的x- - - - - -y坐标与新x- - - - - -y坐标显示它们。每个图添加标签和标题。

图mapshow (x1, y1)包含(“x(米))ylabel (“y”(米)标题(p1.Name)

图包含一个坐标轴对象。坐标轴对象与标题NAD83 /麻萨诸塞州大陆,包含x(米),ylabel y(米)包含一个类型的对象。

图mapshow (x2, y2)包含(“x(米))ylabel (“y”(米)标题(p2.Name)

图包含一个坐标轴对象。坐标轴对象标题NAD83 / UTM区19 n,包含x(米),ylabel y(米)包含一个类型的对象。

可视化是相似的,但沿轴坐标显示统治者对应不同的投影crs。

另请参阅