1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106
| proc sort data = bd;
by commune;
run;
data bdr;
length ID $15.;
set bd;
ID='FR-'||Departement;
informat ID $15.;
format ID $15.;
run;
data nombre_commune;
set bdr;
by commune;
if first.commune then n=1;
else n+1;
if last.commune;
run;
/* Construction de la carte */
proc datasets library =work nolist;
modify nombre_commune;
rename latitude=lat;
rename longitude=long;
quit;
data ville;
set nombre_commune;
keep commune lat Long ID Classe_simple n;
run;
data villee;
set ville;
if lat = -361.1000022 then delete;
run;
data village; set villee;
length function $8 color $8;
xsys='2';
ysys='2';
hsys='3';
when='a';
function='pie'; style='psolid';
if classe_simple='A4' then color='yellow';
if classe_simple='A1' then color='red';
if classe_simple='A2' then color='blue';
if classe_simple='A3' then color='green';
if classe_simple='A5' then color='gray';
rotate=360;
size=1;
run;
proc gproject data=village out=villageois degrees NORANGECHECK latlon eastlong dupok;
id ID;
run;
proc gproject data=nombre_commune out=nombre_communes degrees NORANGECHECK latlon eastlong dupok;
id ID;
run;
data villageoise;
set villageois (obs=5);
run;
PROC SQL;
CREATE VIEW WORK.MAPCHARTMAPPREP AS
SELECT *
FROM MAPSgfk.FRANCE
;
QUIT;
PROC SQL;
CREATE VIEW WORK.MAPCHARTRESPONSEPREP AS
SELECT *
FROM WORK.nombre_commune
;
QUIT;
GOPTIONS CBACK= ;
TITLE;FOOTNOTE;
TITLE1 "Carte";
pattern1 v=s c=white;
PROC GMAP GOUT=MAPCHART DATA=WORK.MAPCHARTRESPONSEPREP MAP=WORK.MAPCHARTMAPPREP ALL anno = villageoise;
ID ID;
CHORO ID/
WOUTLINE=1
legend = legendnn
;
RUN;
QUIT; |
Partager