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-boun...@r-project.org> On Behalf Of SIBYLLE STÖCKLI via R-help Sent: Friday, August 16, 2024 10:36 AM To: 'Ivan Krylov' <ikry...@disroot.org>; 'SIBYLLE STÖCKLI via R-help' <r-help@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:ikry...@disroot.org> ikry...@disroot.org> Sent: Friday, August 16, 2024 9:28 AM To: SIBYLLE STÖCKLI via R-help < <mailto:r-help@r-project.org> r-help@r-project.org> Cc: <mailto:sibylle.stoec...@gmx.ch> sibylle.stoec...@gmx.ch Subject: Re: [R] allequal diff В Fri, 16 Aug 2024 07:19:38 +0200 SIBYLLE STÖCKLI via R-help < <mailto:r-help@r-project.org> r-help@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@r-project.org> R-help@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]] ______________________________________________ R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.
