[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