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

_______________________________________________
grass-user mailing list
[email protected]
https://lists.osgeo.org/mailman/listinfo/grass-user

Reply via email to