[GRASS-user] Scatter plot cannot be added
mega saputra
megasaputra8 at gmail.com
Mon Dec 28 15:31:58 PST 2020
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 at 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 at 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 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
>>
>> --
> 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/20201229/f045fb40/attachment.html>
More information about the grass-user
mailing list