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 <[email protected]> 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 <[email protected]> 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 <[email protected]> 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 <[email protected] >>> > <mailto:[email protected]>> 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 >>> > > [email protected] <mailto:[email protected]> >>> > > 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 >> >>
_______________________________________________ grass-user mailing list [email protected] https://lists.osgeo.org/mailman/listinfo/grass-user
