|
|
| Rivi 4: |
Rivi 4: |
|
| |
|
| <rcode embed=1 label="Näytä pyörätelineet" showcode=-1> | | <rcode embed=1 label="Näytä pyörätelineet" showcode=-1> |
|
| |
|
| |
| # X.rakennustunnus kayttotarkoitus kokonaisala tilavuus ika valmistumisvuosi paaasiallinen_lammitystapa x_koord y_koord
| |
| #1 29700200180021000.00 039 Muut kerrostalot 3962 11800 14 1998 Vesikeskuslämmitys 6973868.64 535256.24
| |
| #2 29700200180022000.00 039 Muut kerrostalot 6685 22110 68 1944 Vesikeskuslämmitys 6973928.61 535228.24
| |
| #3 29700200200006004.00 162 Kulkuneuvojen suoja- ja huoltorakennukset 120 360 4 2008 Ei kiinteää lämmityslaitetta 6973854.64 535102.30
| |
| #4 29700200200006008.00 039 Muut kerrostalot 4735 15380 1 2011 Vesikeskuslämmitys 6973850.64 535120.29
| |
| #5 29700200200006008.00 012 Kahden asunnon talot 213 790 1 2011 Vesikeskuslämmitys 6973826.66 535145.27
| |
| #6 29700200200007008.00 711 Teollisuusvarastot 501 3900 39 1973 Ei kiinteää lämmityslaitetta 6973808.66 535113.29
| |
|
| |
|
| library(OpasnetUtils) | | library(OpasnetUtils) |
| Rivi 22: |
Rivi 13: |
| oprint(dat) | | oprint(dat) |
|
| |
|
| par(mfrow=c(6,1), mar=c(3,1,0,1), cex=1.5)
| | coordinates(out)=c("LO","LA") |
| BOTTOM <- 1
| |
| | |
| par
| |
| | |
| colorstrip <- function(colors, labels)
| |
| {
| |
| count <- length(colors)
| |
| m <- matrix(1:count, count, 1)
| |
| image(m, col=colors, ylab="", axes=FALSE)
| |
| axis(BOTTOM,approx(c(0, 1), n=length(labels))$y, labels)
| |
| }
| |
| | |
| #my_csv <- opasnet.csv("f/fe/Kuo_rak.csv", "heande", sep = ",", header = TRUE)
| |
| | |
| # Create new classes depending on floor area
| |
| | |
| my_csv$class <- cut(my_csv$kokonaisala, breaks = c(0, 2, 5, 10, 20, 50), include.lowest = TRUE, labels = FALSE)
| |
| | |
| coordinates(my_csv) <- c("LA", "LO")
| |
| proj4string(my_csv) <- ("+init=epsg:3067")
| |
| epsg4326String <- CRS("+proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs")
| |
| shp<-spTransform(my_csv,epsg4326String)
| |
| | |
| coordinates
| |
| proj4string
| |
| spTransform
| |
| | |
| #Create blank raster
| |
| rast<-raster()
| |
| | |
| #Set raster extent to that of point data
| |
| extent(rast)<-extent(shp)
| |
| | |
| #Choose number of columns and rows
| |
| ncol(rast)<-200
| |
| nrow(rast)<-300
| |
| | |
| | |
| cat("Pyörätelineet\n")
| |
| | |
| #Rasterize point data
| |
| rast2<-rasterize(shp, rast, shp$class, fun=mean)
| |
|
| |
|
| colors <- rainbow(5)
| | out <- google.point_kml( |
| colorstrip(colors, c('0-100','100-200','200-500','500-1000','> 1000'))
| | obj = dat, |
| | kmlname = "Helsingin pyörätelineet", |
| | kmldescription = "KML-tiedostona", |
| | name = dat$Teline, |
| | description = paste("<b>Year:</b>",dat$Teline) |
| | ) |
|
| |
|
| #Plot data
| | google.show_kml_data_on_maps(out) |
| google.show_raster_on_maps(rast2, col=colors, style="height:600px;") | |
|
| |
|
| </rcode> | | </rcode> |