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

[Ahmadou Dicko]

Hi Erin,

Not sure if understood well the problem, but it looks like a crop can do
the job here.
You'll need to build a rectangular polygon from the desired extent
(sf::st_bbox -> SpatialPolygon), and crop your raster with it.

Best,
Ahmadou





On Tue, Aug 2, 2022 at 3:37 PM Erin Hodgess <erinm.hodgess using gmail.com> wrote:

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


-- 
Ahmadou H. DICKO, PhD

Statistical consultant
Mobile: (+221) 77 123 81 69
Skype: dicko.ahmadou.h
Twitter : @dickoah
Gitlab: gitlab/dickoa
Github: github/dickoa

	[[alternative HTML version deleted]]