Järvisedimenttien metallipitoisuudet

Versio hetkellä 4. toukokuuta 2012 kello 07.59 – tehnyt Marjo (keskustelu | muokkaukset) (→‎Vastaus)

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

Tämä sivu on tietokide alatyypiltään muuttuja. Sivutunniste: Op_fi2953
Moderaattori:Ei ole (katso kaikki) Kuinka ryhtyä moderaattoriksi? 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 * coeff.landuse
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 ison järven pitoisuuden perusteella.

Katso myös

Viitteet

Navigointivalikko