[R-es] Resumen de R-help-es, Vol 80, Envío 26
Gilsanz, Jose Luis
jluis.gilsanz en eu.jll.com
Mie Oct 21 12:36:41 CEST 2015
Gracias por la respuesta Carlos pero sigo sin ver el camino.
Con stat_density2d puedo representar la densidad de los puntos en una zona determinada, pero eso no es lo que yo quiero.
Lo que quiero es mostrar el valor de z (independientemente de que haya mucho o pocos puntos en la zona)
El problema lo tengo a la hora de comprender como definir bien el aes para que el mapa de calor muestre una escala (la escala tim.colors(n = 64)) de los valores de z y no la cantidad de puntos.
> -----Mensaje original-----
> De: R-help-es [mailto:r-help-es-bounces en r-project.org] En nombre de r-
> help-es-request en r-project.org
> Enviado el: martes, 20 de octubre de 2015 22:36
> Para: r-help-es en r-project.org
> Asunto: Resumen de R-help-es, Vol 80, Envío 26
>
> Envíe los mensajes para la lista R-help-es a
> r-help-es en r-project.org
>
> Para subscribirse o anular su subscripción a través de la WEB
> https://stat.ethz.ch/mailman/listinfo/r-help-es
>
> O por correo electrónico, enviando un mensaje con el texto "help" en el
> asunto (subject) o en el cuerpo a:
> r-help-es-request en r-project.org
>
> Puede contactar con el responsable de la lista escribiendo a:
> r-help-es-owner en r-project.org
>
> Si responde a algún contenido de este mensaje, por favor, edite la linea del
> asunto (subject) para que el texto sea mas especifico que:
> "Re: Contents of R-help-es digest...". Además, por favor, incluya en la
> respuesta sólo aquellas partes del mensaje a las que está respondiendo.
>
>
> Asuntos del día:
>
> 1. Re: Mapa de Calor con Google Maps de fondo
> (Carlos J. Gil Bellosta )
>
>
> ----------------------------------------------------------------------
>
> Message: 1
> Date: Tue, 20 Oct 2015 22:35:37 +0200
> From: "Carlos J. Gil Bellosta " <cgb en datanalytics.com>
> To: "Gilsanz, Jose Luis" <jluis.gilsanz en eu.jll.com>
> Cc: "r-help-es en r-project.org" <r-help-es en r-project.org>
> Subject: Re: [R-es] Mapa de Calor con Google Maps de fondo
> Message-ID:
> <CADg83eefbuPD7bkSX08toLi62mo--kFGxWEVfLP6U28V-
> fvxrw en mail.gmail.com>
> Content-Type: text/plain; charset="utf-8"
>
> Hola, ¿qué tal?
>
> Tienes que hacer algo así como
>
> ggmap(mapa) + stat_density2d([...])
>
> Tienes ejemplos completos en el artículo de Wickham
>
> https://journal.r-project.org/archive/2013-1/kahle-wickham.pdf
>
> Un saludo,
>
> Carlos J. Gil Bellosta
> http://www.datanalytics.com
>
> El 20 de octubre de 2015, 17:39, Gilsanz, Jose Luis < jluis.gilsanz en eu.jll.com>
> escribió:
>
> > Hola:
> >
> >
> >
> > Estoy intentando generar un mapa de calor que muestre, dadas unas
> > coordenadas geográficas determinadas, el ?calor? (cuanto mas alto el
> > valor de esta variable mas rojo y cuanto mas bajo mas azul) de una
> > determinada variable.
> >
> >
> >
> > Es decir con unos datos como estos:
> >
> > x y z
> >
> > -3,716575 40,41743 3.169
> >
> > -3,710735 40,42179 4.134
> >
> > -3,704625 40,40333 2.606
> >
> > -3,703929 40,40363 2.563
> >
> > -3,686431 40,36133 1.452
> >
> > -3,682371 40,35542 1.523
> >
> > -3,675325 40,41986 4.122
> >
> > -3,673758 40,41978 4.098
> >
> > -3,645127 40,42780 2.306
> >
> >
> >
> > Donde el data frame datos contiene:
> >
> > x: Longitud
> >
> > y: Latitud
> >
> > z: Variable a estudiar
> >
> >
> >
> > lo que quiero es mostrar sobre un fondo de GoogleMaps, otra
> > capa(semitransparente) con el Calor de la variable.
> >
> >
> >
> > Para ello primero calculo una aproximación de superficies de mis
> > puntos usando MBA
> >
> >
> >
> > ##Cargamos paquetes
> >
> > library(MBA)
> >
> >
> >
> > ##Calculamos aproximacion de superficies segun coordenadas y valores
> >
> > superf <- mba.surf(datos, 2000, 2000,extend=FALSE)$xyz.est
> >
> >
> >
> > ##Obtenemos la capa con los colores a partir de la aproximación de
> > superficies
> >
> > image.plot(superf, bg= NULL, nlevel = param$NIVELES)
> >
> >
> >
> >
> >
> > Por otro lado obtengo el mapa de google maps sobre el cual deberia ir
> > esto
> > asi:
> >
> > #Cargamos paquete
> >
> > library(ggmap)
> >
> >
> >
> > ##Centro del mapa
> >
> > center <- c(mean(datos$x), mean(datos$y))
> >
> > ##Zoom
> >
> > zoom <- min(MaxZoom(range(datos$x), range(datos$x)))
> >
> >
> >
> > ##Obtenenmos mapa base
> >
> > mapa <- get_map(center, zoom = 12)
> >
> >
> >
> > Lo que ya no consigo hacer es poner la imagen con los colores sobre el
> > mapa de forma que este correctamente calibrada.
> >
> >
> >
> > Seguro que todo esto se pude hacer de un tiron con ggmap y/o
> > RgoogleMaps pero por mas que miro no veo la manera de meter el objeto
> superf ?dentro?
> > de un mapa.
> >
> >
> >
> > Muchas Gracias
> >
> >
> >
> >
> >
> > TASACIONES HIPOTECARIAS S.A.
> > Registration number: A-28/806222.
> > Registered Office: Pº de la Castellana, 79 - 1ª ; 28046 Madrid
> >
> > This e-mail is for the use of the intended recipient(s) only. If you
> > have received this e-mail in error, please notify the sender
> > immediately and then delete it. If you are not the intended recipient,
> > you must not use, disclose or distribute this e-mail without the author's
> prior permission.
> > We have taken precautions to minimise the risk of transmitting
> > software viruses, but we advise you to carry out your own virus checks
> > on any attachment to this message. We cannot accept liability for any
> > loss or damage caused by software viruses. If you are the intended
> > recipient and you do not wish to receive similar electronic messages
> > from us in future then please respond to the sender to this effect
> >
> > _______________________________________________
> > R-help-es mailing list
> > R-help-es en r-project.org
> > https://stat.ethz.ch/mailman/listinfo/r-help-es
> >
> ------------ próxima parte ------------
> Se ha borrado un adjunto en formato HTML...
> URL: <https://stat.ethz.ch/pipermail/r-help-
> es/attachments/20151020/7d4339f9/attachment.html>
> ------------ próxima parte ------------
> A non-text attachment was scrubbed...
> Name: image002.png
> Type: image/png
> Size: 265395 bytes
> Desc: no disponible
> URL: <https://stat.ethz.ch/pipermail/r-help-
> es/attachments/20151020/7d4339f9/attachment.png>
> ------------ próxima parte ------------
> A non-text attachment was scrubbed...
> Name: image001.png
> Type: image/png
> Size: 13312 bytes
> Desc: no disponible
> URL: <https://stat.ethz.ch/pipermail/r-help-
> es/attachments/20151020/7d4339f9/attachment-0001.png>
>
> ------------------------------
>
> Subject: Pié de página del digest
>
> _______________________________________________
> R-help-es mailing list
> R-help-es en r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-help-es
>
> ------------------------------
>
> Fin de Resumen de R-help-es, Vol 80, Envío 26
> *********************************************
Más información sobre la lista de distribución R-help-es