The GREDUCE procedure creates a new map data set, REFLIB.CAN2, containing a DENSITY variable. The ID statement specifies the variable in the map data set that defines unit areas.
proc greduce data=maps.canada2 out=reflib.can2;
   id province;
run;