[GRASS-user] unexpected i.pansharpen results

Nikos Alexandris nik at nikosalexandris.net
Sat Jul 27 05:49:09 PDT 2013


[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 at PERMANENT)"
> r.mapcalc "ik_c1r2_f.3 = float(ik_c1r2_mss.3 at PERMANENT)"
> r.mapcalc "ik_c1r2_f.2 = float(ik_c1r2_mss.2 at PERMANENT)"
> r.mapcalc "ik_c1r2_f.p = float(ik_c1r2_pan at 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 at PERMANENT
> ms2=ik_c1r2_f.3 at PERMANENTms1=ik_c1r2_f.2 at PERMANENTpan=ik_c1r2_f.p at 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


More information about the grass-user mailing list