Hi, i tried to convert the following enviroments settings to GRASS g.region and r.mask for the next ARCGIS output (ASCII file):<br><br> fileTIF = r"x:\test\fileTIF.tif"<br> MASK = r"x:\test\mask_file" <br>
env.extent = fileTIF<br> env.snapRaster = MASK<br> env.mask = MASK<br> fileGRD = r"x:\test\fileGRD"<br> outCon = Con(fileTIF, fileTIF, "", "VALUE >= 0")<br> outCon.save(fileGRD)<br>
fileASC = r"x:\test\fileASC.asc"<br> env.workspace = fileGRD<br> RasterToASCII_conversion(fileGRD,fileASC)<br><br>So, i put on the GRASS command-line:<br><br>r.mask input=mask_file@default <br>g.region rast=fileTIF@default res=0.00208333333 align=MASK@default <br>
r.mapcalc fileGRD = if(fileTIF@default < 0 , null(), fileTIF@default) <br>r.out.arc input= fileGRD@default output=x:\test\fileASC_grass.asc<br> <br>However, the ArcGIS and GRASS outputs had a different extend and curiously the GRASS output didn't snap their cells with the MASK. Also, when i observe the GRASS output cellsize, it only has the first 6 decimals (f.e. 0.002083)....<br>
<br>Finally i tried to setup the parameters of g.region with the ArcGIS output but GRASS output definetly didnt have the same properties of ArcGIS output (cellsize, columns and rows, extend)<br><br>If somebody knows how can i fix it, i'll be grateful,<br>
<br>Thanks in advance,<br><br>Alejandro<br>