Hi Carlos, >From the manual [1], the 0-255 rescaling is default. So I believe if you want to use some other rescaling, you have to do it manually. To be honest, I never understood what are the best options to use in the rescale parameter. I imagine the 0-255 is reminiscent from the 8-bit images. So what about higher bit images. Should we use other rescaling factors?
[1]: https://grass.osgeo.org/grass71/manuals/i.atcorr.html On Mon, Mar 21, 2016 at 11:40 AM, Carlos Grohmann <[email protected] > wrote: > Thanks Daniel! > > That was really dumb of not double-checking the input parameters.. > > Now another problem: the resulting image is empty! r.info returns "min = > nan max = nan". > > I tried with radiance and reflectance as inputs. One thing that I realized > is that, even though I'm not setting it to rescale the output range, it's > doing it: > > | Comments: > | > | i.atcorr --overwrite --verbose -r > input="ast_18560.01.toa_reflectanc\ | > | e@lateritas" range=0,255 elevation="srtm_xingu_3sec@lateritas" > param\ | > | > eters="/Volumes/MacintoshHD2/Dropbox/USP/projetosPesquisa/Xingu_late\ | > | ritas/aster/params6s_ast_18560.txt" > output="ast_18560.01.toa_reflect\ | > | ance_atcorr_6s" rescale=0,255 > > Here I didn't set neither range nor rescale (left them blank), but > apparently the command is not respecting my choices. bug detected? > > After this I tried setting the input range and rescale options, but > without success. Still an empty raster... > > best > > Carlos > > > > > > On Mon, Mar 21, 2016 at 8:54 AM, Daniel Victoria < > [email protected]> wrote: > >> 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 >>> >> >> > > > -- > 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
