[GRASS-user] i.atcorr returning empty raster and not respecting input parameters [WAS: i.atcorr - segfault on OSX]

GRASS GRASS user
1

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 <daniel.victoria@gmail.com> 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 <carlos.grohmann@gmail.com> 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 “”, line 1, in
    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
grass-user@lists.osgeo.org
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.

2

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?

···

On Mon, Mar 21, 2016 at 11:40 AM, Carlos Grohmann <carlos.grohmann@gmail.com> 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 <daniel.victoria@gmail.com> 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

–

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.

On Fri, Mar 18, 2016 at 4:49 PM, Carlos Grohmann <carlos.grohmann@gmail.com> 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 “”, line 1, in
    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
grass-user@lists.osgeo.org
http://lists.osgeo.org/mailman/listinfo/grass-user