Message-ID: <CAASFQpShC9gNg2BZewB5ePpu3u7cUHH2eDjrq7fxPmx-xJGEOw@mail.gmail.com>
Date: 2016-02-19T11:09:05Z
From: Chris English
Subject: calculate "regional" slope
In-Reply-To: <CAF1jk_mQ0Jv1u+a7HnxwUDiBSo3RX1QHXgYaFjq+Cur3QX5_hg@mail.gmail.com>
Dominik,
r <- raster(nrows=22, ncols=20, xmn=-58, xmx=-48, ymn=-33, ymx=-22)
vals <- sample.int(1e3,440)
r[ ] <- vals
#raster::terrain
terr_r <- terrain(r, opt='slope', unit='degrees', neighbors=8)
Ah, but it appears you want up sampling to 1.5 degrees rather than 0.5 deg.
so maybe spatial.tools::projectRaster_rigorous then raster:terrain.
I'm inclined to end that last so maybe with a question mark. Sorry for an
essentially inconclusive response but I was happy to find terrain in any
case.
Chris
On Fri, Feb 19, 2016 at 2:59 AM, Dominik Schneider <
Dominik.Schneider at colorado.edu> wrote:
> I need to calculate slope at different scales. In the case below, r is a
> 0.5deg resolution raster and I want the slope for 1.5 deg centered on each
> of those 0.5 deg pixels. I'm trying to estimate which side of mountain
> range each pixel is on. So the resulting raster would have the same number
> of pixels as r. The edges can be NA.
> any suggestions would be appreciated. Thanks
>
>
> r <- raster(nrows=22, ncols=20, xmn=-58, xmx=-48, ymn=-33, ymx=-22)
> setValues(r,rnorm(440))
>
> [[alternative HTML version deleted]]
>
> _______________________________________________
> R-sig-Geo mailing list
> R-sig-Geo at r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-sig-geo
>
[[alternative HTML version deleted]]