Депроектируйте координаты x-y в координаты широта-долгота путем определения проектируемой CRS координат x-y. Затем отобразите координаты широта-долгота на географических осях.

Для этого сначала импортируйте файл shapefile, содержащий координаты x и y дорог в Concord, MA. Получите информацию о файле shapefile как о структуре. Найдите спроецированную CRS для координат путем доступа к CoordinateReferenceSystem поле структуры.

roads = shaperead('concord_roads.shp');
x = [roads.X];
y = [roads.Y];
info = shapeinfo('concord_roads.shp');
proj = info.CoordinateReferenceSystem;

Отключите координаты x-y в координаты широта-долгота.

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

Отображение координат на географических осях.

figure
geoplot(lat,lon)
hold on
geobasemap('streets')

Географическая CRS координат x-y, используемых в этом примере, NAD83. Вы можете найти географическую CRS, которая лежит в основе спроецированной CRS, запросив GeographicCRS свойство.

proj.GeographicCRS.Name

ans = 
"NAD83"

Географическая CRS, лежащая в основе 'streets' basemap WGS84. NAD83 и WGS84 похожи, но не идентичны. Поэтому на высоких уровнях масштаба координаты и основная карта могут показаться несовпадающими.