Carlos, The aerosol model is the 4th line in the parameter file. You have set the aerosol model to 15. But in the i.atcor manual, the values range from 0 to 11. Could that be the problem? i.atcor issues an unknown aerosol model warning and assumes you are trying to enter aerosol models by hand. Hence it sets your visibility to -0.1 (ground altitude) and other parameters wrong.
Cheers Daniel On Fri, Mar 18, 2016 at 4:49 PM, Carlos Grohmann <[email protected]> wrote: > 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. > > _______________________________________________ > grass-user mailing list > [email protected] > http://lists.osgeo.org/mailman/listinfo/grass-user >
_______________________________________________ grass-user mailing list [email protected] http://lists.osgeo.org/mailman/listinfo/grass-user
