[GRASS-dev] [GRASS GIS] #3690: r.stats.quantile crashing with data linked through r.external

[GRASS GIS]

#3690: r.stats.quantile crashing with data linked through r.external
--------------------------+------------------------------------------
  Reporter:  neteler      |      Owner:  grass-dev@…
      Type:  defect       |     Status:  new
  Priority:  normal       |  Milestone:  7.4.3
 Component:  Raster       |    Version:  svn-releasebranch74
Resolution:               |   Keywords:  r.external, r.stats.quantile
       CPU:  Unspecified  |   Platform:  Unspecified
--------------------------+------------------------------------------

Comment (by mmetz):

 Replying to [ticket:3690 neteler]:
 > We discovered that r.stats.quantile crashes with data linked through
 r.external (when importing the same data sets with r.in.gdal it works
 fine).
 >
 > It seems that a corruption occurs with GDAL?
 >
 >
 > {{{
 > GRASS 7.4.3svn (utm30N):  > r.external input="covermap.tif" band=1
 output=tmp15409185587112 --overwrite -o
 > Reading band 1 of 1...
 > r.external komplett. Link to raster map <tmp15409185587112> created.
 >
 > GRASS 7.4.3svn (utm30N):  > r.external input="basemap.tif" band=1
 output=tmp15409185587113 --overwrite -oÜ
 > Reading band 1 of 1...
 > r.external komplett. Link to raster map <tmp15409185587113> created.
 >
 > GRASS 7.4.3svn (utm30N):  > r.stats.quantile -p base="tmp15409185587112"
 cover="tmp15409185587113" bins="2000"
 > Computing histograms
 >  100%
 > Speicherzugriffsfehler (Speicherabzug geschrieben)
 > }}}
 >
 > ... i.e. Core dump.

 The above commands do not work for me because basemap and covermap are
 swapped and the base map must be CELL type.

 These commands

 {{{
 r.external input=covermap.tif band=1 output=covermap_ext
 r.external input=basemap.tif band=1 output=basemap_ext
 r.stats.quantile -p base=basemap_ext cover=covermap_ext bins=2000
 }}}

 work for me without segfault and valgrind does not find any memory
 violations.

 I am using GDAL 2.3.2.

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