В этом примере показано, как составить карту крупных городов мира с улучшенными береговыми линиями и рельефом местности. В нём используется модифицированная азимутальная проекция карты Бриземейстера. Пример включает в себя дополнительный код, который позволяет пользователю в интерактивном режиме выбирать местоположение и получать имя и местоположение ближайшего города. Чтобы увидеть эту часть примера, необходимо выполнить полный пример, вывести последнюю иллюстрацию на отдельный рисунок MATLAB, а затем запустить дополнительный код в командной строке MATLAB.
Создайте объект оси карты.
figure axesm bries text(2.8,-1.8,'Briesemeister projection','HorizontalAlignment','right') framem('FLineWidth',1)
Загрузите растровые данные отметки и объект ссылки на географические ячейки. Отображение данных на карте.
load topo60c geoshow(topo60c,topo60cR,'DisplayType','texturemap')
Примените карту цветов, соответствующую данным отметки. Сделайте дисплей ярче.
demcmap(topo60c) brighten(0.5)
Загрузка координат глобальной береговой линии. Обобщить береговые линии до 0,25-градусного допуска. Затем обведите береговые линии коричневым цветом.
load coastlines [rlat,rlon] = reducem(coastlat,coastlon,0.25); geoshow(rlat,rlon,'Color',[.6 .5 .2],'LineWidth',1.5)
Прочитайте файл формы, содержащий названия городов по всему миру и их координаты в широте и долготе.
cities = shaperead('worldcities','UseGeoCoords',true);
Извлеките точечные широты и долготы с помощью extractfieldи добавьте их на карту.
lats = extractfield(cities,'Lat'); lons = extractfield(cities,'Lon'); geoshow(lats, lons,... 'DisplayType', 'point',... 'Marker', 'o',... 'MarkerEdgeColor', 'r',... 'MarkerFaceColor', 'r',... 'MarkerSize', 3) text(-2.8,-1.8,'Major World Cities')
Теперь, используя созданную карту, можно настроить простой цикл для запроса щелчков на карте и отображения имени и координат ближайшего города. Необходимо вставить последнюю карту, созданную на шаге 4, в отдельное окно фигуры MATLAB с помощью кнопки, которая отображается в верхней части карты. Также, в следующем коде, установить runCitySelectionLoop кому trueи выполните код в командной строке.
Код сначала отображает текстовые инструкции в верхней левой части карты. Затем он входит в цикл, в котором он захватывает выбранные широты и долготы с inputm. Использовать distance для вычисления большого расстояния по окружности между каждой выбранной точкой и каждым городом в базе данных. Определите индекс ближайшего города, измените внешний вид его маркерного символа и отобразите имя города и координаты широты/долготы.
runCitySelectionLoop = false; % Set to true to run optional city selection loop if(runCitySelectionLoop) h1 = text(-2.8, 1.9, 'Click on a dot for its city name. Press ENTER to stop'); h2 = text(-2.8, 1.7, ''); h3 = text(-2.8, 1.5, 'City Coordinates.'); while true [selected_lat,selected_lon] = inputm(1); if isempty(selected_lat) break % User typed ENTER end d = distance(lats, lons, selected_lat, selected_lon); k = find(d == min(d(:)),1); city = cities(k); geoshow(city.Lat, city.Lon, ... 'DisplayType', 'point', ... 'Marker', 'o', ... 'MarkerEdgeColor', 'k', ... 'MarkerFaceColor', 'y', ... 'MarkerSize', 3) h2.String = city.Name; h3.String = num2str([city.Lat, city.Lon],'%10.2f'); end disp('End of input.') end