我有一个带有横向墨卡托投影的shapefile。它还在ArcGIS中具有与其关联的错误北向/东向值。
我试图将多边形映射到欧洲地图上,但它们没有出现在它上面。当我使用worldmap world
时,多边形出现在非洲,即使它们应该在爱尔兰。
如何在正确的位置绘制多边形?他们有与他们相关的纬度和经度,所以我认为这必须是预测。
S = shaperead(polygons.shp','UseGeoCoords', true)
h = axesm('tranmerc','FalseEasting',60000,'FalseNorthing',750000,'Grid','on','Frame','on', 'MlabelParallel',0,'PlabelMeridian',0)
当我试图描绘纬度和经度时,它们最终会出现在非洲:
h = plot([S.Lat],[S.Lon], 'k:');
或者当我尝试使用geoshow
时,他们会用颜色填充地图:
geoshow(S,'FaceColor',[1 1 .5],'EdgeColor',[.6 .6 .6]);