Wojciech Sobala

Wojciech Sobala Redaktor
statystyczny,
biostatystyk,
Instytut Medycyny
Pr...

Temat: Mapa polski

Czy ktoś ma lub potrafi wskazać gdzie można dostać mapę Polski z podziałem na województwa.
Chodziło by mi o możliwość narysowania mapy z różnymi kolorami dla każdego województwa (wypełnienie tła).
Najlepiej jeżeli byłby to dane do rysowania w pakiecie maps.
Mogą to być również dane dotyczące granic województw aby można było je narysować przy pomocy funkcji poly.
Ostatni plan to wykorzystanie danych z formatu SVG. Takie dane posiadam więc potrzebna byłaby pomoc w zaimportowaniu do R tak aby można było skorzystać z poly do narysowania mapy.
Marcin S.

Marcin S. Kierownik ds
Raportowania i
Analiz
Windykacyjnych

Temat: Mapa polski

Tutaj znajdziesz województwa i powiaty w formacie shapefile (do obrobienia w R):

ftp://ump.waw.pl/Incoming/

Dostałem pytanie o prawa autorskie tej mapy. Znajduje się ona na stronie projektu ump.waw.pl i tam należy szukać takich informacji. Ja wykorzystuje ją do celów własnych (edukacyjnych), stąd nie do końca się tym interesowałem.

Inne źródło to:

http://www.geovertical.pl/pl/product_6.htm

Jest tam tylko podział na województwa i znów brak informacji o prawach autorskich. Tym razem jest to jednak strona konkretnej firmy i udostępnia ona te dane jako próbkę swojego produktu (link do mapki na ump.pl znalazłem w sieci i nie mogę znaleźć do niej odnośnika na ich stronie).Marcin S. edytował(a) ten post dnia 06.03.10 o godzinie 08:09

Temat: Mapa polski

Cześć!
A mógłbyś podać jakieś wskazówki jak te mapy wprowadzić pod R? Jak mi się wydaje trzeba tu korzystać z foprmatu ESRI, czy się mylę? Przepraszam za głupie pytania, ale nie mam kompletnie pojęcia jak po ściągnięciu mógłbym wykorzystać mapki w R-rze.... ;(
marcin
Wojciech Sobala

Wojciech Sobala Redaktor
statystyczny,
biostatystyk,
Instytut Medycyny
Pr...

Temat: Mapa polski

Potrzebna jest biblioteka maptools.
Poniżej program do narysowania mapki wraz z legendą i nazwami województw.

woj.shape = readShapePoly("wojew_region.shp",proj4string=CRS("+proj=longlat +ellps=clrk66")) #,IDvar="NAZ_WOJ"

xw = as.numeric(strsplit("14.71769 17.21543 20.19659 14.74992 16.21633 17.84389 22.32370 18.90744 20.35774 22.14644 15.37838 17.39268 18.71407 19.43922 19.98711 21.45352"," ")[[1]])
yw = as.numeric(strsplit("53.63893 54.23526 53.95697 52.37669 52.39657 53.20162 53.26125 51.68097 52.55559 51.37287 51.19397 50.72684 50.61751 50.00130 50.77654 50.08081"," ")[[1]])
nw = c("zachodnio-\npomorskie","pomorskie","warmińsko-\nmazurskie","lubuskie","wielkopolskie","kujawsko-\npomorskie","podlaskie","łódzkie","mazowieckie","lubelskie","dolnośląskie","opolskie","śląskie","małopolskie","świętokrzyskie","podkarpackie")
#
kolory <- c("#f4e3d7","#e9c6af","#d38d5f","#a05a2c")
ind <- c(3,2,2,2,1,3,1,3,1,3,1,4,2,3,2,3)
legtxt <- c("<14","15-24","25-49",">50")
par(mar=c(0,0,0,0))
plot(woj.shape, border="black",lwd=0.5, axes=T,col=kolory[ind])
text(xw,yw,nw,adj=0,cex=0.8)
legend(x=14, y=49.5, legend=legtxt, fill=kolory, bty="n",cex=0.7)

Temat: Mapa polski

Dziękuję.
Czy otwierałeś może ten plik przed zaimplementowaniem w Rrze? Tak z programu (nazwy zmiennych itp.) wnioskuję, że tak. Pan google mówi mi, ze to plik modelowania trójwymiarowego i można go otworzyć np. przy pomocy AutoCad-a, 3D Studio itp. A czy są może jakieś shareware-owe programy do tegoż służące? (chodzi o to, że przy komendzie: woj.shape = readShapePoly("wojew_region.shp",proj4string=CRS("+proj=longlat +ellps=clrk66")) #,IDvar="NAZ_WOJ" z poprawkami na miejsce zapisania pliku dostaję błąd w I linii "error opening shp file"...).
Wojciech Sobala

Wojciech Sobala Redaktor
statystyczny,
biostatystyk,
Instytut Medycyny
Pr...

Temat: Mapa polski

Marcin Pracz:
Dziękuję.
Czy otwierałeś może ten plik przed zaimplementowaniem w Rrze? Tak z programu (nazwy zmiennych itp.) wnioskuję, że tak. Pan google mówi mi, ze to plik modelowania trójwymiarowego i można go otworzyć np. przy pomocy AutoCad-a, 3D Studio itp. A czy są może jakieś shareware-owe programy do tegoż służące? (chodzi o to, że przy komendzie: woj.shape = readShapePoly("wojew_region.shp",proj4string=CRS("+proj=longlat +ellps=clrk66")) #,IDvar="NAZ_WOJ" z poprawkami na miejsce zapisania pliku dostaję błąd w I linii "error opening shp file"...).

Plikowi z roszerzeniem shp towarzyszą trzy inne pilki z roszerzeniami: dbf, prj, shx. Jeżeli ich nie w tym samym katalogu co plik z roszerzeniem shp to pojawia się błąd.
Programem o których słyszałem, że odczytuje takie pliki jest ArcView.

Temat: Mapa polski

Jeszcze raz dzięki. No cóż, są w tym samym katalogu także szukam dalej :)
Wojciech Sobala

Wojciech Sobala Redaktor
statystyczny,
biostatystyk,
Instytut Medycyny
Pr...

Temat: Mapa polski

Sprawdził bym czy są zainstalowane zależności do pakietu maptools, np. pakiet sp (najlepiej instalować komendą install.packages("maptools",dep=T))

Temat: Mapa polski

Dzięki. Wszystko już działa. Po prostu nie rozpakowałem katalogu i stąd to całe zamieszanie..... ;)

Temat: Mapa polski

Panie Wojciechu,
czy mógłbym jeszcze mieć pytanie do programu, który Pan zamieścił? Nie rozumiem na jakiej zasadzie działa tu kolejność kolorowania. Województwa są kolorowane zgodnie ze zmienną "ind", ale kolejność wydaje się być zupełnia inna od tej, której Pan użył przy wpisywaniu nazw (czyli zachodnio-pomorskie, pomorskie, warmińsko-mazurskie, itd.) Być może ma to związek z tym parametrem [[1]], którego też nie rozumiem niestety...
Dawid Pyrdał

Dawid Pyrdał Specjalista ds. GIS

Temat: Mapa polski

Do pracy z plikami formatu *.shp polecam wykorzystanie programu z rodziny Open Source o nazwie QuantumGIS. Jest to program, który za darmo można ściągnąć np. ze strony http://linfiniti.com/downloads/QGIS-1.4.0-1-No-GrassSe...

Można w nim zarówno obejrzeć zawartość plików *.shp jak również prowadzić na nim różne analizy przestrzenne. Jest również możliwość korzystania z funkcji "R"-a po zainstalowaniu specjalnej wtyczki.
Wojciech Sobala

Wojciech Sobala Redaktor
statystyczny,
biostatystyk,
Instytut Medycyny
Pr...

Temat: Mapa polski

Parametr [[1]] wynika z odwołania do pierwszego elementu listy która jest zwracana przez funkcję strsplit, a potem konwertowana na wektor z wartościami numerycznymi.
xw i yw to współrzędne X oraz Y dla nazw województw (1 do 16). Kolejność tych województw nie ma nic wspólnego z kolejnością granic obszarów. Kolejność współrzędnych jest taka sama jak nazw (nw).
W zmiennej ind kolejność elementów odpowiada kolejności granic obszarów w zbiorze danych. Wartości tego wektora służą do odnajdowania indeksów w wektorze z kodami kolorów (nazwanym kolory).

Do generowania nazw na wykresie można użyć nazw, które występują w zbiorze z granicami ale ja chciałem niektóre z nich napisać w dwóch liniach więc skorzystałem z osobnego wektora z odpowiednio sformatowanymi nazwami ( w tym przypadku znak \n powoduje wstawienie znaku nowej lini).
Michał Bojanowski

Michał Bojanowski socjolog, analityk

Temat: Mapa polski

W wolnej chwili zacząłem sklecać pakiet do rysowania map Polski. Link: http://r-forge.r-project.org/projects/mapoland/

Na razie zawiera on tylko jedną mapę dla województw o stosunkowo niskiej rozdzielczości i "nieGISową" (tzn. współrzędne poligonów nie mają żadnej specjalnej interpretacji). Wystarczy jednak do celów prezentacyjnych.

Mam już w miarę przygotowane województwa w wyższej rozdzielczości oraz powiaty. Jak znajdę trochę więcej czasu to je do pakietu dołożę.

Jeżeli ktoś chce się włączyć w pisanie to zapraszam!

~michal

PS. No i oczywiście mile widziane sugestie i informacje o ew. błędachMichał Bojanowski edytował(a) ten post dnia 25.03.10 o godzinie 21:12
Wojciech Sobala

Wojciech Sobala Redaktor
statystyczny,
biostatystyk,
Instytut Medycyny
Pr...

Temat: Mapa polski

Na testy + sugestie możesz liczyć.

Temat: Mapa polski

W porządku, dziękuję za wyjaśnienia. Wszystko co dotyczy nanoszenia nazw województw na mapę bez problemu rozumiem. Chodziło mi tylko o kolejność granic, które jak z tego wynika są ustalone w pliku źródłowym.
Marcin Szymkowiak

Marcin Szymkowiak Specjalista w
Ośrodku Statystyki
Małych Obszarów,
Urząd S...

Temat: Mapa polski

Do robienia wykresów mapowych polecam pliki shapefile, które można ściągnąć ze strony http://www.gadm.org/country. Znajdują się tam pliki do ściągnięcia dla bardzo wielu państw. W przypadku Polski mamy układ województw łącznie z powiatami. Pozdrawiam
Marek R.

Marek R. Software Architect

Temat: Mapa polski

Zapraszam wszystkich do korzystania z darmowego narzędzia do prezentacji danych na mapie:
http://mapy.easysoftware.pl
dane można nanosić na podstawie kodów pocztowych lub kodów TERYT z dokładnością do gmin

Następna dyskusja:

mapa Polski




Wyślij zaproszenie do