#Rstats #ggplot2 #geom_raster() #Calpuff Plotear salidas de Calpuff en R

He estado trabajando para plotear las salidas de concentraciones 24h promedio de Calpuff en R y después de harto trabajo he llegado a muy pocas líneas de programación:

Cargaré el archivo de salida de Calpuff, promedio de 24 h:

Como se ve, tiene un header de 4 líneas innecesarias hasta llegar a la información que tiene el nombre de las variables, por lo que es necesario cargar el archivo, eliminando esas 4 líneas, automáticamente se eliminará la línea 6 que no tiene información, y despues es necesario llenar los vacíos de los nombres de las variables:

	setwd("c:/tesis/emisiones"); #define tu ruta Set Working Directory
	#Eliminar header innecesario
	#skip elimina las primeras 4 filas
	#fill llena con NA en las ?ltimas columnas
	calpuff<-read.table("RANK(ALL)_PM10_24HR_CONC_C.DAT",header=T, skip=4, fill=T)
	attach(calpuff)
	names(calpuff)
	head(calpuff)
	#Se aprecia como los nombres de las variables estan desplazados
	#Es necesario renombrar las variables
	#install.packages("reshape"); #solo si es necesario
	library(reshape)
	names(calpuff) <- c("x","y","mp10", "lat","long")
	names(calpuff)
	attach(calpuff)
	#aca se aprecia como se renombraron las variables
	summary(calpuff)
	#Ahora, a plotear las concentraciones
	library(ggplot2)
	ggplot(calpuff, aes(x,y, fill = mp10))+
	geom_raster()+coord_equal()+
	ggtitle("Concentraciones de MP10")+
	theme(plot.title = element_text(lineheight=.8, face="bold"),
	axis.text.x=element_text(size=12, face="bold"),
	axis.text.y=element_text(size=12, face="bold"),
	legend.position=c(1,0.82),
	legend.justification=c(1,1))+
	ylab("Coordenada Y UTM") +
	xlab("Coordenada X UTM")+
	scale_fill_gradient("Count", low = "white", high = "black")

view raw gistfile1.r hosted with ❤ by GitHub