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

Для этого сначала задайте координаты широты и долготы ориентиров в Бостоне. Задайте координаты в географическом CRS NAD83.

lat = [42.3604 42.3691 42.3469 42.3480 42.3612];
lon = [-71.0580 -71.0710 -71.0623 -71.0968 -71.0941];

Затем импортируйте изображение GeoTIFF Бостона как массив и объект ссылки ячеек карты. Получите информацию о проекции карты путем запроса ProjectedCRS свойство ссылочного объекта. Проверьте, что географический CRS, лежащий в основе спроектированного CRS, является NAD83.

[A,R] = readgeoraster('boston.tif');
proj = R.ProjectedCRS;
proj.GeographicCRS.Name

ans = 
"NAD83"

Спроектируйте координаты долготы широты к координатам x-y с помощью того же спроектированного CRS в качестве изображения GeoTIFF.

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

Отобразите изображение GeoTIFF и спроектированные координаты на той же карте. Измените символ маркера и цвет координат, таким образом, они больше отображаются. Затем добавьте подписи по осям.

mapshow(A,R)
mapshow(x,y,'DisplayType','point','Marker','o', ...
    'MarkerFaceColor','y','MarkerEdgeColor','none')
xlabel('x (survey feet)')
ylabel('y (survey feet)')