GMAP程序和中国地图文件（包括县）

相似文件 换一批

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

今天学习了下GMAP，代码见下，供新手学习！

1. proc mapimport datafile="E:\diji\Bount_poly.shp" out=chinamap;
   
2. run;
   
3. proc sort data=chinamap(where=(SH2="32")) out=response nodupkey;
   
4. by adcode99;
   
5. run;
   
6. data response;
   
7. set response;
   
8. if di2=01 then city="南京";if di2=01 then gdp=9720;
   
9. if di2=02 then city="无锡";if di2=02 then gdp=8518;
   
10. if di2=03 then city="徐州";if di2=03 then gdp=5319;
    
11. if di2=04 then city="常州";if di2=04 then gdp=5273;
    
12. if di2=05 then city="苏州";if di2=05 then gdp=14504;
    
13. if di2=06 then city="南通";if di2=06 then gdp=6148;
    
14. if di2=07 then city="连云港";if di2=07 then gdp=2160;
    
15. if di2=08 then city="淮安";if di2=08 then gdp=2745;
    
16. if di2=09 then city="盐城";if di2=09 then gdp=4213;
    
17. if di2=10 then city="扬州";if di2=10 then gdp=4016;
    
18. if di2=11 then city="镇江";if di2=11 then gdp=3502;
    
19. if di2=12 then city="泰州";if di2=12 then gdp=3655;
    
20. if di2=13 then city="宿迁";if di2=13 then gdp=2110;
    
21. run;
    
22. goptions reset=all;
    
23. legend1 label=("江苏省地级市");
    
24. legend2 down=10 label=("江苏省地级市GDP（亿元）") position=(middle right outside) across=1;
    
25. proc gmap map=chinamap(where=(sh2="32")) data=response all;
    
26. id name99;
    
27. choro name99;
    
28. run;
    
29. proc gmap map=chinamap(where=(sh2="32")) data=response all;
    
30. id name99;
    
31. choro city/legend=legend1;
    
32. run;
    
33. proc gmap map=chinamap(where=(sh2="32")) data=response all;
    
34. id di2;
    
35. choro gdp/legend=legend2 levels=13;
    
36. run;
    
37. proc gmap map=chinamap(where=(sh2="32")) data=response all;
    
38. id di2;
    
39. choro gdp/legend=legend2 discrete levels=13;
    
40. run;
    
41. proc gmap map=chinamap(where=(sh2="32")) data=response all;
    
42. id di2;
    
43. prism gdp/legend=legend2 discrete levels=13;
    
44. run;
    
45. proc gmap map=chinamap(where=(sh2="32")) data=response all;
    
46. id di2;
    
47. area di2/legend=legend2;
    
48. block gdp/legend=legend2;
    
49. run;
    
50. proc print data=chinamap(obs=5);
    
51. run;
    
52. proc print data=response;
    
53. run;
    
54. quit;

复制代码

diji.zip (10.33 MB, 需要: 2 个论坛币) 本附件包括：

• BOUNT_poly.shp
• BOUNT_poly.shx
• BOUNT_poly.dbf

2016-11-21 02:09:52 上传

包括县的中国地图
需要: 2 个论坛币  [购买]

中华人民共和国行政区划代码标准.doc (3.97 MB, 需要: 2 个论坛币)

2016-11-21 02:04:11 上传

中国行政区代码
需要: 2 个论坛币  [购买]

扫码加我 拉你入群

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

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

分享0 收藏2 回帖

关键词：中国地图 gmap 中国地 Map GMA 中国地图 程序

不知道哪位大大知道如何将按县画地图改为按地级市画图？，将地级市内的县的分界去掉？
或者如何利用县的地图画地级市的地图？

mark mark