用sas的gmap也可以作出这种图,可以用颜色代表区域,用柱子代表城市。sas里有现成的例子,代码如下:
goptions reset=all border;
proc format;
value popfmt low-1000000="0-1"
1000001-10000000="1-10"
10000001-100000000="10-100"
100000001-500000000="100-500"
500000001-high="over 500";
run;
proc format;
value $ regionfmt "SEAR" = "South-East Asia"
"EUR" = "Europe"
"EMR" = "Eastern Mediterranean"
"WPR" = "Western Pacific";
run;
title1 "Population Data for Asia (2005)";
pattern1 color=gold value=msolid;
pattern2 color=yellow value=msolid;
pattern3 color=cyan value=msolid;
pattern4 color=light_blue value=msolid;
legend1 label=("Population (Millions)");
proc gmap data=sashelp.demographics(where=(cont=95))
map=maps.asia all;
format pop popfmt.;
format region $regionfmt.;
id id;
area region / midpoints="SEAR" "EUR" "EMR" "WPR";
block pop / legend=legend1
relzero
levels=all;
run;
quit;