Locate the pixels which have the biggest standard deviation in a time-series layer stack raster using r [on hold]Calculating mean and standard deviation for raster values using whole pixels within a polygonR: raster stack (time series) loop pixelwise and write layer number of value occurence into a single output layerHow to store large time series of raster imagesPixel-wize trend analysis of an irregular time series stack of NDVI rasters in RAveraging MODIS 8-day LST for each pixel in a region to generate annual mean LSTUsing spei() function on time series from rasterstack in R?MODIS LST maximum value compositeFormatting h-axis of time series chart in Google Earth EngineDoes the MODIStsp package computes the QA's and eventually corrects the bad quality pixels?
If human space travel is limited by the G force vulnerability, is there a way to counter G forces?
iPad being using in wall mount battery swollen
What's the in-universe reasoning behind sorcerers needing material components?
In 'Revenger,' what does 'cove' come from?
ssTTsSTtRrriinInnnnNNNIiinngg
What is a romance in Latin?
How can saying a song's name be a copyright violation?
How would I stat a creature to be immune to everything but the Magic Missile spell? (just for fun)
Is it inappropriate for a student to attend their mentor's dissertation defense?
Why would the Red Woman birth a shadow if she worshipped the Lord of the Light?
Im going to France and my passport expires June 19th
Forgetting the musical notes while performing in concert
Avoiding direct proof while writing proof by induction
Do UK voters know if their MP will be the Speaker of the House?
What does “the session was packed” mean in this context?
What method can I use to design a dungeon difficult enough that the PCs can't make it through without killing them?
How does a predictive coding aid in lossless compression?
Why is it a bad idea to hire a hitman to eliminate most corrupt politicians?
Why are the 737's rear doors unusable in a water landing?
Alternative to sending password over mail?
How much of data wrangling is a data scientist's job?
How did the Super Star Destroyer Executor get destroyed exactly?
Plagiarism or not?
How to show a landlord what we have in savings?
Locate the pixels which have the biggest standard deviation in a time-series layer stack raster using r [on hold]
Calculating mean and standard deviation for raster values using whole pixels within a polygonR: raster stack (time series) loop pixelwise and write layer number of value occurence into a single output layerHow to store large time series of raster imagesPixel-wize trend analysis of an irregular time series stack of NDVI rasters in RAveraging MODIS 8-day LST for each pixel in a region to generate annual mean LSTUsing spei() function on time series from rasterstack in R?MODIS LST maximum value compositeFormatting h-axis of time series chart in Google Earth EngineDoes the MODIStsp package computes the QA's and eventually corrects the bad quality pixels?
I have 30 monthly tiff images of MOD11A2, and I created a time-series layer stuck using the code
mod = timeStackMODIS(list, pattern = glob2rx("\.tif$"))
using the bfastspatial package. I want to check which pixels have the biggest standard deviation, so I can focus my analysis in those pixels. I have made an attempt using this code
modsd = (list.files(path = "my_dir", full.names = FALSE, pattern = "\.tif$" ))
modsdstack <- stack(modsd)
testmodsd = calc(modsdstack, fun = sd)
and as a result the code returned a single raster, which I was unable to explain his raster values (because I had 30 images and the code returned one). So I think thats not the correct way to proceed. What I can not understand is if I have to create 30 new raster (each for every original image) which would have the sd as pixel values or I can, somehow, calculate the sd using the time-series layer stack which I have.
raster r spatial-statistics modis time-series
put on hold as unclear what you're asking by Spacedman, Erik, csk, Andre Silva, wetland yesterday
Please clarify your specific problem or add additional details to highlight exactly what you need. As it's currently written, it’s hard to tell exactly what you're asking. See the How to Ask page for help clarifying this question. If this question can be reworded to fit the rules in the help center, please edit the question.
|
show 5 more comments
I have 30 monthly tiff images of MOD11A2, and I created a time-series layer stuck using the code
mod = timeStackMODIS(list, pattern = glob2rx("\.tif$"))
using the bfastspatial package. I want to check which pixels have the biggest standard deviation, so I can focus my analysis in those pixels. I have made an attempt using this code
modsd = (list.files(path = "my_dir", full.names = FALSE, pattern = "\.tif$" ))
modsdstack <- stack(modsd)
testmodsd = calc(modsdstack, fun = sd)
and as a result the code returned a single raster, which I was unable to explain his raster values (because I had 30 images and the code returned one). So I think thats not the correct way to proceed. What I can not understand is if I have to create 30 new raster (each for every original image) which would have the sd as pixel values or I can, somehow, calculate the sd using the time-series layer stack which I have.
raster r spatial-statistics modis time-series
put on hold as unclear what you're asking by Spacedman, Erik, csk, Andre Silva, wetland yesterday
Please clarify your specific problem or add additional details to highlight exactly what you need. As it's currently written, it’s hard to tell exactly what you're asking. See the How to Ask page for help clarifying this question. If this question can be reworded to fit the rules in the help center, please edit the question.
Yes, I want to compare each pixel in a raster for a specific year against the same pixel of the other 29 overlapping years/rasters. It's not what I expected because in order to check which pixels have the biggest sd I need to know the mean (average) of the time-series layer stack, no?
– NikosGr
2 days ago
1
When you runsd
on a stack of 30 layers you get a single-band raster of thesd
of each set of 30 values at each pixel. You're not going to get 30 layers. Its like adding the pixel values, or computing their mean. You get one layer out for all the layers put in.
– Spacedman
2 days ago
Thank you for your fast reply. What I can not understand is how will I know which pixels have a big value of sd if I don't know what is the mean (average) of each pixel in the overlapping images? What I mean is in order to locate the pixels with high sd values, firstly, I need to know the mean (average) of each pixel for every year, am I right?
– NikosGr
2 days ago
What I thought is this: from the 8day products of MOD11A2 which I had, I should create monthly products using the sd and after to utilize thecalc
function to create a single image as I did before. Maybe in this way I will be able to understand the values.
– NikosGr
2 days ago
1
One should note that the sd is likely irrelevant in this case because this type of data is quite non-normal. If the skewness is high enough the sd could even take a value outside the observed distribution. Besides, you are not listing to the advice being presented and misunderstanding sd.
– Jeffrey Evans
2 days ago
|
show 5 more comments
I have 30 monthly tiff images of MOD11A2, and I created a time-series layer stuck using the code
mod = timeStackMODIS(list, pattern = glob2rx("\.tif$"))
using the bfastspatial package. I want to check which pixels have the biggest standard deviation, so I can focus my analysis in those pixels. I have made an attempt using this code
modsd = (list.files(path = "my_dir", full.names = FALSE, pattern = "\.tif$" ))
modsdstack <- stack(modsd)
testmodsd = calc(modsdstack, fun = sd)
and as a result the code returned a single raster, which I was unable to explain his raster values (because I had 30 images and the code returned one). So I think thats not the correct way to proceed. What I can not understand is if I have to create 30 new raster (each for every original image) which would have the sd as pixel values or I can, somehow, calculate the sd using the time-series layer stack which I have.
raster r spatial-statistics modis time-series
I have 30 monthly tiff images of MOD11A2, and I created a time-series layer stuck using the code
mod = timeStackMODIS(list, pattern = glob2rx("\.tif$"))
using the bfastspatial package. I want to check which pixels have the biggest standard deviation, so I can focus my analysis in those pixels. I have made an attempt using this code
modsd = (list.files(path = "my_dir", full.names = FALSE, pattern = "\.tif$" ))
modsdstack <- stack(modsd)
testmodsd = calc(modsdstack, fun = sd)
and as a result the code returned a single raster, which I was unable to explain his raster values (because I had 30 images and the code returned one). So I think thats not the correct way to proceed. What I can not understand is if I have to create 30 new raster (each for every original image) which would have the sd as pixel values or I can, somehow, calculate the sd using the time-series layer stack which I have.
raster r spatial-statistics modis time-series
raster r spatial-statistics modis time-series
edited 2 days ago
NikosGr
asked 2 days ago
NikosGrNikosGr
317
317
put on hold as unclear what you're asking by Spacedman, Erik, csk, Andre Silva, wetland yesterday
Please clarify your specific problem or add additional details to highlight exactly what you need. As it's currently written, it’s hard to tell exactly what you're asking. See the How to Ask page for help clarifying this question. If this question can be reworded to fit the rules in the help center, please edit the question.
put on hold as unclear what you're asking by Spacedman, Erik, csk, Andre Silva, wetland yesterday
Please clarify your specific problem or add additional details to highlight exactly what you need. As it's currently written, it’s hard to tell exactly what you're asking. See the How to Ask page for help clarifying this question. If this question can be reworded to fit the rules in the help center, please edit the question.
Yes, I want to compare each pixel in a raster for a specific year against the same pixel of the other 29 overlapping years/rasters. It's not what I expected because in order to check which pixels have the biggest sd I need to know the mean (average) of the time-series layer stack, no?
– NikosGr
2 days ago
1
When you runsd
on a stack of 30 layers you get a single-band raster of thesd
of each set of 30 values at each pixel. You're not going to get 30 layers. Its like adding the pixel values, or computing their mean. You get one layer out for all the layers put in.
– Spacedman
2 days ago
Thank you for your fast reply. What I can not understand is how will I know which pixels have a big value of sd if I don't know what is the mean (average) of each pixel in the overlapping images? What I mean is in order to locate the pixels with high sd values, firstly, I need to know the mean (average) of each pixel for every year, am I right?
– NikosGr
2 days ago
What I thought is this: from the 8day products of MOD11A2 which I had, I should create monthly products using the sd and after to utilize thecalc
function to create a single image as I did before. Maybe in this way I will be able to understand the values.
– NikosGr
2 days ago
1
One should note that the sd is likely irrelevant in this case because this type of data is quite non-normal. If the skewness is high enough the sd could even take a value outside the observed distribution. Besides, you are not listing to the advice being presented and misunderstanding sd.
– Jeffrey Evans
2 days ago
|
show 5 more comments
Yes, I want to compare each pixel in a raster for a specific year against the same pixel of the other 29 overlapping years/rasters. It's not what I expected because in order to check which pixels have the biggest sd I need to know the mean (average) of the time-series layer stack, no?
– NikosGr
2 days ago
1
When you runsd
on a stack of 30 layers you get a single-band raster of thesd
of each set of 30 values at each pixel. You're not going to get 30 layers. Its like adding the pixel values, or computing their mean. You get one layer out for all the layers put in.
– Spacedman
2 days ago
Thank you for your fast reply. What I can not understand is how will I know which pixels have a big value of sd if I don't know what is the mean (average) of each pixel in the overlapping images? What I mean is in order to locate the pixels with high sd values, firstly, I need to know the mean (average) of each pixel for every year, am I right?
– NikosGr
2 days ago
What I thought is this: from the 8day products of MOD11A2 which I had, I should create monthly products using the sd and after to utilize thecalc
function to create a single image as I did before. Maybe in this way I will be able to understand the values.
– NikosGr
2 days ago
1
One should note that the sd is likely irrelevant in this case because this type of data is quite non-normal. If the skewness is high enough the sd could even take a value outside the observed distribution. Besides, you are not listing to the advice being presented and misunderstanding sd.
– Jeffrey Evans
2 days ago
Yes, I want to compare each pixel in a raster for a specific year against the same pixel of the other 29 overlapping years/rasters. It's not what I expected because in order to check which pixels have the biggest sd I need to know the mean (average) of the time-series layer stack, no?
– NikosGr
2 days ago
Yes, I want to compare each pixel in a raster for a specific year against the same pixel of the other 29 overlapping years/rasters. It's not what I expected because in order to check which pixels have the biggest sd I need to know the mean (average) of the time-series layer stack, no?
– NikosGr
2 days ago
1
1
When you run
sd
on a stack of 30 layers you get a single-band raster of the sd
of each set of 30 values at each pixel. You're not going to get 30 layers. Its like adding the pixel values, or computing their mean. You get one layer out for all the layers put in.– Spacedman
2 days ago
When you run
sd
on a stack of 30 layers you get a single-band raster of the sd
of each set of 30 values at each pixel. You're not going to get 30 layers. Its like adding the pixel values, or computing their mean. You get one layer out for all the layers put in.– Spacedman
2 days ago
Thank you for your fast reply. What I can not understand is how will I know which pixels have a big value of sd if I don't know what is the mean (average) of each pixel in the overlapping images? What I mean is in order to locate the pixels with high sd values, firstly, I need to know the mean (average) of each pixel for every year, am I right?
– NikosGr
2 days ago
Thank you for your fast reply. What I can not understand is how will I know which pixels have a big value of sd if I don't know what is the mean (average) of each pixel in the overlapping images? What I mean is in order to locate the pixels with high sd values, firstly, I need to know the mean (average) of each pixel for every year, am I right?
– NikosGr
2 days ago
What I thought is this: from the 8day products of MOD11A2 which I had, I should create monthly products using the sd and after to utilize the
calc
function to create a single image as I did before. Maybe in this way I will be able to understand the values.– NikosGr
2 days ago
What I thought is this: from the 8day products of MOD11A2 which I had, I should create monthly products using the sd and after to utilize the
calc
function to create a single image as I did before. Maybe in this way I will be able to understand the values.– NikosGr
2 days ago
1
1
One should note that the sd is likely irrelevant in this case because this type of data is quite non-normal. If the skewness is high enough the sd could even take a value outside the observed distribution. Besides, you are not listing to the advice being presented and misunderstanding sd.
– Jeffrey Evans
2 days ago
One should note that the sd is likely irrelevant in this case because this type of data is quite non-normal. If the skewness is high enough the sd could even take a value outside the observed distribution. Besides, you are not listing to the advice being presented and misunderstanding sd.
– Jeffrey Evans
2 days ago
|
show 5 more comments
0
active
oldest
votes
0
active
oldest
votes
0
active
oldest
votes
active
oldest
votes
active
oldest
votes
Yes, I want to compare each pixel in a raster for a specific year against the same pixel of the other 29 overlapping years/rasters. It's not what I expected because in order to check which pixels have the biggest sd I need to know the mean (average) of the time-series layer stack, no?
– NikosGr
2 days ago
1
When you run
sd
on a stack of 30 layers you get a single-band raster of thesd
of each set of 30 values at each pixel. You're not going to get 30 layers. Its like adding the pixel values, or computing their mean. You get one layer out for all the layers put in.– Spacedman
2 days ago
Thank you for your fast reply. What I can not understand is how will I know which pixels have a big value of sd if I don't know what is the mean (average) of each pixel in the overlapping images? What I mean is in order to locate the pixels with high sd values, firstly, I need to know the mean (average) of each pixel for every year, am I right?
– NikosGr
2 days ago
What I thought is this: from the 8day products of MOD11A2 which I had, I should create monthly products using the sd and after to utilize the
calc
function to create a single image as I did before. Maybe in this way I will be able to understand the values.– NikosGr
2 days ago
1
One should note that the sd is likely irrelevant in this case because this type of data is quite non-normal. If the skewness is high enough the sd could even take a value outside the observed distribution. Besides, you are not listing to the advice being presented and misunderstanding sd.
– Jeffrey Evans
2 days ago