<div dir="ltr"><div></div><div>I guess Grass GIS can create interactive scatter plot from my data. Because, in other menu on Grass GIS (not in interactive scatter plot), it can create a scatter plot. <br></div><div><br></div><div>I will post new email in new title as your suggestion about interactive scatter plot.<br></div><div><br></div><div>Thank you for your response.<br></div><div><br></div><div>Regards, <br></div><div>mega<br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Mon, Dec 28, 2020 at 3:28 PM Micha Silver <<a href="mailto:tsvibar@gmail.com" target="_blank">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">
<div style="direction:ltr">
<p><br>
</p>
<div>On 12/28/2020 1:38 AM, mega saputra
wrote:<br>
</div>
<blockquote type="cite">
<div dir="ltr">
<div>Here is the output :</div>
<div>g.proj -p
</div>
</div>
</blockquote>
<p><br>
</p>
<p>.....</p>
<p>....<br>
</p>
<blockquote type="cite">
<div dir="ltr">
<div>But, there is no scatter plot in Map Display. Then I click
twice in vector > scatter_1_3 in Data tab. Then there is a
message : <br>
</div>
<div>"Unable to zoom to vector map <scatter_1_3@PERMANENT></div>
<div>Details: Illegal latitude for North: 65528"</div>
<div><br>
</div>
</div>
</blockquote>
<p>Yes, that's correct. You need to understand how r.scatterplot
works. The X-Y coordinates in the newly created point vector layer
are the raster values at each pixel of the original rasters. In
this case, since the original are color bands, the values
represent reflectance for each band. These values can go as high
as 65528, which is way beyond the range in a latitude/longitude
coordinate system (maximum long/lat of 360 X 180).</p>
<p><br>
</p>
<p>If you want to carry on with this (in your next mail you are
beginning with the interactive scatter plot), it might be helpful
to export the X-Y "coordinates" to a text file for plotting in
some other software. Just keep in mind that you will get a list of
12,000,000 rows with two columns of numbers-the reflectances of
both bands.</p>
<p> <br>
</p>
<blockquote type="cite">
<div dir="ltr">
<div>How to solve the problem?<br>
</div>
<div><br>
</div>
<div>Regards,</div>
<div>mega<br>
</div>
</div>
<br>
<div class="gmail_quote">
<div dir="ltr" class="gmail_attr">On Sun, Dec 27, 2020 at 3:22
PM Micha Silver <<a href="mailto:tsvibar@gmail.com" target="_blank">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">
<div style="direction:ltr">
<p>Here's my attempt:</p>
<p>I downloaded the rgb.tif that you linked to, then used it
to create a new, fresh location</p>
<p><br>
</p>
<p># Create a NEW location "mega" using the georeferenced
rgb.tif file to define location<br>
# and import the file</p>
<p># Run at the command line, *before* starting grass:<br>
</p>
<p>grass78 -c rgb.tif /home/micha/grass/mega</p>
<p><br>
</p>
<p># Now at the grass prompt:</p>
<p># Check projection settings<br>
g.proj -p<br>
-PROJ_INFO-------------------------------------------------<br>
name : WGS 84<br>
datum : wgs84<br>
ellps : wgs84<br>
proj : ll<br>
no_defs : defined<br>
-PROJ_EPSG-------------------------------------------------<br>
epsg : 4326<br>
-PROJ_UNITS------------------------------------------------<br>
unit : degree<br>
units : degrees<br>
meters : 1.0<br>
<br>
# Install the scatterplot extension<br>
g.extension r.scatterplot<br>
Downloading precompiled GRASS Addons
<r.scatterplot>...<br>
Updating extensions metadata file...<br>
Updating extension modules metadata file...<br>
Installation of <r.scatterplot> successfully
finished<br>
<br>
# Now run r.scatterplot on two bands<br>
r.scatterplot input=rgb.3,rgb.2 output=scatter_3_2<br>
100%<br>
Building topology for vector map
<scatter_3_2@PERMANENT>...<br>
Registering primitives...<br>
11830000<br>
# Almost 12 million points (!)<br>
</p>
<div><br>
</div>
<div><br>
</div>
<div>The X-Y coordinates of the points are the scatterplot
values.</div>
<div><br>
</div>
<div>HTH,</div>
<div>Micha<br>
</div>
<div><br>
</div>
<div>On 12/26/2020 11:40 PM, mega saputra wrote:<br>
</div>
<blockquote type="cite">
<div dir="ltr">
<div>The result is the same. It cannot add scatterplot.<br>
</div>
<div>Here is the output :</div>
<div><br>
</div>
<div>(Sun Dec 27 05:20:25 2020)
<br>
g.region -ap raster=rgb.1
<br>
projection: 3 (Latitude-Longitude)<br>
zone: 0<br>
datum: wgs84<br>
ellipsoid: wgs84<br>
north: 13:07:42.9N<br>
south: 27:38:15.5S<br>
west: 100:55:15.7E<br>
east: 131:15:38.1E<br>
nsres: 0:00:36.8<br>
ewres: 0:00:36.8<br>
rows: 3988<br>
cols: 2968<br>
cells: 11836384<br>
(Sun Dec 27 05:20:25 2020) Command finished (0 sec)
<br>
G__open(read): Unable to open
'/home/mega/newLocation/PERMAN<br>
ENT/.tmp/mega/vector/trAreas54950/frmt': No such file
or<br>
directory</div>
<div><br>
</div>
<div>If you want to try my data, my data is in : <a href="https://1drv.ms/u/s!AjXAWS5maHffhyH0q9oTos1umy3D?e=bKa7mX" target="_blank">https://1drv.ms/u/s!AjXAWS5maHffhyH0q9oTos1umy3D?e=bKa7mX</a></div>
<div><br>
</div>
<div>Thank you.</div>
<div><br>
</div>
<div>regards,</div>
<div>mega<br>
</div>
<div><br>
</div>
</div>
<br>
<div class="gmail_quote">
<div dir="ltr" class="gmail_attr">On Sat, Dec 26, 2020
at 9:39 PM Micha Silver <<a href="mailto:tsvibar@gmail.com" target="_blank">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"><br>
On 12/25/20 11:46 PM, mega saputra wrote:<br>
> Sorry, there is holiday.<br>
><br>
> I import the raster. Then set the region. This is
my syntax :<br>
><br>
> g.region -up raster=rgb.1<br>
<br>
<br>
The "-u" flag does not set the current computational
region to your <br>
raster. It just reports what would happen if you do
reset the region.<br>
<br>
Please try:<br>
<br>
<br>
g.region -ap rast=rgb.1<br>
<br>
<br>
and report back the output.<br>
<br>
<br>
<br>
> projection: 3 (Latitude-Longitude)<br>
> zone: 0<br>
> datum: wgs84<br>
> ellipsoid: wgs84<br>
> north: 13:07:42.9N<br>
> south: 27:38:15.5S<br>
> west: 100:55:15.7E<br>
> east: 131:15:38.1E<br>
> nsres: 0:00:36.8<br>
> ewres: 0:00:36.8<br>
> rows: 3988<br>
> cols: 2968<br>
> cells: 11836384<br>
> (Sat Dec 26 05:24:05 2020) Command finished (0
sec)<br>
><br>
> Then, from menu Imagery > Classify image >
Interactive input for <br>
> supervised classification. Then the output is
like this :<br>
><br>
> G__open(read): Unable to open
'/home/mega/newLocation/PERMAN<br>
> ENT/.tmp/mega/vector/trAreas114130/frmt': No such
file or<br>
> directory<br>
><br>
> I don't know why there is that message.<br>
> Then in Plots Menu, I choose Scatter Plots. I
click Add scatter plot. <br>
> In the Name of imagery group, I choose the
raster. In the Name of <br>
> imagery subgroup, I type the name. I click
Create/edit group. In x and <br>
> y axis, I choose a band. I click Add. Then there
is a message that :<br>
> "Scatter plot cannot be added.<br>
> Multiple of bands ranges
<x53.1@PERMANENT:65529 * <br>
> x53.1@PERMANENT:65529 = 4294049841 > is higher
than maximum limit <br>
> <16810000> "<br>
><br>
> In this case, I have defined the region. Thank
you for your attention.<br>
><br>
> Regards,<br>
> mega<br>
><br>
> On Wed, Dec 23, 2020 at 4:30 PM Micha Silver <<a href="mailto:tsvibar@gmail.com" target="_blank">tsvibar@gmail.com</a> <br>
> <mailto:<a href="mailto:tsvibar@gmail.com" target="_blank">tsvibar@gmail.com</a>>>
wrote:<br>
><br>
><br>
> On 12/23/20 5:09 AM, mega saputra wrote:<br>
> > Hello guys,<br>
> ><br>
> > I have MODIS data. It is sixteen bit
unsigned integer. When I<br>
> want to<br>
> > add scatter plot, there is a warning :<br>
> ><br>
> > "Scatter plot cannot be added.<br>
> > Multiple of bands ranges
<x53.1@PERMANENT:65529 *<br>
> > x53.1@PERMANENT:65529 = 4294049841 >
is higher than maximum limit<br>
> > <16810000> "<br>
><br>
><br>
> Please post your computation region (i.e. the
output of g.region<br>
> -p) and<br>
> the command you used.<br>
><br>
><br>
> ><br>
> > So, how to make big maximum limit on
Grass GIS, so that I can add<br>
> > scatter plot of MODIS data?<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>
<mailto:<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>
> <<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>
-- <br>
Micha Silver<br>
Ben Gurion Univ.<br>
Sde Boker, Remote Sensing Lab<br>
cell: +972-523-665918<br>
<br>
</blockquote>
</div>
</blockquote>
<pre cols="72">--
Micha Silver
Remote Sensing Lab, Sde Boker
Ben Gurion University
cell: +972-523-665918</pre>
</div>
</blockquote>
</div>
</blockquote>
<pre cols="72">--
Micha Silver
Remote Sensing Lab, Sde Boker
Ben Gurion University
cell: +972-523-665918</pre>
</div>
</blockquote></div>