[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