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 <tsvi...@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@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 <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 >> >> -- > 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