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@default
g.region rast=fileTIF@default res=0.00208333333 align=MASK@default
r.mapcalc fileGRD = if(fileTIF@default < 0 , null(), fileTIF@default)
r.out.arc input= fileGRD@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
_______________________________________________
grass-user mailing list
[email protected]
http://lists.osgeo.org/mailman/listinfo/grass-user
