全部版块 我的主页
论坛 计量经济学与统计论坛 五区 计量经济学与统计软件 Stata专版
9333 5
2016-01-01
    大家好,我现在需要绘制中国地级行政区的等值图(choropleth maps),网上有很多教程教授如何在stata上完成这一任务,大体不出“先用--shp2dta-- 将shp 和 dbf文件转成stata能够读取的.dta文件,在和需要分析的数据merge起来”这一思路,我用http://www.stata.com/support/faqs/graphics/spmap-and-maps/上的数据跑了一遍,大体也没有问题,但是现在问题出在中国找不到合适的,地级市级别的shp文件


  美国shp转成.dta文件后的内容: (下载地址: http://www.nws.noaa.gov/geodata/catalog/national/html/us_state.htm )              
---------------------------------------------------------------------------------------------------------
              storage   display    value
variable name   type    format     label      variable label
---------------------------------------------------------------------------------------------------------
OBJECTID        int     %10.0g                OBJECTID
STATE            str2    %9s                   STATE
NAME             str20   %20s                  NAME
FIPS               str2    %9s                   FIPS
LON               float   %9.0g                 LON
LAT               float   %9.0g                 LAT
Shape_Leng    float   %9.0g                 Shape_Leng
Shape_Area     float   %9.0g                 Shape_Area
id                  byte    %12.0g               
---------------------------------------------------------------------------------------------------------
可以看到包括了州名,经纬度,还有生成的id

我在论坛上找到了标注 最新地级市shp文件的帖子(https://bbs.pinggu.org/thread-3956363-1-1.html?mobile=no),发现里面的若干shp文件用--shp2dta--转换后和上图的变量很不一样,比较类似一个是


---------------------------------------------------------------------------------------------------------
              storage   display    value
variable name   type    format     label      variable label
---------------------------------------------------------------------------------------------------------
AREA            double  %10.0g                AREA
PERIMETER       double  %10.0g                PERIMETER
BOU2_4M_        int     %10.0g                BOU2_4M_
BOU2_4M_ID      int     %10.0g                BOU2_4M_ID
ADCODE93        long    %10.0g                ADCODE93
ADCODE99        long    %10.0g                ADCODE99
NAME            str16   %16s                  NAME
id              int     %12.0g               
---------------------------------------------------------------------------------------------------------

   AREA   PERIME~R   BOU2_4M_   BOU2_4~D   ADCODE93   ADCODE99               NAME   id |
     |----------------------------------------------------------------------------------------|
  1. |  54.447     68.489          2         23     230000     230000           黑龙江省    1 |
  2. | 129.113    129.933          3         15     150000     150000       内蒙古自治区    2 |
  3. | 175.591     84.905          4         65     650000     650000   新疆维吾尔自治区    3 |
  4. |  21.315     41.186          5         22     220000     220000             吉林省    4 |
  5. |  15.603     38.379          6         21     210000     210000             辽宁省    5 |
     +----------------------------------------------------------------------------------------+


没有经纬度,没有label, 而且有的省份有好几个观察值,有的只有一个。这种情况怎么处理,有没有适合 --shp2dta-- 和--spmap-- 使用的 shp文件??















二维码

扫码加我 拉你入群

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

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

全部回复
2016-1-2 19:44:14
  有谁用过这个命令吗?
二维码

扫码加我 拉你入群

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

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

2017-3-12 21:31:38
同求,论坛上标价很高的帖子买了不少,没一个好用的
二维码

扫码加我 拉你入群

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

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

2019-5-5 14:47:05
https://gadm.org/download_country_v3.html
这里可以下载中国国家级、省级、地市级、县级四个层级的shp 数据
二维码

扫码加我 拉你入群

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

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

2019-7-15 23:56:00
kangshulong 发表于 2019-5-5 14:47
https://gadm.org/download_country_v3.html
这里可以下载中国国家级、省级、地市级、县级四个层级的shp 数 ...
注意,这个地图数据不能用,中印边界那里是错的
二维码

扫码加我 拉你入群

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

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

2021-2-6 00:59:18
yishuchanglq 发表于 2019-7-15 23:56
注意,这个地图数据不能用,中印边界那里是错的
哇,你这个可以的
二维码

扫码加我 拉你入群

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

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

相关推荐
栏目导航
热门文章
推荐文章

说点什么

分享

扫码加好友,拉您进群
各岗位、行业、专业交流群