Possibly wrong statistics calculation

[rmorales]

Summary:
Possibly wrong statistics calculation

gnuastro version:
0.11

test image set:
https://cloud.iaa.csic.es/public.php?service=files&t=9344cdb85c640c7dd74446d8c266591d

image pixel sequence as csv file:
https://cloud.iaa.csic.es/public.php?service=files&t=56b7ade952ede68f496b313c52887aa5

command:
aststatistics 019_00.fits --hdu=0

result :
  Number of elements:                      223111
  Minimum:                                 0
  Maximum:                                 65184
  Median:                                  1295
  Mean:                                    1340.290403
  Standard deviation:                      524.5495077

question:

Using my own statistics some values are different: average, stdDev and 
max:
average = 1340.57812668077
variance = 486126.1640894976
stdDev = 697.2274837450813
mode = 0
median = 1295
min = 0
max = 65535

Writing all pixels in to a calc sheet (see download link at beginning of 
this message) :
Mean    1340.57812668077
Mode    0
Median    1295
Variance    486126.164089498
Standard Deviation    697.227483745081
Minimum    0
Maximum    65535

suggestion:
Review how the pixel values are obtained from FITS file applying the 
linear transformation with BSCALE and BZERO keywords.