天涯在何方 发表于 2010-8-31 12:04 谢谢老兄，可是这些文件都有，需要的东西是如何根据我的主帖的附件，做出一份中国地图的完整的do命令

天涯在何方 发表于 2010-8-31 20:46 不知道这个是否可以满足STATA作图的条件

sungmoo 发表于 2010-9-1 07:03

天涯在何方 发表于 2010-8-31 20:46 不知道这个是否可以满足STATA作图的条件

*简单地话，可以直接借用2楼链接中的数据uscoord.dta与usdb.dta。再结合你的数据，可操作（设所有文件都在默认目录中）：

use usdb,clear
g Market=.
replace Market=8.619 if name=="北京市"
replace Market=8.342 if name=="天津市"
replace Market=6.405 if name=="河北省"
replace Market=5.262 if name=="山西省"
replace Market=5.521 if name=="内蒙古自治区"
replace Market=7.835 if name=="辽宁省"
replace Market=5.890 if name=="吉林省"
replace Market=5.263 if name=="黑龙江省"
replace Market=10.407 if name=="上海市"
replace Market=9.065  if name=="江苏省"
replace Market=9.896  if name=="浙江省"
replace Market=6.556  if name=="安徽省"
replace Market=8.624 if name=="福建省"
replace Market=6.225 if name=="江西省"
replace Market=8.210 if name=="山东省"
replace Market=6.198 if name=="河南省"
replace Market=6.652 if name=="湖北省"
replace Market=6.546 if name=="湖南省"
replace Market=10.057 if name=="广东省"
replace Market=5.818 if name=="广西壮族自治区"
replace Market=5.542 if name=="海南省"
replace Market=7.234 if name=="重庆市"
replace Market=6.860 if name=="四川省"
replace Market=4.572 if name=="贵州省"
replace Market=5.153 if name=="云南省"
replace Market=2.500 if name=="西藏自治区"
replace Market=4.797 if name=="陕西省"
replace Market=4.445 if name=="甘肃省"
replace Market=3.838 if name=="青海省"
replace Market=4.850 if name=="宁夏回族自治区"
replace Market=5.024 if name== "新疆维吾尔自治区"
spmap Market using "uscoord.dta",id(id) fc(Greys2) legs(2) cln(4) ndf(none)

*另外，你给的数据中还缺少label数据库，该数据库表明在地图上哪个位置（坐标）显示什么内容。

天涯在何方 发表于 2010-9-1 11:01 label文件确实是我没有找到，这个文件也不知道何处寻找

sungmoo 发表于 2010-9-1 11:12

天涯在何方 发表于 2010-9-1 11:01 label文件确实是我没有找到，这个文件也不知道何处寻找

这个文件可以根据自己的需要自己设计。里面列明你想要显示的内容以及显示的位置。

你想在地图上显示什么内容（label）？

sungmoo 发表于 2010-9-1 16:54
*结合上面的数据粗略做（需要细化，但基本原理不变）：

use uscoord,clear
collapse _X _Y,by(_ID)
sort _ID
save label,replace

use usdb,clear
ren id _ID
sort _ID
mer _ID using label
drop _merge

g Market=.
replace Market=8.619 if name=="北京市"
……
replace Market=5.024 if name=="新疆维吾尔自治区"

g obs=.
replace obs=160 if name=="北京市"
……
replace obs=400  if name=="新疆维吾尔自治区"

g label=substr(name,1,4)+" "+string(obs)+" "+string(Market) ///
   if (_ID16&_ID21&_ID

macross509 发表于 2010-9-1 18:34 搞什么，这么复杂，用Arc GIS两下就完了，搞了这么多命令，晕

zyjbs 发表于 2011-3-16 20:36 tmap choropleth pop1990, id(id) map(us_coord) palette(Blues)

zyjbs 发表于 2011-3-17 16:33 您说的“这一步可以实现吗”我没看懂您的意思，是我的命令输入不对吗

zyjbs 发表于 2011-3-16 20:36 以下是我的程序：
sshp2dta using s_01au07, database(us_db) coordinates(us_coord) genid(id)
use us_db, clear
describe
list id NAME
use stats
merge scode using trans, sort unique
tabulate _merge
drop _merge
merge id using us_db, sort unique
tabulate _merge
drop _merge
tmap choropleth pop1990, id(id) map(us_coord) palette(Blues)
前面都行 但是到最后制图不行 总是出现这样的提示：

. tmap choropleth pop1990, id(id) map(us_db) palette(Blues)
file us_db specified in option map() does not exist
r(498);

linglinboy 发表于 2011-3-17 20:02
. tmap choropleth pop1990, id(id) map(us_coord.dta) palette(Blues)
invalid file specification
r(198);
. run "C:\Documents and Settings\hp\桌面\10090113404cbace5894657b4d\china\sungmoo.do"
unrecognized command:  spmap
r(199);

linglinboy 发表于 2011-3-18 21:35 下了这个帖子里的数据，别的就没了。http://www.pinggu.org/bbs/thread-337213-1-1.html

sungmoo 发表于 2011-3-19 07:02

linglinboy 发表于 2011-3-18 21:35 下了这个帖子里的数据，别的就没了。http://www.pinggu.org/bbs/thread-337213-1-1.html

http://www.pinggu.org/bbs/thread-331184-2-1.html中下到usdb与uscoord两个文件。

http://www.pinggu.org/bbs/thread-337213-2-1.html

做图前要准备三份文件：

（1）各地区底层（二维）坐标数据库：里面列明各地区对应的（二维）坐标序列——（uscoord）；

（2）各地区label数据库：里面列明各地区对应的label的内容及其显示的位置（二维坐标）；

（3）各地区属性数据库：里面列明各地区对应的属性（比如各地区的某一经济指标——usdb。

三份数据库的地区变量的各值的意义必须相同（变量名可以不完全相同）。

先在usdb文件生成属性变量aaa（它就是要画图所表现的指标），然后画关于aaa的地图：

spmap aaa using uscoord, id(id)