<!DOCTYPE html>
<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=UTF-8">
<style id="bidiui-paragraph-margins" type="text/css">body p { margin-bottom: 0cm; margin-top: 0pt; } </style>
</head>
<body bidimailui-detected-decoding-type="UTF-8"
bidimailui-charset-is-forced="true" style="direction: ltr;"
text="#933131" bgcolor="#FFFFFF">
<p><font size="4">Hi Sibylle:</font></p>
<p><font size="4">I'm not sure I fully understand the problem. If
the below example is not relevant, just ignore...</font></p>
<p><font size="4"><br>
</font></p>
<p><font size="4">I did a simple test using the nc_spm_08 location.
The PERMANENT mapset includes an elevation raster
'elev_srtm_30m' and a categorical landcover raster
'landclass96'. I used the landcover raster as zones, and
elevation for the values in the boxplots.</font></p>
<p><font size="5"><br>
</font></p>
<p><font size="5"><span style="font-family:monospace"><span
style="color:#000000;background-color:#ffffff;">g.region -ap
rast=elev_srtm_30m</span></span></font></p>
<p><font size="5"><span style="font-family:monospace"><span
style="color:#000000;background-color:#ffffff;"># All
landcover zones<br>
</span></span></font></p>
<p><font size="5"><span style="font-family:monospace"><span
style="color:#000000;background-color:#ffffff;"> r.boxplot
input=elev_srtm_30m zones=landclass96
raster_statistics=median output=~/work/tmp/boxplot.png --o</span></span></font></p>
<p><font size="5"><span style="font-family:monospace"><span
style="color:#000000;background-color:#ffffff;"># See
attached boxplot image<br>
</span></span></font></p>
<p><font size="5"><span style="font-family:monospace"><span
style="color:#000000;background-color:#ffffff;"># Now limit
to only one zone, the forests, class 5 in this raster:</span></span></font></p>
<p><font size="5"><span style="font-family:monospace"><span
style="color:#000000;background-color:#ffffff;">r.mapcalc
"MASK = if(landclass96 == 5, 1, null())"</span></span></font></p>
<p><font size="5"><span style="font-family:monospace"><span
style="color:#000000;background-color:#ffffff;">r.boxplot
input=elev_srtm_30m zones=landclass96
raster_statistics=median
output=~/work/tmp/boxplot_forest.png</span></span></font></p>
<p><font size="5"><span style="font-family:monospace"><span
style="color:#000000;background-color:#ffffff;"># See
attached boxplot_forest<br>
</span></span></font></p>
<p><font size="5"><span style="font-family:monospace"><span
style="color:#000000;background-color:#ffffff;">r.mask -r</span><br>
</span></font></p>
<p><span style="font-family:monospace"><span
style="color:#000000;background-color:#ffffff;"><br>
</span></span></p>
<p><font size="4">Hope this is somehow helpful.<br>
</font></p>
<p></p>
<p><span style="font-family:monospace"><span
style="color:#000000;background-color:#ffffff;"><br>
</span></span></p>
<p><span style="font-family:monospace"><span
style="color:#000000;background-color:#ffffff;"><br>
</span></span></p>
<div class="moz-cite-prefix">On 31/05/2024 10:35, sibylle via
grass-user wrote:<br>
</div>
<blockquote type="cite"
cite="mid:022501dab32d$20d23c00$6276b400$@gmx.ch">
<meta http-equiv="Content-Type" content="text/html; charset=UTF-8">
<meta name="Generator"
content="Microsoft Word 15 (filtered medium)">
<!--[if !mso]><style>v\:* {behavior:url(#default#VML);}
o\:* {behavior:url(#default#VML);}
w\:* {behavior:url(#default#VML);}
.shape {behavior:url(#default#VML);}
</style><![endif]-->
<style>@font-face
{font-family:Wingdings;
panose-1:5 0 0 0 0 0 0 0 0 0;}@font-face
{font-family:"Cambria Math";
panose-1:2 4 5 3 5 4 6 3 2 4;}@font-face
{font-family:Calibri;
panose-1:2 15 5 2 2 2 4 3 2 4;}@font-face
{font-family:Aptos;}p.MsoNormal, li.MsoNormal, div.MsoNormal
{margin:0cm;
font-size:12.0pt;
font-family:"Aptos",sans-serif;}a:link, span.MsoHyperlink
{mso-style-priority:99;
color:#0563C1;
text-decoration:underline;}p.MsoPlainText, li.MsoPlainText, div.MsoPlainText
{mso-style-priority:99;
mso-style-link:"Plain Text Char";
margin:0cm;
font-size:11.0pt;
font-family:"Calibri",sans-serif;
mso-ligatures:standardcontextual;
mso-fareast-language:EN-US;}p.MsoListParagraph, li.MsoListParagraph, div.MsoListParagraph
{mso-style-priority:34;
margin-top:0cm;
margin-right:0cm;
margin-bottom:0cm;
margin-left:36.0pt;
font-size:12.0pt;
font-family:"Aptos",sans-serif;}span.PlainTextChar
{mso-style-name:"Plain Text Char";
mso-style-priority:99;
mso-style-link:"Plain Text";
font-family:"Calibri",sans-serif;}span.EmailStyle22
{mso-style-type:personal-reply;
font-family:"Calibri",sans-serif;
color:windowtext;}.MsoChpDefault
{mso-style-type:export-only;
font-size:10.0pt;
mso-ligatures:none;}div.WordSection1
{page:WordSection1;}ol
{margin-bottom:0cm;}ul
{margin-bottom:0cm;}</style><!--[if gte mso 9]><xml>
<o:shapedefaults v:ext="edit" spidmax="1026" />
</xml><![endif]--><!--[if gte mso 9]><xml>
<o:shapelayout v:ext="edit">
<o:idmap v:ext="edit" data="1" />
</o:shapelayout></xml><![endif]-->
<div class="WordSection1">
<p class="MsoNormal"><span
style="font-size:11.0pt;font-family:"Calibri",sans-serif;mso-fareast-language:EN-US">Dear
community<o:p></o:p></span></p>
<p class="MsoNormal"><span
style="font-size:11.0pt;font-family:"Calibri",sans-serif;mso-fareast-language:EN-US"><o:p> </o:p></span></p>
<p class="MsoNormal"><span
style="font-size:11.0pt;font-family:"Calibri",sans-serif;mso-fareast-language:EN-US">I
tested the r.boxplot addon and it is definitively possible
to run the addon with more than 20
classes/regions/sites/squares.<o:p></o:p></span></p>
<p class="MsoNormal"><span
style="font-size:11.0pt;font-family:"Calibri",sans-serif;mso-fareast-language:EN-US">Furthermore
I realised the reason for the wrong assignment of class to
raster, but I am not yet able to fix it.<o:p></o:p></span></p>
<ul style="margin-top:0cm" type="disc">
<li class="MsoListParagraph"
style="margin-left:0cm;mso-list:l1 level1 lfo4"><span
style="font-size:11.0pt;font-family:"Calibri",sans-serif;mso-fareast-language:EN-US">The
challenge is, that (at the moment) I have raster data for
one class (one square), but in total I have 150 squares.
So at the moment it just selects the first square (either
by cat oder ID_Quadrat).<o:p></o:p></span></li>
<li class="MsoListParagraph"
style="margin-left:0cm;mso-list:l1 level1 lfo4"><span
style="font-size:11.0pt;font-family:"Calibri",sans-serif;mso-fareast-language:EN-US">It
seems however, if I have raster data for all squares,
r.boxplot is working fine.<o:p></o:p></span></li>
<li class="MsoListParagraph"
style="margin-left:0cm;mso-list:l1 level1 lfo4"><span
style="font-size:11.0pt;font-family:"Calibri",sans-serif;mso-fareast-language:EN-US">My
question: Is it possible to assign raster data for one
individual class/squares, of you have missing data for
other classes/squares?<o:p></o:p></span></li>
</ul>
<p class="MsoNormal"><span
style="font-size:11.0pt;font-family:"Calibri",sans-serif;mso-fareast-language:EN-US"><o:p> </o:p></span></p>
<p class="MsoNormal"><span
style="font-size:11.0pt;font-family:"Calibri",sans-serif;mso-fareast-language:EN-US">Kind
regards<o:p></o:p></span></p>
<p class="MsoNormal"><span
style="font-size:11.0pt;font-family:"Calibri",sans-serif;mso-fareast-language:EN-US">Sibylle<o:p></o:p></span></p>
<p class="MsoNormal"><span
style="font-size:11.0pt;font-family:"Calibri",sans-serif;mso-fareast-language:EN-US"><o:p> </o:p></span></p>
<p class="MsoNormal"><span
style="font-size:11.0pt;font-family:"Calibri",sans-serif;mso-fareast-language:EN-US">First
using this code here:<o:p></o:p></span></p>
<p class="MsoNormal">v.to.rast
input=Untersuchungsquadrate_2020_2024@PERMANENT
output=Quadrate use=attr <span
style="background:yellow;mso-highlight:yellow">attribute_column="ID_Quadrat"</span>
label_column="ID_Quadrat"<span
style="font-family:"Calibri",sans-serif"><o:p></o:p></span></p>
<p class="MsoNormal"><span lang="DE-CH">r.boxplot -o
input=U803174_abgeschnitten zones=Quadrate bx_sort=ascending
raster_statistics=median<o:p></o:p></span></p>
<p class="MsoNormal"><span
style="font-size:11.0pt;font-family:"Calibri",sans-serif;mso-fareast-language:EN-US"
lang="DE-CH"><o:p> </o:p></span></p>
<p class="MsoNormal"><img style="width:8.2812in;height:3.7083in"
id="Picture_x0020_5"
src="cid:part1.qn0UOO3Z.olFn5dKL@gmail.com" class=""
width="795" height="356"><span
style="font-size:11.0pt;font-family:"Calibri",sans-serif;mso-fareast-language:EN-US"
lang="DE-CH"><o:p></o:p></span></p>
<p class="MsoNormal"><span
style="font-size:11.0pt;font-family:"Calibri",sans-serif;mso-fareast-language:EN-US"><o:p> </o:p></span></p>
<p class="MsoNormal"><span
style="font-size:11.0pt;font-family:"Calibri",sans-serif;mso-fareast-language:EN-US"><o:p> </o:p></span></p>
<p class="MsoNormal"><span
style="font-size:11.0pt;font-family:"Calibri",sans-serif;mso-fareast-language:EN-US"><o:p> </o:p></span></p>
<p class="MsoNormal"><span
style="font-size:11.0pt;font-family:"Calibri",sans-serif;mso-fareast-language:EN-US">Second
using this code here:<o:p></o:p></span></p>
<p class="MsoNormal">v.to.rast
input=Untersuchungsquadrate_2020_2024@PERMANENT
output=Quadrate use=attr <span
style="background:yellow;mso-highlight:yellow">attribute_column=cat</span>
label_column="ID_Quadrat"<span
style="font-family:"Calibri",sans-serif"><o:p></o:p></span></p>
<p class="MsoNormal"><span lang="DE-CH">r.boxplot -o
input=U803174_abgeschnitten zones=Quadrate bx_sort=ascending
raster_statistics=median<o:p></o:p></span></p>
<p class="MsoNormal"><span
style="font-size:11.0pt;font-family:"Calibri",sans-serif;mso-fareast-language:EN-US"
lang="DE-CH"><o:p> </o:p></span></p>
<p class="MsoNormal"><span
style="font-size:11.0pt;font-family:"Calibri",sans-serif;mso-fareast-language:EN-US"
lang="DE-CH"><o:p> </o:p></span></p>
<p class="MsoNormal"><img style="width:10.375in;height:4.8854in"
id="Picture_x0020_6"
src="cid:part2.QIFRSvdP.1G6LywKC@gmail.com" class=""
width="996" height="469"><span
style="font-size:11.0pt;font-family:"Calibri",sans-serif;mso-fareast-language:EN-US"><o:p></o:p></span></p>
<div>
<div
style="border:none;border-top:solid #E1E1E1 1.0pt;padding:3.0pt 0cm 0cm 0cm">
<p class="MsoNormal"><b><span
style="font-size:11.0pt;font-family:"Calibri",sans-serif">From:</span></b><span
style="font-size:11.0pt;font-family:"Calibri",sans-serif">
<a class="moz-txt-link-abbreviated" href="mailto:sibylle.stoeckli@gmx.ch">sibylle.stoeckli@gmx.ch</a> <a class="moz-txt-link-rfc2396E" href="mailto:sibylle.stoeckli@gmx.ch"><sibylle.stoeckli@gmx.ch></a>
<br>
<b>Sent:</b> Thursday, May 30, 2024 9:27 AM<br>
<b>To:</b> '<a class="moz-txt-link-abbreviated" href="mailto:grass-user@lists.osgeo.org">grass-user@lists.osgeo.org</a>'
<a class="moz-txt-link-rfc2396E" href="mailto:grass-user@lists.osgeo.org"><grass-user@lists.osgeo.org></a><br>
<b>Subject:</b> RE: raster boxplots wrong assignment<o:p></o:p></span></p>
</div>
</div>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal"><span
style="font-size:11.0pt;font-family:"Calibri",sans-serif;mso-fareast-language:EN-US">Dear
community<o:p></o:p></span></p>
<p class="MsoNormal"><span
style="font-size:11.0pt;font-family:"Calibri",sans-serif;mso-fareast-language:EN-US"><o:p> </o:p></span></p>
<p class="MsoNormal"><span
style="font-size:11.0pt;font-family:"Calibri",sans-serif;mso-fareast-language:EN-US">I
am still wondering, if the r.boxplot addon is working with
just a limited number of areass/sites/regions?<o:p></o:p></span></p>
<p class="MsoNormal"><span
style="font-size:11.0pt;font-family:"Calibri",sans-serif;mso-fareast-language:EN-US"><o:p> </o:p></span></p>
<p class="MsoNormal"><span
style="font-size:11.0pt;font-family:"Calibri",sans-serif;mso-fareast-language:EN-US">Kind
regards<o:p></o:p></span></p>
<p class="MsoNormal"><span
style="font-size:11.0pt;font-family:"Calibri",sans-serif;mso-fareast-language:EN-US">Sibylle<o:p></o:p></span></p>
<p class="MsoNormal"><span
style="font-size:11.0pt;font-family:"Calibri",sans-serif;mso-fareast-language:EN-US"><o:p> </o:p></span></p>
<div>
<div
style="border:none;border-top:solid #E1E1E1 1.0pt;padding:3.0pt 0cm 0cm 0cm">
<p class="MsoNormal"><b><span
style="font-size:11.0pt;font-family:"Calibri",sans-serif">From:</span></b><span
style="font-size:11.0pt;font-family:"Calibri",sans-serif"> <a
href="mailto:sibylle.stoeckli@gmx.ch"
moz-do-not-send="true" class="moz-txt-link-freetext">sibylle.stoeckli@gmx.ch</a>
<<a href="mailto:sibylle.stoeckli@gmx.ch"
moz-do-not-send="true" class="moz-txt-link-freetext">sibylle.stoeckli@gmx.ch</a>>
<br>
<b>Sent:</b> Saturday, May 25, 2024 3:48 PM<br>
<b>To:</b> '<a class="moz-txt-link-abbreviated" href="mailto:grass-user@lists.osgeo.org">grass-user@lists.osgeo.org</a>' <<a
href="mailto:grass-user@lists.osgeo.org"
moz-do-not-send="true" class="moz-txt-link-freetext">grass-user@lists.osgeo.org</a>><br>
<b>Subject:</b> raster boxplots wrong assignment<o:p></o:p></span></p>
</div>
</div>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoPlainText"><span lang="DE-CH">Dear community<o:p></o:p></span></p>
<p class="MsoPlainText"><span lang="DE-CH"><o:p> </o:p></span></p>
<p class="MsoPlainText">I would like to plot boxplots from a
raster _<i>pres_apple_spring</i>_merged for different 1 km2
squares (38 squares, Fig. 1, layer “Untersuchungsquadrate”). I
was using the same code as used before for 12 regions:<o:p></o:p></p>
<p class="MsoPlainText"><a
href="https://ecodiv.earth/post/drawing-boxplots-of-raster-values/"
moz-do-not-send="true" class="moz-txt-link-freetext">https://ecodiv.earth/post/drawing-boxplots-of-raster-values/</a>
(Author Paulo van Breugel). All layers have the same
projection (EPSG:2056, CH1903+/LV95). However, when loading
the layer squares (Untersuchungsquadrate) GRASS GIS was asking
for reprojection.<o:p></o:p></p>
<p class="MsoPlainText"><o:p> </o:p></p>
<ul style="margin-top:0cm" type="disc">
<li class="MsoPlainText" style="mso-list:l0 level1 lfo3">I got
now a warning message, that I have a lot of categories,
however, the modelling time was fast (Fig 2)<o:p></o:p></li>
<li class="MsoPlainText" style="mso-list:l0 level1 lfo3">The
error is now, that the boxplot is not showing all squares
(e.g. ID_Quadrat 803174) and the boxplots are assigned to
the wrong square (ID_Quadrat) (Fig. 3)<o:p></o:p></li>
</ul>
<p class="MsoPlainText"><o:p> </o:p></p>
<p class="MsoPlainText">Kind regards<o:p></o:p></p>
<p class="MsoPlainText">Sibylle Stöckli<o:p></o:p></p>
<p class="MsoPlainText"><o:p> </o:p></p>
<p class="MsoPlainText">Fig. 1<o:p></o:p></p>
<p class="MsoPlainText"><o:p> </o:p></p>
<p class="MsoPlainText"><span style="mso-ligatures:none"><img
style="width:6.9166in;height:3.625in" id="Picture_x0020_1"
src="cid:part3.WV1wgxW8.WIVULOdS@gmail.com" class=""
width="664" height="348" border="0"></span><o:p></o:p></p>
<p class="MsoPlainText"><o:p> </o:p></p>
<p class="MsoPlainText">Fig. 2<o:p></o:p></p>
<p class="MsoPlainText"><o:p> </o:p></p>
<p class="MsoPlainText"><span style="mso-ligatures:none"><img
style="width:7.0416in;height:3.5in" id="Picture_x0020_3"
src="cid:part4.tbo7LG9I.TIObktR5@gmail.com" class=""
width="676" height="336" border="0"></span><span
lang="DE-CH"><o:p></o:p></span></p>
<p class="MsoPlainText"><span lang="DE-CH"><o:p> </o:p></span></p>
<p class="MsoPlainText"><span lang="DE-CH">Figure 3<o:p></o:p></span></p>
<p class="MsoPlainText"><span lang="DE-CH"><o:p> </o:p></span></p>
<p class="MsoPlainText"><span style="mso-ligatures:none"><img
style="width:9.6354in;height:2.8645in"
id="Picture_x0020_2"
src="cid:part5.Zr0deI1s.UQqkGoZF@gmail.com" class=""
width="925" height="275" border="0"></span><span
lang="DE-CH"><o:p></o:p></span></p>
<p class="MsoPlainText"><span lang="DE-CH"><o:p> </o:p></span></p>
<p class="MsoPlainText"><span lang="DE-CH"><o:p> </o:p></span></p>
</div>
<br>
<fieldset class="moz-mime-attachment-header"></fieldset>
<pre class="moz-quote-pre" wrap="">_______________________________________________
grass-user mailing list
<a class="moz-txt-link-abbreviated" href="mailto:grass-user@lists.osgeo.org">grass-user@lists.osgeo.org</a>
<a class="moz-txt-link-freetext" href="https://lists.osgeo.org/mailman/listinfo/grass-user">https://lists.osgeo.org/mailman/listinfo/grass-user</a>
</pre>
</blockquote>
<pre class="moz-signature" cols="72">--
Micha Silver
Ben Gurion Univ.
Sde Boker, Remote Sensing Lab
cell: +972-523-665918</pre>
</body>
</html>