[GRASS-user] Help - Output difference between ARCGIS enviroment settings and GRASS g.region function

Alejandro Coca Castro acocac at gmail.com
Tue May 8 12:09:21 EDT 2012


Hi, i tried to convert the following enviroments settings to GRASS g.region
and r.mask for the next ARCGIS output (ASCII file):

    fileTIF = r"x:\test\fileTIF.tif"
    MASK = r"x:\test\mask_file"
    env.extent = fileTIF
    env.snapRaster = MASK
    env.mask = MASK
    fileGRD = r"x:\test\fileGRD"
    outCon = Con(fileTIF, fileTIF, "", "VALUE >= 0")
    outCon.save(fileGRD)
    fileASC = r"x:\test\fileASC.asc"
    env.workspace = fileGRD
    RasterToASCII_conversion(fileGRD,fileASC)

So, i put on the GRASS command-line:

r.mask input=mask_file at default
g.region rast=fileTIF at default res=0.00208333333 align=MASK at default
r.mapcalc fileGRD = if(fileTIF at default  < 0 , null(), fileTIF at default)
r.out.arc input= fileGRD at default output=x:\test\fileASC_grass.asc

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)....

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)

If somebody knows how can i fix it, i'll be grateful,

Thanks in advance,

Alejandro
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://lists.osgeo.org/pipermail/grass-user/attachments/20120508/e514780a/attachment.html


More information about the grass-user mailing list