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

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

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)')