[R-sig-Geo] Problem to retrieve/plot selected bands in a stars proxy object

[Loïc Valéry]

Thank you so much Edzer for fixing this so quickly. Everything is working fine now.
Cheers,
Loïc


De : R-sig-Geo <r-sig-geo-bounces using r-project.org> de la part de Edzer Pebesma <edzer.pebesma using uni-muenster.de>
Envoyé : mardi 14 septembre 2021 17:45
À : r-sig-geo using r-project.org <r-sig-geo using r-project.org>
Objet : Re: [R-sig-Geo] Problem to retrieve/plot selected bands in a stars proxy object 
 
Thanks for the clear report, Loïc; this should now be fixed in the 
version on github.

On 14/09/2021 14:55, Loïc Valéry wrote:
> Dear list members,
> 
> I hope this message finds you well.
> I submit to you a problem with which I have been struggling for several days without finding the beginning of a solution.
> 
> I would like to select some bands of a 'stars proxy' object that I have previously created by merging two images with the function read_stars(c(image1, image2), proxy = true, along = "band") in order to plot them.
> 
> My problem is that R returns the following error message when trying to plot the selected bands:
> Error in x[[i]][, , 10:12, , drop = FALSE] : subscript out of bounds
> 
> Here is a short REPREX:
> library(stars)
> tif_1 <- system.file("tif/L7_ETMs.tif", package = "stars")
> tif_2 <- system.file("tif/L7_ETMs.tif", package = "stars")
> 
> tif_merge <- read_stars(c(tif_1,tif_2), proxy = TRUE, along = "band")
> 
> plot(tif_merge[,,,10:12], rgb = 1:3)
> #> Error in x[[i]][, , 10:12, , drop = FALSE]: indice hors limites
> 
> 
> 
> I guess the problem is that the "tif_merge" object has different dimensions before and after the promise evaluation:
> 
> dim(tif_merge)
> #>    x    y band
> #>  349  352   12
> 
> dim(st_as_stars(tif_merge))
> #>       x       y    band new_dim
> #>     349     352       6       2
> 
> However, this did not help me to find any solution!
> 
> Thank you in advance for your help.
> Cheers,
> Loïc
> 
> _______________________________________________
> R-sig-Geo mailing list
> R-sig-Geo using r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-sig-geo
> 

-- 
Edzer Pebesma
Institute for Geoinformatics
Heisenbergstrasse 2, 48151 Muenster, Germany
Phone: +49 251 8333081

_______________________________________________
R-sig-Geo mailing list
R-sig-Geo using r-project.org
https://stat.ethz.ch/mailman/listinfo/r-sig-geo