中国地图与世界地图的绘制 |
世界地图 |
- ssc install wbopendta //世界银行数据
- ssc install spmap // 制图命令
将数据画在地图上,首先,将研究数据和地图数据绑在一起
*准备地图数据
- use "http://fmwww.bc.edu/repec/bocode/w/world-c.dta",clear
- save world-c
- use "http://fmwww.bc.edu/repec/bocode/w/world-d.dta",clear
- save world-d
*拼接数据,绘制地图
- merge 1:1 countrycode using world-d.dta,keepusing(_ID)
- rename _ID id
- spmap ny_gdp_pcap_kd using world-c,id(id) fcolor(Heat)
*最为简洁的命令
- worldstat world,stat(GDP) year(2012) cname
*中国地图的绘制(该部分参考连玉君老师讲义,感谢连老师) |
- findit china map
-
- use china_label,clear
- tab name
- replace name = subinstr(name, "省", "", .)
- replace name = subinstr(name, "市", "", .)
- replace name = subinstr(name, "回族自治区", "", .)
- replace name = subinstr(name, "壮族自治区", "", .)
- replace name = subinstr(name, "特别行政区", "", .)
- replace name = subinstr(name, "自治区", "", .)
- replace name = subinstr(name, "维吾尔", "", .)
- tab name
- gen x = uniform()
- format x %9.3g
-
- spmap x using "china_map.dta", id(id) ///
- label(label(name) ///
- xcoord(x_coord) ycoord(y_coord) size(*.9)) ///
- plotregion(icolor(stone)) graphregion(icolor(stone)) ///
- clnumber(8) fc(Greens2) oc(white ..) osize(medthin ..)
“一个酒鬼眼中的中国地图”* Dr. Chuntao Li(中南财经政法大学)绘图讲义 |
|
- //命令安装
- ssc install spgrid
- ssc install spkde
- ssc install spmap
- //数据集
- net get spkde, from(http://fmwww.bc.edu/RePEc/bocode/s)
- net get spmap, from(http://fmwww.bc.edu/RePEc/bocode/s)
- net get spgrid, from(http://fmwww.bc.edu/RePEc/bocode/s)
- //修改变量值单位
- use "italy-outlinecoordinates",clear
- replace _X=_X/1000
- replace _Y=_Y/1000
- save,replace
- //略去(kilometers)
- spgrid using "italy-outlinecoordinates", ///
- resolution(w10) ///
- cells("GridCells.dta") ///
- points("GridPoints.dta") ///
- replace compress dots
- use "Italy-DataPoints.dta", clear
- spkde dcvd95 pop95 using "GridPoints.dta", ///
- xcoord(xcoord) ycoord(ycoord) ///
- bandwidth(fbw) fbw(100) dots ///
- saving("Kde.dta", replace)