<div dir="ltr"><div>Finally, I succeed to display scatter plot. I display scatter plot from Map Display > Create bivariate scatterplot of raster maps. Then I choose two band. You see in the attachment. <br></div><div>But, in this scatter plot, I can not create a polygon like in <a href="https://grass.osgeo.org/grass78/manuals/wxGUI.iscatt.html">https://grass.osgeo.org/grass78/manuals/wxGUI.iscatt.html</a>. <br></div><div>I want to create region of interest in the scatter plot link with the raster like in the website page. <br></div><div>Creating scatter plot from the website page will return error. It is about maximum limit. <br></div><div>How to solve this?</div><div><br></div><div>Regards,</div><div>mega<br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Mon, Dec 28, 2020 at 7:38 AM mega saputra <<a href="mailto:megasaputra8@gmail.com">megasaputra8@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 dir="ltr"><div>Here is the output :</div><div>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>(Mon Dec 28 07:25:33 2020) Command finished (0 sec) <br>(Mon Dec 28 07:25:42 2020) <br>g.extension r.scatterplot <br>Fetching <r.scatterplot> from GRASS GIS Addons repository (be patient)...<br>Compiling...<br>Installing...<br>Updating addons metadata file...<br>Installation of <r.scatterplot> successfully finished<br>(Mon Dec 28 07:25:50 2020) Command finished (7 sec) <br>(Mon Dec 28 07:26:18 2020) <br>r.scatterplot input=rgb.1,rgb.3 output=scatter_1_3 <br>Building topology for vector map <scatter_1_3@PERMANENT>...<br>Registering primitives...<br>(Mon Dec 28 07:28:40 2020) Command finished (2 min 21 sec) <br>(Mon Dec 28 07:30:53 2020) <br><a href="http://v.info" target="_blank">v.info</a> map=scatter_1_3@PERMANENT <br> +----------------------------------------------------------------------------+<br> | Name: scatter_1_3 |<br> | Mapset: PERMANENT |<br> | Location: newLocation |<br> | Database: /home/mega |<br> | Title: |<br> | Map scale: 1:1 |<br> | Name of creator: mega |<br> | Organization: |<br> | Source date: Mon Dec 28 07:26:19 2020 |<br> | Timestamp (first layer): none |<br> |----------------------------------------------------------------------------|<br> | Map format: native |<br> |----------------------------------------------------------------------------|<br> | Type of map: vector (level: 2) |<br> | |<br> | Number of points: 11836384 Number of centroids: 0 |<br> | Number of lines: 0 Number of boundaries: 0 |<br> | Number of areas: 0 Number of islands: 0 |<br> | Number of faces: 0 Number of kernels: 0 |<br> | Number of volumes: 0 Number of holes: 0 |<br> | |<br> | Map is 3D: Yes |<br> | Number of dblinks: 0 |<br> | |<br> | Projection: Latitude-Longitude |<br> | |<br> | N: 65528N S: 0 |<br> | E: 65528E W: 0 |<br> | B: 0 T: 0 |<br> | |<br> | Digitization threshold: 0 |<br> | Comment: |<br> | |<br> +----------------------------------------------------------------------------+<br>(Mon Dec 28 07:30:53 2020) Command finished (0 sec) <br></div><div><br></div><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>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></div>