Файл:Mis3 france warm kriged tjuly 1.svg
Повна роздільність (SVG-файл, номінально 1650 × 1275 пікселів, розмір файлу: 1,02 МБ)
Відомості про цей файл містяться на Вікісховищі — централізованому сховищі вільних файлів мультимедіа для використання у проектах Фонду Вікімедіа. |
Опис файлу
ОписMis3 france warm kriged tjuly 1.svg |
English: Temperature of July in France during warm instage on MIS3. |
|||
Час створення | ||||
Джерело | Власна робота | |||
Автор | Merikanto | |||
SVG розвиток InfoField | Вихідний код цього SVG-файлу правильний.
|
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")
Ліцензування
- Ви можете вільно:
- ділитися – копіювати, поширювати і передавати твір
- модифікувати – переробляти твір
- При дотриманні таких умов:
- зазначення авторства – Ви повинні вказати авторство, надати посилання на ліцензію і вказати, чи якісь зміни було внесено до оригінального твору. Ви можете зробити це в будь-який розсудливий спосіб, але так, щоб він жодним чином не натякав на те, наче ліцензіар підтримує Вас чи Ваш спосіб використання твору.
- поширення на тих же умовах – Якщо ви змінюєте, перетворюєте або створюєте іншу похідну роботу на основі цього твору, ви можете поширювати отриманий у результаті твір тільки на умовах такої ж або сумісної ліцензії.
Об'єкти, показані на цьому файлі
зображує
Якесь значення без елемента на сайті Вікідані
30 вересня 2019
image/svg+xml
Історія файлу
Клацніть на дату/час, щоб переглянути, як тоді виглядав файл.
Дата/час | Мініатюра | Розмір об'єкта | Користувач | Коментар | |
---|---|---|---|---|---|
поточний | 08:52, 30 вересня 2019 | 1650 × 1275 (1,02 МБ) | Merikanto | User created page with UploadWizard |
Використання файлу
Така сторінка використовує цей файл:
Метадані
Файл містить додаткові дані, які зазвичай додаються цифровими камерами чи сканерами. Якщо файл редагувався після створення, то деякі параметри можуть не відповідати цьому зображенню.
Ширина | 1650px |
---|---|
Висота | 1275px |