On 27/11/2016 20:06, Markus Neteler wrote:
On Sun, Nov 27, 2016 at 7:47 PM, Yann <[email protected]> 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
_______________________________________________
grass-dev mailing list
[email protected]
http://lists.osgeo.org/mailman/listinfo/grass-dev
