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 <tsvi...@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@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 <tsvi...@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@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 > >
_______________________________________________ grass-user mailing list grass-user@lists.osgeo.org https://lists.osgeo.org/mailman/listinfo/grass-user