Hallo Hans,
ich würde mit PROC GREMOVE und PROC GMAP arbeiten.
GREMOVE bietet die Möglichkeit, aus einer vorhandenen Koordinatendatei interne Grenzen zu beseitigen und so neue - aggregierte - Polygone zu schaffen.
Voraussetzung hierfür ist, das neben der ID-Variablen, die die vorhandenen Polygone definiert, eine weitere ID-Variable vorhanden ist
(oder mit Hilfe eines Daten-/Prozedurschrittes berechnet wird), die die gewünschte Aggregationsebene darstellt.
GREMOVE berechnet dann die neuen Polygone und speichert die neue Koordinatendatei in einem SAS-Dataset ab.
PROC GREMOVE data = alte_koordinaten out = neue_koordinaten;
ID alte_id;
BY neue_id;
run;
Beispiel:
PROC GREMOVE DATA=MAPS.GERMANY OUT=BUNDESLAENDER;
by DISTRICT;
id ID;
WHERE STATE in ('05','04');
RUN;
Im Anhang habe ich Dir ein Excel-Sheet zur Verfügung gestellt, da sind dem DISTRICT (PLZ) der Name der Stadt zugeordnet.
Diese Datei ich habe ich GERMANY_KREIS genannt. Kann dann im PROC GMAP verwendet werden z.B.
PROC GMAP data = GERMANY_KREIS
map = BUNDESLAENDER /* MAP = Koordinatengrundlage die verwendet werden soll */
......
Eine kleinere Einheit als PLZ habe ich nicht gefunden.
Viele Grüße
Dito
... View more
