[Stata]地图绘制（一）：中国地图与世界地图

相似文件 换一批

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

中国地图与世界地图的绘制

世界地图

注：（该部分主要参考中南财经政法大学Stata协会讲座，在此向中南财经政法大学Stata协会致敬感谢）

1. ssc install wbopendta //世界银行数据
   
2. ssc install spmap   // 制图命令

复制代码

然后输入命令db wbopendata打开下载数据的向导

将数据画在地图上，首先，将研究数据和地图数据绑在一起
*准备地图数据

1. use  "http://fmwww.bc.edu/repec/bocode/w/world-c.dta",clear        
   
2. save world-c
   
3. use  "http://fmwww.bc.edu/repec/bocode/w/world-d.dta",clear        
   
4. save world-d

复制代码

*拼接数据，绘制地图

1. merge 1:1 countrycode using world-d.dta,keepusing(_ID)
   
2. rename _ID id
   
3. spmap ny_gdp_pcap_kd using world-c,id(id) fcolor(Heat)

复制代码

*最为简洁的命令

1. worldstat world,stat(GDP) year(2012) cname

复制代码

*中国地图的绘制（该部分参考连玉君老师讲义，感谢连老师）

1. findit china map
   
2.   
   
3.     use china_label,clear
   
4.     tab name
   
5.     replace name = subinstr(name, "省", "", .)
   
6.     replace name = subinstr(name, "市", "", .)
   
7.     replace name = subinstr(name, "回族自治区", "", .)
   
8.     replace name = subinstr(name, "壮族自治区", "", .)
   
9.     replace name = subinstr(name, "特别行政区", "", .)
   
10.     replace name = subinstr(name, "自治区", "", .)
    
11.     replace name = subinstr(name, "维吾尔", "", .)
    
12.     tab name
    
13.     gen x = uniform()
    
14.     format x %9.3g
    
15.   
    
16.     spmap x using "china_map.dta", id(id)                  ///
    
17.       label(label(name)                                    ///
    
18.         xcoord(x_coord) ycoord(y_coord) size(*.9))         ///
    
19.           plotregion(icolor(stone)) graphregion(icolor(stone)) ///
    
20.       clnumber(8) fc(Greens2) oc(white ..) osize(medthin ..)

复制代码

---

“一个酒鬼眼中的中国地图”* Dr. Chuntao Li（中南财经政法大学）绘图讲义

use china_label ,clear replace name="北大荒" if name=="黑龙江省" replace name="马奶酒" if name=="内蒙古自治区" replace name="葡萄美酒" if name=="新疆维吾尔自治区" replace name="宝丰" if name=="河南省" replace name="大泉源酒" if name=="吉林省" replace name="道光廿五酒" if name=="辽宁省" replace name="古河州" if name=="甘肃省" replace name="老白干" if name=="河北省" replace name="二锅头" if name=="北京市" replace name="汾酒" if name=="山西省" replace name="青稞酒" if name=="青海省" replace name="西凤酒" if name=="陕西省" replace name="老银川" if name=="宁夏回族自治区" replace name="青稞酒" if name=="西藏自治区" replace name="青岛啤酒" if name=="山东省" replace name="洋河" if name=="江苏省" replace name="古井" if name=="安徽省" replace name="五粮液" if name=="四川省" replace name="劲酒" if name=="湖北省" replace name="" if name=="上海市" replace name="江津" if name=="重庆市" replace name="花雕" if name=="浙江省" replace name="四特" if name=="江西省" replace name="酒鬼" if name=="湖南省" replace name="云南白" if name=="云南省" replace name="茅台" if name=="贵州省" replace name="刘三姐" if name=="广西壮族自治区" replace name="珠江" if name=="广东省" replace name="椰岛" if name=="海南省" replace name="金门高粱" if name=="台湾省" replace name="津酒" if name=="天津市" replace name="干红" if name=="香港特别行政区" replace name="金门高粱" if name=="福建省" format y %9.3g spmap y using "china_map.dta", /// id(id)title("用Stata13.1绘制的酒鬼眼中的中国地图", size(*0.8)) label(label(name) ///         xcoord(x_coord) ycoord(y_coord) size(*.76))fcolor(Oranges) legenda(off) 复制代码

1. //命令安装
   
2. ssc install spgrid
   
3. ssc install spkde
   
4. ssc install spmap
   
5. 
   
6. //数据集
   
7. net get spkde, from(http://fmwww.bc.edu/RePEc/bocode/s)
   
8. net get spmap, from(http://fmwww.bc.edu/RePEc/bocode/s)
   
9. net get spgrid, from(http://fmwww.bc.edu/RePEc/bocode/s)
   
10. 
    
11. //修改变量值单位
    
12. use "italy-outlinecoordinates",clear
    
13. replace _X=_X/1000
    
14. replace _Y=_Y/1000
    
15. save,replace
    
16. 
    
17. //略去(kilometers)
    
18. spgrid using "italy-outlinecoordinates",   ///
    
19. resolution(w10)  ///
    
20. cells("GridCells.dta")      ///
    
21. points("GridPoints.dta")   ///
    
22. replace compress dots
    
23. 
    
24. use "Italy-DataPoints.dta", clear
    
25. spkde dcvd95 pop95 using "GridPoints.dta",   ///
    
26. xcoord(xcoord) ycoord(ycoord)              ///
    
27. bandwidth(fbw) fbw(100) dots               ///
    
28. saving("Kde.dta", replace)

复制代码

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

关键词：Stata 中国地图 世界地图 tata 中国地 中国地图 世界地图

niuniuyiwan 发表于 2015-7-4 19:53
*世界地图（该部分主要参考中南财经政法大学Stata协会讲座）            

在安装wbopendata后使用命令 ...

:)

谢谢分享。

这些方法的用途会越来越广

赞！我要试一下 画地图 很有实用性的命令

非常详细的操作过程，楼主太棒了

可是如何用自己的数据画出中国地图呢，如何更改示例数据为自己的面板数据呢。

谢谢楼主u分享 好开心