|
|
| (28 välissä olevaa versiota 3 käyttäjän tekeminä ei näytetä) |
| Rivi 1: |
Rivi 1: |
| [[Luokka:Minera]]
| | #REDIRECT [[Sedimenttien metallipitoisuuksien arviointi]] |
| {{edistymisluokitus|edistyminen = Luonnos}}
| |
| {{muuttuja|moderator=Marjo}}
| |
| | |
| == Kysymys ==
| |
| | |
| Miten lasketaan suomalaisen järven sedimentin metallipitoisuus?
| |
| | |
| == Vastaus ==
| |
| | |
| === 1) Lasketaan järven sedimentin metallipitoisuus perustuen ympäröivän, ympyrän säteen määräämän alueen järvien keskipitoisuuteen. ===
| |
| | |
| <rcode graphics="1" include="page:OpasnetBaseUtils|name:generic"
| |
| variables="
| |
| name:metalli|description:Valitse metallit|type:checkbox|options:'Co';Co;'Cr';Cr;'Cu';Cu;'Fe';Fe;'Mn';Mn;'Ni';Ni;'Pb';Pb;'S';S;'V';V;'Zn';Zn|
| |
| name:northing|description:Northing-koordinaatti (KKJ)|default:6995650|
| |
| name:easting|description:Easting-koordinaatti (KKJ)|default:3653220|
| |
| name:säde|description:Määritä säde (km; järvet säteen määräämän ympyränmuotoisen alueen sisäpuolella huomioidaan)|default:100
| |
| ">
| |
| | |
| library(OpasnetBaseUtils)
| |
| library(xtable)
| |
| library(ggplot2)
| |
| library(reshape)
| |
| | |
| data <- (op_baseGetData("opasnet_base", "Op_en5462"))
| |
| data <- tidy(data)
| |
| | |
| data$Result <- as.numeric(data$Result) # This replaces all as.numeric rows below.
| |
| data <- reshape(data, timevar = "Observation", idvar = "ID", direction = "wide")
| |
| colnames(data) <- gsub("Result.", "", colnames(data))
| |
| | |
| data$distance <- sqrt((as.numeric(data$Easting) - easting)^2 + (as.numeric(data$Northing) - northing)^2)
| |
| | |
| data <- melt(data, measure.vars = c("Co", "Cr", "Cu", "Fe", "Mn", "Ni", "Pb", "S", "V", "Zn"))
| |
| | |
| out <- data
| |
| out$distance <- data$distance <= säde*1000
| |
| out <- as.data.frame(as.table(tapply(out$value, out[c("variable", "distance")], mean)))
| |
| colnames(out) <- gsub("variable", "metalli", colnames(out))
| |
| colnames(out) <- gsub("Freq", "pitoisuus", colnames(out))
| |
| cat("Metallipitoisuuden keskiarvo ", säde, " km säteellä pisteestä ", easting, "E, ", northing, " N.\n")
| |
| print(xtable(out[out$distance == TRUE, ]), type = 'html')
| |
| | |
| cat("Pitoisuudet metallille", metalli, ".\n")
| |
| concentration <- data[data$distance <= säde*1000 & data$variable == metalli, "value"]
| |
| plot(concentration) # simple plot about concentrations.
| |
| qplot(Easting, Northing, data = data[data$variable == metalli, ], size = value) # Concentrations on a "map".
| |
| </rcode>
| |
| | |
| === 2) Lasketaan järven sedimentin metallipitoisuus läheisen/sopivan suuren järven pitoisuuden perusteella. ===
| |
| | |
| <rcode graphics="1" include="page:OpasnetBaseUtils|name:generic|page:Olio-ohjelmointi_Opasnetissä|name:answer"
| |
| variables="
| |
| name:metalli|description:Valitse metallit|type:checkbox|options:'Co';Co;'Cr';Cr;'Cu';Cu;'Fe';Fe;'Mn';Mn;'Ni';Ni;'Pb';Pb;'S';S;'V';V;'Zn';Zn|
| |
| name:Järvi|description:Valitse järvi|type:checkbox|options:'Pielinen';Pielinen, Lieksa;'Oulujärvi';Oulujärvi, Vaala;'Kallavesi;Kallavesi, Kuopio;'Suvasvesi';Suvasvesi, Vehmersalmi;'Yli-Kitka;Ylikitka, Kuusamo;'Nilakka';Nilakka, Keitele;'Koitere';Koitere, Ilomantsi;'Kivijärvi';Kivijärvi, Kivijärvi;'Viinijärvi';Viinijärvi, Polvijärvi;'Keurusselkä';Keurusselkä, Keuruu;'Ontojärvi';Ontojärvi, Kuhmo;'Vuohijärvi';Vuohijärvi, Vaala;'Kukkia';Kukkia, Luopioinen;'Kostonjärvi';Kostonjärvi, Taivalkoski;'Sääksjärvi';Sääksjärvi, Kokemäki;'Pieksänjärvi';Pieksänjärvi, Pieksämäki;'Pihlajavesi';Pihlajavesi, Keuruu;'Keyritty';Keyritty, Rautavaara;'Nerkoonjärvi';Nerkoonjärvi, Lapinlahti;'Pyhäjärvi';Pyhäjärvi, Artjärvi;'Uuranjärvi';Uuranjärvi, Multia;'Sääjärvi';Sääjärvi, Janakkala;'Haukijärvi';Haukijärvi, Maaninka;'Kinnasjärvi';Kinnasjärvi, Tuupovaara;'Tervajärvi';Tervajärvi, Luumäki;'Kalliojärvi';Kalliojärvi, Juuka;'Pyhäluoma';Pyhäluoma, Pieksämäki;'Puntasjärvi';Puntasjärvi, Kuru;'Ylimmäinen';Ylimmäinen, Jaala;'Hanhijärvi';Hanhijärvi, Jämsänkoski;'Ruuhilampi';Ruuhilampi, Virtasalmi;'Ala-Saimaa';Ala-Saimaa, Taipalsaari;'Pihlajavesi';Pihlajavesi, Puumala;'Orivesi';Orivesi, Kitee;'Puulavesi';Puulavesi, Hirvensalmi;'Höytiäinen';Höytiäinen, Kontiolahti;'Pyhäjärvi';Pyhäjärvi, Säkylä;'Kiantajärvi';Kiantajärvi, Suomussalmi;'Lappajärvi';Lappajärvi, Lappajärvi;'Pyhäjärvi';Pyhäjärvi, Pyhäsalmi;'Vesijärvi';Vesijärvi, Asikkala;'Unnukka';Unnukka, Leppävirta;'Kolimajärvi';Kolimajärvi, Viitasaari;'Lestijärvi';Lestijärvi, Lestijärvi;'Muuratjärvi';Muuratjärvi, Muurame;'Karijärvi';Karijärvi, Jaala;'Kolkonjärvi';Kolkonjärvi, Rantasalmi;'Ätäskö';Ätäskö, Kitee;'Nurmijärvi';Nurmijärvi, Rautjärvi;'Kevätön';Kevätön, Siilinjärvi;'Kotajärvi';Kotajärvi, Lemi;'Ihalanjärvi';Ihalanjärvi, Ruokolahti;'Vahermanjärvi';Vahermanjärvi, Nummi-Pusula;'Johdasjärvi';Johdasjärvi, Jaala;'Saarijärvi';Saarijärvi, Saarijärvi;'Muntsurinjärvi';Muntsurinjärvi, Lieksa;'Kypärijärvi';Kypärijärvi, Vihti;'Vehmasjärvi';Vehmasjärvi, Suonenjoki;'Ristijärvi';Ristijärvi, Valtimo;'Mäkijärvi';Mäkijärvi, Kuopio;'Keskijärvi';Keskijärvi, Luhanka;'Piilampi';Piilampi, Ruokolahti;'Tammalammi';Tammalammi, Tohmajärvi;'Tuomaslampi';Tuomaslampi, Savitaipale;'Postilampi';Postilampi, Nilsiä;'Päijänne';Päijänne, Asikkala;'Puruvesi';Puruvesi, Kerimäki;'Kemijärvi';Kemijärvi, Kemijärvi;'Näsijärvi';Näsijärvi, Ylöjärvi;'Pyhäjärvi';Pyhäjärvi, Kitee;'Juojärvi';Juojärvi, Tuusniemi;'Konnevesi';Konnevesi, Konnevesi;'Längelmävesi';Längelmävesi, Kangasala;'Vanajavesi';Vanajavesi, Valkeakoski;'Onkivesi';Onkivesi, Lapinlahti;'Pielavesi';Pielavesi, Pielavesi;'Luonteri';Luonteri, Anttola;'Kuonanjärvi';Kuonanjärvi, Kerimäki;'Hossanlampi';Hossanlampi, Hossanlampi;'Pintamojärvi';Pintamojärvi, Pudasjärvi;'Kätkänjärvi';Kätkänjärvi, Lehtimäki;'Ylä-Siika';Ylä-Siika, Nilsiä;'Havanganjärvi';Havanganjärvi, Virrat;|
| |
| | |
| name:name|description:Suuren järven nimi|type:text|
| |
| name:coeff.pH|description:pH:sta johtuva kerroin: ''Ohje: http://intra.eracnet.fi/main/index.php/Minera:Pintavesi''|default:1|
| |
| name:coeff.landuse|description:Maankäytöstä johtuva kerroin: ''Ohje: http://intra.eracnet.fi/main/index.php/Minera:Pintavesi''|default:1
| |
| ">
| |
| library(OpasnetBaseUtils)
| |
| n <- 1
| |
| data <- tidy(op_baseGetData("opasnet_base", "Op_en5462"))
| |
| temp <- data[data$Observation == "Järvi", ]
| |
| temp <- temp[temp$Result == name, "ID"]
| |
| data <- data[data$ID == temp & data$Observation %in% metalli, ]
| |
| conc <- make.ovariable(data)
| |
| conc <- (conc * coeff.pH + conc * coeff.landuse)/2
| |
| conc@sample
| |
| colnames(conc@sample) <- gsub("Observation", "Metalli", colnames(conc@sample))
| |
| colnames(conc@sample) <- gsub("Result", "Pitoisuus", colnames(conc@sample))
| |
| library(xtable)
| |
| print(xtable(conc@sample), type = 'html')
| |
| | |
| </rcode>
| |
| | |
| ==Perustelut==
| |
| | |
| == Katso myös ==
| |
| | |
| Perustuu GTK:n dataan. [http://en.gtk.fi GTK].
| |
| | |
| {{Minera}}
| |
| | |
| == Viitteet ==
| |
| | |
| Lähdeluettelo tulostuu tänne automaattisesti kun viittauksissa käytetään ''<nowiki><ref></ref></nowiki>'' -tagia.
| |
| Lähteet voi myös lisätä halutessaan manuaalisesti.
| |
| | |
| <references/>
| |