<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>