File:Maya area ptopet july 1990 1.svg
Original file (SVG file, nominally 630 × 630 pixels, file size: 1.3 MB)
Captions
Captions
Summary[edit]
| Description |
English: Maya area P/PET of July 1990. This means aridity of summer. Peten area neat Tikal is anccording this quite dry area. |
| Date | |
| Source | Own work |
| Author | Merikanto |
This image is based on Chelsa bioclim+ dataset.
Brun, P., Zimmermann, N.E., Hari, C., Pellissier, L., Karger, D. (2022): Data from: CHELSA-BIOCLIM+ A novel set of global climate-related predictors at kilometre-resolution. EnviDat. https://doi.org/10.16904/envidat.332
Karger D.N., Conrad, O., Böhner, J., Kawohl, T., Kreft, H., Soria-Auza, R.W., Zimmermann, N.E, Linder, H.P., Kessler, M. (2018): Data from: Climatologies at high resolution for the earth’s land surface areas. EnviDat. https://doi.org/10.16904/envidat.228.v2.1
R code to download data and visualize it, valid at January 2024:
library(raster)
library(ggplot2)
library(dplyr)
library(RColorBrewer)
- download jult pet, precip from CHELSA
download_from_chelsa=0
if(download_from_chelsa==1)
{
download.file(url = 'https://os.zhdk.cloud.switch.ch/envicloud/chelsa/chelsa_V2/GLOBAL/monthly/pet/CHELSA_pet_penman_07_1990_V.2.1.tif',
destfile = 'chelsa_1990_07_pet.tif')
download.file(url = 'https://os.zhdk.cloud.switch.ch/envicloud/chelsa/chelsa_V2/GLOBAL/monthly/pet/CHELSA_pet_penman_03_1990_V.2.1.tif',
destfile = 'chelsa_1990_03_pet.tif')
download.file(url = 'https://os.zhdk.cloud.switch.ch/envicloud/chelsa/chelsa_V2/GLOBAL/monthly/pr/CHELSA_pr_07_1990_V.2.1.tif',
destfile = 'chelsa_1990_07_ppt.tif')
download.file(url = 'https://os.zhdk.cloud.switch.ch/envicloud/chelsa/chelsa_V2/GLOBAL/monthly/pr/CHELSA_pr_03_1990_V.2.1.tif',
destfile = 'chelsa_1990_03_ppt.tif')
download.file(url = 'https://os.zhdk.cloud.switch.ch/envicloud/chelsa/chelsa_V1/chelsa_trace/orog/CHELSA_TraCE21k_dem_0_V1.0.tif',
destfile = 'chelsa_dem_00.tif')
}
pet07<-raster('chelsa_1990_07_pet.tif')
prec07<-raster('chelsa_1990_07_ppt.tif')
pet03<-raster('chelsa_1990_03_pet.tif')
prec03<-raster('chelsa_1990_03_ppt.tif')
- pet0<-raster('chelsa_2010_07_pet.tif')
- prec0<-raster('chelsa_2010_07_ppt.tif')
- pet03<-raster('chelsa_2010_03_pet.tif')
- prec03<-raster('chelsa_2010_03_ppt.tif')
dem0<-raster('chelsa_dem_00.tif')
ext <- extent(c(xmin = -94, xmax = -86, ymin = 14, ymax = 22))
prec <- crop(x = prec07, y = ext)
pet<- crop(x = pet07,y = ext)
dem <- crop(x = dem0, y = ext)
- plot(prec)
ptopet<-(prec*100)/pet
ptopet[ptopet>200]<-200
- plot(ptopet)
pal <- colorRampPalette(c("brown","yellow","green","blue"))
plot(ptopet,
col=pal(256),
main="Maya area July 1990 P/PET percent",
xlab="longitude",
ylab="latitude")
#sub="P/PET percent July 2010")
contour(ptopet, add=T, alpha=50)
points(c(-89.622778 ),c(17.221944 ), col = "red",pch=2, cex=2, lwd=3, add=T)
text(c(-89.622778 ),c(17.221944 ), labels=c("Tikal"), cex= 2, pos=3)
dev.off()
system("pdf2svg Rplots.pdf maya_area_ptopet_july_1990_1.svg")
writeRaster(ptopet , filename="ptopet.nc", bandorder='BSQ',
format="NetCDF", overwrite=TRUE)
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 | 11:10, 7 January 2024 | | 630 × 630 (1.3 MB) | Merikanto (talk | contribs) | Uploaded own work with UploadWizard |
You cannot overwrite this file.
File usage on Commons
There are no pages that use this file.
