Rainfall Amount Flooding Quezon City Philippines

The rainfall received by Quezon City, Philippines was almost double of what the city normally receive for the entire month of August, causing flooding and land slide to various villages in the area.

August 6-7 Rainfall on Metro Manila

Majority of Metro Manila is affected by floods. Looking at the hourly data from PAGASA weather  station located at Bicutan, Taguig, the graph below will gave the viewers of this blog on the rainfall situation in Manila from August 6 (12:00am)-7(11:15 pm), 2012.

Provincial Monthly Rainfall of the Philippines from WORLDCLIM

Preparing for a future conference on climate change, I downloaded and extracted average monthly rainfall in the Philippines from worldclim.org. Using maptools, raster, and animation package in R, I produced an animation of average monthly rainfall of the country.

Provincial Map using GADM

This blog demonstrates how to produce political/provincial boundary map (below) using R maptools and raster packages.


 

	## Load required packages
	library(maptools)
	library(raster)
	## Download data from gadm.org
	adm <- getData('GADM', country='PHL', level=2)
	mar<-(adm[adm$NAME_1=="Marinduque",])plot(mar, bg="dodgerblue", axes=T)
	##Plot downloaded data
	plot(mar, lwd=10, border="skyblue", add=T)
	plot(mar,col="green4", add=T)
	grid()
	box()
	invisible(text(getSpPPolygonsLabptSlots(mar), labels=as.character(mar$NAME_2), cex=1.1, col="white", font=2))
	mtext(side=3, line=1, "Provincial Map of Marinduque", cex=2)
	mtext(side=1, "Longitude", line=2.5, cex=1.1)
	mtext(side=2, "Latitude", line=2.5, cex=1.1)
	text(122.08,13.22, "Projection: Geographic\nCoordinate System: WGS 1984\nData Source: GADM.org\nCreated by: ARSsalvacion", adj=c(0,0), cex=0.7, col="grey20")

view raw Marinduque.r hosted with ❤ by GitHub