用--shp2dta-- 绘制中国地级市行政区图

zlqs1985

9333

收藏 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文件？？