[GRASS-dev] [GRASS GIS] #1689: weird import with r.in.gdal

[GRASS GIS]

#1689: weird import with r.in.gdal
-----------------------+----------------------------------------------------
 Reporter:  madi       |       Owner:  grass-dev@…              
     Type:  defect     |      Status:  new                      
 Priority:  normal     |   Milestone:  6.4.3                    
Component:  Raster     |     Version:  svn-releasebranch64      
 Keywords:  r.in.gdal  |    Platform:  Linux                    
      Cpu:  x86-64     |  
-----------------------+----------------------------------------------------
 Hi all,

 I pointed this out already in ML, but in the same thread there were many
 topics so that i better separate this one in a dedicated ticket.
 I have generated a vrt file by this procedure
 http://grass.osgeo.org/wiki/Large_raster_data_processing and imported it
 into grass with the following commands:

 r.in.gdal input="/forest/ASTER-GDEM/Dem_lzw/mosaic.vrt" output="mosaic"

 and

 r.external -r input="/forest/ASTER-GDEM/Dem_lzw/mosaic.vrt"
 output="mosaic2"

 Then, comparing mosaic and mosaic2, i get:

 GRASS 6.4.3svn (ASTER_GDEM):~ > r.info map=mosaic
 +----------------------------------------------------------------------------+
  | Layer:    mosaic                         Date: Tue Jul 17 17:57:31 2012
 |
  | Mapset:   PERMANENT                      Login of Creator: leomarg
 |
  | Location: ASTER_GDEM
 |
  | DataBase: /forest/GRASS/grassdata
 |
  | Title:     ( mosaic )
 |
  | Timestamp: none
 |
 |----------------------------------------------------------------------------|
  |
 |
  |   Type of Map:  raster               Number of Categories: 6518
 |
  |   Data Type:    CELL
 |
  |   Rows:         140401
 |
  |   Columns:      198001
 |
  |   Total Cells:  27799538401
 |
  |        Projection: Latitude-Longitude
 |
  |            N: 72:00:00.5N    S: 32:59:59.5N   Res: 0:00:01
 |
  |            E: 44:00:00.5E    W: 11:00:00.5W   Res: 0:00:01
 |
  |   Range of data:    min = -9999  max = 6518
 |
  |
 |
  |   Data Description:
 |
  |    generated by r.in.gdal
 |
  |
 |
  |   Comments:
 |
  |    r.in.gdal input="/forest/ASTER-GDEM/Dem_lzw/mosaic.vrt"
 output="mosa\   |
  |    ic"
 |
  |
 |
 +----------------------------------------------------------------------------+

 GRASS 6.4.3svn (ASTER_GDEM):~ > r.info map=mosaic2
 +----------------------------------------------------------------------------+
  | Layer:    mosaic2                        Date: Wed Jul 18 10:54:37 2012
 |
  | Mapset:   PERMANENT                      Login of Creator: leomarg
 |
  | Location: ASTER_GDEM
 |
  | DataBase: /forest/GRASS/grassdata
 |
  | Title:     ( mosaic2 )
 |
  | Timestamp: none
 |
 |----------------------------------------------------------------------------|
  |
 |
  |   Type of Map:  raster               Number of Categories: 0
 |
  |   Data Type:    CELL
 |
  |   Rows:         140401
 |
  |   Columns:      198001
 |
  |   Total Cells:  27799538401
 |
  |        Projection: Latitude-Longitude
 |
  |            N: 72:00:00.5N    S: 32:59:59.5N   Res: 0:00:01
 |
  |            E: 44:00:00.5E    W: 11:00:00.5W   Res: 0:00:01
 |
  |   Range of data:    min = -9999  max = 6518
 |
  |
 |
  |   Data Description:
 |
  |    generated by r.external
 |
  |
 |
  |   Comments:
 |
  |    r.external -r input="/forest/ASTER-GDEM/Dem_lzw/mosaic.vrt"
 output="\   |
  |    mosaic2"
 |
  |
 |
 +----------------------------------------------------------------------------+

 The problem is that with r.in.gdal, the map is shrinked, while the map
 produced by r.external is correct. See the attached pictures.

-- 
Ticket URL: <https://trac.osgeo.org/grass/ticket/1689>
GRASS GIS <http://grass.osgeo.org>