请选择 进入手机版 | 继续访问电脑版
楼主: 匿名
58498 69

[作图问题求助] [Stata]地图绘制(一):中国地图与世界地图   [推广有奖]

匿名网友
匿名网友  发表于 2015-7-4 19:53:14
相似文件 换一批

+2 论坛币
k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群
赵安豆老师微信:zhaoandou666

经管之家联合CDA

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

感谢您参与论坛问题回答

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

+2 论坛币

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


世界地图

注:(该部分主要参考中南财经政法大学Stata协会讲座,在此向中南财经政法大学Stata协会致敬感谢)
  1. ssc install wbopendta //世界银行数据
  2. ssc install spmap   // 制图命令
复制代码
然后输入命令db wbopendata打开下载数据的向导

1.JPG


将数据画在地图上,首先,将研究数据和地图数据绑在一起
*准备地图数据
  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
复制代码

2.JPG


*拼接数据,绘制地图
  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)
复制代码

4.JPG




5.png


*最为简洁的命令

  1. worldstat world,stat(GDP) year(2012) cname
复制代码

6.png




*中国地图的绘制(该部分参考连玉君老师讲义,感谢连老师)

  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 ..)
复制代码

7.JPG


10.png





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

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

Graph.png


  1. //命令安装
  2. ssc install spgrid
  3. ssc install spkde
  4. ssc install spmap

  5. //数据集
  6. net get spkde, from(http://fmwww.bc.edu/RePEc/bocode/s)
  7. net get spmap, from(http://fmwww.bc.edu/RePEc/bocode/s)
  8. net get spgrid, from(http://fmwww.bc.edu/RePEc/bocode/s)

  9. //修改变量值单位
  10. use "italy-outlinecoordinates",clear
  11. replace _X=_X/1000
  12. replace _Y=_Y/1000
  13. save,replace

  14. //略去(kilometers)
  15. spgrid using "italy-outlinecoordinates",   ///
  16. resolution(w10)  ///
  17. cells("GridCells.dta")      ///
  18. points("GridPoints.dta")   ///
  19. replace compress dots

  20. use "Italy-DataPoints.dta", clear
  21. spkde dcvd95 pop95 using "GridPoints.dta",   ///
  22. xcoord(xcoord) ycoord(ycoord)              ///
  23. bandwidth(fbw) fbw(100) dots               ///
  24. saving("Kde.dta", replace)
复制代码

Graph.png


二维码

扫码加我 拉你入群

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

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

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

已有 2 人评分经验 论坛币 学术水平 热心指数 信用等级 收起 理由
うめし + 5 精彩帖子
xddlovejiao1314 + 100 + 100 + 5 + 5 + 5 精彩帖子

总评分: 经验 + 100  论坛币 + 105  学术水平 + 5  热心指数 + 5  信用等级 + 5   查看全部评分

本帖被以下文库推荐

刺客王朝 学生认证  发表于 2015-7-4 19:58:01 |显示全部楼层
已有 1 人评分论坛币 收起 理由
xddlovejiao1314 + 5 鼓励积极发帖讨论

总评分: 论坛币 + 5   查看全部评分

使用道具

gaojianwqjk 发表于 2015-7-4 22:13:21 来自手机 |显示全部楼层
niuniuyiwan 发表于 2015-7-4 19:53
*世界地图(该部分主要参考中南财经政法大学Stata协会讲座)            

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

使用道具

xddlovejiao1314 学生认证  发表于 2015-7-5 10:44:58 |显示全部楼层
谢谢分享。
已有 1 人评分经验 论坛币 收起 理由
niuniuyiwan + 30 + 3 感谢对楼主及本版的支持

总评分: 经验 + 30  论坛币 + 3   查看全部评分

使用道具

qitelata 发表于 2015-7-11 21:21:42 |显示全部楼层
这些方法的用途会越来越广
已有 2 人评分经验 论坛币 热心指数 收起 理由
niuniuyiwan + 30 + 3 + 1 感谢对楼主及本版的支持
xddlovejiao1314 + 5 鼓励积极发帖讨论

总评分: 经验 + 30  论坛币 + 8  热心指数 + 1   查看全部评分

使用道具

byjr211 发表于 2015-7-28 11:01:25 |显示全部楼层
赞!我要试一下 画地图 很有实用性的命令
已有 1 人评分经验 论坛币 热心指数 收起 理由
niuniuyiwan + 60 + 3 + 1 感谢对楼主及本版的支持

总评分: 经验 + 60  论坛币 + 3  热心指数 + 1   查看全部评分

使用道具

soulin27 发表于 2015-8-13 15:23:10 |显示全部楼层
非常详细的操作过程,楼主太棒了

使用道具

nc3722 发表于 2015-8-15 00:12:23 |显示全部楼层
可是如何用自己的数据画出中国地图呢,如何更改示例数据为自己的面板数据呢。

使用道具

sunnyyyyy123 发表于 2015-8-15 00:14:48 |显示全部楼层
谢谢楼主u分享 好开心

使用道具

niuniuyiwan 在职认证  发表于 2015-9-8 07:55:14 |显示全部楼层

graph.gif



使用道具

您需要登录后才可以回帖 登录 | 我要注册

京ICP备16021002-2号 京B2-20170662号 京公网安备 11010802022788号 论坛法律顾问:王进律师 知识产权保护声明   免责及隐私声明

GMT+8, 2023-1-31 22:41