% This script is use to demo the map function in the matlab
% It will use the default dataset 'geoid' and 'coast' to generate several
% maps in different projection
clc; close all; clear all

% load dataset
load geoid;
load coast;

figure(1)
% define map area
worldmap('world');

% plot the map
geoshow(geoid, geoidrefvec,'DisplayType','texturemap');
geoshow('landareas.shp', 'FaceColor', 'black');
contourm(geoid, geoidrefvec, 'b');

%%

figure(2)

worldmap('sa')
geoshow(geoid, geoidrefvec,'DisplayType','surface');
light; material dull; lightangle(90,45)
plotm(lat, long,'color','black');

%%

figure(3)
mapshow(geoid, geoidrefvec,'DisplayType','texturemap');
mapshow('landareas.shp', 'FaceColor', 'green');
mapshow(long, lat,'color','red');

%%

figure(4)
axesm sinusoid;
framem; gridm;
geoshow(lat, long);

%%
figure(5)
worldmap([-10,45],[60,160]);
geoshow('landareas.shp', 'FaceColor', 'green');
geoshow(lat, long);