<div dir="ltr"><div>The answer for your question : <br></div><div>1. I really input band1 for x axis and band3 for y axis. I don't know the message tell rgb1*rgb1. May be it just gives information about dimension. You can see in screenshot-1. <br></div><div>2. Actually I don't want to downscale. Because the digital number will vague.</div><div><br></div><div>Btw I have tried and I rescale the data with recommendation of GRASS GIS. Recommendation limit is 4000000. So, I rescale to 2000.</div><div><br></div><div>(Thu Dec 31 13:02:58 2020) <br>r.rescale input=rgb.5 output=rgb.5.rescale to=0,1999 --o <br>Rescale rgb.5[0,65528] to rgb.5.rescale[0,1999]<br>(Thu Dec 31 13:02:59 2020) Command finished (0 sec) <br>(Thu Dec 31 13:03:09 2020) <br>r.rescale input=rgb.4 output=rgb.4.rescale to=0,1999 --o <br>Rescale rgb.4[0,65528] to rgb.4.rescale[0,1999]<br>(Thu Dec 31 13:03:10 2020) Command finished (0 sec) <br>(Thu Dec 31 13:03:20 2020) <br>r.rescale input=rgb.3 output=rgb.3.rescale to=0,1999 --o <br>Rescale rgb.3[0,65528] to rgb.3.rescale[0,1999]<br>(Thu Dec 31 13:03:21 2020) Command finished (0 sec) <br>(Thu Dec 31 13:03:41 2020) <br>r.rescale input=rgb.2 output=rgb.2.rescale to=0,1999 --o <br>Rescale rgb.2[0,65528] to rgb.2.rescale[0,1999]<br>(Thu Dec 31 13:03:42 2020) Command finished (0 sec) <br>(Thu Dec 31 13:03:54 2020) <br>r.rescale input=rgb.1 output=rgb.1.rescale to=0,1999 --o <br>Rescale rgb.1[0,65528] to rgb.1.rescale[0,1999]<br>(Thu Dec 31 13:03:55 2020) Command finished (0 sec</div><div><br></div><div>Then I group the result with command i.group.</div><div><br></div><div>Then, from menu Imagery > Classify image > Interactive input for supervised classification. Then I succeed to create interactive scatter plot (screenshot-2).</div><div><br></div><div>In this window I can create polygon in scatter plot, but no effect with the raster like in <a href="https://grass.osgeo.org/grass78/manuals/wxGUI.iscatt.html">https://grass.osgeo.org/grass78/manuals/wxGUI.iscatt.html</a>. May be I will ask in new topic in new email.</div><div><br></div><div>Thank you.<br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Tue, Dec 29, 2020 at 6:49 PM Micha Silver <<a href="mailto:tsvibar@gmail.com">tsvibar@gmail.com</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">Hello Mega<br>
<br>
<br>
On 12/29/20 1:51 AM, mega saputra wrote:<br>
> Hello guys,<br>
><br>
> I have tried to add (passive) scatter plot with my data, and it is <br>
> succeed.<br>
> But, I want to add in interactive scatter plot, and not succeed. The <br>
> message :<br>
> "Scatter plot cannot be added. Multiple of bands ranges <br>
> <rgb.1@PERMANENT:65529 * rgb.1@PERMANENT:65529 = 4294049841 > is <br>
> higher than maximum limit <16810000> "<br>
><br>
<br>
Before the interactive scatterplot, can you clarify two questions?<br>
<br>
Why are you trying to plot a band against itself (rgb.1, rgb.1 above) ? <br>
This will give you 12,000,000 rows where each row contains two identical <br>
values. So the scatterplot will be essentially a line along 45º.<br>
<br>
<br>
Second, with rasters of about 3000x4000 you will get a huge number of <br>
points in the scatterplot. Would you consider to downscale to a lower <br>
resolution to make the resulting scatterplot more manageable?<br>
<br>
<br>
Now some suggestions:<br>
<br>
<br>
The range of reflectance values in your RGB image are up from 0 to <br>
65528. But the possible X-Y values in the Long/lat location need to be <br>
between 0-90 maximum in order to display the scatterplot. So here's a <br>
way to rescale the reflectance values to display the scatterplot on the <br>
Map Display.<br>
<br>
<br>
# ----------------------------------------------------------<br>
# 1- Create new location based on the CRS of the input GeoTiff<br>
# and import to grass<br>
# ----------------------------------------------------------<br>
grass -c rgb.tif /home/micha/GIS/grass/mega<br>
r.import input=rgb.tif output=rgb --overwrite<br>
# imports 5 bands<br>
<br>
# Check univariate statistics of one of the bands<br>
r.univar rgb.1<br>
100%<br>
total null and non-null cells: 11836384<br>
total null cells: 0<br>
<br>
Of the non-null cells:<br>
----------------------<br>
n: 11836384<br>
minimum: 0<br>
maximum: 65528<br>
range: 65528<br>
mean: 3982.64<br>
mean of absolute values: 3982.64<br>
standard deviation: 4159.48<br>
variance: 1.73012e+07<br>
variation coefficient: 104.44 %<br>
sum: 47140019279<br>
<br>
# Notice that the maximum value can be as high as 65528<br>
# (maximum of unsigned short integer)<br>
<br>
<br>
# ----------------------------------------------------------<br>
# 2- In order to display the scatterplot in a Lon/Lat coordinate system<br>
# the max value cannot be higher than 90 (maximum latitude)<br>
# So create rescaled rasters that "squeeze" reflectance values to the <br>
range 0,90<br>
# ----------------------------------------------------------<br>
r.rescale input=rgb.1 output=rgb.1.rescale to=0,90 --o<br>
r.rescale input=rgb.2 output=rgb.2.rescale to=0,90 --o<br>
r.rescale input=rgb.3 output=rgb.3.rescale to=0,90 --o<br>
r.rescale input=rgb.4 output=rgb.4.rescale to=0,90 --o<br>
r.rescale input=rgb.5 output=rgb.5.rescale to=0,90 --o<br>
<br>
# ----------------------------------------------------------<br>
# 3- Now prepare the scatterplot vector using the rescaled rasters, and <br>
display<br>
# 12 million points takes a while...<br>
<br>
# Here I choose the red and NIR bands for the scatterplot<br>
<br>
# ----------------------------------------------------------<br>
r.scatterplot input=rgb.3.rescale,rgb.4.rescale out=scatter_3_4 --o<br>
<br>
# See the attached image<br>
<br>
<br>
HTH<br>
<br>
<br>
> So, how to add in interactive scatter plot?<br>
><br>
> My data in : <a href="https://1drv.ms/u/s!AjXAWS5maHffhyH0q9oTos1umy3D?e=bKa7mX" rel="noreferrer" target="_blank">https://1drv.ms/u/s!AjXAWS5maHffhyH0q9oTos1umy3D?e=bKa7mX</a> <br>
> <<a href="https://1drv.ms/u/s!AjXAWS5maHffhyH0q9oTos1umy3D?e=bKa7mX" rel="noreferrer" target="_blank">https://1drv.ms/u/s!AjXAWS5maHffhyH0q9oTos1umy3D?e=bKa7mX</a>><br>
><br>
> regards,<br>
> mega<br>
><br>
> _______________________________________________<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>
<br>
-- <br>
Micha Silver<br>
Ben Gurion Univ.<br>
Sde Boker, Remote Sensing Lab<br>
cell: +972-523-665918<br>
<br>
</blockquote></div>