[GRASS-user] Scatter plot cannot be added
mega saputra
megasaputra8 at gmail.com
Sun Dec 27 16:52:18 PST 2020
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 at 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 at 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 at 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 at 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 at 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 at 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 at 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 at PERMANENT:65529 *
>>> > x53.1 at 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 at gmail.com
>>> > <mailto:tsvibar at 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 at PERMANENT:65529 *
>>> > > x53.1 at 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 at lists.osgeo.org <mailto:grass-user at 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
>>
>>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/grass-user/attachments/20201228/d9d94b7a/attachment-0001.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: Screenshot_2020-12-28_08-44-48.png
Type: image/png
Size: 31575 bytes
Desc: not available
URL: <http://lists.osgeo.org/pipermail/grass-user/attachments/20201228/d9d94b7a/attachment-0001.png>
More information about the grass-user
mailing list