filename=gunzip('gshhs_c.b.gz',tempdir);
world=gshhs(filename{1});
figure
worldmap world
geoshow([world.Lat],[world.Lon])
levels=[world.Level];
land=(levels==1);
lake=(levels==2);
island=(levels==3);
figure
worldmap world
geoshow(world(land),'facecolor',[0 1 0])
geoshow(world(lake),'facecolor',[0 0 1])
geoshow(world(island),'facecolor',[1 1 0])
figure
worldmap Italy
projection=gcm;
latlim=projection.maplatlimit;
lonlim=projection.maplonlimit;
italy=gshhs(filename{1},latlim,lonlim);
[,ix]=sort([italy.Level],'descend');
italy=italy(ix);
geoshow(italy,'Facecolor','green')
setm(gca,FFaceColor','cyan')