source("~/SparkleShare/Rprofile/R/Rsettings/phdRsettings.R")
### Make study Boundary
jp1.SPDF <- getData('GADM', country='JPN', level=1, path = "~/Dropbox/2data//dataRaw/gadm2")
basemap.r <- readRDS("~/Dropbox/2data/dataProduct/hkd/hkd_google_satellite_142.5_43.5_zoom6_141018_2339.Rds")
jpVolA.spdf <- readRDS("~/Dropbox/2data/dataProduct/jpVolcanoes/jpVol110_140812_174525.Rds")
jpVolQ.spdf <- readRDS("~/Dropbox/2data/dataProduct/jpVolcanoes/jpVol455_140812_172148.Rds")
sap.spdf <- readRDS("~/Dropbox/2data/data/Sapporo_140817_162919.Rds")
xmin <- 139
xmax <- 146
ymin <- 41.4
ymax <- 45.8
bbox.SPDF <- ge.xy2bboxSPDF(xmin,xmax,ymin,ymax,wgs84GRS)
bh <- readRDS("~/Dropbox/2data/dataProduct/hkd/hkd_profiles_140806_164333.Rds")
bh_xy <- bh[!duplicated(bh$ID),]
bh_xy$grp <- cut(bh_xy$TD, breaks = c(0,500,1000,1500,2000,2200),
labels = c("0-500","500-1000","1000-1500","1500-2000","2000-2200"))
proj4string(jpVolQ.spdf) <- proj4string(bbox.SPDF)
volQ <- jpVolQ.spdf[bbox.SPDF,]
volA <- jpVolA.spdf[bbox.SPDF,]
hkdLand <- ge.LargestPolys(jp1.SPDF, Polygon=T)
plot(hkdLand)
volQ2 <- jpVolQ.spdf[hkdLand,]
limitsX <- c(138,147)
breaksX <- seq(limitsX[1], limitsX[2],1)
labelsX=parse(text=paste(breaksX, "^o ", "*E", sep=""))
limitsY <- c(41,47)
breaksY <- seq(limitsY[1],limitsY[2],1)
labelsY=parse(text=paste(breaksY, "^o ", "*N", sep=""))
## Layer0: Base map
ggBH <- ggmap(basemap.r, extent = "panel") +
### Layers
geom_point(data = bh_xy, aes(Lon, Lat,fill = grp, size = grp),
shape = 21, alpha = 0.9) +
### X
xlab("Lontitude") +
scale_x_continuous(breaks=breaksX,
labels=labelsX,
limits=limitsX,
expand = c(0.01,0.01)) +
theme(axis.text.y = element_text(angle = 90, hjust = 0.5, vjust = 0)) +
### Y
ylab("Latitude") +
scale_y_continuous(breaks=breaksY,
labels=labelsY,
limits=limitsY,
expand = c(0.01,0.01)) +
###Legend
### Size
labs(size = "Borehole Depth (m)") +
scale_size_manual(values=c(1,1.5,2,3,4)) +
### fill
scale_fill_brewer("Borehole Depth (m)", palette="Blues")
ggBH
### Color
ggVol <- ggBH +
geom_point(data = volQ@data,
aes(as.numeric(lon), as.numeric(lat),
color="blue"), shape = 17, alpha = 0.7) +
geom_point(data = volA@data,
aes(as.numeric(lon), as.numeric(lat),
color="red"), shape = 17, size = 3) +
scale_color_manual(name = "Volcanoes", values = c("orange","red"), labels = c("Quaternary Volcanoes","Active Volcanoes")) +
geom_path(data = volQ2@data, aes(as.numeric(lon), as.numeric(lat)),size = 12, alpha = 0.3, colour = "yellow",lineend = "round")
ggVol
ggSap <- ggVol + geom_point(data = sap.spdf, aes(x = lon, y = lat), colour = "White") + geom_text(data = sap.spdf, aes(x = lon, y = lat, label = name), hjust = -0.1,family="Times", face="italic", colour="white")
# ggSap
ggBar <- ggSap +scaleBar(lon = 145.5, lat = 41, distanceLon = 50, distanceLegend = 30,distanceLat = 15, dist.unit = "km", arrow.length = 60, arrow.distance = 500, arrow.North.size = 4,legend.colour = "white", arrow.North.color = "white", arrow.colour = "blue")
ggFont <- ggBar +
#coord_equal() +
theme_bw(base_family = "Times", base_size = 12)
ggFont
f01_hkdStudyArea <- ggFont
##ge.ggsave(hkdStudyArea)
## ggWRS2
library(wrspathrow)
wrs2.SPDF <- pathrow_num(x = hkdLand, as_polys = TRUE)
plot(wrs2.SPDF,col = "red")
wrs2.SPDF@data
wrs2.df <- fortify(wrs2.SPDF)
ggWRS <- ggFont + geom_polygon(aes(long,lat,group=group),
color = "grey", alpha = 0.3, fill = NA,
linetype = 3,
data=wrs2.df) +
geom_point(data = wrs2.SPDF@data, color = "black",
aes(x = centroid_x, y = centroid_y, shape = MODE ))
ggWRS2 <- ggWRS + geom_text(data = wrs2.SPDF@data,
aes(x = centroid_x, y = centroid_y),
label = paste(wrs2.SPDF@data$PATH, wrs2.SPDF@data$ROW," "),
family="Times", face = "Italic", colour="black") +
scale_shape_manual(name = " WRS 2", values = 20 , labels = c("Path and Row"))
ggWRS2
f01_hkdStudyArea <-ggWRS2
ge.ggsave(f01_hkdStudyArea)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.