[GRASS-dev] [GRASS GIS] #2048: i.pansharpen limited to 8-bit imagery

[GRASS GIS]

#2048: i.pansharpen limited to 8-bit imagery
-------------------------+-------------------------------------------------
  Reporter:  Nikos       |      Owner:  grass-dev@…
  Alexandris             |
      Type:  defect      |     Status:  new
  Priority:  normal      |  Milestone:  7.6.1
 Component:  Imagery     |    Version:  svn-trunk
Resolution:              |   Keywords:  i.pansharpen, sharpening, fusion,
                         |  brovey, ihs, pca, 8-bit
       CPU:  All         |   Platform:  Unspecified
-------------------------+-------------------------------------------------

Comment (by Nikos Alexandris):

 Michael,

 it works for me, i.e. 16-bit Landsat 8 bands used as input, are processed
 using the default method, and output 8-bit pan-sharpened bands.

 Here,

 {{{
 # input
 for RASTER in B8 B4 B5 B1 ;do echo $RASTER ;r.info -gr $RASTER ;echo ;done

 B8
 north=4423507.5
 south=4188892.5
 east=704707.5
 west=474292.5
 nsres=15
 ewres=15
 rows=15641
 cols=15361
 cells=240261401
 datatype=CELL
 ncats=0
 min=0
 max=65535

 B4
 north=4423515
 south=4188885
 east=704715
 west=474285
 nsres=30
 ewres=30
 rows=7821
 cols=7681
 cells=60073101
 datatype=CELL
 ncats=0
 min=0
 max=65535

 B5
 north=4423515
 south=4188885
 east=704715
 west=474285
 nsres=30
 ewres=30
 rows=7821
 cols=7681
 cells=60073101
 datatype=CELL
 ncats=0
 min=0
 max=65535

 B1
 north=4423515
 south=4188885
 east=704715
 west=474285
 nsres=30
 ewres=30
 rows=7821
 cols=7681
 cells=60073101
 datatype=CELL
 ncats=0
 min=0
 max=65535
 }}}

 and

 {{{
 #process
 i.pansharpen pan=B8 red=B4 green=B5 blue=B1 bitdepth=16 output=pan_ihs
 method=ihs --o
 i.pansharpen pan=B8 red=B4 green=B5 blue=B1 bitdepth=16 output=pan_brovey
 method=brovey --o
 i.pansharpen pan=B8 red=B4 green=B5 blue=B1 bitdepth=16 output=pan_pca
 method=pca --o
 }}}

 and
 {{{
 # output
 for RASTER in $(glr pattern=pan*) ;do echo $RASTER ;r.info -gr $RASTER
 ;echo ;done

 pan_brovey_blue
 north=4315237.5
 south=4302412.5
 east=524647.5
 west=512302.5
 nsres=15
 ewres=15
 rows=855
 cols=823
 cells=703665
 datatype=DCELL
 ncats=0
 min=0
 max=21.8425196850394

 pan_brovey_green
 north=4315237.5
 south=4302412.5
 east=524647.5
 west=512302.5
 nsres=15
 ewres=15
 rows=855
 cols=823
 cells=703665
 datatype=DCELL
 ncats=0
 min=0
 max=65.8663101604278

 pan_brovey_red
 north=4315237.5
 south=4302412.5
 east=524647.5
 west=512302.5
 nsres=15
 ewres=15
 rows=855
 cols=823
 cells=703665
 datatype=DCELL
 ncats=0
 min=0
 max=28.5573122529644

 pan_ihs_blue
 north=4315237.5
 south=4302412.5
 east=524647.5
 west=512302.5
 nsres=15
 ewres=15
 rows=855
 cols=823
 cells=703665
 datatype=CELL
 ncats=0
 min=0
 max=77

 pan_ihs_green
 north=4315237.5
 south=4302412.5
 east=524647.5
 west=512302.5
 nsres=15
 ewres=15
 rows=855
 cols=823
 cells=703665
 datatype=CELL
 ncats=0
 min=0
 max=109

 pan_ihs_red
 north=4315237.5
 south=4302412.5
 east=524647.5
 west=512302.5
 nsres=15
 ewres=15
 rows=855
 cols=823
 cells=703665
 datatype=CELL
 ncats=0
 min=0
 max=87

 pan_pca_blue
 north=4315237.5
 south=4302412.5
 east=524647.5
 west=512302.5
 nsres=15
 ewres=15
 rows=855
 cols=823
 cells=703665
 datatype=DCELL
 ncats=0
 min=31.40914130045
 max=96.916420771218

 pan_pca_green
 north=4315237.5
 south=4302412.5
 east=524647.5
 west=512302.5
 nsres=15
 ewres=15
 rows=855
 cols=823
 cells=703665
 datatype=DCELL
 ncats=0
 min=44.8984732158751
 max=71.6383695727307

 pan_pca_red
 north=4315237.5
 south=4302412.5
 east=524647.5
 west=512302.5
 nsres=15
 ewres=15
 rows=855
 cols=823
 cells=703665
 datatype=DCELL
 ncats=0
 min=-41.5285779366908
 max=33.7977309887422
 }}}

-- 
Ticket URL: <https://trac.osgeo.org/grass/ticket/2048#comment:41>
GRASS GIS <https://grass.osgeo.org>