[R-sig-Geo] stplot - legend/classes of categorical vairable

[richard.redweik at posteo.de]

Thank you, Edzer!

That was exactly what I was looking for. It works as expected after 
coercing my variable to a factor!

Best,
Richard

Am 04.03.2015 19:12 schrieb Edzer Pebesma:
> Hi Richard, stplot as it is on CRAN does not support factor variables. 
> I
> added this support to the dev version on github, e.g. try this
> (reproducible!) example:
> 
> # install spacetime from github:
> devtools::install_github("edzer/spacetime")
> 
> library(spacetime)
> example(STFDF)
> set.seed(42)
> gridded(stfdf at sp) = TRUE
> stfdf$f = factor(sample(c("a", "b", "c"), 12, replace = TRUE))
> library(RColorBrewer)
> stplot(stfdf[,,"f"], col.regions = brewer.pal(3, "Accent"))
> 
> 
> Make sure that your variable IS a factor, e.g. by
> 
> fdf at crop_id = as.factor(fdf at crop_id)
> 
> Your data contains many points. If they are laid out on a grid, the 
> plot
> improves if you do
> 
> gridded(fdf at sp) = TRUE
> 
> so that stplot knows, too.
> 
> On 03/04/2015 11:45 AM, richard.redweik at posteo.de wrote:
>> Hello,
>> 
>> I have an object of class STFDF named 'fdf'. I want to create a stplot
>> in which the categorical variable crop_id resp. crop_short is 
>> visualized
>> per year.
>> 
>> Currently I am only able to assign a color/legend entry to a range of
>> crop_ids, e.g. green --> [1, 2], yellow --> (2, 3]:
>> 
>>> col <- brewer.pal(8, "Accent")
>>> color <- colorRampPalette(col)(18)
>>> stplot(fdf[,, "crop_id"], col.regions=color, cuts=18)
>> 
>> However, this does mean, that the crop_ids 1 and 2 are assigned to 
>> green.
>> When changing the cuts parameter to cuts=19, I get an unambigious but
>> not very readable assignment, e.g. green --> [1, 1.947], yellow -->
>> (1.947, 2.895].
>> 
>> Now I am wondering, what is the best way to get one class resp. legend
>> entry per crop_id, e.g. green --> 1, yellow --> 2 etc.
>> 
>> Does someone know how to do this?
>> 
>> Thanks in advance,
>> Richard
>> 
>> PS:
>> Summary of my fdf object:
>> 
>>> summary(fdf)
>> Object of class STFDF
>>  with Dimensions (s, t, attr): (17511, 11, 6)
>> [[Spatial:]]
>> Object of class SpatialPoints
>> Coordinates:
>>       min     max
>> x 3464500 3654500
>> y 2828500 3026500
>> Is projected: NA
>> proj4string : [NA]
>> Number of points: 17511
>> [[Temporal:]]
>>      Index                       timeIndex
>>  Min.   :1995-01-01 00:00:00   Min.   : 1.0
>>  1st Qu.:1997-07-02 12:00:00   1st Qu.: 3.5
>>  Median :2000-01-01 00:00:00   Median : 6.0
>>  Mean   :2000-01-01 08:43:38   Mean   : 6.0
>>  3rd Qu.:2002-07-02 12:00:00   3rd Qu.: 8.5
>>  Max.   :2005-01-01 00:00:00   Max.   :11.0
>> [[Data attributes:]]
>>     objectid          crop_id        crop_short         year           
>> prob
>>  Min.   :3222754   Min.   : 1.00   WSWH   :79136   Min.   :1995   Min. 
>>   :1
>>  1st Qu.:3255362   1st Qu.: 3.00   MAIF   :62594   1st Qu.:1997   1st 
>> Qu.:1
>>  Median :3287732   Median :10.00   RAPE   :13881   Median :2000   
>> Median :1
>>  Mean   :3289614   Mean   :11.56   PULS   :10785   Mean   :2000   Mean 
>>   :1
>>  3rd Qu.:3322247   3rd Qu.:19.00   WBAR   :10095   3rd Qu.:2003   3rd 
>> Qu.:1
>>  Max.   :3367519   Max.   :19.00   MAIZ   : 5069   Max.   :2005   Max. 
>>   :1
>>                                    (Other):11061
>>  nuts_code
>>  FR25:192621
>> 
>> _______________________________________________
>> R-sig-Geo mailing list
>> R-sig-Geo at r-project.org
>> https://stat.ethz.ch/mailman/listinfo/r-sig-geo
> 
> _______________________________________________
> R-sig-Geo mailing list
> R-sig-Geo at r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-sig-geo