Questions tagged [stack]

stack in this context is intended for questions relating to the stacking of rasters, it is not intended for questions relating to a GIS software stack

Filter by
Sorted by
Tagged with
20
votes
4answers
16k views

Layer stacking rasters in QGIS?

I have some multi-spectral data, where the individual bands are separate Raster files. How do I layer stack these rasters so that I get one raster with 3 bands instead of 3 different rasters? This ...
9
votes
1answer
4k views

How to stack bands in Google Earth Engine?

I have created an Image collection in GEE and with the help of a function I have computed NDVI index and mapped it to create another collection with NDVI as a band. Now I want the create a stacked ...
7
votes
1answer
343 views

Stacking Vector Layers in TileStache

I am serving several vector tile layers of different geometries (point, line & polygon) in TileStache using VecTiles:Provider, with this cfgfile: { "cache": { "name": "DISK", "path":"...
6
votes
2answers
3k views

Mean by month on R stacked raster

I work with MODIS NDVI rasters in 2016. I have 23 rasters stacked in one object. I have 2 raster by month. I would like the average by months and conserve a raster for each month. ndvi.stack <- ...
5
votes
1answer
2k views

how to rename the band-names of a layer stack?

I created a layer stack of 20 images (all GTiff-files) using the gdal_merge tool in QGIS, available at: raster > miscellaneous > merge The result is a GTiff layer stack of the 20 images, but ...
5
votes
0answers
198 views

Stacking and masking Landsat at the same time with rasterio

Based on answers from previous questions 1,2 and the Masking a raster using a shapefile guide from rasterio cookbook, I wrote the code below. The code runs correctly. I would like to improve it. 1- ...
4
votes
1answer
484 views

How to create raster stack from names on R

I have a dataframe containing the best predictors by their contribution to a model: bestpredictors var contrib 2 Pfallwinter 8.1141 15 Tofall 58.0000 23 ai 5.3573 26 ...
4
votes
1answer
413 views

Where are the “band names” of an ENVI file once it is imported in R?

I have a .dat file (ENVI format) with specific band names, as you can see in the snapshot. I've imported this file in R, but I cannot find the band names in the R object. The command names(...
4
votes
1answer
1k views

R package gapfill: How to convert R raster stack to 4-dimensional array and then back to raster stack?

I'm currently analyzing MODIS EVI time series (MOD13Q1) for a small region between 2000 and 2015. Now, I would like to fill missing values using the new R package gapfill. What I have done so far: ...
3
votes
1answer
2k views

Import multi-band TIFF in R

I am trying to import a TIFF file that contains 2 bands. When using the following code in R, it seems only the first band is being recognised. S1<-"my/path/" S1<-list.files(S1, full.names = ...
3
votes
1answer
1k views

Raster calculation loop - R

I have a raster stack of 108 elements The stack 'images' is created using 4 bands (RGB+NIR) of 27 satellite images. The sequence of the rasters in the stack is: Position [[1]]: blue band image 1 ...
3
votes
2answers
5k views

Handling multiple-extent-problem to create raster stack in R?

I am storing ~ 30 raster files in a folder and want to create a raster stack. I know that not all of these rasters have completely the same extent. They do have the same resolution and CRS. Thus I ...
3
votes
2answers
3k views

Extract raster layer data by attributes in RasterStacks R

I have a RasterStack with 2 raster layers of the same area, one from a SpatialPolygonsDataFrame (layer1) and another NDVI satellite image (NDVI). layer1 has a series of polygons, each with various ...
3
votes
3answers
4k views

Stack rasters with different origin, dimensions and resolution in R

I have three raster files with the same extent and CRS, but with different origin, dimensions and resolution. r1 class : RasterLayer dimensions : 3020, 3043, 9189860 (nrow, ncol, ncell) ...
3
votes
2answers
3k views

How to remove the Na values from a stack of rasters?

In R I have a stack (nsw.stack) which contains 6 variables. I have clipped the extent. However, When I look at the summary it appears that there is a lot of NA values (see below). I think this is ...
3
votes
3answers
991 views

Unstack band layers from raster

I have stacked layer (from ERDAS) of Landsat 1-5 and 7, bands. I want each of the bands in a separate layer. How can I do that?
3
votes
1answer
174 views

Stack sat images and run operation on all

I have about 300 MODIS satellite with different acqusition dates, covering the same area. I am not used to work with so many tiles at once so I ask for advice. I am working with ArcMap10.1 What I ...
3
votes
0answers
510 views

How to calculate 90% confidence interval in a raster stack

I have a stack of grid files (23) in ASC format from a MODIS sensor and I need to calculate a 90% confidence interval, mean, min, etc, for each one in R, how can I do that? specially the interval. #...
2
votes
1answer
3k views

R: extract values from raster stack and aggregate results

I'd like to extract values from a rasterstack, cell by cell, and create a list of the values: ext <- extent(0,100,0,100) r1 <- raster(nrows=100, ncols=100,ext) r1[] <- sample(seq(from = 1, ...
2
votes
1answer
592 views

Calculating quantile values over a stack of images with grass gis. Help needed on r.series

I need to calculate the values for the 25%, 50%, 75% and 90% quantile over a series of rasters (time series of ndvi values from landsat). I am exploring the r.series command in Grass GIS, but the ...
2
votes
1answer
2k views

Stacking overlapping rasters, keeping cell values on top?

I have a land cover classification that consists of three files. One file is the LCC for the western half of my study area, another covers the eastern half (these two overlap in the middle), and the ...
2
votes
0answers
43 views

Biweekly time series for raster stack in r

I have a netcdf file with daily data and I calculate for each week some statistics, using the "apply.weekly" command in the "xts" library. However, I would like to calculate bi-weekly statistics, ...
2
votes
0answers
134 views

Layer Stack 6 Landsat8 bands

I have 6 images from Landsat 8 (Green,Red,NIR,SWIR1,SWIR2,NDVI). When I want to layerstack them with "Merge" I always get a result tif which is just black and white, however I can see that all 6 bands ...
2
votes
1answer
76 views

How to determine the appropriate size of a 3 dimesional raster block/chunk

I am processing large MODIS time series (12 years) dataset with an array dimension of 4800 x 4800 with a data type of Int16-Sixteen bit signed integer. I want to build a 3 dimensional array (row,col,...
2
votes
0answers
381 views

R raster stack 4 dimensions

I am trying to create a raster stack / brick from a list of ncdf files with the dimensions 160x160x24x1). r = stack(files) This gives me a raster stack but without the 4th dimensions which is ...
1
vote
2answers
9k views

Can ArcGIS Desktop stack layers?

I have 7 bands of Landsat and I would like to combine all layers into one layer. Can ArcMap 10 do this like Erdas software can?
1
vote
1answer
159 views

stack 10 modis ndvi geotiffs and scale them by 0.0001

I have 10 geotiffs Modis ndvi named 2005.tif ... 2014.tif I want to stack them all and create timeseries. i know how to do that, but the stacked images still need scaling by 0.0001. how can i bulk ...
1
vote
1answer
149 views

Matching Decedal Average NDVI and Worldclim Data

Ive never used NDVI (or EVI) datasets in R. Does anyone have any resource I could use to access a decadal average NDVI dataset for the globe. It would be preferable if this dataset was also comparable ...
1
vote
1answer
3k views

stack rasters with different extent and resolution in a loop with R package raster

I would like to stack rasters with the R package "raster". Layers are in ascii format. I want to crop them to have the same extent and project them to have the same resolution. My code gives raster ...
1
vote
1answer
2k views

Using GDAL/Python to stack georeferenced images of different sizes

I am in the process of porting a code I wrote in IDL (interactive data language) to Python but am running into a bit of a problem that I am hoping someone can help me with. The code goes like this: ...
1
vote
2answers
46 views

stack() in R does not include all .tiff files in my directory

I am trying to perform stack() in R using enhanced vegetation index (EVI) time-series maps but it doesn't include all the .tiff files in my directory. Here's the list of files when I dir(files.evi): ...
1
vote
1answer
40 views

Determining which layer in raster stack is first to meet condition using R?

I want to analyze a raster stack and determine the first layer in the stack that meets a condition for every pixel. Here's a reproducible code example generating rasters and making a stack. # ...
1
vote
1answer
1k views

R Stack. Extent and CRS

I'm new on R. I was listing .tif files from a folder to R. My final product is going to be a value mean of all the files in only one image. So I start stacking all of them, then brick the stack, and ...
1
vote
2answers
247 views

How to convert a raster output in a stack in R into ArcGIS

I am modelling the response of a bird (skylark) to various climatic and topographical variables in the program R. Using the predict function I predicted the abundance per pixel (see below for code). I ...
1
vote
1answer
367 views

Running Erdas command in MS-DOS

I have generated the following command line in ERDAS, and want to run it in MS-DOS. Getting an error- "system cannot find the path specified". I have checked the path by opening it in windows ...
1
vote
1answer
88 views

RasterStack extract() for duplicated dates

I have a a Large RasterStack of temperatures and a 10570 points in a Large SpatialPointsDataFrame and I'm using extract() to assign the temperatures to each point in R. sub<-stack(temp@layers[...
1
vote
0answers
20 views

Merging satellite images error

I am trying to merge 3 satellite images of 3 times. The images are median images. I have used the following code but there is an error. var maskL8 = function(image_m) { var qa = image_m.select('BQA'...
1
vote
0answers
22 views

How do I compare multiple years worth of geomorph changes in a river while drawing one profile line?

I have many years worth of DEMs of a river I'm studying. I'm very, very new to ArcGIS (think one intro class), and learned how to use the 3D analyst tool. What I'd like to be able to do, is to draw ...
1
vote
0answers
56 views

How to save statistical parameters of a raster stack in a csv using R

I have a reference raster and a stack of other rasters. I wanted to perform Mean, Std Dev and Bias of the reference with respect to the remaining. So I calculated the Cell Stats of all but when I try ...
1
vote
0answers
263 views

Is it possible to list SpatialGridDataFrame files to be used in a loop?

I have a lot (>50) of SpatialGridDataFrame files that I need to turn into raster files, stack with a set of raster files I already have and then extract data from random points from each layer of all ...
1
vote
0answers
71 views

Stacking/Aligning High-Res UAV Mosaics from Different Dates in QGIS?

I have two mosaics taken in different band sets respectively, on two different flights- of a crop field. The flights were taken a day apart, so the vegetation cover has not changed. The only issue is ...
1
vote
0answers
115 views

Determining Highest Position in QGIS for raster stacks with 10+ images?

I want to follow up on this question from last year Emulating ArcGIS Highest Position tool in QGIS? Here the question is if QGIS has a tool like the "highest position" tool in ArcGIS, where, in a ...
1
vote
0answers
221 views

Raster package limitations

I am very new to R and to raster package and I was hoping that someone could provide some guidance. I am trying to run a function that averages a series of rasters, 300 to be exact, they are in geotif ...
1
vote
0answers
352 views

gdal_merge shifts rasters after processing, layer stack is slightly off

I have experienced that gdal shifts rasters. The outputs of gdal_warp are often skewed, have different pixel sizes and so forth. A known error which has also been discussed here. But now, I had to ...
1
vote
0answers
181 views

write rasters from estudm2spixdf with significant digits

I have estimated my kernels using kernelUD from the adehabitatHS package which produces information from multiple animals as a class estUDm. I am trying to write out from R to get into ArcGIS the ...
0
votes
2answers
2k views

R -plot 12 images per plot or save 12 images per image file in R

Here is R script and my question is How can I plot 12 images per plot or save 12 images per image file because I have total 864 images. RAD1998 <- raster(paste(getwd(), "/1998bil /...
0
votes
1answer
1k views

Is there anyway to stack or merge .hdf files in R Studio?

I would like to be able to merge, mosaic, or stack several .hdf files. Does anyone know the package or library i should be looking for the tool to do this?
0
votes
1answer
131 views

Using raster stack function time in R?

I successfully ran a function over a smaller test raster stack (48 layers) in about 30 minutes but cannot run it over my large stack (1462 layers). The last one I tried over the large stack ran from ...
0
votes
1answer
37 views

Plot a series of subsetted pixels from a raster stack in R

I have created two raster stacks and would like to plot a row or column of pixels through the stack, but am not sure how to approach this. I know I can subset by layer, but how do I subset ...
0
votes
2answers
62 views

Why does the Extent of NIR Band 8 Landsat image not match other bands, even when cropping to a specified extent?

I'm trying to load and then stack four bands 2,3,4 and 8 RGB and NIR from Landsat images. I have downloaded and extracted them and they are in .tif format. # loading the bands B2 <- raster("./...