GGPLOT2     HEALTH     R     EMISSIONS     STATISTICS     GIS   WRF Open Air Air Pollution  

miércoles, 29 de abril de 2015

Error If you REALLY want Grib2 output from WRF, modify the arch/Config_new.pl - WRF 3.7 in Ubuntu (Elementary)

Well, i was installing WRF Chem in a server and also in my local machine, following this great tutorial, and when i did ./configure, i had the following error. I had all the libraries installed and environmental variables declared.

If you REALLY want Grib2 output from WRF, modify the arch/Config_new.pl script.
Right now you are not getting the Jasper lib, from the environment, compiled into WRF.


I searched and searched and lots of dudes had the same problem. The solution was ridiculous, was change  one line in VRFV3/arch/Config_new.pl, the line 223, from
 $I_really_want_to_output_grib2_from_WRF = "FALSE" ;

to

 $I_really_want_to_output_grib2_from_WRF = "TRUE" ;


and ./configure in WRFV3 and, voilá

The funny was, when i was searching for the answer, lots of people had the same problem. I hope it helps



jueves, 16 de abril de 2015

instalar netcdf, rgdal y rgeos en R con Elementary (Ubuntu 14)

Creé este post por que me dio mas problemas de lo que me imaginaría

hay que instalar las librerías netcdf
Luego instalar los paquetes en R y listo

martes, 10 de marzo de 2015

Open Street Map data em Qgis

Bom, eu trabalho com +OpenStreetMap  e +QGIS  faz tempo, e eu sempre baixaba dados de Geofabrik. Mas hj aprendi como fazer o proceso direitamente em qgis..... e é tão fácil!!!! Nossa, mais rápido e melhor que baixar um shapefile...
Preciso dados OSM de varias cidades de America do Sul, então,
primeiro seleciono o complemento Web -> Open Layer -> Open Street Map.

Fazer sum na cidade, Lima neste caso

ir pra Vetor -> Open Street Map e baixar dados do OpenStreetMap com Extensão a partir do mapa na tela.  O nome de arquivo de saida tem que ter a estensão .osm, lima.osm neste caso.

Agora, ir para Vetor -> Open Street Map -> importar Topologia de XML e procurar o arquivo lima.osm

ir para Vetor -> Open Street Map -> Exportar Topologia para Spatialite e procurar o arquivo lima.osm.db, marcar arquivos de polylinhas (neste caso estou procurando as ruas), e clicar Carregar do DB, e logo selecionar os campos de cada linha.
PRONTO!



jueves, 18 de diciembre de 2014

#QGIS Criar Mapa Simple com QGIS

1._ Baixar shapefiles de http://www.mapcruzin.com/free-mozambique-maps.htm ou http://www.mapmakerdata.co.uk.s3-website-eu-west-1.amazonaws.com/library/stacks/Africa/Mozambique/index.htm ou http://www.diva-gis.org/gdata

2._ Abrir os shapefiles com QGIS


3._ Selecionar o SRC,  neste caso é Lat Lon WGS 84, EPSG: 4326 (mas em outro caso pode ser coordenadas projetadas em UTM 36 S)

4._ Baixar complemento de Open Layers


5._ Selecionar Open Street Map / OCM LandScape e logo fazer sum na camada de MOZ_boundaries com botom dereito acima da camada

6._ Selecionar escala 1:9.000.000 



7._ Clicar um novo compositor de impressão e escrever nome de mapa


8._ Clicar em OK, isto abre nova janela de compositor de impressão

9._ Selecionar em Adicionar mapa

10._ Ajustar o mapa




11._ Adicionar o título

12._ Adicionar escala


13._ Adicionar seta

14._ Exportar o mapa  jpg, png ou otro formato

15._ Olhar mapa

17._ Melhorar mapa e pesquisar perguntando em  Gis.StackExchange 

martes, 25 de noviembre de 2014

#ggplot2 e #RColorBrewer. Grafico com escala invertida


Para uma apresentação tinha que mostrar uma imagem de fatores de emissão e queria mostrar uma escala com os fatores mais poluentes a os menos com ggplot2 e para a escala usei RColorBrewer.


#ggplot2 #R Simple histograma com media de grupos


Soo uma forma simple de fazer um histograma agrupado em ggplot2 para cegos de cores (Cookbook).

martes, 25 de febrero de 2014

R, ggplot2: graficar Intervalos de confianza

Es normal querer graficar intervalos de confianza en las ciencias, y siempre se anda preguntando como hacerlo, sobre todo en R. Pues bien, acá va un tutorial de como hacerlo. Yo lo tengo a partir de las mismas regresiones, pero sa supone puede hacer hacer de cualquier modo. Lo que se necesita es:

Límite inferior del intervalo (Yo lo nombre Yi)
yi
Límite superior del intervalo (Yo lo nombre Yf que viene de final)
yf
y punto central
y
Luego necesitamos el valor de X
x
Por último, necesitamos que todas estas variables estén en una dataframe con columnas y, yi, yf, x , que indican los nombres de las variables de los intervalos de confianza de alguna regresión o datos de error.

El dataframe se vería algo así:




Luego se puede graficar los intervalos de confianza con ggplot2:

Y el gráfico asi



martes, 5 de noviembre de 2013

Ubuntu: Como cambiar de Nautilus a Thunar

Me demoré mas de lo habitual, por eso comparto este post de Ubuntu Geek para olvidarse del lento de Nautilus y pasar a algo mas ligero, como thunar.

Primero, instala thunar... para ello vaya simplemente al Ubuntu Sofware Center y escriba thunar


Después viene la parte de dejar thunar por defecto. Busqué en varias partes, pero lo más efectivo fue:


Después busca


Quedando en mi caso asi


Después logoff, login, se demora un poco y ya está

jueves, 10 de octubre de 2013

#GRASS: Lo básico, importar shapefiles 200 Mb y hacer intersección, Visualizar con Qgis

DATO: Al empezar grass pide definir directorio, que yo llamo SIG y localización, que es una subcarpeta con los datos grass de ese proyecto. Es mejor llamar a esa subcarpeta con iniciales como SP y agregar la proyección, como WGS84LatLon, quedando SPWGS84LatLon .... para evitar posibles errores tontos que hacen perder mucho tiempo.
Mis primeros comandos con grass


Tenía que cargar un shapefile de 200 Mb en qgis y tardaba demasiado. Entonces pregunté en Gis.Stackexchange y me recomendaron usar grass. :)

Las dos primeras líneas del código importan un vector que están en shapefile y lo guardan en formato de grass como roads. Lo mismo con LandUse.
Despues la intersección, que recoge los mapas de grass, y con el operador 'and' realiza la intersección.

Aun no puedo mostrar el mapa en grass, no se por que, pero sí lo hice cargando los mapas de grass en QGIS