<div dir="ltr">Hi Micha<div><br></div><div>  Â  Â Many thanks for your reply. </div><div>  Â  Â Here is the command I am using: </div><div><br></div><div>  Â  Â if(float(cat1_acc_riv) == float(cat1_minacc), str_r, null()) </div><div>  Â Â </div><div>  Â  Â  The str_r is a CELL raster. the result is different when I change it to: </div><div>  Â  Â 

 Â if(int(cat1_acc_riv) == int(cat1_minacc), str_r, null()) 

</div><div>  Â  Â  Â </div><div>  Â  Here is output of <a href="http://r.info">r.info</a> for first DCELL raster </div><div><br></div><div>+----------------------------------------------------------------------------+<br> | Map: Â  Â  Â cat1_acc_riv@PERMANENT Â  Â  Â  Â  Date: Sat Jan 23 22:58:42 2021 Â  Â |<br> | Mapset: Â  PERMANENT Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â Login of Creator: m43han Â  Â  Â  Â  Â |<br> | Location: main_working_location Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â |<br> | DataBase: C:\Users\m43han\Documents\Routing_Prod\Prod000001\grassdb Â  Â  Â  Â |<br> | Title: Â  Â cat1_acc_riv Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  |<br> | Timestamp: none Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â |<br> |----------------------------------------------------------------------------|<br> | Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â |<br> | Â  Type of Map: Â raster Â  Â  Â  Â  Â  Â  Â  Number of Categories: 19432056 Â  Â  Â  Â |<br> | Â  Data Type: Â  Â DCELL Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â |<br> | Â  Rows: Â  Â  Â  Â  4239 Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  |<br> | Â  Columns: Â  Â  Â 9254 Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  |<br> | Â  Total Cells: Â 39227706 Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  |<br> | Â  Â  Â  Â Projection: Latitude-Longitude Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â |<br> | Â  Â  Â  Â  Â  Â N: Â 50:52:39N Â  Â S: Â 40:16:48N Â  Res: 0:00:09 Â  Â  Â  Â  Â  Â  Â  Â  Â  |<br> | Â  Â  Â  Â  Â  Â E: Â 70:10:39W Â  Â W: Â 93:18:45W Â  Res: 0:00:09 Â  Â  Â  Â  Â  Â  Â  Â  Â  |<br> | Â  Range of data: Â  Â min = 250752 Â max = 19432056 Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  |<br> | Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â |<br> | Â  Data Description: Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â |<br> | Â  Â generated by r.mapcalc Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â |<br> | Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â |<br> | Â  Comments: Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â |<br> | Â  Â if(isnull(str_r), null(), acc) Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â |<br> | Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â |<br> +----------------------------------------------------------------------------+<br>(Sun Jan 24 04:45:38 2021) Command finished (0 sec)  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â <br></div><div><br></div><div><br></div><div>Here is <a href="http://r.info">r.info</a> output for second raster</div><div> +----------------------------------------------------------------------------+<br> | Map: Â  Â  Â cat1_minacc@PERMANENT Â  Â  Â  Â  Â Date: Sat Jan 23 22:58:48 2021 Â  Â |<br> | Mapset: Â  PERMANENT Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â Login of Creator: m43han Â  Â  Â  Â  Â |<br> | Location: main_working_location Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â |<br> | DataBase: C:\Users\m43han\Documents\Routing_Prod\Prod000001\grassdb Â  Â  Â  Â |<br> | Title: Â  Â cat1_minacc Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â |<br> | Timestamp: none Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â |<br> |----------------------------------------------------------------------------|<br> | Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â |<br> | Â  Type of Map: Â raster Â  Â  Â  Â  Â  Â  Â  Number of Categories: 0 Â  Â  Â  Â  Â  Â  Â  |<br> | Â  Data Type: Â  Â FCELL Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â |<br> | Â  Rows: Â  Â  Â  Â  4239 Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  |<br> | Â  Columns: Â  Â  Â 9254 Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  |<br> | Â  Total Cells: Â 39227706 Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  |<br> | Â  Â  Â  Â Projection: Latitude-Longitude Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â |<br> | Â  Â  Â  Â  Â  Â N: Â 50:52:39N Â  Â S: Â 40:16:48N Â  Res: 0:00:09 Â  Â  Â  Â  Â  Â  Â  Â  Â  |<br> | Â  Â  Â  Â  Â  Â E: Â 70:10:39W Â  Â W: Â 93:18:45W Â  Res: 0:00:09 Â  Â  Â  Â  Â  Â  Â  Â  Â  |<br> | Â  Range of data: Â  Â min = 250752 Â max = 1.817303e+007 Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â |<br> | Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â |<br> | Â  Data Description: Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â |<br> | Â  Â generated by r.stats.zonal Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â |<br> | Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â |<br> | Â  Comments: Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â |<br> | Â  Â r.stats.zonal --overwrite base="str_r" cover="cat1_acc_riv" method="\ Â  |<br> | Â  Â min" output="cat1_minacc" Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  |<br> | Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â  Â |<br> +----------------------------------------------------------------------------+<br>(Sun Jan 24 04:46:50 2021) Command finished (0 sec) <br></div><div><br></div><div>Thanks</div><div>Ming </div><div><br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">Micha Silver <<a href="mailto:tsvibar@gmail.com">tsvibar@gmail.com</a>> äºŽ2021å¹´1月24日周日 ä¸Šåˆ3:29写道:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">Is there some reason that you expect the rasters to be the same? Maybe<br>
begin by posting the outputs of `<a href="http://r.info" rel="noreferrer" target="_blank">r.info</a>` for both rasters, and what<br>
command you used to compare them?<br>
<br>
On Sun, Jan 24, 2021 at 6:13 AM ming han <<a href="mailto:dustming@gmail.com" target="_blank">dustming@gmail.com</a>> wrote:<br>
><br>
> Hi Everyone<br>
><br>
>  Â  Â I tried to compare if grids in two rasters are the same, one raster is FCELL and another raster is DCELL. I got different result  when I int both raster first before comparing them; and when I float() both raster first before I comparing them<br>
><br>
> Is there any reason for this?<br>
><br>
> Thanks<br>
> Ming<br>
> _______________________________________________<br>
> grass-user mailing list<br>
> <a href="mailto:grass-user@lists.osgeo.org" target="_blank">grass-user@lists.osgeo.org</a><br>
> <a href="https://lists.osgeo.org/mailman/listinfo/grass-user" rel="noreferrer" target="_blank">https://lists.osgeo.org/mailman/listinfo/grass-user</a><br>
<br>
<br>
<br>
-- <br>
Micha Silver<br>
Ben Gurion Univ<br>
Sde-Boker Remote Sensing Lab<br>
cell: +972 (52) 3665918<br>
</blockquote></div>