Ahh - I forgot one thing in the gproject. I had specified the parmin= dataset, but I forgot to tell gproject which entry in that dataset to use.  Here's the new/fixed code:   data coordinates; input lat long; datalines; 28.9787 -95.1672 29.1468 -94.8645 29.2313 -89.8375 29.4695 -94.0852 29.5513 -92.899 29.5282 -93.9293 29.016 -95.098 29.8865 -88.1793 29.5088 -93.9617 29.4177 -92.587 29.1007 -94.02 29.5347 -93.5135 25.2593 -81.9585 28.9765 -91.7018 28.1155 -93.6067 28.3387 -94.4748 28.9765 -91.7018 ; run; data us48; set mapsgfk.us_states; if state not in (2 15 72); run; proc gproject data=us48 out=gulf latlon degrees eastlong longmin=-98 longmax=-68 latmin=25 latmax=36.5 parmout=projparm; where density<5; id state; run; /* project the coordinates using the same projection parameters as were used to project the map */ proc gproject data=coordinates out=coordinates latlon degrees eastlong longmin=-98 longmax=-68 latmin=25 latmax=36.5 parmin=projparm parmentry=gulf; id; run; /* convert the projected coordinates into annotate dataset */ data coordinates; set coordinates; xsys='2'; ysys='2'; hsys='3'; when='a'; function='pie'; rotate=360; size=1.0; color='red'; run; title "US Southern Atlantic and Gulf of Mexico"; pattern value=mempty color=blue; proc gmap map=gulf data=gulf all anno=coordinates; id state; choro state / statistic=first nolegend levels=1; ... View more