[GRASS-user] question on i.nightlights.intercalibration

[Gabriel Cotlier]

Dear Markus I fowled your instructions and worked :

- first start GRASS in the PERMANENT mapset, then run
'g.region -p n=75:00:15N s=65:00:15S w=179:59:45W e=180:00:15E res=00:00:30
-s'

- close grass,

open grass again in my mapsert where my data is and run the command:
'g.region -p'
 and and it worked !!!!


[image: image.png]

Thanks a lot!



On Wed, Aug 28, 2019 at 8:46 AM Markus Metz <markus.metz.giswork at gmail.com>
wrote:

>
>
> On Wed, Aug 28, 2019 at 12:54 PM Gabriel Cotlier <gabiklm01 at gmail.com>
> wrote:
> >
> > Dear Markus,
> > Thanks a lot for your replay,
> > I have tried the command to change default region with g.region -p
> n=75:00:15N s=65:00:15S w=179:59:45W e=180:00:15E res=00:00:30 -s
> > and did not allow me because according to the error message bellow my
> data is in a mapset that is not <PERMANENT> . I have two mapsets one called
> "gabriel" where all the data is and another <PERMANENT>....
> > How can I solve this problem?
>
> start GRASS in the PERMANENT mapset, then run
> g.region -p n=75:00:15N s=65:00:15S w=179:59:45W e=180:00:15E res=00:00:30
> -s
>
> after that start GRASS again in the mapset with your data.
>
> Markus M
>
> > Thanks a lot again for your help.
> > Regards,
> > Gabriel
> >
> >
> >
> >
> >
> >
> >
> >
> > On Tue, Aug 27, 2019 at 4:02 PM Markus Metz <
> markus.metz.giswork at gmail.com> wrote:
> >>
> >>
> >>
> >> On Tue, Aug 27, 2019 at 2:25 PM Gabriel Cotlier <gabiklm01 at gmail.com>
> wrote:
> >> >
> >> > Dear Markus,
> >> >
> >> > here is my log of the command, "g.region -d -u -p":
> >> >
> >> >
> >> > (Tue Aug 27 09:23:26 2019)
>
> >> > g.region -d -u -p
>
> >> > 360 degree EW extent is exceeded by 1.99983 cells
> >> > 360 degree EW extent is exceeded by 1.99983 cells
> >> > projection: 3 (Latitude-Longitude)
> >> > zone:       0
> >> > datum:      wgs84
> >> > ellipsoid:  wgs84
> >> > north:      75:00:29.99892N
> >> > south:      65:00:29.999064S
> >> > west:       180:00:29.997408W
> >> > east:       180:00:29.997408E
> >> > nsres:      0:00:30
> >> > ewres:      0:00:30
> >> > rows:       16802
> >> > cols:       43202
> >> > cells:      725880004
> >> > (Tue Aug 27 09:23:26 2019) Command finished (0 sec)
>
> >>
> >> if you want to avoid this message, change our default region with
> >>
> >> g.region -p n=75:00:15N s=65:00:15S w=179:59:45W e=180:00:15E
> res=00:00:30 -s
> >>
> >> QGIS uses GDAL to fetch metadata of raster maps, therefore the output
> of gdalinfo would be more informative, but you know the relevant parts
> already: these raster maps need to be imported with r.in.gdal -a
> >>
> >> Markus M
> >>
> >> >
> >> >
> >> >
> >> > On Mon, Aug 26, 2019 at 4:01 PM Markus Metz <
> markus.metz.giswork at gmail.com> wrote:
> >> >>
> >> >>
> >> >>
> >> >> On Mon, Aug 26, 2019 at 1:14 AM Gabriel Cotlier <gabiklm01 at gmail.com>
> wrote:
> >> >> >
> >> >> > Dear Markus,
> >> >> >
> >> >> > I uninstall and reinstall grass 7.6.1 64 bits and still my log for
> the commands
> >> >> >
> >> >> > g.region -p
> >> >> >
> >> >> > and
> >> >> >
> >> >> > g.region -p n=75:00:15N s=65:00:15S w=179:59:45W e=180:00:15E
> res=00:00:30
> >> >> >
> >> >> > are respectively :
> >> >> >
> >> >> > (Sun Aug 25 20:01:47 2019)
>
> >> >> > g.region -p
>
> >> >> > 360 degree EW extent is exceeded by 1.99983 cells
> >> >> > projection: 3 (Latitude-Longitude)
> >> >> > zone:       0
> >> >> > datum:      wgs84
> >> >> > ellipsoid:  wgs84
> >> >> > north:      75:00:15N
> >> >> > south:      65:00:15S
> >> >> > west:       179:59:45W
> >> >> > east:       180:00:15E
> >> >> > nsres:      0:00:30
> >> >> > ewres:      0:00:30
> >> >> > rows:       16801
> >> >> > cols:       43200
> >> >> > cells:      725803200
> >> >> > (Sun Aug 25 20:01:47 2019) Command finished (0 sec)
>
> >> >> > (Sun Aug 25 20:02:10 2019)
> >> >> >
> >> >> > g.region -p n=75:00:15N s=65:00:15S w=179:59:45W e=180:00:15E
> res=00:00:30
> >> >> > 360 degree EW extent is exceeded by 1.99983 cells
> >> >>
> >> >> Strange, the region defined with "n=75:00:15N s=65:00:15S
> w=179:59:45W e=180:00:15E res=00:00:30" should not cause a message like
> >> >> "360 degree EW extent is exceeded by 1.99983 cells"
> >> >>
> >> >> Maybe the default region has strange settings. What is the output of
> "g.region -d -u -p" ?
> >> >>
> >> >> Markus M
> >> >>
> >> >> > projection: 3 (Latitude-Longitude)
> >> >> > zone:       0
> >> >> > datum:      wgs84
> >> >> > ellipsoid:  wgs84
> >> >> > north:      75:00:15N
> >> >> > south:      65:00:15S
> >> >> > west:       179:59:45W
> >> >> > east:       180:00:15E
> >> >> > nsres:      0:00:30
> >> >> > ewres:      0:00:30
> >> >> > rows:       16801
> >> >> > cols:       43200
> >> >> > cells:      725803200
> >> >> > (Sun Aug 25 20:02:10 2019) Command finished (0 sec)
> >> >> >
> >> >> > Are the resulting intercallibrated layers ready for further use
> outside grass exported as GeoTIFF or I need to get inside grass without any
> messages that 360 degree EW extent is exceeded for the layer to be OK?
> >> >> >
> >> >> > further more I would like to ask if maybe the max and min for
> intercallibtated layer are somewhere already available to check whether the
> max and min I got are correct?
> >> >> > Thanks a lot for your help.
> >> >> > Regards,
> >> >> > Gabriel
> >> >> >
> >> >> >
> >> >> >
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/grass-user/attachments/20190901/0aaf5bf4/attachment-0001.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: image.png
Type: image/png
Size: 21207 bytes
Desc: not available
URL: <http://lists.osgeo.org/pipermail/grass-user/attachments/20190901/0aaf5bf4/attachment-0001.png>