[R-sig-Geo] selecting a subset of a raster layer
Ahmadou Dicko
d|cko@@hm@dou @end|ng |rom gm@||@com
Tue Aug 2 17:46:56 CEST 2022
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]]
More information about the R-sig-Geo
mailing list