32 votes
Accepted

Splitting .tif image into several tiles?

I propose two solutions: the first one using QGIS, the second one using Python (GDAL). Solution using QGIS In QGIS you may create a VRT mosaic. Please follow this procedure (see the image below): ...
mgri's user avatar
  • 16.1k
14 votes
Accepted

completing satellite imagery tiles when data coverage is less than 100%

For images of the same location but different dates, I would rather talk about compositing than mosaicing (which combines images from different extents into a larger image). You will find a lot of ...
radouxju's user avatar
  • 49.5k
11 votes
Accepted

"Copy" Image coordinates to another image that is nd.array

Save result array (rf_array) as in the following lines: rf_array = df_join['Prediction'].values # returns numpy.ndarray rf_array = rf_array.reshape(img.shape[0], img.shape[1]) # rf_array = rf_array....
Kadir Şahbaz's user avatar
11 votes

Applying HDR (high dynamic range) rendering in QGIS

First of all, what you are trying to achieve is not HDR. To attain (real) HDR one needs to have excess bit depth, and to discard in a meaningful fashion some bits. For example, if you have a 3 channel ...
RafDouglas C. Tommasi's user avatar
10 votes

Is it scientifically correct to pansharpen landsat reflectance product with pan band

Fundamentally the question here is "what does 'scientifically valid' mean". If you are looking to do spectral modelling on the data, then the answer is possibly different than if you are looking at ...
Mikkel Lydholm Rasmussen's user avatar
10 votes

Counting rocks of different sizes on shoreline using ArcGIS Desktop?

For the most part not an ArcGIS answer but you could try it anyway since it is completely free software. You could try using scikit-image. You will get it if you install Anaconda (with Anaconda you ...
BERA's user avatar
  • 71.2k
10 votes
Accepted

Should I always choose Sentinel-2 atmospheric-corrected imagery?

Atmospheric correction is particularly usefull if you need consistency between dates and/or adjacent tiles. Some studies showed that the use of atmospherically corrected images improve the subsequent ...
radouxju's user avatar
  • 49.5k
7 votes

Calibrating Landsat Level-1 Precision and Terrain (L1TP) corrected data?

The three sensors are all slightly different. However the OLI/TIRs setup is a marked departure from the TM/ETM+ sensors. The changes are succintly summarised by Li et al. 2013 as the: replacing of ...
RoperMaps's user avatar
  • 2,166
6 votes

Splitting .tif image into several tiles?

Don't split You gdalbuildvrt, you can create virtual tiles that will only use a few bytes on your disk. Then you can use most softwares that will take your vrt's as input to perform your processing. ...
radouxju's user avatar
  • 49.5k
6 votes

Splitting .tif image into several tiles?

Another solution using GDAL is the gdal_retile.py tool: mkdir image_tiles gdal_retile.py -v -r bilinear -levels 1 -ps 20000 20000 -co "TILED=YES" -co "COMPRESS=JPEG" -targetDir image_tiles ...
Saleika's user avatar
  • 1,320
6 votes
Accepted

Create pandas DataFrame from raster image - one row per pixel with bands as columns

Quick solution pd.DataFrame(array.reshape([3,-1]).T) Explanation Take array of shape (3, x, y) and flatten out the 2nd and 3rd dimension. From the numpy docs: One shape dimension can be -1. In this ...
StefanBrand_EOX's user avatar
6 votes
Accepted

What is lost when converting 12 bit imagery to 8 bit?

I know that you have tried to frame this as a general question, rather than specific to SPOT 6/7, but it's really a little of both. The naive answer to the "general" question "what is ...
rs_burner's user avatar
5 votes
Accepted

Filling in missing values in raster using R

I am the author of the R package gapfill, which is a flexible tool to predict missing values in spatio-temporal remote sensing data sets. https://CRAN.R-project.org/package=gapfill It could be helpful ...
Nairolf's user avatar
  • 211
5 votes

Pansharpen landsat mosaic in Google Earth Engine

You can mosaic a collection of pan-sharpened images as follows: // Function to mask clouds using the quality band of Landsat 8. var maskL8 = function(image) { var qa = image.select('BQA'); /// ...
Tyler Erickson's user avatar
5 votes

What does cumulative count cut actually do?

I think, it does the same than numpy.nanpercentile. It takes the 2nd and the 98th percentile of your grey values and sets it as the colour-map minimum and maximum. Like colour stretch it doesn't alter ...
Rebecca's user avatar
  • 51
5 votes
Accepted

Applying HDR (high dynamic range) rendering in QGIS

You can use Orfeo Toolbox, an open source libray specialized in image processing for remote sensing images. It can be used in QGIS - see here how to make it work. For QGIS versions older than QGIS 3.8,...
Babel's user avatar
  • 70.4k
4 votes

How to create a mosaic in QGIS with cutline and feathering for Landsat-8 imagery

I am a developer of the open-source GIS Whitebox GAT, which contains several image processing tools including a tool called Mosaic With Feathering. I doubt that it is as sophisticated as ERDAS or ENVI ...
WhiteboxDev's user avatar
4 votes

Calculate NDVI image from R and NIR channels and rescale it to 0-255

the general equation for rescaling is (pixel_value - min_value)/(max_pixel_value - min_pixel_value)*(max_rescale_value-min_rescale_value)+min_rescale_value in your example, this would be (...
radouxju's user avatar
  • 49.5k
4 votes

completing satellite imagery tiles when data coverage is less than 100%

I think that what you describe is still part of what is called mosaicing (or image stitching). Mosaicing involves indeed joining adjacent tiles, but usually tiles do have some overlap. Here you are ...
Matifou's user avatar
  • 2,001
4 votes

R: Display single band raster in true color

I'm a little late to the party, but figured I could contribute a more complete answer to the question than the bits left in comments here. raster::raster() only reads a single band when handed a multi-...
Mike Mahoney's user avatar
4 votes
Accepted

Is image registration necessary among temporal Landsat images?

Using Collection 1 is not necessary: In 2016, the USGS reorganized the Landsat archive into a tiered collection. This structure ensures that Landsat Level-1 products provide a consistent ...
aldo_tapia's user avatar
  • 13.4k
4 votes

Pansharpen landsat mosaic in Google Earth Engine

Here's a variant that you may find useful/different (continuing Tyler's example): var pansharpen1 = function(scene, kernel) { // Compute the per-pixel means of the unsharpened bands. var bgr = ...
Nicholas Clinton's user avatar
4 votes

Reading large image using GDAL?

To read an image using gdal as a numpy array you can use the following two lines: ds = gdal.Open('path_to_img', 0) arr = ds.ReadAsArray() However, you can run into a MemoryError if your image is too ...
Marcelo Villa's user avatar
4 votes

Create pandas DataFrame from raster image - one row per pixel with bands as columns

Or another simple solution with numpy ravel(): import rasterio as rio src= rio.open('myraster.tif') # number of bands src.count 3 # read bands array = src.read() # convert to a DataFrame import ...
gene's user avatar
  • 54.7k
4 votes

"Copy" Image coordinates to another image that is nd.array

You must define 2 elements in order to have a geolocated image The first is the geotransform that converts the row/column coordinates into X/Y coordinates. In your case this will be done using ...
radouxju's user avatar
  • 49.5k
4 votes

Removing dropped line using fourier transformation?

It is worth noting that the strip (dropped line) you are observing is likely an error created during the collection process and is not likely to have happened down stream. It is worth checking though. ...
AWGIS's user avatar
  • 3,210
3 votes

Auto digitizing features on map?

You can use ArcScan extension of ArcGIS for this. This is pretty much easy for auto digitizing such simple map with ArcScan For more information you can check http://desktop.arcgis.com/en/arcmap/...
Ahsan Mukhtar's user avatar
3 votes

Pattern search in satellite imagery

This isn't really an answer, but it is funny that I just saw an article today about this problem. Maybe it will lead you in the right direction. http://www.nextgov.com/big-data/2016/05/promise-...
LanderTome's user avatar
3 votes
Accepted

Make Image Pyramid from JP2 Raster

please, follow iant suggestion about the jpep compressed (big)tiff with jpeg compressed overviews. I would not go for serving jp2 directly unless you can use kakadu or ecw sdk. I have not had time to ...
simogeo's user avatar
  • 1,441

Only top scored, non community-wiki answers of a minimum length are eligible