..
Yann:
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
Nikos:
Back in August, I thoroughly checked running the (modified) python script
against each and every sensor's spectral response data (in form of CSVs)
we have in our repo.
I'll give it a quick check -- I can only work on it after December the
1st though.
The uploaded file is wrong in term of units. Please see the attached,
here-in, spectral response values (+ the output template) set in the
correct unit. It works fine for me. I think we had uploaded, back in
August, the "correct" one.
Best, Nikos
nm,pan,blue,green,red,nir
300,0,0,0,0,0
305,0,0,0,0,0
310,0,0,0,0,0
315,0,0,0,0,0
320,0,0,0,0,0
325,0,0.0010102,0,0,0.0014092
330,0,0.002051,0,0,0.0028608
335,0,0.0031232,0,0,0.0043557
340,0,0.0042264,0,0,0.005895
345,0,0.0045373,0,0,0.006329
350,0.0014033,0.0053218,0.0030856,0.0037148,0.0060022
355,0.001179,0.004057,0.0040267,0.002443,0.0053844
360,0.0011042,0.0032654,0.0021878,0.0016479,0.0035169
365,0.0009189,0.0027029,0.0023189,0.0019254,0.0021084
370,0.0011777,0.0023032,0.0014787,0.001181,0.0024355
375,0.0020335,0.0012599,0.0014276,0.0010728,0.0017002
380,0.0068664,0.0018015,0.0008263,0.0009032,0.0022714
385,0.0108822,0.0022036,0.0011403,0.0009703,0.0018699
390,0.0167311,0.0023949,0.0012422,0.0009224,0.0016207
395,0.0259192,0.0025521,0.0011512,0.0013128,0.0019495
400,0.0361999,0.0024719,0.0012394,0.000677,0.001831
405,0.0506315,0.0029075,0.0011408,0.0012157,0.0015098
410,0.0669468,0.0040027,0.0011869,0.0012567,0.0014986
415,0.0845848,0.0065597,0.0013912,0.0013207,0.0019432
420,0.1031961,0.0131362,0.0018941,0.0011522,0.0017185
425,0.121867,0.0283802,0.0023351,0.0009955,0.0024463
430,0.1449896,0.068135,0.0028135,0.001519,0.0025377
435,0.1683812,0.1551567,0.0037063,0.0020873,0.003128
440,0.1895051,0.2775976,0.0055455,0.0028027,0.00405
445,0.21272,0.4312732,0.0098195,0.0030937,0.0047331
450,0.2198261,0.6299497,0.0178466,0.0028875,0.0047484
455,0.2585699,0.7240532,0.0234254,0.0028213,0.0049004
460,0.2752517,0.791362,0.0300974,0.0029808,0.0041819
465,0.2955355,0.8438716,0.0411627,0.0034599,0.0044438
470,0.3153968,0.8895741,0.0596683,0.0036881,0.0047788
475,0.3364159,0.9177772,0.0981647,0.0031561,0.0052347
480,0.3456366,0.9293006,0.1566285,0.0020766,0.0055566
485,0.3546477,0.9601538,0.228192,0.0021746,0.0052475
490,0.3611199,1,0.3196465,0.0023823,0.0041908
495,0.3684864,0.9916359,0.4174751,0.0026856,0.0029904
500,0.3781246,0.9178985,0.5276728,0.0030709,0.0026452
505,0.3969294,0.7724284,0.6578772,0.0032914,0.0025776
510,0.4141351,0.5734303,0.7697888,0.0035025,0.0025379
515,0.4410028,0.3840363,0.8311977,0.0037279,0.0028625
520,0.4676101,0.2653482,0.8444885,0.0045702,0.0025801
525,0.4981964,0.1591609,0.8704719,0.0053096,0.0033657
530,0.5351925,0.0816323,0.9081774,0.0059502,0.0052825
535,0.5711692,0.0692367,0.9459902,0.006323,0.0078701
540,0.6137484,0.0633006,0.98194,0.0065926,0.0098718
545,0.6627288,0.0469325,1,0.0069161,0.0102328
550,0.7050662,0.0293462,0.9933748,0.0077701,0.0078634
555,0.7499391,0.0176386,0.9818187,0.0089839,0.0048711
560,0.7909615,0.0127288,0.9691108,0.0102276,0.0048416
565,0.8301127,0.0124551,0.9774484,0.0117892,0.0063449
570,0.8628948,0.0128916,0.9868932,0.0138696,0.0076518
575,0.891069,0.0126216,0.9652411,0.0176781,0.0079247
580,0.9250918,0.0111146,0.8979307,0.0237713,0.0064695
585,0.9334622,0.0100924,0.7909738,0.0334155,0.0043683
590,0.9518124,0.0095605,0.6550701,0.0464361,0.0037961
595,0.9620912,0.0092837,0.4822322,0.0630683,0.00374
600,0.9598312,0.0097672,0.3124203,0.0871887,0.0035986
605,0.990737,0.010639,0.2033228,0.1281145,0.0036339
610,0.9974814,0.0103876,0.1378029,0.2050211,0.0039231
615,0.9869457,0.0078934,0.0900831,0.3300581,0.0045064
620,1,0.0053712,0.0557136,0.4959814,0.0048696
625,0.9875694,0.0035571,0.0364832,0.6620421,0.005671
630,0.9878027,0.0026321,0.0278317,0.8053931,0.0060849
635,0.9662325,0.0025031,0.0261576,0.8839703,0.0067965
640,0.94831,0.0024856,0.027381,0.9112235,0.0072915
645,0.9197402,0.0025446,0.0275229,0.9418237,0.008561
650,0.9183001,0.0027851,0.0241831,0.9798105,0.0100233
655,0.8943635,0.0026625,0.0185783,0.9958087,0.0112326
660,0.8880265,0.0027425,0.0147736,1,0.0122499
665,0.8781178,0.003288,0.0131941,0.9922692,0.0129368
670,0.8949128,0.0044459,0.0132372,0.9523635,0.0135234
675,0.8892688,0.0088555,0.013176,0.872592,0.0143963
680,0.8822929,0.0161798,0.011635,0.7530428,0.0155814
685,0.8915592,0.0241669,0.0091548,0.6004284,0.0172409
690,0.8780062,0.0270878,0.0060522,0.423465,0.0195995
695,0.897761,0.0190696,0.0043438,0.2806013,0.0240677
700,0.8888378,0.008303,0.0032056,0.1650672,0.0294353
705,0.8713531,0.0056106,0.002855,0.0957075,0.0343437
710,0.8667318,0.0069981,0.0028712,0.0531459,0.0406824
715,0.8687264,0.0101687,0.00405,0.0302447,0.0489796
720,0.8711503,0.0120955,0.0061748,0.0219501,0.059924
725,0.8624942,0.0091411,0.007627,0.0200762,0.0728935
730,0.8810285,0.0036613,0.0076368,0.020729,0.0905184
735,0.881845,0.0026953,0.0056925,0.0214752,0.1194788
740,0.8976004,0.002929,0.0037252,0.0196225,0.1624776
745,0.9071913,0.0026597,0.0026429,0.014756,0.2423822
750,0.9295069,0.0030559,0.0021649,0.0095228,0.360557
755,0.9326408,0.0032472,0.0026131,0.006632,0.4930289
760,0.906958,0.0035759,0.0030395,0.0050248,0.6349644
765,0.9191317,0.0032778,0.002827,0.0045135,0.7983151
770,0.9155077,0.003613,0.0030903,0.0041632,0.9469408
775,0.9038986,0.0036439,0.0033719,0.004527,1
780,0.8936451,0.0030898,0.0039527,0.0050636,0.9971452
785,0.8805316,0.0033737,0.0038226,0.006018,0.9701849
790,0.8717182,0.0031601,0.0043911,0.0066649,0.9349833
795,0.8492843,0.003139,0.0041374,0.006201,0.9008774
800,0.8457008,0.0029854,0.0038086,0.004455,0.8697601
805,0.8553407,0.003293,0.0034867,0.003063,0.8515358
810,0.8581889,0.0033344,0.0032214,0.0027053,0.8382033
815,0.8582515,0.0035658,0.0030436,0.0024349,0.8322304
820,0.8671206,0.0037282,0.0030211,0.002471,0.8256435
825,0.8339193,0.0039858,0.00312,0.0024841,0.8203544
830,0.844024,0.0041173,0.003213,0.0023666,0.8069696
835,0.8499554,0.0042771,0.0033022,0.0026162,0.784373
840,0.8553289,0.004175,0.0036992,0.0024444,0.7494696
845,0.8341391,0.004954,0.0044282,0.0026118,0.704013
850,0.8393621,0.0047827,0.0050802,0.0025035,0.6782487
855,0.8075909,0.0058532,0.0068986,0.0025965,0.6493954
860,0.8105692,0.0057233,0.0077976,0.0025047,0.6131506
865,0.7594286,0.0056168,0.0090757,0.0023465,0.5668747
870,0.7638184,0.0053032,0.0092009,0.0023744,0.5316688
875,0.7341904,0.0052007,0.0090439,0.0022134,0.4916918
880,0.7065537,0.0051649,0.0086195,0.0021206,0.4371382
885,0.694113,0.0051542,0.0079798,0.0019849,0.3515431
890,0.6634336,0.0052812,0.0076004,0.0018715,0.2753125
895,0.654453,0.0056987,0.0074959,0.002009,0.2215915
900,0.6185557,0.0068015,0.0076005,0.0021329,0.1604435
905,0.6002698,0.0082558,0.0078219,0.0024529,0.1165635
910,0.5918959,0.0097551,0.0081012,0.0027978,0.0839996
915,0.5744231,0.0105366,0.0081533,0.0027873,0.0671801
920,0.5309092,0.0109482,0.0077694,0.0025152,0.040855
925,0.4968949,0.0093236,0.0067075,0.0015319,0.020482
930,0.4797736,0.0076207,0.0055921,0.0011564,0.0135345
935,0.4463188,0.0061257,0.0046417,0.0010802,0.010032
940,0.4147031,0.0050368,0.003859,0.0009794,0.0077714
945,0.4047691,0.0042263,0.0030915,0.0009748,0.0063828
950,0.3738311,0.0035359,0.0026681,0.0010349,0.005366
955,0.344805,0.0029168,0.0020248,0.0008131,0.004555
960,0.3117545,0.0021153,0.0016715,0.0008282,0.0034837
965,0.2951745,0.0017263,0.0014821,0.0007164,0.0028074
970,0.2728967,0.0015684,0.001242,0.0006983,0.0024346
975,0.2567351,0.0013848,0.0012525,0.0007062,0.0022279
980,0.2427786,0.0012033,0.0010517,0.0006728,0.0020881
985,0.2216402,0.0012069,0.0009425,0.0005846,0.0018583
990,0.2107166,0.0012223,0.0009661,0.000629,0.0018199
995,0.197115,0.0011194,0.0008666,0.0006789,0.0019029
1000,0.177585,0.0011202,0.0007896,0.0007596,0.0016151
1005,0.1629801,0.0010435,0.0008717,0.0006784,0.0017322
1010,0.1484601,0.0010853,0.0008884,0.0006441,0.001725
1015,0.1357785,0.0012269,0.0009619,0.0008511,0.0020431
1020,0.1239647,0.0012386,0.0011177,0.0009305,0.0021422
1025,0.1085287,0.0014607,0.0012735,0.0007567,0.0020229
1030,0.0954432,0.0018346,0.0010511,0.0008329,0.0019033
1035,0.0831004,0.0014991,0.0010623,0.0009111,0.0017056
1040,0.0715127,0.0015231,0.0008597,0.0008798,0.0018068
1045,0.0634132,0.0014061,0.0009766,0.0006438,0.0016984
1050,0.0548069,0.001219,0.0008612,0.0008323,0.0020005
1055,0.0475781,0,0,0,0
1060,0.0100588,0,0,0,0
1065,0,0,0,0,0
1070,0,0,0,0,0
1075,0,0,0,0,0
1080,0,0,0,0,0
1085,0,0,0,0,0
1090,0,0,0,0,0
1095,0,0,0,0,0
1100,0,0,0,0,0
/* Following filter function created using create_iwave.py */
void IWave::quickbird2(int iwa)
{
/* pan of quickbird2 */
static const float sr1[321] = {
.0000, .0000, .0000, .0000, .0000, .0000, .0000, .0000,
.0000, .0000, .0000, .0000, .0000, .0000, .0000, .0000,
.0000, .0000, .0000, .0007, .0014, .0013, .0012, .0011,
.0011, .0010, .0009, .0010, .0012, .0016, .0020, .0044,
.0069, .0089, .0109, .0138, .0167, .0213, .0259, .0311,
.0362, .0434, .0506, .0588, .0669, .0758, .0846, .0939,
.1032, .1125, .1219, .1334, .1450, .1567, .1684, .1789,
.1895, .2011, .2127, .2163, .2198, .2392, .2586, .2669,
.2753, .2854, .2955, .3055, .3154, .3259, .3364, .3410,
.3456, .3501, .3546, .3579, .3611, .3648, .3685, .3733,
.3781, .3875, .3969, .4055, .4141, .4276, .4410, .4543,
.4676, .4829, .4982, .5167, .5352, .5532, .5712, .5925,
.6137, .6382, .6627, .6839, .7051, .7275, .7499, .7705,
.7910, .8105, .8301, .8465, .8629, .8770, .8911, .9081,
.9251, .9293, .9335, .9426, .9518, .9570, .9621, .9610,
.9598, .9753, .9907, .9941, .9975, .9922, .9869, .9935,
1.0000, .9938, .9876, .9877, .9878, .9770, .9662, .9573,
.9483, .9340, .9197, .9190, .9183, .9063, .8944, .8912,
.8880, .8831, .8781, .8865, .8949, .8921, .8893, .8858,
.8823, .8869, .8916, .8848, .8780, .8879, .8978, .8933,
.8888, .8801, .8714, .8690, .8667, .8677, .8687, .8699,
.8712, .8668, .8625, .8718, .8810, .8814, .8818, .8897,
.8976, .9024, .9072, .9183, .9295, .9311, .9326, .9198,
.9070, .9130, .9191, .9173, .9155, .9097, .9039, .8988,
.8936, .8871, .8805, .8761, .8717, .8605, .8493, .8475,
.8457, .8505, .8553, .8568, .8582, .8582, .8583, .8627,
.8671, .8505, .8339, .8390, .8440, .8470, .8500, .8526,
.8553, .8447, .8341, .8368, .8394, .8235, .8076, .8091,
.8106, .7850, .7594, .7616, .7638, .7490, .7342, .7204,
.7066, .7003, .6941, .6788, .6634, .6589, .6545, .6365,
.6186, .6094, .6003, .5961, .5919, .5832, .5744, .5527,
.5309, .5139, .4969, .4883, .4798, .4630, .4463, .4305,
.4147, .4097, .4048, .3893, .3738, .3593, .3448, .3283,
.3118, .3035, .2952, .2840, .2729, .2648, .2567, .2498,
.2428, .2322, .2216, .2162, .2107, .2039, .1971, .1873,
.1776, .1703, .1630, .1557, .1485, .1421, .1358, .1299,
.1240, .1162, .1085, .1020, .0954, .0893, .0831, .0773,
.0715, .0675, .0634, .0591, .0548, .0512, .0476, .0288,
.0101, .0050, .0000, .0000, .0000, .0000, .0000, .0000,
.0000, .0000, .0000, .0000, .0000, .0000, .0000, .0000,
.0000
};
/* blue of quickbird2 */
static const float sr2[321] = {
.0000, .0000, .0000, .0000, .0000, .0000, .0000, .0000,
.0000, .0005, .0010, .0015, .0021, .0026, .0031, .0037,
.0042, .0044, .0045, .0049, .0053, .0047, .0041, .0037,
.0033, .0030, .0027, .0025, .0023, .0018, .0013, .0015,
.0018, .0020, .0022, .0023, .0024, .0025, .0026, .0025,
.0025, .0027, .0029, .0035, .0040, .0053, .0066, .0098,
.0131, .0208, .0284, .0483, .0681, .1116, .1552, .2164,
.2776, .3544, .4313, .5306, .6299, .6770, .7241, .7577,
.7914, .8176, .8439, .8667, .8896, .9037, .9178, .9235,
.9293, .9447, .9602, .9801, 1.0000, .9958, .9916, .9548,
.9179, .8452, .7724, .6729, .5734, .4787, .3840, .3247,
.2653, .2123, .1592, .1204, .0816, .0754, .0692, .0663,
.0633, .0551, .0469, .0381, .0293, .0235, .0176, .0152,
.0127, .0126, .0125, .0127, .0129, .0128, .0126, .0119,
.0111, .0106, .0101, .0098, .0096, .0094, .0093, .0095,
.0098, .0102, .0106, .0105, .0104, .0091, .0079, .0066,
.0054, .0045, .0036, .0031, .0026, .0026, .0025, .0025,
.0025, .0025, .0025, .0027, .0028, .0027, .0027, .0027,
.0027, .0030, .0033, .0039, .0044, .0067, .0089, .0125,
.0162, .0202, .0242, .0256, .0271, .0231, .0191, .0137,
.0083, .0070, .0056, .0063, .0070, .0086, .0102, .0111,
.0121, .0106, .0091, .0064, .0037, .0032, .0027, .0028,
.0029, .0028, .0027, .0029, .0031, .0032, .0032, .0034,
.0036, .0034, .0033, .0034, .0036, .0036, .0036, .0034,
.0031, .0032, .0034, .0033, .0032, .0031, .0031, .0031,
.0030, .0031, .0033, .0033, .0033, .0035, .0036, .0036,
.0037, .0039, .0040, .0041, .0041, .0042, .0043, .0042,
.0042, .0046, .0050, .0049, .0048, .0053, .0059, .0058,
.0057, .0057, .0056, .0055, .0053, .0053, .0052, .0052,
.0052, .0052, .0052, .0052, .0053, .0055, .0057, .0063,
.0068, .0075, .0083, .0090, .0098, .0101, .0105, .0107,
.0109, .0101, .0093, .0085, .0076, .0069, .0061, .0056,
.0050, .0046, .0042, .0039, .0035, .0032, .0029, .0025,
.0021, .0019, .0017, .0016, .0016, .0015, .0014, .0013,
.0012, .0012, .0012, .0012, .0012, .0012, .0011, .0011,
.0011, .0011, .0010, .0011, .0011, .0012, .0012, .0012,
.0012, .0013, .0015, .0016, .0018, .0017, .0015, .0015,
.0015, .0015, .0014, .0013, .0012, .0006, .0000, .0000,
.0000, .0000, .0000, .0000, .0000, .0000, .0000, .0000,
.0000, .0000, .0000, .0000, .0000, .0000, .0000, .0000,
.0000
};
/* green of quickbird2 */
static const float sr3[321] = {
.0000, .0000, .0000, .0000, .0000, .0000, .0000, .0000,
.0000, .0000, .0000, .0000, .0000, .0000, .0000, .0000,
.0000, .0000, .0000, .0015, .0031, .0036, .0040, .0031,
.0022, .0023, .0023, .0019, .0015, .0015, .0014, .0011,
.0008, .0010, .0011, .0012, .0012, .0012, .0012, .0012,
.0012, .0012, .0011, .0012, .0012, .0013, .0014, .0016,
.0019, .0021, .0023, .0026, .0028, .0033, .0037, .0046,
.0055, .0077, .0098, .0138, .0178, .0206, .0234, .0268,
.0301, .0356, .0412, .0504, .0597, .0789, .0982, .1274,
.1566, .1924, .2282, .2739, .3196, .3686, .4175, .4726,
.5277, .5928, .6579, .7138, .7698, .8005, .8312, .8378,
.8445, .8575, .8705, .8893, .9082, .9271, .9460, .9640,
.9819, .9910, 1.0000, .9967, .9934, .9876, .9818, .9755,
.9691, .9733, .9774, .9822, .9869, .9761, .9652, .9316,
.8979, .8445, .7910, .7230, .6551, .5687, .4822, .3973,
.3124, .2579, .2033, .1706, .1378, .1139, .0901, .0729,
.0557, .0461, .0365, .0322, .0278, .0270, .0262, .0268,
.0274, .0275, .0275, .0259, .0242, .0214, .0186, .0167,
.0148, .0140, .0132, .0132, .0132, .0132, .0132, .0124,
.0116, .0104, .0092, .0076, .0061, .0052, .0043, .0038,
.0032, .0030, .0029, .0029, .0029, .0035, .0040, .0051,
.0062, .0069, .0076, .0076, .0076, .0067, .0057, .0047,
.0037, .0032, .0026, .0024, .0022, .0024, .0026, .0028,
.0030, .0029, .0028, .0030, .0031, .0032, .0034, .0037,
.0040, .0039, .0038, .0041, .0044, .0043, .0041, .0040,
.0038, .0036, .0035, .0034, .0032, .0031, .0030, .0030,
.0030, .0031, .0031, .0032, .0032, .0033, .0033, .0035,
.0037, .0041, .0044, .0048, .0051, .0060, .0069, .0073,
.0078, .0084, .0091, .0091, .0092, .0091, .0090, .0088,
.0086, .0083, .0080, .0078, .0076, .0075, .0075, .0075,
.0076, .0077, .0078, .0080, .0081, .0081, .0082, .0080,
.0078, .0072, .0067, .0061, .0056, .0051, .0046, .0043,
.0039, .0035, .0031, .0029, .0027, .0023, .0020, .0018,
.0017, .0016, .0015, .0014, .0012, .0012, .0013, .0012,
.0011, .0010, .0009, .0010, .0010, .0009, .0009, .0008,
.0008, .0008, .0009, .0009, .0009, .0009, .0010, .0010,
.0011, .0012, .0013, .0012, .0011, .0011, .0011, .0010,
.0009, .0009, .0010, .0009, .0009, .0004, .0000, .0000,
.0000, .0000, .0000, .0000, .0000, .0000, .0000, .0000,
.0000, .0000, .0000, .0000, .0000, .0000, .0000, .0000,
.0000
};
/* red of quickbird2 */
static const float sr4[321] = {
.0000, .0000, .0000, .0000, .0000, .0000, .0000, .0000,
.0000, .0000, .0000, .0000, .0000, .0000, .0000, .0000,
.0000, .0000, .0000, .0019, .0037, .0031, .0024, .0020,
.0016, .0018, .0019, .0016, .0012, .0011, .0011, .0010,
.0009, .0009, .0010, .0009, .0009, .0011, .0013, .0010,
.0007, .0009, .0012, .0012, .0013, .0013, .0013, .0012,
.0012, .0011, .0010, .0013, .0015, .0018, .0021, .0024,
.0028, .0029, .0031, .0030, .0029, .0029, .0028, .0029,
.0030, .0032, .0035, .0036, .0037, .0034, .0032, .0026,
.0021, .0021, .0022, .0023, .0024, .0025, .0027, .0029,
.0031, .0032, .0033, .0034, .0035, .0036, .0037, .0041,
.0046, .0049, .0053, .0056, .0060, .0061, .0063, .0065,
.0066, .0068, .0069, .0073, .0078, .0084, .0090, .0096,
.0102, .0110, .0118, .0128, .0139, .0158, .0177, .0207,
.0238, .0286, .0334, .0399, .0464, .0548, .0631, .0751,
.0872, .1077, .1281, .1666, .2050, .2675, .3301, .4130,
.4960, .5790, .6620, .7337, .8054, .8447, .8840, .8976,
.9112, .9265, .9418, .9608, .9798, .9878, .9958, .9979,
1.0000, .9961, .9923, .9723, .9524, .9125, .8726, .8128,
.7530, .6767, .6004, .5119, .4235, .3520, .2806, .2228,
.1651, .1304, .0957, .0744, .0531, .0417, .0302, .0261,
.0220, .0210, .0201, .0204, .0207, .0211, .0215, .0205,
.0196, .0172, .0148, .0121, .0095, .0081, .0066, .0058,
.0050, .0048, .0045, .0043, .0042, .0043, .0045, .0048,
.0051, .0055, .0060, .0063, .0067, .0064, .0062, .0053,
.0045, .0038, .0031, .0029, .0027, .0026, .0024, .0025,
.0025, .0025, .0025, .0024, .0024, .0025, .0026, .0025,
.0024, .0025, .0026, .0026, .0025, .0026, .0026, .0026,
.0025, .0024, .0023, .0024, .0024, .0023, .0022, .0022,
.0021, .0021, .0020, .0019, .0019, .0019, .0020, .0021,
.0021, .0023, .0025, .0026, .0028, .0028, .0028, .0027,
.0025, .0020, .0015, .0013, .0012, .0011, .0011, .0010,
.0010, .0010, .0010, .0010, .0010, .0009, .0008, .0008,
.0008, .0008, .0007, .0007, .0007, .0007, .0007, .0007,
.0007, .0006, .0006, .0006, .0006, .0007, .0007, .0007,
.0008, .0007, .0007, .0007, .0006, .0007, .0009, .0009,
.0009, .0008, .0008, .0008, .0008, .0009, .0009, .0009,
.0009, .0008, .0006, .0007, .0008, .0004, .0000, .0000,
.0000, .0000, .0000, .0000, .0000, .0000, .0000, .0000,
.0000, .0000, .0000, .0000, .0000, .0000, .0000, .0000,
.0000
};
/* nir of quickbird2 */
static const float sr5[321] = {
.0000, .0000, .0000, .0000, .0000, .0000, .0000, .0000,
.0000, .0007, .0014, .0021, .0029, .0036, .0044, .0051,
.0059, .0061, .0063, .0062, .0060, .0057, .0054, .0045,
.0035, .0028, .0021, .0023, .0024, .0021, .0017, .0020,
.0023, .0021, .0019, .0017, .0016, .0018, .0019, .0019,
.0018, .0017, .0015, .0015, .0015, .0017, .0019, .0018,
.0017, .0021, .0024, .0025, .0025, .0028, .0031, .0036,
.0040, .0044, .0047, .0047, .0047, .0048, .0049, .0045,
.0042, .0043, .0044, .0046, .0048, .0050, .0052, .0054,
.0056, .0054, .0052, .0047, .0042, .0036, .0030, .0028,
.0026, .0026, .0026, .0026, .0025, .0027, .0029, .0027,
.0026, .0030, .0034, .0043, .0053, .0066, .0079, .0089,
.0099, .0101, .0102, .0090, .0079, .0064, .0049, .0049,
.0048, .0056, .0063, .0070, .0077, .0078, .0079, .0072,
.0065, .0054, .0044, .0041, .0038, .0038, .0037, .0037,
.0036, .0036, .0036, .0038, .0039, .0042, .0045, .0047,
.0049, .0053, .0057, .0059, .0061, .0064, .0068, .0070,
.0073, .0079, .0086, .0093, .0100, .0106, .0112, .0117,
.0122, .0126, .0129, .0132, .0135, .0140, .0144, .0150,
.0156, .0164, .0172, .0184, .0196, .0218, .0241, .0268,
.0294, .0319, .0343, .0375, .0407, .0448, .0490, .0545,
.0599, .0664, .0729, .0817, .0905, .1050, .1195, .1410,
.1625, .2024, .2424, .3015, .3606, .4268, .4930, .5640,
.6350, .7166, .7983, .8726, .9469, .9735, 1.0000, .9986,
.9971, .9837, .9702, .9526, .9350, .9179, .9009, .8853,
.8698, .8606, .8515, .8449, .8382, .8352, .8322, .8289,
.8256, .8230, .8204, .8137, .8070, .7957, .7844, .7669,
.7495, .7267, .7040, .6911, .6782, .6638, .6494, .6313,
.6132, .5900, .5669, .5493, .5317, .5117, .4917, .4644,
.4371, .3943, .3515, .3134, .2753, .2485, .2216, .1910,
.1604, .1385, .1166, .1003, .0840, .0756, .0672, .0540,
.0409, .0307, .0205, .0170, .0135, .0118, .0100, .0089,
.0078, .0071, .0064, .0059, .0054, .0050, .0046, .0040,
.0035, .0031, .0028, .0026, .0024, .0023, .0022, .0022,
.0021, .0020, .0019, .0018, .0018, .0019, .0019, .0018,
.0016, .0017, .0017, .0017, .0017, .0019, .0020, .0021,
.0021, .0021, .0020, .0020, .0019, .0018, .0017, .0018,
.0018, .0018, .0017, .0018, .0020, .0010, .0000, .0000,
.0000, .0000, .0000, .0000, .0000, .0000, .0000, .0000,
.0000, .0000, .0000, .0000, .0000, .0000, .0000, .0000,
.0000
};
static const float wli[5] = {0.300, 0.300, 0.300, 0.300, 0.300};
static const float wls[5] = {1.100, 1.100, 1.100, 1.100, 1.100};
ffu.wlinf = (float)wli[iwa-1];
ffu.wlsup = (float)wls[iwa-1];
int i;
for(i = 0; i < 1501; i++) ffu.s[i] = 0;
switch(iwa)
{
case 1: for(i = 0; i < 321; i++) ffu.s[20+i] = sr1[i];
break;
case 2: for(i = 0; i < 321; i++) ffu.s[20+i] = sr2[i];
break;
case 3: for(i = 0; i < 321; i++) ffu.s[20+i] = sr3[i];
break;
case 4: for(i = 0; i < 321; i++) ffu.s[20+i] = sr4[i];
break;
case 5: for(i = 0; i < 321; i++) ffu.s[20+i] = sr5[i];
break;
}
}
_______________________________________________
grass-dev mailing list
[email protected]
http://lists.osgeo.org/mailman/listinfo/grass-dev
