File:Mis3 france warm kriged tjuly 1.svg
Original file (SVG file, nominally 1,650 × 1,275 pixels, file size: 1.02 MB)
Captions
Summary
[edit]DescriptionMis3 france warm kriged tjuly 1.svg |
English: Temperature of July in France during warm instage on MIS3. |
|||
Date | ||||
Source | Own work | |||
Author | Merikanto | |||
SVG development InfoField |
|
Tjuly data is from Davies& Gollop 2003
Davies, William & Gollop, Piers. (2003). The Human Presence in Europe during the Last Glacial Period II: Climate Tolerance and Climate Preferences of Mid-and Late Glacial Hominids.
Original Tjuly data is from Pollard et al OIS3 project.
Barron, E., & Pollard, D. (2002). High-resolution climate simulations of oxygen isotope stage 3 in Europe. Quaternary Research, 58(3), 296-309. https://doi.org/10.1006/qres.2002.2374
Some points data is intuitively manually selected from plottings.
Kriging with etopo1 dataset and r. Visualization with NASA Panoply. Conversion to SVG on net.
"ID","lon","lat","TS","Elev"
"Calais",1.85,50.95,10.0,0
"Marse",5.0,41.0,25.0,0.0
"Marse2",4.0,43.0,23.0,0.0
"P5",0.0,50.0,10.0,0.0
"Up1",15.0,47.5,12,1000.0
"Up2",15.0,46.5,8,2000.0
"Mare1",-10.0,46.5,10,00.0
"Aqvit1",0.5,44.5,17,200
"LesEyzies",1.018,44.936,17.0,75
"Loire",0,47,17.0,250
library(raster)
require(rgdal)
library(sp)
library(gstat)
library(automap)
library(RColorBrewer)
library(viridis)
elevname0<-"./eurodata/area_glacial.nc"
outname1="base_tk7.nc"
elev1<-raster(elevname0)
names(elev1)<-"Elev"
print(names(elev1))
print(elev1)
pointable1 <- read.csv("./mis3warm4.txt")
projstr1<-"+proj=longlat +datum=WGS84 +no_defs +ellps=WGS84 +towgs84=0,0,0"
print (projstr1)
points1 <- SpatialPointsDataFrame(coords = pointable1[,c("lon","lat")],
data = pointable1,
proj4string = CRS(projstr1))
names(points1)<-c("ID","lon","lat","TS","Elev")
knitr::kable(head(points1, n=10))
print("SP")
points1$extracted <- extract(x=elev1, y=points1)
newgrid <- as(elev1, "SpatialGridDataFrame")
names(newgrid) <- "Elev"
print("Variogram ...")
tsvar1 <- variogram(TS~Elev, points1)
print(".")
sph.model <- vgm(psill=15, model="Exp", range=500, nugget=20)
print(".")
sph.fit <- fit.variogram(object = tsvar1, model = sph.model)
print("Kriging ...")
rk1 <- krige(TS~Elev, points1, newgrid, model = sph.fit)
rout0 <- as(rk1[1], "RasterLayer")
plot(rout0)
plot(rout0, col=viridis(100))
crs(rout0) <- "+proj=longlat +datum=WGS84 +no_defs +ellps=WGS84 +towgs84=0,0,0"
writeRaster(rout0, filename=outname1, overwrite=TRUE, format="CDF", varname="T_July", varunit="Celsius",
longname="Temperature of July", xname="lon", yname="lat")
Licensing
[edit]- You are free:
- to share – to copy, distribute and transmit the work
- to remix – to adapt the work
- Under the following conditions:
- attribution – You must give appropriate credit, provide a link to the license, and indicate if changes were made. You may do so in any reasonable manner, but not in any way that suggests the licensor endorses you or your use.
- share alike – If you remix, transform, or build upon the material, you must distribute your contributions under the same or compatible license as the original.
File history
Click on a date/time to view the file as it appeared at that time.
Date/Time | Thumbnail | Dimensions | User | Comment | |
---|---|---|---|---|---|
current | 08:52, 30 September 2019 | 1,650 × 1,275 (1.02 MB) | Merikanto (talk | contribs) | User created page with UploadWizard |
You cannot overwrite this file.
File usage on Commons
There are no pages that use this file.
File usage on other wikis
The following other wikis use this file:
- Usage on uk.wikipedia.org
Metadata
This file contains additional information such as Exif metadata which may have been added by the digital camera, scanner, or software program used to create or digitize it. If the file has been modified from its original state, some details such as the timestamp may not fully reflect those of the original file. The timestamp is only as accurate as the clock in the camera, and it may be completely wrong.
Width | 1650px |
---|---|
Height | 1275px |