[GRASS-user] i.atcorr - segfault on OSX

[Carlos Grohmann]

Hi all.

I'm trying to run some atmospheric correction of ASTER images. First I
calculated radiance and reflectance (with my own script, not
i.aster.toar).  I've tried with radiance and reflectance as input, but both
end up in a segfault.. the parameters file has this content:

10
6 22 13.99 -51.931 -3.355
1
15
-0.100
-1000
72


and this is the call to i.atcorr and the error:
(region is defined to the raster)

 i.atcorr --overwrite --verbose input=ast_18560.01.radiance
elevation=srtm_xingu_3sec parameters=params6s_ast_18560.txt
output=ast_18560.01.radiance_atcorr

WARNING: Unknown aerosol model!
 wavelength  less  than  0.25  micron:
 let's take s(l)=s(0.25)



* ****************************** 6s version 4.2b
****************************** *
*                        geometrical conditions identity
     *
*                        -------------------------------
     *
*                        ASTER    observation
    *
*
    *
*    month: 6 day: 22
   *
*    solar zenith angle:   34.61 deg  solar azimuthal angle:       38.12
deg   *
*    view zenith angle:     0.00 deg  view azimuthal angle:         0.00
deg   *
*    scattering angle:    145.39 deg  azimuthal angle difference:  38.12
deg   *
*
    *
*                        atmospheric model description
     *
*                        -----------------------------
     *
*            atmospheric model identity :
    *
*                tropical            (uh2o=4.12g/cm2,uo3=.247cm-atm)
     *
*
    *
*            aerosols type identity :
    *
*                 user defined aerosols model
    *
*
    *
*            optical condition identity :
    *
*                 visibility :   -0.10 km  opt. thick. 550nm :   0.0000
  *
*
    *
*                        spectral condition
    *
*                        ------------------
    *
*
    *
*                        target type
     *
*                        -----------
     *
*            homogeneous ground
    *
*              constant reflectance over the spectra     0.000
    *
*
    *
*                        target elevation description
    *
*                        ----------------------------
    *
*            ground pressure  [mb]          0.00
    *
*            ground altitude  [km]       99.990
   *
*                 gaseous content at target level:
     *
*                 uh2o=    0.000 g/cm2        uo3=    0.000 cm-atm
    *
*
    *
*                        plane simulation description
    *
*                        ----------------------------
    *
*            plane  pressure          [mb]      0.00
    *
*            plane  altitude absolute [km] 396884967424.000
   *
*                 atmosphere under plane description:
    *
*                 ozone content                0.000
    *
*                 h2o   content                0.000
    *
*                aerosol opt. thick. 550nm     0.000
    *
*
    *
*                         atmospheric correction activated
     *
*                         --------------------------------
     *
Atmospheric correction...
Segmentation fault: 11


I've tried with 7.0.2 and 6.0.4 versions, both end up in the same error.

If I try to run it from a python script I a bit more information:

grass.run_command('i.atcorr', input=ast_band, elevation='srtm_xingu_3sec',
parameters=fname, output=ast_band+'_atcorr_6s', verbose=True)
WARNING: Unknown aerosol model!
WARNING: Unsupported iwave value: 809775095
 wavelength  less  than  0.25  micron:
 let's take s(l)=s(0.25)



* ****************************** 6s version 4.2b
****************************** *
*                        geometrical conditions identity
     *
*                        -------------------------------
     *
*                        ASTER    observation
    *
*
    *
*    month: 6 day: 22
   *
*    solar zenith angle:   34.65 deg  solar azimuthal angle:       38.18
deg   *
*    view zenith angle:     0.00 deg  view azimuthal angle:         0.00
deg   *
*    scattering angle:    145.35 deg  azimuthal angle difference:  38.18
deg   *
*
    *
*                        atmospheric model description
     *
*                        -----------------------------
     *
*            atmospheric model identity :
    *
*                tropical            (uh2o=4.12g/cm2,uo3=.247cm-atm)
     *
*
    *
*            aerosols type identity :
    *
*                 user defined aerosols model
    *
*
    *
*            optical condition identity :
    *
*                 visibility :   -0.10 km  opt. thick. 550nm :   0.0000
  *
*
    *
*                        spectral condition
    *
*                        ------------------
    *
*           Traceback (most recent call last):
  File "<stdin>", line 1, in <module>
  File
"/Applications/GRASS-7.0.app/Contents/MacOS/etc/python/grass/script/core.py",
line 394, in run_command
    return handle_errors(returncode, returncode, args, kwargs)
  File
"/Applications/GRASS-7.0.app/Contents/MacOS/etc/python/grass/script/core.py",
line 312, in handle_errors
    returncode=returncode)
  File
"/Applications/GRASS-7.0.app/Contents/MacOS/etc/python/grass/exceptions/__init__.py",
line 68, in __init__
    msg = _("Module run %s %s ended with error") % (module, code)
TypeError: 'str' object is not callable


any help is appreciated

Carlos
-- 
Prof. Carlos Henrique Grohmann
Institute of Energy and Environment - Univ. of São Paulo, Brazil
- Digital Terrain Analysis | GIS | Remote Sensing -

http://carlosgrohmann.com
http://orcid.org/0000-0001-5073-5572
________________
Can’t stop the signal.
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/grass-user/attachments/20160318/e2eb6617/attachment.html>