<div dir="ltr"><div dir="ltr"><br></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Wed, Aug 7, 2024 at 11:00 AM sibylle via grass-user <<a href="mailto:grass-user@lists.osgeo.org">grass-user@lists.osgeo.org</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div class="msg-5273633622732490447"><div lang="EN-GB" style="overflow-wrap: break-word;"><div class="m_-5273633622732490447WordSection1"><p class="m_-5273633622732490447MsoPlainText"><span lang="DE-CH">Dear coummunity<u></u><u></u></span></p><p class="m_-5273633622732490447MsoPlainText"><span lang="DE-CH"><u></u> <u></u></span></p><p class="m_-5273633622732490447MsoPlainText">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.<u></u><u></u></p><p class="m_-5273633622732490447MsoPlainText"><u></u> <u></u></p><p class="m_-5273633622732490447MsoPlainText"><span lang="DE-CH" style="background:yellow">v.extract quadrat out=quadrat_3 where="ID_Quadrat=527190 and ID_Quadrat=569228 and ID_Quadrat=581166" --o</span><span lang="DE-CH"><u></u><u></u></span></p><p class="m_-5273633622732490447MsoPlainText"><span lang="DE-CH"><u></u> </span></p></div></div></div></blockquote><div>The where statement will result in no quadrats selected. Use 'or' instead of 'and':</div><div><br></div><div><span lang="DE-CH" style="background:yellow">v.extract quadrat out=quadrat_3 where="ID_Quadrat=527190 or ID_Quadrat=569228 or ID_Quadrat=581166" --o</span></div><div><br></div><div>This will give you the three quadrats (or the quadrats that have those IDs, in case there are more quadrats with the same ID). <br></div><div><br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div class="msg-5273633622732490447"><div lang="EN-GB" style="overflow-wrap: break-word;"><div class="m_-5273633622732490447WordSection1"><p class="m_-5273633622732490447MsoPlainText"><span lang="DE-CH"><u></u></span></p><p class="m_-5273633622732490447MsoPlainText">The resason is that I get an error on the next line<u></u><u></u></p><p class="m_-5273633622732490447MsoPlainText"><u></u> <u></u></p><p class="m_-5273633622732490447MsoPlainText">v.to.rast quadrat_3 out=quadrat_3 use=attr attribute_column=ID_Quadrat --o      <u></u><u></u></p><p class="m_-5273633622732490447MsoPlainText"><span style="background:yellow">ERROR: Database connection not defined for layer <1></span><u></u><u></u></p><p class="m_-5273633622732490447MsoPlainText"><u></u> <u></u></p><p class="m_-5273633622732490447MsoPlainText">Full Code:<u></u><u></u></p><p class="m_-5273633622732490447MsoPlainText"><span lang="DE-CH">r.import input=allema_LZall_UAA_pres_merged.tif output=al_pres --o<u></u><u></u></span></p><p class="m_-5273633622732490447MsoPlainText"><span lang="DE-CH">g.region -ap rast=al_pres<u></u><u></u></span></p><p class="m_-5273633622732490447MsoPlainText"><span lang="DE-CH">v.import input=untersuchungsquadrate_2020_2024.gpkg output=quadrat --o<u></u><u></u></span></p><p class="m_-5273633622732490447MsoPlainText"><span lang="DE-CH"><a href="http://v.info" target="_blank">v.info</a> -g quadrat<u></u><u></u></span></p><p class="m_-5273633622732490447MsoPlainText"><span lang="DE-CH"><u></u> <u></u></span></p><p class="m_-5273633622732490447MsoPlainText"><span lang="DE-CH">v.extract quadrat out=quadrat_3 where="ID_Quadrat=527190 and ID_Quadrat=569228 and ID_Quadrat=581166" --o <u></u><u></u></span></p><p class="m_-5273633622732490447MsoPlainText">v.to.rast quadrat_3 out=quadrat_3 use=attr attribute_column=ID_Quadrat --o<u></u><u></u></p><p class="m_-5273633622732490447MsoPlainText">r.boxplot -o input=al_pres zones=quadrat_3 output="al_pres_boxplot_3.png" raster_statistics=median –o<u></u><u></u></p><p class="m_-5273633622732490447MsoPlainText"><u></u> <u></u></p><p class="m_-5273633622732490447MsoPlainText">Thanks a lot<u></u><u></u></p><p class="m_-5273633622732490447MsoPlainText">Sibylle<u></u><u></u></p></div></div>_______________________________________________<br>
grass-user mailing list<br>
<a href="mailto:grass-user@lists.osgeo.org" target="_blank">grass-user@lists.osgeo.org</a><br>
<a href="https://lists.osgeo.org/mailman/listinfo/grass-user" rel="noreferrer" target="_blank">https://lists.osgeo.org/mailman/listinfo/grass-user</a><br>
</div></blockquote></div></div>