[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 veroandreo):

 Replying to [comment:44 Nikos Alexandris]:
 > Update(d) with latest
 https://trac.osgeo.org/grass/attachment/ticket/2048/i.pansharpen.py and
 repeat for PCA. Here https://i.imgur.com/lkEoBRw.png, I think it is still
 not what we should expect color-wise. If someone else sees the same fuzzy
 colors, maybe we should go again through the steps for the PCA method.

 > {{{
 > # process
 > i.pansharpen pan=B8 red=B4 green=B5 blue=B1 bitdepth=16 output=pan_pca
 method=pca --o
 >
 > # in & out
 > for RASTER in $(echo B8 B4 B5 B1 $(glr pattern=pan_pca*)) ;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
 >
 > 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=37.4635679507703
 > max=73.8477809701631
 >
 > 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=-60.815057189378
 > max=71.127075165179
 > }}}

 Maybe a silly question, but why does pan_pca_green has a minimum of -60?
 Shouldn't that be >= 0?

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