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

Markus Metz markus.metz.giswork at gmail.com
Wed Aug 28 04:46:31 PDT 2019


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/20190828/a46b3fe4/attachment.html>


More information about the grass-user mailing list