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


martes, 17 de septiembre de 2013

GIS.StackExchange: Identificar el mayor polígono dentro de cada celda de una grilla

POSTED ON GIS.StackExchange:
I have a 9x9 grid and a polygon of municipalities of an extensive area of Brazil. My goal is to have a field indicating the biggest municipality per grid cell.
I've been trying with intercept without good results. I know there are similar questions but i cant figure it out how to handle it in QGIS.
Here is a screenshot of what i'm talking about:
ANSWER by Micha:


Resulted in:
Which is the id of the municipality with the largest intersection within each grid cell. Continue...

I've also been trying with Kosmo, according to Juán José del Toro

1.- Create the grid 10x10 km and concatenate Row and Column into a key field

I thid the concatenation with the "calculadora option"



And picture

martes, 10 de septiembre de 2013

Criar clusters em R

Venho fazendo análise de dados con clusters em R e estes são os meus códigos y dataset.

Aqui estou identificando o número de clusters

8

Então códigos gráficos jerárquiso são:


E o gráfico é, e eu posso ver, especificamente, como meus dados são classificados

miércoles, 4 de septiembre de 2013

Seleccionar atributos de un campo y agregarlos a otro campo QGIS Field Calculator

En un vector Polyline en QGIS tengo un campo que se llama TYPE y 'secondary' y 'secondary_link' y otro campo con la longitud de cada arco, entonces me gustaría crear un campo llamado SECONDARY con la longitud LONG de cada arco del campo TYPE para los casos  'secondary' y 'secondary_link'. En QGIS el código es

y la figura es
POSTERIORMENTE  SE RECOMIENDA BORRAR LOS CAMPOS QUE MANTENGAN VALORES NULL, O REEMPLAZAR NULL POR 0

martes, 3 de septiembre de 2013

Creación de grillas de polígono con QGIS / Polygon Grid creation with QGIS

Me pidieron un trabajo en la USP que es generar grillas que sirvan para despues hacer inventarios de emisiones y modelaciones. Estuve buscando y lo hice en QGIS

Me dieron dos puntos, el izquierdo inferior y el derecho superior para crear una grilla con 129 celdas horizontales y 69 verticales. Datos

                        inicial       final
west_east    : -45.02942  -41.26178  -> 129 pontos
south_north  : -23.57699  -21.69774  -> 69  pontos

Luego, en esta página los convertí a UTM WGS 84 23 S acá, y además creé un punto con estas coordenadas geográficas en QGIS y el pluggin MMQGIS, cuargando las coordenadas desde Libre Office a Excel guardadas como .csv:




Despues con MMQGIS importe el punto


Y creé el vector de punto con las coordenadas dadas.

Posteriormente guardé el archivo con el nombre de pUTM.shp con la proyección UTM WGS84 23S, y con este nuevo punto creé el punto para crear la grilla con MMQGIS. Para crear las coordenadas X e Y en la tabla de atributos edite la tabla y actualicé el campo x e y como $x y $y, y automáticamente escribe las coordenadas:



Una vez con el punto en coordenadas UTM, creé la grilla con MMQGIS:


Con las coordenadas x e y. H  V son las dimensiones de cada celda. Width y Height es la altura de la grilla> 129*1000 y 69 * 1000

Y así obtuve la grilla: