[GRASS-user] Pansharpening Landsat 8 Imagery

Michele_Toma · July 14, 2016, 12:39am

Hi everyone,

I am having difficulty with pansharpening Landsat 8 imagery using this tutorial: http://planet.qgis.org/planet/tag/landsat/. I am able to create a natural color image, but the pansharpening results in images with distorted colors.

Here are the steps that I have followed:

Import all Landsat 8 bands to GRASS 7.0.4

a. r.in.gdal input=C:\Users\Michele\Desktop\LC81100362016082LGN00\LC81100362016082LGN00_B1.TIF output=LC81100362016082LGN00_B1 -k -e

Converted digital number data to reflectance using i.landsat.toar

a. i.landsat.toar input=LC81100362016082LGN00_B output=LC81100362016082LGN00_refl metfile=C:\Users\Michele\Desktop\LC81100362016082LGN00\LC81100362016082LGN00_MTL.txt sensor=oli8 date=2016-03-22

Color balanced using i.colors.enhance for bands 4, 3, 2 as RGB respectively

a. i.colors.enhance red=LC81100362016082LGN00_refl4@Osaka green=LC81100362016082LGN00_refl3@Osaka blue=LC81100362016082LGN00_refl2@Osaka

Displayed RGB map using bands 4, 3, 2 (I get a really great natural color image)
Pansharpened bands 4, 3, 2 using band 8 as the panchromatic band with the IHS method

a. i.pansharpen red=LC81100362016082LGN00_refl4@Osaka green=LC81100362016082LGN00_refl3@Osaka blue=LC81100362016082LGN00_refl2@Osaka pan=LC81100362016082LGN00_refl8@Osaka output=LC81100362016082LGN00_ihs method=ihs

Displayed RGB map with the “pansharpened” bands for 4, 3, 2
Result: An image that is entirely gray (I found this link that says to use i.fusion.hpf, but following those instructions, I get “‘make’ is not recognized as an internal or external command, operable program or batch file.” http://gis.stackexchange.com/questions/111763/pan-sharpen-landsat-8-using-grass-7 Am I supposed to use the GUI for this?)

From searching online, I have found that the i.pansharpen module doesn’t take 16 bit imagery (which Landsat imagery is) but rather 8 bit. I found this link (http://gis.stackexchange.com/questions/39174/how-to-color-balance-false-composites-of-pan-sharpened-landsat-images-in-grass/67319#67319) that says to rescale the image. I attempted this twice:

When I run it on the bands produced after running the i.landsat.toar module, I get numbers such as min=-0.00813278402972343 and max=1.58301384751301. After rescaling the images to 0-255 and proceeding with step 3 I end up with an entirely white image.
However, if I rescale the original bands which all have ranges of 0-65535, to 0-255, then pansharpen, the result is a washed out version of the natural color image. I believe this has been pansharpened, but the colors are off. I also attempted to color balance the pansharpened bands but was unsuccessful as well.

Am I rescaling the images incorrectly or is there another step I am getting wrong? Any guidance would be greatly appreciated.

Thank you,

Michele Toma

michele.toma@airborne.aero

www.airborne.aero

Michele_Toma · July 18, 2016, 3:00pm

Hi everyone,

I was wondering if anyone had any suggestions or workarounds to the issue that I am having with pansharpening Landsat imagery? My email below has the steps that I took to attempt to pansharpen, but once I create a natural color image and run the pansharpen tool, I end up with distorted colors. Are there any other methods that I should try?

Thanks in advance,

Michele

(attachments)

···

From: grass-user [mailto:grass-user-bounces@lists.osgeo.org] On Behalf Of Michele Toma
Sent: Wednesday, July 13, 2016 5:39 PM
To: grass-user@lists.osgeo.org
Subject: [GRASS-user] Pansharpening Landsat 8 Imagery

Hi everyone,

I am having difficulty with pansharpening Landsat 8 imagery using this tutorial: http://planet.qgis.org/planet/tag/landsat/. I am able to create a natural color image, but the pansharpening results in images with distorted colors.

Here are the steps that I have followed:

Import all Landsat 8 bands to GRASS 7.0.4

a. r.in.gdal input=C:\Users\Michele\Desktop\LC81100362016082LGN00\LC81100362016082LGN00_B1.TIF output=LC81100362016082LGN00_B1 -k -e

Converted digital number data to reflectance using i.landsat.toar

a. i.landsat.toar input=LC81100362016082LGN00_B output=LC81100362016082LGN00_refl metfile=C:\Users\Michele\Desktop\LC81100362016082LGN00\LC81100362016082LGN00_MTL.txt sensor=oli8 date=2016-03-22

Color balanced using i.colors.enhance for bands 4, 3, 2 as RGB respectively

a. i.colors.enhance red=LC81100362016082LGN00_refl4@Osaka green=LC81100362016082LGN00_refl3@Osaka blue=LC81100362016082LGN00_refl2@Osaka

Displayed RGB map using bands 4, 3, 2 (I get a really great natural color image)
Pansharpened bands 4, 3, 2 using band 8 as the panchromatic band with the IHS method

a. i.pansharpen red=LC81100362016082LGN00_refl4@Osaka green=LC81100362016082LGN00_refl3@Osaka blue=LC81100362016082LGN00_refl2@Osaka pan=LC81100362016082LGN00_refl8@Osaka output=LC81100362016082LGN00_ihs method=ihs

Displayed RGB map with the “pansharpened” bands for 4, 3, 2
Result: An image that is entirely gray (I found this link that says to use i.fusion.hpf, but following those instructions, I get “‘make’ is not recognized as an internal or external command, operable program or batch file.” http://gis.stackexchange.com/questions/111763/pan-sharpen-landsat-8-using-grass-7 Am I supposed to use the GUI for this?)

From searching online, I have found that the i.pansharpen module doesn’t take 16 bit imagery (which Landsat imagery is) but rather 8 bit. I found this link (http://gis.stackexchange.com/questions/39174/how-to-color-balance-false-composites-of-pan-sharpened-landsat-images-in-grass/67319#67319) that says to rescale the image. I attempted this twice:

When I run it on the bands produced after running the i.landsat.toar module, I get numbers such as min=-0.00813278402972343 and max=1.58301384751301. After rescaling the images to 0-255 and proceeding with step 3 I end up with an entirely white image.
However, if I rescale the original bands which all have ranges of 0-65535, to 0-255, then pansharpen, the result is a washed out version of the natural color image. I believe this has been pansharpened, but the colors are off. I also attempted to color balance the pansharpened bands but was unsuccessful as well.

Am I rescaling the images incorrectly or is there another step I am getting wrong? Any guidance would be greatly appreciated.

Thank you,

Michele Toma

michele.toma@airborne.aero

www.airborne.aero

veroandreo · July 18, 2016, 7:40pm

Hello Michele,

I haven’t tried myself, but have you seen Markus’s blog from some time ago??

http://courses.neteler.org/processing-landsat-8-data-in-grass-gis-7-rgb-composites-and-pan-sharpening/

Maybe useful

Best,
Vero

(attachments)

···

2016-07-18 12:00 GMT-03:00 Michele Toma <michele.toma@airborne.aero>:

Hi everyone,

I was wondering if anyone had any suggestions or workarounds to the issue that I am having with pansharpening Landsat imagery? My email below has the steps that I took to attempt to pansharpen, but once I create a natural color image and run the pansharpen tool, I end up with distorted colors. Are there any other methods that I should try?

Thanks in advance,

Michele

From: grass-user [mailto:grass-user-bounces@lists.osgeo.org] On Behalf Of Michele Toma
Sent: Wednesday, July 13, 2016 5:39 PM
To: grass-user@lists.osgeo.org
Subject: [GRASS-user] Pansharpening Landsat 8 Imagery

Hi everyone,

I am having difficulty with pansharpening Landsat 8 imagery using this tutorial: http://planet.qgis.org/planet/tag/landsat/. I am able to create a natural color image, but the pansharpening results in images with distorted colors.

Here are the steps that I have followed:

Import all Landsat 8 bands to GRASS 7.0.4

a. r.in.gdal input=C:\Users\Michele\Desktop\LC81100362016082LGN00\LC81100362016082LGN00_B1.TIF output=LC81100362016082LGN00_B1 -k -e

Converted digital number data to reflectance using i.landsat.toar

a. i.landsat.toar input=LC81100362016082LGN00_B output=LC81100362016082LGN00_refl metfile=C:\Users\Michele\Desktop\LC81100362016082LGN00\LC81100362016082LGN00_MTL.txt sensor=oli8 date=2016-03-22

Color balanced using i.colors.enhance for bands 4, 3, 2 as RGB respectively

a. i.colors.enhance red=LC81100362016082LGN00_refl4@Osaka green=LC81100362016082LGN00_refl3@Osaka blue=LC81100362016082LGN00_refl2@Osaka

Displayed RGB map using bands 4, 3, 2 (I get a really great natural color image)

Pansharpened bands 4, 3, 2 using band 8 as the panchromatic band with the IHS method

a. i.pansharpen red=LC81100362016082LGN00_refl4@Osaka green=LC81100362016082LGN00_refl3@Osaka blue=LC81100362016082LGN00_refl2@Osaka pan=LC81100362016082LGN00_refl8@Osaka output=LC81100362016082LGN00_ihs method=ihs

Displayed RGB map with the “pansharpened” bands for 4, 3, 2

Result: An image that is entirely gray (I found this link that says to use i.fusion.hpf, but following those instructions, I get “‘make’ is not recognized as an internal or external command, operable program or batch file.” http://gis.stackexchange.com/questions/111763/pan-sharpen-landsat-8-using-grass-7 Am I supposed to use the GUI for this?)

From searching online, I have found that the i.pansharpen module doesn’t take 16 bit imagery (which Landsat imagery is) but rather 8 bit. I found this link (http://gis.stackexchange.com/questions/39174/how-to-color-balance-false-composites-of-pan-sharpened-landsat-images-in-grass/67319#67319) that says to rescale the image. I attempted this twice:

When I run it on the bands produced after running the i.landsat.toar module, I get numbers such as min=-0.00813278402972343 and max=1.58301384751301. After rescaling the images to 0-255 and proceeding with step 3 I end up with an entirely white image.

However, if I rescale the original bands which all have ranges of 0-65535, to 0-255, then pansharpen, the result is a washed out version of the natural color image. I believe this has been pansharpened, but the colors are off. I also attempted to color balance the pansharpened bands but was unsuccessful as well.

Am I rescaling the images incorrectly or is there another step I am getting wrong? Any guidance would be greatly appreciated.

Thank you,

Michele Toma

michele.toma@airborne.aero

www.airborne.aero

grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

Michele_Toma · July 18, 2016, 11:49pm

Hi Vero,

Thanks for the link. I believe the tutorial that I was originally following was the same as what you listed – just on a different site.

Paul Shapley was kind enough to share his tutorial (http://www.paulshapley.com/2014/04/processing-landsat-8-bands-in-grass-gis.html) which has worked for me. (Thanks Paul!)

For anyone that is interested, I followed Paul’s tutorial and was able to get a pansharpened image. I did end up color balancing (i.colors.enhance) as a final step because urban areas were extremely washed out. This corrected the problem and I now have an image that has been pansharpened with natural colors.

(Note: I used version 7.0.4 and installed the i.fusion.brovey extension as I believe this was replaced by i.pansharpen. To do this – Settings → Addons extensions → Install extension from addons. Then search for i.fusion.brovey in the “imagery” drop down.)

Thanks again for all the help!

-Michele

(attachments)

···

Hello Michele,

I haven’t tried myself, but have you seen Markus’s blog from some time ago??

http://courses.neteler.org/processing-landsat-8-data-in-grass-gis-7-rgb-composites-and-pan-sharpening/

Maybe useful

Best,

Vero

2016-07-18 12:00 GMT-03:00 Michele Toma <michele.toma@airborne.aero>:

Hi everyone,

I was wondering if anyone had any suggestions or workarounds to the issue that I am having with pansharpening Landsat imagery? My email below has the steps that I took to attempt to pansharpen, but once I create a natural color image and run the pansharpen tool, I end up with distorted colors. Are there any other methods that I should try?

Thanks in advance,

Michele

From: grass-user [mailto:grass-user-bounces@lists.osgeo.org] On Behalf Of Michele Toma
Sent: Wednesday, July 13, 2016 5:39 PM
To: grass-user@lists.osgeo.org
Subject: [GRASS-user] Pansharpening Landsat 8 Imagery

Hi everyone,

I am having difficulty with pansharpening Landsat 8 imagery using this tutorial: http://planet.qgis.org/planet/tag/landsat/. I am able to create a natural color image, but the pansharpening results in images with distorted colors.

Here are the steps that I have followed:

Import all Landsat 8 bands to GRASS 7.0.4

a. r.in.gdal input=C:\Users\Michele\Desktop\LC81100362016082LGN00\LC81100362016082LGN00_B1.TIF output=LC81100362016082LGN00_B1 -k -e

Converted digital number data to reflectance using i.landsat.toar

a. i.landsat.toar input=LC81100362016082LGN00_B output=LC81100362016082LGN00_refl metfile=C:\Users\Michele\Desktop\LC81100362016082LGN00\LC81100362016082LGN00_MTL.txt sensor=oli8 date=2016-03-22

Color balanced using i.colors.enhance for bands 4, 3, 2 as RGB respectively

a. i.colors.enhance red=LC81100362016082LGN00_refl4@Osaka green=LC81100362016082LGN00_refl3@Osaka blue=LC81100362016082LGN00_refl2@Osaka

Displayed RGB map using bands 4, 3, 2 (I get a really great natural color image)
Pansharpened bands 4, 3, 2 using band 8 as the panchromatic band with the IHS method

a. i.pansharpen red=LC81100362016082LGN00_refl4@Osaka green=LC81100362016082LGN00_refl3@Osaka blue=LC81100362016082LGN00_refl2@Osaka pan=LC81100362016082LGN00_refl8@Osaka output=LC81100362016082LGN00_ihs method=ihs

Displayed RGB map with the “pansharpened” bands for 4, 3, 2
Result: An image that is entirely gray (I found this link that says to use i.fusion.hpf, but following those instructions, I get “‘make’ is not recognized as an internal or external command, operable program or batch file.” http://gis.stackexchange.com/questions/111763/pan-sharpen-landsat-8-using-grass-7 Am I supposed to use the GUI for this?)

From searching online, I have found that the i.pansharpen module doesn’t take 16 bit imagery (which Landsat imagery is) but rather 8 bit. I found this link (http://gis.stackexchange.com/questions/39174/how-to-color-balance-false-composites-of-pan-sharpened-landsat-images-in-grass/67319#67319) that says to rescale the image. I attempted this twice:

When I run it on the bands produced after running the i.landsat.toar module, I get numbers such as min=-0.00813278402972343 and max=1.58301384751301. After rescaling the images to 0-255 and proceeding with step 3 I end up with an entirely white image.
However, if I rescale the original bands which all have ranges of 0-65535, to 0-255, then pansharpen, the result is a washed out version of the natural color image. I believe this has been pansharpened, but the colors are off. I also attempted to color balance the pansharpened bands but was unsuccessful as well.

Am I rescaling the images incorrectly or is there another step I am getting wrong? Any guidance would be greatly appreciated.

Thank you,

Michele Toma

michele.toma@airborne.aero

www.airborne.aero

grass-user mailing list
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user

NikosAlexandris · July 19, 2016, 7:40am

Hello Michele,

* Michele Toma:
..

on i.fusion.hpf,

7. Result: An image that is entirely gray (I found this link that
says to use i.fusion.hpf, but following those instructions, I get
"'make' is not recognized as an internal or external command, operable
program or batch file."
http://gis.stackexchange.com/questions/111763/pan-sharpen-landsat-8-using-grass-7
Am I supposed to use the GUI for this?)

You can use the module from either the GUI or the TUI, within from a
running GRASS-GIS session.

I know that people have installed and used it in Windows too. I have
never had the opportunity to test it on a Windows system myself though.

How did you try to install it?

I (still!) have to look and update (ie merge a pull request from
Panagiotis Mavrogiorgos which might fix path name issues in Windows).

Nikos

mlennert · July 19, 2016, 9:56am

I've seen that you've found a solution, but some comments anyhow:

On 18/07/16 17:00, Michele Toma wrote:

5. Pansharpened bands 4, 3, 2 using band 8 as the panchromatic
band with the IHS method

Have you seen that brovey is also available as a method in i.pansharpen ?

7. Result: An image that is entirely gray (I found this link that
says to use i.fusion.hpf, but following those instructions, I get
“'make' is not recognized as an internal or external command, operable
program or batch file.”

This means that you do not have make installed. How did you install GRASS GIS ? Maybe you need to install a specific development package ?

From searching online, I have found that the i.pansharpen module doesn’t
take 16 bit imagery (which Landsat imagery is) but rather 8 bit. I found
this link
(http://gis.stackexchange.com/questions/39174/how-to-color-balance-false-composites-of-pan-sharpened-landsat-images-in-grass/67319#67319)
that says to rescale the image. I attempted this twice:

1. When I run it on the bands produced after running the
i.landsat.toar module, I get numbers such as min=-0.00813278402972343
and max=1.58301384751301. After rescaling the images to 0-255 and
proceeding with step 3 I end up with an entirely white image.

How did you rescale ? And what does r.info give you on the rescaled maps ?

Moritz

mlennert · July 19, 2016, 9:56am

On 19/07/16 09:40, Nikos Alexandris wrote:

Hello Michele,

* Michele Toma:
..

on i.fusion.hpf,

7. Result: An image that is entirely gray (I found this link that
says to use i.fusion.hpf, but following those instructions, I get
"'make' is not recognized as an internal or external command, operable
program or batch file."
http://gis.stackexchange.com/questions/111763/pan-sharpen-landsat-8-using-grass-7
Am I supposed to use the GUI for this?)

You can use the module from either the GUI or the TUI, within from a
running GRASS-GIS session.

I know that people have installed and used it in Windows too. I have
never had the opportunity to test it on a Windows system myself though.

How did you try to install it?

I (still!) have to look and update (ie merge a pull request from
Panagiotis Mavrogiorgos which might fix path name issues in Windows).

How about putting it into the addon repository ?

Moritz

Ken_Mankoff · July 19, 2016, 9:59am

Hi Michele (and Paul),

On 2016-07-18 at 23:49, Michele Toma <michele.toma@airborne.aero> wrote:

Paul Shapley was kind enough to share his tutorial
(http://www.paulshapley.com/2014/04/processing-landsat-8-bands-in-grass-gis.html)
which has worked for me. (Thanks Paul!)

Thanks for this link! It is nice to finally have pan-sharpened images without reducing them to 8-bit images.

One change I suggest is that instead of using r.colors for histogram equalization, to use i.colors.enhance with the "-p" flag. I think this considers RGB bands together, which preserves "true color" a bit better, than histogram-equalizing the RGB bands independently of each other.

-k.

Ken_Mankoff · July 19, 2016, 10:11am

On 2016-07-19 at 09:56, Moritz Lennert <mlennert@club.worldonline.be> wrote:

Have you seen that brovey is also available as a method in
i.pansharpen ?

But i.pansharpen in GRASS 7 seems to only works on 8 bit images. This is a significant data loss compared to the full bandwidth of LS8.

-k.

neteler · July 29, 2016, 6:10am

On Tue, Jul 19, 2016 at 9:40 AM, Nikos Alexandris
<nik@nikosalexandris.net> wrote:
...

You can use the module from either the GUI or the TUI, within from a
running GRASS-GIS session.

I know that people have installed and used it in Windows too. I have
never had the opportunity to test it on a Windows system myself though.

Now it is also available for Windows (the Windows Addons are compiled
on a weekly base)

https://wingrass.fsv.cvut.cz/grass70/x86_64/addons/latest/

Hence, g.extension or the extension menu entry will fetch it for you.

Markus

hellik · July 29, 2016, 5:36pm

Markus Neteler wrote

On Tue, Jul 19, 2016 at 9:40 AM, Nikos Alexandris
<

nik@

> wrote:
...

You can use the module from either the GUI or the TUI, within from a
running GRASS-GIS session.

I know that people have installed and used it in Windows too. I have
never had the opportunity to test it on a Windows system myself though.

Now it is also available for Windows (the Windows Addons are compiled
on a weekly base)

https://wingrass.fsv.cvut.cz/grass70/x86_64/addons/latest/

Hence, g.extension or the extension menu entry will fetch it for you.

installing by g.extension in winGRASS7 works, but starting the script:

---------
C:\>i.fusion.hpf --help

Description:
Fusing high resolution Panchromatic and low resolution Multi-Spectral data
base
d on the High-Pass Filter Addition technique (Gangkofner, 2008)

Keywords:
imagery, fusion, sharpening, high pass filter, HPFA

Usage:
i.fusion.hpf.py [-l2c] pan=filename msx=filename(s)[,filename(s),...]
   suffix=suffix string [ratio=rational number] [center=string]
   [center2=string] [modulation=string] [modulation2=string]
   [trim=rational number] [--overwrite] [--help] [--verbose] [--quiet]
   [--ui]

Flags:
  -l Linearly match histogram of Pan-sharpened output to Multi-Spectral
input
  -2 2-Pass Processing (recommended) for large resolution ratio (>=5.5)
  -c Match color table of Pan-Sharpened output to Multi-Spectral input
--o Allow output files to overwrite existing files
--h Print usage summary
--v Verbose module output
--q Quiet module output
--qq Super quiet module output
--ui Force launching GUI dialog

Parameters:
          pan High resolution Panchromatic image
          msx Low resolution Multi-Spectral image(s)
       suffix Suffix for output image(s)
                 Names of Pan-Sharpened image(s) will end with this suffix
                default: hpf
        ratio Custom ratio
                 Custom ratio overriding standard calculation
                options: 1.0-10.0
       center Center cell value
                 Center cell value of the High-Pass-Filter
                options: low,mid,high
                default: low
      center2 2nd Pass center cell value
                 Center cell value for the second High-Pass-Filter (use -2
flag)

                options: low,mid,high
                default: low
   modulation Modulation level
                 Modulation level weighting the HPF image determining
crispness
                options: min,mid,max
                default: mid
  modulation2 2nd Pass modulation level (use -2 flag)
                 Modulation level weighting the second HPF image determining
cri
spness (use -2 flag)
                options: min,mid,max
                default: mid
                 min: Minimum: 0.25
                 mid: Mid: 0.35
                 max: Maximum: 0.5
         trim Trimming factor
                 Trim output border pixels by a factor of the pixel size of
the
low resolution image. A factor of 1.0 may suffice.

This application has requested the Runtime to terminate it in an unusual
way.
Please contact the application's support team for more information.
--------

terminating in an unusal way.

any idea?

-----
best regards
Helmut
--
View this message in context: http://osgeo-org.1560.x6.nabble.com/Pansharpening-Landsat-8-Imagery-tp5276266p5278737.html
Sent from the Grass - Users mailing list archive at Nabble.com.

hellik · July 29, 2016, 6:00pm

Helmut Kudrnovsky wrote

Markus Neteler wrote

On Tue, Jul 19, 2016 at 9:40 AM, Nikos Alexandris
<

nik@

> wrote:
...

You can use the module from either the GUI or the TUI, within from a
running GRASS-GIS session.

I know that people have installed and used it in Windows too. I have
never had the opportunity to test it on a Windows system myself though.

Now it is also available for Windows (the Windows Addons are compiled
on a weekly base)

https://wingrass.fsv.cvut.cz/grass70/x86_64/addons/latest/

Hence, g.extension or the extension menu entry will fetch it for you.

installing by g.extension in winGRASS7 works, but starting the script:

---------
C:\>i.fusion.hpf --help

[...]
terminating in an unusal way.

any idea?

tested on the command line in winGRASS73

-----------
g.gisenv set=DEBUG=3
---
i.fusion.hpf pan=lsat7_2000_80a
msx=lsat7_2000_30a,lsat7_2000_20a,lsat7_2000_10a suffix=my

[...]
|! Original Region restored
D1/3: G_set_program_name(): g.message
D2/3: G_file_name(): path = J:\gisdata\nc_spm_08_grass7\user1
D1/3: g.message message=

Hint, rebalancing colors (via i.colors.enhance) may improve appearance
of RG

B composites! debug=1 -i
D1/3: G_set_program_name(): g.message
D2/3: G_file_name(): path = J:\gisdata\nc_spm_08_grass7\user1

Hint, rebalancing colors (via i.colors.enhance) may improve appearance

of RGB composites!
D1/3: grass.script.core.start_command(): g.remove --q -f pattern=tmp.5792*
type=
raster
D1/3: G_set_program_name(): g.remove
D2/3: G_file_name(): path = J:\gisdata\nc_spm_08_grass7\user1
D2/3: G_file_name(): path = J:\gisdata\nc_spm_08_grass7\user1\cell
WARNING: No data base element files found

This application has requested the Runtime to terminate it in an unusual
way.
Please contact the application's support team for more information
---
g.list type=raster mapset=user1
lsat7_2000_10a
lsat7_2000_10a.my
lsat7_2000_20a
lsat7_2000_20a.my
lsat7_2000_30a
lsat7_2000_30a.my
lsat7_2000_80a
---
r.info map=lsat7_2000_10a.my@user1

|----------------------------------------------------------------------------|
|
|
| Type of Map: raster Number of Categories: 0
|
| Data Type: DCELL
|
| Rows: 950
|
| Columns: 1054
|
| Total Cells: 1001300
|
| Projection: Lambert Conformal Conic
|
| N: 228513 S: 214975.5 Res: 14.25
|
| E: 645012 W: 629992.5 Res: 14.25
|
| Range of data: min = 46.8253308406372 max = 292.357619939526
|
|
|
| Data Description:
|
| generated by r.mapcalc
|
|
|
| Comments:
|
| tmp.5792.4_msx_blnr + tmp.5792.4_pan_hpf * 0.019623261
|
| Weigthing applied: 13.976 / 178.050 * 0.250
|
|
|

+----------------------------------------------------------------------------+
-----------

so it seems the script calculates the images with higher resolution, but
terminates at some stage before finishing correctly

-----
best regards
Helmut
--
View this message in context: http://osgeo-org.1560.x6.nabble.com/Pansharpening-Landsat-8-Imagery-tp5276266p5278738.html
Sent from the Grass - Users mailing list archive at Nabble.com.