改变技术的投影图像

使用投影存储在更改投影axesm的地图

这个例子展示了如何项目技术图像通过使用存储在一个投影方法axesm的地图。关于预测支持的更多信息金宝appaxesm的地图,看到总结和指导预测。

加载数据

阅读欧洲国家首都的纬度和经度坐标从一个文本文件。

(纬度、经度)= readvars (“european_capitals.txt”);

阅读基础图图像

指定经度和纬度限制的技术形象。找到包含大写字母的限制的地区使用geoquadline函数。扩展的限制10学位由使用纬度和经度bufgeoquad函数。

[latlim, lonlim] = geoquadline(纬度、经度);[latlim, lonlim] = bufgeoquad (latlim lonlim 10、10);

读的区域的图像“bluegreen”技术,使用的缩放级别3引用对象,数组和地图细胞在Web墨卡托投影坐标。

(RA) = readBasemapImage (“bluegreen”,latlim lonlim 3);

Unproject基础图图像

提取Web墨卡托投影坐标的基础图图像对象的引用。然后,unproject地理坐标的坐标。

[xGrid, yGrid] = worldGrid (RA);[latGrid, lonGrid] = projinv (RA.ProjectedCRS、xGrid yGrid);

欧洲各国首都的坐标已经在地理坐标。如果你是这个例子适应您自己的代码和数据投影坐标,你必须unproject地理坐标数据。

创建地图

创建一个axesm的地图投影,适合该地区使用worldmap函数。或者,您可以创建一个axesm与指定的地图投影ID使用axesm函数。

图h = worldmap (latlim lonlim);

查看存储在地图投影方法。

getm (h,“mapprojection”)

ans = ' eqdconic '

结果“eqdconic”意味着地图使用等距圆锥投影。

在地图上显示技术和欧洲国家。的geoshow功能项目和显示经度坐标的基础图形象和首都使用存储在地图投影方法。

geoshow latGrid, lonGrid geoshow(纬度,经度,DisplayType =“点”标志=“五角星形”)标题(“欧洲首都”)副标题(“技术归因:自然的地球”)

投影技术图像中的文本会出现扭曲。

改变投影使用投影CRS对象

这个例子展示了如何使用预计CRS读基础图图像对象,然后显示坐标轴上的投影图像。

加载数据

一个包含世界城市的坐标shapefile解读工作区作为地理空间表。提取的纬度和经度坐标。

GT = readgeotable (“worldcities.shp”);lat = GT.Shape.Latitude;朗= GT.Shape.Longitude;

指定使用的CRS,找到投影面积

指定投影CRS您想要使用。对于这个示例,创建一个projcrs对象为WGS 84 /等于地球亚太EPSG代码8859年。

equalEarth = projcrs (8859);

不同的投影crs不同区域的使用有效。对于许多预计CRS,你可以找到使用的面积显示CRS的著名的文本(WKT)字符串和寻找BBOX属性。这个WKT指定BBOX在倒数第二个。

wktstring (equalEarth格式=“格式化”)

ans = " PROJCRS [“WGS 84 /等于地球亚太”,BASEGEOGCRS [“WGS 84”,基准面(“1984年世界大地系统”,椭球(“WGS 84”, 6378137, 298.257223563, LENGTHUNIT[“米”,1]]],PRIMEM[“格林威治”,0,ANGLEUNIT[“学位”,0.0174532925199433]],ID (EPSG, 4326]],转换(“平等地球亚太”、方法(地球“平等”,ID (EPSG, 1078]],参数(“经度的自然起源”,150年,ANGLEUNIT[“学位”,0.0174532925199433],ID (EPSG, 8802]],参数(“假以东”,0,LENGTHUNIT(“米”,1)ID (EPSG, 8806]],参数(“假北航”,0,LENGTHUNIT(“米”,1)ID (EPSG, 8807]]], CS(笛卡尔,2)轴(“(E)”东,[1],LENGTHUNIT[“米”,1]],轴[" (N)”,[2],北LENGTHUNIT[“米”,1]],使用[[范围”等积很小的规模——Asia-Pacific-centred映射。”,“世界区域(集中在亚太地区。”),BBOX [-90, -29.99, 90, -30.01]], ID [“EPSG”, 8859]]”

指定的经度和纬度限制技术图像通过使用信息BBOX。预计CRS,打开最大经度通过添加360年度最小经度。

latlim = (-90 - 90);lonmin = -29.99;lonlim = [lonmin lonmin + 360];

阅读基础图图像

阅读的图像“卫星”技术是一个数组,Web墨卡托投影坐标映射细胞参考对象,和一个属性字符串。

(RA,鲜明)= readBasemapImage (“卫星”、latlim lonlim);

Reproject数据和基础图形象

Reproject基础图图像的坐标等于地球坐标。

[xGrid, yGrid] = worldGrid (RA);[latGrid, lonGrid] = projinv (RA.ProjectedCRS、xGrid yGrid);[xEqualEarth, yEqualEarth] = projfwd (equalEarth、latGrid lonGrid);

项目世界城市的地理坐标等于地球坐标。如果你是这个例子适应您自己的代码和数据投影坐标,你必须首先unproject地理坐标数据。

[xData, yData] = projfwd (equalEarth、纬度、经度);

创建地图

显示投影技术形象和世界城市在地图上没有轴标签。

图mapshow (xEqualEarth, yEqualEarth)在mapshow (xData yData DisplayType =“点”标志=“。”)标题(“世界城市”)副标题(基础图归因:“+鲜明)轴从

投影技术图像中的文本会出现扭曲。