<html style="direction: ltr;">
  <head>
    <meta http-equiv="Content-Type" content="text/html; charset=UTF-8">
    <style id="bidiui-paragraph-margins" type="text/css">body p { margin-bottom: 0cm; margin-top: 0pt; } </style>
  </head>
  <body bidimailui-charset-is-forced="true" style="direction: ltr;">
    <p><br>
    </p>
    <div class="moz-cite-prefix">On 12/28/2020 1:38 AM, mega saputra
      wrote:<br>
    </div>
    <blockquote type="cite"
cite="mid:CAPCs1vcpytN73_GwFtA1LOEf8DCDUL0JqzY0XcaZx+UPZ8vd3w@mail.gmail.com">
      <meta http-equiv="content-type" content="text/html; charset=UTF-8">
      <div dir="ltr">
        <div>Here is the output :</div>
        <div>g.proj -p                                                  
                              </div>
      </div>
    </blockquote>
    <p><br>
    </p>
    <p>.....</p>
    <p>....<br>
    </p>
    <blockquote type="cite"
cite="mid:CAPCs1vcpytN73_GwFtA1LOEf8DCDUL0JqzY0XcaZx+UPZ8vd3w@mail.gmail.com">
      <div dir="ltr">
        <div>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 : <br>
        </div>
        <div>"Unable to zoom to vector map <scatter_1_3@PERMANENT></div>
        <div>Details: Illegal latitude for North: 65528"</div>
        <div><br>
        </div>
      </div>
    </blockquote>
    <p>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).</p>
    <p><br>
    </p>
    <p>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.</p>
    <p> <br>
    </p>
    <blockquote type="cite"
cite="mid:CAPCs1vcpytN73_GwFtA1LOEf8DCDUL0JqzY0XcaZx+UPZ8vd3w@mail.gmail.com">
      <div dir="ltr">
        <div>How to solve the problem?<br>
        </div>
        <div><br>
        </div>
        <div>Regards,</div>
        <div>mega<br>
        </div>
      </div>
      <br>
      <div class="gmail_quote">
        <div dir="ltr" class="gmail_attr">On Sun, Dec 27, 2020 at 3:22
          PM Micha Silver <<a href="mailto:tsvibar@gmail.com"
            moz-do-not-send="true">tsvibar@gmail.com</a>> wrote:<br>
        </div>
        <blockquote class="gmail_quote" style="margin:0px 0px 0px
          0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
          <div style="direction:ltr">
            <p>Here's my attempt:</p>
            <p>I downloaded the rgb.tif that you linked to, then used it
              to create a new, fresh location</p>
            <p><br>
            </p>
            <p># Create a NEW location "mega" using the georeferenced
              rgb.tif file to define location<br>
              # and import the file</p>
            <p># Run at the command line, *before* starting grass:<br>
            </p>
            <p>grass78 -c rgb.tif /home/micha/grass/mega</p>
            <p><br>
            </p>
            <p># Now at the grass prompt:</p>
            <p># Check projection settings<br>
              g.proj -p<br>
-PROJ_INFO-------------------------------------------------<br>
              name       : WGS 84<br>
              datum      : wgs84<br>
              ellps      : wgs84<br>
              proj       : ll<br>
              no_defs    : defined<br>
-PROJ_EPSG-------------------------------------------------<br>
              epsg       : 4326<br>
-PROJ_UNITS------------------------------------------------<br>
              unit       : degree<br>
              units      : degrees<br>
              meters     : 1.0<br>
              <br>
              # Install the scatterplot extension<br>
              g.extension r.scatterplot<br>
              Downloading precompiled GRASS Addons
              <r.scatterplot>...<br>
              Updating extensions metadata file...<br>
              Updating extension modules metadata file...<br>
              Installation of <r.scatterplot> successfully
              finished<br>
              <br>
              # Now run r.scatterplot on two bands<br>
              r.scatterplot input=rgb.3,rgb.2 output=scatter_3_2<br>
               100%<br>
              Building topology for vector map
              <scatter_3_2@PERMANENT>...<br>
              Registering primitives...<br>
                11830000<br>
              # Almost 12 million points (!)<br>
            </p>
            <div><br>
            </div>
            <div><br>
            </div>
            <div>The X-Y coordinates of the points are the scatterplot
              values.</div>
            <div><br>
            </div>
            <div>HTH,</div>
            <div>Micha<br>
            </div>
            <div><br>
            </div>
            <div>On 12/26/2020 11:40 PM, mega saputra wrote:<br>
            </div>
            <blockquote type="cite">
              <div dir="ltr">
                <div>The result is the same. It cannot add scatterplot.<br>
                </div>
                <div>Here is the output :</div>
                <div><br>
                </div>
                <div>(Sun Dec 27 05:20:25 2020)                        
                                               <br>
                  g.region -ap raster=rgb.1                            
                                            <br>
                  projection: 3 (Latitude-Longitude)<br>
                  zone:       0<br>
                  datum:      wgs84<br>
                  ellipsoid:  wgs84<br>
                  north:      13:07:42.9N<br>
                  south:      27:38:15.5S<br>
                  west:       100:55:15.7E<br>
                  east:       131:15:38.1E<br>
                  nsres:      0:00:36.8<br>
                  ewres:      0:00:36.8<br>
                  rows:       3988<br>
                  cols:       2968<br>
                  cells:      11836384<br>
                  (Sun Dec 27 05:20:25 2020) Command finished (0 sec)  
                                            <br>
                  G__open(read): Unable to open
                  '/home/mega/newLocation/PERMAN<br>
                  ENT/.tmp/mega/vector/trAreas54950/frmt': No such file
                  or<br>
                  directory</div>
                <div><br>
                </div>
                <div>If you want to try my data, my data is in : <a
                    href="https://1drv.ms/u/s!AjXAWS5maHffhyH0q9oTos1umy3D?e=bKa7mX"
                    target="_blank" moz-do-not-send="true">https://1drv.ms/u/s!AjXAWS5maHffhyH0q9oTos1umy3D?e=bKa7mX</a></div>
                <div><br>
                </div>
                <div>Thank you.</div>
                <div><br>
                </div>
                <div>regards,</div>
                <div>mega<br>
                </div>
                <div><br>
                </div>
              </div>
              <br>
              <div class="gmail_quote">
                <div dir="ltr" class="gmail_attr">On Sat, Dec 26, 2020
                  at 9:39 PM Micha Silver <<a
                    href="mailto:tsvibar@gmail.com" target="_blank"
                    moz-do-not-send="true">tsvibar@gmail.com</a>>
                  wrote:<br>
                </div>
                <blockquote class="gmail_quote" style="margin:0px 0px
                  0px 0.8ex;border-left:1px solid
                  rgb(204,204,204);padding-left:1ex"><br>
                  On 12/25/20 11:46 PM, mega saputra wrote:<br>
                  > Sorry, there is holiday.<br>
                  ><br>
                  > I import the raster. Then set the region. This is
                  my syntax :<br>
                  ><br>
                  > g.region -up raster=rgb.1<br>
                  <br>
                  <br>
                  The "-u" flag does not set the current computational
                  region to your <br>
                  raster. It just reports what would happen if you do
                  reset the region.<br>
                  <br>
                  Please try:<br>
                  <br>
                  <br>
                  g.region -ap rast=rgb.1<br>
                  <br>
                  <br>
                  and report back the output.<br>
                  <br>
                  <br>
                  <br>
                  > projection: 3 (Latitude-Longitude)<br>
                  > zone:       0<br>
                  > datum:      wgs84<br>
                  > ellipsoid:  wgs84<br>
                  > north:      13:07:42.9N<br>
                  > south:      27:38:15.5S<br>
                  > west:       100:55:15.7E<br>
                  > east:       131:15:38.1E<br>
                  > nsres:      0:00:36.8<br>
                  > ewres:      0:00:36.8<br>
                  > rows:       3988<br>
                  > cols:       2968<br>
                  > cells:      11836384<br>
                  > (Sat Dec 26 05:24:05 2020) Command finished (0
                  sec)<br>
                  ><br>
                  > Then, from menu Imagery > Classify image >
                  Interactive input for <br>
                  > supervised classification. Then the output is
                  like this :<br>
                  ><br>
                  > G__open(read): Unable to open
                  '/home/mega/newLocation/PERMAN<br>
                  > ENT/.tmp/mega/vector/trAreas114130/frmt': No such
                  file or<br>
                  > directory<br>
                  ><br>
                  > I don't know why there is that message.<br>
                  > Then in Plots Menu, I choose Scatter Plots. I
                  click Add scatter plot. <br>
                  > In the Name of imagery group, I choose the
                  raster. In the Name of <br>
                  > imagery subgroup, I type the name. I click
                  Create/edit group. In x and <br>
                  > y axis, I choose a band. I click Add. Then there
                  is a message that :<br>
                  > "Scatter plot cannot be added.<br>
                  > Multiple of bands ranges
                  <x53.1@PERMANENT:65529 * <br>
                  > x53.1@PERMANENT:65529 = 4294049841 > is higher
                  than maximum limit <br>
                  > <16810000> "<br>
                  ><br>
                  > In this case, I have defined the region. Thank
                  you for your attention.<br>
                  ><br>
                  > Regards,<br>
                  > mega<br>
                  ><br>
                  > On Wed, Dec 23, 2020 at 4:30 PM Micha Silver <<a
                    href="mailto:tsvibar@gmail.com" target="_blank"
                    moz-do-not-send="true">tsvibar@gmail.com</a> <br>
                  > <mailto:<a href="mailto:tsvibar@gmail.com"
                    target="_blank" moz-do-not-send="true">tsvibar@gmail.com</a>>>
                  wrote:<br>
                  ><br>
                  ><br>
                  >     On 12/23/20 5:09 AM, mega saputra wrote:<br>
                  >     > Hello guys,<br>
                  >     ><br>
                  >     > I have MODIS data. It is sixteen bit
                  unsigned integer. When I<br>
                  >     want to<br>
                  >     > add scatter plot, there is a warning :<br>
                  >     ><br>
                  >     > "Scatter plot cannot be added.<br>
                  >     > Multiple of bands ranges
                  <x53.1@PERMANENT:65529 *<br>
                  >     > x53.1@PERMANENT:65529 = 4294049841 >
                  is higher than maximum limit<br>
                  >     > <16810000> "<br>
                  ><br>
                  ><br>
                  >     Please post your computation region (i.e. the
                  output of g.region<br>
                  >     -p) and<br>
                  >     the command you used.<br>
                  ><br>
                  ><br>
                  >     ><br>
                  >     > So, how to make big maximum limit on
                  Grass GIS, so that I can add<br>
                  >     > scatter plot of MODIS data?<br>
                  >     ><br>
                  >     > Regards,<br>
                  >     > mega<br>
                  >     ><br>
                  >     >
                  _______________________________________________<br>
                  >     > grass-user mailing list<br>
                  >     > <a
                    href="mailto:grass-user@lists.osgeo.org"
                    target="_blank" moz-do-not-send="true">grass-user@lists.osgeo.org</a>
                  <mailto:<a href="mailto:grass-user@lists.osgeo.org"
                    target="_blank" moz-do-not-send="true">grass-user@lists.osgeo.org</a>><br>
                  >     > <a
                    href="https://lists.osgeo.org/mailman/listinfo/grass-user"
                    rel="noreferrer" target="_blank"
                    moz-do-not-send="true">https://lists.osgeo.org/mailman/listinfo/grass-user</a><br>
                  >     <<a
                    href="https://lists.osgeo.org/mailman/listinfo/grass-user"
                    rel="noreferrer" target="_blank"
                    moz-do-not-send="true">https://lists.osgeo.org/mailman/listinfo/grass-user</a>><br>
                  ><br>
                  >     -- <br>
                  >     Micha Silver<br>
                  >     Ben Gurion Univ.<br>
                  >     Sde Boker, Remote Sensing Lab<br>
                  >     cell: +972-523-665918<br>
                  ><br>
                  -- <br>
                  Micha Silver<br>
                  Ben Gurion Univ.<br>
                  Sde Boker, Remote Sensing Lab<br>
                  cell: +972-523-665918<br>
                  <br>
                </blockquote>
              </div>
            </blockquote>
            <pre cols="72">-- 
Micha Silver
Remote Sensing Lab, Sde Boker
Ben Gurion University
cell: +972-523-665918</pre>
          </div>
        </blockquote>
      </div>
    </blockquote>
    <pre class="moz-signature" cols="72">-- 
Micha Silver
Remote Sensing Lab, Sde Boker
Ben Gurion University
cell: +972-523-665918</pre>
  </body>
</html>