making sense of r.reclass.area

Dear GRASS list,

I'm having trouble understanding how to properly use ~r.reclass.area~. It seems to give inconsistent results. I'm using the following MWE to help me figure this out.

Several questions... such as:

Why are the units in hectare instead of m^2? 

Why is -c and -d mutually exclusive for "mode=reclass" but not "mode=rmarea"?

Why does mode=rmarea appear to ignore the -d flag

Why does mode=reclass say it cannot find areas <=1 ha, but only if -c is used, and the module never seems to have trouble finding areas of 1 ha otherwise?

Thanks,

  -k.
  

grass72 -c EPSG:3413 ./Gclump.debug
# frink "1 hectare -> m^2" 10,000 m^2
g.region w=0 e=500 s=0 n=500 res=100 -p

debug() { r.out.xyz -i input=$@  | cut -d"|" -f3 | xargs -n5; }

r.mapcalc "foo = if(row() == 2 && ((col() > 1) && (col() < 4)) , 1, null())" --o
r.mapcalc "foo = if(row() == 3 && col() == 4, 1, foo)" --o
r.mapcalc "foo = if(row() == 4 && col() == 2, 1, foo)" --o

debug foo
# * * * * *
# * 1 1 * *
# * * * 1 *
# * 1 * * *
# * * * * *

# RMAREA

# simple run with no args
r.reclass.area input=foo output=o value=1 mode=lesser method=rmarea --o --q && debug o
r.reclass.area -c input=foo output=o value=1 mode=lesser method=rmarea --o --q && debug o
# WANT: include diagonal
r.reclass.area -d input=foo output=o value=1 mode=lesser method=rmarea --o --q && debug o
r.reclass.area -d -c input=foo output=o value=1 mode=lesser method=rmarea --o --q && debug o

# All produce the same output:
# * * * * *
# * 1 1 * *
# * * * * *
# * * * * *
# * * * * *

# RECLASS

r.reclass.area input=foo output=o value=1 mode=lesser method=reclass --o --q && debug o
# GOOD: flags areas <=1. Does not consider diagonal:
# * * * * *
# * * * * *
# * * * 1 *
# * 1 * * *
# * * * * *

r.reclass.area -c input=foo output=o value=1 mode=lesser method=reclass --o --q && debug o
# ERROR: "No areas <= 1.000 hectar

r.reclass.area -d input=foo output=o value=1 mode=lesser method=reclass --o --q && debug o
# Appears to work, only areas <=1 flagged, including diagonal
# * * * * *
# * * * * *
# * * * * *
# * 1 * * *
# * * * * *

r.reclass.area -c -d input=foo output=o value=1 mode=lesser method=reclass --o --q && debug o
# ERROR: flags c and d are mutually exclusive

exit
rm -fR  Gclump.debug
_______________________________________________
grass-user mailing list
[email protected]
https://lists.osgeo.org/mailman/listinfo/grass-user