This demo illustrates using WORLDDATAMAP to draw maps of the world using data from a variety of shapefiles and regular data grids.
The following examples draw maps of the World using data from landareas.shp.
Draw a map of the World with coastlines.
figure worlddatamap world title('World')
Draw a map of Antarctica as a patch.
figure worlddatamap('antarctica','patch') title('Antartica')
Draw a map of Africa and India and major lakes and rivers.
figure
worlddatamap({'Africa','India'},{'lakes', 'rivers'})
title('Africa and India with Major Lakes and Rivers')
Draw a map of New Zealand's major cities.
figure worlddatamap('New Zealand','citiesonly') title('Cities of New Zealand')
Draw a map of the geoid surface.
figure load geoid worlddatamap(geoid, geoidrefvec) title('Geoid Surface')
Draw a map of terrain elevations in Korea.
load korea figure worlddatamap(map, refvec, 'dem'); title('Terrain Elevation Map of Korea')
Draw a line map of the region of Korea.
load korea figure worlddatamap(map, refvec, 'line'); title('Line Map of Korea')
Mask the World's land surface.
figure worlddatamap worldonly patch h=handlem('world'); set(h,'facecolor','black') title('Land Areas of the World')
Display the geoid and mask the land surface.
figure worlddatamap worldonly patch h=handlem('world'); set(h,'facecolor','black') load geoid geoshow(geoid, geoidrefvec, 'display', 'texture') zdatam(h,max(geoid(:))) title('Geoid Surface with Masked Land Areas')
The following examples draw a map of the United States using data from landareas.shp or usastatehi.shp.
Draw the United States and add rivers and lakes.
figure worlddatamap('usa',{'rivers','lakes'}); title('USA Rivers Lakes')
Draw a map of the United States using the usastatehi shapefile.
figure worlddatamap usastatehi.shp usa title('USA States')
Draw a map of the United States as a patch.
figure worlddatamap usastatehi.shp usa patch title('USA Patch States');
Draw a map of the states near Mexico with labels.
figure worlddatamap usastatehi.shp mexico patch labels title('States near Mexico with Labels')
Draw a map of the states near Mexico as patches only (notice that labels are not displayed).
figure worlddatamap usastatehi.shp mexico patchonly labels title('States near Mexico as Patch Only')
Draw a map of the state of Missouri with a state label.
figure
worlddatamap({'usastatehi.shp','name'},'missouri','patch','labels')
title('Missouri')
Draw a map of only the state of Missouri with all features.
figure
worlddatamap({'usastatehi.shp','name'},'missourionly','patch','all')
title('Missouri Only with Major Rivers, Lakes, and Cities')
A variety of world vector data in shapefile format can be downloaded from the Internet.
These examples download and use data from three different world vector shapefiles, cntry98, cntry02, and world_borders.
For details on locating other geospatial data for download over the Internet, see the following documentation at The MathWorks web site: http://www.mathworks.com/support/tech-notes/2100/2101.html
Download the cntry98 shapefiles The cntry98 shapefiles can be downloaded from:
http://www.zmuc.dk/public/GIS-intro/GIS-data/Lektion1/Temaer/
- or -
http://giscenter.isu.edu/spatial_library/Other/World/data/world/
url = 'http://www.zmuc.dk/public/GIS-intro/GIS-data/Lektion1/Temaer/'; try urlwrite([url '/cntry98.shp'],'cntry98.shp'); urlwrite([url '/cntry98.dbf'],'cntry98.dbf'); urlwrite([url '/cntry98.shx'],'cntry98.shx'); catch error('Unable to download cntry98 shapefiles.'); end
Download the cntry02 shapefiles. The cntry02 shapefiles can be downloaded from:
http://openmap.bbn.com/data/shape/cntry02.tar.gz
- or -
http://www.zmuc.dk/public/GIS-intro/GIS-data/Lektion1/Temaer/
url = 'http://openmap.bbn.com/data/shape/cntry02.tar.gz'; try files = untar(url) catch error('Unable to download %s.', url); end
files =
Columns 1 through 4
'cntry02/' 'cntry02/.DS_Store' 'cntry02/cntry02.avl' 'cntry02/cntry02.dbf'
Columns 5 through 7
'cntry02/cntry02.prj' 'cntry02/cntry02.sbn' 'cntry02/cntry02.sbx'
Columns 8 through 11
'cntry02/cntry02.shp' [1x23 char] 'cntry02/cntry02.shx' 'cntry02/cntry02.ssx'
Download the world_borders shapefiles. The world_borders shapefiles can be downloaded from:
http://mappinghacks.com/data/world_borders.zip
url = 'http://mappinghacks.com/data/world_borders.zip'; try files = unzip(url) catch error('Unable to download %s.', url); end
files =
'world_borders.dbf' 'world_borders.shp' 'world_borders.shx'
The following examples illustrate making maps from political boundaries defined by data in cntry98.shp.
Draw a map of the World using the countries defined by cntry98.shp.
figure worlddatamap cntry98.shp world title({'World Countries','Using cntry98.shp'})
Draw a map of the countries of Europe as a line.
figure worlddatamap('cntry98.shp','europe','line') title({'European Countries','Using cntry98.shp'})
Draw a map of North and South America as a patch and include the major rivers and lakes.
figure worlddatamap('cntry98.shp',{'North America','South America'}, ... 'patch',{'lakes', 'rivers'}); title({'North and South America','Using cntry98.shp'})
These examples use data from cntry02.shp.
Draw a map of the World using the countries defined by cntry02.shp.
figure worlddatamap cntry02.shp world title({'World Countries','Using cntry02.shp'})
Draw a map of the countries of Europe.
figure worlddatamap('cntry02.shp','europe','patch') title({'European Countries','Using cntry02.shp'})
Draw a map of the cities of South America.
figure worlddatamap('cntry02.shp','south america','cities') title({'Cities of South America','Using cntry02.shp'})
Draw a map of Africa with all cities, rivers, and lakes. Set the width of the Nile river to be larger than the others.
figure worlddatamap('cntry02.shp','africa',{'rivers','lakes'}) nile = handlem('Nile'); whiteNile = handlem('Nile - White Nile'); set([nile; whiteNile],'linewidth',4.0) title({'African Countries with Major Rivers and Lakes','Using cntry02.shp'})
Draw a map of all landlocked countries. List the attribute names of the file cntry02.shp. Notice that the dataset contains a LANDLOCKED field name. Use the LANDLOCKED attribute (Y) to display only the countries that are landlocked.
info = shapeinfo('cntry02.shp') {info.Attributes.Name} figure worlddatamap({'cntry02.shp','landlocked'},'yonly','patch') load coast plotm(lat,long) title({'Landlocked Countries','Using cntry02.shp'})
info =
Filename: [3x52 char]
ShapeType: 'Polygon'
BoundingBox: [2x2 double]
NumFeatures: 251
Attributes: [14x1 struct]
ans =
Columns 1 through 5
'FIPS_CNTRY' 'GMI_CNTRY' 'ISO_2DIGIT' 'ISO_3DIGIT' 'CNTRY_NAME'
Columns 6 through 11
'LONG_NAME' 'SOVEREIGN' 'POP_CNTRY' 'CURR_TYPE' 'CURR_CODE' 'LANDLOCKED'
Columns 12 through 14
'SQKM' 'SQMI' 'COLOR_MAP'
The following example illustrates making a map from political boundaries defined by data in world_borders.shp
Draw a map of Australia including major cities.
figure worlddatamap world_borders.shp Australia patch cities title({'Australia','Using world\_borders.shp'})
The following examples illustrate making maps of the World with data from the Global Self-consistant Hierarchical High-resolution Shorelines dataset.
Download GSHHS low resolution data.
url = 'http://www.ngdc.noaa.gov/mgg/shorelines/data/gshhs/gshhs_l.b.gz';
files = gunzip(url)
files =
'gshhs_l.b'
Create a shapefile of the GSHHS data.
Since the GSHHS data does not have a default REGIONATTRIBUTE field name of 'name', 'cntry_name', 'country_name', or 'long_name', add a Name field name with 'World' as the attribute value.
s = gshhs('gshhs_l.b'); [s.Name] = deal('World'); shapewrite(s, 'gshhs_l');
Create a map of global shorelines using GSHHS low resolution data.
figure worlddatamap gshhs_l.shp world title({'Global Shorelines','Using gshhs\_l.shp'})
Create a map of the Pacific region's shorelines.
figure worlddatamap gshhs_l.shp pacific title({'Pacific Region Shorelines','Using gshhs\_l.shp'})
