filename=gunzip('gshhs_c.b.gz',tempdir);
worldmap australia
projection=gcm;
latlim=projection.maplatlimit;
lonlim=projection.maplonlimit;
australia=gshhs(filename{1},latlim,lonlim);
geoshow(australia,'Facecolor','green');
levels=[australia.Level];
unique(levels)
L1 = australia(levels == 1);
L2 = australia(levels == 2);
figure
axesm('mercator', 'MapLatLimit', latlim, 'MapLonLimit', lonlim)
gridm; mlabel; plabel
geoshow([L1.Lat], [L1.Lon], 'Color', 'blue')
geoshow([L2.Lat], [L2.Lon], 'Color', 'red')
tightmap

