[GRASS-user] R: Re: r.hazard.flood white map flood output

Sure,

g.region -p

projection: 1 (UTM)

zone: 32

datum: wgs84

ellipsoid: wgs84

north: 5110010

south: 5103990

west: 621990

east: 634010

nsres: 1.204

ewres: 1.202

rows: 5000

cols: 10000

cells: 50000000

r.info bacino_dtm

±---------------------------------------------------------------------------+

|

|----------------------------------------------------------------------------|

| |

| Type of Map: raster Number of Categories: 255 |

| Data Type: FCELL |

| Rows: 5000 |

| Columns: 10000 |

| Total Cells: 50000000 |

| Projection: UTM (zone 32) |

| N: 5110010 S: 5103990 Res: 1.204 |

| E: 634010 W: 621990 Res: 1.202 |

| Range of data: min = 1026.483 max = 3465.358 |

| |

| Data Description: |

| generated by r.mapcalc |

| |

| Comments: |

| |

| |

±---------------------------------------------------------------------------+

r.hazard.flood.PY --overwrite map=zzzzzzz flood=zzzzz_flood mti=zzzzzzzzzz_MTI
Cellsize : 1.203
SEZIONE 1a (di 4): Inizializando Memoria.
SECTION 1b (of 4): Determining Offmap Flow.
SEZIONE 2: A * Search.
SECTION 3: Accumulating Surface Flow with MFD.
SEZIONE 4: Chiudendo Mappe.
WARNING: Writing out only positive flow accumulation values.
WARNING: Cells with a likely underestimate for flow accumulation can no longer be identified.
Flow accumulation done.
Mappa raster delle pendenze <r_slope> completata
Slope map done.
Exponent : 0.0174197601602
MTI threshold : 2.47170118814
Calculating mti raster map…
Cleaning up…
Removing raster <r_accumulation>
Removing raster <r_slope>
Calculating flood raster map…
Running r.clump…
Pass 1…
Pass 2…
r.clump completo. 38808 clumps.
Deleting areas of less than 690985 cells…
Reading…
Writing…
Done!
Cleaning up…
Removing raster <r_clump>
Removing raster <r_flood_th>
Removing raster <r_flood>
Done!

thanks…

m

Thank you Mattia, and what does

r.info zzzzzzzzzz_MTI

say?

···

Best regards,

Margherita DI LEO
Postdoctoral Researcher

European Commission - DG JRC
Institute for Environment and Sustainability (IES)
Via Fermi, 2749
I-21027 Ispra (VA) - Italy - TP 261

Tel. +39 0332 78 3600
margherita.di-leo@jrc.ec.europa.eu

Disclaimer: The views expressed are purely those of the writer and may not in any circumstance be regarded as stating an official position of the European Commission.

On Fri, Jun 14, 2013 at 10:57 AM, M. Riccadonna <Mattia_r@libero.it> wrote:

Sure,

g.region -p

projection: 1 (UTM)

...

nsres: 1.204
ewres: 1.202
rows: 5000
cols: 10000
cells: 50000000

Are you sure that you want
- non-square cells
- 1.2 meter resolution
?

Please see
http://grasswiki.osgeo.org/wiki/Computational_region

r.info bacino_dtm

...

| Total Cells: 50000000

...

| N: 5110010 S: 5103990 Res: 1.204

...

| E: 634010 W: 621990 Res: 1.202

... same here.
It could well be, just these 1m DEMs are stil rare :slight_smile:

| Range of data: min = 1026.483 max = 3465.358

... looks like an Alpine region... does flooding risk occur there?

ciao
Markus