Here is a whole example so you can see MAPSGFK and GPROJECT: data schools (label='Wake County, NC public schools'); infile datalines dlm=','; length school $64 address $32 city $24 state $2 zip 5 type $12 color $10; input school address city state zip type; type = upcase(type); /* Set point color based on school type. */ if type = 'HIGH' then color = 'CXff0000'; else if type = 'MIDDLE' then color = 'CX00ff00'; else if type = 'ELEMENTARY' then color = 'CX0000ff'; else delete; /* Drop multi-type schools */ datalines; Adams Elementary, 805 Cary Towne Blvd., Cary, NC, 27511, Elementary Apex High, 1501 Laura Duncan Road, Apex, NC, 27502, High Apex Middle, 400 E. Moore Street, Apex, NC, 27502, Middle Athens Drive High, 1420 Athens Drive, Raleigh, NC, 27606, High ; run; proc geocode method=street data=work.schools out=work.geocoded lookupstreet=sashelp.geoexm; run; %make_dots(anno,work.geocoded,color,1); data anno; set anno; flag=1; run; data us_states(drop=state); set mapsgfk.us_states; run; data us_states; set us_states; state=statecode; x=long; y=lat; if state="AK" or state="HI" then; else output; run; data all; set us_states anno; run; proc sort data=all out=allsorted; by state; run; proc gproject data=allsorted out=annop project=albers degrees eastlong; id state; run; data annodata usmap; set annop; if flag=1 then output annodata; else output usmap; run; proc gmap data=usmap map=usmap anno=annodata; id state; choro state/ nolegend; run; quit;
... View more