|
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
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
Thank you.
regards,
mega
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 <tsvi...@gmail.com
> <mailto:tsvi...@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>
>
> --
> 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
|
