Mapping Puerto Rico at the county (Municipio) level using mapsgfk

[mitchell_keener]

I apologize if this is the wrong board, but since I am using PROC GMAP, I posted here. I am using SAS 9.4 TS1M5. 

I am new to mapping and I am trying to map Puerto Rico using the mapsgfk.us_all dataset supplied with SAS. I need the map to be at the county (municipio) level. When I create the map using this dataset, it maps it at the county subdivision level. My data is still mapping correctly (filling up each sub division to fill the entire county), but considering my data is at the county level I would rather not give the impression my data is more granular. Does anyone know of a way to remove the subdivision lines and only show the major county lines? The odd thing is that it seems Puerto Rico has more granular information than the rest of the US in this file. When using the ID field in mapsgfk.us_all for the id statement in PROC GMAP, it maps all of the US correctly using a code such as 'US-15001' and the ID2 field is empty, but for Puerto Rico data, the ID field looks more granular such as 'PR-7200100401' and the ID2 field looks more like the other identifiers 'PR-72001'. When I try to use ID2 to map the data, it looks very distored. I included code below which maps Puerto Rico at the subdivision level and I included images of what the non-distorted map looks like (using ID) and what the distorted map looks like (using ID2 which seems to be at the county level). I have tried using PROC GREDUCE, but it seems like this is not solving the problem but I wonder if it is needed here. If anyone knows how to only map the county lines and not the subdivision lines, I would greatly appreciate it. Thank you!

proc gmap map=mapsgfk.us_all data=mapsgfk.us_all ;
   id id;
   choro state / nolegend ;
   where statecode="PR";
run;
	quit;

[mitchell_keener · Accepted solution]

Well FWIW I believe I answered my own question. This requires PROC GREMOVE and not PROC GREDUCE. I will leave this up in case it happens to help someone down the road. 

Fairly simply done using this code. 

data PR;
	set mapsgfk.us_all;
	where statecode="PR";
	run;

proc gremove data=PR out=PR2;
	by id2;
	id id;
run;



proc gmap map=PR2 data=PR2 all ;
   id id2;
   choro state / nolegend ;
   where statecode="PR";
run;
	quit;

[mkeintz]

@mitchell_keener 

Glad you found the answer to your question.    Go ahead and mark your own answer as the solution.

[GraphGuy]

Nice job of self-help, Mitchell!

Here are a couple of suggestions you might want to add...

When you subset your map to PR, I recommend dropping the resolution variable. The mapping procs (such as gmap) automatically apply reduce the number of points used in drawing the map, based on the size of the output ... but when you subset just a piece of a larger map, this throws off the calculations, and the results might or might not look good. If you delete the resolution variable, then the automated algorithms will not be applied.

After you get rid of the resolution variable and disable the automated system, then I recommend using the density values to just use the number/level of pixels detail as looks good to you. I tried several different settings, and I think density<=3 looks pretty good (not too sparse, and not too crowded borders).

And when I create a new map, I always like to use some test data with the map area names in it, and use that as mouse-over text (with gmap's html= option). You can see a sample mouse-over in the screen-capture below:

data PR; set mapsgfk.us_all
(where=(statecode="PR" and density<=3) drop=resolution);
run;

proc gremove data=PR out=PR2;
by id2;
id id;
run;

proc sql noprint;
create table pr2_data as
select unique id2, id2name, 'title='||quote(trim(left(id2name))) as my_html
from mapsgfk.us_all_attr;
quit; run;

title1 ls=1.5 "Puerto Rico Municipios";

goptions xpixels=800 ypixels=350;

proc gmap map=PR2 data=pr2_data all ;
id id2;
choro id2name / nolegend html=my_html;
run;