I gave up on that pansharpening attempt. I can test some of my imagery
tomorrow and let you know if I still have issues with it.
On Jul 27, 2013 3:13 PM, "Nikos Alexandris" <[email protected]> wrote:
> [Back to this old thread!]
>
> Eric Goddard wrote:
> > Alright, I'm at home on my linux box running grass7, still with the same
> > issue. I converted the inputs for my test area to float:
> >
> > r.mapcalc "ik_c1r2_f.4 = float(ik_c1r2_mss.4@PERMANENT)"
> > r.mapcalc "ik_c1r2_f.3 = float(ik_c1r2_mss.3@PERMANENT)"
> > r.mapcalc "ik_c1r2_f.2 = float(ik_c1r2_mss.2@PERMANENT)"
> > r.mapcalc "ik_c1r2_f.p = float(ik_c1r2_pan@PERMANENT)"
> >
> > and set the color table to grey for each one. (I also remembered to
> change
> > region to reflect the higher resolution of the pan image). I then ran
> > i.pansharpen:
> > i.pansharpen sharpen=ihs ms3=ik_c1r2_f.4@PERMANENT
> > ms2=ik_c1r2_f.3@PERMANENTms1=ik_c1r2_f.2@PERMANENTpan
> =ik_c1r2_f.p@PERMANENTo
> > utput_prefix=ik_ihs_f
> >
> > running r.info on the output bands from i.pansharpen shows that their
> range
> > is still 0-38.
> >
> > The output image is attached.
> >
> > Also, there is one small bug in the output. Line 330: grass.message(_("To
> > visualize output, run: g.region -p rast=%s.red" % out)) -- the %s.red
> > should be %s_red.
> >
> > I've tried using the PCA algorithm and the output looks the same, though
> > the range is different. Later on I'll try with different imagery and see
> > what happens...
>
> Is this still an issue? I get akward outputs using i.pansharpen on IKONOS
> imagery.
>
> Nikos
>
_______________________________________________
grass-user mailing list
[email protected]
http://lists.osgeo.org/mailman/listinfo/grass-user