{"id":6125,"date":"2015-08-14T09:00:58","date_gmt":"2015-08-14T13:00:58","guid":{"rendered":"https:\/\/blogs.mathworks.com\/pick\/?p=6125"},"modified":"2015-08-12T19:49:38","modified_gmt":"2015-08-12T23:49:38","slug":"geom-2d","status":"publish","type":"post","link":"https:\/\/blogs.mathworks.com\/pick\/2015\/08\/14\/geom-2d\/","title":{"rendered":"Geom 2d"},"content":{"rendered":"
\r\n \r\n

Sean<\/a>'s pick this week is geom2d<\/a> by David Legland<\/a>.\r\n <\/p>\r\n <\/introduction>\r\n\r\n

Working at MathWorks, I get to use MATLAB every day, learn a lot about it, and get to try everything. But one of the things\r\n I forget is the experience of someone who doesn't use it every day but occasionally needs it for a specific task. Thus, when\r\n I get the opportunity to solve a problem like this, it reminds me of my early days and reminds me what a privilege it is to\r\n have access to professional quality software.\r\n <\/p>\r\n

Yesterday, one of our program managers came to my office and wanted me to create a map of the United States for her with different\r\n regions color coded. Seems straight forward right? Here's the story of the steps I took to create the map for her.\r\n <\/p>\r\n

I knew off the top of my head, that usamap<\/tt><\/a> in the Mapping Toolbox<\/a> was the function that I needed to build the map. So I went to the doc page linked above by running:\r\n <\/p>

doc usamap<\/span><\/pre>
The third example at the bottom of the doc page was exactly what I was looking for.  I almost always start with examples because\r\n      it's easier to modify or delete than to add.  As an added bonus, it already even excluded Alaska and Hawaii which she didn't\r\n      want on her map!\r\n   <\/p>\r\n   
 <\/p>\r\n   
I copied this example in, and ran it.  I quickly realized that the line I really needed to modify was the following one; it\r\n      even had a nice comment noting the current randomness.\r\n   <\/p>
faceColors = makesymbolspec('Polygon',...\r\n   {'INDEX', [1 numel(states)], 'FaceColor', ...\r\n   polcmap(numel(states))}); %NOTE - colors are random<\/pre>
I figured polcmap<\/tt> probably worked like any of the other colormap functions and takes the number of unique colors to make.  My task involved\r\n      six regions and so I needed six colors.\r\n   <\/p>
figure\r\npolcmap6 = polcmap(6);\r\ncolormap(polcmap6)\r\ncolorbar<\/pre> 
Minor victory!<\/p>\r\n   
Now I needed to index into this colormap each state's region.  I asked for an Excel sheet containing the states and the region\r\n      indices.\r\n   <\/p>\r\n   
 <\/p>\r\n   
This can be easily read in with readtable<\/tt><\/a>.\r\n   <\/p>
Tstateregions = readtable('StateRegionsLookup.xlsx'<\/span>);<\/pre>
Warning: Variable names were modified to make\r\nthem valid MATLAB identifiers. \r\n<\/pre>
With the region indices, I should be able to plot the map color coded by indexing into the colormap:<\/p>
figure\r\nax = usamap('conus'<\/span>);\r\n\r\nstates = shaperead('usastatelo'<\/span>, 'UseGeoCoords'<\/span>, true,...<\/span>\r\n  'Selector'<\/span>,...<\/span>\r\n  {@(name) ~any(strcmp(name,{'Alaska'<\/span>,'Hawaii'<\/span>})), 'Name'<\/span>});\r\n\r\nfaceColors = makesymbolspec('Polygon'<\/span>,...<\/span>\r\n    {'INDEX'<\/span>, [1 numel(states)], 'FaceColor'<\/span>, ...<\/span>\r\n    polcmap6(Tstateregions.Region,:)});\r\n\r\ngeoshow(ax, states, 'DisplayType'<\/span>, 'polygon'<\/span>, ...<\/span>\r\n   'SymbolSpec'<\/span>, faceColors)\r\n\r\n\r\nframem off<\/span>\r\ngridm off<\/span>\r\nmlabel off<\/span>\r\nplabel off<\/span><\/pre> 
But that doesn't look right!  So what went wrong, I took a look at the states and compared them:<\/p>
table({states(:).Name}.', Tstateregions.State,'VariableNames'<\/span>,{'MATLAB'<\/span>,'Excel'<\/span>})<\/pre>
ans = \r\n            MATLAB                     Excel          \r\n    ______________________    ________________________\r\n    'Alabama'                 ''Alabama''             \r\n    'Arizona'                 ''Arizona''             \r\n    'Arkansas'                ''Arkansas''            \r\n    'California'              ''California''          \r\n    'Colorado'                ''Colorado''            \r\n    'Connecticut'             ''Connecticut''         \r\n    'Delaware'                ''Delaware''            \r\n    'Florida'                 ''District of Columbia''\r\n    'Georgia'                 ''Florida''             \r\n    'Idaho'                   ''Georgia''             \r\n    'Illinois'                ''Idaho''               \r\n    'Indiana'                 ''Illinois''            \r\n    'Iowa'                    ''Indiana''             \r\n    'Kansas'                  ''Iowa''                \r\n    'Kentucky'                ''Kansas''              \r\n    'Louisiana'               ''Kentucky''            \r\n    'Maine'                   ''Louisiana''           \r\n    'Maryland'                ''Maine''               \r\n    'Massachusetts'           ''Maryland''            \r\n    'Michigan'                ''Massachusetts''       \r\n    'Minnesota'               ''Michigan''            \r\n    'Mississippi'             ''Minnesota''           \r\n    'Missouri'                ''Mississippi''         \r\n    'Montana'                 ''Missouri''            \r\n    'Nebraska'                ''Montana''             \r\n    'Nevada'                  ''Nebraska''            \r\n    'New Hampshire'           ''Nevada''              \r\n    'New Jersey'              ''New Hampshire''       \r\n    'New Mexico'              ''New Jersey''          \r\n    'New York'                ''New Mexico''          \r\n    'North Carolina'          ''New York''            \r\n    'North Dakota'            ''North Carolina''      \r\n    'Ohio'                    ''North Dakota''        \r\n    'Oklahoma'                ''Ohio''                \r\n    'Oregon'                  ''Oklahoma''            \r\n    'Pennsylvania'            ''Oregon''              \r\n    'Rhode Island'            ''Pennsylvania''        \r\n    'South Carolina'          ''Rhode Island''        \r\n    'South Dakota'            ''South Carolina''      \r\n    'Tennessee'               ''South Dakota''        \r\n    'Texas'                   ''Tennessee''           \r\n    'Utah'                    ''Texas''               \r\n    'Vermont'                 ''Utah''                \r\n    'Virginia'                ''Vermont''             \r\n    'Washington'              ''Virginia''            \r\n    'West Virginia'           ''Washington''          \r\n    'Wisconsin'               ''West Virginia''       \r\n    'Wyoming'                 ''Wisconsin''           \r\n    'District of Columbia'    ''Wyoming''             \r\n<\/pre>
Ahh ha!  Washington DC is appended to the end of MATLAB's shipping state border file where as it's alphabetical in the Excel\r\n      sheet causing an offset.  Simple enough manual fix:\r\n   <\/p>
Tstateregions = [Tstateregions([1:7 9:end],:); Tstateregions(8,:)];<\/pre>
figure\r\nax = usamap('conus'<\/span>);\r\n\r\nfaceColors = makesymbolspec('Polygon'<\/span>,...<\/span>\r\n    {'INDEX'<\/span>, [1 numel(states)], 'FaceColor'<\/span>, ...<\/span>\r\n    polcmap6(Tstateregions.Region,:)});\r\n\r\ngeoshow(ax, states, 'DisplayType'<\/span>, 'polygon'<\/span>, ...<\/span>\r\n   'SymbolSpec'<\/span>, faceColors)\r\n\r\n\r\nframem off<\/span>\r\ngridm off<\/span>\r\nmlabel off<\/span>\r\nplabel off<\/span><\/pre> 
By now you're probably wondering what this has to do with geom2d<\/tt>.  Well I still had to put the two letter code on each state.  To do this I needed a way to identify the centroid of the polygon.\r\n       I went to the File Exchange and searched for \"polygon centroid\", and sure enough geom2d<\/tt> has a function named exactly that!\r\n   <\/p>
for<\/span> ii = 1:numel(states)\r\n    latii = states(ii).Lat; % Remove nans indicating end of border<\/span>\r\n    lonii = states(ii).Lon;\r\n    cent = polygonCentroid(latii(~isnan(latii)),lonii(~isnan(lonii)));\r\n    textm(cent(1),cent(2),Tstateregions.x2letter(ii))\r\nend<\/span><\/pre> 
And there you have it.<\/p>\r\n   
Comments<\/a><\/h3>\r\n   
Do you ever get a random one-off task or activity that is quickly solvable in MATLAB and that makes you feel good?<\/p>\r\n   
Give it a try and let us know what you think here<\/a> or leave a comment<\/a> for David.\r\n   <\/p>
The polygon centroid function is just one of dozens in this well documented and organized toolbox.  There are tools for all\r\n      sorts of math, aggregation and plotting of geometric primitives, e.g. polygons, lines, rays, and polynomial curves.  Each\r\n      group of functions has its own Contents.m<\/i> file so one can quickly see all of the functions and a quick description of what they do.  The functions have a full help\r\n      description, many with an example.  There are also a bunch of example files that have been published so you can see the example\r\n      without having to run it. Overall this is a great package to have around and now it can join its older brother, geom3d<\/tt><\/a> with a Pick of the Week banner!\r\n   <\/p>\r\n