Проецируйте координаты широта-долгота на x - y сопоставьте координаты
Проецируйте координаты широта-долгота на координаты 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.Nameans = "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)')
Если географическая CRS lat и lon не соответствует географической CRS proj, затем значения x и y может быть неточным. Когда proj является projcrs объект, вы можете найти его географический CRS, запросив его GeographicCRS свойство. Для примера этот код показов, как создать projcrs объект из кода EPSG 32610 и найти связанный географический CRS.
proj = projcrs(32610); proj.GeographicCRS.Name
ans =
"WGS 84"