So what is in the three bands? Paul
On 11/07/2010 02:06 PM, Manuel Spínola wrote: > Species richness (from 3 to 20 more or less) > > On 07/11/2010 06:46 a.m., Paul Hiemstra wrote: >> What is the data captured by the tif file? >> >> Paul >> >> On 11/07/2010 12:45 PM, Manuel Spínola wrote: >>> Thank you very much Paul. >>> >>> This is the data frame that I created >>> >>> > head(riquezakrig_df) >>> band1 band2 band3 x y >>> 269 1.760396 1.620212 1.272875 460590.1 1216391 >>> 1686 1.760396 1.620212 1.272875 460540.1 1216341 >>> 1687 1.760396 1.620212 1.272875 460590.1 1216341 >>> 3103 1.760396 1.620212 1.272875 460490.1 1216291 >>> 3104 1.760396 1.620212 1.272875 460540.1 1216291 >>> 3105 1.760396 1.620212 1.272875 460590.1 1216291 >>> >>> >>> You suggest to make a new column with the color, but I don't >>> understand which information should I provide to the new column. >>> >>> Best, >>> >>> Manuel >>> >>> >>> On 07/11/2010 04:00 a.m., Paul Hiemstra wrote: >>>> Hi Manuel, >>>> >>>> Please do not respond to me alone, include the mailing list in the >>>> cc. You waste the option of someone else answering your question >>>> and the answers are not stored in the archives for other people. >>>> >>>> But back to your question. The tiff file has 3 bands, rgb bands. >>>> Probably the contents in those bands is a number between 0 and 255, >>>> but you do not specify this in your mail. In this case the tiff >>>> file does not contain the 'real' values but is a picture. So >>>> drawing individual bands does not make sense. What you could try is >>>> to construct a new column in the dataset with the color and try to >>>> use in the aes() argument (col = col, where col is the name of the >>>> column with the color). >>>> >>>> cheers, >>>> Paul >>>> >>>> >>>> On 11/06/2010 03:03 PM, Paul Hiemstra wrote: >>>>> Gpplot cannot plot SpatialPixelsDataFrame's in one go. You need to >>>>> convert the grid to a data frame. >>>>> >>>>> An example using the meuse dataset: >>>>> >>>>> library(sp) >>>>> library(ggplot2) >>>>> data(meuse.grid) >>>>> gridded(meuse.grid) = ~x+y >>>>> summary(meuse.grid) >>>>> >>>>> meuse.grid_df = as(meuse.grid, "data.frame") >>>>> >>>>> # Now make a plot with ggplot >>>>> p = ggplot(aes(x = x, y = y, fill = dist), data = meuse.grid_df) + >>>>> geom_tile() >>>>> print(p) >>>>> >>>>> # Add some tricks specially suitable for >>>>> # spatial plots >>>>> print(p + >>>>> scale_x_continuous('', labels = NA, breaks = NA) + >>>>> scale_y_continuous('', labels = NA, breaks = NA) + >>>>> opts(aspect.ratio = 1, alpha = 0.5)) >>>>> >>>>> cheers, >>>>> Paul >>>>> >>>>> On 11/06/2010 02:28 PM, Manuel Spínola wrote: >>>>> >>>>>> Thank you very much Paul. >>>>>> >>>>>> With plot I got a black polygon and with spplpot I got: >>>>>> >>>>>> >>>>>>> spplot(riquezakrig) >>>>>>> >>>>>> Error: cannot allocate vector of size 22.6 Mb >>>>>> >>>>>> with image show me what I want but no legend. >>>>>> >>>>>> Any specific geom to plot with ggplot? >>>>>> >>>>>> Best, >>>>>> >>>>>> Manuel >>>>>> >>>>>> On 06/11/2010 06:46 a.m., Paul Hiemstra wrote: >>>>>> >>>>>>> On 11/06/2010 11:15 AM, Manuel Spínola wrote: >>>>>>> >>>>>>>> Dear list members, >>>>>>>> >>>>>>>> I read a .tif file with readGDAL. >>>>>>>> What is the best way to display it?. >>>>>>>> >>>>>>>> >>>>>>> spplot or ggplot >>>>>>> >>>>>>>> I tried function image but there is no legend. >>>>>>>> >>>>>>>> >>>>>>> true >>>>>>> >>>>>>>> I tried plot and spplot but I couldn't get them to work. >>>>>>>> >>>>>>>> >>>>>>> why not? What error? What version of R and sp? You're not giving us >>>>>>> any information that could give us a clue to what's going wrong. >>>>>>> >>>>>>> cheers, >>>>>>> Paul >>>>>>> >>>>>>>> Any help will be appreciated. >>>>>>>> >>>>>>>> Best, >>>>>>>> >>>>>>>> Manuel >>>>>>>> >>>>>>>> >>>>>>>> >>>>>>>> >>>>>>>> >>>>>>>> _______________________________________________ >>>>>>>> R-sig-Geo mailing list >>>>>>>> R-sig-Geo@stat.math.ethz.ch >>>>>>>> https://stat.ethz.ch/mailman/listinfo/r-sig-geo >>>>>>>> >>>>>>>> >>>>>>> -- >>>>>>> Drs. Paul Hiemstra >>>>>>> Department of Physical Geography >>>>>>> Faculty of Geosciences >>>>>>> University of Utrecht >>>>>>> Heidelberglaan 2 >>>>>>> P.O. Box 80.115 >>>>>>> 3508 TC Utrecht >>>>>>> Phone: +3130 253 5773 >>>>>>> http://intamap.geo.uu.nl/~paul >>>>>>> http://nl.linkedin.com/pub/paul-hiemstra/20/30b/770 >>>>>>> >>>>>> -- >>>>>> *Manuel Spínola, Ph.D.* >>>>>> Instituto Internacional en Conservación y Manejo de Vida Silvestre >>>>>> Universidad Nacional >>>>>> Apartado 1350-3000 >>>>>> Heredia >>>>>> COSTA RICA >>>>>> mspin...@una.ac.cr >>>>>> mspinol...@gmail.com >>>>>> Teléfono: (506) 2277-3598 >>>>>> Fax: (506) 2237-7036 >>>>>> Personal website: Lobito de río >>>>>> <https://sites.google.com/site/lobitoderio/> >>>>>> Institutional website: ICOMVIS<http://www.icomvis.una.ac.cr/> >>>>>> >>>>> >>>>> >>>>> >>>>> _______________________________________________ >>>>> R-sig-Geo mailing list >>>>> R-sig-Geo@stat.math.ethz.ch >>>>> https://stat.ethz.ch/mailman/listinfo/r-sig-geo >>>>> >>>> >>>> >>>> -- >>>> Drs. Paul Hiemstra >>>> Department of Physical Geography >>>> Faculty of Geosciences >>>> University of Utrecht >>>> Heidelberglaan 2 >>>> P.O. Box 80.115 >>>> 3508 TC Utrecht >>>> Phone: +3130 253 5773 >>>> http://intamap.geo.uu.nl/~paul >>>> http://nl.linkedin.com/pub/paul-hiemstra/20/30b/770 >>> >>> >>> -- >>> *Manuel Spínola, Ph.D.* >>> Instituto Internacional en Conservación y Manejo de Vida Silvestre >>> Universidad Nacional >>> Apartado 1350-3000 >>> Heredia >>> COSTA RICA >>> mspin...@una.ac.cr >>> mspinol...@gmail.com >>> Teléfono: (506) 2277-3598 >>> Fax: (506) 2237-7036 >>> Personal website: Lobito de río >>> <https://sites.google.com/site/lobitoderio/> >>> Institutional website: ICOMVIS <http://www.icomvis.una.ac.cr/> >> >> >> -- >> Drs. Paul Hiemstra >> Department of Physical Geography >> Faculty of Geosciences >> University of Utrecht >> Heidelberglaan 2 >> P.O. Box 80.115 >> 3508 TC Utrecht >> Phone: +3130 253 5773 >> http://intamap.geo.uu.nl/~paul >> http://nl.linkedin.com/pub/paul-hiemstra/20/30b/770 > > > -- > *Manuel Spínola, Ph.D.* > Instituto Internacional en Conservación y Manejo de Vida Silvestre > Universidad Nacional > Apartado 1350-3000 > Heredia > COSTA RICA > mspin...@una.ac.cr > mspinol...@gmail.com > Teléfono: (506) 2277-3598 > Fax: (506) 2237-7036 > Personal website: Lobito de río > <https://sites.google.com/site/lobitoderio/> > Institutional website: ICOMVIS <http://www.icomvis.una.ac.cr/> -- Drs. Paul Hiemstra Department of Physical Geography Faculty of Geosciences University of Utrecht Heidelberglaan 2 P.O. Box 80.115 3508 TC Utrecht Phone: +3130 253 5773 http://intamap.geo.uu.nl/~paul http://nl.linkedin.com/pub/paul-hiemstra/20/30b/770 [[alternative HTML version deleted]]
_______________________________________________ R-sig-Geo mailing list R-sig-Geo@stat.math.ethz.ch https://stat.ethz.ch/mailman/listinfo/r-sig-geo