R help - Aug 2024 - allequal diff

Many thanks Ivan


Use is.na() on getValues() outputs, combine the two masks using the | operator
to get a mask of values that are missing in either raster, then negate the mask
to choose the non-missing values:

all.equal(getValues(r1)[!mask], getValues(r2)[!mask])

--> what do you mean by use is.na() in getValues(). So I need to call
getValues a second time? I suppose you mean to first prepare a mask using is.na
without getValues and then in the second step your code?

Kind regards
Sibylle

-----Original Message-----
From: Ivan Krylov <ikrylov at disroot.org> 
Sent: Friday, August 16, 2024 9:28 AM
To: SIBYLLE ST?CKLI via R-help <r-help at r-project.org>
Cc: sibylle.stoeckli at gmx.ch
Subject: Re: [R] allequal diff

? Fri, 16 Aug 2024 07:19:38 +0200
SIBYLLE ST?CKLI via R-help <r-help at r-project.org> ?????:
> Is it possible to consider na.rm=TRUE?
> > all.equal(getValues(r1), getValues(r2_resampled), tolerance = 0)
> 
> [1] "'is.NA' value mismatch: 9544032 in current 66532795 in
target"
Use is.na() on getValues() outputs, combine the two masks using the | operator
to get a mask of values that are missing in either raster, then negate the mask
to choose the non-missing values:

all.equal(getValues(r1)[!mask], getValues(r2)[!mask])

--
Best regards,
Ivan

? Fri, 16 Aug 2024 10:35:35 +0200
<sibylle.stoeckli at gmx.ch> ?????:
> what do you mean by use is.na() in getValues(). So I need to call
> getValues a second time?
Not necessarily, but it's one of the options. I was thinking along the
lines of:

values1 <- getValues(r1)
mask1 <- is.na(values1)
# Do the same for r2
# Combine the masks
all.equal(values1[!combined_mask], values2[!combined_mask])

Unlike compareRaster(), this assumes that the coordinate grid of r1 and
r2 is already the same and that only some of the values may differ.
> I suppose you mean to first prepare a mask using is.na without
> getValues and then in the second step your code?  
'raster' documentation says that is.na() works on raster objects, so it
should work.

Even if it didn't work, since you already access the underlying data
using getValues() and then compare the resulting vectors using
all.equal(), using is.na(getValues(...)) should definitely work.

-- 
Best regards,
Ivan

Here my idea including the error:

 
> m1=r1[r1="NA",]
> m2=r2_resampled[r2_resampled="NA",]
> 
> 
> all.equal(getValues(r1)[!m1], getValues(r2_resampled)[!m2])
[1] "Numeric: lengths (80706867, 65806339) differ"

 

-----Original Message-----
From: R-help <r-help-bounces at r-project.org> On Behalf Of SIBYLLE
ST?CKLI via R-help
Sent: Friday, August 16, 2024 10:36 AM
To: 'Ivan Krylov' <ikrylov at disroot.org>; 'SIBYLLE ST?CKLI
via R-help' <r-help at r-project.org>
Subject: Re: [R] allequal diff

 

Many thanks Ivan

 

 

Use is.na() on getValues() outputs, combine the two masks using the | operator
to get a mask of values that are missing in either raster, then negate the mask
to choose the non-missing values:

 

all.equal(getValues(r1)[!mask], getValues(r2)[!mask])

 

--> what do you mean by use is.na() in getValues(). So I need to call
getValues a second time? I suppose you mean to first prepare a mask using is.na
without getValues and then in the second step your code?

 

Kind regards

Sibylle

 

-----Original Message-----

From: Ivan Krylov < <mailto:ikrylov at disroot.org> ikrylov at
disroot.org>

Sent: Friday, August 16, 2024 9:28 AM

To: SIBYLLE ST?CKLI via R-help < <mailto:r-help at r-project.org>
r-help at r-project.org>

Cc:  <mailto:sibylle.stoeckli at gmx.ch> sibylle.stoeckli at gmx.ch

Subject: Re: [R] allequal diff

 

? Fri, 16 Aug 2024 07:19:38 +0200

SIBYLLE ST?CKLI via R-help < <mailto:r-help at r-project.org> r-help at
r-project.org> ?????:

 
> Is it possible to consider na.rm=TRUE?
 
> > all.equal(getValues(r1), getValues(r2_resampled), tolerance = 0)
> 
> [1] "'is.NA' value mismatch: 9544032 in current 66532795 in
target"
 

Use is.na() on getValues() outputs, combine the two masks using the | operator
to get a mask of values that are missing in either raster, then negate the mask
to choose the non-missing values:

 

all.equal(getValues(r1)[!mask], getValues(r2)[!mask])

 

--

Best regards,

Ivan

 

______________________________________________

 <mailto:R-help at r-project.org> R-help at r-project.org mailing list --
To UNSUBSCRIBE and more, see 
<https://stat.ethz.ch/mailman/listinfo/r-help>
https://stat.ethz.ch/mailman/listinfo/r-help

PLEASE do read the posting guide 
<http://www.R-project.org/posting-guide.html>
http://www.R-project.org/posting-guide.html

and provide commented, minimal, self-contained, reproducible code.


	[[alternative HTML version deleted]]