[GRASS-dev] [GRASS GIS] #3622: r.out.gdal export in USGSDEM format takes very long time and finally crashes

[GRASS GIS]

#3622: r.out.gdal export in USGSDEM format takes very long time and finally
crashes
-------------------------------------------------+-------------------------
 Reporter:  dido                                 |      Owner:  grass-dev@…
     Type:  defect                               |     Status:  new
 Priority:  critical                             |  Milestone:
Component:  Raster                               |    Version:  7.4.0
 Keywords:  r.out.gdal export USGSDEM long time  |        CPU:  x86-64
  crashes                                        |
 Platform:  MSWindows 7                          |
-------------------------------------------------+-------------------------
 I cannot attach any test raster because the one that manifests the problem
 is very large. Here I apply raster metadata (I suppose any raster of this
 size will do):

 {{{
 r.info map=BGMdem_8x8m_holes_filled_feathered at PERMANENT
 +----------------------------------------------------------------------------+
  | Map:      BGMdem_8x8m_holes_filled_feat  Date: Fri Aug 17 13:30:15 2018
 |
  | Mapset:   PERMANENT                      Login of Creator: a
 |
  | Location: BG_DEM
 |
  | DataBase: E:\mapping\GRASS7DB
 |
  | Title:    BGMdem_8x8m_holes_filled_feathered
 |
  | Timestamp: none
 |
 |----------------------------------------------------------------------------|
  |
 |
  |   Type of Map:  raster               Number of Categories: 0
 |
  |   Data Type:    DCELL
 |
  |   Rows:         52101
 |
  |   Columns:      77901
 |
  |   Total Cells:  4058720001
 |
  |        Projection: UTM (zone 35)
 |
  |            N: 4928368.63825695    S: 4511560.45143113   Res: 8.00000359
 |
  |            E: 663749.15806402    W: 40538.62472039   Res: 8.00003252
 |
  |   Range of data:    min = -52.070711097653  max = 2925
 |
  |
 |
  |   Data Description:
 |
  |    generated by r.mapcalc
 |
  |
 |
  |   Comments:
 |
  |    if(isnull(defekti_rast_buf_recoded), BGMdem_8x8m_holes_filled,
 |
  |    if(defekti_rast_buf_recoded < 100, BGMdem_8x8m_holes_filled * (100 -
 |
  |    defekti_rast_buf_recoded) / 100 + ALOSdem_interp *
 |
  |    defekti_rast_buf_recoded / 100, BGMdem_8x8m_holes_filled))
 |
  |
 |
 +----------------------------------------------------------------------------+
 }}}

 At crash time progress indicator has not advanced at all. Tool output
 follows:

 {{{
 r.out.gdal -f --overwrite
 input=BGMdem_8x8m_holes_filled_feathered at PERMANENT output=E:\tmp\BGM-dem-
 f5p-csovp-wbf-cas-clgeo-s0_UTM_blended_20180816.dem format=USGSDEM
 type=Int16
 Driver <USGSDEM> does not support direct writing. Using MEM driver for
 intermediate dataset.
 WARNING: Precision loss: Raster map
 <BGMdem_8x8m_holes_filled_feathered at PERMANENT> of type DCELL to be
 exported as Int16. This can be avoided by using Float64.
 WARNING: Forcing raster export
 Checking GDAL data type and nodata value...
 Using GDAL data type <Int16>
 Input raster map contains cells with NULL-value (no-data). The value
 -32768 will be used to represent no-data values in the input map. You can
 specify a nodata value with the nodata option.
 Exporting raster data to USGSDEM format...
 }}}

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