[GRASS-dev] Implementation of the High Pass Filter Additive Fusion technique (i.fusion.hpf)

[Michael Barton]

Attached is a new version of i.pansharpen (i.pansharpen2) to test.

I've changed it so that the user can set the number of grey levels (defaults to 256). One thing that worries me is this warning I got. It comes from r.stat:

WARNING: Raster map <p034r032_7dp20010924_z13_80> is reading as integer
         map! Flag '-C' and/or 'nsteps' option will be ignored.

I'd been running r.stats with the -i flag to force it to read map values as integer. IHS pan sharpening requires integer input and most image data are integers anyway. But the warning suggests that if the map is read as integer, it can only have 256 levels. Not sure what to make of that. I've dropped the -i flag in r.stats. It seems to work OK, but I'm not sure whether it has problematic effects or not.

If anyone has test data at 8 bit and higher, please test and see if this solves the 256 level limit.

Michael



____________________
C. Michael Barton
Director, Center for Social Dynamics & Complexity
Professor of Anthropology, School of Human Evolution & Social Change
Arizona State University

voice:   480-965-6262 (SHESC), 480-965-8130/727-9746 (CSDC)
fax:          480-965-7671 (SHESC),  480-727-0709 (CSDC)
www: http://www.public.asu.edu/~cmbarton, http://csdc.asu.edu












On Nov 17, 2013, at 2:08 PM, Markus Neteler <neteler at osgeo.org>
 wrote:

> On Sun, Nov 17, 2013 at 8:47 PM, Michael Barton <Michael.Barton at asu.edu> wrote:
>> I do not have an image set that is >8bit and floating point for testing. But someone on the list did not too long ago.
>
> This may be of interest:
>
> http://landsat.usgs.gov/L8_12_bit.php
> "Landsat 8 products are delivered as 16-bit images."
>
> Markus

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/grass-dev/attachments/20131118/8cade6ea/attachment-0001.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: i.pansharpen2
Type: application/octet-stream
Size: 17865 bytes
Desc: i.pansharpen2
URL: <http://lists.osgeo.org/pipermail/grass-dev/attachments/20131118/8cade6ea/attachment-0001.obj>