R se está convirtiendo en una herramienta bastante fuerte para manejar y analizar datos espaciales. Aprendí algunas cosas útiles a través de preguntas como estas en SO y pensé que podría ser útil tener algo similar, pero más orientado 'espacialmente'.
¿Puedes compartir algunos consejos y trucos espaciales de R que te resultaron útiles?
tips-and-tricks
r
revs radek
fuente
fuente
Respuestas:
Esto no es tanto un truco como lo es
spplot()
la ingeniosa funcionalidad incorporada.spplot()
La capacidad de escalar muestras de leyenda (para que coincida con los rangos de ruptura de clasificación) sirve como una herramienta pedagógica útil cuando se discute la distribución de datos de atributos y los tipos de clasificación. La combinación de parcelas de distribución acumulativa con los mapas ayuda en este esfuerzo.Los estudiantes solo necesitan modificar algunos parámetros de script para explorar los tipos de clasificación y los efectos de transformación de datos. Esta suele ser su primera incursión en R en un curso centrado principalmente en ArcGIS.
Aquí hay un fragmento de código:
Ref: Análisis de datos espaciales aplicados con R (R. Bivand, E Pebesma y V. Gomez-Rubio)
fuente
EDITAR: tenga en cuenta que esto ya no funciona 2018-10-24, debido a los nuevos requisitos para las fuentes de mapas de Google.
Estaba bastante feliz de encontrar el paquete de dismo con geocodificación y descarga de google maps:
Eso está en R 2.12.0 en Windows, es trivial instalar dismo y sus dependencias allí, no estoy seguro en otros sistemas.
fuente
e <- extent(x[4:7] + c(-0.001, 0.001, -0.001, 0.001))
anuncio de línea y recibo un mensaje de errorError: c("x", "y") %in% names(x) is not all TRUE
.x[4:7]
parece bien sin embargo; alguna idea sobre cuál podría ser el problema?x <- geocode('110 George Street, Bathurst, NSW, Australia')
regresaZERO_RESULTS
por ejemplo, y cuando se utiliza un ejemplo que devuelve una latitud / longitud, la funcióne <- extent(x[4:7] + c(-0.001, 0.001, -0.001, 0.001)) also fails.
extent
requiere un vector de números. Entonces esto funcionae <- extent(c(x[,4], x[,5], x[,6], x[,7]) + c(-0.001, 0.001, -0.001, 0.001))
.e <- extent(as.numeric(x[4:7]) + c(-0.001, 0.001, -0.001, 0.001))
Aquí hay algunos que he recogido.
cran.r espacial , blog gmane , análisis de series de tiempo , geodatacenter
fuente
Tampoco es un truco, pero aquí hay algunos recursos / ejemplos que he recopilado
Un ejemplo de trazar pequeños mapas múltiples de datos de Areal en R usando el paquete de celosía.
Hay algunas preguntas sobre StackOverflow sobre mapeo y R, y aquí hay una con un buen ejemplo. Me gustaría ver las otras respuestas y los recursos que ofrecen (así como buscar algunos ejemplos más) en SO también.
Un enlace diferente al mismo grupo r-sig-geo que Brad ya dio. Es muy activo, y Roger Bivand responde preguntas prácticamente todos los días en el grupo. Ambos relacionados con la programación y el análisis estadístico.
Además de revisar la página espacial del cráneo , también sugeriría que revise específicamente la página Spatstat mantenida por Adrian Baddeley. Muchos ejemplos, un curso y un próximo libro electrónico. (En este momento he estado siguiendo el curso de statstat , y creo que es una introducción mucho más amable que el libro Bivand).
No es un recurso gratuito, pero para cualquier persona interesada en RI, sugeriría que eche un vistazo a Use R! serie de Springer. Tienen un libro Análisis de datos espaciales aplicados con R directamente pertinente (también el libro A Beginner's Guide to R es mi libro de aprendizaje sugerido de R).
Un libro electrónico gratuito, Una guía práctica para el mapeo geoestadístico (Hengl 2009), tiene ejemplos de geoestatos aplicados en R, GRASS y Google Earth (KML).
Si encuentro más ejemplos buenos, continuaré actualizando (¡espero que otras personas publiquen buenos ejemplos también!)
fuente
Para el análisis ráster, el paquete ráster es extremadamente poderoso. Además del manual estándar, hay algunas viñetas para comenzar.
fuente
No soy un usuario de PostGIS, pero después de sugerir polígonos de Voronoi para una pregunta del vecino más cercano , hice un poco de búsqueda. Descubrí que con R, puedes crear polígonos Voronoi para PostGIS . Estoy impresionado.
fuente
tess
objeto en unsp
objeto utilizando esta función proporcionada por Adrian Baddeley. Desde unsp
objeto, puede exportarlo a un shapefile si lo desea.Me topé con Spatial-Analyst.net . Muy informativo, completo y útil. Más específico a esta pregunta y en línea con algunas de las respuestas anteriores, vea esta página .
fuente
Vea también aquí cómo disfrutar del análisis estadístico de alta calidad en GRASS: http://grass.osgeo.org/wiki/R_statistics
fuente
Con esta función, puede hacer fácilmente uniones espaciales, pero solo si todas las áreas están llenas de polígonos.
fuente
Ejemplo de análisis de patrón de puntos:
Crea un patrón de puntos y lo representa. El paquete spatstat tiene varias funciones para analizar datos geográficos. Aquí hay algunos tutoriales de statstat :
fuente
No estoy seguro de si esto califica como un "truco", pero soy un gran admirador de la combinación del
acs
paquete (para seleccionar datos del Censo de EE. UU.) Y elleaflet
paquete (para hacer mapas interactivos de JavaScript que se pueden alojar en línea).Este tutorial hace un excelente trabajo ilustrando el beneficio de usar estos dos paquetes juntos.
fuente