Error in geos_op2_geom("intersection", x, y, ...) : st_crs(x) == st_crs(y) is not TRUE

General
2

Without a reprex (see the FAQ), I haven't tested this, but it should work.

To conform the Coordinate Reference System (CRS) of a raster layer to another raster layer in R, you can use the projectRaster() function from the raster package[1][2]. Here's a step-by-step guide on how to do this:

  1. Load the raster package:
library(raster)

  1. Load or create the two raster objects, for example, raster1 and raster2. You want to conform the CRS of raster2 to match the CRS of raster1.

  2. Use the crs() function to obtain the CRS of raster1:

raster1_crs <- crs(raster1)
  1. Use the projectRaster() function to reproject raster2 to match the CRS of raster1:
raster2_conformed <- projectRaster(raster2, crs = raster1_crs)

Now, raster2_conformed has the same CRS as raster1, and you can perform further analysis or visualization with these raster layers[1][2].

Citations:
[1] Introduction to Geospatial Raster and Vector Data with R: Reproject Raster Data in R
[2] Raster 02: When Rasters Don't Line Up - Reproject Raster Data in R | NSF NEON | Open Data to Understand our Ecosystems

By Perplexity at For raster object layer in the R programming language, what is the function to set or change crs?

1 Like