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

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
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

On Mon, Aug 26, 2019 at 1:14 AM Gabriel Cotlier <gabiklm01@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

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)

On Mon, Aug 26, 2019 at 4:01 PM Markus Metz <markus.metz.giswork@gmail.com> wrote:

On Mon, Aug 26, 2019 at 1:14 AM Gabriel Cotlier <gabiklm01@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

What should I do in this case?

On Tue, Aug 27, 2019 at 9:24 AM Gabriel Cotlier <gabiklm01@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)

On Mon, Aug 26, 2019 at 4:01 PM Markus Metz <markus.metz.giswork@gmail.com> wrote:

On Mon, Aug 26, 2019 at 1:14 AM Gabriel Cotlier <gabiklm01@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

while for individual raster maps I get:

image.png

On Tue, Aug 27, 2019 at 9:25 AM Gabriel Cotlier <gabiklm01@gmail.com> wrote:

What should I do in this case?

On Tue, Aug 27, 2019 at 9:24 AM Gabriel Cotlier <gabiklm01@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)

On Mon, Aug 26, 2019 at 4:01 PM Markus Metz <markus.metz.giswork@gmail.com> wrote:

On Mon, Aug 26, 2019 at 1:14 AM Gabriel Cotlier <gabiklm01@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

imported into QGIS layer visible and the of the layer properties are:

image.png

On Tue, Aug 27, 2019 at 9:27 AM Gabriel Cotlier <gabiklm01@gmail.com> wrote:

while for individual raster maps I get:

image.png

On Tue, Aug 27, 2019 at 9:25 AM Gabriel Cotlier <gabiklm01@gmail.com> wrote:

What should I do in this case?

On Tue, Aug 27, 2019 at 9:24 AM Gabriel Cotlier <gabiklm01@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)

On Mon, Aug 26, 2019 at 4:01 PM Markus Metz <markus.metz.giswork@gmail.com> wrote:

On Mon, Aug 26, 2019 at 1:14 AM Gabriel Cotlier <gabiklm01@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

On Tue, Aug 27, 2019 at 2:25 PM Gabriel Cotlier <gabiklm01@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@gmail.com> wrote:

On Mon, Aug 26, 2019 at 1:14 AM Gabriel Cotlier <gabiklm01@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

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 . I have two mapsets one called “gabriel” where all the data is and another …

How can I solve this problem?
Thanks a lot again for your help.
Regards,
Gabriel

image.png

On Tue, Aug 27, 2019 at 4:02 PM Markus Metz <markus.metz.giswork@gmail.com> wrote:

On Tue, Aug 27, 2019 at 2:25 PM Gabriel Cotlier <gabiklm01@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@gmail.com> wrote:

On Mon, Aug 26, 2019 at 1:14 AM Gabriel Cotlier <gabiklm01@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

On Wed, Aug 28, 2019 at 12:54 PM Gabriel Cotlier <gabiklm01@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 . I have two mapsets one called “gabriel” where all the data is and another …

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@gmail.com> wrote:

On Tue, Aug 27, 2019 at 2:25 PM Gabriel Cotlier <gabiklm01@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@gmail.com> wrote:

On Mon, Aug 26, 2019 at 1:14 AM Gabriel Cotlier <gabiklm01@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

Dear Markus I fowled your instructions and worked :

  • firststart 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.png

Thanks a lot!

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

On Wed, Aug 28, 2019 at 12:54 PM Gabriel Cotlier <gabiklm01@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 . I have two mapsets one called “gabriel” where all the data is and another …

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@gmail.com> wrote:

On Tue, Aug 27, 2019 at 2:25 PM Gabriel Cotlier <gabiklm01@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@gmail.com> wrote:

On Mon, Aug 26, 2019 at 1:14 AM Gabriel Cotlier <gabiklm01@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