[GRASS-user] Raster boxplots selection

sibylle.stoeckli at gmx.ch sibylle.stoeckli at gmx.ch
Wed Aug 7 01:23:54 PDT 2024


Dear coummunity



I am wondering if my code for selecting specific raster points (e.g.
squares, Quadrate) has a bug. I do not want to plot all squares, because I
still have not yet data for all squares.



v.extract quadrat out=quadrat_3 where="ID_Quadrat=527190 and
ID_Quadrat=569228 and ID_Quadrat=581166" --o



The resason is that I get an error on the next line



v.to.rast quadrat_3 out=quadrat_3 use=attr attribute_column=ID_Quadrat --o


ERROR: Database connection not defined for layer <1>



Full Code:

r.import input=allema_LZall_UAA_pres_merged.tif output=al_pres --o

g.region -ap rast=al_pres

v.import input=untersuchungsquadrate_2020_2024.gpkg output=quadrat --o

v.info -g quadrat



v.extract quadrat out=quadrat_3 where="ID_Quadrat=527190 and
ID_Quadrat=569228 and ID_Quadrat=581166" --o

v.to.rast quadrat_3 out=quadrat_3 use=attr attribute_column=ID_Quadrat --o

r.boxplot -o input=al_pres zones=quadrat_3 output="al_pres_boxplot_3.png"
raster_statistics=median -o



Thanks a lot

Sibylle

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/grass-user/attachments/20240807/59f4481c/attachment.htm>


More information about the grass-user mailing list