我试图通过各县与之相关的价值来绘制美国大陆的县。我正在使用spmap
包。我的形状文件来自:https://www.census.gov/geo/maps-data/data/cbf/cbf_counties.html。我使用的是500k版本。
capture ssc install spmap
capture ssc install shp2dta
我解压缩文件。然后我按如下方式加载地图:
shp2dta using "cb_2015_us_county_500k.shp", database(county2015db) coordinates(county2015coord) replace
use county2015db, clear
destring STATEFP, replace
destring GEOID, replace
如果我不想遮蔽各个县,一切都很完美:
spmap using county2015coord if STATEFP<60 & STATEFP!=2 & STATEFP!=15, id(GEOID)
但是,如果我执行以下操作(仅使用_ID作为此玩具示例的值),则只会填充少数几个县(并且地图会小得多并且会显示阿拉斯加):
spmap _ID using county2015coord if STATEFP<60 & STATEFP!=2 & STATEFP!=15, id(GEOID) fcolor(Blues)
怎么办?
答案 0 :(得分:2)
我应该在_ID
字段中使用id( )
而不是GEOID
。这解决了这个问题。