如何使用sas map 绘制地图

如何使用sasmap繪製地圖以台灣地圖為例想一想如何作畫材料1.一張台灣地圖2.一份各地區觀測值資料工具

GMAPprocedure表示方式---平面圖或長方圖Sasmapdatasetsdefaultsaslibrary:MapsTaiwanTaiwan2縣市為最小單位Taiwan:Mapdata-set座標Taiwan2:Responsedataset才能畫界線以鄉鎮為最小單位GMAPprocedureUsethePROCGMAPstatementtoidentifythemapandresponsedatasets.UsetheIDstatementtonametheidentificationvariable(s).UseaBLOCK,CHORO,PRISM,orSURFACEstatementtoidentifytheresponsevariableandgeneratethemapProcedureSyntaxPROCGMAP

MAP=map-data-set

<DATA=response-data-set>

<ALL>

<ANNOTATE=Annotate-data-set>

<GOUT=<libref.>output-catalog>

<IMAGEMAP=output-data-set>

IDid-variable(s);

CHORO

response-variable(s)</option(s)>;

BLOCK

response-variable(s)</option(s)>;

PRISM

response-variable(s)</option(s)>;SURFACE

response-variable(s)</option(s)>;ProducingaSimpleSurfaceMapBasicprocedureSimpleandEasychoroCc

gmapdata=maps.taiwan2map=maps.taiwan

idid;

choro

idname/discrete

coutline=black;run;ResponsedatasetMapdata-setResponsevar.?blockCreatesthree-dimensionalblockmapsonwhichlevelsofmagnitudeofthespecifiedresponsevariablesarerepresentedbyblocksofvaryingheight,pattern,andcolor.

datamap;setmaps.taiwan2;ifidin(5

10

15

20)thenno=5;elseifidin(3

6

9

12

15

18)thenno=3;elseifidin(7

14

21)thenno=7;elseno=0;run;

proc

gmapdata=mapmap=maps.taiwanidid;blockno/discretecoutline=blackCBLKOUT=black;run;prismCreatesthree-dimensionalprismmapsinwhichlevelsofmagnitudeofthespecifiedresponsevariablesarerepresentedbypolyhedrons(raisedpolygons)ofvaryingheight,pattern,andcolor.

proc

gmapdata=mapmap=maps.taiwaidid;prismno/discretecoutline=black;run;Advanceprocedure---

toenhancethemapComplicatedandOptionallegendAnnotateByvarpatternexamplegoptionsreset=all;legend1position=(middlerightinside)label=none shape=bar(1,1)pct mode=share across=1;Datacitystar;……;run;proc

gmapdata=strata8map=d.taiwan;bystrata;idid;chorocity1city/discrete

coutline=black

legend=legend1

annotate=citystar;formatcity1citycity.stratast.;run;MethodGoptions---SetthegraphicsenvironmentPatternLegendAnnotatedatasetChorostatement</options>--appearanceoptions--legendoptions--mappingoptionsGoptions---

setthegraphicsenvironmentoptionsthataffecttheappearanceofthedisplayareaandthegraphicsoutputoptionsthataffectcolor

optionsthatcontrolfontselectionortextappearancegoptions

reset=global

gunit=pctborder

cback=whitecolors=(blackbluegreenred)

ctext=black

ftext=swisshtitle=6

htext=3;

patternpattern1v=msc=LIB;pattern2v=msc=LIBG;V=map/plot-pattern

PlotpatternlegendTheLEGENDstatementcontrolsthelocationandappearanceoflegendsontwo-dimensionalplots,contourplots,maps,andcharts.

legend1position=(middlerightinside)label=none shape=bar(1,1)pct mode=share across=1;Annotatedatasetdata

citystar;lengthfunctionstylecolor$8position$1text$20;retainfunction'label'xsysysys'2'hsys'3'when'a';setxy(keep=ididnamexy);ifidin(3051101111201630710154011327120312141510

1003

1702

619);/*city為縣市代碼*/

ifid>999thencity=substr(id,9,2)*1;

elseifid>99thencity=substr(id,10,1)*1;

elsecity=substr(id,11,2)*1;

ifid^=0;/*strata分層*/

if(idin(63

101

102

103

104

105

106

110

111))or(cityin(63))thenstrata=1;if(idin(107

108

109

112

113

114

115

116

117

118

119

120

121

122

123

124

125

126

127

128

129

17))or(cityin(17))thenstrata=2;if(cityin(3

4

5

18))thenstrata=3;if(cityin(6

7

8

19))thenstrata=4;if(cityin(9

20

10

21

11))thenstrata=5;if(cityin(64

12

13

16))thenstrata=6;if(cityin(2

14

15))thenstrata=7;style='special';text='M';color='red';size=7;position='5';output;style='KAI';text=idname;color='yellow';size=5;position='8';output;run;Chorostatement</options>;proc

gmapdata=strata8map=d.taiwan;bystrata;idid;chorocity1/discretecoutline=black

legend=legend1

annotate=citystar;formatcity1city.stratast.;run;appearanceoptionsANNOTATE=Annotate-data-setCEMPTY=empty-area-outline-color