Subsetting rasterstack for subtraction between elements using R?Correlation between rasters using RStep-by-step: How do I extract Raster values from Polygon overlay with Q-GIS or R?For loops + boxplot for rasterGrid vs. Polygon Overlay in R gives different results than in QGISChanging raster calculation values - shorten this script?Error while comparing two raster layer in RTransforming geostationary satellite image to lon/latMonitoring Progress of projectRaster in RSubset polygon data of SpatialPolygonsDataFrame by Distance to a Set of Coordinates in R?Using spei() function on time series from rasterstack in R?
Is it legal for company to use my work email to pretend I still work there?
Collect Fourier series terms
Approximately how much travel time was saved by the opening of the Suez Canal in 1869?
US citizen flying to France today and my passport expires in less than 2 months
What would happen to a modern skyscraper if it rains micro blackholes?
A newer friend of my brother's gave him a load of baseball cards that are supposedly extremely valuable. Is this a scam?
Mathematical cryptic clues
Pattern match does not work in bash script
Why don't electron-positron collisions release infinite energy?
What do you call a Matrix-like slowdown and camera movement effect?
LaTeX closing $ signs makes cursor jump
Arthur Somervell: 1000 Exercises - Meaning of this notation
The use of multiple foreign keys on same column in SQL Server
Is a tag line useful on a cover?
"You are your self first supporter", a more proper way to say it
can i play a electric guitar through a bass amp?
Have astronauts in space suits ever taken selfies? If so, how?
Dragon forelimb placement
Why doesn't Newton's third law mean a person bounces back to where they started when they hit the ground?
What's the point of deactivating Num Lock on login screens?
Today is the Center
How do I create uniquely male characters?
How to find program name(s) of an installed package?
Accidentally leaked the solution to an assignment, what to do now? (I'm the prof)
Subsetting rasterstack for subtraction between elements using R?
Correlation between rasters using RStep-by-step: How do I extract Raster values from Polygon overlay with Q-GIS or R?For loops + boxplot for rasterGrid vs. Polygon Overlay in R gives different results than in QGISChanging raster calculation values - shorten this script?Error while comparing two raster layer in RTransforming geostationary satellite image to lon/latMonitoring Progress of projectRaster in RSubset polygon data of SpatialPolygonsDataFrame by Distance to a Set of Coordinates in R?Using spei() function on time series from rasterstack in R?
.everyoneloves__top-leaderboard:empty,.everyoneloves__mid-leaderboard:empty,.everyoneloves__bot-mid-leaderboard:empty margin-bottom:0;
I want ro calculate the differences between rasters within a raster stack with a certain step, and I am unsure how to do so.
What I want is to subtract every element from its 5th next one, something like layer1-layer5, layer5-layer9, layer9 -layer13 and so on.
I tried to use the following to avoid loops but it does not work.
r = raster(matrix(1:30,5,6))
mstack<-stack(r,r,r,r,r,r,r,r,r,r,r,r,r,r,r,r,r,r,r,r,r,r,r,r)
msub<-subset(mstack, 1:nlayers(mstack)) - subset(mstack, 1:(nlayers(mstack)+5))
r element rasterstack
add a comment |
I want ro calculate the differences between rasters within a raster stack with a certain step, and I am unsure how to do so.
What I want is to subtract every element from its 5th next one, something like layer1-layer5, layer5-layer9, layer9 -layer13 and so on.
I tried to use the following to avoid loops but it does not work.
r = raster(matrix(1:30,5,6))
mstack<-stack(r,r,r,r,r,r,r,r,r,r,r,r,r,r,r,r,r,r,r,r,r,r,r,r)
msub<-subset(mstack, 1:nlayers(mstack)) - subset(mstack, 1:(nlayers(mstack)+5))
r element rasterstack
1
I think you want this:msub < -subset(mstack, 1:(nlayers(mstack)-4)) - subset(mstack, 5:nlayers(mstack))
In your second subset you are adding five to the total number of layers which is what is causing the invalid subset error.
– Jacob F
Apr 2 at 13:42
1
Do you also want layer2-layer6? All layer[n]-layer[n+4] for all n such that [n+4]<= total number of layers
– Spacedman
Apr 2 at 14:32
@Spacedman yes, apologies for my english, this is what I try to do but I messed up the subsets.
– foo
Apr 2 at 14:34
1
Do you need to usesubset
? Can you directly access layers thus:nl=nlayers(mstack); mstack[[1:(nl-4)]] - mstack[[5:nl]]
? You can see the pairs of layers being subtracted withcbind(1:(nl-4), 5:nl)
– Spacedman
Apr 2 at 14:39
Thank you @Spacedman, it took me some hours but it worked! Can you please post this as answer to accept it?
– foo
Apr 3 at 10:20
add a comment |
I want ro calculate the differences between rasters within a raster stack with a certain step, and I am unsure how to do so.
What I want is to subtract every element from its 5th next one, something like layer1-layer5, layer5-layer9, layer9 -layer13 and so on.
I tried to use the following to avoid loops but it does not work.
r = raster(matrix(1:30,5,6))
mstack<-stack(r,r,r,r,r,r,r,r,r,r,r,r,r,r,r,r,r,r,r,r,r,r,r,r)
msub<-subset(mstack, 1:nlayers(mstack)) - subset(mstack, 1:(nlayers(mstack)+5))
r element rasterstack
I want ro calculate the differences between rasters within a raster stack with a certain step, and I am unsure how to do so.
What I want is to subtract every element from its 5th next one, something like layer1-layer5, layer5-layer9, layer9 -layer13 and so on.
I tried to use the following to avoid loops but it does not work.
r = raster(matrix(1:30,5,6))
mstack<-stack(r,r,r,r,r,r,r,r,r,r,r,r,r,r,r,r,r,r,r,r,r,r,r,r)
msub<-subset(mstack, 1:nlayers(mstack)) - subset(mstack, 1:(nlayers(mstack)+5))
r element rasterstack
r element rasterstack
edited Apr 2 at 21:49
PolyGeo♦
53.9k1781245
53.9k1781245
asked Apr 2 at 13:31
foofoo
555
555
1
I think you want this:msub < -subset(mstack, 1:(nlayers(mstack)-4)) - subset(mstack, 5:nlayers(mstack))
In your second subset you are adding five to the total number of layers which is what is causing the invalid subset error.
– Jacob F
Apr 2 at 13:42
1
Do you also want layer2-layer6? All layer[n]-layer[n+4] for all n such that [n+4]<= total number of layers
– Spacedman
Apr 2 at 14:32
@Spacedman yes, apologies for my english, this is what I try to do but I messed up the subsets.
– foo
Apr 2 at 14:34
1
Do you need to usesubset
? Can you directly access layers thus:nl=nlayers(mstack); mstack[[1:(nl-4)]] - mstack[[5:nl]]
? You can see the pairs of layers being subtracted withcbind(1:(nl-4), 5:nl)
– Spacedman
Apr 2 at 14:39
Thank you @Spacedman, it took me some hours but it worked! Can you please post this as answer to accept it?
– foo
Apr 3 at 10:20
add a comment |
1
I think you want this:msub < -subset(mstack, 1:(nlayers(mstack)-4)) - subset(mstack, 5:nlayers(mstack))
In your second subset you are adding five to the total number of layers which is what is causing the invalid subset error.
– Jacob F
Apr 2 at 13:42
1
Do you also want layer2-layer6? All layer[n]-layer[n+4] for all n such that [n+4]<= total number of layers
– Spacedman
Apr 2 at 14:32
@Spacedman yes, apologies for my english, this is what I try to do but I messed up the subsets.
– foo
Apr 2 at 14:34
1
Do you need to usesubset
? Can you directly access layers thus:nl=nlayers(mstack); mstack[[1:(nl-4)]] - mstack[[5:nl]]
? You can see the pairs of layers being subtracted withcbind(1:(nl-4), 5:nl)
– Spacedman
Apr 2 at 14:39
Thank you @Spacedman, it took me some hours but it worked! Can you please post this as answer to accept it?
– foo
Apr 3 at 10:20
1
1
I think you want this:
msub < -subset(mstack, 1:(nlayers(mstack)-4)) - subset(mstack, 5:nlayers(mstack))
In your second subset you are adding five to the total number of layers which is what is causing the invalid subset error.– Jacob F
Apr 2 at 13:42
I think you want this:
msub < -subset(mstack, 1:(nlayers(mstack)-4)) - subset(mstack, 5:nlayers(mstack))
In your second subset you are adding five to the total number of layers which is what is causing the invalid subset error.– Jacob F
Apr 2 at 13:42
1
1
Do you also want layer2-layer6? All layer[n]-layer[n+4] for all n such that [n+4]<= total number of layers
– Spacedman
Apr 2 at 14:32
Do you also want layer2-layer6? All layer[n]-layer[n+4] for all n such that [n+4]<= total number of layers
– Spacedman
Apr 2 at 14:32
@Spacedman yes, apologies for my english, this is what I try to do but I messed up the subsets.
– foo
Apr 2 at 14:34
@Spacedman yes, apologies for my english, this is what I try to do but I messed up the subsets.
– foo
Apr 2 at 14:34
1
1
Do you need to use
subset
? Can you directly access layers thus: nl=nlayers(mstack); mstack[[1:(nl-4)]] - mstack[[5:nl]]
? You can see the pairs of layers being subtracted with cbind(1:(nl-4), 5:nl)
– Spacedman
Apr 2 at 14:39
Do you need to use
subset
? Can you directly access layers thus: nl=nlayers(mstack); mstack[[1:(nl-4)]] - mstack[[5:nl]]
? You can see the pairs of layers being subtracted with cbind(1:(nl-4), 5:nl)
– Spacedman
Apr 2 at 14:39
Thank you @Spacedman, it took me some hours but it worked! Can you please post this as answer to accept it?
– foo
Apr 3 at 10:20
Thank you @Spacedman, it took me some hours but it worked! Can you please post this as answer to accept it?
– foo
Apr 3 at 10:20
add a comment |
0
active
oldest
votes
Your Answer
StackExchange.ready(function()
var channelOptions =
tags: "".split(" "),
id: "79"
;
initTagRenderer("".split(" "), "".split(" "), channelOptions);
StackExchange.using("externalEditor", function()
// Have to fire editor after snippets, if snippets enabled
if (StackExchange.settings.snippets.snippetsEnabled)
StackExchange.using("snippets", function()
createEditor();
);
else
createEditor();
);
function createEditor()
StackExchange.prepareEditor(
heartbeatType: 'answer',
autoActivateHeartbeat: false,
convertImagesToLinks: false,
noModals: true,
showLowRepImageUploadWarning: true,
reputationToPostImages: null,
bindNavPrevention: true,
postfix: "",
imageUploader:
brandingHtml: "Powered by u003ca class="icon-imgur-white" href="https://imgur.com/"u003eu003c/au003e",
contentPolicyHtml: "User contributions licensed under u003ca href="https://creativecommons.org/licenses/by-sa/3.0/"u003ecc by-sa 3.0 with attribution requiredu003c/au003e u003ca href="https://stackoverflow.com/legal/content-policy"u003e(content policy)u003c/au003e",
allowUrls: true
,
onDemand: true,
discardSelector: ".discard-answer"
,immediatelyShowMarkdownHelp:true
);
);
Sign up or log in
StackExchange.ready(function ()
StackExchange.helpers.onClickDraftSave('#login-link');
);
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
Post as a guest
Required, but never shown
StackExchange.ready(
function ()
StackExchange.openid.initPostLogin('.new-post-login', 'https%3a%2f%2fgis.stackexchange.com%2fquestions%2f317500%2fsubsetting-rasterstack-for-subtraction-between-elements-using-r%23new-answer', 'question_page');
);
Post as a guest
Required, but never shown
0
active
oldest
votes
0
active
oldest
votes
active
oldest
votes
active
oldest
votes
Thanks for contributing an answer to Geographic Information Systems Stack Exchange!
- Please be sure to answer the question. Provide details and share your research!
But avoid …
- Asking for help, clarification, or responding to other answers.
- Making statements based on opinion; back them up with references or personal experience.
To learn more, see our tips on writing great answers.
Sign up or log in
StackExchange.ready(function ()
StackExchange.helpers.onClickDraftSave('#login-link');
);
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
Post as a guest
Required, but never shown
StackExchange.ready(
function ()
StackExchange.openid.initPostLogin('.new-post-login', 'https%3a%2f%2fgis.stackexchange.com%2fquestions%2f317500%2fsubsetting-rasterstack-for-subtraction-between-elements-using-r%23new-answer', 'question_page');
);
Post as a guest
Required, but never shown
Sign up or log in
StackExchange.ready(function ()
StackExchange.helpers.onClickDraftSave('#login-link');
);
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
Post as a guest
Required, but never shown
Sign up or log in
StackExchange.ready(function ()
StackExchange.helpers.onClickDraftSave('#login-link');
);
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
Post as a guest
Required, but never shown
Sign up or log in
StackExchange.ready(function ()
StackExchange.helpers.onClickDraftSave('#login-link');
);
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
Post as a guest
Required, but never shown
Required, but never shown
Required, but never shown
Required, but never shown
Required, but never shown
Required, but never shown
Required, but never shown
Required, but never shown
Required, but never shown
1
I think you want this:
msub < -subset(mstack, 1:(nlayers(mstack)-4)) - subset(mstack, 5:nlayers(mstack))
In your second subset you are adding five to the total number of layers which is what is causing the invalid subset error.– Jacob F
Apr 2 at 13:42
1
Do you also want layer2-layer6? All layer[n]-layer[n+4] for all n such that [n+4]<= total number of layers
– Spacedman
Apr 2 at 14:32
@Spacedman yes, apologies for my english, this is what I try to do but I messed up the subsets.
– foo
Apr 2 at 14:34
1
Do you need to use
subset
? Can you directly access layers thus:nl=nlayers(mstack); mstack[[1:(nl-4)]] - mstack[[5:nl]]
? You can see the pairs of layers being subtracted withcbind(1:(nl-4), 5:nl)
– Spacedman
Apr 2 at 14:39
Thank you @Spacedman, it took me some hours but it worked! Can you please post this as answer to accept it?
– foo
Apr 3 at 10:20