Vous n'êtes pas identifié.
Vous êtes sur le forum du master ESA !
Le site du master ESA - description de la formation, notes de cours, contacts... vient de déménager !!!
Venez visiter notre nouveau site : www.master-esa.fr
Bonjour à tous !
J'ai fait une petite découverte en fouinant dans la doc SAS.. du coté SAS/GRAPH, et donc je suis tombé sur la PROC GMAP qui permet de s'amuser avec les .shp ( cartes ).
Tout d'abord il faut se munir d'un .shp ( exemple : le « south.shp » dans les samples fournis avec GeoDa )
On commence par une PROC MAPIMPORT pour charger la carte dans SAS :
proc mapimport out=south datafile='C:\Program Files\GeoDa\Sample Data\south.shp' ; run;
NOTE: 41163 observations copiées de la table WORK.SOUTH.
NOTE: La table WORK.SOUTH a 41163 observations et 72 variables.
NOTE: La procédure SORT a utilisé (Durée totale du processus) :
temps réel 3.67 secondes
temps processeur 0.26 secondes
Pas bon !! ça le fait avec toutes les cartes je ne sais pas pourquoi...
astuce :
proc sort data=south ;by fipsno;run; data south1 ; set south; by fipsno ; if first.fipsno then output;run;
NOTE: 41163 observations copiées de la table WORK.SOUTH.
NOTE: La table WORK.SOUTH1 a 1412 observations et 72 variables.
NOTE: L'étape DATA a utilisé (Durée totale du processus) :
temps réel 1.15 secondes
temps processeur 0.04 secondes
On dispose d'un SOUTH1.WORK avec 1412 obs ( plus logique car il n'y a que 1412 contés dans le sud des usa )
Et pour afficher... :
proc gmap data=south1 map=south; id fipsno; choro HR80 / levels=5 coutline=black woutline=3 ; run; quit;
id : clé primaire de la table (unique)
choro : pour afficher la carte en couleur selon les niveaux de la variable HR80
on peut choisir l'option « PRISM » au lieu de choro pour afficher la carte en... prism lol ( ça fait une carte en 3D )
« BLOCK » pour avoir des barres verticales..
ex :
Block HR80 / levels=5 shape=cylinder
Dans le cas des variables discrètes on rajoute l'option « discrete » :
title1 'elections 2000'; footnote1 'usa'; pattern1 value=solid color=lime; pattern2 value=solid color=cyan; proc gmap data=elections map=usa; id fipsno; choro winner / discrete; run; quit;
...etc ! ça marche à peu prêt comme une PROC GPLOT.
voili voualou ! Plus d'infos « The GMAP Procedure » dans l'aide SAS.
oulaaa 15h25... sauve-qui-peut !
Hors ligne
arf !! désolé je me suis tromper de rubrique ! je voulais mettre ça dans trucs et astuces !
Hors ligne