{"id":630,"date":"2018-05-03T16:10:04","date_gmt":"2018-05-03T16:10:04","guid":{"rendered":"http:\/\/ipn.paleofire.org\/?p=630"},"modified":"2018-08-11T21:10:36","modified_gmt":"2018-08-11T21:10:36","slug":"gcd-4-0-2-and-paleofire-1-2-2","status":"publish","type":"post","link":"https:\/\/ipn.paleofire.org\/?p=630","title":{"rendered":"GCD 4.0.2 and paleofire 1.2.2"},"content":{"rendered":"

\"\"<\/a><\/p>\n

We are please to announce that the GCD and paleofire R packages have been updated to their 4.0.2 and 1.2.2 versions, respectively. GCD major update from 3.X.X to 4.X.X version number is associated to extensive changes in the way data is added into the package. GCD 4.X.X is now mirroring the online GCD SQL database on a monthly basis or whenever a significant number of new charcoal sites are added to the database.<\/p>\n

Three new fields in the paleofiresites<\/em> table (see ?paleofiresites<\/em> for details) enable to perform analyses with specific database versions or enable to select sites according to the date at which they have been added into the database. See the quick example below for explanations:<\/p>\n

[code language=”R”]
\npar(mfrow=c(2,2))
\nplot( pfSiteSel(num_version &amp;lt; 400 ) ) # All sites in GCD version before 4.0.0
\nplot( pfSiteSel(gcd_version == "GCD1" ) ) # All sites in GCDv1
\nplot( pfSiteSel(update_date &amp;lt; "2016-01-01" ) ) # Sites included before 2016-01-01
\nplot( pfSiteSel(update_date &amp;gt; "2018-01-01" ) ) # Sites included since 2018-01-01&lt;\/code&gt;
\n[\/code]<\/p>\n

\"\"<\/p>\n

The map above that is used to display all GCD charcoal sites has been produced by adapting the code available at\u00a0https:\/\/seethedatablog.wordpress.com\/2016\/12\/23\/r-simple-world-map-robinson-ggplot\/<\/a>\u00a0please follow that link for additional references and explanations. The modified code is copied below:<\/p>\n

[code language=”R”]<\/pre>
\n# ======================================================================================
\n# Create a simple world map in Robinson projection with labeled graticules using ggplot
\n# ======================================================================================<\/p>\n

# Set a working directory with setwd() or work with an RStudio project<\/p>\n

# __________ Set libraries
\nlibrary(rgdal) # for spTransform() & project()
\nlibrary(ggplot2) # for ggplot()
\nlibrary(paleofire) #<\/p>\n

setwd("\/Users\/OlivierB\/Documents\/R_in_progress\/MAP_GCD\/")<\/p>\n

# __________ Load ready to use data from GitHub
\nload(url("https:\/\/github.com\/valentinitnelav\/RandomScripts\/blob\/master\/NaturalEarth.RData?raw=true"))
\n# This will load 6 objects:
\n# xbl.X & lbl.Y are two data.frames that contain labels for graticule lines
\n# They can be created with the code at this link:
\n# https:\/\/gist.github.com\/valentinitnelav\/8992f09b4c7e206d39d00e813d2bddb1
\n# NE_box is a SpatialPolygonsDataFrame object and represents a bounding box for Earth
\n# NE_countries is a SpatialPolygonsDataFrame object representing countries
\n# NE_graticules is a SpatialLinesDataFrame object that represents 10 dg latitude lines and 20 dg longitude lines
\n# (for creating graticules check also the graticule package or gridlines fun. from sp package)
\n# (or check this gist: https:\/\/gist.github.com\/valentinitnelav\/a7871128d58097e9d227f7a04e00134f)
\n# NE_places – SpatialPointsDataFrame with city and town points
\n# NOTE: data downloaded from http:\/\/www.naturalearthdata.com\/
\n# here is a sample script how to download, unzip and read such shapefiles:
\n# https:\/\/gist.github.com\/valentinitnelav\/a415f3fbfd90f72ea06b5411fb16df16<\/p>\n

# __________ Project from long-lat (unprojected data) to Robinson projection
\n# spTransform() is used for shapefiles and project() in the case of data frames
\n# for more PROJ.4 strings check the followings
\n# http:\/\/proj4.org\/projections\/index.html
\n# https:\/\/epsg.io\/<\/p>\n

PROJ <- "+proj=robin +lon_0=0 +x_0=0 +y_0=0 +ellps=WGS84 +datum=WGS84 +units=m +no_defs"
\n# or use the short form "+proj=robin"
\nNE_countries_rob <- spTransform(NE_countries, CRSobj = PROJ)
\nNE_graticules_rob <- spTransform(NE_graticules, CRSobj = PROJ)
\nNE_box_rob <- spTransform(NE_box, CRSobj = PROJ)<\/p>\n

# Add lakes http:\/\/www.naturalearthdata.com\/http\/\/www.naturalearthdata.com\/download\/10m\/physical\/ne_10m_lakes.zip<\/p>\n

download.file("http:\/\/www.naturalearthdata.com\/http\/\/www.naturalearthdata.com\/download\/10m\/physical\/ne_10m_lakes.zip", destfile="ne_10m_lakes.zip")
\nunzip("ne_10m_lakes.zip")
\nNE_lakes <- readOGR(‘ne_10m_lakes.shp’,
\n‘ne_10m_lakes’)
\nNE_lakes <- spTransform(NE_lakes, CRSobj = PROJ)
\nNE_lakes <- NE_lakes[NE_lakes@data$scalerank==0,]<\/p>\n

# project long-lat coordinates for graticule label data frames
\n# (two extra columns with projected XY are created)
\nprj.coord <- project(cbind(lbl.Y$lon, lbl.Y$lat), proj=PROJ)
\nlbl.Y.prj <- cbind(prj.coord, lbl.Y)
\nnames(lbl.Y.prj)[1:2] <- c("X.prj","Y.prj")<\/p>\n

prj.coord <- project(cbind(lbl.X$lon, lbl.X$lat), proj=PROJ)
\nlbl.X.prj <- cbind(prj.coord, lbl.X)
\nnames(lbl.X.prj)[1:2] <- c("X.prj","Y.prj")
\ndata("paleofiresites")
\n# Project GCD sites<\/p>\n

gcd_rob=project(cbind(paleofiresites$long,paleofiresites$lat),PROJ)
\ngcd=data.frame(long=gcd_rob[,1],lat=gcd_rob[,2],GCD_version=paleofiresites$gcd_version)<\/p>\n

# __________ Plot layers
\nggplot() +
\n# add Natural Earth countries projected to Robinson, give black border and fill with gray
\ngeom_polygon(data=NE_box_rob, aes(x=long, y=lat), colour="black", fill="grey90", size = 0.25) +
\ngeom_polygon(data=NE_countries_rob, aes(long,lat, group=group), colour="black", fill="white", size = 0.25) +
\ngeom_polygon(data=NE_lakes, aes(long,lat, group=group), colour="black", fill="white", size = 0.25) +
\n# Note: "Regions defined for each Polygons" warning has to do with fortify transformation. Might get deprecated in future!
\n# alternatively, use use map_data(NE_countries) to transform to data frame and then use project() to change to desired projection.
\n# add Natural Earth box projected to Robinson
\ngeom_polygon(data=NE_box_rob, aes(x=long, y=lat), colour="black", fill="transparent", size = 0.25) +
\n# add graticules projected to Robinson
\ngeom_path(data=NE_graticules_rob, aes(long, lat, group=group), linetype="dotted", color="grey50", size = 0.25) +
\n# add graticule labels – latitude and longitude
\ngeom_text(data = lbl.Y.prj, aes(x = X.prj, y = Y.prj, label = lbl), color="grey50", size=2) +
\ngeom_text(data = lbl.X.prj, aes(x = X.prj, y = Y.prj, label = lbl), color="grey50", size=2) +
\n# the default, ratio = 1 in coord_fixed ensures that one unit on the x-axis is the same length as one unit on the y-axis
\ngeom_point(data=gcd, aes(long,lat, col= GCD_version),pch=2)+
\ncoord_fixed(ratio = 1) +
\n# remove the background and default gridlines
\ntheme_void()+
\ntheme(legend.title = element_text(colour="black", size=8, face="bold"), # adjust legend title
\nlegend.position = c(0.1, 0.2), # relative position of legend
\nplot.margin = unit(c(t=0, r=0, b=0, l=0), unit="cm"),
\nlegend.background = element_rect(fill="white",
\nsize=0.5, linetype="solid",
\ncolour ="black")) # adjust margins<\/p>\n

# save to pdf and png file
\nggsave("map_draft_1.pdf", width=28, height=13.5, units="cm")
\nggsave("map_1.png", width=26, height=13.5, units="cm", dpi=300)<\/p>\n

# REFERENCES:
\n# This link was useful for graticule idea
\n# http:\/\/stackoverflow.com\/questions\/38532070\/how-to-add-lines-of-longitude-and-latitude-on-a-map-using-ggplot2
\n# Working with shapefiles, projections and world maps in ggplot
\n# http:\/\/rpsychologist.com\/working-with-shapefiles-projections-and-world-maps-in-ggplot
\n<pre>[\/code]<\/p>\n","protected":false},"excerpt":{"rendered":"

We are please to announce that the GCD and paleofire R packages have been updated to their 4.0.2 and 1.2.2 versions, respectively. GCD major update from 3.X.X to 4.X.X version…<\/p>\n","protected":false},"author":1,"featured_media":632,"comment_status":"closed","ping_status":"open","sticky":false,"template":"","format":"standard","meta":{"footnotes":""},"categories":[3],"tags":[],"class_list":["post-630","post","type-post","status-publish","format-standard","has-post-thumbnail","hentry","category-analyses"],"_links":{"self":[{"href":"https:\/\/ipn.paleofire.org\/index.php?rest_route=\/wp\/v2\/posts\/630","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/ipn.paleofire.org\/index.php?rest_route=\/wp\/v2\/posts"}],"about":[{"href":"https:\/\/ipn.paleofire.org\/index.php?rest_route=\/wp\/v2\/types\/post"}],"author":[{"embeddable":true,"href":"https:\/\/ipn.paleofire.org\/index.php?rest_route=\/wp\/v2\/users\/1"}],"replies":[{"embeddable":true,"href":"https:\/\/ipn.paleofire.org\/index.php?rest_route=%2Fwp%2Fv2%2Fcomments&post=630"}],"version-history":[{"count":10,"href":"https:\/\/ipn.paleofire.org\/index.php?rest_route=\/wp\/v2\/posts\/630\/revisions"}],"predecessor-version":[{"id":677,"href":"https:\/\/ipn.paleofire.org\/index.php?rest_route=\/wp\/v2\/posts\/630\/revisions\/677"}],"wp:featuredmedia":[{"embeddable":true,"href":"https:\/\/ipn.paleofire.org\/index.php?rest_route=\/wp\/v2\/media\/632"}],"wp:attachment":[{"href":"https:\/\/ipn.paleofire.org\/index.php?rest_route=%2Fwp%2Fv2%2Fmedia&parent=630"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/ipn.paleofire.org\/index.php?rest_route=%2Fwp%2Fv2%2Fcategories&post=630"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/ipn.paleofire.org\/index.php?rest_route=%2Fwp%2Fv2%2Ftags&post=630"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}