Skip to main content
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
5 votes
Accepted

Calling gdal_pansharpen.py from another Python script?

For older versions of GDAL: Ensure gdal_pansharpen.py is in your python path. You may need to add it by setting the PYTHONPATH env var before running your script, or within your script, i.e. sys.path....
user2856's user avatar
  • 66.8k
5 votes

Pansharpening Sentinel-2 imagery in Google Earth Engine

Thanks for sharing the article. I haven't seen this before, and I'm sure I'll get some use of this at one point or another. Here's my stab at getting the code running, and cleaning up a bit: var ...
Daniel Wiell's user avatar
  • 14.3k
5 votes

Pansharpen a Landsat TM 4-5 image in QGIS

No, Landsat 4 and 5 did not carry a panchromatic sensor with 15m resolution. Sources: https://www.usgs.gov/landsat-missions/landsat-4 https://earth.esa.int/eogateway/instruments/thematic-mapper-tm-
bugmenot123's user avatar
  • 10.9k
4 votes

Blank Image - Pansharpening Sentinel-2 imagery in Google Earth Engine

You cannot use resample on a composite image. ee.Image.resample(mode) An algorithm that returns an image identical to its argument, but which uses bilinear or bicubic interpolation (rather than the ...
Noel Gorelick's user avatar
2 votes

Which preprocessing comes first for pansharpened imagery- atmospheric correction or orthorectification?

In most processing chains, orthorectification or reprojection is performed as late as possible in the workflow, to avoid altering the spatial relationships between neighbouring pixels. In your case, ...
Benjamin's user avatar
  • 911
2 votes

How to pansharpen LANDSAT 8 in GRASS?

After struggling with this issue I found a simple solution that lets you pan-sharpen in QGIS (2.18.21 is my version) using SCP. The update to SCP incorporating this method is described here https://...
Mark Thompson's user avatar
2 votes
Accepted

Landsat 8 composite image using different resolutions?

Yes. Differing resolutions will likely cause issues for you later on. I'd suggest either resampling the 30m bands to 15m using a bilinear interpolation, or even pansharpening those bands as well. ...
Mikkel Lydholm Rasmussen's user avatar
2 votes
Accepted

export grass-gis pansharpened image

Please be sure to set the computational region right prior to the export, see https://grasswiki.osgeo.org/wiki/Computational_region In your case: g.region raster=SPOT6_MS_20151113.1._hpf@PERMANENT # ...
markusN's user avatar
  • 13k
2 votes

Pansharpening using GDAL tools

I had the same question when using the gdal_pansharpen.py script. I believed GDAL developped their own algorithm... But this post (here) answered your question... GDAL pan sharpening algorithm = ...
rdmato33's user avatar
  • 159
2 votes

Google Earth Engine pan sharpening of Landsat 7 image is not producing the correct true color image

Your code is correct, but this does not work for Landsat 7 as its panchromatic band includes a range of wavelengths that is both wider and shifted longer (less blue, more infrared) than visible light. ...
jscarto's user avatar
  • 21
2 votes

Calling gdal_pansharpen.py from another Python script?

ERROR 4: 3: No such file or directory I think code should be: from osgeo.scripts.gdal_pansharpen import gdal_pansharpen gdal_pansharpen(['', '-b', '1', '-b', '2', '-b', '3', 'input_pan.tif', '...
诸葛村夫's user avatar
1 vote
Accepted

Pan sharpening Landsat 8 composites in Google Earth Engine (Python API)

rgbtohsv and hsvtorgb expect normalized inputs (0 to 1). Your simpleComposite isn't normalized. Add asFloat=True to simpleComposite to get normalized values.
Noel Gorelick's user avatar
1 vote

Sentinel-2 SWIR (band 11) pansharpening

You could run a Principal Component Analysis (PCA) on the 10m bands, take PCA 1 and then perform image fusion with the 20m band(s). Not sure about ArcGIS, but in GRASS GIS rather easily done (with i....
markusN's user avatar
  • 13k
1 vote

gdal_pansharpen.py raise module errors

Use from osgeo_utils.gdal_pansharpen import gdal_pansharpen
Irin Bush's user avatar
1 vote

Pansharpening Sentinel-2 Image Collection in Google Earth Engine

You just put what you want to do in a function that takes 1 image argument and map it over the collection. (Don't do this to the entire S2 collection; filter it down to something reasonable first): ...
Noel Gorelick's user avatar
1 vote
Accepted

How to pan sharpen Skysat Multi spectral image in Google Earth Engine

ee.Image.rgbToHsv expects a 3 band image in the range [0, 1]. The Skysat imagery is not in this range (see estimated min/max values). Use ee.Image.unitScale to rescale/normalise to a 0-1 range.
user2856's user avatar
  • 66.8k
1 vote
Accepted

Landsat 8 Imagery not Creating Composite RGB Image with Pansharpening Correctly

All the zeros in your image are throwing out your image statistics so your image appears washed out when stretched for display by your GIS. Set your NoData value to 0. Alternatively, in your GIS ...
user2856's user avatar
  • 66.8k
1 vote

Pansharpening problem in ArcGIS: failure when using TOA calibrated Landsat

TOA reflectance is a float value and the tool that you are using only works with integers. May I suggest that you use the tools available in QGis instead, which are capable of dealing with both ...
Mikkel Lydholm Rasmussen's user avatar
1 vote
Accepted

Gram-Schmidt pan-sharpening in QGIS?

To summarize the comments above at the moment the is no such functionality in QGIS. You could try to implement your own plug-in to perform this task by using this scientific paper: How to pan-sharpen ...
Comrade Che's user avatar
  • 7,251
1 vote

NITF raster (Worldview) failing to load correctly in QGIS

Seems it was a change to default settings in QGIS 3.10 for CRS (see below) which meant newly loaded assets had no projection - classed as 'invalid projected' unless you specifically set it. So I ...
Rossoe's user avatar
  • 353
1 vote

Pansharpening ALL(8) worldview 3 and worldview 2 bands

Think I solved it. If you add the WV2 data to a mosaicked dataset and use "Pansharpen and Multispectral" Processing Template it looks like it sharpens all of the bands
Emtomp's user avatar
  • 543
1 vote
Accepted

Speeding up pan-sharpening with OTB plugin for QGIS (orfeo-toolbox)

Perhaps this is a bug because execution BundleToPerfectSensor with OTB's mapla.bat took only 359 seconds. Here is my issue for this bug: https://gitlab.orfeo-toolbox.org/orfeotoolbox/qgis-otb-plugin/...
Comrade Che's user avatar
  • 7,251
1 vote

Speeding up pan-sharpening with OTB plugin for QGIS (orfeo-toolbox)

RCS method (default) is fast, but poor result. Bayes method is better, but the price to pay is the processing time. The pansharpening application requires that the xs image is resampled on the same ...
Rémi's user avatar
  • 96
1 vote

Pansharpening Landsat 8 16 bit raster images

Since you have QGIS you'll have access to Python and the various GDAL tools and scripts. As of GDAL 2.1 there is a pansharpening Python script available called gdal_pansharpen.py. If you don't have it ...
Ali's user avatar
  • 4,065
1 vote

Retaining bands when pansharpening

The Create Pan-sharpened Raster Dataset tool ONLY takes a maximum of four bands (R/G/B/IR) so I think that answers your question. How to Pan-sharpen Landsat Imagery at the Earth Matters blog is ...
Hornbydd's user avatar
  • 43.7k
1 vote

Pansharpening many L8 images in folder?

Since you already know about gdal_pansharpen, what you're asking is more of a general Python question than anything. That said, here's how you could build a lookup dict with regex: import os import ...
mikewatt's user avatar
  • 5,093
1 vote

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

First of all - unless you REALLY know what you're doing and what you're experimenting with - you can not correctly convert PAN from DN to TOA reflectance. This data is made solely for the purpose of ...
Mouad Alami's user avatar
1 vote

Pansharpening of atmospheric corrected image

PanSharpening is a process of merging high-resolution panchromatic and lower resolution multispectral imagery to create a single high-resolution color image. Google Maps and nearly every map creating ...
Oualid Fouad's user avatar
1 vote

Accessing "old" ("previous date") Pansharpened Landsat data in ArcGIS

You can download Landsat imagery directly for free from a number of sources online. The GLOVIS server (http://glovis.usgs.gov/) is excellent as you can specify what the acceptable cloud cover ...
7LeagueBoots's user avatar

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