<div dir="ltr"><div>Hi, <br></div><div><br></div><div>I had the same problem last time I used i.fusion.hpf for pansharpening a SPOT image (6m to 1.5m). Question is in my view if the resulting values mean anything (i.e., they can be used for estimation of spectral indices) or it is just for visual purposes? Sorry for my ignorance, but I haven't found a clear-cut answer to this... <br></div><div><br></div><div>On the other hand, i.pansharpen does include a re-scaling step, but whatever bit depth is passed as input it will be re-scaled to 8-bit. If the input data is >8-bit, it is a bit of a loss of radiometric resolution (if the resulting values can be used for indices estimation)...</div><div><br></div><div>I'd appreciate your insights</div><div><br></div><div>Cheers,</div><div>Vero<br></div><div><br></div><div><br></div><div><br></div><div><br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">El lun., 28 oct. 2019 a las 17:40, Markus Neteler (<<a href="mailto:neteler@osgeo.org">neteler@osgeo.org</a>>) escribió:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">On Mon, Oct 28, 2019 at 4:44 PM Moritz Lennert<br>
<<a href="mailto:mlennert@club.worldonline.be" target="_blank">mlennert@club.worldonline.be</a>> wrote:<br>
> On 28/10/19 12:33, Markus Neteler wrote:<br>
> > On Mon, Oct 28, 2019 at 10:46 AM Moritz Lennert<br>
> > <<a href="mailto:mlennert@club.worldonline.be" target="_blank">mlennert@club.worldonline.be</a>> wrote:<br>
> >> On 27/10/19 18:52, Markus Neteler wrote:<br>
> >>> Hi,<br>
> >>><br>
> >>> I am using i.fusion.hpf on Worldview-2 data and see that the resulting<br>
> >>> range is completely different:<br>
> >>><br>
> >>> # orig<br>
> >>> <a href="http://r.info" rel="noreferrer" target="_blank">r.info</a> -r wv2_17OCT08182034_M2AS_058891334010_01_P001.5<br>
> >>> min=0<br>
> >>> max=1509<br>
> >>><br>
> >>> # HPF<br>
> >>> <a href="http://r.info" rel="noreferrer" target="_blank">r.info</a> -r wv2_17OCT08182034_M2AS_058891334010_01_P001.5.hpf<br>
> >>> min=-1045.79783524995<br>
> >>> max=3416.06793724612<br>
> >><br>
> >> Just guessing here, but did you try the -c flag ?<br>
> ><br>
> > Now yes:<br>
> ><br>
> > i.fusion.hpf pan=$PAN msx=$(g.list raster pattern="$COLORPREFIX.?"<br>
> > sep=comma) -c --o<br>
> ><br>
> > # band 8<br>
> > <a href="http://r.info" rel="noreferrer" target="_blank">r.info</a> -r wv2_17OCT08182034_M2AS_058891334010_01_P001.8<br>
> > min=0<br>
> > max=1605<br>
> ><br>
> > <a href="http://r.info" rel="noreferrer" target="_blank">r.info</a> -r wv2_17OCT08182034_M2AS_058891334010_01_P001.8.hpf<br>
> > min=-1163.98434592277<br>
> > max=3879.99538075697<br>
> ><br>
> > So, -c doesn't seem to have any effect here.<br>
> ><br>
><br>
> Sorry, wrong flag, should have been -l ("Linearly match histogram of<br>
> Pan-sharpened output to Multi-Spectral input")<br>
<br>
This comes with -l:<br>
<br>
<a href="http://r.info" rel="noreferrer" target="_blank">r.info</a> -r wv2_17OCT08182034_M2AS_058891334010_01_P001.8.hpf<br>
min=-932.700879728634<br>
max=3456.94617091076<br>
<br>
> And if you rescale the output manually to the input value range, do you<br>
> still see a good fusion result ? If yes, then this could be added as a<br>
> last step to the module.<br>
<br>
In lack of time I cannot test that right now.<br>
<br>
Nikos has a subset of my data, maybe he's able to check with that.<br>
<br>
Markus<br>
_______________________________________________<br>
grass-dev mailing list<br>
<a href="mailto:grass-dev@lists.osgeo.org" target="_blank">grass-dev@lists.osgeo.org</a><br>
<a href="https://lists.osgeo.org/mailman/listinfo/grass-dev" rel="noreferrer" target="_blank">https://lists.osgeo.org/mailman/listinfo/grass-dev</a></blockquote></div>