> ## Estimer postnummersoner basert på koordinatar til postnummer
>
> library(rgdal) # For lasting av kartfiler
> library(rgeos) # For geometriske operasjonar
> library(plyr) # For enkel datamanipulering
> library(spatstat) # For Voronoi-tesselering og formatkonvertering
> library(spdep) # For utrekning av polygon-naboskap
> library(maptools) # For formatkonvertering og ymse
> library(RColorBrewer) # For fine fargepalettar
>
>
>
> # Last inn postnummerdatabasen. Last først alt som tekst
> # (for å bevara leiane nullar i postnummer), og gjer så
> # lengde- og breiddegradar til tal.
> pnr=read.delim("postnummer.csv", colClasses="character")
> pnr=transform(pnr, LON=as.numeric(LON), LAT=as.numeric(LAT))
>
> # Fjern postnummer som ikkje er i bruk
> pnr.ok=pnr[!grepl("ikkje i bruk", pnr$POSTSTAD),]
>
> # Nokre posisjonar er dei same for ulike postnummer. Dette
> # fører til problem, så slå dei derfor saman til eitt nummer.
> pnr.unik=ddply(pnr.ok, ~LON+LAT, summarise, POSTNR=paste(POSTNR, collapse=", "))
>
>
> uproj=CRS("+init=epsg:32633") # Uprojiserte data
>
> proj4string(pnr.unik)=uproj
Error in function (classes, fdef, mtable) :
unable to find an inherited method for function "proj4string pnr.proj=spTransform(pnr.unik,uproj) # Transformer til UTM sone 33
Error in function (classes, fdef, mtable) :
unable to find an inherited method for function "spTransform", for signature "data.frame", "CRS"
>

For berre å gjera om koordinatane til UTM-koordinatar treng du ikkje lasta inn kartet i det heile. Berre bruk denne proj4string-en i spTransform()-funksjonen: CRS("+init=epsg:32633"). Dette er koden for UTM sone 33, som ofte er brukt for å representera heile Norge. Byt ut dei to siste siffera for andre sonar. For delar av Norge er andre projeksjonar betre, frå UTM sone 31 i vest til UTM sone 36 heilt aust (i Finnmark), pluss andre projeksjonar for endå meir lokale forhold.

Jeg skjønte ikke syntaksen på
fylke=readOGR(«N5000 shape», «N5000_ArealdekkeFlate»)
men fant ut etterhvert at første felt er katalogen til fila. Da gikk alt bra tilslutt.

Det jeg primært trengte var å generer en CSV-fil med UTM-xy koordinater til postnummer. Det gikk etter litt plundring.

Takk for svar!

> # Vis kor mange færre «postnummer» me no fekk
> nrow(pnr)
[1] 4665
> nrow(pnr.unik)
[1] 4331
>
>
> # Norgeskartet er i UTM-projeksjon sone 33, mens postnummera
> # er uprojiserte (dvs. lengde- og breiddegradar). Gjer derfor
> # koordinatane til postnummera om til UTM sone 33-koordinatar.
> proj=CRS(proj4string(norge)) # UTM sone 33
Error in proj4string(norge) :
error in evaluating the argument 'obj' in selecting a method for function 'proj4string': Error: object 'norge' not found

No har eg endeleg lært meg å generere filer, så det er berre å kome med ynske om andre format.

Det er for øvrig ein liten feil i filformatet. CSV-fila er ikkje ei ekte CSV-fil; ho er ei TSV-fil. Skildringa «Tabseparert CSV» gjev ikkje heilt meining, då C-en i CSV nettopp står for «comma». Sjølv om nokre folk «misbrukar» namnet CSV til å tyda DSV, bør ein unngå dette. Det fører til forvirring.

Postnummerlistene er tilgjengelege i ei rekkje ulike format, ikkje berre HTML, men CSV, TXT, XML og KML.

Sjekk ut http://www.erikbolstad.no/geo/noreg/postnummer/last-ned/.

Ha ein fin dag!
Erik