作者paranoia5201 (社会黑暗、行路难)
看板Statistics
标题[程式] SAS画地图,但结果不够完美
时间Fri Jan 12 18:02:23 2018
[软体程式类别]:
SAS
[程式问题]:
地图绘制
[软体熟悉度]:
熟悉
[问题叙述]:
我正在用SAS练习画一个疾病发生率的县市地图,地图是跑出来了,但还是有几个问题:
1. 有些县市的边界不见了(後来有跑出来,但是不晓得是修正了哪部分)。
2. 有些县市在地图上呈现空白,可是它们的资料各有内容,以下附上相关的log讯息:
3. 原本县市边界会消失的县市、与呈现空白的县市,我发现是同一群县市。
4. 希望能更改legand的名称,现在是ADJ_R、即是变项名称。
5. 希望颜色是改为渐层ex浅绿→深绿,下列的自动分五层、会是五种不同的颜色。
感谢~
[程式范例]:
proc mapimport
datafile="C:\Users\COUNTY_MOI_1060525.shp"
out=map.townmap;
run;
data adjr;
input countyname $ ADJ_R @@;
datalines;
彰化县 305.55
嘉义县 306.91
基隆市 294.2
新竹县 253.79
花莲县 290.73
宜兰县 292.74
苗栗县 278.47
南投县 271.14
澎湖县 216.92
屏东县 302.71
新北市 306.03
台北市 285.11
台东县 293.06
桃园市 288.52
云林县 300.12
台中市 320.16
台南市 323.06
高雄市 327.39
新竹市 294.1
嘉义市 335.72
金门县 140.13
连江县 132.19
; run;
/* Map01: legand自动分五个等级scale,部分县市无边界 & 皆空白*/
proc gmap data=adjr map=map.townmap;
id countyname;
choro ADJ_R / statistic=sum levels=5;
format ADJ_R degfmt.;
run;
log内容:
WARNING: ID variable countyname is of length 8 in response data set and
length 12 in map data set. This might result in ID values not matching or
multiple regions to becombined.
NOTE: 前景颜色 WHITE 与背景相同。您的图形可能有部分看不见。
/* Map02: legand各县市一个scale,颜色皆不同,仍是部分县市无边界*/
proc gmap data=adjr map=map.townmap;
id countyname;
choro ADJ_R / discrete coutline = black;
format ADJ_R degfmt.;
run;
log内容:
WARNING: ID variable countyname is of length 8 in response data set and
length 12 in map data set. This might result in ID values not matching or
multiple regions to becombined.
NOTE: 前景颜色 WHITE 与背景相同。您的图形可能有部分看不见。
/* Map03: legand各县市一个scale,颜色皆不同,各县市皆有边界,仍有空白县市*/
TITLE;FOOTNOTE;
TITLE1 "2015县市别发生率";
FOOTNOTE"资料来源:中心";
PROC GMAP GOUT=MAPCHART DATA=adjr MAP=map.townmap ALL ;
ID COUNTYNAME;
CHORO ADJ_R /
discrete
WOUTLINE=1
levels=5
statistic = sum
LEGEND= legend1
ANNOTATE=adjr;
RUN;
QUIT; TITLE;FOOTNOTE;GOPTIONS RESET=LEGEND;
log内容:
WARNING: ID variable countyname is of length 8 in response data set and
length 12 in map data set. This might result in ID values not matching or
multiple regions to becombined.
WARNING: Legend statement 1 could not be found. The default legend
description will be used.
NOTE: 在 ANNOTATE= 资料集 WORK.ADJR 中侦测到错误。
NOTE: 观测值 1 中有问题 - 已要求标签函数,但资料集中没有文字变数
NOTE: 已达到注解程序中的错误上限。处理已终止。
NOTE: 处理已依据个别错误计数而终止。
NOTE: 1 错误总计。
NOTE: 前景颜色 WHITE 与背景相同。您的图形可能有部分看不见。
--
※ 发信站: 批踢踢实业坊(ptt.cc), 来自: 210.241.78.121
※ 文章网址: https://webptt.com/cn.aspx?n=bbs/Statistics/M.1515751345.A.58A.html
1F:→ jasonfun44: 如果你部介意用gis用qgis可能会符合你的需求 02/02 19:17