[GRASS-user] Scatter plot cannot be added

GRASS GRASS user
1

Hello guys,

I have MODIS data. It is sixteen bit unsigned integer. When I want to add scatter plot, there is a warning :

"Scatter plot cannot be added.
Multiple of bands ranges <x53.1@PERMANENT:65529 * x53.1@PERMANENT:65529 = 4294049841 > is higher than maximum limit <16810000> "

So, how to make big maximum limit on Grass GIS, so that I can add scatter plot of MODIS data?

Regards,
mega

2

On 12/23/20 5:09 AM, mega saputra wrote:

Hello guys,

I have MODIS data. It is sixteen bit unsigned integer. When I want to add scatter plot, there is a warning :

"Scatter plot cannot be added.
Multiple of bands ranges <x53.1@PERMANENT:65529 * x53.1@PERMANENT:65529 = 4294049841 > is higher than maximum limit <16810000> "

Please post your computation region (i.e. the output of g.region -p) and the command you used.

So, how to make big maximum limit on Grass GIS, so that I can add scatter plot of MODIS data?

Regards,
mega

_______________________________________________
grass-user mailing list
grass-user@lists.osgeo.org
https://lists.osgeo.org/mailman/listinfo/grass-user

--
Micha Silver
Ben Gurion Univ.
Sde Boker, Remote Sensing Lab
cell: +972-523-665918

3

Sorry, there is holiday.

I import the raster. Then set the region. This is my syntax :

g.region -up raster=rgb.1
projection: 3 (Latitude-Longitude)
zone: 0
datum: wgs84
ellipsoid: wgs84
north: 13:07:42.9N
south: 27:38:15.5S
west: 100:55:15.7E
east: 131:15:38.1E
nsres: 0:00:36.8
ewres: 0:00:36.8
rows: 3988
cols: 2968
cells: 11836384
(Sat Dec 26 05:24:05 2020) Command finished (0 sec)

Then, from menu Imagery > Classify image > Interactive input for supervised classification. Then the output is like this :

G__open(read): Unable to open ‘/home/mega/newLocation/PERMAN
ENT/.tmp/mega/vector/trAreas114130/frmt’: No such file or
directory

I don’t know why there is that message.

Then in Plots Menu, I choose Scatter Plots. I click Add scatter plot. In the Name of imagery group, I choose the raster. In the Name of imagery subgroup, I type the name. I click Create/edit group. In x and y axis, I choose a band. I click Add. Then there is a message that :

"Scatter plot cannot be added.
Multiple of bands ranges <x53.1@PERMANENT:65529 * x53.1@PERMANENT:65529 = 4294049841 > is higher than maximum limit <16810000> "

In this case, I have defined the region. Thank you for your attention.

Regards,
mega

On Wed, Dec 23, 2020 at 4:30 PM Micha Silver <tsvibar@gmail.com> wrote:

On 12/23/20 5:09 AM, mega saputra wrote:

Hello guys,

I have MODIS data. It is sixteen bit unsigned integer. When I want to
add scatter plot, there is a warning :

"Scatter plot cannot be added.
Multiple of bands ranges <x53.1@PERMANENT:65529 *
x53.1@PERMANENT:65529 = 4294049841 > is higher than maximum limit
<16810000> "

Please post your computation region (i.e. the output of g.region -p) and
the command you used.

So, how to make big maximum limit on Grass GIS, so that I can add
scatter plot of MODIS data?

Regards,
mega

grass-user mailing list
grass-user@lists.osgeo.org
https://lists.osgeo.org/mailman/listinfo/grass-user

–
Micha Silver
Ben Gurion Univ.
Sde Boker, Remote Sensing Lab
cell: +972-523-665918

4

On 12/25/20 11:46 PM, mega saputra wrote:

Sorry, there is holiday.

I import the raster. Then set the region. This is my syntax :

g.region -up raster=rgb.1

The "-u" flag does not set the current computational region to your raster. It just reports what would happen if you do reset the region.

Please try:

g.region -ap rast=rgb.1

and report back the output.

projection: 3 (Latitude-Longitude)
zone: 0
datum: wgs84
ellipsoid: wgs84
north: 13:07:42.9N
south: 27:38:15.5S
west: 100:55:15.7E
east: 131:15:38.1E
nsres: 0:00:36.8
ewres: 0:00:36.8
rows: 3988
cols: 2968
cells: 11836384
(Sat Dec 26 05:24:05 2020) Command finished (0 sec)

Then, from menu Imagery > Classify image > Interactive input for supervised classification. Then the output is like this :

G__open(read): Unable to open '/home/mega/newLocation/PERMAN
ENT/.tmp/mega/vector/trAreas114130/frmt': No such file or
directory

I don't know why there is that message.
Then in Plots Menu, I choose Scatter Plots. I click Add scatter plot. In the Name of imagery group, I choose the raster. In the Name of imagery subgroup, I type the name. I click Create/edit group. In x and y axis, I choose a band. I click Add. Then there is a message that :
"Scatter plot cannot be added.
Multiple of bands ranges <x53.1@PERMANENT:65529 * x53.1@PERMANENT:65529 = 4294049841 > is higher than maximum limit <16810000> "

In this case, I have defined the region. Thank you for your attention.

Regards,
mega

On Wed, Dec 23, 2020 at 4:30 PM Micha Silver <tsvibar@gmail.com <mailto:tsvibar@gmail.com>> wrote:

    On 12/23/20 5:09 AM, mega saputra wrote:
    > Hello guys,
    >
    > I have MODIS data. It is sixteen bit unsigned integer. When I
    want to
    > add scatter plot, there is a warning :
    >
    > "Scatter plot cannot be added.
    > Multiple of bands ranges <x53.1@PERMANENT:65529 *
    > x53.1@PERMANENT:65529 = 4294049841 > is higher than maximum limit
    > <16810000> "

    Please post your computation region (i.e. the output of g.region
    -p) and
    the command you used.

    >
    > So, how to make big maximum limit on Grass GIS, so that I can add
    > scatter plot of MODIS data?
    >
    > Regards,
    > mega
    >
    > _______________________________________________
    > grass-user mailing list
    > grass-user@lists.osgeo.org <mailto:grass-user@lists.osgeo.org>
    > https://lists.osgeo.org/mailman/listinfo/grass-user
    <https://lists.osgeo.org/mailman/listinfo/grass-user&gt;

    -- Micha Silver
    Ben Gurion Univ.
    Sde Boker, Remote Sensing Lab
    cell: +972-523-665918

--
Micha Silver
Ben Gurion Univ.
Sde Boker, Remote Sensing Lab
cell: +972-523-665918

5

The result is the same. It cannot add scatterplot.

Here is the output :

(Sun Dec 27 05:20:25 2020)
g.region -ap raster=rgb.1
projection: 3 (Latitude-Longitude)
zone: 0
datum: wgs84
ellipsoid: wgs84
north: 13:07:42.9N
south: 27:38:15.5S
west: 100:55:15.7E
east: 131:15:38.1E
nsres: 0:00:36.8
ewres: 0:00:36.8
rows: 3988
cols: 2968
cells: 11836384
(Sun Dec 27 05:20:25 2020) Command finished (0 sec)
G__open(read): Unable to open ‘/home/mega/newLocation/PERMAN
ENT/.tmp/mega/vector/trAreas54950/frmt’: No such file or
directory

If you want to try my data, my data is in : https://1drv.ms/u/s!AjXAWS5maHffhyH0q9oTos1umy3D?e=bKa7mX

Thank you.

regards,
mega

On Sat, Dec 26, 2020 at 9:39 PM Micha Silver <tsvibar@gmail.com> wrote:

On 12/25/20 11:46 PM, mega saputra wrote:

Sorry, there is holiday.

I import the raster. Then set the region. This is my syntax :

g.region -up raster=rgb.1

The “-u” flag does not set the current computational region to your
raster. It just reports what would happen if you do reset the region.

Please try:

g.region -ap rast=rgb.1

and report back the output.

projection: 3 (Latitude-Longitude)
zone: 0
datum: wgs84
ellipsoid: wgs84
north: 13:07:42.9N
south: 27:38:15.5S
west: 100:55:15.7E
east: 131:15:38.1E
nsres: 0:00:36.8
ewres: 0:00:36.8
rows: 3988
cols: 2968
cells: 11836384
(Sat Dec 26 05:24:05 2020) Command finished (0 sec)

Then, from menu Imagery > Classify image > Interactive input for
supervised classification. Then the output is like this :

G__open(read): Unable to open ‘/home/mega/newLocation/PERMAN
ENT/.tmp/mega/vector/trAreas114130/frmt’: No such file or
directory

I don’t know why there is that message.
Then in Plots Menu, I choose Scatter Plots. I click Add scatter plot.
In the Name of imagery group, I choose the raster. In the Name of
imagery subgroup, I type the name. I click Create/edit group. In x and
y axis, I choose a band. I click Add. Then there is a message that :
"Scatter plot cannot be added.
Multiple of bands ranges <x53.1@PERMANENT:65529 *
x53.1@PERMANENT:65529 = 4294049841 > is higher than maximum limit
<16810000> "

In this case, I have defined the region. Thank you for your attention.

Regards,
mega

On Wed, Dec 23, 2020 at 4:30 PM Micha Silver <tsvibar@gmail.com
mailto:[tsvibar@gmail.com](mailto:tsvibar@gmail.com)> wrote:

On 12/23/20 5:09 AM, mega saputra wrote:

Hello guys,

I have MODIS data. It is sixteen bit unsigned integer. When I
want to
add scatter plot, there is a warning :

"Scatter plot cannot be added.
Multiple of bands ranges <x53.1@PERMANENT:65529 *
x53.1@PERMANENT:65529 = 4294049841 > is higher than maximum limit
<16810000> "

Please post your computation region (i.e. the output of g.region
-p) and
the command you used.

So, how to make big maximum limit on Grass GIS, so that I can add
scatter plot of MODIS data?

Regards,
mega

grass-user mailing list
grass-user@lists.osgeo.org mailto:[grass-user@lists.osgeo.org](mailto:grass-user@lists.osgeo.org)
https://lists.osgeo.org/mailman/listinfo/grass-user
<https://lists.osgeo.org/mailman/listinfo/grass-user>

–
Micha Silver
Ben Gurion Univ.
Sde Boker, Remote Sensing Lab
cell: +972-523-665918

–
Micha Silver
Ben Gurion Univ.
Sde Boker, Remote Sensing Lab
cell: +972-523-665918

6

Here’s my attempt:

I downloaded the rgb.tif that you linked to, then used it to create a new, fresh location

Create a NEW location “mega” using the georeferenced rgb.tif file to define location

and import the file

Run at the command line, before starting grass:

grass78 -c rgb.tif /home/micha/grass/mega

Now at the grass prompt:

Check projection settings

g.proj -p
-PROJ_INFO-------------------------------------------------
name : WGS 84
datum : wgs84
ellps : wgs84
proj : ll
no_defs : defined
-PROJ_EPSG-------------------------------------------------
epsg : 4326
-PROJ_UNITS------------------------------------------------
unit : degree
units : degrees
meters : 1.0

Install the scatterplot extension

g.extension r.scatterplot
Downloading precompiled GRASS Addons <r.scatterplot>…
Updating extensions metadata file…
Updating extension modules metadata file…
Installation of <r.scatterplot> successfully finished

Now run r.scatterplot on two bands

r.scatterplot input=rgb.3,rgb.2 output=scatter_3_2
100%
Building topology for vector map scatter_3_2@PERMANENT…
Registering primitives…
11830000

Almost 12 million points (!)

···

The X-Y coordinates of the points are the scatterplot values.

HTH,
Micha

On 12/26/2020 11:40 PM, mega saputra wrote:

The result is the same. It cannot add scatterplot.

Here is the output :

(Sun Dec 27 05:20:25 2020)
g.region -ap raster=rgb.1
projection: 3 (Latitude-Longitude)
zone: 0
datum: wgs84
ellipsoid: wgs84
north: 13:07:42.9N
south: 27:38:15.5S
west: 100:55:15.7E
east: 131:15:38.1E
nsres: 0:00:36.8
ewres: 0:00:36.8
rows: 3988
cols: 2968
cells: 11836384
(Sun Dec 27 05:20:25 2020) Command finished (0 sec)
G__open(read): Unable to open ‘/home/mega/newLocation/PERMAN
ENT/.tmp/mega/vector/trAreas54950/frmt’: No such file or
directory

If you want to try my data, my data is in : https://1drv.ms/u/s!AjXAWS5maHffhyH0q9oTos1umy3D?e=bKa7mX

Thank you.

regards,
mega

On Sat, Dec 26, 2020 at 9:39 PM Micha Silver <tsvibar@gmail.com> wrote:

On 12/25/20 11:46 PM, mega saputra wrote:

Sorry, there is holiday.

I import the raster. Then set the region. This is my syntax :

g.region -up raster=rgb.1

The “-u” flag does not set the current computational region to your
raster. It just reports what would happen if you do reset the region.

Please try:

g.region -ap rast=rgb.1

and report back the output.

projection: 3 (Latitude-Longitude)
zone: 0
datum: wgs84
ellipsoid: wgs84
north: 13:07:42.9N
south: 27:38:15.5S
west: 100:55:15.7E
east: 131:15:38.1E
nsres: 0:00:36.8
ewres: 0:00:36.8
rows: 3988
cols: 2968
cells: 11836384
(Sat Dec 26 05:24:05 2020) Command finished (0 sec)

Then, from menu Imagery > Classify image > Interactive input for
supervised classification. Then the output is like this :

G__open(read): Unable to open ‘/home/mega/newLocation/PERMAN
ENT/.tmp/mega/vector/trAreas114130/frmt’: No such file or
directory

I don’t know why there is that message.
Then in Plots Menu, I choose Scatter Plots. I click Add scatter plot.
In the Name of imagery group, I choose the raster. In the Name of
imagery subgroup, I type the name. I click Create/edit group. In x and
y axis, I choose a band. I click Add. Then there is a message that :
"Scatter plot cannot be added.
Multiple of bands ranges <x53.1@PERMANENT:65529 *
x53.1@PERMANENT:65529 = 4294049841 > is higher than maximum limit
<16810000> "

In this case, I have defined the region. Thank you for your attention.

Regards,
mega

On Wed, Dec 23, 2020 at 4:30 PM Micha Silver <tsvibar@gmail.com
mailto:[tsvibar@gmail.com](mailto:tsvibar@gmail.com)> wrote:

On 12/23/20 5:09 AM, mega saputra wrote:

Hello guys,

I have MODIS data. It is sixteen bit unsigned integer. When I
want to
add scatter plot, there is a warning :

"Scatter plot cannot be added.
Multiple of bands ranges <x53.1@PERMANENT:65529 *
x53.1@PERMANENT:65529 = 4294049841 > is higher than maximum limit
<16810000> "

Please post your computation region (i.e. the output of g.region
-p) and
the command you used.

So, how to make big maximum limit on Grass GIS, so that I can add
scatter plot of MODIS data?

Regards,
mega

grass-user mailing list
grass-user@lists.osgeo.org mailto:[grass-user@lists.osgeo.org](mailto:grass-user@lists.osgeo.org)
https://lists.osgeo.org/mailman/listinfo/grass-user
<https://lists.osgeo.org/mailman/listinfo/grass-user>

–
Micha Silver
Ben Gurion Univ.
Sde Boker, Remote Sensing Lab
cell: +972-523-665918

–
Micha Silver
Ben Gurion Univ.
Sde Boker, Remote Sensing Lab
cell: +972-523-665918

-- 
Micha Silver
Remote Sensing Lab, Sde Boker
Ben Gurion University
cell: +972-523-665918
7

Here is the output :
g.proj -p
-PROJ_INFO-------------------------------------------------
name : WGS 84
datum : wgs84
ellps : wgs84
proj : ll
no_defs : defined
-PROJ_EPSG-------------------------------------------------
epsg : 4326
-PROJ_UNITS------------------------------------------------
unit : degree
units : degrees
meters : 1.0
(Mon Dec 28 07:25:33 2020) Command finished (0 sec)
(Mon Dec 28 07:25:42 2020)
g.extension r.scatterplot
Fetching <r.scatterplot> from GRASS GIS Addons repository (be patient)…
Compiling…
Installing…
Updating addons metadata file…
Installation of <r.scatterplot> successfully finished
(Mon Dec 28 07:25:50 2020) Command finished (7 sec)
(Mon Dec 28 07:26:18 2020)
r.scatterplot input=rgb.1,rgb.3 output=scatter_1_3
Building topology for vector map scatter_1_3@PERMANENT…
Registering primitives…
(Mon Dec 28 07:28:40 2020) Command finished (2 min 21 sec)
(Mon Dec 28 07:30:53 2020)
v.info map=scatter_1_3@PERMANENT
±---------------------------------------------------------------------------+
| Name: scatter_1_3 |
| Mapset: PERMANENT |
| Location: newLocation |
| Database: /home/mega |
| Title: |
| Map scale: 1:1 |
| Name of creator: mega |
| Organization: |
| Source date: Mon Dec 28 07:26:19 2020 |

Timestamp (first layer): none
Map format: native
----------------------------------------------------------------------------
Type of map: vector (level: 2)
Number of points: 11836384 Number of centroids: 0
Number of lines: 0 Number of boundaries: 0
Number of areas: 0 Number of islands: 0
Number of faces: 0 Number of kernels: 0
Number of volumes: 0 Number of holes: 0
Map is 3D: Yes
Number of dblinks: 0
Projection: Latitude-Longitude
N: 65528N S: 0
E: 65528E W: 0
B: 0 T: 0
Digitization threshold: 0
Comment:
±---------------------------------------------------------------------------+
(Mon Dec 28 07:30:53 2020) Command finished (0 sec)

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 :

“Unable to zoom to vector map scatter_1_3@PERMANENT
Details: Illegal latitude for North: 65528”

How to solve the problem?

Regards,
mega

On Sun, Dec 27, 2020 at 3:22 PM Micha Silver <tsvibar@gmail.com> wrote:

Here’s my attempt:

I downloaded the rgb.tif that you linked to, then used it to create a new, fresh location

Create a NEW location “mega” using the georeferenced rgb.tif file to define location

and import the file

Run at the command line, before starting grass:

grass78 -c rgb.tif /home/micha/grass/mega

Now at the grass prompt:

Check projection settings

g.proj -p
-PROJ_INFO-------------------------------------------------
name : WGS 84
datum : wgs84
ellps : wgs84
proj : ll
no_defs : defined
-PROJ_EPSG-------------------------------------------------
epsg : 4326
-PROJ_UNITS------------------------------------------------
unit : degree
units : degrees
meters : 1.0

Install the scatterplot extension

g.extension r.scatterplot
Downloading precompiled GRASS Addons <r.scatterplot>…
Updating extensions metadata file…
Updating extension modules metadata file…
Installation of <r.scatterplot> successfully finished

Now run r.scatterplot on two bands

r.scatterplot input=rgb.3,rgb.2 output=scatter_3_2
100%
Building topology for vector map scatter_3_2@PERMANENT…
Registering primitives…
11830000

Almost 12 million points (!)

The X-Y coordinates of the points are the scatterplot values.

HTH,
Micha

On 12/26/2020 11:40 PM, mega saputra wrote:

The result is the same. It cannot add scatterplot.

Here is the output :

(Sun Dec 27 05:20:25 2020)
g.region -ap raster=rgb.1
projection: 3 (Latitude-Longitude)
zone: 0
datum: wgs84
ellipsoid: wgs84
north: 13:07:42.9N
south: 27:38:15.5S
west: 100:55:15.7E
east: 131:15:38.1E
nsres: 0:00:36.8
ewres: 0:00:36.8
rows: 3988
cols: 2968
cells: 11836384
(Sun Dec 27 05:20:25 2020) Command finished (0 sec)
G__open(read): Unable to open ‘/home/mega/newLocation/PERMAN
ENT/.tmp/mega/vector/trAreas54950/frmt’: No such file or
directory

If you want to try my data, my data is in : https://1drv.ms/u/s!AjXAWS5maHffhyH0q9oTos1umy3D?e=bKa7mX

Thank you.

regards,
mega

On Sat, Dec 26, 2020 at 9:39 PM Micha Silver <tsvibar@gmail.com> wrote:

On 12/25/20 11:46 PM, mega saputra wrote:

Sorry, there is holiday.

I import the raster. Then set the region. This is my syntax :

g.region -up raster=rgb.1

The “-u” flag does not set the current computational region to your
raster. It just reports what would happen if you do reset the region.

Please try:

g.region -ap rast=rgb.1

and report back the output.

projection: 3 (Latitude-Longitude)
zone: 0
datum: wgs84
ellipsoid: wgs84
north: 13:07:42.9N
south: 27:38:15.5S
west: 100:55:15.7E
east: 131:15:38.1E
nsres: 0:00:36.8
ewres: 0:00:36.8
rows: 3988
cols: 2968
cells: 11836384
(Sat Dec 26 05:24:05 2020) Command finished (0 sec)

Then, from menu Imagery > Classify image > Interactive input for
supervised classification. Then the output is like this :

G__open(read): Unable to open ‘/home/mega/newLocation/PERMAN
ENT/.tmp/mega/vector/trAreas114130/frmt’: No such file or
directory

I don’t know why there is that message.
Then in Plots Menu, I choose Scatter Plots. I click Add scatter plot.
In the Name of imagery group, I choose the raster. In the Name of
imagery subgroup, I type the name. I click Create/edit group. In x and
y axis, I choose a band. I click Add. Then there is a message that :
"Scatter plot cannot be added.
Multiple of bands ranges <x53.1@PERMANENT:65529 *
x53.1@PERMANENT:65529 = 4294049841 > is higher than maximum limit
<16810000> "

In this case, I have defined the region. Thank you for your attention.

Regards,
mega

On Wed, Dec 23, 2020 at 4:30 PM Micha Silver <tsvibar@gmail.com
mailto:[tsvibar@gmail.com](mailto:tsvibar@gmail.com)> wrote:

On 12/23/20 5:09 AM, mega saputra wrote:

Hello guys,

I have MODIS data. It is sixteen bit unsigned integer. When I
want to
add scatter plot, there is a warning :

"Scatter plot cannot be added.
Multiple of bands ranges <x53.1@PERMANENT:65529 *
x53.1@PERMANENT:65529 = 4294049841 > is higher than maximum limit
<16810000> "

Please post your computation region (i.e. the output of g.region
-p) and
the command you used.

So, how to make big maximum limit on Grass GIS, so that I can add
scatter plot of MODIS data?

Regards,
mega

grass-user mailing list
grass-user@lists.osgeo.org mailto:[grass-user@lists.osgeo.org](mailto:grass-user@lists.osgeo.org)
https://lists.osgeo.org/mailman/listinfo/grass-user
<https://lists.osgeo.org/mailman/listinfo/grass-user>

–
Micha Silver
Ben Gurion Univ.
Sde Boker, Remote Sensing Lab
cell: +972-523-665918

–
Micha Silver
Ben Gurion Univ.
Sde Boker, Remote Sensing Lab
cell: +972-523-665918

-- 
Micha Silver
Remote Sensing Lab, Sde Boker
Ben Gurion University
cell: +972-523-665918
8

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.

But, in this scatter plot, I can not create a polygon like in https://grass.osgeo.org/grass78/manuals/wxGUI.iscatt.html.

I want to create region of interest in the scatter plot link with the raster like in the website page.

Creating scatter plot from the website page will return error. It is about maximum limit.

How to solve this?

Regards,
mega

On Mon, Dec 28, 2020 at 7:38 AM mega saputra <megasaputra8@gmail.com> wrote:

Here is the output :
g.proj -p
-PROJ_INFO-------------------------------------------------
name : WGS 84
datum : wgs84
ellps : wgs84
proj : ll
no_defs : defined
-PROJ_EPSG-------------------------------------------------
epsg : 4326
-PROJ_UNITS------------------------------------------------
unit : degree
units : degrees
meters : 1.0
(Mon Dec 28 07:25:33 2020) Command finished (0 sec)
(Mon Dec 28 07:25:42 2020)
g.extension r.scatterplot
Fetching <r.scatterplot> from GRASS GIS Addons repository (be patient)…
Compiling…
Installing…
Updating addons metadata file…
Installation of <r.scatterplot> successfully finished
(Mon Dec 28 07:25:50 2020) Command finished (7 sec)
(Mon Dec 28 07:26:18 2020)
r.scatterplot input=rgb.1,rgb.3 output=scatter_1_3
Building topology for vector map scatter_1_3@PERMANENT…
Registering primitives…
(Mon Dec 28 07:28:40 2020) Command finished (2 min 21 sec)
(Mon Dec 28 07:30:53 2020)
v.info map=scatter_1_3@PERMANENT
±---------------------------------------------------------------------------+
| Name: scatter_1_3 |
| Mapset: PERMANENT |
| Location: newLocation |
| Database: /home/mega |
| Title: |
| Map scale: 1:1 |
| Name of creator: mega |
| Organization: |
| Source date: Mon Dec 28 07:26:19 2020 |

Timestamp (first layer): none
Map format: native
----------------------------------------------------------------------------
Type of map: vector (level: 2)
Number of points: 11836384 Number of centroids: 0
Number of lines: 0 Number of boundaries: 0
Number of areas: 0 Number of islands: 0
Number of faces: 0 Number of kernels: 0
Number of volumes: 0 Number of holes: 0
Map is 3D: Yes
Number of dblinks: 0
Projection: Latitude-Longitude
N: 65528N S: 0
E: 65528E W: 0
B: 0 T: 0
Digitization threshold: 0
Comment:
±---------------------------------------------------------------------------+
(Mon Dec 28 07:30:53 2020) Command finished (0 sec)

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 :

“Unable to zoom to vector map scatter_1_3@PERMANENT
Details: Illegal latitude for North: 65528”

How to solve the problem?

Regards,
mega

On Sun, Dec 27, 2020 at 3:22 PM Micha Silver <tsvibar@gmail.com> wrote:

Here’s my attempt:

I downloaded the rgb.tif that you linked to, then used it to create a new, fresh location

Create a NEW location “mega” using the georeferenced rgb.tif file to define location

and import the file

Run at the command line, before starting grass:

grass78 -c rgb.tif /home/micha/grass/mega

Now at the grass prompt:

Check projection settings

g.proj -p
-PROJ_INFO-------------------------------------------------
name : WGS 84
datum : wgs84
ellps : wgs84
proj : ll
no_defs : defined
-PROJ_EPSG-------------------------------------------------
epsg : 4326
-PROJ_UNITS------------------------------------------------
unit : degree
units : degrees
meters : 1.0

Install the scatterplot extension

g.extension r.scatterplot
Downloading precompiled GRASS Addons <r.scatterplot>…
Updating extensions metadata file…
Updating extension modules metadata file…
Installation of <r.scatterplot> successfully finished

Now run r.scatterplot on two bands

r.scatterplot input=rgb.3,rgb.2 output=scatter_3_2
100%
Building topology for vector map scatter_3_2@PERMANENT…
Registering primitives…
11830000

Almost 12 million points (!)

The X-Y coordinates of the points are the scatterplot values.

HTH,
Micha

On 12/26/2020 11:40 PM, mega saputra wrote:

The result is the same. It cannot add scatterplot.

Here is the output :

(Sun Dec 27 05:20:25 2020)
g.region -ap raster=rgb.1
projection: 3 (Latitude-Longitude)
zone: 0
datum: wgs84
ellipsoid: wgs84
north: 13:07:42.9N
south: 27:38:15.5S
west: 100:55:15.7E
east: 131:15:38.1E
nsres: 0:00:36.8
ewres: 0:00:36.8
rows: 3988
cols: 2968
cells: 11836384
(Sun Dec 27 05:20:25 2020) Command finished (0 sec)
G__open(read): Unable to open ‘/home/mega/newLocation/PERMAN
ENT/.tmp/mega/vector/trAreas54950/frmt’: No such file or
directory

If you want to try my data, my data is in : https://1drv.ms/u/s!AjXAWS5maHffhyH0q9oTos1umy3D?e=bKa7mX

Thank you.

regards,
mega

On Sat, Dec 26, 2020 at 9:39 PM Micha Silver <tsvibar@gmail.com> wrote:

On 12/25/20 11:46 PM, mega saputra wrote:

Sorry, there is holiday.

I import the raster. Then set the region. This is my syntax :

g.region -up raster=rgb.1

The “-u” flag does not set the current computational region to your
raster. It just reports what would happen if you do reset the region.

Please try:

g.region -ap rast=rgb.1

and report back the output.

projection: 3 (Latitude-Longitude)
zone: 0
datum: wgs84
ellipsoid: wgs84
north: 13:07:42.9N
south: 27:38:15.5S
west: 100:55:15.7E
east: 131:15:38.1E
nsres: 0:00:36.8
ewres: 0:00:36.8
rows: 3988
cols: 2968
cells: 11836384
(Sat Dec 26 05:24:05 2020) Command finished (0 sec)

Then, from menu Imagery > Classify image > Interactive input for
supervised classification. Then the output is like this :

G__open(read): Unable to open ‘/home/mega/newLocation/PERMAN
ENT/.tmp/mega/vector/trAreas114130/frmt’: No such file or
directory

I don’t know why there is that message.
Then in Plots Menu, I choose Scatter Plots. I click Add scatter plot.
In the Name of imagery group, I choose the raster. In the Name of
imagery subgroup, I type the name. I click Create/edit group. In x and
y axis, I choose a band. I click Add. Then there is a message that :
"Scatter plot cannot be added.
Multiple of bands ranges <x53.1@PERMANENT:65529 *
x53.1@PERMANENT:65529 = 4294049841 > is higher than maximum limit
<16810000> "

In this case, I have defined the region. Thank you for your attention.

Regards,
mega

On Wed, Dec 23, 2020 at 4:30 PM Micha Silver <tsvibar@gmail.com
mailto:[tsvibar@gmail.com](mailto:tsvibar@gmail.com)> wrote:

On 12/23/20 5:09 AM, mega saputra wrote:

Hello guys,

I have MODIS data. It is sixteen bit unsigned integer. When I
want to
add scatter plot, there is a warning :

"Scatter plot cannot be added.
Multiple of bands ranges <x53.1@PERMANENT:65529 *
x53.1@PERMANENT:65529 = 4294049841 > is higher than maximum limit
<16810000> "

Please post your computation region (i.e. the output of g.region
-p) and
the command you used.

So, how to make big maximum limit on Grass GIS, so that I can add
scatter plot of MODIS data?

Regards,
mega

grass-user mailing list
grass-user@lists.osgeo.org mailto:[grass-user@lists.osgeo.org](mailto:grass-user@lists.osgeo.org)
https://lists.osgeo.org/mailman/listinfo/grass-user
<https://lists.osgeo.org/mailman/listinfo/grass-user>

–
Micha Silver
Ben Gurion Univ.
Sde Boker, Remote Sensing Lab
cell: +972-523-665918

–
Micha Silver
Ben Gurion Univ.
Sde Boker, Remote Sensing Lab
cell: +972-523-665918

-- 
Micha Silver
Remote Sensing Lab, Sde Boker
Ben Gurion University
cell: +972-523-665918
(attachments)

Screenshot_2020-12-28_08-44-48.png

9
···

On 12/28/2020 1:38 AM, mega saputra wrote:

Here is the output :
g.proj -p

…

…

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 :

“Unable to zoom to vector map scatter_1_3@PERMANENT
Details: Illegal latitude for North: 65528”

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).

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.

How to solve the problem?

Regards,
mega

On Sun, Dec 27, 2020 at 3:22 PM Micha Silver <tsvibar@gmail.com> wrote:

Here’s my attempt:

I downloaded the rgb.tif that you linked to, then used it to create a new, fresh location

Create a NEW location “mega” using the georeferenced rgb.tif file to define location

and import the file

Run at the command line, before starting grass:

grass78 -c rgb.tif /home/micha/grass/mega

Now at the grass prompt:

Check projection settings

g.proj -p
-PROJ_INFO-------------------------------------------------
name : WGS 84
datum : wgs84
ellps : wgs84
proj : ll
no_defs : defined
-PROJ_EPSG-------------------------------------------------
epsg : 4326
-PROJ_UNITS------------------------------------------------
unit : degree
units : degrees
meters : 1.0

Install the scatterplot extension

g.extension r.scatterplot
Downloading precompiled GRASS Addons <r.scatterplot>…
Updating extensions metadata file…
Updating extension modules metadata file…
Installation of <r.scatterplot> successfully finished

Now run r.scatterplot on two bands

r.scatterplot input=rgb.3,rgb.2 output=scatter_3_2
100%
Building topology for vector map scatter_3_2@PERMANENT…
Registering primitives…
11830000

Almost 12 million points (!)

The X-Y coordinates of the points are the scatterplot values.

HTH,
Micha

On 12/26/2020 11:40 PM, mega saputra wrote:

The result is the same. It cannot add scatterplot.

Here is the output :

(Sun Dec 27 05:20:25 2020)
g.region -ap raster=rgb.1
projection: 3 (Latitude-Longitude)
zone: 0
datum: wgs84
ellipsoid: wgs84
north: 13:07:42.9N
south: 27:38:15.5S
west: 100:55:15.7E
east: 131:15:38.1E
nsres: 0:00:36.8
ewres: 0:00:36.8
rows: 3988
cols: 2968
cells: 11836384
(Sun Dec 27 05:20:25 2020) Command finished (0 sec)
G__open(read): Unable to open ‘/home/mega/newLocation/PERMAN
ENT/.tmp/mega/vector/trAreas54950/frmt’: No such file or
directory

If you want to try my data, my data is in : https://1drv.ms/u/s!AjXAWS5maHffhyH0q9oTos1umy3D?e=bKa7mX

Thank you.

regards,
mega

On Sat, Dec 26, 2020 at 9:39 PM Micha Silver <tsvibar@gmail.com> wrote:

On 12/25/20 11:46 PM, mega saputra wrote:

Sorry, there is holiday.

I import the raster. Then set the region. This is my syntax :

g.region -up raster=rgb.1

The “-u” flag does not set the current computational region to your
raster. It just reports what would happen if you do reset the region.

Please try:

g.region -ap rast=rgb.1

and report back the output.

projection: 3 (Latitude-Longitude)
zone: 0
datum: wgs84
ellipsoid: wgs84
north: 13:07:42.9N
south: 27:38:15.5S
west: 100:55:15.7E
east: 131:15:38.1E
nsres: 0:00:36.8
ewres: 0:00:36.8
rows: 3988
cols: 2968
cells: 11836384
(Sat Dec 26 05:24:05 2020) Command finished (0 sec)

Then, from menu Imagery > Classify image > Interactive input for
supervised classification. Then the output is like this :

G__open(read): Unable to open ‘/home/mega/newLocation/PERMAN
ENT/.tmp/mega/vector/trAreas114130/frmt’: No such file or
directory

I don’t know why there is that message.
Then in Plots Menu, I choose Scatter Plots. I click Add scatter plot.
In the Name of imagery group, I choose the raster. In the Name of
imagery subgroup, I type the name. I click Create/edit group. In x and
y axis, I choose a band. I click Add. Then there is a message that :
"Scatter plot cannot be added.
Multiple of bands ranges <x53.1@PERMANENT:65529 *
x53.1@PERMANENT:65529 = 4294049841 > is higher than maximum limit
<16810000> "

In this case, I have defined the region. Thank you for your attention.

Regards,
mega

On Wed, Dec 23, 2020 at 4:30 PM Micha Silver <tsvibar@gmail.com
mailto:[tsvibar@gmail.com](mailto:tsvibar@gmail.com)> wrote:

On 12/23/20 5:09 AM, mega saputra wrote:

Hello guys,

I have MODIS data. It is sixteen bit unsigned integer. When I
want to
add scatter plot, there is a warning :

"Scatter plot cannot be added.
Multiple of bands ranges <x53.1@PERMANENT:65529 *
x53.1@PERMANENT:65529 = 4294049841 > is higher than maximum limit
<16810000> "

Please post your computation region (i.e. the output of g.region
-p) and
the command you used.

So, how to make big maximum limit on Grass GIS, so that I can add
scatter plot of MODIS data?

Regards,
mega

grass-user mailing list
grass-user@lists.osgeo.org mailto:[grass-user@lists.osgeo.org](mailto:grass-user@lists.osgeo.org)
https://lists.osgeo.org/mailman/listinfo/grass-user
<https://lists.osgeo.org/mailman/listinfo/grass-user>

–
Micha Silver
Ben Gurion Univ.
Sde Boker, Remote Sensing Lab
cell: +972-523-665918

–
Micha Silver
Ben Gurion Univ.
Sde Boker, Remote Sensing Lab
cell: +972-523-665918

-- 
Micha Silver
Remote Sensing Lab, Sde Boker
Ben Gurion University
cell: +972-523-665918

-- 
Micha Silver
Remote Sensing Lab, Sde Boker
Ben Gurion University
cell: +972-523-665918
10

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.

I will post new email in new title as your suggestion about interactive scatter plot.

Thank you for your response.

Regards,

mega

On Mon, Dec 28, 2020 at 3:28 PM Micha Silver <tsvibar@gmail.com> wrote:

On 12/28/2020 1:38 AM, mega saputra wrote:

Here is the output :
g.proj -p

…

…

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 :

“Unable to zoom to vector map scatter_1_3@PERMANENT
Details: Illegal latitude for North: 65528”

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).

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.

How to solve the problem?

Regards,
mega

On Sun, Dec 27, 2020 at 3:22 PM Micha Silver <tsvibar@gmail.com> wrote:

Here’s my attempt:

I downloaded the rgb.tif that you linked to, then used it to create a new, fresh location

Create a NEW location “mega” using the georeferenced rgb.tif file to define location

and import the file

Run at the command line, before starting grass:

grass78 -c rgb.tif /home/micha/grass/mega

Now at the grass prompt:

Check projection settings

g.proj -p
-PROJ_INFO-------------------------------------------------
name : WGS 84
datum : wgs84
ellps : wgs84
proj : ll
no_defs : defined
-PROJ_EPSG-------------------------------------------------
epsg : 4326
-PROJ_UNITS------------------------------------------------
unit : degree
units : degrees
meters : 1.0

Install the scatterplot extension

g.extension r.scatterplot
Downloading precompiled GRASS Addons <r.scatterplot>…
Updating extensions metadata file…
Updating extension modules metadata file…
Installation of <r.scatterplot> successfully finished

Now run r.scatterplot on two bands

r.scatterplot input=rgb.3,rgb.2 output=scatter_3_2
100%
Building topology for vector map scatter_3_2@PERMANENT…
Registering primitives…
11830000

Almost 12 million points (!)

The X-Y coordinates of the points are the scatterplot values.

HTH,
Micha

On 12/26/2020 11:40 PM, mega saputra wrote:

The result is the same. It cannot add scatterplot.

Here is the output :

(Sun Dec 27 05:20:25 2020)
g.region -ap raster=rgb.1
projection: 3 (Latitude-Longitude)
zone: 0
datum: wgs84
ellipsoid: wgs84
north: 13:07:42.9N
south: 27:38:15.5S
west: 100:55:15.7E
east: 131:15:38.1E
nsres: 0:00:36.8
ewres: 0:00:36.8
rows: 3988
cols: 2968
cells: 11836384
(Sun Dec 27 05:20:25 2020) Command finished (0 sec)
G__open(read): Unable to open ‘/home/mega/newLocation/PERMAN
ENT/.tmp/mega/vector/trAreas54950/frmt’: No such file or
directory

If you want to try my data, my data is in : https://1drv.ms/u/s!AjXAWS5maHffhyH0q9oTos1umy3D?e=bKa7mX

Thank you.

regards,
mega

On Sat, Dec 26, 2020 at 9:39 PM Micha Silver <tsvibar@gmail.com> wrote:

On 12/25/20 11:46 PM, mega saputra wrote:

Sorry, there is holiday.

I import the raster. Then set the region. This is my syntax :

g.region -up raster=rgb.1

The “-u” flag does not set the current computational region to your
raster. It just reports what would happen if you do reset the region.

Please try:

g.region -ap rast=rgb.1

and report back the output.

projection: 3 (Latitude-Longitude)
zone: 0
datum: wgs84
ellipsoid: wgs84
north: 13:07:42.9N
south: 27:38:15.5S
west: 100:55:15.7E
east: 131:15:38.1E
nsres: 0:00:36.8
ewres: 0:00:36.8
rows: 3988
cols: 2968
cells: 11836384
(Sat Dec 26 05:24:05 2020) Command finished (0 sec)

Then, from menu Imagery > Classify image > Interactive input for
supervised classification. Then the output is like this :

G__open(read): Unable to open ‘/home/mega/newLocation/PERMAN
ENT/.tmp/mega/vector/trAreas114130/frmt’: No such file or
directory

I don’t know why there is that message.
Then in Plots Menu, I choose Scatter Plots. I click Add scatter plot.
In the Name of imagery group, I choose the raster. In the Name of
imagery subgroup, I type the name. I click Create/edit group. In x and
y axis, I choose a band. I click Add. Then there is a message that :
"Scatter plot cannot be added.
Multiple of bands ranges <x53.1@PERMANENT:65529 *
x53.1@PERMANENT:65529 = 4294049841 > is higher than maximum limit
<16810000> "

In this case, I have defined the region. Thank you for your attention.

Regards,
mega

On Wed, Dec 23, 2020 at 4:30 PM Micha Silver <tsvibar@gmail.com
mailto:[tsvibar@gmail.com](mailto:tsvibar@gmail.com)> wrote:

On 12/23/20 5:09 AM, mega saputra wrote:

Hello guys,

I have MODIS data. It is sixteen bit unsigned integer. When I
want to
add scatter plot, there is a warning :

"Scatter plot cannot be added.
Multiple of bands ranges <x53.1@PERMANENT:65529 *
x53.1@PERMANENT:65529 = 4294049841 > is higher than maximum limit
<16810000> "

Please post your computation region (i.e. the output of g.region
-p) and
the command you used.

So, how to make big maximum limit on Grass GIS, so that I can add
scatter plot of MODIS data?

Regards,
mega

grass-user mailing list
grass-user@lists.osgeo.org mailto:[grass-user@lists.osgeo.org](mailto:grass-user@lists.osgeo.org)
https://lists.osgeo.org/mailman/listinfo/grass-user
<https://lists.osgeo.org/mailman/listinfo/grass-user>

–
Micha Silver
Ben Gurion Univ.
Sde Boker, Remote Sensing Lab
cell: +972-523-665918

–
Micha Silver
Ben Gurion Univ.
Sde Boker, Remote Sensing Lab
cell: +972-523-665918

-- 
Micha Silver
Remote Sensing Lab, Sde Boker
Ben Gurion University
cell: +972-523-665918

-- 
Micha Silver
Remote Sensing Lab, Sde Boker
Ben Gurion University
cell: +972-523-665918