[GRASS-dev] seg fault of i.atcorr on quickbird

[Markus Neteler]

On Sun, Nov 27, 2016 at 7:47 PM, Yann <dr.yann.chemin at gmail.com> wrote:
> trying to run i.atcorr on an image of Colombo.
>
> i.atcorr --overwrite input=radxs.1 parameters=6s_b1.txt output=refxs.1
...
> Segmentation fault (core dumped)
>
> I looked into the iwave.cpp, it seems OK. What is the wavelength being less than 0.25 micron, never seen that before...

I wonder about this:

grep 'ffu.s\[-' imagery/i.atcorr/iwave.cpp
    case 1: for(i = 0; i < 65; i++)  ffu.s[-99+i] = sr1[i];
    case 2: for(i = 0; i < 65; i++)  ffu.s[-99+i] = sr2[i];
    case 3: for(i = 0; i < 65; i++)  ffu.s[-99+i] = sr3[i];
    case 4: for(i = 0; i < 65; i++)  ffu.s[-99+i] = sr4[i];
    case 5: for(i = 0; i < 65; i++)  ffu.s[-99+i] = sr5[i];

Negative value (-99)? If I remember correctly the Quickbird support
was added still without using the script.

In r69930 (just now) I have added the Quickbird2 filter functions as
CSV which I got from Nikos in August to that you may re-run the
create_iwave.py script and compare the resulting iwave.cpp etc code to
that present in the source code at time.

Markus