writeRaster(r2,  names(list_raster)[[i]], bylayer = F, suffix='names',
format="GTiff", overwrite = T)
}
# rasterien summaaminen
# luetaan pisteytetyt rasteripinnat kansiosta
raster_files_uusi <- list.files("D:/HYMY/gis/demo/R/", pattern=".tif$")
list_raster_uusi <- lapply(raster_files_uusi, function(fl){
raster(paste0("D:/HYMY/gis/demo/R/", fl))
})
# summataan rasterit
rasteri_summattu<-do.call("sum",list_raster_uusi)
# normalisoidaan rasteri välille [0,1] jakamalla pikselin maksimiarvolla
tulosrasteri<-rasteri_summattu/maxValue(rasteri_summattu)
# plottaan lopullinen rasteripinta
plot(tulosrasteri)
# kirjoitetaan aineiston kansioon
writeRaster(tulosrasteri,  names(list_raster)[[i]], bylayer = F, suffix='names',
format="GTiff", overwrite = T)
writeRaster(tulosrasteri,  summarasteri.tif, bylayer = F, suffix='names',
format="GTiff", overwrite = T)
writeRaster(tulosrasteri,  "D:/HYMY/gis/demo/output/summarasteri.tif", bylayer = F, suffix='names',
format="GTiff", overwrite = T)
# summataan rasterit
rasteri_summattu<-do.call("sum",list_raster_uusi)
# normalisoidaan rasteri välille [0,1] jakamalla pikselin maksimiarvolla
tulosrasteri<-rasteri_summattu/maxValue(rasteri_summattu)
# plottaan lopullinen rasteripinta
plot(tulosrasteri)
# kirjoitetaan lopulliset tulosaineistot
# summattu rasteri
writeRaster(rasteri_summattu,  "D:/HYMY/gis/demo/output/summarasteri.tif", bylayer = F, suffix='names',
format="GTiff", overwrite = T)
# summattu ja normalisoitu rasteri
writeRaster(tulosrasteri,  "D:/HYMY/gis/demo/output/summarasteri_n.tif", bylayer = F, suffix='names',
format="GTiff", overwrite = T)
# plottaan lopullinen rasteripinta
plot(rasteri_summattu)
plot(tulosrasteri)
# plottaan lopullinen rasteripinta
plot(rasteri_summattu)
# pisteytetään aineistoa vastaavat rasterit
foreach (i = 1:length(list_raster)) %do%{
r1=list_raster[[i]] # valitsee listalta pisteytettävän rasterin
id=names(list_raster[[i]]) # muodostaa uuden [id] muuttujan rasterin nimestä
p1<-pisteet[pisteet$AineistoNimi == id,] # valitsee rasteria vastaavan pistemäärä taulukosta
# pisteyttää rasterin
r2<-(r1*(p1[,2]))
# kirjoittaa ja nimeää pisteytetyn rasterin ulos
writeRaster(r2,  names(list_raster)[[i]], bylayer = F, suffix='names',
format="GTiff", overwrite = T)
}
# luetaan pisteytetyt rasteripinnat kansiosta
raster_files_uusi <- list.files("D:/HYMY/gis/demo/R/", pattern=".tif$")
list_raster_uusi <- lapply(raster_files_uusi, function(fl){
raster(paste0("D:/HYMY/gis/demo/R/", fl))
})
# summataan rasterit
rasteri_summattu<-do.call("sum",list_raster_uusi)
# plottaan lopullinen rasteripinta
plot(rasteri_summattu)
raster_files <- list.files("D:/HYMY/gis/demo/input/", pattern=".tif$") # vaihda tiedostopolku kansioon, missä aineistot ovat
list_raster <- lapply(raster_files, function(fl){
raster(paste0("D:/HYMY/gis/demo/input/", fl)) # vaihda tiedostopolku kansioon, missä aineistot ovat
})
names(list_raster)
# nimetään rasterit samoin kuin kansiossa
names(list_raster) <- (raster_files)
names(list_raster)
# luetaan pistetaulukko
# pistetaulukon tulee sisältää kolumnit [AineistoNimi], [Pisteet]
# lisäksi taulukossa on hyvä olla lyhyt aineistoKuvaus, mutta tämä ei ole välttämätöntä mallin kannalta
pisteet<-read.table("D:/HYMY/gis/demo/pisteytysDEMO.csv", header = TRUE, sep =";", dec=".") # vaihda tiedostopolku kansioon
foreach (i = 1:length(list_raster)) %do%{
r1=list_raster[[i]] # valitsee listalta pisteytettävän rasterin
id=names(list_raster[[i]]) # muodostaa uuden [id] muuttujan rasterin nimestä
p1<-pisteet[pisteet$AineistoNimi == id,] # valitsee rasteria vastaavan pistemäärä taulukosta
# pisteyttää rasterin
r2<-(r1*(p1[,2]))
# kirjoittaa ja nimeää pisteytetyn rasterin ulos
writeRaster(r2, "D:/HYMY/gis/demo/testi/",paste0(names(list_raster))[[i]], bylayer = F, suffix='names',
format="GTiff", overwrite = T)
}
foreach (i = 1:length(list_raster)) %do%{
r1=list_raster[[i]] # valitsee listalta pisteytettävän rasterin
id=names(list_raster[[i]]) # muodostaa uuden [id] muuttujan rasterin nimestä
p1<-pisteet[pisteet$AineistoNimi == id,] # valitsee rasteria vastaavan pistemäärä taulukosta
# pisteyttää rasterin
r2<-(r1*(p1[,2]))
# kirjoittaa ja nimeää pisteytetyn rasterin ulos
writeRaster(r2, "D:/HYMY/gis/demo/testi/",names(list_raster)[[i]], bylayer = F, suffix='names',
format="GTiff", overwrite = T)
}
raster_files <- list.files("D:/HYMY/gis/demo/input/", pattern=".tif$") # vaihda tiedostopolku kansioon, missä aineistot sijaitsevat
list_raster <- lapply(raster_files, function(fl){
raster(paste0("D:/HYMY/gis/demo/input/", fl)) # vaihda tiedostopolku kansioon, missä aineistot ovat
})
# nimetään rasterit samoin kuin kansiossa
names(list_raster) <- (raster_files)
# luetaan pistetaulukko
# pistetaulukon tulee sisältää kolumnit [AineistoNimi], [Pisteet]
# lisäksi taulukossa on hyvä olla lyhyt aineistoKuvaus, mutta tämä ei ole välttämätöntä mallin kannalta
pisteet<-read.table("D:/HYMY/gis/demo/pisteytysDEMO.csv", header = TRUE, sep =";", dec=".") # vaihda tiedostopolku kansioon
# pisteytetään aineistoa vastaavat rasterit
foreach (i = 1:length(list_raster)) %do%{
r1=list_raster[[i]] # valitsee listalta pisteytettävän rasterin
id=names(list_raster[[i]]) # muodostaa uuden [id] muuttujan rasterin nimestä
p1<-pisteet[pisteet$AineistoNimi == id,] # valitsee rasteria vastaavan pistemäärä taulukosta
# pisteyttää rasterin
r2<-(r1*(p1[,2]))
# kirjoittaa ja nimeää pisteytetyn rasterin ulos
writeRaster(r2, "D:/HYMY/gis/demo/testi/",names(list_raster)[[i]], bylayer = F, suffix='names',
format="GTiff", overwrite = T)
}
foreach (i = 1:length(list_raster)) %do%{
r1=list_raster[[i]] # valitsee listalta pisteytettävän rasterin
id=names(list_raster[[i]]) # muodostaa uuden [id] muuttujan rasterin nimestä
p1<-pisteet[pisteet$AineistoNimi == id,] # valitsee rasteria vastaavan pistemäärä taulukosta
# pisteyttää rasterin
r2<-(r1*(p1[,2]))
# kirjoittaa ja nimeää pisteytetyn rasterin ulos
writeRaster(r2, file.path("D:/HYMY/gis/demo/testi/",names(list_raster)[[i]]), bylayer = F, suffix='names',
format="GTiff", overwrite = T)
}
foreach (i = 1:length(list_raster)) %do%{
r1=list_raster[[i]] # valitsee listalta pisteytettävän rasterin
id=names(list_raster[[i]]) # muodostaa uuden [id] muuttujan rasterin nimestä
p1<-pisteet[pisteet$AineistoNimi == id,] # valitsee rasteria vastaavan pistemäärä taulukosta
# pisteyttää rasterin
r2<-(r1*(p1[,2]))
# kirjoittaa ja nimeää pisteytetyn rasterin ulos
writeRaster(r2, names(list_raster)[[i]], bylayer = F, suffix='names',
format="GTiff", overwrite = T)
}
require(raster)
require(rgdal)
require(sp)
require(spdep)
require(foreach)
# luetaan rasterit kansiosta
raster_files <- list.files("D:/HYMY/gis/demo/input/", pattern=".tif$") # vaihda tiedostopolku kansioon, missä aineistot sijaitsevat
list_raster <- lapply(raster_files, function(fl){
raster(paste0("D:/HYMY/gis/demo/input/", fl)) # vaihda tiedostopolku kansioon, missä aineistot ovat
})
# nimetään rasterit samoin kuin kansiossa
names(list_raster) <- (raster_files)
# plotataan haluttu rasteri testiksi
plot(list_raster[[1]]) # [[1]] vastaa kansion ensimmäistä rasteria
# luetaan pistetaulukko
# pistetaulukon tulee sisältää kolumnit [AineistoNimi], [Pisteet]
# lisäksi taulukossa on hyvä olla lyhyt aineistoKuvaus, mutta tämä ei ole välttämätöntä mallin kannalta
pisteet<-read.table("D:/HYMY/gis/demo/pisteytysDEMO.csv", header = TRUE, sep =";", dec=".") # vaihda tiedostopolku kansioon
# pisteytetään aineistoa vastaavat rasterit
foreach (i = 1:length(list_raster)) %do%{
r1=list_raster[[i]] # valitsee listalta pisteytettävän rasterin
id=names(list_raster[[i]]) # muodostaa uuden [id] muuttujan rasterin nimestä
p1<-pisteet[pisteet$AineistoNimi == id,] # valitsee rasteria vastaavan pistemäärä taulukosta
# pisteyttää rasterin
r2<-(r1*(p1[,2]))
# kirjoittaa ja nimeää pisteytetyn rasterin ulos
writeRaster(r2, names(list_raster)[[i]], bylayer = F, suffix='names',
format="GTiff", overwrite = T)
}
pisteet
# pisteytetään aineistoa vastaavat rasterit
foreach (i = 1:length(list_raster)) %do%{
r1=list_raster[[i]] # valitsee listalta pisteytettävän rasterin
id=names(list_raster[[i]]) # muodostaa uuden [id] muuttujan rasterin nimestä
p1<-pisteet[pisteet$AineistoNimi == id,] # valitsee rasteria vastaavan pistemäärä taulukosta
# pisteyttää rasterin
r2<-(r1*(p1[,2]))
# kirjoittaa ja nimeää pisteytetyn rasterin ulos
writeRaster(r2, names(list_raster)[[i]], bylayer = F, suffix='names',
format="GTiff", overwrite = T)
}
raster_files_uusi <- list.files("D:/HYMY/gis/demo/R/", pattern=".tif$") # vaihda tiedostopolku kansioon, missä aineistot sijaitsevat
list_raster_uusi <- lapply(raster_files_uusi, function(fl){
raster(paste0("D:/HYMY/gis/demo/R/", fl)) # vaihda tiedostopolku kansioon, missä aineistot sijaitsevat
})
raster_files_uusi
# nimetään rasterit samoin kuin kansiossa
names(list_raster) <- (raster_files)
names(list_raster)
# luetaan pistetaulukko
# pistetaulukon tulee sisältää kolumnit [AineistoNimi], [Pisteet]
# lisäksi taulukossa on hyvä olla lyhyt aineistoKuvaus, mutta tämä ei ole välttämätöntä mallin kannalta
pisteet<-read.table("D:/HYMY/gis/demo/pisteytysDEMO.csv", header = TRUE, sep =";", dec=".") # vaihda tiedostopolku kansioon
# pisteytetään aineistoa vastaavat rasterit
foreach (i = 1:length(list_raster)) %do%{
r1=list_raster[[i]] # valitsee listalta pisteytettävän rasterin
id=names(list_raster[[i]]) # muodostaa uuden [id] muuttujan rasterin nimestä
p1<-pisteet[pisteet$AineistoNimi == id,] # valitsee rasteria vastaavan pistemäärä taulukosta
# pisteyttää rasterin
r2<-(r1*(p1[,2]))
# kirjoittaa ja nimeää pisteytetyn rasterin ulos
writeRaster(r2, names(list_raster)[[i]], bylayer = F, suffix='names',
format="GTiff", overwrite = T)
}
# summataan rasterit
rasteri_summattu<-do.call("sum",list_raster_uusi)
rasteri_summattu
rasteri_summattu<-do.call("sum",list_raster_uusi)
# normalisoidaan rasteri välille [0,1] jakamalla pikselin maksimiarvolla
tulosrasteri<-rasteri_summattu/maxValue(rasteri_summattu)
# plottaan lopullinen rasteripinta
plot(rasteri_summattu)
plot(tulosrasteri)
# kirjoitetaan lopulliset tulosaineistot
# summattu rasteri
writeRaster(rasteri_summattu,  "D:/HYMY/gis/demo/output/summarasteri.tif", bylayer = F, suffix='names',
format="GTiff", overwrite = T) # vaihda tiedostopolku kansioon, mihin aineistot kirjoitetaan
# summattu ja normalisoitu rasteri
writeRaster(tulosrasteri,  "D:/HYMY/gis/demo/output/summarasteri_n.tif", bylayer = F, suffix='names',
format="GTiff", overwrite = T) # vaihda tiedostopolku kansioon, mihin aineistot kirjoitetaan
# otetaan paketit käyttöön
require(raster)
require(rgdal)
require(sp)
require(spdep)
require(foreach)
# luetaan rasterit kansiosta
raster_files <- list.files("D:/HYMY/gis/demo/input/", pattern=".tif$") # vaihda tiedostopolku kansioon, missä aineistot sijaitsevat
list_raster <- lapply(raster_files, function(fl){
raster(paste0("D:/HYMY/gis/demo/input/", fl)) # vaihda tiedostopolku kansioon, missä aineistot ovat
})
# nimetään rasterit samoin kuin kansiossa
names(list_raster) <- (raster_files)
# plotataan haluttu rasteri testiksi
plot(list_raster[[1]]) # [[1]] vastaa kansion ensimmäistä rasteria
# luetaan pistetaulukko
# pistetaulukon tulee sisältää kolumnit [AineistoNimi], [Pisteet]
# lisäksi taulukossa on hyvä olla lyhyt aineistoKuvaus, mutta tämä ei ole välttämätöntä mallin kannalta
pisteet<-read.table("D:/HYMY/gis/demo/pisteytysDEMO.csv", header = TRUE, sep =";", dec=".") # vaihda tiedostopolku kansioon
# missä pistetaulukko sijaitsee
# pisteytetään aineistoa vastaavat rasterit
foreach (i = 1:length(list_raster)) %do%{
r1=list_raster[[i]] # valitsee listalta pisteytettävän rasterin
id=names(list_raster[[i]]) # muodostaa uuden [id] muuttujan rasterin nimestä
p1<-pisteet[pisteet$AineistoNimi == id,] # valitsee rasteria vastaavan pistemäärä taulukosta
# pisteyttää rasterin
r2<-(r1*(p1[,2]))
# kirjoittaa ja nimeää pisteytetyn rasterin ulos
writeRaster(r2, names(list_raster)[[i]], bylayer = F, suffix='names',
format="GTiff", overwrite = T)
}
length(list_raster)
require(raster)
require(rgdal)
require(sp)
require(spdep)
require(foreach)
raster_files <- list.files("D:/HYMY/gis/demo/input/", pattern=".tif$") # vaihda tiedostopolku kansioon, missä aineistot sijaitsevat
list_raster <- lapply(raster_files, function(fl){
raster(paste0("D:/HYMY/gis/demo/input/", fl)) # vaihda tiedostopolku kansioon, missä aineistot ovat
})
raster_files
names(list_raster)
# nimetään rasterit samoin kuin kansiossa
names(list_raster) <- (raster_files)
names(list_raster)
pisteet<-read.table("D:/HYMY/gis/demo/pisteytysDEMO.csv", header = TRUE, sep =";", dec=".") # vaihda tiedostopolku kansioon
# missä pistetaulukko sijaitsee
pisteet
# pisteytetään aineistoa vastaavat rasterit
foreach (i = 1:length(list_raster)) %do%{
r1=list_raster[[i]] # valitsee listalta pisteytettävän rasterin
id=names(list_raster[[i]]) # muodostaa uuden [id] muuttujan rasterin nimestä
p1<-pisteet[pisteet$AineistoNimi == id,] # valitsee rasteria vastaavan pistemäärä taulukosta
# pisteyttää rasterin
r2<-(r1*(p1[,2]))
# kirjoittaa ja nimeää pisteytetyn rasterin ulos
writeRaster(r2, names(list_raster)[[i]], bylayer = F, suffix='names',
format="GTiff", overwrite = T)
}
# pisteytetään aineistoa vastaavat rasterit
foreach (i = 1:length(list_raster)) %do%{
r1=list_raster[[i]] # valitsee listalta pisteytettävän rasterin
id=names(list_raster[[i]]) # muodostaa uuden [id] muuttujan rasterin nimestä
p1<-pisteet[pisteet$AineistoNimi == id,] # valitsee rasteria vastaavan pistemäärä taulukosta
# pisteyttää rasterin
r2<-(r1*(p1[,2]))
}
pisteet$AineistoNimi
(p1[,2]))
p1[,2]
names(r1)
foreach (i = 1:length(list_raster)) %do%{
r1=list_raster[[i]] # valitsee listalta pisteytettävän rasterin
id=names(r1) # muodostaa uuden [id] muuttujan rasterin nimestä
p1<-pisteet[pisteet$AineistoNimi == id,] # valitsee rasteria vastaavan pistemäärä taulukosta
}
foreach (i = 1:length(list_raster)) %do%{
r1=list_raster[[i]] # valitsee listalta pisteytettävän rasterin
id=names(r1) # muodostaa uuden [id] muuttujan rasterin nimestä
p1<-pisteet[pisteet$AineistoNimi == id,] # valitsee rasteria vastaavan pistemäärä taulukosta
# pisteyttää rasterin
r2<-(r1*(p1[,2]))
}
# pisteyttää rasterin
r2<-r1*(p1[,2])
r1
p1
p1[,2]
id
# nimetään rasterit samoin kuin kansiossa
names(list_raster) <- (raster_files)
names(list_raster)
r1
length(list_raster)
names(list_raster)
names(r1)
pisteet$AineistoNimi
# pisteytetään aineistoa vastaavat rasterit
foreach (i = 1:length(list_raster)) %do%{
r1=list_raster[[i]] # valitsee listalta pisteytettävän rasterin
id=names(list_raster)[[i]] # muodostaa uuden [id] muuttujan rasterin nimestä
p1<-pisteet[pisteet$AineistoNimi == id,] # valitsee rasteria vastaavan pistemäärä taulukosta
# pisteyttää rasterin
r2<-r1*(p1[,2])
# kirjoittaa ja nimeää pisteytetyn rasterin ulos
writeRaster(r2, names(list_raster)[[i]], bylayer = F, suffix='names',
format="GTiff", overwrite = T)
}
# rasterien summaaminen
# luetaan pisteytetyt rasteripinnat kansiosta
raster_files_uusi <- list.files("D:/HYMY/gis/demo/R/", pattern=".tif$") # vaihda tiedostopolku kansioon, missä aineistot sijaitsevat
list_raster_uusi <- lapply(raster_files_uusi, function(fl){
raster(paste0("D:/HYMY/gis/demo/R/", fl)) # vaihda tiedostopolku kansioon, missä aineistot sijaitsevat
})
# summataan rasterit
rasteri_summattu<-do.call("sum",list_raster_uusi)
# plottaan lopullinen rasteripinta
plot(rasteri_summattu)
# plottaan lopullinen rasteripinta
plot(rasteri_summattu)
plot(tulosrasteri)
# pisteytetään aineistoa vastaavat rasterit
foreach (i = 1:length(list_raster)) %do%{
r1=list_raster[[i]] # valitsee listalta pisteytettävän rasterin
id=names(list_raster)[[i]] # muodostaa uuden [id] muuttujan rasterin nimestä
p1<-pisteet[pisteet$AineistoNimi == id,] # valitsee rasteria vastaavan pistemäärä taulukosta
# pisteyttää rasterin
r2<-r1*(p1[,2])
# kirjoittaa ja nimeää pisteytetyn rasterin ulos
writeRaster(r2, file.path("D:/HYMY/gis/demo/testi",names(list_raster)[[i]]), bylayer = F, suffix='names',
format="GTiff", overwrite = T)
}
rm(list = ls())
lapply(paste('package:',names(sessionInfo()$otherPkgs),sep=""),detach,character.only=TRUE,unload=TRUE)
.rs.restartR()
require(raster)
require(rgdal)
require(sp)
require(spdep)
require(foreach)
raster_files <- list.files("D:/HYMY/gis/demo/input/", pattern=".tif$") # vaihda tiedostopolku kansioon, missä aineistot sijaitsevat
list_raster <- lapply(raster_files, function(fl){
raster(paste0("D:/HYMY/gis/demo/input/", fl)) # vaihda tiedostopolku kansioon, missä aineistot ovat
})
# nimetään rasterit samoin kuin kansiossa
names(list_raster) <- (raster_files)
# plotataan haluttu rasteri testiksi
plot(list_raster[[1]]) # [[1]] vastaa kansion ensimmäistä rasteria
# luetaan pistetaulukko
# pistetaulukon tulee sisältää kolumnit [AineistoNimi], [Pisteet]
# lisäksi taulukossa on hyvä olla lyhyt aineistoKuvaus, mutta tämä ei ole välttämätöntä mallin kannalta
pisteet<-read.table("D:/HYMY/gis/demo/pisteytysDEMO.csv", header = TRUE, sep =";", dec=".") # vaihda tiedostopolku kansioon
# pisteytetään aineistoa vastaavat rasterit
foreach (i = 1:length(list_raster)) %do%{
r1=list_raster[[i]] # valitsee listalta pisteytettävän rasterin
id=names(list_raster)[[i]] # muodostaa uuden [id] muuttujan rasterin nimestä
p1<-pisteet[pisteet$AineistoNimi == id,] # valitsee rasteria vastaavan pistemäärä taulukosta
# pisteyttää rasterin
r2<-r1*(p1[,2])
# kirjoittaa ja nimeää pisteytetyn rasterin ulos
writeRaster(r2, file.path("D:/HYMY/gis/demo/input/pisteytetyt/",names(list_raster)[[i]]), bylayer = F, suffix='names',
format="GTiff", overwrite = T)
}
# luetaan pisteytetyt rasteripinnat kansiosta
raster_files_uusi <- list.files("D:/HYMY/gis/demo/R/", pattern=".tif$") # vaihda tiedostopolku kansioon, missä aineistot sijaitsevat
list_raster_uusi <- lapply(raster_files_uusi, function(fl){
raster(paste0("D:/HYMY/gis/demo/R/", fl)) # vaihda tiedostopolku kansioon, missä aineistot sijaitsevat
})
# summataan rasterit
rasteri_summattu<-do.call("sum",list_raster_uusi)
# normalisoidaan rasteri välille [0,1] jakamalla pikselin maksimiarvolla
tulosrasteri<-rasteri_summattu/maxValue(rasteri_summattu)
# plottaan lopullinen rasteripinta
plot(rasteri_summattu)
plot(tulosrasteri)
# kirjoitetaan lopulliset tulosaineistot
# summattu rasteri
writeRaster(rasteri_summattu,  "D:/HYMY/gis/demo/output/summarasteri.tif", bylayer = F, suffix='names',
format="GTiff", overwrite = T) # vaihda tiedostopolku kansioon, mihin aineistot kirjoitetaan
# summattu ja normalisoitu rasteri
writeRaster(tulosrasteri,  "D:/HYMY/gis/demo/output/summarasteri_n.tif", bylayer = F, suffix='names',
format="GTiff", overwrite = T) # vaihda tiedostopolku kansioon, mihin aineistot kirjoitetaan
# luetaan pisteytetyt rasteripinnat kansiosta
raster_files_uusi <- list.files("D:/HYMY/gis/demo/input/pisteytetyt/", pattern=".tif$") # vaihda tiedostopolku kansioon, missä aineistot sijaitsevat
list_raster_uusi <- lapply(raster_files_uusi, function(fl){
raster(paste0("D:/HYMY/gis/demo/input/pisteytetyt/", fl)) # vaihda tiedostopolku kansioon, missä aineistot sijaitsevat
})
# summataan rasterit
rasteri_summattu<-do.call("sum",list_raster_uusi)
# normalisoidaan rasteri välille [0,1] jakamalla pikselin maksimiarvolla
tulosrasteri<-rasteri_summattu/maxValue(rasteri_summattu)
plot(rasteri_summattu)
plot(tulosrasteri)
# summattu rasteri
writeRaster(rasteri_summattu,  "D:/HYMY/gis/demo/output/summarasteri.tif", bylayer = F, suffix='names',
format="GTiff", overwrite = T) # vaihda tiedostopolku kansioon, mihin aineistot kirjoitetaan
# summattu ja normalisoitu rasteri
writeRaster(tulosrasteri,  "D:/HYMY/gis/demo/output/summarasteri_n.tif", bylayer = F, suffix='names',
format="GTiff", overwrite = T) # vaihda tiedostopolku kansioon, mihin aineistot kirjoitetaan
# luetaan rasterit kansiosta
raster_files <- list.files("D:/HYMY/gis/demo/input/", pattern=".tif$") # vaihda tiedostopolku kansioon, missä aineistot sijaitsevat
list_raster <- lapply(raster_files, function(fl){
raster(paste0("D:/HYMY/gis/demo/input/", fl)) # vaihda tiedostopolku kansioon, missä aineistot sijaitsevat
})
# nimetään rasterit samoin kuin kansiossa
names(list_raster) <- (raster_files)
# plotataan haluttu rasteri testiksi
plot(list_raster[[1]]) # [[1]] vastaa kansion ensimmäistä rasteria
# luetaan pistetaulukko
# pistetaulukon tulee sisältää kolumnit [AineistoNimi], [Pisteet]
# lisäksi taulukossa on hyvä olla lyhyt aineistoKuvaus, mutta tämä ei ole välttämätöntä mallin kannalta
pisteet<-read.table("D:/HYMY/gis/demo/pisteytysDEMO.csv", header = TRUE, sep =";", dec=".") # vaihda tiedostopolku kansioon
# missä pistetaulukko sijaitsee
# pisteytetään aineistoa vastaavat rasterit
foreach (i = 1:length(list_raster)) %do%{
r1=list_raster[[i]] # valitsee listalta pisteytettävän rasterin
id=names(list_raster)[[i]] # muodostaa uuden [id] muuttujan rasterin nimestä
p1<-pisteet[pisteet$AineistoNimi == id,] # valitsee rasteria vastaavan pistemäärä taulukosta
# pisteyttää rasterin kertomalla jokaisen solun taulukon pistemäärällä
r2<-r1*(p1[,2])
# nimeää rasterin ja exporttaa sen ulos
writeRaster(r2, file.path("D:/HYMY/gis/demo/pisteytetyt/",names(list_raster)[[i]]), bylayer = F, suffix='names',
format="GTiff", overwrite = T)
}
# rasterien summaaminen
# luetaan pisteytetyt rasteripinnat kansiosta
raster_files_uusi <- list.files("D:/HYMY/gis/demo/pisteytetyt/", pattern=".tif$") # vaihda tiedostopolku kansioon, missä aineistot sijaitsevat
list_raster_uusi <- lapply(raster_files_uusi, function(fl){
raster(paste0("D:/HYMY/gis/demo/pisteytetyt/", fl)) # vaihda tiedostopolku kansioon, missä aineistot sijaitsevat
})
# summataan rasterit
rasteri_summattu<-do.call("sum",list_raster_uusi)
# normalisoidaan rasteri välille [0,1] jakamalla pikselin maksimiarvolla
rasteri_summattu_n<-rasteri_summattu/maxValue(rasteri_summattu)
# plottaan lopullinen rasteripinta
plot(rasteri_summattu) # Pisteytyksien perusteella muodostettu loppukartta, joka kuvastaa viherrakenteen arvoja
# ekosysteemipalveluiden näkökulmasta. Pikselin korkeampi arvo tarkoittaa viherrakenteen kannalta tärkeämpää aluetta.
plot(rasteri_summattu_n) # Pisteytyksien perusteella muodostettu loppukartta, joka kuvastaa viherrakenteen arvoja ekosysteemipalveluiden näkökulmasta.
# Aineisto on normalisoitu välille [0-1]. Pikselin korkeampi arvo tarkoittaa viherrakenteen kannalta tärkeämpää aluetta.
# exportataan lopulliset tulosaineistot
# summattu rasteri
writeRaster(rasteri_summattu,  "D:/HYMY/gis/demo/output/summarasteri.tif", bylayer = F, suffix='names',
format="GTiff", overwrite = T) # vaihda tiedostopolku kansioon, mihin aineistot kirjoitetaan
# summattu ja normalisoitu rasteri
writeRaster(rasteri_summattu_n,  "D:/HYMY/gis/demo/output/summarasteri_norm.tif", bylayer = F, suffix='names',
format="GTiff", overwrite = T) # vaihda tiedostopolku kansioon, mihin aineistot kirjoitetaan
require(raster)
require(rgdal)
require(sp)
require(spdep)
require(foreach)
# luetaan rasterit kansiosta
raster_files <- list.files("D:/HYMY/gis/demo/input/", pattern=".tif$") # vaihda tiedostopolku kansioon, missä aineistot sijaitsevat
list_raster <- lapply(raster_files, function(fl){
raster(paste0("D:/HYMY/gis/demo/input/", fl)) # vaihda tiedostopolku kansioon, missä aineistot sijaitsevat
})
# nimetään rasterit samoin kuin kansiossa
names(list_raster) <- (raster_files)
# plotataan haluttu rasteri testiksi
plot(list_raster[[1]]) # [[1]] vastaa kansion ensimmäistä rasteria
# luetaan pistetaulukko
# pistetaulukon tulee sisältää kolumnit [AineistoNimi], [Pisteet]
# lisäksi taulukossa on hyvä olla lyhyt aineistoKuvaus, mutta tämä ei ole välttämätöntä mallin kannalta
pisteet<-read.table("D:/HYMY/gis/demo/pisteytysDEMO.csv", header = TRUE, sep =";", dec=".") # vaihda tiedostopolku kansioon
# missä pistetaulukko sijaitsee
# pisteytetään aineistoa vastaavat rasterit
foreach (i = 1:length(list_raster)) %do%{
r1=list_raster[[i]] # valitsee listalta pisteytettävän rasterin
id=names(list_raster)[[i]] # muodostaa uuden [id] muuttujan rasterin nimestä
p1<-pisteet[pisteet$AineistoNimi == id,] # valitsee rasteria vastaavan pistemäärä taulukosta
# pisteyttää rasterin kertomalla jokaisen solun taulukon pistemäärällä
r2<-r1*(p1[,2])
# nimeää rasterin ja exporttaa sen ulos
writeRaster(r2, file.path("D:/HYMY/gis/demo/pisteytetyt/",names(list_raster)[[i]]), bylayer = F, suffix='names',
format="GTiff", overwrite = T)
}
# rasterien summaaminen
# luetaan pisteytetyt rasteripinnat kansiosta
raster_files_uusi <- list.files("D:/HYMY/gis/demo/pisteytetyt/", pattern=".tif$") # vaihda tiedostopolku kansioon, missä aineistot sijaitsevat
list_raster_uusi <- lapply(raster_files_uusi, function(fl){
raster(paste0("D:/HYMY/gis/demo/pisteytetyt/", fl)) # vaihda tiedostopolku kansioon, missä aineistot sijaitsevat
})
# summataan rasterit
rasteri_summattu<-do.call("sum",list_raster_uusi)
# normalisoidaan rasteri välille [0,1] jakamalla pikselin maksimiarvolla
rasteri_summattu_n<-rasteri_summattu/maxValue(rasteri_summattu)
# plottaan lopullinen rasteripinta
plot(rasteri_summattu) # Pisteytyksien perusteella muodostettu loppukartta, joka kuvastaa viherrakenteen arvoja
# ekosysteemipalveluiden näkökulmasta. Pikselin korkeampi arvo tarkoittaa viherrakenteen kannalta tärkeämpää aluetta.
plot(rasteri_summattu_n) # Pisteytyksien perusteella muodostettu loppukartta, joka kuvastaa viherrakenteen arvoja ekosysteemipalveluiden näkökulmasta.
# Aineisto on normalisoitu välille [0-1]. Pikselin korkeampi arvo tarkoittaa viherrakenteen kannalta tärkeämpää aluetta.
# exportataan lopulliset tulosaineistot
# summattu rasteri
writeRaster(rasteri_summattu,  "D:/HYMY/gis/demo/output/summarasteri.tif", bylayer = F, suffix='names',
format="GTiff", overwrite = T) # vaihda tiedostopolku kansioon, mihin aineistot kirjoitetaan
# summattu ja normalisoitu rasteri
writeRaster(rasteri_summattu_n,  "D:/HYMY/gis/demo/output/summarasteri_norm.tif", bylayer = F, suffix='names',
format="GTiff", overwrite = T) # vaihda tiedostopolku kansioon, mihin aineistot kirjoitetaan
