[GRASS-dev] [GRASS GIS] #2089: i.landsat.toar doesn't take into account the resolution difference between panchromatic and the other bands.

[GRASS GIS]

#2089: i.landsat.toar doesn't take into account the resolution difference between
panchromatic and the other bands.
----------------------------+-----------------------------------------------
 Reporter:  vesnikos        |       Owner:  grass-dev@…              
     Type:  defect          |      Status:  new                      
 Priority:  major           |   Milestone:  7.0.0                    
Component:  Imagery         |     Version:  svn-trunk                
 Keywords:  i.landsat.toar  |    Platform:  All                      
      Cpu:  Unspecified     |  
----------------------------+-----------------------------------------------

Comment(by neteler):

 Replying to [ticket:2089 vesnikos]:
 > i.landsat.toar transforms the sensors DN to top of the atmosphere
 radiance/reflectance.
 >
 > All bands except band 8 (panchromatic) have a pixel size of 30m. Band 8
 has pixel size of 15.

 Did you check the input resolution? I tested it with the Raleigh scene
 LC80160352013134LGN03.tar.gz:

 {{{
 for i in `ls LC*.TIF` ; do echo -n "$i: " ; gdalinfo $i | grep Pixel ;
 done
 LC80160352013134LGN03_B10.TIF: Pixel Size =
 (30.000000000000000,-30.000000000000000)
 LC80160352013134LGN03_B11.TIF: Pixel Size =
 (30.000000000000000,-30.000000000000000)
 LC80160352013134LGN03_B1.TIF: Pixel Size =
 (30.000000000000000,-30.000000000000000)
 LC80160352013134LGN03_B2.TIF: Pixel Size =
 (30.000000000000000,-30.000000000000000)
 LC80160352013134LGN03_B3.TIF: Pixel Size =
 (30.000000000000000,-30.000000000000000)
 LC80160352013134LGN03_B4.TIF: Pixel Size =
 (30.000000000000000,-30.000000000000000)
 LC80160352013134LGN03_B5.TIF: Pixel Size =
 (30.000000000000000,-30.000000000000000)
 LC80160352013134LGN03_B6.TIF: Pixel Size =
 (30.000000000000000,-30.000000000000000)
 LC80160352013134LGN03_B7.TIF: Pixel Size =
 (30.000000000000000,-30.000000000000000)
 LC80160352013134LGN03_B8.TIF: Pixel Size =
 (15.000000000000000,-15.000000000000000)
 LC80160352013134LGN03_B9.TIF: Pixel Size =
 (30.000000000000000,-30.000000000000000)
 LC80160352013134LGN03_BQA.TIF: Pixel Size =
 (30.000000000000000,-30.000000000000000)
 }}}

 The TIFF files are delivered as such (only B8 is 15m).

 I also see that even the PAN B8 band is set to 30m in GRASS 7:
 {{{
 i.landsat.toar input_prefix=LC80160352013134LGN03_B output_prefix=toar \
   metfile=LC80160352013134LGN03_MTL.txt sensor=oli8 date=2013-05-14  -r
 ...

 for i in `g.mlist rast pattern="toar*"` ; do echo "$i: " ; r.info -g $i |
 grep -i res ; done
 toar1:
 nsres=30
 ewres=30
 toar10:
 nsres=30
 ewres=30
 toar11:
 nsres=30
 ewres=30
 toar2:
 nsres=30
 ewres=30
 toar3:
 nsres=30
 ewres=30
 toar4:
 nsres=30
 ewres=30
 toar5:
 nsres=30
 ewres=30
 toar6:
 nsres=30
 ewres=30
 toar7:
 nsres=30
 ewres=30
 toar8:
 nsres=30  <<-- should be 15
 ewres=30  <<-- should be 15
 toar9:
 nsres=30
 ewres=30
 }}}

 Using

 {{{
 g.version -g
 version=7.0.svn
 date=2013
 revision=57610
 build_date=2013-09-08
 }}}

 NOTE: In GRASS 6 i.landsat.toar preserves the resolution. Must be a bug in
 the
 usage of Rast_get_cellhd() in the G7 version?

-- 
Ticket URL: <http://trac.osgeo.org/grass/ticket/2089#comment:3>
GRASS GIS <http://grass.osgeo.org>