<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">
</head>
<body>
<div class="BodyFragment"><font size="2"><span style="font-size:10pt;">
<div class="PlainText">Attached is a new version of i.pansharpen (i.pansharpen2) to test.
<br>
<br>
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:<br>
<br>
WARNING: Raster map <p034r032_7dp20010924_z13_80> is reading as integer<br>
map! Flag '-C' and/or 'nsteps' option will be ignored.<br>
<br>
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.<br>
<br>
If anyone has test data at 8 bit and higher, please test and see if this solves the 256 level limit.<br>
<br>
Michael<br>
<br>
</div>
</span></font></div>
<div class="BodyFragment"><font size="2"><span style="font-size:10pt;">
<div class="PlainText"><br>
<br>
____________________<br>
C. Michael Barton<br>
Director, Center for Social Dynamics & Complexity <br>
Professor of Anthropology, School of Human Evolution & Social Change<br>
Arizona State University<br>
<br>
voice: 480-965-6262 (SHESC), 480-965-8130/727-9746 (CSDC)<br>
fax: 480-965-7671 (SHESC), 480-727-0709 (CSDC)<br>
www: <a href="http://www.public.asu.edu/~cmbarton">http://www.public.asu.edu/~cmbarton</a>,
<a href="http://csdc.asu.edu">http://csdc.asu.edu</a><br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
<br>
On Nov 17, 2013, at 2:08 PM, Markus Neteler <neteler@osgeo.org><br>
wrote:<br>
<br>
> On Sun, Nov 17, 2013 at 8:47 PM, Michael Barton <Michael.Barton@asu.edu> wrote:<br>
>> 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.<br>
> <br>
> This may be of interest:<br>
> <br>
> <a href="http://landsat.usgs.gov/L8_12_bit.php">http://landsat.usgs.gov/L8_12_bit.php</a><br>
> "Landsat 8 products are delivered as 16-bit images."<br>
> <br>
> Markus<br>
<br>
</div>
</span></font></div>
</body>
</html>