<div dir="ltr"><br><br>On Tue, Jun 26, 2018 at 9:53 AM, Nikos Alexandris <<a href="mailto:nik@nikosalexandris.net">nik@nikosalexandris.net</a>> wrote:<br>><br>> * Markus Metz <<a href="mailto:markus.metz.giswork@gmail.com">markus.metz.giswork@gmail.com</a>> [2018-06-25 08:29:45 +0200]:<br>><br>> [..]<br>><br>>> The resolution is a bit wrong, it is 0.008333333300000 but should be<br>>> 0.008333333333333, i.e. exactly 30 arc-seconds. This can be solved with the<br>>> -a flag of r.in.gdal, or after import with r.region -a.<br>>><br>>> The message<br>>><br>>>> 360 degree EW extent is exceeded by 0.999827 cells<br>>><br>>><br>>> will change to<br>>><br>>>> 360 degree EW extent is exceeded by 1 cells<br>>><br>>><br>>> but will not go away, because 360 degree EW extent is exceeded in the input<br>>> data, the first and last column cover the same geographical area. You can<br>>> change your current region to chop of e.g. the first column: set the region<br>>> to the raster, then modify the current region with g.region w=179:59:45W -p<br>>> and use this region for further processing.<br>><br>><br><div>> I guess this is worth being documented in the manual of the add-on.</div><div><br></div><div>This is a universal problem applying to various raster data in latlong. The first issue, 30 sec represented as 0.008333333300000 instead of 0.008333333333333 is solved by r.in.gdal -a. The second problem, this extra column responsible that 360 degree EW extent is exceeded by 1 cell can be solved by setting the current region accordingly. This is also a universal problem. Maybe the manual of r.in.gdal could include a hint about how to do this. Generally, users are encouraged to inspect the output of <a href="http://r.info">r.info</a> after importing raster data to check if everything is as expected.<br></div><div><br></div><div>> Would</div>> it also make sense to let the module attempt to perform this "correction"?<br><div><br></div><div>If you refer to i.nightlights.intercalibration, I would say no, because it is a more general issue not restricted to DMSP-OLS nightlight data. Even more general, the current region as set by the user is used for raster processing (with a very few exceptions).<br></div><div><br></div>Markus M<br></div>