<div dir="ltr"><div>Hi all. </div><div><br></div><div>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:</div><div><br></div><div><div>10</div><div>6 22 13.99 -51.931 -3.355</div><div>1</div><div>15</div><div>-0.100</div><div>-1000</div><div>72</div></div><div><br></div><div><br></div><div>and this is the call to i.atcorr and the error:</div><div>(region is defined to the raster)</div><div><br></div><div> i.atcorr --overwrite --verbose input=ast_18560.01.radiance elevation=srtm_xingu_3sec parameters=params6s_ast_18560.txt output=ast_18560.01.radiance_atcorr</div><div><br></div><div>WARNING: Unknown aerosol model!</div><div> wavelength less than 0.25 micron: </div><div> let's take s(l)=s(0.25)</div><div><br></div><div><br></div><div><br></div><div>* ****************************** 6s version 4.2b ****************************** *</div><div>* geometrical conditions identity *</div><div>* ------------------------------- *</div><div>* ASTER observation *</div><div>* *</div><div>* month: 6 day: 22 *</div><div>* solar zenith angle: 34.61 deg solar azimuthal angle: 38.12 deg *</div><div>* view zenith angle: 0.00 deg view azimuthal angle: 0.00 deg *</div><div>* scattering angle: 145.39 deg azimuthal angle difference: 38.12 deg *</div><div>* *</div><div>* atmospheric model description *</div><div>* ----------------------------- *</div><div>* atmospheric model identity : *</div><div>* tropical (uh2o=4.12g/cm2,uo3=.247cm-atm) *</div><div>* *</div><div>* aerosols type identity : *</div><div>* user defined aerosols model *</div><div>* *</div><div>* optical condition identity : *</div><div>* visibility : -0.10 km opt. thick. 550nm : 0.0000 *</div><div>* *</div><div>* spectral condition *</div><div>* ------------------ *</div><div>* *</div><div>* target type *</div><div>* ----------- *</div><div>* homogeneous ground *</div><div>* constant reflectance over the spectra 0.000 *</div><div>* *</div><div>* target elevation description *</div><div>* ---------------------------- *</div><div>* ground pressure [mb] 0.00 *</div><div>* ground altitude [km] 99.990 *</div><div>* gaseous content at target level: *</div><div>* uh2o= 0.000 g/cm2 uo3= 0.000 cm-atm *</div><div>* *</div><div>* plane simulation description *</div><div>* ---------------------------- *</div><div>* plane pressure [mb] 0.00 *</div><div>* plane altitude absolute [km] 396884967424.000 *</div><div>* atmosphere under plane description: *</div><div>* ozone content 0.000 *</div><div>* h2o content 0.000 *</div><div>* aerosol opt. thick. 550nm 0.000 *</div><div>* *</div><div>* atmospheric correction activated *</div><div>* -------------------------------- *</div><div>Atmospheric correction...</div><div>Segmentation fault: 11</div><div><br></div><div><br></div><div>I've tried with 7.0.2 and 6.0.4 versions, both end up in the same error.</div><div><br></div><div>If I try to run it from a python script I a bit more information:</div><div><br></div><div><div>grass.run_command('i.atcorr', input=ast_band, elevation='srtm_xingu_3sec', parameters=fname, output=ast_band+'_atcorr_6s', verbose=True)</div><div>WARNING: Unknown aerosol model!</div><div>WARNING: Unsupported iwave value: 809775095</div><div> wavelength less than 0.25 micron: </div><div> let's take s(l)=s(0.25)</div><div><br></div><div><br></div><div><br></div><div>* ****************************** 6s version 4.2b ****************************** *</div><div>* geometrical conditions identity *</div><div>* ------------------------------- *</div><div>* ASTER observation *</div><div>* *</div><div>* month: 6 day: 22 *</div><div>* solar zenith angle: 34.65 deg solar azimuthal angle: 38.18 deg *</div><div>* view zenith angle: 0.00 deg view azimuthal angle: 0.00 deg *</div><div>* scattering angle: 145.35 deg azimuthal angle difference: 38.18 deg *</div><div>* *</div><div>* atmospheric model description *</div><div>* ----------------------------- *</div><div>* atmospheric model identity : *</div><div>* tropical (uh2o=4.12g/cm2,uo3=.247cm-atm) *</div><div>* *</div><div>* aerosols type identity : *</div><div>* user defined aerosols model *</div><div>* *</div><div>* optical condition identity : *</div><div>* visibility : -0.10 km opt. thick. 550nm : 0.0000 *</div><div>* *</div><div>* spectral condition *</div><div>* ------------------ *</div><div>* Traceback (most recent call last):</div><div> File "<stdin>", line 1, in <module></div><div> File "/Applications/GRASS-7.0.app/Contents/MacOS/etc/python/grass/script/core.py", line 394, in run_command</div><div> return handle_errors(returncode, returncode, args, kwargs)</div><div> File "/Applications/GRASS-7.0.app/Contents/MacOS/etc/python/grass/script/core.py", line 312, in handle_errors</div><div> returncode=returncode)</div><div> File "/Applications/GRASS-7.0.app/Contents/MacOS/etc/python/grass/exceptions/__init__.py", line 68, in __init__</div><div> msg = _("Module run %s %s ended with error") % (module, code)</div><div>TypeError: 'str' object is not callable</div></div><div><br></div><div><br></div><div>any help is appreciated</div><div><br></div><div>Carlos</div>-- <br><div class="gmail_signature"><div dir="ltr">Prof. Carlos Henrique Grohmann<br>Institute of Energy and Environment - Univ. of São Paulo, Brazil<div>- Digital Terrain Analysis | GIS | Remote Sensing - </div><div><br></div><div><a href="http://carlosgrohmann.com/" target="_blank">http://carlosgrohmann.com</a></div><div><a href="http://orcid.org/0000-0001-5073-5572" style="font-size:13px;color:rgb(17,85,204);font-family:arial,sans-serif;background-color:rgb(255,255,255)" target="_blank">http://orcid.org/0000-0001-5073-5572</a><br><div>________________<br>Can’t stop the signal.</div></div></div></div>
</div>