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?
What is the word for reserving something for yourself before others do?
How much RAM could one put in a typical 80386 setup?
How does one intimidate enemies without having the capacity for violence?
I'm flying to France today and my passport expires in less than 2 months
How do I deal with an unproductive colleague in a small company?
Does an object always see its latest internal state irrespective of thread?
Does detail obscure or enhance action?
High voltage LED indicator 40-1000 VDC without additional power supply
What does it mean to describe someone as a butt steak?
Could an aircraft fly or hover using only jets of compressed air?
Why "Having chlorophyll without photosynthesis is actually very dangerous" and "like living with a bomb"?
How is the claim "I am in New York only if I am in America" the same as "If I am in New York, then I am in America?
What typically incentivizes a professor to change jobs to a lower ranking university?
Theorems that impeded progress
Cross compiling for RPi - error while loading shared libraries
Replacing matching entries in one column of a file by another column from a different file
Is it tax fraud for an individual to declare non-taxable revenue as taxable income? (US tax laws)
Has there ever been an airliner design involving reducing generator load by installing solar panels?
Perform and show arithmetic with LuaLaTeX
If human space travel is limited by the G force vulnerability, is there a way to counter G forces?
How to determine what difficulty is right for the game?
Why can't I see bouncing of a switch on an oscilloscope?
What are these boxed doors outside store fronts in New York?
Modeling an IP Address
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
edited Apr 2 at 21:49
PolyGeo♦
53.9k1781245
asked Apr 2 at 13:31
foofoo
555
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
edited Apr 2 at 21:49
PolyGeo♦
53.9k1781245
asked Apr 2 at 13:31
foofoo
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
2 days ago
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
edited Apr 2 at 21:49
PolyGeo♦
53.9k1781245
asked Apr 2 at 13:31
foofoo
555
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
asked Apr 2 at 13:31
foofoo
555
edited Apr 2 at 21:49
PolyGeo♦
53.9k1781245
asked Apr 2 at 13:31
foofoo
555
edited Apr 2 at 21:49
PolyGeo♦
53.9k1781245
edited Apr 2 at 21:49
PolyGeo♦
53.9k1781245
edited Apr 2 at 21:49
PolyGeo♦
53.9k1781245
53.9k1781245
asked Apr 2 at 13:31
foofoo
555
asked Apr 2 at 13:31
foofoo
555
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
2 days ago
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
2 days ago
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
2 days ago
Thank you @Spacedman, it took me some hours but it worked! Can you please post this as answer to accept it?
– foo
2 days ago
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
2 days ago