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
[[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
