<div dir="ltr">Thanks a lot..!!!!!!!..Nikos<div>We will look into it</div><div>If we couldn't do it we will definitely ask your help</div><div><br></div><div>Thanks again for the help.</div><div><br></div><div>Best Regards</div><div>Kalindu Perera</div></div><br><div class="gmail_quote"><div dir="ltr">On Tue, 16 Oct 2018 at 09:39, Nikos Alexandris <<a href="mailto:nik@nikosalexandris.net">nik@nikosalexandris.net</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Kalindu Perera:<br>
<br>
>Commands are as follows<br>
><br>
>*r.ros,*<br>
>r.ros --overwrite model=fuel_model@demomapset<br>
>moisture_1h=1hour_moisture@demomapset moisture_live=live_moisture@demomapset<br>
>velocity=wind_speed@demomapset direction=wind_direction@demomapset<br>
>slope=slope@demomapset aspect=aspect@demomapset<br>
>elevation=elevation@demomapset base_ros=my_ros.base max_ros=my_ros.max<br>
>direction_ros=my_ros.maxdir spotting_distance=my_ros.spotdist<br>
><br>
>*r.spread,*<br>
>r.spread --overwrite base_ros=my_ros.base@demomapset<br>
>max_ros=my_ros.max@demomapset direction_ros=my_ros.maxdir@demomapset<br>
>start=fire_origin@demomapset wind_speed=wind_speed@demomapset<br>
>output=spread_time_observed<br>
><br>
>*After entering; r.category spread_time_observed,*<br>
>[image: 1.png]<br>
><br>
>*After entering; r.report spread_time_observed units=c,p,*<br>
>[image: Screenshot from 2018-10-15 21-40-46.png]<br>
<br>
Dear Kalindu,<br>
<br>
thank you for the above.<br>
<br>
If all is right, you might just need to properly prepare the input data<br>
and/or the output of `r.ros`, for example, to label raster categories<br>
appropriately.  Then, modules like `r.report` will show also the label<br>
for each raster category.<br>
<br>
When drawing the `spread_time_observed` map, can you identify visually the<br>
raster category of your interest?  If yes, pick up the the "Query<br>
raster/vector map(s)" tool (4th button in wxGUI's Map Display window),<br>
and select any pixel of the category of your interest.  This will tell<br>
the value (of the pixel you selected) and a label (if any is given).<br>
<br>
Knowing this value, you can extract all pixels of this value using<br>
`r.mapcalc`. For example, say the raster category (or value) of<br>
interest, is 3. Then, something like<br>
```<br>
r.mapcalc "category_3 = if(spread_time_observed == 3, 3, null())"<br>
```<br>
will create a map with all pixels of value 3.<br>
<br>
I guess using this "category_3" map as an input to `r.surf.area` would<br>
give the area estimation you require. Please, take care to read the<br>
NOTES in `r.surf`area`'s manual about "edge effects" though.<br>
<br>
<br>
Else, I think the easiest/fastest for me would to be to have access to sample<br>
data (even via a private message). If possible, cutting out (setting via<br>
`g.region`) a small fragment from your study area, and sharing would<br>
allow me to "resolve" this.<br>
<br>
Also, please don't post screenshots (if it is not about part of maps or other<br>
visuals). For textual output, please always post the output of your commands<br>
(as in copy and paste) or attach a text file.<br>
<br>
For example, I would copy-and-paste the complete output of (or put in a<br>
text file if a specific output of a command is really too long):<br>
<br>
# in Bash<br>
for MAP in fuel_model 1hour_moisture live_moisture wind_speed wind_direction slope aspect elevation my_ros.base my_ros.max my_ros.maxdir my_ros.spotdist fire_origin spread_time_observed ;do r.category $MAP ;done<br>
<br>
# in Windows<br>
see example in: <a href="https://lists.osgeo.org/pipermail/grass-user/2018-June/078526.html" rel="noreferrer" target="_blank">https://lists.osgeo.org/pipermail/grass-user/2018-June/078526.html</a><br>
<br>
Nikos<br>
</blockquote></div>