[GRASS-user] Scatter plot cannot be added
mega saputra
megasaputra8 at gmail.com
Sun Dec 27 15:38:30 PST 2020
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/dd5f061a/attachment-0001.html>
More information about the grass-user
mailing list