Hi, i forgot... below I have included the gdalinfo from both ouputs (ArcMap and GRASS). Any ideas as to
what else I should check would be appreciated. <br><br>#### ArcMap output gdalinfo ####<br>gdalinfo x:\test\fileASC_ARCMAP.asc <br>Driver: AAIGrid/Arc/Info ASCII Grid<br>Files: x:\test\fileASC_ARCMAP.asc <br> x:\test\fileASC_ARCMAP.asc.aux.xml<br>
Size is 5245, 4800<br>Coordinate System is `'<br>Origin = (-60.925000056119998,-0.000000095999999)<br>Pixel Size = (0.002083333330000,-0.002083333330000)<br>Metadata:<br> PyramidResamplingType=NEAREST<br>Corner Coordinates:<br>
Upper Left ( -60.9250001, -0.0000001) <br>Lower Left ( -60.9250001, -10.0000001) <br>Upper Right ( -49.9979167, -0.0000001) <br>Lower Right ( -49.9979167, -10.0000001) <br>Center ( -55.4614584, -5.0000001) <br>Band 1 Block=5245x1 Type=Int32, ColorInterp=Undefined<br>
NoData Value=-9999<br>(Tue May 08 11:21:09 2012) Command finished (10 sec) <br>(Tue May 08 11:21:42 2012) <br><br>#### GRASS output gdalinfo #### <br>
gdalinfo x:\test\fileASC_GRASS.asc <br>Driver: AAIGrid/Arc/Info ASCII Grid<br>Files: x:\test\fileASC_GRASS.asc <br> x:\test\fileASC_GRASS.asc.aux.xml<br>Size is 5246, 4801<br>Coordinate System is `'<br>
Origin = (-60.927083000000003,0.000483000000001)<br>Pixel Size = (0.002083000000000,-0.002083000000000)<br>Metadata:<br> PyramidResamplingType=NEAREST<br>Corner Coordinates:<br>Upper Left ( -60.9270830, 0.0004830) <br>
Lower Left ( -60.9270830, -10.0000000) <br>Upper Right ( -49.9996650, 0.0004830) <br>Lower Right ( -49.9996650, -10.0000000) <br>Center ( -55.4633740, -4.9997585) <br>Band 1 Block=5246x1 Type=Int32, ColorInterp=Undefined<br>
NoData Value=-9999<br><br><br><br><br><div class="gmail_quote">2012/5/8 Alejandro Coca Castro <span dir="ltr"><<a href="mailto:acocac@gmail.com" target="_blank">acocac@gmail.com</a>></span><br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
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>
</blockquote></div><br><br clear="all"><br>-- <br><b><i><span style="font-size:9.0pt;font-family:"Arial","sans-serif";color:#404040">Alejandro Coca Castro</span></i></b><b><i><span style="font-size:9.0pt;font-family:"Arial","sans-serif";color:black"></span></i></b><b><i><span style="font-size:9.0pt;font-family:"Arial","sans-serif";color:#00502e"></span></i></b><span style="font-size:9.0pt;font-family:"Arial","sans-serif";color:#404040"></span>
<b><i><span style="font-size:8.0pt;font-family:"Arial","sans-serif";color:black" lang="ES-CO"></span></i></b><i><span style="font-size:9.0pt;font-family:"Arial","sans-serif";color:#404040"></span></i><br>