Järvisedimenttien metallipitoisuudet

Versio hetkellä 3. kesäkuuta 2012 kello 12.03 – tehnyt Jouni (keskustelu | muokkaukset) (luokiteltu)

(ero) ← Vanhempi versio | Nykyinen versio (ero) | Uudempi versio → (ero)

Tämä sivu on tietokide alatyypiltään muuttuja. Sivutunniste: Op_fi2953
Moderaattori:Marjo (katso kaikki) Sivun edistymistä ei ole arvioitu. Arvostuksen määrää ei ole arvioitu (ks. peer review).
Lisää dataa {{#opasnet_base_link:Op_fi2953}}

Tämä sivu on tietokide alatyypiltään muuttuja.
Sivutunniste: Op_fi2953


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".

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')

Järvisedimenttien metallipitoisuudet

Sisällys

Kysymys

Vastaus

1) Lasketaan järven sedimentin metallipitoisuus perustuen ympäröivän, ympyrän säteen määräämän alueen järvien keskipitoisuuteen.

2) Lasketaan järven sedimentin metallipitoisuus läheisen/sopivan suuren järven pitoisuuden perusteella.

Katso myös

Viitteet

Navigointivalikko