作者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/m.aspx?n=bbs/Statistics/M.1515751345.A.58A.html
1F:→ jasonfun44: 如果你部介意用gis用qgis可能會符合你的需求 02/02 19:17