%% Interactive Global City Finder 
% This demo constructs a map of major world cities enhanced with coastlines
% and terrain. It uses the modified azimuthal Briesemeister map projection.
% Then it allows a user to interactively pick a location and automatically
% highlights the nearest city, displaying its name and coordinates.

% Copyright 2003-2004 The MathWorks, Inc.
% $Revision: 1.1.6.2 $  $Date: 2004/12/18 07:42:19 $

%% Step 1: Set up a map axes object and render a global elevation grid
%
% Construct the axes.
figure; axesm bries
text(1.2, -1.8, 'Briesemeister projection', 'FontSize', 8)
framem('FLineWidth',1)

%%
% Load and display a 1-by-1-degree elevation grid.
load topo
geoshow(topo, topolegend, 'DisplayType', 'texturemap')

%% Step 2: Improve the terrain display
%
% Get a colormap appropriate for elevation.
demcmap(topo)

%%
% Make it brighter.
brighten(.5)

%% Step 3: Add simplified coastlines
%
% Load global coastline coordinates.
load coast

%%
% Generalize the coastlines to 0.25-degree tolerance.
[rlat, rlon] = reducem(lat, long, 0.25);

%%
% Plot the coastlines in brown.
geoshow(rlat, rlon, 'Color', [.6 .5 .2], 'LineWidth', 1.5)

%% Step 4: Plot city locations with red point markers
%
% Read a shapefile containing names of cities worldwide and their
% coordinates in latitude and longitude.
cities = shaperead('worldcities', 'UseGeoCoords', true);

%%
% Extract the point latitudes and longitudes with |extractfield|, and add
% them to the map.
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', 'FontSize',8)

%% Step 5: Loop forever while selecting cities
%
% Prompt the user to click on or near a city, capturing the selected
% latitude and longitude with |inputm|.  Use |distance| to calculate the
% great circle distance betweeen the selected point and every city in the
% database.  Determine index of the closest city, change the appearance of
% its marker symbol, and display the city's name and coordinates.

h1 = text(-2.8, 1.7, 'Click on a dot for its city name');
h2 = text(-2.8, 1.5, 'Press ENTER to stop');
while true
    [selected_lat,selected_lon] = inputm(1);
    if isempty(selected_lat) 
        break % User typed ENTER
    end
    d = distance(lats, lons, repmat(selected_lat,size(lats)),...
                             repmat(selected_lon,size(lons)));
    k = find(d == min(d(:)));
    where = [cities(k).Lat, cities(k).Lon];
    plotm(where, 'o',...
          'MarkerEdgeColor', 'k',...
          'MarkerFaceColor', 'y',...
          'MarkerSize', 3)
    delete(h1,h2)
    h1 = text(-2.8, 1.7, cities(k(1)).Name);
    h2 = text(-2.8, 1.5, num2str(where,'%10.2f'));
end
disp('End of input.')