[R-sig-Geo] selecting a subset of a raster layer

[Erin Hodgess]

Thank you so much!

Sincerely,
Erin


On Tue, Aug 2, 2022 at 12:15 PM Bede-Fazekas Ákos <bfalevlist using gmail.com>
wrote:

> Dear Erin,
> I think that projectRaster(from = import_raster, to = r) or crop(x =
> import_raster, y = r) is what you need. crop() simply crop x with the
> extent of y, while projectRaster(from, to) create a raster from from
> that has exactly the same structure (resolution, origin, nrow, ncol,
> crs) as to.
> HTH,
> Ákos
> --------
> Ákos Bede-Fazekas
> Centre for Ecological Research, Hungary
>
> 2022.08.02. 17:37 keltezéssel, Erin Hodgess írta:
> > Hello!
> >
> > I have a raster with the following bbox:
> >
> >> bbox(import_raster)
> >            min       max
> > s1 -180.00417 180.00417
> > s2  -65.00417  75.00417
> >
> > I would like to create a subset of that raster layer with the following
> > dimensions:
> >> r
> > class      : RasterLayer
> > dimensions : 180, 360, 64800  (nrow, ncol, ncell)
> > resolution : 0.02222147, 0.003172611  (x, y)
> > extent     : -90.4184, -82.41867, 41.69613, 42.2672  (xmin, xmax, ymin,
> > ymax)
> > crs        : +proj=longlat +datum=WGS84 +no_defs
> >
> > Is there an easy way to do this, please?
> >
> > I did a "which" subset, but it produced an incredibly long vector.
> >
> > Any suggestions much appreciated.
> >
> > Sincerely,
> > Erin
> >
> >
> > Erin Hodgess, PhD
> > mailto: erinm.hodgess using gmail.com
> >
> >       [[alternative HTML version deleted]]
> >
> > _______________________________________________
> > R-sig-Geo mailing list
> > R-sig-Geo using r-project.org
> > https://stat.ethz.ch/mailman/listinfo/r-sig-geo
>
> _______________________________________________
> R-sig-Geo mailing list
> R-sig-Geo using r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-sig-geo
>
-- 
Erin Hodgess, PhD
mailto: erinm.hodgess using gmail.com

	[[alternative HTML version deleted]]