<div dir="ltr">Dear Vero and Markus,<div><br><div>Thanks a lot for the guidance and help.</div><div><br></div><div>If I understood correctly:</div><div><br></div><div>I first have to import one layer when I select a location setting the coordinate system of that layer to the mapset of GRASS (at the start).</div><div><div><br></div></div><div>Then I have to import all other layers to mapset using command "r.in.gdal -a" </div><div><br></div><div>r.in.gdal -a input=C:\Users\Gabriel\Documents\layers\F121995.tif output=F121995<br></div><div><br></div><div><span style="font-size:small;background-color:rgb(255,255,255);text-decoration-style:initial;text-decoration-color:initial;float:none;display:inline">Then I have to running "g.region" to set up the same extent to all layers as follows:</span>
<br></div><div><br></div><div>
<span style="font-size:12.8px;background-color:rgb(255,255,255);text-decoration-style:initial;text-decoration-color:initial;float:none;display:inline">g.region raster=F121995</span><br></div><div><br></div><div>and I got the layer processed and finished but the error did not went way....</div><div><br></div><div><div>|i Inter-satellite calibration of DMSP-OLS Nighttime Stable Lights</div><div>WARNING: Operating on current region</div><div><br></div><div>|> Calibrating average visible Digital Number values</div><div>Regression coefficients: (-0.0515, 1.2293, -0.0038) | Associated R^2: 0.9178</div><div><br></div><div>360 degree EW extent is exceeded by 0.999827 cells</div><div>360 degree EW extent is exceeded by 1 cells</div><div><br></div><div>WARNING: As requested, timestamp transferring not attempted.</div><div>WARNING: No data base element files found</div><div><br></div><div>(Sun Jun 24 17:05:31 2018) Command finished (1 min 2 sec) </div></div><div><br></div><div>What could be doing wrong?</div><div>In addition I'm doing it for one layer here, is there a way to do it for many layers at onces with command line or GUI dialog box?</div><div><br></div><div>Thanks a lot again.</div><div><br></div><div>Gabriel </div></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Sat, Jun 23, 2018 at 6:31 AM, Markus Metz <span dir="ltr"><<a href="mailto:markus.metz.giswork@gmail.com" target="_blank">markus.metz.giswork@gmail.com</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div dir="ltr"><span class=""><br><br>On Sat, Jun 23, 2018 at 11:10 AM, Veronica Andreo <<a href="mailto:veroandreo@gmail.com" target="_blank">veroandreo@gmail.com</a>> wrote:<br>><br>> Hi,<br>><br>> El sáb., 23 jun. 2018 4:35, Markus Metz <<a href="mailto:markus.metz.giswork@gmail.com" target="_blank">markus.metz.giswork@gmail.com</a><wbr>> escribió:<br>>><br>>><br>>><br>>> On Fri, Jun 22, 2018 at 10:34 PM, Veronica Andreo <<a href="mailto:veroandreo@gmail.com" target="_blank">veroandreo@gmail.com</a>> wrote:<br>>> ><br>>> > Hi Gabriel,<br>>> ><br>>> > What you could do is import with r.in.gdal -a that adjusts resolution for lat long maps [0]<br>>><br>>> that will help to fix the resolution from 0.008333333300000 to 0.008333333333333, i.e. exactly 30 arc-seconds. The software used to create the raster data has stored the resolution with limited precision.<br>><br>><br>> Right, I overlooked this <br>><br>>> > and then (before the intercalibration step), set the region to one of the imported maps with g.region raster=yourmap<br>>><br>>> you will then get a message like<br>>> 360 degree EW extent is exceeded by 1 cells<br>>><br>>> which is correct because the first and last column are duplicates. The cell centers cover -180, 180, and the EW extents regarding cell borders are E: 180:00:15E, W: 180:00:15W, grown by half a cell, i.e. 15 arc-seconds.<br>><br>><br>> So, solution is to just use the data with the extra 15 arc-seconds in each side? <br><div><br></div></span><div>yes, or chop off the first or last column: set the region to the raster, then modify the current region with g.region w=179:59:45W -p<br></div><span class=""><div>></div><div>> Of I want data to fit 180/-180, would r.region mymap e=180 w=180 help or will it change the data?</div><div><br></div></span><div>this would modify the data because</div><div>1) the raster will be shifted by half a cell to the east</div><div>2) the cell size (ew resolution) will be changed<br></div><div><br></div><div>Markus M</div><div><br></div>><br>> Vero <br><br></div>
</blockquote></div><br></div>