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 <[email protected]>
wrote:
> On Sun, Nov 17, 2013 at 8:47 PM, Michael Barton <[email protected]>
> 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
i.pansharpen2
Description: i.pansharpen2
_______________________________________________ grass-dev mailing list [email protected] http://lists.osgeo.org/mailman/listinfo/grass-dev
