The GREDUCE procedure creates a new map data set, 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=can2;
   id province;
run;