proc gmap map=reflib.can2(where=(density<3)) data=reflib.can2 all; id province; choro province / nolegend; run; quit;