question on i.nightlights.intercalibration

classic Classic list List threaded Threaded
10 messages Options
Reply | Threaded
Open this post in threaded view
|

question on i.nightlights.intercalibration

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

                


_______________________________________________
grass-user mailing list
[hidden email]
https://lists.osgeo.org/mailman/listinfo/grass-user
Reply | Threaded
Open this post in threaded view
|

Re: question on i.nightlights.intercalibration

Markus Metz-3


On Mon, Aug 26, 2019 at 1:14 AM Gabriel Cotlier <[hidden email]> 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
>
>                
>

_______________________________________________
grass-user mailing list
[hidden email]
https://lists.osgeo.org/mailman/listinfo/grass-user
Reply | Threaded
Open this post in threaded view
|

Re: question on i.nightlights.intercalibration

Gabriel Cotlier
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 <[hidden email]> wrote:


On Mon, Aug 26, 2019 at 1:14 AM Gabriel Cotlier <[hidden email]> 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
>
>                
>

_______________________________________________
grass-user mailing list
[hidden email]
https://lists.osgeo.org/mailman/listinfo/grass-user
Reply | Threaded
Open this post in threaded view
|

Re: question on i.nightlights.intercalibration

Gabriel Cotlier
What should I do in this case?

On Tue, Aug 27, 2019 at 9:24 AM Gabriel Cotlier <[hidden email]> 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 <[hidden email]> wrote:


On Mon, Aug 26, 2019 at 1:14 AM Gabriel Cotlier <[hidden email]> 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
>
>                
>

_______________________________________________
grass-user mailing list
[hidden email]
https://lists.osgeo.org/mailman/listinfo/grass-user
Reply | Threaded
Open this post in threaded view
|

Re: question on i.nightlights.intercalibration

Gabriel Cotlier
while for individual raster maps I get:

image.png

On Tue, Aug 27, 2019 at 9:25 AM Gabriel Cotlier <[hidden email]> wrote:
What should I do in this case?

On Tue, Aug 27, 2019 at 9:24 AM Gabriel Cotlier <[hidden email]> 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 <[hidden email]> wrote:


On Mon, Aug 26, 2019 at 1:14 AM Gabriel Cotlier <[hidden email]> 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
>
>                
>

_______________________________________________
grass-user mailing list
[hidden email]
https://lists.osgeo.org/mailman/listinfo/grass-user
Reply | Threaded
Open this post in threaded view
|

Re: question on i.nightlights.intercalibration

Gabriel Cotlier
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 <[hidden email]> wrote:
while for individual raster maps I get:

image.png

On Tue, Aug 27, 2019 at 9:25 AM Gabriel Cotlier <[hidden email]> wrote:
What should I do in this case?

On Tue, Aug 27, 2019 at 9:24 AM Gabriel Cotlier <[hidden email]> 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 <[hidden email]> wrote:


On Mon, Aug 26, 2019 at 1:14 AM Gabriel Cotlier <[hidden email]> 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
>
>                
>

_______________________________________________
grass-user mailing list
[hidden email]
https://lists.osgeo.org/mailman/listinfo/grass-user
Reply | Threaded
Open this post in threaded view
|

Re: question on i.nightlights.intercalibration

Markus Metz-3
In reply to this post by Gabriel Cotlier


On Tue, Aug 27, 2019 at 2:25 PM Gabriel Cotlier <[hidden email]> 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 <[hidden email]> wrote:
>>
>>
>>
>> On Mon, Aug 26, 2019 at 1:14 AM Gabriel Cotlier <[hidden email]> 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
>> >
>> >                
>> >
_______________________________________________
grass-user mailing list
[hidden email]
https://lists.osgeo.org/mailman/listinfo/grass-user
Reply | Threaded
Open this post in threaded view
|

Re: question on i.nightlights.intercalibration

Gabriel Cotlier
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?
Thanks a lot again for your help.
Regards,
Gabriel


image.png





On Tue, Aug 27, 2019 at 4:02 PM Markus Metz <[hidden email]> wrote:


On Tue, Aug 27, 2019 at 2:25 PM Gabriel Cotlier <[hidden email]> 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 <[hidden email]> wrote:
>>
>>
>>
>> On Mon, Aug 26, 2019 at 1:14 AM Gabriel Cotlier <[hidden email]> 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
>> >
>> >                
>> >
_______________________________________________
grass-user mailing list
[hidden email]
https://lists.osgeo.org/mailman/listinfo/grass-user
Reply | Threaded
Open this post in threaded view
|

Re: question on i.nightlights.intercalibration

Markus Metz-3


On Wed, Aug 28, 2019 at 12:54 PM Gabriel Cotlier <[hidden email]> 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 <[hidden email]> wrote:
>>
>>
>>
>> On Tue, Aug 27, 2019 at 2:25 PM Gabriel Cotlier <[hidden email]> 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 <[hidden email]> wrote:
>> >>
>> >>
>> >>
>> >> On Mon, Aug 26, 2019 at 1:14 AM Gabriel Cotlier <[hidden email]> 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
>> >> >
>> >> >                
>> >> >
_______________________________________________
grass-user mailing list
[hidden email]
https://lists.osgeo.org/mailman/listinfo/grass-user
Reply | Threaded
Open this post in threaded view
|

Re: 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.png

Thanks a lot!

 

On Wed, Aug 28, 2019 at 8:46 AM Markus Metz <[hidden email]> wrote:


On Wed, Aug 28, 2019 at 12:54 PM Gabriel Cotlier <[hidden email]> 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 <[hidden email]> wrote:
>>
>>
>>
>> On Tue, Aug 27, 2019 at 2:25 PM Gabriel Cotlier <[hidden email]> 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 <[hidden email]> wrote:
>> >>
>> >>
>> >>
>> >> On Mon, Aug 26, 2019 at 1:14 AM Gabriel Cotlier <[hidden email]> 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
>> >> >
>> >> >                
>> >> >
_______________________________________________
grass-user mailing list
[hidden email]
https://lists.osgeo.org/mailman/listinfo/grass-user