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