Przeskocz nawigację.
Strona główna

Mapy w R

Tworzenie map w R
Żeby rozpocząć pracę z mapami trzeba się najpierw zaopatrzyć w ich wersje cyfrowe.
Za pomocą linku:
 
http://sdi.geoportal.gov.pl/wfs_prg/wfservice.aspx?request=getfeature&ve...
 
pobrać można mapę powiatów a zamieniając wartość parametru typename na WOJEWODZTWA
 
http://sdi.geoportal.gov.pl/wfs_prg/wfservice.aspx?request=getfeature&ve...
 
pobrać można mapę województw w formacie gml. Zgodnie z regulaminem, zasoby geoportalu używane mogą być jedynie do celów niekomercyjnych.
 
W programie R potrzebujemy formatu shapefile. W ramach tego formatu niezbędne są 3 pliki: shp (kontury obszarów), dbf (baza z danymi), shx (plik wiażący bazę z obszarami).
za pomocą formularza 
 
http://mygeodata.eu/apps/converter/main_EN.html?dataType=vector
 
w prosty sposób można przekonwertować mapy na pożądany format. Trzeba jedynie spakować pojedyncze pliki gml za pomocą zip-a. Następnie umieścić w formularzu i przesłać.
 
Do wyświetlania map służy przede wszystkim pakiet maps. Na początek warto zaznajomić się z funkcją readShapePoly
 
http://bm2.genes.nig.ac.jp/RGM2/R_current/library/maptools/man/readShape...
 
W dalszej części zakładam, że skrypt który będzie rysował mapy znajduje się w katalogu zawierającym pliki POWIATY.shp, POWIATY.shx, POWIATY.dbf, WOJ.shp, WOJ.shx, WOJ.dbf.
 
Rysowanie map województw
library("maptools")
powiaty <- readShapePoly("POWIATY.shp")
plot(powiaty)
 
Rysowanie map województw
library("maptools")
woj <- readShapePoly("WOJ.shp")
plot(woj)
 
Rysowanie map województw i powiatów
library("maptools")
powiaty <- readShapePoly("POWIATY.shp")
woj <- readShapePoly("WOJ.shp")
plot(powiaty) plot(woj, border="red", add=TRUE)

ZałącznikWielkość
POWIATY.gml3.18 MB
POWIATY.dbf9.71 KB
POWIATY.shp1.78 MB
POWIATY.shx3.07 KB
WOJEWODZTWA.gml658.47 KB
WOJ.dbf817 bajtów
WOJ.shp352.88 KB
WOJ.shx228 bajtów