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

[Yann]

On 27/11/2016 20:06, Markus Neteler wrote:
> 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
Thanks Markus, will look into that.

Also I believe the "less than 0.25" thing comes from either (or both) of 
these:
   static const float wli[5] = {0.001, 0.001, 0.001, 0.001, 0.001};
   static const float wls[5] = {0.161, 0.161, 0.161, 0.161, 0.161};

line 2942 in iwave.cpp