[GRASS-SVN] r61314 - in grass/branches/releasebranch_7_0: . imagery/i.atcorr imagery/i.atcorr/sensors_csv
svn_grass at osgeo.org
svn_grass at osgeo.org
Tue Jul 22 03:57:04 PDT 2014
Author: neteler
Date: 2014-07-22 03:57:03 -0700 (Tue, 22 Jul 2014)
New Revision: 61314
Added:
grass/branches/releasebranch_7_0/imagery/i.atcorr/sensors_csv/landsat_8.csv
Modified:
grass/branches/releasebranch_7_0/
grass/branches/releasebranch_7_0/imagery/i.atcorr/create_iwave.py
grass/branches/releasebranch_7_0/imagery/i.atcorr/geomcond.cpp
grass/branches/releasebranch_7_0/imagery/i.atcorr/geomcond.h
grass/branches/releasebranch_7_0/imagery/i.atcorr/i.atcorr.html
grass/branches/releasebranch_7_0/imagery/i.atcorr/iwave.cpp
grass/branches/releasebranch_7_0/imagery/i.atcorr/iwave.h
Log:
i.atcorr: backport of Landsat8 support (r60610 + r60761, trac #2305)
Property changes on: grass/branches/releasebranch_7_0
___________________________________________________________________
Modified: svn:mergeinfo
- /grass/trunk:60219,60807,61275,61288,61290,61292,61294,61301
+ /grass/trunk:60219,60610,60807,61275,61288,61290,61292,61294,61301
Modified: grass/branches/releasebranch_7_0/imagery/i.atcorr/create_iwave.py
===================================================================
--- grass/branches/releasebranch_7_0/imagery/i.atcorr/create_iwave.py 2014-07-22 09:48:50 UTC (rev 61313)
+++ grass/branches/releasebranch_7_0/imagery/i.atcorr/create_iwave.py 2014-07-22 10:57:03 UTC (rev 61314)
@@ -87,6 +87,7 @@
w = values[:,1] >= 0
response = values[w]
+
# interpolating
f = interpolate.interp1d(response[:,0],response[:,1])
Modified: grass/branches/releasebranch_7_0/imagery/i.atcorr/geomcond.cpp
===================================================================
--- grass/branches/releasebranch_7_0/imagery/i.atcorr/geomcond.cpp 2014-07-22 09:48:50 UTC (rev 61313)
+++ grass/branches/releasebranch_7_0/imagery/i.atcorr/geomcond.cpp 2014-07-22 10:57:03 UTC (rev 61314)
@@ -369,6 +369,9 @@
case 13: /* rapideye * enter month,day,hh.ddd,long.,lat. */
case 14: /* vgt1_spot4 * enter month,day,hh.ddd,long.,lat. */
case 15: /* vgt2_spot5 * enter month,day,hh.ddd,long.,lat. */
+ case 16: /* worldview2 * enter month,day,hh.ddd,long.,lat. */
+ case 17: /* quickbird2 * enter month,day,hh.ddd,long.,lat. */
+ case 18: /* Landsat 8 * enter month,day,hh.ddd,long.,lat. */
{
cin >> month;
cin >> jday;
@@ -414,24 +417,27 @@
/* ---- print geometrical conditions ---- */
void GeomCond::print()
{
- static const string etiq1[16] = {
+ static const string etiq1[] = {
string(" user defined conditions "),
- string(" meteosat observation "),
- string(" goes east observation "),
- string(" goes west observation "),
- string(" avhrr (AM noaa) observation "),
- string(" avhrr (PM noaa) observation "),
- string(" h.r.v. observation "),
- string(" t.m. observation "),
- string(" etm+ observation "),
- string(" liss observation "),
- string(" aster observation "),
- string(" avnir observation "),
- string(" ikonos observation "),
- string(" rapideye observation "),
- string(" vgt1_spot4 observation "),
- string(" vgt2_spot5 observation ")
- };
+ string(" Meteosat observation "),
+ string(" GOES east observation "),
+ string(" GOES west observation "),
+ string(" AVHRR (AM noaa) observation "),
+ string(" AVHRR (PM noaa) observation "),
+ string(" H.R.V. observation "),
+ string(" T.M. observation "),
+ string(" ETM+ observation "),
+ string(" LISS observation "),
+ string(" ASTER observation "),
+ string(" AVNIR observation "),
+ string(" Ikonos observation "),
+ string(" Rapideye observation "),
+ string(" VGT1-SPOT4 observation "),
+ string(" VGT2-SPOT5 observation "),
+ string(" Worldview2 observation "),
+ string(" Quickbird2 observation "),
+ string(" Landsat 8 observation ")
+ };
static const string head(" geometrical conditions identity ");
static const string line(" ------------------------------- ");
Modified: grass/branches/releasebranch_7_0/imagery/i.atcorr/geomcond.h
===================================================================
--- grass/branches/releasebranch_7_0/imagery/i.atcorr/geomcond.h 2014-07-22 09:48:50 UTC (rev 61313)
+++ grass/branches/releasebranch_7_0/imagery/i.atcorr/geomcond.h 2014-07-22 10:57:03 UTC (rev 61314)
@@ -59,15 +59,20 @@
/* c */
/* 15 vgt2_spot5 * enter month,day,hh.ddd,long.,lat.c */
/* c */
+/* 16 worldview2 * enter month,day,hh.ddd,long.,lat.c */
/* c */
+/* 17 quickbird * enter month,day,hh.ddd,long.,lat. c */
+/* c */
+/* 18 Landsat 8 * enter month,day,hh.ddd,long.,lat. c */
+/* c */
/* note: for hrv and tm experiments long. and lat. are the c */
/* coordinates of the scene center. c */
/* lat. must be > 0 for north lat., < 0 for south lat. c */
-/* long. must be > 0 for east long., <0 for west long. c */
+/* long. must be > 0 for east long., < 0 for west long. c */
/* c */
/* solar and viewing positions are computed c */
/* c */
-/* **********************************************************************c */
+/* *********************************************************************c */
struct GeomCond
{
Modified: grass/branches/releasebranch_7_0/imagery/i.atcorr/i.atcorr.html
===================================================================
--- grass/branches/releasebranch_7_0/imagery/i.atcorr/i.atcorr.html 2014-07-22 09:48:50 UTC (rev 61313)
+++ grass/branches/releasebranch_7_0/imagery/i.atcorr/i.atcorr.html 2014-07-22 10:57:03 UTC (rev 61314)
@@ -161,6 +161,25 @@
<td><b>VGT2 (SPOT5)</b></td>
<td>enter month,day,hh.ddd,long.,lat. *</td>
</tr>
+
+<tr>
+<td>16</td>
+<td><b>WorldView 2</b></td>
+<td>enter month,day,hh.ddd,long.,lat. *</td>
+</tr>
+
+<tr>
+<td>17</td>
+<td><b>QuickBird</b></td>
+<td>enter month,day,hh.ddd,long.,lat. *</td>
+</tr>
+
+<tr>
+<td>17</td>
+<td><b>LandSat 8</b></td>
+<td>enter month,day,hh.ddd,long.,lat. *</td>
+</tr>
+
</table>
<blockquote>
@@ -557,6 +576,36 @@
<tr><td>98</td><td>VGT2 (SPOT5) band 2 (0.580-0.780)</td></tr>
<tr><td>99</td><td>VGT2 (SPOT5) band 3 (0.700-1.000)</td></tr>
<tr><td>100</td><td>VGT2 (SPOT5) MIR band (1.450-1.800)</td></tr>
+
+<tr><td>101</td><td>WorldView 2 Panchromatic band (0.447-0.808)</td></tr>
+<tr><td>102</td><td>WorldView 2 Coastal Blue band (0.396-0.458)</td></tr>
+<tr><td>103</td><td>WorldView 2 Blue band (0.442-0.515)</td></tr>
+<tr><td>104</td><td>WorldView 2 Green band (0.506-0.586)</td></tr>
+<tr><td>105</td><td>WorldView 2 Yellow band (0.584-0.632)</td></tr>
+<tr><td>106</td><td>WorldView 2 Red band (0.624-0.694)</td></tr>
+<tr><td>107</td><td>WorldView 2 Red Edge band (0.699-0.749)</td></tr>
+<tr><td>108</td><td>WorldView 2 NIR1 band (0.765-0.901)</td></tr>
+<tr><td>109</td><td>WorldView 2 NIR2 band (0.856-0.1043)</td></tr>
+
+<tr><td>110</td><td><b>QuickBird</b> Panchromatic band (0.405-1.053)</td></tr>
+<tr><td>111</td><td>QuickBird Blue band (0.430-0.545)</td></tr>
+<tr><td>112</td><td>QuickBird Green band (0.466-0.620)</td></tr>
+<tr><td>113</td><td>QuickBird Red band (0.590-0.710)</td></tr>
+<tr><td>114</td><td>QuickBird NIR1 band (0.715-0.918)</td></tr>
+
+<tr><td>115</td><td><b>Landsat 8 </b> Coastal Aerosol Band (0.427nm - 0.459nm)</td></tr>
+<tr><td>116</td><td>Landsat 8 Blue Band (436nm - 527nm)</td></tr>
+<tr><td>117</td><td>Landsat 8 Green Band (512nm-610nm)</td></tr>
+<tr><td>118</td><td>Landsat 8 Red Band (625nm-691nm)</td></tr>
+<tr><td>119</td><td>Landsat 8 Panchromatic Band (488nm-692nm)</td></tr>
+<tr><td>120</td><td>Landsat 8 NIR Band (829nm-900nm)</td></tr>
+<tr><td>121</td><td>Landsat 8 Cirrus Band (1340nm-1409nm)</td></tr>
+<tr><td>122</td><td>Landsat 8 SWIR1 Band (1515nm - 1697nm)</td></tr>
+<tr><td>123</td><td>Landsat 8 SWIR2 Band (2037nm - 2355nm)</td></tr>
+
+
+
+
</table>
<h2>EXAMPLES</h2>
@@ -696,6 +745,7 @@
<a href="http://www.rsgis.ait.ac.th/~honda/textbooks/advrs/6smanv2.0_P2.pdf">PDF2</a>,
and <a href="http://www.rsgis.ait.ac.th/~honda/textbooks/advrs/6smanv2.0_P3.pdf">PDF3</a>
<li>RapidEye sensors have been provided by <a href="http://www.rapideye.com/">RapidEye AG, Germany</a>
+<li> Julia A. Barsi, Brian L. Markham and Jeffrey A. Pedelty "The operational land imager: spectral response and spectral uniformity", Proc. SPIE 8153, 81530G (2011); doi:10.1117/12.895438
</ul>
@@ -722,5 +772,10 @@
<p><em>VGT1 and VGT2 sensors addition from <a href="http://6s.ltdri.org/">6SV-1.1 sources</a>, addition 07/2011:</em>
<br>Alfredo Alessandrini, Anne Ghisla
-<p><i>Last changed: $Date$</i>
+<p><em>Added Landsat 8 from <a href="http://landsat.gsfc.nasa.gov/?p=5779">NASA sources</a>, addition 05/2014:</em>
+<br>Nikolaos Ves
+
+<p>
+<i>Last changed: $Date$</i>
+
Modified: grass/branches/releasebranch_7_0/imagery/i.atcorr/iwave.cpp
===================================================================
--- grass/branches/releasebranch_7_0/imagery/i.atcorr/iwave.cpp 2014-07-22 09:48:50 UTC (rev 61313)
+++ grass/branches/releasebranch_7_0/imagery/i.atcorr/iwave.cpp 2014-07-22 10:57:03 UTC (rev 61314)
@@ -2457,6 +2457,644 @@
}
}
+void IWave::worldview2(int iwa)
+{
+
+ /* WV2 Pan of wv2_spectral_response */
+ static const float sr1[301] = {
+ .0055, .0062, .0068, .0045, .0022, .0024, .0027, .0034,
+ .0041, .0038, .0030, .0039, .0049, .0042, .0034, .0037,
+ .0040, .0050, .0060, .0050, .0040, .0036, .0033, .0037,
+ .0042, .0052, .0063, .0048, .0034, .0061, .0088, .0105,
+ .0121, .0499, .0881, .2723, .4599, 1.6988, 2.9515, 55.9008,
+ 109.7660, 271.3440, 436.1630, 461.8220, 487.5330, 494.8560, 504.9920, 515.2150,
+ 525.5960, 529.5890, 535.5830, 529.7560, 522.2330, 527.2070, 533.6600, 528.6020,
+ 523.0890, 541.1120, 557.2800, 558.1700, 558.2320, 580.8990, 605.5330, 616.2080,
+ 626.7920, 612.1180, 595.4140, 612.0080, 630.7120, 649.1530, 671.5410, 663.7040,
+ 654.7330, 672.5320, 693.0980, 691.3350, 688.2370, 678.2900, 664.9730, 683.3230,
+ 701.5600, 705.9920, 713.5760, 709.8480, 709.7860, 738.4220, 765.7520, 762.7560,
+ 757.8270, 759.1620, 759.8180, 779.9830, 794.0850, 797.9350, 800.7150, 811.3510,
+ 819.2900, 820.8140, 823.0300, 829.7940, 838.8690, 848.9130, 860.8020, 858.0370,
+ 858.4620, 865.1970, 870.2850, 877.9320, 889.2780, 896.4750, 909.2530, 913.0800,
+ 917.8850, 914.2210, 910.7620, 916.9790, 924.8900, 936.8300, 950.7910, 943.9980,
+ 936.9110, 932.3110, 929.4170, 938.2740, 951.2320, 960.0120, 962.7810, 957.5010,
+ 953.8650, 954.7520, 957.3760, 968.1550, 977.7350, 961.3320, 949.2170, 946.9550,
+ 953.3340, 946.6560, 952.1930, 962.4700, 969.5950, 983.9680, 1000.0000, 994.3200,
+ 994.0030, 980.9740, 969.5630, 968.2530, 968.3970, 973.6740, 978.8550, 966.8470,
+ 963.7410, 952.2520, 946.6540, 945.9200, 952.9380, 954.7590, 962.5830, 962.1420,
+ 956.6370, 950.3810, 940.9620, 936.2240, 924.5960, 923.3790, 921.6840, 933.5680,
+ 940.8310, 957.9770, 963.8320, 958.1430, 941.2970, 924.2600, 903.0000, 895.8780,
+ 887.1520, 910.5230, 929.4020, 744.0190, 553.8670, 320.2460, 85.8164, 51.4292,
+ 17.2983, 12.0102, 6.5998, 4.3680, 2.1313, 1.5356, .9222, .7218,
+ .5287, .4897, .4497, .4318, .4090, .4071, .3964, .4203,
+ .4465, .4033, .3613, .3303, .2980, .3011, .3085, .2862,
+ .2615, .2536, .2423, .2746, .3070, .2950, .2837, .2615,
+ .2375, .2483, .2590, .2701, .2825, .2617, .2408, .2486,
+ .2580, .2735, .2876, .2828, .2762, .2753, .2731, .2568,
+ .2422, .2392, .2343, .2234, .2135, .2157, .2186, .2048,
+ .1911, .1799, .1690, .1771, .1832, .1676, .1513, .1388,
+ .1258, .1297, .1337, .1225, .1114, .1029, .0942, .1047,
+ .1152, .0995, .0839, .0845, .0853, .0776, .0701, .0707,
+ .0711, .0657, .0600, .0621, .0642, .0557, .0473, .0433,
+ .0392, .0361, .0329, .0341, .0352, .0314, .0276, .0255,
+ .0234, .0211, .0188, .0180, .0171, .0151, .0131, .0121,
+ .0111, .0100, .0088, .0088, .0087, .0080, .0072, .0071,
+ .0069, .0067, .0064, .0064, .0063
+ };
+
+ /* WV2 Coastal Blue of wv2_spectral_response */
+ static const float sr2[301] = {
+ .0019, .0033, .0033, .0028, .0017, .0033, .0036, .0044,
+ .0029, .0062, .0054, .0050, .0074, .0135, .0454, .1731,
+ .6347, 3.4150, 27.4284, 132.9030, 355.9880, 592.5800, 700.7190, 739.6250,
+ 765.3070, 787.8420, 817.0880, 838.6870, 861.7800, 883.5130, 905.0550, 917.0720,
+ 927.2080, 947.6580, 965.0880, 979.3040, 985.9210, 989.0940, 997.2360, 974.6130,
+ 853.9400, 588.0320, 281.5770, 80.5215, 8.9461, 1.3331, .3609, .1454,
+ .0712, .0380, .0115, .0096, .0051, .0112, .0110, .0108,
+ .0091, .0101, .0117, .0064, .0013, .0014, .0009, .0006,
+ .0007, .0009, .0009, .0011, .0010, .0007, .0009, .0012,
+ .0036, .0013, .0004, .0005, .0003, .0005, .0004, .0126,
+ .1417, .0024, .0003, .0003, .0005, .0005, .0003, .0003,
+ .0001, .0000, .0002, .0016, .0017, .0004, .0001, .0001,
+ .0005, .0029, .0045, .0028, .0010, .0017, .0013, .0014,
+ .0023, .0031, .0061, .0116, .0101, .0056, .0039, .0032,
+ .0027, .0031, .0025, .0026, .0027, .0027, .0025, .0021,
+ .0022, .0018, .0017, .0018, .0018, .0021, .0024, .0030,
+ .0044, .0093, .0172, .0163, .0078, .0042, .0044, .0042,
+ .0036, .0042, .0066, .0185, .0264, .0124, .0081, .0076,
+ .0085, .0096, .0122, .0292, .0976, .0936, .0333, .0238,
+ .0298, .0442, .0669, .1526, .4227, .2563, .0653, .0414,
+ .0578, .2646, .7112, .4411, .0894, .0305, .0131, .0081,
+ .0060, .0059, .0070, .0112, .0193, .0331, .0553, .1062,
+ .2462, .2458, .1162, .0691, .1216, .1248, .1605, .1959,
+ .1637, .1168, .0928, .0797, .0730, .0625, .0716, .0692,
+ .0699, .0721, .0593, .0474, .0512, .0388, .0287, .0282,
+ .0423, .0284, .0276, .0314, .0321, .0255, .0155, .0168,
+ .0152, .0173, .0183, .0163, .0163, .0211, .0168, .0132,
+ .0099, .0096, .0092, .0160, .0189, .0195, .0168, .0153,
+ .0145, .0112, .0095, .0074, .0062, .0047, .0040, .0040,
+ .0030, .0033, .0025, .0030, .0030, .0034, .0036, .0042,
+ .0041, .0050, .0054, .0060, .0054, .0064, .0054, .0049,
+ .0044, .0039, .0037, .0041, .0039, .0040, .0046, .0053,
+ .0049, .0046, .0042, .0035, .0032, .0026, .0024, .0023,
+ .0024, .0023, .0026, .0030, .0031, .0037, .0037, .0037,
+ .0038, .0034, .0031, .0027, .0025, .0021, .0019, .0016,
+ .0015, .0015, .0015, .0017, .0019, .0025, .0029, .0034,
+ .0041, .0043, .0046, .0048, .0047, .0050, .0050, .0056,
+ .0058, .0063, .0068, .0078, .0081
+ };
+
+ /* WV2 Blue of wv2_spectral_response */
+ static const float sr3[301] = {
+ .0008, .0004, .0003, .0002, .0003, .0001, .0002, .0003,
+ .0004, .0005, .0004, .0005, .0004, .0003, .0005, .0006,
+ .0007, .0005, .0006, .0007, .0028, .0024, .0017, .0010,
+ .0026, .0062, .0135, .0275, .0412, .0515, .0701, .1333,
+ .3374, .9579, 2.2771, 5.4329, 19.7821, 81.8229, 245.2000, 503.3090,
+ 700.9700, 757.6050, 778.2550, 788.6260, 800.0000, 801.7370, 809.8590, 821.0700,
+ 826.0450, 830.9250, 842.2620, 862.7060, 878.4820, 876.9430, 878.2670, 887.6130,
+ 907.5760, 930.5590, 950.1920, 976.8010, 1000.0000, 978.9400, 855.7790, 593.8730,
+ 297.6290, 125.6610, 50.3684, 18.1990, 5.7341, 2.9735, 1.6560, .9803,
+ .6049, .3915, .2592, .1681, .1089, .0763, .0514, .0520,
+ .0667, .0623, .0607, .0594, .0626, .0549, .0587, .0488,
+ .0494, .0415, .0310, .0297, .0174, .0175, .0143, .0109,
+ .0128, .0092, .0086, .0092, .0081, .0085, .0090, .0076,
+ .0081, .0091, .0074, .0072, .0080, .0061, .0055, .0061,
+ .0042, .0035, .0041, .0038, .0028, .0038, .0038, .0032,
+ .0042, .0054, .0057, .0057, .0100, .0162, .0175, .0215,
+ .0359, .0457, .0343, .0239, .0243, .0255, .0195, .0134,
+ .0126, .0141, .0120, .0086, .0076, .0088, .0101, .0090,
+ .0070, .0081, .0109, .0134, .0124, .0120, .0149, .0224,
+ .0310, .0326, .0310, .0340, .0454, .0611, .0653, .0523,
+ .0415, .0420, .0510, .0638, .0645, .0492, .0354, .0300,
+ .0320, .0400, .0444, .0393, .0299, .0248, .0245, .0294,
+ .0405, .0507, .0515, .0447, .0595, .0543, .0609, .0773,
+ .0848, .0927, .0950, .0805, .0570, .0509, .0510, .0429,
+ .0566, .0641, .0635, .0645, .0686, .0662, .0780, .0785,
+ .0748, .0557, .0449, .0383, .0324, .0264, .0202, .0219,
+ .0214, .0236, .0270, .0351, .0351, .0308, .0297, .0280,
+ .0215, .0221, .0192, .0175, .0185, .0174, .0181, .0187,
+ .0178, .0187, .0184, .0172, .0158, .0142, .0148, .0135,
+ .0130, .0131, .0118, .0109, .0106, .0099, .0088, .0079,
+ .0070, .0061, .0064, .0069, .0074, .0090, .0105, .0120,
+ .0136, .0146, .0152, .0146, .0128, .0114, .0094, .0078,
+ .0063, .0056, .0053, .0053, .0052, .0054, .0053, .0048,
+ .0038, .0032, .0027, .0022, .0020, .0021, .0022, .0024,
+ .0029, .0029, .0028, .0025, .0021, .0019, .0015, .0014,
+ .0014, .0018, .0022, .0027, .0030, .0032, .0030, .0027,
+ .0025, .0021, .0018, .0021, .0027, .0038, .0047, .0056,
+ .0063, .0059, .0051, .0043, .0033
+ };
+
+ /* WV2 Green of wv2_spectral_response */
+ static const float sr4[301] = {
+ .0004, .0003, .0001, .0002, .0002, .0001, .0002, .0002,
+ .0003, .0003, .0002, .0003, .0004, .0002, .0004, .0005,
+ .0004, .0004, .0004, .0005, .0004, .0005, .0005, .0005,
+ .0005, .0004, .0004, .0005, .0003, .0005, .0004, .0004,
+ .0003, .0004, .0005, .0003, .0002, .0003, .0005, .0004,
+ .0006, .0006, .0005, .0005, .0023, .0044, .0085, .0128,
+ .0141, .0159, .0198, .0212, .0226, .0266, .0356, .0435,
+ .0664, .1438, .3963, 1.1357, 3.0919, 11.1042, 40.5198, 133.7880,
+ 351.5100, 606.0390, 723.0070, 760.0260, 779.3890, 788.8110, 806.8500, 828.0350,
+ 840.7560, 842.0390, 844.3050, 846.8570, 850.3350, 859.7230, 863.5160, 869.1430,
+ 882.5710, 900.7200, 924.4690, 942.1680, 962.2850, 976.6010, 985.7260, 994.7640,
+ 1000.0000, 993.0530, 979.1210, 907.7450, 664.0400, 301.4360, 92.9241, 25.2385,
+ 5.5627, 2.4767, 1.1294, .5085, .2229, .1072, .0539, .0257,
+ .0069, .0035, .0017, .0003, .0018, .0017, .0062, .0086,
+ .0152, .0161, .0257, .0325, .0292, .0270, .0241, .0141,
+ .0002, .0001, .0002, .0002, .0003, .0005, .0007, .0008,
+ .0010, .0008, .0006, .0010, .0011, .0006, .0009, .0011,
+ .0012, .0011, .0009, .0008, .0011, .0011, .0010, .0009,
+ .0011, .0009, .0010, .0011, .0011, .0011, .0010, .0013,
+ .0009, .0011, .0008, .0008, .0010, .0012, .0008, .0010,
+ .0011, .0013, .0011, .0010, .0011, .0012, .0014, .0012,
+ .0017, .0035, .0040, .0027, .0014, .0010, .0008, .0007,
+ .0003, .0007, .0008, .0011, .0010, .0026, .0009, .0002,
+ .0000, .0000, .0000, .0000, .0000, .0000, .0000, .0006,
+ .0051, .0042, .0018, .0005, .0003, .0004, .0008, .0003,
+ .0001, .0023, .0086, .0047, .0017, .0002, .0001, .0017,
+ .0013, .0017, .0010, .0002, .0117, .0095, .0112, .0091,
+ .0090, .0085, .0077, .0067, .0077, .0070, .0060, .0055,
+ .0061, .0069, .0062, .0066, .0066, .0058, .0060, .0054,
+ .0054, .0049, .0040, .0036, .0035, .0035, .0029, .0026,
+ .0023, .0021, .0029, .0024, .0022, .0022, .0020, .0018,
+ .0018, .0017, .0015, .0015, .0012, .0014, .0014, .0013,
+ .0012, .0011, .0013, .0013, .0011, .0013, .0014, .0016,
+ .0016, .0020, .0019, .0018, .0019, .0019, .0017, .0016,
+ .0017, .0014, .0013, .0013, .0012, .0013, .0013, .0015,
+ .0017, .0019, .0023, .0030, .0038, .0052, .0072, .0097,
+ .0132, .0170, .0211, .0255, .0300, .0351, .0407, .0474,
+ .0551, .0653, .0811, .1055, .1401
+ };
+
+ /* WV2 Yellow of wv2_spectral_response */
+ static const float sr5[301] = {
+ .0003, .0001, .0001, .0001, .0000, .0000, .0001, .0001,
+ .0001, .0001, .0001, .0001, .0001, .0001, .0001, .0003,
+ .0006, .0211, .0273, .0005, .0000, .0001, .0000, .0001,
+ .0002, .0002, .0001, .0002, .0002, .0003, .0003, .0002,
+ .0003, .0004, .0004, .0004, .0010, .0050, .0032, .0012,
+ .0005, .0003, .0004, .0003, .0007, .0006, .0005, .0005,
+ .0004, .0007, .0007, .0007, .0007, .0007, .0007, .0006,
+ .0007, .0006, .0003, .0004, .0002, .0003, .0002, .0002,
+ .0002, .0002, .0001, .0003, .0001, .0002, .0004, .0003,
+ .0005, .0002, .0002, .0002, .0001, .0004, .0004, .0038,
+ .0070, .0062, .0044, .0049, .0038, .0111, .0267, .0575,
+ .1271, .3445, .8950, 2.2247, 5.8789, 26.5572, 134.6520, 389.2200,
+ 663.5500, 848.1980, 899.8060, 920.0230, 932.6340, 940.1810, 947.0260, 952.4410,
+ 966.4120, 981.4690, 989.7230, 992.7060, 1000.0000, 969.3680, 759.4330, 436.6510,
+ 134.3220, 30.8525, 7.3852, 2.2367, .6906, .2851, .1576, .0963,
+ .0443, .0231, .0129, .0075, .0050, .0033, .0022, .0014,
+ .0008, .0010, .0008, .0007, .0009, .0005, .0008, .0009,
+ .0006, .0009, .0005, .0007, .0010, .0005, .0009, .0008,
+ .0007, .0011, .0009, .0008, .0012, .0011, .0015, .0013,
+ .0012, .0011, .0013, .0013, .0010, .0009, .0011, .0011,
+ .0010, .0008, .0009, .0008, .0008, .0009, .0008, .0008,
+ .0006, .0005, .0005, .0006, .0007, .0002, .0004, .0006,
+ .0003, .0004, .0006, .0006, .0234, .0159, .0189, .0213,
+ .0113, .0140, .0189, .0200, .0194, .0099, .0183, .0115,
+ .0093, .0135, .0106, .0113, .0119, .0076, .0120, .0127,
+ .0171, .0114, .0067, .0086, .0089, .0068, .0044, .0064,
+ .0049, .0058, .0060, .0049, .0019, .0039, .0036, .0029,
+ .0033, .0030, .0031, .0047, .0039, .0043, .0039, .0038,
+ .0033, .0028, .0025, .0021, .0020, .0017, .0018, .0018,
+ .0010, .0014, .0013, .0013, .0007, .0008, .0009, .0009,
+ .0008, .0008, .0009, .0009, .0006, .0010, .0009, .0007,
+ .0009, .0007, .0007, .0008, .0009, .0008, .0010, .0014,
+ .0015, .0016, .0018, .0018, .0015, .0014, .0012, .0010,
+ .0005, .0004, .0005, .0004, .0004, .0005, .0004, .0003,
+ .0003, .0004, .0004, .0004, .0003, .0004, .0003, .0003,
+ .0002, .0002, .0002, .0001, .0001, .0001, .0001, .0001,
+ .0001, .0001, .0001, .0001, .0001, .0002, .0002, .0003,
+ .0003, .0004, .0005, .0005, .0006
+ };
+
+ /* WV2 Red of wv2_spectral_response */
+ static const float sr6[301] = {
+ .0003, .0001, .0001, .0000, .0001, .0000, .0001, .0001,
+ .0002, .0001, .0001, .0002, .0002, .0001, .0002, .0003,
+ .0004, .0002, .0003, .0004, .0003, .0003, .0003, .0004,
+ .0004, .0004, .0006, .0006, .0003, .0003, .0003, .0003,
+ .0003, .0003, .0002, .0002, .0003, .0002, .0004, .0002,
+ .0003, .0003, .0003, .0002, .0019, .0034, .0008, .0006,
+ .0012, .0024, .0010, .0009, .0006, .0013, .0013, .0014,
+ .0029, .0019, .0029, .0022, .0015, .0023, .0008, .0002,
+ .0003, .0002, .0002, .0002, .0002, .0002, .0001, .0000,
+ .0002, .0000, .0000, .0000, .0001, .0001, .0001, .0001,
+ .0001, .0000, .0001, .0002, .0001, .0001, .0001, .0001,
+ .0002, .0001, .0001, .0003, .0002, .0004, .0005, .0006,
+ .0008, .0009, .0013, .0091, .0126, .0268, .0410, .0704,
+ .1229, .2687, .6281, 1.6492, 4.8674, 22.7608, 95.5659, 294.3530,
+ 596.0710, 822.2890, 872.7760, 885.3350, 927.9280, 955.6950, 965.0420, 967.0950,
+ 965.5880, 957.7760, 955.5000, 952.8310, 955.8010, 965.2530, 970.3580, 976.4930,
+ 987.5260, 996.5750, 999.0850, 992.1070, 983.6940, 978.3110, 885.2800, 632.5950,
+ 291.8230, 83.2547, 20.4412, 5.6355, 1.9137, .8432, .4036, .2037,
+ .1018, .0569, .0272, .0162, .0142, .0168, .0133, .0144,
+ .0132, .0125, .0106, .0061, .0031, .0034, .0015, .0014,
+ .0012, .0004, .0006, .0012, .0030, .0050, .0095, .0073,
+ .0067, .0089, .0069, .0049, .0055, .0022, .0015, .0005,
+ .0022, .0007, .0006, .0009, .0138, .0140, .0134, .0163,
+ .0106, .0124, .0108, .0077, .0089, .0090, .0079, .0043,
+ .0108, .0045, .0070, .0065, .0062, .0064, .0084, .0069,
+ .0027, .0065, .0063, .0065, .0051, .0071, .0092, .0108,
+ .0067, .0079, .0041, .0043, .0086, .0062, .0072, .0061,
+ .0039, .0055, .0048, .0050, .0057, .0051, .0042, .0044,
+ .0040, .0035, .0035, .0040, .0027, .0031, .0032, .0024,
+ .0031, .0027, .0021, .0017, .0017, .0016, .0013, .0011,
+ .0012, .0010, .0013, .0010, .0011, .0009, .0008, .0007,
+ .0007, .0007, .0008, .0007, .0007, .0006, .0006, .0006,
+ .0005, .0005, .0005, .0006, .0005, .0006, .0005, .0005,
+ .0002, .0006, .0005, .0004, .0004, .0004, .0003, .0003,
+ .0004, .0005, .0005, .0006, .0006, .0008, .0007, .0006,
+ .0005, .0003, .0003, .0002, .0001, .0002, .0001, .0001,
+ .0002, .0002, .0002, .0002, .0002, .0002, .0002, .0002,
+ .0001, .0001, .0001, .0001, .0001
+ };
+
+ /* WV2 RedEdge of wv2_spectral_response */
+ static const float sr7[301] = {
+ .0001, .0002, .0003, .0002, .0002, .0001, .0001, .0002,
+ .0002, .0002, .0002, .0002, .0003, .0003, .0003, .0004,
+ .0008, .0023, .0013, .0013, .0020, .0011, .0016, .0010,
+ .0010, .0011, .0009, .0020, .0011, .0009, .0007, .0010,
+ .0097, .0509, .0025, .0003, .0008, .0239, .0026, .0002,
+ .0002, .0003, .0001, .0002, .0002, .0002, .0002, .0001,
+ .0001, .0005, .0008, .0017, .0009, .0058, .0138, .0079,
+ .0039, .0057, .0153, .0099, .0012, .0015, .0042, .0014,
+ .0002, .0002, .0001, .0001, .0002, .0001, .0001, .0001,
+ .0002, .0001, .0002, .0002, .0001, .0001, .0001, .0000,
+ .0000, .0001, .0002, .0001, .0004, .0020, .0021, .0008,
+ .0009, .0020, .0031, .0031, .0045, .0067, .0061, .0050,
+ .0071, .0113, .0138, .0161, .0240, .0344, .0375, .0398,
+ .0538, .0794, .0884, .0759, .0904, .1744, .1848, .0517,
+ .0088, .0045, .0027, .0015, .0011, .0007, .0014, .0028,
+ .0050, .0064, .0081, .0106, .0116, .0115, .0118, .0105,
+ .0113, .0084, .0082, .0089, .0097, .0153, .0328, .0787,
+ .2073, .7591, 3.1275, 18.8368, 99.6101, 346.8350, 622.9080, 853.7850,
+ 952.7920, 978.4220, 987.6340, 989.0260, 993.9010, 999.8580, 999.2800, 989.4170,
+ 986.8130, 980.9620, 977.1160, 944.1910, 846.3990, 615.8370, 359.7630, 138.8130,
+ 30.5442, 7.0917, 2.1585, .7300, .2699, .1188, .0553, .0240,
+ .0096, .0069, .0053, .0033, .0043, .0040, .0051, .0077,
+ .0083, .0094, .0064, .0050, .0001, .0003, .0003, .0003,
+ .0002, .0003, .0002, .0002, .0002, .0001, .0001, .0003,
+ .0002, .0002, .0002, .0001, .0003, .0003, .0002, .0004,
+ .0004, .0002, .0008, .0008, .0019, .0039, .0088, .0112,
+ .0114, .0077, .0047, .0043, .0104, .0084, .0076, .0074,
+ .0089, .0118, .0153, .0190, .0174, .0171, .0183, .0187,
+ .0180, .0190, .0157, .0145, .0127, .0102, .0081, .0062,
+ .0045, .0039, .0032, .0025, .0020, .0016, .0015, .0011,
+ .0010, .0010, .0008, .0007, .0008, .0006, .0007, .0005,
+ .0006, .0006, .0007, .0006, .0007, .0004, .0004, .0006,
+ .0004, .0005, .0005, .0005, .0005, .0004, .0004, .0004,
+ .0004, .0004, .0004, .0004, .0003, .0003, .0002, .0002,
+ .0002, .0002, .0002, .0002, .0002, .0001, .0001, .0001,
+ .0001, .0002, .0002, .0001, .0002, .0002, .0002, .0002,
+ .0002, .0001, .0001, .0001, .0001, .0001, .0001, .0000,
+ .0001, .0001, .0001, .0003, .0005
+ };
+
+ /* WV2 NIR1 of wv2_spectral_response */
+ static const float sr8[301] = {
+ .0003, .0000, .0000, .0001, .0001, .0000, .0000, .0001,
+ .0001, .0000, .0000, .0000, .0000, .0000, .0000, .0000,
+ .0001, .0000, .0000, .0001, .0002, .0000, .0001, .0001,
+ .0001, .0002, .0001, .0001, .0001, .0001, .0001, .0001,
+ .0001, .0002, .0001, .0001, .0016, .0061, .0004, .0001,
+ .0004, .0003, .0003, .0003, .0003, .0004, .0005, .0006,
+ .0004, .0004, .0007, .0014, .0006, .0009, .0016, .0006,
+ .0005, .0006, .0003, .0002, .0003, .0001, .0002, .0004,
+ .0002, .0003, .0003, .0004, .0003, .0004, .0002, .0001,
+ .0001, .0000, .0000, .0000, .0000, .0000, .0001, .0001,
+ .0000, .0001, .0001, .0002, .0001, .0000, .0001, .0001,
+ .0000, .0000, .0001, .0000, .0000, .0000, .0000, .0001,
+ .0001, .0000, .0000, .0000, .0002, .0000, .0000, .0002,
+ .0000, .0001, .0000, .0001, .0000, .0000, .0000, .0001,
+ .0000, .0000, .0000, .0001, .0000, .0001, .0001, .0004,
+ .0005, .0005, .0006, .0005, .0007, .0007, .0010, .0008,
+ .0011, .0008, .0009, .0010, .0010, .0009, .0009, .0012,
+ .0013, .0010, .0009, .0021, .0017, .0035, .0048, .0064,
+ .0052, .0043, .0039, .0044, .0050, .0060, .0066, .0111,
+ .0162, .0284, .0440, .0705, .1061, .1653, .2562, .4123,
+ .7130, 1.3994, 2.8857, 5.9533, 12.6868, 27.6938, 58.8599, 123.1170,
+ 260.2040, 506.7740, 759.9670, 940.5800, 995.8050, 997.6220, 988.9070, 983.0910,
+ 972.9340, 958.3120, 941.1390, 931.9070, 919.3820, 909.0800, 897.6000, 882.1200,
+ 871.1910, 863.0540, 844.2840, 836.4960, 838.0360, 841.1890, 835.6570, 833.8190,
+ 842.5480, 830.7770, 818.6040, 816.8370, 806.7110, 810.7600, 794.1610, 779.8840,
+ 767.3130, 749.6020, 739.2710, 738.7140, 734.1450, 711.8330, 697.1410, 680.4620,
+ 658.6350, 649.7810, 632.4710, 619.9080, 609.1200, 598.9920, 590.9140, 566.1070,
+ 507.1120, 402.7990, 275.1740, 157.6670, 73.8047, 32.7870, 13.8147, 5.3908,
+ 2.6246, 1.8086, 1.3063, 1.0328, .9172, .8912, .7983, .7457,
+ .7272, .7833, .7011, .6553, .6028, .5579, .2812, .2742,
+ .0044, .0032, .0029, .0022, .0020, .0015, .0011, .0010,
+ .0011, .0009, .0009, .0007, .0007, .0007, .0007, .0007,
+ .0004, .0006, .0007, .0006, .0003, .0005, .0005, .0005,
+ .0003, .0004, .0004, .0003, .0002, .0005, .0003, .0003,
+ .0003, .0004, .0003, .0003, .0002, .0002, .0002, .0002,
+ .0001, .0002, .0002, .0001, .0001, .0001, .0001, .0001,
+ .0001, .0001, .0001, .0001, .0001, .0001, .0001, .0001,
+ .0001, .0001, .0001, .0000, .0001
+ };
+
+ /* WV2 NIR2 of wv2_spectral_response */
+ static const float sr9[301] = {
+ .0003, .0070, .0199, .0874, .0748, .0052, .0005, .0006,
+ .0003, .0001, .0002, .0003, .0001, .0002, .0002, .0002,
+ .0002, .0002, .0002, .0002, .0002, .0001, .0002, .0002,
+ .0004, .0002, .0002, .0002, .0001, .0002, .0002, .0002,
+ .0001, .0002, .0001, .0008, .0012, .0001, .0002, .0019,
+ .0012, .0008, .0006, .0002, .0000, .0001, .0004, .0002,
+ .0000, .0009, .0037, .0012, .0002, .0001, .0007, .0009,
+ .0000, .0003, .0001, .0000, .0007, .0008, .0000, .0001,
+ .0002, .0004, .0008, .0009, .0004, .0002, .0002, .0027,
+ .0088, .0007, .0004, .0040, .0047, .0001, .0003, .0021,
+ .0010, .0001, .0002, .0001, .0001, .0000, .0000, .0000,
+ .0000, .0000, .0000, .0001, .0013, .0001, .0003, .0002,
+ .0002, .0005, .0006, .0002, .0006, .0005, .0005, .0006,
+ .0012, .0005, .0003, .0002, .0005, .0006, .0006, .0011,
+ .0011, .0003, .0004, .0005, .0001, .0006, .0007, .0008,
+ .0012, .0019, .0024, .0010, .0005, .0015, .0005, .0000,
+ .0003, .0002, .0003, .0005, .0008, .0002, .0002, .0005,
+ .0008, .0002, .0004, .0003, .0001, .0002, .0001, .0004,
+ .0005, .0007, .0007, .0009, .0009, .0006, .0010, .0009,
+ .0007, .0007, .0009, .0011, .0006, .0009, .0010, .0010,
+ .0012, .0009, .0008, .0005, .0008, .0007, .0009, .0006,
+ .0005, .0011, .0010, .0005, .0006, .0009, .0009, .0008,
+ .0012, .1141, .4534, .5047, .7796, 1.0461, .9763, 1.1133,
+ 1.0687, 1.0909, 1.0955, 1.1311, 1.0945, 1.0891, 1.0416, 1.1403,
+ 1.1270, 1.1706, 1.1509, 1.2061, 1.3547, 1.3002, 1.4918, 2.2310,
+ 4.3284, 10.5836, 30.5957, 107.8530, 313.7460, 571.8840, 789.2690, 951.5410,
+ 987.1440, 991.6770, 993.3630, 999.6910, 996.6910, 983.3270, 970.8560, 952.3600,
+ 933.0200, 922.3760, 910.9220, 899.7670, 892.1120, 878.5650, 864.6450, 850.5250,
+ 841.2980, 827.2570, 808.1800, 794.0050, 774.3770, 760.2820, 746.9060, 725.8480,
+ 708.9670, 689.1840, 661.7200, 640.6300, 620.9290, 596.4570, 574.6850, 551.8710,
+ 530.9050, 511.6300, 493.6670, 474.1180, 454.3110, 437.7660, 421.4490, 406.0470,
+ 389.7410, 376.6220, 363.5170, 349.8550, 336.6970, 326.9620, 317.0800, 305.8060,
+ 294.9200, 282.2030, 271.4390, 260.8070, 250.9120, 243.3730, 235.7340, 227.6580,
+ 217.5590, 208.7220, 198.5760, 188.1520, 178.7050, 168.4410, 159.0320, 148.1250,
+ 137.7700, 124.5570, 109.6120, 91.7406, 71.5068, 50.3125, 32.8064, 19.8817,
+ 9.7957, 4.7076, 2.1849, 1.0292, .5237, .2985, .2000, .1460,
+ .1143, .1020, .0801, .0766, .0801, .0645, .0701, .0646,
+ .0538, .0298, .0300, .0277, .0007
+ };
+
+ static const float wli[9] = {0.350, 0.350, 0.350, 0.350, 0.350, 0.350, 0.350, 0.350, 0.350};
+ static const float wls[9] = {1.100, 1.100, 1.100, 1.100, 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 < 301; i++) ffu.s[40+i] = sr1[i];
+ break;
+ case 2: for(i = 0; i < 301; i++) ffu.s[40+i] = sr2[i];
+ break;
+ case 3: for(i = 0; i < 301; i++) ffu.s[40+i] = sr3[i];
+ break;
+ case 4: for(i = 0; i < 301; i++) ffu.s[40+i] = sr4[i];
+ break;
+ case 5: for(i = 0; i < 301; i++) ffu.s[40+i] = sr5[i];
+ break;
+ case 6: for(i = 0; i < 301; i++) ffu.s[40+i] = sr6[i];
+ break;
+ case 7: for(i = 0; i < 301; i++) ffu.s[40+i] = sr7[i];
+ break;
+ case 8: for(i = 0; i < 301; i++) ffu.s[40+i] = sr8[i];
+ break;
+ case 9: for(i = 0; i < 301; i++) ffu.s[40+i] = sr9[i];
+ break;
+ }
+}
+
+void IWave::quickbird2(int iwa)
+{
+
+ /* QuickBird Pan of quickbird2 */
+ static const float sr1[65] = {
+ 300.0000, 312.5000, 325.0000, 337.5000, 350.0000, 362.5000, 375.0000, 387.5000,
+ 400.0000, 412.5000, 425.0000, 437.5000, 450.0000, 462.5000, 475.0000, 487.5000,
+ 500.0000, 512.5000, 525.0000, 537.5000, 550.0000, 562.5000, 575.0000, 587.5000,
+ 600.0000, 612.5000, 625.0000, 637.5000, 650.0000, 662.5000, 675.0000, 687.5000,
+ 700.0000, 712.5000, 725.0000, 737.5000, 750.0000, 762.5000, 775.0000, 787.5000,
+ 800.0000, 812.5000, 825.0000, 837.5000, 850.0000, 862.5000, 875.0000, 887.5000,
+ 900.0000, 912.5000, 925.0000, 937.5000, 950.0000, 962.5000, 975.0000, 987.5000,
+ 1000.0000, 1012.5000, 1025.0000, 1037.5000, 1050.0000, 1062.5000, 1075.0000, 1087.5000,
+ 1100.0000
+ };
+
+ /* QuickBird Blue of quickbird2 */
+ static const float sr2[65] = {
+ .0000, .0000, .0000, .0000, .0014, .0010, .0020, .0138,
+ .0362, .0758, .1219, .1789, .2198, .2854, .3364, .3579,
+ .3781, .4276, .4982, .5925, .7051, .8105, .8911, .9426,
+ .9598, .9922, .9876, .9573, .9183, .8831, .8893, .8848,
+ .8888, .8677, .8625, .8897, .9295, .9130, .9039, .8761,
+ .8457, .8582, .8339, .8526, .8394, .7850, .7342, .6788,
+ .6186, .5832, .4969, .4305, .3738, .3035, .2567, .2162,
+ .1776, .1421, .1085, .0773, .0548, .0050, .0000, .0000,
+ .0000
+ };
+
+ /* QuickBird Green of quickbird2 */
+ static const float sr3[65] = {
+ .0000, .0000, .0010, .0037, .0053, .0030, .0013, .0023,
+ .0025, .0053, .0284, .2164, .6299, .8176, .9178, .9801,
+ .9179, .4787, .1592, .0663, .0293, .0126, .0126, .0098,
+ .0098, .0091, .0036, .0025, .0028, .0030, .0089, .0256,
+ .0083, .0086, .0091, .0028, .0031, .0034, .0036, .0033,
+ .0030, .0035, .0040, .0042, .0048, .0057, .0052, .0052,
+ .0068, .0101, .0093, .0056, .0035, .0019, .0014, .0012,
+ .0011, .0012, .0015, .0015, .0012, .0000, .0000, .0000,
+ .0000
+ };
+
+ /* QuickBird Red of quickbird2 */
+ static const float sr4[65] = {
+ .0000, .0000, .0000, .0000, .0031, .0023, .0014, .0012,
+ .0012, .0013, .0023, .0046, .0178, .0356, .0982, .2739,
+ .5277, .8005, .8705, .9640, .9934, .9733, .9652, .7230,
+ .3124, .1139, .0365, .0268, .0242, .0140, .0132, .0076,
+ .0032, .0035, .0076, .0047, .0022, .0029, .0034, .0041,
+ .0038, .0031, .0031, .0035, .0051, .0084, .0090, .0078,
+ .0076, .0081, .0067, .0043, .0027, .0016, .0013, .0010,
+ .0008, .0009, .0013, .0010, .0009, .0000, .0000, .0000,
+ .0000
+ };
+
+ /* QuickBird NIR of quickbird2 */
+ static const float sr5[65] = {
+ .0000, .0000, .0000, .0000, .0037, .0018, .0011, .0009,
+ .0007, .0013, .0010, .0024, .0029, .0032, .0032, .0023,
+ .0031, .0036, .0053, .0065, .0078, .0110, .0177, .0399,
+ .0872, .2675, .6620, .8976, .9798, .9961, .8726, .5119,
+ .1651, .0417, .0201, .0205, .0095, .0048, .0045, .0063,
+ .0045, .0026, .0025, .0025, .0025, .0024, .0022, .0019,
+ .0021, .0028, .0015, .0010, .0010, .0008, .0007, .0006,
+ .0008, .0007, .0008, .0009, .0008, .0000, .0000, .0000,
+ .0000
+ };
+
+ 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};
+
+ 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 < 65; i++) ffu.s[-99+i] = sr1[i];
+ break;
+ case 2: for(i = 0; i < 65; i++) ffu.s[-99+i] = sr2[i];
+ break;
+ case 3: for(i = 0; i < 65; i++) ffu.s[-99+i] = sr3[i];
+ break;
+ case 4: for(i = 0; i < 65; i++) ffu.s[-99+i] = sr4[i];
+ break;
+ case 5: for(i = 0; i < 65; i++) ffu.s[-99+i] = sr5[i];
+ break;
+ }
+}
+
+void IWave::landsat_8(int iwa)
+{
+
+ /* Coastal Aerosol of landsat_8 */
+ static const float sr1[13] = {
+ .0001, .0025, .0248, .3860, .9087, .9806, .9867, .9966,
+ .9828, .8257, .2264, .0256, .0024
+ };
+
+ /* Blue of landsat_8 */
+ static const float sr2[37] = {
+ .0000, .0002, .0005, .0016, .0069, .0429, .2714, .7907,
+ .9030, .9047, .8897, .8792, .8797, .8898, .8485, .8363,
+ .8685, .9115, .9317, .9549, .9564, .9838, .9895, .9681,
+ .9887, .9611, .9661, .9821, .9631, .9982, .8449, .1195,
+ .0053, .0013, .0005, .0001, .0000
+ };
+
+ /* Green of landsat_8 */
+ static const float sr3[35] = {
+ .0000, .0003, .0009, .0022, .0049, .0133, .0415, .1676,
+ .5459, .8963, .9542, .9631, .9525, .9779, .9809, .9961,
+ .9838, .9612, .9479, .9721, .9534, .9743, .9810, .9733,
+ .9758, .9744, .9678, .9717, .9837, .9603, .8093, .4353,
+ .0878, .0100, .0010
+ };
+
+ /* Red of landsat_8 */
+ static const float sr4[23] = {
+ .0003, .0027, .0145, .1077, .5268, .9267, .9474, .9689,
+ .9836, .9783, .9556, .9550, 1.0000, .9831, .9730, .9945,
+ .9640, .9690, .9499, .4629, .0460, .0053, .0006
+ };
+
+ /* Pan of landsat_8 */
+ static const float sr5[82] = {
+ .0002, .0013, .0038, .0123, .0427, .1601, .4725, .7454,
+ .8319, .8553, .8596, .8577, .8585, .8583, .8502, .8582,
+ .8615, .8577, .8792, .8917, .9063, .9129, .9029, .9207,
+ .9130, .8857, .8794, .8742, .8754, .8917, .8741, .8869,
+ .9035, .9110, .9132, .9202, .9244, .9298, .9489, .9405,
+ .9457, .9394, .9467, .9340, .9408, .9580, .9682, .9665,
+ .9572, .9477, .9525, .9575, .9642, .9674, .9770, .9760,
+ .9696, .9728, .9658, .9667, .9721, .9793, .9711, .9534,
+ .9639, .9671, .9706, .9800, .9883, .9918, 1.0000, .9985,
+ .9926, .9858, .9139, .5244, .1673, .0462, .0152, .0067,
+ .0032, .0012
+ };
+
+ /* Nir of landsat_8 */
+ static const float sr6[27] = {
+ .0000, .0002, .0005, .0012, .0032, .0090, .0283, .1149,
+ .4035, .8177, .9868, .9884, .9898, .9740, .9512, .9523,
+ .9401, .9613, .9366, .7120, .2888, .0793, .0207, .0062,
+ .0019, .0006, .0001
+ };
+
+ /* Cirrus of landsat_8 */
+ static const float sr7[25] = {
+ .0001, .0004, .0009, .0016, .0030, .0063, .0170, .0560,
+ .1828, .5269, .8624, .9465, .9387, .9933, .9832, .9638,
+ .8727, .5295, .1910, .0538, .0162, .0056, .0018, .0007,
+ .0001
+ };
+
+ /* SWIR1 of landsat_8 */
+ static const float sr8[73] = {
+ .0001, .0003, .0006, .0010, .0016, .0023, .0033, .0046,
+ .0064, .0090, .0124, .0176, .0259, .0378, .0554, .0826,
+ .1202, .1755, .2529, .3528, .4736, .6052, .7213, .8047,
+ .8734, .9035, .9187, .9270, .9251, .9241, .9214, .9237,
+ .9346, .9451, .9464, .9489, .9513, .9548, .9615, .9591,
+ .9609, .9611, .9679, .9722, .9791, .9813, .9959, .9999,
+ .9966, .9826, .9551, .8983, .7965, .6654, .5099, .3684,
+ .2513, .1657, .1067, .0694, .0450, .0291, .0187, .0124,
+ .0082, .0054, .0036, .0024, .0015, .0009, .0005, .0003,
+ .0000
+ };
+
+ /* SWIR2 of landsat_8 */
+ static const float sr9[125] = {
+ .0000, .0002, .0003, .0005, .0007, .0010, .0014, .0019,
+ .0024, .0032, .0041, .0054, .0071, .0094, .0123, .0159,
+ .0210, .0274, .0354, .0482, .0633, .0841, .1128, .1534,
+ .2028, .2673, .3425, .4315, .5225, .6134, .6897, .7665,
+ .8137, .8500, .8836, .8907, .9097, .9221, .9322, .9420,
+ .9431, .9464, .9497, .9510, .9531, .9509, .9475, .9494,
+ .9533, .9494, .9515, .9437, .9434, .9426, .9381, .9468,
+ .9474, .9473, .9341, .9390, .9347, .9264, .9369, .9295,
+ .9358, .9364, .9566, .9543, .9644, .9630, .9637, .9600,
+ .9578, .9502, .9588, .9576, .9520, .9607, .9493, .9619,
+ .9643, .9694, .9809, .9844, .9915, .9954, .9927, .9996,
+ .9944, .9860, .9891, .9765, .9742, .9757, .9737, .9584,
+ .9551, .9419, .9071, .8695, .7849, .6751, .5526, .4283,
+ .3157, .2284, .1617, .1143, .0803, .0567, .0404, .0294,
+ .0213, .0156, .0114, .0085, .0063, .0047, .0034, .0025,
+ .0018, .0012, .0008, .0005, .0002
+ };
+
+ static const float wli[9] = {0.427, 0.436, 0.513, 0.626, 0.488, 0.830, 1.341, 1.516, 2.038};
+ static const float wls[9] = {0.459, 0.527, 0.600, 0.682, 0.692, 0.896, 1.402, 1.696, 2.350};
+
+ 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 < 13; i++) ffu.s[70+i] = sr1[i];
+ break;
+ case 2: for(i = 0; i < 37; i++) ffu.s[74+i] = sr2[i];
+ break;
+ case 3: for(i = 0; i < 35; i++) ffu.s[105+i] = sr3[i];
+ break;
+ case 4: for(i = 0; i < 23; i++) ffu.s[150+i] = sr4[i];
+ break;
+ case 5: for(i = 0; i < 82; i++) ffu.s[95+i] = sr5[i];
+ break;
+ case 6: for(i = 0; i < 27; i++) ffu.s[232+i] = sr6[i];
+ break;
+ case 7: for(i = 0; i < 25; i++) ffu.s[436+i] = sr7[i];
+ break;
+ case 8: for(i = 0; i < 73; i++) ffu.s[506+i] = sr8[i];
+ break;
+ case 9: for(i = 0; i < 125; i++) ffu.s[715+i] = sr9[i];
+ break;
+ }
+}
+
+
+
float IWave::equivwl() const
{
float seb = 0;
@@ -2529,6 +3167,9 @@
else if(iwave <= 92) rapideye(iwave - 87);
else if(iwave <= 96) vgt1_spot4(iwave - 92);
else if(iwave <= 100) vgt2_spot5(iwave - 96);
+ else if(iwave <= 109) worldview2(iwave - 100);
+ else if(iwave <= 114) quickbird2(iwave - 109);
+ else if(iwave <= 123) landsat_8(iwave - 114);
else G_warning(_("Unsupported iwave value: %d"), iwave);
}
@@ -2545,7 +3186,8 @@
/* --- spectral condition ---- */
void IWave::print()
{
- static const string nsat[102] = {
+/* The compiler will deduce the dimension from the initializer list. */
+ static const string nsat[] = {
string(" constant "), string(" user s "),
string(" meteosat "), string(" goes east "), string(" goes west "),
string(" avhrr 1 (noaa6) "), string(" avhrr 2 (noaa6) "),
@@ -2588,7 +3230,15 @@
string(" vgt1 spot4 0 "), string(" vgt1 spot4 2 "), string(" vgt1 spot4 3 "),
string(" vgt1 spot4 MIR "),
string(" vgt2 spot5 0 "), string(" vgt2 spot5 2 "), string(" vgt2 spot5 3 "),
- string(" vgt2 spot5 MIR "),
+ string(" vgt2 spot5 MIR "), string(" worldview 2 pan "), string(" worldview 2 coastal blue"),
+ string(" worldview 2 blue "), string(" worldview 2 blue "), string(" worldview 2 green"),
+ string(" worldview 2 yellow"), string(" worldview 2 red"), string(" worldview 2 red edge"),
+ string(" worldview 2 nir1"), string(" worldview 2 nir2"), string(" quickbird 2 pan"),
+ string(" quickbird 2 blue"), string(" quickbird 2 green"), string(" quickbird 2 red"),
+ string(" quickbird 2 nir"), string(" Landsat 8 Coastal Aerosol"), string(" Landsat 8 Blue"),
+ string(" Landsat 8 Green"),string(" Landsat 8 Red"), string(" Landsat 8 Panchromatic"),
+ string(" Landsat 8 NIR"),string(" Landsat 8 Cirrus"),string(" Landsat 8 SWIR1"),
+ string(" Landsat 8 SWIR2")
};
Modified: grass/branches/releasebranch_7_0/imagery/i.atcorr/iwave.h
===================================================================
--- grass/branches/releasebranch_7_0/imagery/i.atcorr/iwave.h 2014-07-22 09:48:50 UTC (rev 61313)
+++ grass/branches/releasebranch_7_0/imagery/i.atcorr/iwave.h 2014-07-22 10:57:03 UTC (rev 61314)
@@ -122,6 +122,30 @@
c 98 B2 " ( 0.580-0.780 ) c
c 99 B3 " ( 0.700-1.000 ) c
c 100 MIR " ( 1.450-1.800 ) c
+c 101 Pan band of Worldview 2 ( 0.450-0.800 ) c
+c 102 Coastal " ( 0.396-0.458 ) c
+c 103 Blue " ( 0.442-0.515 ) c
+c 104 Green " ( 0.506-0.586 ) c
+c 105 Yellow " ( 0.584-0.632 ) c
+c 106 Red " ( 0.624-0.694 ) c
+c 107 RedEdge " ( 0.699-0.749 ) c
+c 108 NIR1 " ( 0.765-0.901 ) c
+c 109 NIR2 " ( 0.856-1.043 ) c
+c 110 Pan band of QuickBird 2 ( 0.405-1.053 ) c
+c 111 Blue " ( 0.430-0.545 ) c
+c 112 Green " ( 0.466-0.620 ) c
+c 113 Red " ( 0.590-0.710 ) c
+c 114 NIR " ( 0.715-0.918 ) c
+c 115 Coastal Aerosol of LS8 ( 0.427-0.459 ) c
+c 116 Blue " ( 0.436-0.527 ) c
+c 117 Green " ( 0.512-0.610 ) c
+c 118 Red " ( 0.625-0.691 ) c
+c 119 Pan " ( 0.488-0.692 ) c
+c 120 NIR " ( 0.829-0.900 ) c
+c 121 Cirrus " ( 1.340-1.409 ) c
+c 122 SWIR1 " ( 1.515-1.697 ) c
+c 123 SWIR2 " ( 2.037-2.037 ) c
+c c
c note: wl has to be in micrometer c
c**********************************************************************/
@@ -163,8 +187,10 @@
void rapideye(int iwa);
void vgt1_spot4(int iwa);
void vgt2_spot5(int iwa);
+ void worldview2(int iwa);
+ void quickbird2(int iwa);
+ void landsat_8(int iwa);
-
public:
/* To compute the equivalent wavelength needed for the calculation of the
downward radiation field used in the computation of the non lambertian
Added: grass/branches/releasebranch_7_0/imagery/i.atcorr/sensors_csv/landsat_8.csv
===================================================================
--- grass/branches/releasebranch_7_0/imagery/i.atcorr/sensors_csv/landsat_8.csv (rev 0)
+++ grass/branches/releasebranch_7_0/imagery/i.atcorr/sensors_csv/landsat_8.csv 2014-07-22 10:57:03 UTC (rev 61314)
@@ -0,0 +1,1930 @@
+nm,Coastal Aerosol,Blue,Green,Red,Pan,Nir,Cirrus,SWIR1,SWIR2
+427,0.000073,,,,,,,,
+428,0.000609,,,,,,,,
+429,0.001628,,,,,,,,
+430,0.003421,,,,,,,,
+431,0.008019,,,,,,,,
+432,0.024767,,,,,,,,
+433,0.085688,,,,,,,,
+434,0.254149,,,,,,,,
+435,0.517821,,,,,,,,
+436,0.765117,0.00001,,,,,,,
+437,0.908749,0.000061,,,,,,,
+438,0.958204,0.000117,,,,,,,
+439,0.977393,0.000241,,,,,,,
+440,0.98379,0.000349,,,,,,,
+441,0.989052,0.000455,,,,,,,
+442,0.986713,0.000756,,,,,,,
+443,0.993683,0.001197,,,,,,,
+444,0.993137,0.00207,,,,,,,
+445,1,0.003712,,,,,,,
+446,0.996969,0.006869,,,,,,,
+447,0.98278,0.013212,,,,,,,
+448,0.972692,0.02717,,,,,,,
+449,0.905808,0.058606,,,,,,,
+450,0.745606,0.130876,,,,,,,
+451,0.471329,0.27137,,,,,,,
+452,0.226412,0.493542,,,,,,,
+453,0.09286,0.723971,,,,,,,
+454,0.036603,0.85751,,,,,,,
+455,0.014537,0.894222,,,,,,,
+456,0.005829,0.903034,,,,,,,
+457,0.002414,0.910928,,,,,,,
+458,0.000984,0.90988,,,,,,,
+459,0.000255,0.899475,,,,,,,
+460,,0.897977,,,,,,,
+461,,0.889667,,,,,,,
+462,,0.883542,,,,,,,
+463,,0.877453,,,,,,,
+464,,0.881011,,,,,,,
+465,,0.874721,,,,,,,
+466,,0.879688,,,,,,,
+467,,0.886569,,,,,,,
+468,,0.891913,,,,,,,
+469,,0.88768,,,,,,,
+470,,0.861157,,,,,,,
+471,,0.848533,,,,,,,
+472,,0.840828,,,,,,,
+473,,0.828339,,,,,,,
+474,,0.844202,,,,,,,
+475,,0.865864,,,,,,,
+476,,0.868497,,,,,,,
+477,,0.890253,,,,,,,
+478,,0.912538,,,,,,,
+479,,0.910385,,,,,,,
+480,,0.918822,,,,,,,
+481,,0.931726,,,,,,,
+482,,0.931813,,,,,,,
+483,,0.954248,,,,,,,
+484,,0.955545,,,,,,,
+485,,0.96242,,,,,,,
+486,,0.956424,,,,,,,
+487,,0.953352,,,,,,,
+488,,0.978564,,,0.000216,,,,
+489,,0.989104,,,0.000514,,,,
+490,,0.985615,,,0.001013,,,,
+491,,0.989469,,,0.001587,,,,
+492,,0.982262,,,0.002619,,,,
+493,,0.968801,,,0.003841,,,,
+494,,0.967332,,,0.006195,,,,
+495,,0.976836,,,0.009003,,,,
+496,,0.988729,,,0.015515,,,,
+497,,0.980826,,,0.02327,,,,
+498,,0.967361,,,0.042723,,,,
+499,,0.954754,,,0.066556,,,,
+500,,0.964132,,,0.124269,,,,
+501,,0.966125,,,0.196006,,,,
+502,,0.966772,,,0.323232,,,,
+503,,0.981834,,,0.472496,,,,
+504,,0.98232,,,0.598302,,,,
+505,,0.965685,,,0.714816,,,,
+506,,0.963135,,,0.776009,,,,
+507,,0.972261,,,0.816555,,,,
+508,,0.996498,,,0.831881,,,,
+509,,1,,,0.837994,,,,
+510,,0.9556,,,0.848975,,,,
+511,,0.844893,,,0.861668,,,,
+512,,0.534592,-0.000046,,0.862705,,,,
+513,,0.190738,0.000016,,0.85964,,,,
+514,,0.048329,0.00011,,0.858532,,,,
+515,,0.013894,0.000247,,0.858244,,,,
+516,,0.005328,0.000362,,0.857148,,,,
+517,,0.002611,0.000648,,0.855759,,,,
+518,,0.001557,0.000935,,0.858455,,,,
+519,,0.0011,0.001332,,0.862834,,,,
+520,,0.000785,0.001816,,0.860609,,,,
+521,,0.000516,0.002515,,0.855993,,,,
+522,,0.000321,0.003446,,0.852785,,,,
+523,,0.000162,0.00488,,0.850183,,,,
+524,,0.000072,0.007024,,0.852788,,,,
+525,,0.000057,0.010441,,0.857317,,,,
+526,,0.000023,0.016247,,0.85913,,,,
+527,,0.000032,0.025513,,0.859876,,,,
+528,,-0.000016,0.041451,,0.861508,,,,
+529,,,0.070551,,0.863494,,,,
+530,,,0.123444,,0.860265,,,,
+531,,,0.21168,,0.855102,,,,
+532,,,0.353885,,0.864423,,,,
+533,,,0.545856,,0.879249,,,,
+534,,,0.741205,,0.885593,,,,
+535,,,0.865225,,0.888663,,,,
+536,,,0.927396,,0.894758,,,,
+537,,,0.954627,,0.902017,,,,
+538,,,0.954163,,0.906294,,,,
+539,,,0.959215,,0.90929,,,,
+540,,,0.961328,,0.911734,,,,
+541,,,0.964902,,0.914,,,,
+542,,,0.969873,,0.909776,,,,
+543,,,0.952489,,0.902939,,,,
+544,,,0.961397,,0.908919,,,,
+545,,,0.97827,,0.919859,,,,
+546,,,0.977533,,0.92162,,,,
+547,,,0.977001,,0.919456,,,,
+548,,,0.980884,,0.91302,,,,
+549,,,0.990784,,0.905134,,,,
+550,,,1,,0.892702,,,,
+551,,,0.992264,,0.878599,,,,
+552,,,0.982642,,0.876838,,,,
+553,,,0.983832,,0.879443,,,,
+554,,,0.977765,,0.878303,,,,
+555,,,0.965081,,0.875692,,,,
+556,,,0.957314,,0.872666,,,,
+557,,,0.946245,,0.869478,,,,
+558,,,0.947871,,0.875361,,,,
+559,,,0.959038,,0.884759,,,,
+560,,,0.966534,,0.889903,,,,
+561,,,0.977656,,0.893427,,,,
+562,,,0.966447,,0.885948,,,,
+563,,,0.953399,,0.874097,,,,
+564,,,0.958314,,0.875674,,,,
+565,,,0.970039,,0.882402,,,,
+566,,,0.978607,,0.891375,,,,
+567,,,0.983397,,0.901194,,,,
+568,,,0.98096,,0.903528,,,,
+569,,,0.974522,,0.902826,,,,
+570,,,0.967229,,0.907551,,,,
+571,,,0.979406,,0.91435,,,,
+572,,,0.978208,,0.914932,,,,
+573,,,0.975818,,0.91319,,,,
+574,,,0.974392,,0.915508,,,,
+575,,,0.979973,,0.919554,,,,
+576,,,0.968827,,0.920802,,,,
+577,,,0.969181,,0.92104,,,,
+578,,,0.967838,,0.924431,,,,
+579,,,0.982956,,0.929105,,,,
+580,,,0.979598,,0.930068,,,,
+581,,,0.963811,,0.929551,,,,
+582,,,0.968886,,0.937549,,,,
+583,,,0.983655,,0.948863,,,,
+584,,,0.986657,,0.947084,,,,
+585,,,0.974207,,0.940273,,,,
+586,,,0.946407,,0.940813,,,,
+587,,,0.904478,,0.944136,,,,
+588,,,0.809275,,0.945674,,,,
+589,,,0.684974,,0.946566,,,,
+590,,,0.525304,,0.942462,,,,
+591,,,0.345364,,0.936298,,,,
+592,,,0.190467,,0.939621,,,,
+593,,,0.087833,,0.946659,,,,
+594,,,0.035393,,0.942374,,,,
+595,,,0.014077,,0.933839,,,,
+596,,,0.005944,,0.93425,,,,
+597,,,0.002574,,0.93784,,,,
+598,,,0.001046,,0.940838,,,,
+599,,,0.000394,,0.943626,,,,
+600,,,0.000085,,0.952415,,,,
+601,,,-0.000084,,0.963664,,,,
+602,,,-0.000194,,0.967381,,,,
+603,,,-0.000234,,0.968241,,,,
+604,,,-0.000292,,0.968233,,,,
+605,,,-0.000307,,0.967878,,,,
+606,,,-0.000343,,0.965083,,,,
+607,,,-0.000348,,0.961367,,,,
+608,,,-0.000329,,0.957232,,,,
+609,,,-0.000351,,0.952939,,,,
+610,,,-0.000317,,0.949345,,,,
+611,,,,,0.946006,,,,
+612,,,,,0.948081,,,,
+613,,,,,0.952465,,,,
+614,,,,,0.954526,,,,
+615,,,,,0.955814,,,,
+616,,,,,0.959149,,,,
+617,,,,,0.963216,,,,
+618,,,,,0.964158,,,,
+619,,,,,0.963953,,,,
+620,,,,,0.965943,,,,
+621,,,,,0.968789,,,,
+622,,,,,0.972688,,,,
+623,,,,,0.977026,,,,
+624,,,,,0.978028,,,,
+625,,,,-0.000342,0.977901,,,,
+626,,,,0.00027,0.974158,,,,
+627,,,,0.000895,0.969078,,,,
+628,,,,0.00185,0.969583,,,,
+629,,,,0.003648,0.972297,,,,
+630,,,,0.007197,0.972906,,,,
+631,,,,0.014515,0.972708,,,,
+632,,,,0.030432,0.969732,,,,
+633,,,,0.066861,0.96578,,,,
+634,,,,0.148518,0.965561,,,,
+635,,,,0.299778,0.966835,,,,
+636,,,,0.526812,0.966641,,,,
+637,,,,0.764443,0.965899,,,,
+638,,,,0.905473,0.972067,,,,
+639,,,,0.947949,0.981077,,,,
+640,,,,0.950823,0.980971,,,,
+641,,,,0.947418,0.977722,,,,
+642,,,,0.951831,0.974405,,,,
+643,,,,0.962705,0.971123,,,,
+644,,,,0.975075,0.963139,,,,
+645,,,,0.984173,0.95326,,,,
+646,,,,0.983613,0.953494,,,,
+647,,,,0.983434,0.957536,,,,
+648,,,,0.982911,0.963851,,,,
+649,,,,0.973636,0.970928,,,,
+650,,,,0.959441,0.969453,,,,
+651,,,,0.955641,0.96475,,,,
+652,,,,0.955548,0.966541,,,,
+653,,,,0.953337,0.970613,,,,
+654,,,,0.956628,0.973941,,,,
+655,,,,0.981688,0.977055,,,,
+656,,,,1,0.982894,,,,
+657,,,,0.992388,0.98977,,,,
+658,,,,0.984615,0.988302,,,,
+659,,,,0.981568,0.983542,,,,
+660,,,,0.97696,0.987968,,,,
+661,,,,0.97298,0.995538,,,,
+662,,,,0.98108,0.998579,,,,
+663,,,,0.996804,1,,,,
+664,,,,0.992142,0.99981,,,,
+665,,,,0.980678,0.999061,,,,
+666,,,,0.964002,0.997891,,,,
+667,,,,0.962154,0.99658,,,,
+668,,,,0.970778,0.992555,,,,
+669,,,,0.96718,0.987369,,,,
+670,,,,0.966928,0.985862,,,,
+671,,,,0.949928,0.985761,,,,
+672,,,,0.848855,0.955679,,,,
+673,,,,0.609359,0.913945,,,,
+674,,,,0.31635,0.78014,,,,
+675,,,,0.123946,0.610197,,,,
+676,,,,0.046033,0.438556,,,,
+677,,,,0.017702,0.266714,,,,
+678,,,,0.007333,0.167313,,,,
+679,,,,0.003205,0.094013,,,,
+680,,,,0.001402,0.057564,,,,
+681,,,,0.000554,0.034787,,,,
+682,,,,0.000117,0.022882,,,,
+683,,,,-0.000122,0.015178,,,,
+684,,,,-0.000289,0.010755,,,,
+685,,,,-0.000376,0.007663,,,,
+686,,,,-0.000396,0.005756,,,,
+687,,,,-0.000458,0.004292,,,,
+688,,,,-0.000488,0.00322,,,,
+689,,,,-0.000457,0.002283,,,,
+690,,,,-0.000429,0.001545,,,,
+691,,,,-0.000417,0.000879,,,,
+692,,,,,0.000214,,,,
+693,,,,,,,,,
+694,,,,,,,,,
+695,,,,,,,,,
+696,,,,,,,,,
+697,,,,,,,,,
+698,,,,,,,,,
+699,,,,,,,,,
+700,,,,,,,,,
+701,,,,,,,,,
+702,,,,,,,,,
+703,,,,,,,,,
+704,,,,,,,,,
+705,,,,,,,,,
+706,,,,,,,,,
+707,,,,,,,,,
+708,,,,,,,,,
+709,,,,,,,,,
+710,,,,,,,,,
+711,,,,,,,,,
+712,,,,,,,,,
+713,,,,,,,,,
+714,,,,,,,,,
+715,,,,,,,,,
+716,,,,,,,,,
+717,,,,,,,,,
+718,,,,,,,,,
+719,,,,,,,,,
+720,,,,,,,,,
+721,,,,,,,,,
+722,,,,,,,,,
+723,,,,,,,,,
+724,,,,,,,,,
+725,,,,,,,,,
+726,,,,,,,,,
+727,,,,,,,,,
+728,,,,,,,,,
+729,,,,,,,,,
+730,,,,,,,,,
+731,,,,,,,,,
+732,,,,,,,,,
+733,,,,,,,,,
+734,,,,,,,,,
+735,,,,,,,,,
+736,,,,,,,,,
+737,,,,,,,,,
+738,,,,,,,,,
+739,,,,,,,,,
+740,,,,,,,,,
+741,,,,,,,,,
+742,,,,,,,,,
+743,,,,,,,,,
+744,,,,,,,,,
+745,,,,,,,,,
+746,,,,,,,,,
+747,,,,,,,,,
+748,,,,,,,,,
+749,,,,,,,,,
+750,,,,,,,,,
+751,,,,,,,,,
+752,,,,,,,,,
+753,,,,,,,,,
+754,,,,,,,,,
+755,,,,,,,,,
+756,,,,,,,,,
+757,,,,,,,,,
+758,,,,,,,,,
+759,,,,,,,,,
+760,,,,,,,,,
+761,,,,,,,,,
+762,,,,,,,,,
+763,,,,,,,,,
+764,,,,,,,,,
+765,,,,,,,,,
+766,,,,,,,,,
+767,,,,,,,,,
+768,,,,,,,,,
+769,,,,,,,,,
+770,,,,,,,,,
+771,,,,,,,,,
+772,,,,,,,,,
+773,,,,,,,,,
+774,,,,,,,,,
+775,,,,,,,,,
+776,,,,,,,,,
+777,,,,,,,,,
+778,,,,,,,,,
+779,,,,,,,,,
+780,,,,,,,,,
+781,,,,,,,,,
+782,,,,,,,,,
+783,,,,,,,,,
+784,,,,,,,,,
+785,,,,,,,,,
+786,,,,,,,,,
+787,,,,,,,,,
+788,,,,,,,,,
+789,,,,,,,,,
+790,,,,,,,,,
+791,,,,,,,,,
+792,,,,,,,,,
+793,,,,,,,,,
+794,,,,,,,,,
+795,,,,,,,,,
+796,,,,,,,,,
+797,,,,,,,,,
+798,,,,,,,,,
+799,,,,,,,,,
+800,,,,,,,,,
+801,,,,,,,,,
+802,,,,,,,,,
+803,,,,,,,,,
+804,,,,,,,,,
+805,,,,,,,,,
+806,,,,,,,,,
+807,,,,,,,,,
+808,,,,,,,,,
+809,,,,,,,,,
+810,,,,,,,,,
+811,,,,,,,,,
+812,,,,,,,,,
+813,,,,,,,,,
+814,,,,,,,,,
+815,,,,,,,,,
+816,,,,,,,,,
+817,,,,,,,,,
+818,,,,,,,,,
+819,,,,,,,,,
+820,,,,,,,,,
+821,,,,,,,,,
+822,,,,,,,,,
+823,,,,,,,,,
+824,,,,,,,,,
+825,,,,,,,,,
+826,,,,,,,,,
+827,,,,,,,,,
+828,,,,,,,,,
+829,,,,,,-0.000034,,,
+830,,,,,,0.000011,,,
+831,,,,,,0.00005,,,
+832,,,,,,0.0001,,,
+833,,,,,,0.000239,,,
+834,,,,,,0.000314,,,
+835,,,,,,0.000495,,,
+836,,,,,,0.000719,,,
+837,,,,,,0.000986,,,
+838,,,,,,0.001445,,,
+839,,,,,,0.002107,,,
+840,,,,,,0.00316,,,
+841,,,,,,0.004744,,,
+842,,,,,,0.007059,,,
+843,,,,,,0.0109,,,
+844,,,,,,0.017346,,,
+845,,,,,,0.028332,,,
+846,,,,,,0.048191,,,
+847,,,,,,0.084363,,,
+848,,,,,,0.145365,,,
+849,,,,,,0.249733,,,
+850,,,,,,0.403526,,,
+851,,,,,,0.582623,,,
+852,,,,,,0.745037,,,
+853,,,,,,0.890315,,,
+854,,,,,,0.960215,,,
+855,,,,,,0.986833,,,
+856,,,,,,0.973133,,,
+857,,,,,,0.980606,,,
+858,,,,,,0.99612,,,
+859,,,,,,1,,,
+860,,,,,,0.989777,,,
+861,,,,,,0.980733,,,
+862,,,,,,0.975935,,,
+863,,,,,,0.972043,,,
+864,,,,,,0.957357,,,
+865,,,,,,0.951209,,,
+866,,,,,,0.947044,,,
+867,,,,,,0.953162,,,
+868,,,,,,0.951499,,,
+869,,,,,,0.94845,,,
+870,,,,,,0.940094,,,
+871,,,,,,0.950632,,,
+872,,,,,,0.956079,,,
+873,,,,,,0.96646,,,
+874,,,,,,0.969821,,,
+875,,,,,,0.93661,,,
+876,,,,,,0.891066,,,
+877,,,,,,0.788733,,,
+878,,,,,,0.63532,,,
+879,,,,,,0.448364,,,
+880,,,,,,0.288847,,,
+881,,,,,,0.174619,,,
+882,,,,,,0.100343,,,
+883,,,,,,0.058265,,,
+884,,,,,,0.034532,,,
+885,,,,,,0.02072,,,
+886,,,,,,0.01244,,,
+887,,,,,,0.007601,,,
+888,,,,,,0.004702,,,
+889,,,,,,0.002944,,,
+890,,,,,,0.00187,,,
+891,,,,,,0.001192,,,
+892,,,,,,0.000743,,,
+893,,,,,,0.000423,,,
+894,,,,,,0.000241,,,
+895,,,,,,0.000116,,,
+896,,,,,,0.000044,,,
+897,,,,,,-0.000013,,,
+898,,,,,,-0.00005,,,
+899,,,,,,-0.000084,,,
+900,,,,,,-0.0001,,,
+901,,,,,,,,,
+902,,,,,,,,,
+903,,,,,,,,,
+904,,,,,,,,,
+905,,,,,,,,,
+906,,,,,,,,,
+907,,,,,,,,,
+908,,,,,,,,,
+909,,,,,,,,,
+910,,,,,,,,,
+911,,,,,,,,,
+912,,,,,,,,,
+913,,,,,,,,,
+914,,,,,,,,,
+915,,,,,,,,,
+916,,,,,,,,,
+917,,,,,,,,,
+918,,,,,,,,,
+919,,,,,,,,,
+920,,,,,,,,,
+921,,,,,,,,,
+922,,,,,,,,,
+923,,,,,,,,,
+924,,,,,,,,,
+925,,,,,,,,,
+926,,,,,,,,,
+927,,,,,,,,,
+928,,,,,,,,,
+929,,,,,,,,,
+930,,,,,,,,,
+931,,,,,,,,,
+932,,,,,,,,,
+933,,,,,,,,,
+934,,,,,,,,,
+935,,,,,,,,,
+936,,,,,,,,,
+937,,,,,,,,,
+938,,,,,,,,,
+939,,,,,,,,,
+940,,,,,,,,,
+941,,,,,,,,,
+942,,,,,,,,,
+943,,,,,,,,,
+944,,,,,,,,,
+945,,,,,,,,,
+946,,,,,,,,,
+947,,,,,,,,,
+948,,,,,,,,,
+949,,,,,,,,,
+950,,,,,,,,,
+951,,,,,,,,,
+952,,,,,,,,,
+953,,,,,,,,,
+954,,,,,,,,,
+955,,,,,,,,,
+956,,,,,,,,,
+957,,,,,,,,,
+958,,,,,,,,,
+959,,,,,,,,,
+960,,,,,,,,,
+961,,,,,,,,,
+962,,,,,,,,,
+963,,,,,,,,,
+964,,,,,,,,,
+965,,,,,,,,,
+966,,,,,,,,,
+967,,,,,,,,,
+968,,,,,,,,,
+969,,,,,,,,,
+970,,,,,,,,,
+971,,,,,,,,,
+972,,,,,,,,,
+973,,,,,,,,,
+974,,,,,,,,,
+975,,,,,,,,,
+976,,,,,,,,,
+977,,,,,,,,,
+978,,,,,,,,,
+979,,,,,,,,,
+980,,,,,,,,,
+981,,,,,,,,,
+982,,,,,,,,,
+983,,,,,,,,,
+984,,,,,,,,,
+985,,,,,,,,,
+986,,,,,,,,,
+987,,,,,,,,,
+988,,,,,,,,,
+989,,,,,,,,,
+990,,,,,,,,,
+991,,,,,,,,,
+992,,,,,,,,,
+993,,,,,,,,,
+994,,,,,,,,,
+995,,,,,,,,,
+996,,,,,,,,,
+997,,,,,,,,,
+998,,,,,,,,,
+999,,,,,,,,,
+1000,,,,,,,,,
+1001,,,,,,,,,
+1002,,,,,,,,,
+1003,,,,,,,,,
+1004,,,,,,,,,
+1005,,,,,,,,,
+1006,,,,,,,,,
+1007,,,,,,,,,
+1008,,,,,,,,,
+1009,,,,,,,,,
+1010,,,,,,,,,
+1011,,,,,,,,,
+1012,,,,,,,,,
+1013,,,,,,,,,
+1014,,,,,,,,,
+1015,,,,,,,,,
+1016,,,,,,,,,
+1017,,,,,,,,,
+1018,,,,,,,,,
+1019,,,,,,,,,
+1020,,,,,,,,,
+1021,,,,,,,,,
+1022,,,,,,,,,
+1023,,,,,,,,,
+1024,,,,,,,,,
+1025,,,,,,,,,
+1026,,,,,,,,,
+1027,,,,,,,,,
+1028,,,,,,,,,
+1029,,,,,,,,,
+1030,,,,,,,,,
+1031,,,,,,,,,
+1032,,,,,,,,,
+1033,,,,,,,,,
+1034,,,,,,,,,
+1035,,,,,,,,,
+1036,,,,,,,,,
+1037,,,,,,,,,
+1038,,,,,,,,,
+1039,,,,,,,,,
+1040,,,,,,,,,
+1041,,,,,,,,,
+1042,,,,,,,,,
+1043,,,,,,,,,
+1044,,,,,,,,,
+1045,,,,,,,,,
+1046,,,,,,,,,
+1047,,,,,,,,,
+1048,,,,,,,,,
+1049,,,,,,,,,
+1050,,,,,,,,,
+1051,,,,,,,,,
+1052,,,,,,,,,
+1053,,,,,,,,,
+1054,,,,,,,,,
+1055,,,,,,,,,
+1056,,,,,,,,,
+1057,,,,,,,,,
+1058,,,,,,,,,
+1059,,,,,,,,,
+1060,,,,,,,,,
+1061,,,,,,,,,
+1062,,,,,,,,,
+1063,,,,,,,,,
+1064,,,,,,,,,
+1065,,,,,,,,,
+1066,,,,,,,,,
+1067,,,,,,,,,
+1068,,,,,,,,,
+1069,,,,,,,,,
+1070,,,,,,,,,
+1071,,,,,,,,,
+1072,,,,,,,,,
+1073,,,,,,,,,
+1074,,,,,,,,,
+1075,,,,,,,,,
+1076,,,,,,,,,
+1077,,,,,,,,,
+1078,,,,,,,,,
+1079,,,,,,,,,
+1080,,,,,,,,,
+1081,,,,,,,,,
+1082,,,,,,,,,
+1083,,,,,,,,,
+1084,,,,,,,,,
+1085,,,,,,,,,
+1086,,,,,,,,,
+1087,,,,,,,,,
+1088,,,,,,,,,
+1089,,,,,,,,,
+1090,,,,,,,,,
+1091,,,,,,,,,
+1092,,,,,,,,,
+1093,,,,,,,,,
+1094,,,,,,,,,
+1095,,,,,,,,,
+1096,,,,,,,,,
+1097,,,,,,,,,
+1098,,,,,,,,,
+1099,,,,,,,,,
+1100,,,,,,,,,
+1101,,,,,,,,,
+1102,,,,,,,,,
+1103,,,,,,,,,
+1104,,,,,,,,,
+1105,,,,,,,,,
+1106,,,,,,,,,
+1107,,,,,,,,,
+1108,,,,,,,,,
+1109,,,,,,,,,
+1110,,,,,,,,,
+1111,,,,,,,,,
+1112,,,,,,,,,
+1113,,,,,,,,,
+1114,,,,,,,,,
+1115,,,,,,,,,
+1116,,,,,,,,,
+1117,,,,,,,,,
+1118,,,,,,,,,
+1119,,,,,,,,,
+1120,,,,,,,,,
+1121,,,,,,,,,
+1122,,,,,,,,,
+1123,,,,,,,,,
+1124,,,,,,,,,
+1125,,,,,,,,,
+1126,,,,,,,,,
+1127,,,,,,,,,
+1128,,,,,,,,,
+1129,,,,,,,,,
+1130,,,,,,,,,
+1131,,,,,,,,,
+1132,,,,,,,,,
+1133,,,,,,,,,
+1134,,,,,,,,,
+1135,,,,,,,,,
+1136,,,,,,,,,
+1137,,,,,,,,,
+1138,,,,,,,,,
+1139,,,,,,,,,
+1140,,,,,,,,,
+1141,,,,,,,,,
+1142,,,,,,,,,
+1143,,,,,,,,,
+1144,,,,,,,,,
+1145,,,,,,,,,
+1146,,,,,,,,,
+1147,,,,,,,,,
+1148,,,,,,,,,
+1149,,,,,,,,,
+1150,,,,,,,,,
+1151,,,,,,,,,
+1152,,,,,,,,,
+1153,,,,,,,,,
+1154,,,,,,,,,
+1155,,,,,,,,,
+1156,,,,,,,,,
+1157,,,,,,,,,
+1158,,,,,,,,,
+1159,,,,,,,,,
+1160,,,,,,,,,
+1161,,,,,,,,,
+1162,,,,,,,,,
+1163,,,,,,,,,
+1164,,,,,,,,,
+1165,,,,,,,,,
+1166,,,,,,,,,
+1167,,,,,,,,,
+1168,,,,,,,,,
+1169,,,,,,,,,
+1170,,,,,,,,,
+1171,,,,,,,,,
+1172,,,,,,,,,
+1173,,,,,,,,,
+1174,,,,,,,,,
+1175,,,,,,,,,
+1176,,,,,,,,,
+1177,,,,,,,,,
+1178,,,,,,,,,
+1179,,,,,,,,,
+1180,,,,,,,,,
+1181,,,,,,,,,
+1182,,,,,,,,,
+1183,,,,,,,,,
+1184,,,,,,,,,
+1185,,,,,,,,,
+1186,,,,,,,,,
+1187,,,,,,,,,
+1188,,,,,,,,,
+1189,,,,,,,,,
+1190,,,,,,,,,
+1191,,,,,,,,,
+1192,,,,,,,,,
+1193,,,,,,,,,
+1194,,,,,,,,,
+1195,,,,,,,,,
+1196,,,,,,,,,
+1197,,,,,,,,,
+1198,,,,,,,,,
+1199,,,,,,,,,
+1200,,,,,,,,,
+1201,,,,,,,,,
+1202,,,,,,,,,
+1203,,,,,,,,,
+1204,,,,,,,,,
+1205,,,,,,,,,
+1206,,,,,,,,,
+1207,,,,,,,,,
+1208,,,,,,,,,
+1209,,,,,,,,,
+1210,,,,,,,,,
+1211,,,,,,,,,
+1212,,,,,,,,,
+1213,,,,,,,,,
+1214,,,,,,,,,
+1215,,,,,,,,,
+1216,,,,,,,,,
+1217,,,,,,,,,
+1218,,,,,,,,,
+1219,,,,,,,,,
+1220,,,,,,,,,
+1221,,,,,,,,,
+1222,,,,,,,,,
+1223,,,,,,,,,
+1224,,,,,,,,,
+1225,,,,,,,,,
+1226,,,,,,,,,
+1227,,,,,,,,,
+1228,,,,,,,,,
+1229,,,,,,,,,
+1230,,,,,,,,,
+1231,,,,,,,,,
+1232,,,,,,,,,
+1233,,,,,,,,,
+1234,,,,,,,,,
+1235,,,,,,,,,
+1236,,,,,,,,,
+1237,,,,,,,,,
+1238,,,,,,,,,
+1239,,,,,,,,,
+1240,,,,,,,,,
+1241,,,,,,,,,
+1242,,,,,,,,,
+1243,,,,,,,,,
+1244,,,,,,,,,
+1245,,,,,,,,,
+1246,,,,,,,,,
+1247,,,,,,,,,
+1248,,,,,,,,,
+1249,,,,,,,,,
+1250,,,,,,,,,
+1251,,,,,,,,,
+1252,,,,,,,,,
+1253,,,,,,,,,
+1254,,,,,,,,,
+1255,,,,,,,,,
+1256,,,,,,,,,
+1257,,,,,,,,,
+1258,,,,,,,,,
+1259,,,,,,,,,
+1260,,,,,,,,,
+1261,,,,,,,,,
+1262,,,,,,,,,
+1263,,,,,,,,,
+1264,,,,,,,,,
+1265,,,,,,,,,
+1266,,,,,,,,,
+1267,,,,,,,,,
+1268,,,,,,,,,
+1269,,,,,,,,,
+1270,,,,,,,,,
+1271,,,,,,,,,
+1272,,,,,,,,,
+1273,,,,,,,,,
+1274,,,,,,,,,
+1275,,,,,,,,,
+1276,,,,,,,,,
+1277,,,,,,,,,
+1278,,,,,,,,,
+1279,,,,,,,,,
+1280,,,,,,,,,
+1281,,,,,,,,,
+1282,,,,,,,,,
+1283,,,,,,,,,
+1284,,,,,,,,,
+1285,,,,,,,,,
+1286,,,,,,,,,
+1287,,,,,,,,,
+1288,,,,,,,,,
+1289,,,,,,,,,
+1290,,,,,,,,,
+1291,,,,,,,,,
+1292,,,,,,,,,
+1293,,,,,,,,,
+1294,,,,,,,,,
+1295,,,,,,,,,
+1296,,,,,,,,,
+1297,,,,,,,,,
+1298,,,,,,,,,
+1299,,,,,,,,,
+1300,,,,,,,,,
+1301,,,,,,,,,
+1302,,,,,,,,,
+1303,,,,,,,,,
+1304,,,,,,,,,
+1305,,,,,,,,,
+1306,,,,,,,,,
+1307,,,,,,,,,
+1308,,,,,,,,,
+1309,,,,,,,,,
+1310,,,,,,,,,
+1311,,,,,,,,,
+1312,,,,,,,,,
+1313,,,,,,,,,
+1314,,,,,,,,,
+1315,,,,,,,,,
+1316,,,,,,,,,
+1317,,,,,,,,,
+1318,,,,,,,,,
+1319,,,,,,,,,
+1320,,,,,,,,,
+1321,,,,,,,,,
+1322,,,,,,,,,
+1323,,,,,,,,,
+1324,,,,,,,,,
+1325,,,,,,,,,
+1326,,,,,,,,,
+1327,,,,,,,,,
+1328,,,,,,,,,
+1329,,,,,,,,,
+1330,,,,,,,,,
+1331,,,,,,,,,
+1332,,,,,,,,,
+1333,,,,,,,,,
+1334,,,,,,,,,
+1335,,,,,,,,,
+1336,,,,,,,,,
+1337,,,,,,,,,
+1338,,,,,,,,,
+1339,,,,,,,,,
+1340,,,,,,,-0.000003,,
+1341,,,,,,,0.000086,,
+1342,,,,,,,0.000194,,
+1343,,,,,,,0.000321,,
+1344,,,,,,,0.000469,,
+1345,,,,,,,0.000647,,
+1346,,,,,,,0.000856,,
+1347,,,,,,,0.001119,,
+1348,,,,,,,0.001429,,
+1349,,,,,,,0.001828,,
+1350,,,,,,,0.002318,,
+1351,,,,,,,0.00303,,
+1352,,,,,,,0.004011,,
+1353,,,,,,,0.005398,,
+1354,,,,,,,0.007194,,
+1355,,,,,,,0.011124,,
+1356,,,,,,,0.016985,,
+1357,,,,,,,0.026321,,
+1358,,,,,,,0.042756,,
+1359,,,,,,,0.06922,,
+1360,,,,,,,0.115351,,
+1361,,,,,,,0.182833,,
+1362,,,,,,,0.316333,,
+1363,,,,,,,0.45703,,
+1364,,,,,,,0.596852,,
+1365,,,,,,,0.772118,,
+1366,,,,,,,0.862394,,
+1367,,,,,,,0.871575,,
+1368,,,,,,,0.930308,,
+1369,,,,,,,0.962664,,
+1370,,,,,,,0.931247,,
+1371,,,,,,,0.938676,,
+1372,,,,,,,0.983748,,
+1373,,,,,,,0.999626,,
+1374,,,,,,,0.986968,,
+1375,,,,,,,1,,
+1376,,,,,,,0.983239,,
+1377,,,,,,,0.990577,,
+1378,,,,,,,0.963584,,
+1379,,,,,,,0.964054,,
+1380,,,,,,,0.871343,,
+1381,,,,,,,0.872689,,
+1382,,,,,,,0.703198,,
+1383,,,,,,,0.606578,,
+1384,,,,,,,0.452434,,
+1385,,,,,,,0.29792,,
+1386,,,,,,,0.190978,,
+1387,,,,,,,0.113055,,
+1388,,,,,,,0.066174,,
+1389,,,,,,,0.041513,,
+1390,,,,,,,0.026084,,
+1391,,,,,,,0.016218,,
+1392,,,,,,,0.009774,,
+1393,,,,,,,0.006839,,
+1394,,,,,,,0.004358,,
+1395,,,,,,,0.00278,,
+1396,,,,,,,0.001809,,
+1397,,,,,,,0.001351,,
+1398,,,,,,,0.000887,,
+1399,,,,,,,0.000481,,
+1400,,,,,,,0.00022,,
+1401,,,,,,,0.000114,,
+1402,,,,,,,0.000029,,
+1403,,,,,,,-0.00007,,
+1404,,,,,,,-0.000178,,
+1405,,,,,,,-0.000182,,
+1406,,,,,,,-0.000238,,
+1407,,,,,,,-0.000272,,
+1408,,,,,,,-0.000296,,
+1409,,,,,,,-0.000308,,
+1410,,,,,,,,,
+1411,,,,,,,,,
+1412,,,,,,,,,
+1413,,,,,,,,,
+1414,,,,,,,,,
+1415,,,,,,,,,
+1416,,,,,,,,,
+1417,,,,,,,,,
+1418,,,,,,,,,
+1419,,,,,,,,,
+1420,,,,,,,,,
+1421,,,,,,,,,
+1422,,,,,,,,,
+1423,,,,,,,,,
+1424,,,,,,,,,
+1425,,,,,,,,,
+1426,,,,,,,,,
+1427,,,,,,,,,
+1428,,,,,,,,,
+1429,,,,,,,,,
+1430,,,,,,,,,
+1431,,,,,,,,,
+1432,,,,,,,,,
+1433,,,,,,,,,
+1434,,,,,,,,,
+1435,,,,,,,,,
+1436,,,,,,,,,
+1437,,,,,,,,,
+1438,,,,,,,,,
+1439,,,,,,,,,
+1440,,,,,,,,,
+1441,,,,,,,,,
+1442,,,,,,,,,
+1443,,,,,,,,,
+1444,,,,,,,,,
+1445,,,,,,,,,
+1446,,,,,,,,,
+1447,,,,,,,,,
+1448,,,,,,,,,
+1449,,,,,,,,,
+1450,,,,,,,,,
+1451,,,,,,,,,
+1452,,,,,,,,,
+1453,,,,,,,,,
+1454,,,,,,,,,
+1455,,,,,,,,,
+1456,,,,,,,,,
+1457,,,,,,,,,
+1458,,,,,,,,,
+1459,,,,,,,,,
+1460,,,,,,,,,
+1461,,,,,,,,,
+1462,,,,,,,,,
+1463,,,,,,,,,
+1464,,,,,,,,,
+1465,,,,,,,,,
+1466,,,,,,,,,
+1467,,,,,,,,,
+1468,,,,,,,,,
+1469,,,,,,,,,
+1470,,,,,,,,,
+1471,,,,,,,,,
+1472,,,,,,,,,
+1473,,,,,,,,,
+1474,,,,,,,,,
+1475,,,,,,,,,
+1476,,,,,,,,,
+1477,,,,,,,,,
+1478,,,,,,,,,
+1479,,,,,,,,,
+1480,,,,,,,,,
+1481,,,,,,,,,
+1482,,,,,,,,,
+1483,,,,,,,,,
+1484,,,,,,,,,
+1485,,,,,,,,,
+1486,,,,,,,,,
+1487,,,,,,,,,
+1488,,,,,,,,,
+1489,,,,,,,,,
+1490,,,,,,,,,
+1491,,,,,,,,,
+1492,,,,,,,,,
+1493,,,,,,,,,
+1494,,,,,,,,,
+1495,,,,,,,,,
+1496,,,,,,,,,
+1497,,,,,,,,,
+1498,,,,,,,,,
+1499,,,,,,,,,
+1500,,,,,,,,,
+1501,,,,,,,,,
+1502,,,,,,,,,
+1503,,,,,,,,,
+1504,,,,,,,,,
+1505,,,,,,,,,
+1506,,,,,,,,,
+1507,,,,,,,,,
+1508,,,,,,,,,
+1509,,,,,,,,,
+1510,,,,,,,,,
+1511,,,,,,,,,
+1512,,,,,,,,,
+1513,,,,,,,,,
+1514,,,,,,,,,
+1515,,,,,,,,-0.000016,
+1516,,,,,,,,0.000067,
+1517,,,,,,,,0.000151,
+1518,,,,,,,,0.000249,
+1519,,,,,,,,0.000348,
+1520,,,,,,,,0.000466,
+1521,,,,,,,,0.000585,
+1522,,,,,,,,0.000758,
+1523,,,,,,,,0.000932,
+1524,,,,,,,,0.00115,
+1525,,,,,,,,0.001369,
+1526,,,,,,,,0.001613,
+1527,,,,,,,,0.001859,
+1528,,,,,,,,0.002172,
+1529,,,,,,,,0.002488,
+1530,,,,,,,,0.002881,
+1531,,,,,,,,0.003277,
+1532,,,,,,,,0.003772,
+1533,,,,,,,,0.004271,
+1534,,,,,,,,0.004898,
+1535,,,,,,,,0.005528,
+1536,,,,,,,,0.006421,
+1537,,,,,,,,0.007319,
+1538,,,,,,,,0.008459,
+1539,,,,,,,,0.009606,
+1540,,,,,,,,0.010989,
+1541,,,,,,,,0.012372,
+1542,,,,,,,,0.014303,
+1543,,,,,,,,0.016248,
+1544,,,,,,,,0.019029,
+1545,,,,,,,,0.021831,
+1546,,,,,,,,0.02589,
+1547,,,,,,,,0.029952,
+1548,,,,,,,,0.035171,
+1549,,,,,,,,0.040434,
+1550,,,,,,,,0.047864,
+1551,,,,,,,,0.055352,
+1552,,,,,,,,0.065732,
+1553,,,,,,,,0.076166,
+1554,,,,,,,,0.089024,
+1555,,,,,,,,0.101893,
+1556,,,,,,,,0.12015,
+1557,,,,,,,,0.138642,
+1558,,,,,,,,0.163127,
+1559,,,,,,,,0.187803,
+1560,,,,,,,,0.220261,
+1561,,,,,,,,0.252894,
+1562,,,,,,,,0.291359,
+1563,,,,,,,,0.329939,
+1564,,,,,,,,0.375648,
+1565,,,,,,,,0.42147,
+1566,,,,,,,,0.47356,
+1567,,,,,,,,0.525681,
+1568,,,,,,,,0.578787,
+1569,,,,,,,,0.631645,
+1570,,,,,,,,0.676683,
+1571,,,,,,,,0.721282,
+1572,,,,,,,,0.75477,
+1573,,,,,,,,0.788248,
+1574,,,,,,,,0.821155,
+1575,,,,,,,,0.854065,
+1576,,,,,,,,0.87337,
+1577,,,,,,,,0.89183,
+1578,,,,,,,,0.899817,
+1579,,,,,,,,0.907252,
+1580,,,,,,,,0.913009,
+1581,,,,,,,,0.918685,
+1582,,,,,,,,0.922953,
+1583,,,,,,,,0.927163,
+1584,,,,,,,,0.92686,
+1585,,,,,,,,0.926413,
+1586,,,,,,,,0.925059,
+1587,,,,,,,,0.923683,
+1588,,,,,,,,0.923953,
+1589,,,,,,,,0.924259,
+1590,,,,,,,,0.922828,
+1591,,,,,,,,0.921383,
+1592,,,,,,,,0.922061,
+1593,,,,,,,,0.922756,
+1594,,,,,,,,0.924648,
+1595,,,,,,,,0.926605,
+1596,,,,,,,,0.934552,
+1597,,,,,,,,0.942525,
+1598,,,,,,,,0.944351,
+1599,,,,,,,,0.945872,
+1600,,,,,,,,0.946175,
+1601,,,,,,,,0.946432,
+1602,,,,,,,,0.947006,
+1603,,,,,,,,0.947589,
+1604,,,,,,,,0.950194,
+1605,,,,,,,,0.952859,
+1606,,,,,,,,0.951303,
+1607,,,,,,,,0.94967,
+1608,,,,,,,,0.953047,
+1609,,,,,,,,0.956494,
+1610,,,,,,,,0.959047,
+1611,,,,,,,,0.961545,
+1612,,,,,,,,0.960048,
+1613,,,,,,,,0.958335,
+1614,,,,,,,,0.959835,
+1615,,,,,,,,0.96147,
+1616,,,,,,,,0.960857,
+1617,,,,,,,,0.960175,
+1618,,,,,,,,0.960813,
+1619,,,,,,,,0.961486,
+1620,,,,,,,,0.964703,
+1621,,,,,,,,0.967943,
+1622,,,,,,,,0.969314,
+1623,,,,,,,,0.970589,
+1624,,,,,,,,0.973713,
+1625,,,,,,,,0.976906,
+1626,,,,,,,,0.9791,
+1627,,,,,,,,0.981258,
+1628,,,,,,,,0.981285,
+1629,,,,,,,,0.981372,
+1630,,,,,,,,0.988609,
+1631,,,,,,,,0.995869,
+1632,,,,,,,,0.998021,
+1633,,,,,,,,1,
+1634,,,,,,,,0.999848,
+1635,,,,,,,,0.999642,
+1636,,,,,,,,0.99659,
+1637,,,,,,,,0.993439,
+1638,,,,,,,,0.986217,
+1639,,,,,,,,0.978989,
+1640,,,,,,,,0.967125,
+1641,,,,,,,,0.95512,
+1642,,,,,,,,0.936199,
+1643,,,,,,,,0.917239,
+1644,,,,,,,,0.879424,
+1645,,,,,,,,0.840967,
+1646,,,,,,,,0.796545,
+1647,,,,,,,,0.751893,
+1648,,,,,,,,0.694313,
+1649,,,,,,,,0.636542,
+1650,,,,,,,,0.573232,
+1651,,,,,,,,0.509946,
+1652,,,,,,,,0.451966,
+1653,,,,,,,,0.39403,
+1654,,,,,,,,0.34275,
+1655,,,,,,,,0.291752,
+1656,,,,,,,,0.251309,
+1657,,,,,,,,0.211153,
+1658,,,,,,,,0.180823,
+1659,,,,,,,,0.15054,
+1660,,,,,,,,0.128463,
+1661,,,,,,,,0.106664,
+1662,,,,,,,,0.090735,
+1663,,,,,,,,0.074941,
+1664,,,,,,,,0.06379,
+1665,,,,,,,,0.052752,
+1666,,,,,,,,0.045028,
+1667,,,,,,,,0.03731,
+1668,,,,,,,,0.031821,
+1669,,,,,,,,0.02635,
+1670,,,,,,,,0.022504,
+1671,,,,,,,,0.018724,
+1672,,,,,,,,0.016045,
+1673,,,,,,,,0.013394,
+1674,,,,,,,,0.011483,
+1675,,,,,,,,0.009587,
+1676,,,,,,,,0.008227,
+1677,,,,,,,,0.006882,
+1678,,,,,,,,0.00591,
+1679,,,,,,,,0.00494,
+1680,,,,,,,,0.004257,
+1681,,,,,,,,0.003576,
+1682,,,,,,,,0.003055,
+1683,,,,,,,,0.002541,
+1684,,,,,,,,0.00216,
+1685,,,,,,,,0.001781,
+1686,,,,,,,,0.001512,
+1687,,,,,,,,0.001244,
+1688,,,,,,,,0.00104,
+1689,,,,,,,,0.000837,
+1690,,,,,,,,0.000677,
+1691,,,,,,,,0.000517,
+1692,,,,,,,,0.000409,
+1693,,,,,,,,0.000301,
+1694,,,,,,,,0.000206,
+1695,,,,,,,,0.000112,
+1696,,,,,,,,0.00004,
+1697,,,,,,,,-0.000031,
+1698,,,,,,,,,
+1699,,,,,,,,,
+1700,,,,,,,,,
+1701,,,,,,,,,
+1702,,,,,,,,,
+1703,,,,,,,,,
+1704,,,,,,,,,
+1705,,,,,,,,,
+1706,,,,,,,,,
+1707,,,,,,,,,
+1708,,,,,,,,,
+1709,,,,,,,,,
+1710,,,,,,,,,
+1711,,,,,,,,,
+1712,,,,,,,,,
+1713,,,,,,,,,
+1714,,,,,,,,,
+1715,,,,,,,,,
+1716,,,,,,,,,
+1717,,,,,,,,,
+1718,,,,,,,,,
+1719,,,,,,,,,
+1720,,,,,,,,,
+1721,,,,,,,,,
+1722,,,,,,,,,
+1723,,,,,,,,,
+1724,,,,,,,,,
+1725,,,,,,,,,
+1726,,,,,,,,,
+1727,,,,,,,,,
+1728,,,,,,,,,
+1729,,,,,,,,,
+1730,,,,,,,,,
+1731,,,,,,,,,
+1732,,,,,,,,,
+1733,,,,,,,,,
+1734,,,,,,,,,
+1735,,,,,,,,,
+1736,,,,,,,,,
+1737,,,,,,,,,
+1738,,,,,,,,,
+1739,,,,,,,,,
+1740,,,,,,,,,
+1741,,,,,,,,,
+1742,,,,,,,,,
+1743,,,,,,,,,
+1744,,,,,,,,,
+1745,,,,,,,,,
+1746,,,,,,,,,
+1747,,,,,,,,,
+1748,,,,,,,,,
+1749,,,,,,,,,
+1750,,,,,,,,,
+1751,,,,,,,,,
+1752,,,,,,,,,
+1753,,,,,,,,,
+1754,,,,,,,,,
+1755,,,,,,,,,
+1756,,,,,,,,,
+1757,,,,,,,,,
+1758,,,,,,,,,
+1759,,,,,,,,,
+1760,,,,,,,,,
+1761,,,,,,,,,
+1762,,,,,,,,,
+1763,,,,,,,,,
+1764,,,,,,,,,
+1765,,,,,,,,,
+1766,,,,,,,,,
+1767,,,,,,,,,
+1768,,,,,,,,,
+1769,,,,,,,,,
+1770,,,,,,,,,
+1771,,,,,,,,,
+1772,,,,,,,,,
+1773,,,,,,,,,
+1774,,,,,,,,,
+1775,,,,,,,,,
+1776,,,,,,,,,
+1777,,,,,,,,,
+1778,,,,,,,,,
+1779,,,,,,,,,
+1780,,,,,,,,,
+1781,,,,,,,,,
+1782,,,,,,,,,
+1783,,,,,,,,,
+1784,,,,,,,,,
+1785,,,,,,,,,
+1786,,,,,,,,,
+1787,,,,,,,,,
+1788,,,,,,,,,
+1789,,,,,,,,,
+1790,,,,,,,,,
+1791,,,,,,,,,
+1792,,,,,,,,,
+1793,,,,,,,,,
+1794,,,,,,,,,
+1795,,,,,,,,,
+1796,,,,,,,,,
+1797,,,,,,,,,
+1798,,,,,,,,,
+1799,,,,,,,,,
+1800,,,,,,,,,
+1801,,,,,,,,,
+1802,,,,,,,,,
+1803,,,,,,,,,
+1804,,,,,,,,,
+1805,,,,,,,,,
+1806,,,,,,,,,
+1807,,,,,,,,,
+1808,,,,,,,,,
+1809,,,,,,,,,
+1810,,,,,,,,,
+1811,,,,,,,,,
+1812,,,,,,,,,
+1813,,,,,,,,,
+1814,,,,,,,,,
+1815,,,,,,,,,
+1816,,,,,,,,,
+1817,,,,,,,,,
+1818,,,,,,,,,
+1819,,,,,,,,,
+1820,,,,,,,,,
+1821,,,,,,,,,
+1822,,,,,,,,,
+1823,,,,,,,,,
+1824,,,,,,,,,
+1825,,,,,,,,,
+1826,,,,,,,,,
+1827,,,,,,,,,
+1828,,,,,,,,,
+1829,,,,,,,,,
+1830,,,,,,,,,
+1831,,,,,,,,,
+1832,,,,,,,,,
+1833,,,,,,,,,
+1834,,,,,,,,,
+1835,,,,,,,,,
+1836,,,,,,,,,
+1837,,,,,,,,,
+1838,,,,,,,,,
+1839,,,,,,,,,
+1840,,,,,,,,,
+1841,,,,,,,,,
+1842,,,,,,,,,
+1843,,,,,,,,,
+1844,,,,,,,,,
+1845,,,,,,,,,
+1846,,,,,,,,,
+1847,,,,,,,,,
+1848,,,,,,,,,
+1849,,,,,,,,,
+1850,,,,,,,,,
+1851,,,,,,,,,
+1852,,,,,,,,,
+1853,,,,,,,,,
+1854,,,,,,,,,
+1855,,,,,,,,,
+1856,,,,,,,,,
+1857,,,,,,,,,
+1858,,,,,,,,,
+1859,,,,,,,,,
+1860,,,,,,,,,
+1861,,,,,,,,,
+1862,,,,,,,,,
+1863,,,,,,,,,
+1864,,,,,,,,,
+1865,,,,,,,,,
+1866,,,,,,,,,
+1867,,,,,,,,,
+1868,,,,,,,,,
+1869,,,,,,,,,
+1870,,,,,,,,,
+1871,,,,,,,,,
+1872,,,,,,,,,
+1873,,,,,,,,,
+1874,,,,,,,,,
+1875,,,,,,,,,
+1876,,,,,,,,,
+1877,,,,,,,,,
+1878,,,,,,,,,
+1879,,,,,,,,,
+1880,,,,,,,,,
+1881,,,,,,,,,
+1882,,,,,,,,,
+1883,,,,,,,,,
+1884,,,,,,,,,
+1885,,,,,,,,,
+1886,,,,,,,,,
+1887,,,,,,,,,
+1888,,,,,,,,,
+1889,,,,,,,,,
+1890,,,,,,,,,
+1891,,,,,,,,,
+1892,,,,,,,,,
+1893,,,,,,,,,
+1894,,,,,,,,,
+1895,,,,,,,,,
+1896,,,,,,,,,
+1897,,,,,,,,,
+1898,,,,,,,,,
+1899,,,,,,,,,
+1900,,,,,,,,,
+1901,,,,,,,,,
+1902,,,,,,,,,
+1903,,,,,,,,,
+1904,,,,,,,,,
+1905,,,,,,,,,
+1906,,,,,,,,,
+1907,,,,,,,,,
+1908,,,,,,,,,
+1909,,,,,,,,,
+1910,,,,,,,,,
+1911,,,,,,,,,
+1912,,,,,,,,,
+1913,,,,,,,,,
+1914,,,,,,,,,
+1915,,,,,,,,,
+1916,,,,,,,,,
+1917,,,,,,,,,
+1918,,,,,,,,,
+1919,,,,,,,,,
+1920,,,,,,,,,
+1921,,,,,,,,,
+1922,,,,,,,,,
+1923,,,,,,,,,
+1924,,,,,,,,,
+1925,,,,,,,,,
+1926,,,,,,,,,
+1927,,,,,,,,,
+1928,,,,,,,,,
+1929,,,,,,,,,
+1930,,,,,,,,,
+1931,,,,,,,,,
+1932,,,,,,,,,
+1933,,,,,,,,,
+1934,,,,,,,,,
+1935,,,,,,,,,
+1936,,,,,,,,,
+1937,,,,,,,,,
+1938,,,,,,,,,
+1939,,,,,,,,,
+1940,,,,,,,,,
+1941,,,,,,,,,
+1942,,,,,,,,,
+1943,,,,,,,,,
+1944,,,,,,,,,
+1945,,,,,,,,,
+1946,,,,,,,,,
+1947,,,,,,,,,
+1948,,,,,,,,,
+1949,,,,,,,,,
+1950,,,,,,,,,
+1951,,,,,,,,,
+1952,,,,,,,,,
+1953,,,,,,,,,
+1954,,,,,,,,,
+1955,,,,,,,,,
+1956,,,,,,,,,
+1957,,,,,,,,,
+1958,,,,,,,,,
+1959,,,,,,,,,
+1960,,,,,,,,,
+1961,,,,,,,,,
+1962,,,,,,,,,
+1963,,,,,,,,,
+1964,,,,,,,,,
+1965,,,,,,,,,
+1966,,,,,,,,,
+1967,,,,,,,,,
+1968,,,,,,,,,
+1969,,,,,,,,,
+1970,,,,,,,,,
+1971,,,,,,,,,
+1972,,,,,,,,,
+1973,,,,,,,,,
+1974,,,,,,,,,
+1975,,,,,,,,,
+1976,,,,,,,,,
+1977,,,,,,,,,
+1978,,,,,,,,,
+1979,,,,,,,,,
+1980,,,,,,,,,
+1981,,,,,,,,,
+1982,,,,,,,,,
+1983,,,,,,,,,
+1984,,,,,,,,,
+1985,,,,,,,,,
+1986,,,,,,,,,
+1987,,,,,,,,,
+1988,,,,,,,,,
+1989,,,,,,,,,
+1990,,,,,,,,,
+1991,,,,,,,,,
+1992,,,,,,,,,
+1993,,,,,,,,,
+1994,,,,,,,,,
+1995,,,,,,,,,
+1996,,,,,,,,,
+1997,,,,,,,,,
+1998,,,,,,,,,
+1999,,,,,,,,,
+2000,,,,,,,,,
+2001,,,,,,,,,
+2002,,,,,,,,,
+2003,,,,,,,,,
+2004,,,,,,,,,
+2005,,,,,,,,,
+2006,,,,,,,,,
+2007,,,,,,,,,
+2008,,,,,,,,,
+2009,,,,,,,,,
+2010,,,,,,,,,
+2011,,,,,,,,,
+2012,,,,,,,,,
+2013,,,,,,,,,
+2014,,,,,,,,,
+2015,,,,,,,,,
+2016,,,,,,,,,
+2017,,,,,,,,,
+2018,,,,,,,,,
+2019,,,,,,,,,
+2020,,,,,,,,,
+2021,,,,,,,,,
+2022,,,,,,,,,
+2023,,,,,,,,,
+2024,,,,,,,,,
+2025,,,,,,,,,
+2026,,,,,,,,,
+2027,,,,,,,,,
+2028,,,,,,,,,
+2029,,,,,,,,,
+2030,,,,,,,,,
+2031,,,,,,,,,
+2032,,,,,,,,,
+2033,,,,,,,,,
+2034,,,,,,,,,
+2035,,,,,,,,,
+2036,,,,,,,,,
+2037,,,,,,,,,-0.00001
+2038,,,,,,,,,0.000037
+2039,,,,,,,,,0.000083
+2040,,,,,,,,,0.000131
+2041,,,,,,,,,0.000179
+2042,,,,,,,,,0.00024
+2043,,,,,,,,,0.000305
+2044,,,,,,,,,0.000368
+2045,,,,,,,,,0.00043
+2046,,,,,,,,,0.000512
+2047,,,,,,,,,0.000599
+2048,,,,,,,,,0.000704
+2049,,,,,,,,,0.000814
+2050,,,,,,,,,0.000947
+2051,,,,,,,,,0.001085
+2052,,,,,,,,,0.001222
+2053,,,,,,,,,0.00136
+2054,,,,,,,,,0.001546
+2055,,,,,,,,,0.001745
+2056,,,,,,,,,0.001964
+2057,,,,,,,,,0.002187
+2058,,,,,,,,,0.002439
+2059,,,,,,,,,0.002696
+2060,,,,,,,,,0.00301
+2061,,,,,,,,,0.003339
+2062,,,,,,,,,0.003733
+2063,,,,,,,,,0.004141
+2064,,,,,,,,,0.004627
+2065,,,,,,,,,0.005137
+2066,,,,,,,,,0.005728
+2067,,,,,,,,,0.006337
+2068,,,,,,,,,0.007139
+2069,,,,,,,,,0.007996
+2070,,,,,,,,,0.008903
+2071,,,,,,,,,0.009824
+2072,,,,,,,,,0.011005
+2073,,,,,,,,,0.012261
+2074,,,,,,,,,0.01361
+2075,,,,,,,,,0.014987
+2076,,,,,,,,,0.016872
+2077,,,,,,,,,0.018899
+2078,,,,,,,,,0.021
+2079,,,,,,,,,0.023121
+2080,,,,,,,,,0.025897
+2081,,,,,,,,,0.028847
+2082,,,,,,,,,0.032071
+2083,,,,,,,,,0.035363
+2084,,,,,,,,,0.040206
+2085,,,,,,,,,0.045432
+2086,,,,,,,,,0.050903
+2087,,,,,,,,,0.056429
+2088,,,,,,,,,0.06327
+2089,,,,,,,,,0.070409
+2090,,,,,,,,,0.079382
+2091,,,,,,,,,0.088907
+2092,,,,,,,,,0.10064
+2093,,,,,,,,,0.11283
+2094,,,,,,,,,0.128292
+2095,,,,,,,,,0.144668
+2096,,,,,,,,,0.162055
+2097,,,,,,,,,0.179714
+2098,,,,,,,,,0.20278
+2099,,,,,,,,,0.227234
+2100,,,,,,,,,0.253732
+2101,,,,,,,,,0.280925
+2102,,,,,,,,,0.311347
+2103,,,,,,,,,0.342526
+2104,,,,,,,,,0.377044
+2105,,,,,,,,,0.412621
+2106,,,,,,,,,0.45047
+2107,,,,,,,,,0.488816
+2108,,,,,,,,,0.522458
+2109,,,,,,,,,0.554715
+2110,,,,,,,,,0.593227
+2111,,,,,,,,,0.633521
+2112,,,,,,,,,0.663067
+2113,,,,,,,,,0.689664
+2114,,,,,,,,,0.722284
+2115,,,,,,,,,0.756529
+2116,,,,,,,,,0.776463
+2117,,,,,,,,,0.792667
+2118,,,,,,,,,0.813716
+2119,,,,,,,,,0.836001
+2120,,,,,,,,,0.846344
+2121,,,,,,,,,0.853714
+2122,,,,,,,,,0.867845
+2123,,,,,,,,,0.883615
+2124,,,,,,,,,0.886411
+2125,,,,,,,,,0.886127
+2126,,,,,,,,,0.895232
+2127,,,,,,,,,0.906527
+2128,,,,,,,,,0.909739
+2129,,,,,,,,,0.911091
+2130,,,,,,,,,0.917985
+2131,,,,,,,,,0.926131
+2132,,,,,,,,,0.929693
+2133,,,,,,,,,0.932227
+2134,,,,,,,,,0.936544
+2135,,,,,,,,,0.941406
+2136,,,,,,,,,0.942571
+2137,,,,,,,,,0.942952
+2138,,,,,,,,,0.943112
+2139,,,,,,,,,0.943194
+2140,,,,,,,,,0.945168
+2141,,,,,,,,,0.947537
+2142,,,,,,,,,0.948776
+2143,,,,,,,,,0.949694
+2144,,,,,,,,,0.949643
+2145,,,,,,,,,0.9494
+2146,,,,,,,,,0.952551
+2147,,,,,,,,,0.956635
+2148,,,,,,,,,0.953083
+2149,,,,,,,,,0.947423
+2150,,,,,,,,,0.949094
+2151,,,,,,,,,0.952773
+2152,,,,,,,,,0.950874
+2153,,,,,,,,,0.947477
+2154,,,,,,,,,0.947014
+2155,,,,,,,,,0.947085
+2156,,,,,,,,,0.951812
+2157,,,,,,,,,0.957717
+2158,,,,,,,,,0.953332
+2159,,,,,,,,,0.946412
+2160,,,,,,,,,0.947778
+2161,,,,,,,,,0.951119
+2162,,,,,,,,,0.951641
+2163,,,,,,,,,0.951518
+2164,,,,,,,,,0.948644
+2165,,,,,,,,,0.944956
+2166,,,,,,,,,0.942515
+2167,,,,,,,,,0.940311
+2168,,,,,,,,,0.9434
+2169,,,,,,,,,0.947923
+2170,,,,,,,,,0.94501
+2171,,,,,,,,,0.940213
+2172,,,,,,,,,0.938737
+2173,,,,,,,,,0.93806
+2174,,,,,,,,,0.941859
+2175,,,,,,,,,0.947019
+2176,,,,,,,,,0.946554
+2177,,,,,,,,,0.944482
+2178,,,,,,,,,0.947414
+2179,,,,,,,,,0.951661
+2180,,,,,,,,,0.949283
+2181,,,,,,,,,0.945318
+2182,,,,,,,,,0.939939
+2183,,,,,,,,,0.934142
+2184,,,,,,,,,0.935493
+2185,,,,,,,,,0.93882
+2186,,,,,,,,,0.939253
+2187,,,,,,,,,0.938955
+2188,,,,,,,,,0.934675
+2189,,,,,,,,,0.929162
+2190,,,,,,,,,0.927085
+2191,,,,,,,,,0.925692
+2192,,,,,,,,,0.930508
+2193,,,,,,,,,0.936899
+2194,,,,,,,,,0.933908
+2195,,,,,,,,,0.927984
+2196,,,,,,,,,0.930981
+2197,,,,,,,,,0.936472
+2198,,,,,,,,,0.935776
+2199,,,,,,,,,0.933523
+2200,,,,,,,,,0.935132
+2201,,,,,,,,,0.937592
+2202,,,,,,,,,0.946217
+2203,,,,,,,,,0.956567
+2204,,,,,,,,,0.955661
+2205,,,,,,,,,0.951991
+2206,,,,,,,,,0.956666
+2207,,,,,,,,,0.963135
+2208,,,,,,,,,0.964442
+2209,,,,,,,,,0.964365
+2210,,,,,,,,,0.963523
+2211,,,,,,,,,0.962434
+2212,,,,,,,,,0.962905
+2213,,,,,,,,,0.963685
+2214,,,,,,,,,0.962473
+2215,,,,,,,,,0.960741
+2216,,,,,,,,,0.959239
+2217,,,,,,,,,0.957814
+2218,,,,,,,,,0.957781
+2219,,,,,,,,,0.958041
+2220,,,,,,,,,0.953274
+2221,,,,,,,,,0.94716
+2222,,,,,,,,,0.951706
+2223,,,,,,,,,0.958833
+2224,,,,,,,,,0.960212
+2225,,,,,,,,,0.960339
+2226,,,,,,,,,0.954785
+2227,,,,,,,,,0.947696
+2228,,,,,,,,,0.951965
+2229,,,,,,,,,0.95906
+2230,,,,,,,,,0.960554
+2231,,,,,,,,,0.960764
+2232,,,,,,,,,0.95575
+2233,,,,,,,,,0.949261
+2234,,,,,,,,,0.953245
+2235,,,,,,,,,0.95997
+2236,,,,,,,,,0.963736
+2237,,,,,,,,,0.966786
+2238,,,,,,,,,0.964315
+2239,,,,,,,,,0.960173
+2240,,,,,,,,,0.965473
+2241,,,,,,,,,0.973416
+2242,,,,,,,,,0.977637
+2243,,,,,,,,,0.980904
+2244,,,,,,,,,0.98276
+2245,,,,,,,,,0.984155
+2246,,,,,,,,,0.984693
+2247,,,,,,,,,0.985056
+2248,,,,,,,,,0.991486
+2249,,,,,,,,,0.9996
+2250,,,,,,,,,0.997654
+2251,,,,,,,,,0.993153
+2252,,,,,,,,,0.992469
+2253,,,,,,,,,0.992675
+2254,,,,,,,,,0.995894
+2255,,,,,,,,,1
+2256,,,,,,,,,0.999279
+2257,,,,,,,,,0.997261
+2258,,,,,,,,,0.994356
+2259,,,,,,,,,0.991127
+2260,,,,,,,,,0.987747
+2261,,,,,,,,,0.984349
+2262,,,,,,,,,0.986037
+2263,,,,,,,,,0.989132
+2264,,,,,,,,,0.984536
+2265,,,,,,,,,0.978024
+2266,,,,,,,,,0.975019
+2267,,,,,,,,,0.972794
+2268,,,,,,,,,0.974168
+2269,,,,,,,,,0.97654
+2270,,,,,,,,,0.976199
+2271,,,,,,,,,0.975206
+2272,,,,,,,,,0.974409
+2273,,,,,,,,,0.973662
+2274,,,,,,,,,0.967502
+2275,,,,,,,,,0.959895
+2276,,,,,,,,,0.956943
+2277,,,,,,,,,0.955095
+2278,,,,,,,,,0.955085
+2279,,,,,,,,,0.955588
+2280,,,,,,,,,0.947195
+2281,,,,,,,,,0.93665
+2282,,,,,,,,,0.922405
+2283,,,,,,,,,0.907109
+2284,,,,,,,,,0.89494
+2285,,,,,,,,,0.883588
+2286,,,,,,,,,0.855489
+2287,,,,,,,,,0.823876
+2288,,,,,,,,,0.78488
+2289,,,,,,,,,0.744025
+2290,,,,,,,,,0.69852
+2291,,,,,,,,,0.651673
+2292,,,,,,,,,0.602539
+2293,,,,,,,,,0.552647
+2294,,,,,,,,,0.502693
+2295,,,,,,,,,0.452698
+2296,,,,,,,,,0.403993
+2297,,,,,,,,,0.355569
+2298,,,,,,,,,0.315712
+2299,,,,,,,,,0.27826
+2300,,,,,,,,,0.244645
+2301,,,,,,,,,0.212213
+2302,,,,,,,,,0.186151
+2303,,,,,,,,,0.161749
+2304,,,,,,,,,0.141435
+2305,,,,,,,,,0.122015
+2306,,,,,,,,,0.106531
+2307,,,,,,,,,0.092029
+2308,,,,,,,,,0.080268
+2309,,,,,,,,,0.069276
+2310,,,,,,,,,0.060703
+2311,,,,,,,,,0.052702
+2312,,,,,,,,,0.046332
+2313,,,,,,,,,0.040405
+2314,,,,,,,,,0.035634
+2315,,,,,,,,,0.031213
+2316,,,,,,,,,0.027516
+2317,,,,,,,,,0.024
+2318,,,,,,,,,0.021262
+2319,,,,,,,,,0.018688
+2320,,,,,,,,,0.016562
+2321,,,,,,,,,0.014545
+2322,,,,,,,,,0.01293
+2323,,,,,,,,,0.011426
+2324,,,,,,,,,0.010155
+2325,,,,,,,,,0.008959
+2326,,,,,,,,,0.007992
+2327,,,,,,,,,0.007088
+2328,,,,,,,,,0.006348
+2329,,,,,,,,,0.005643
+2330,,,,,,,,,0.005019
+2331,,,,,,,,,0.004415
+2332,,,,,,,,,0.003903
+2333,,,,,,,,,0.003421
+2334,,,,,,,,,0.003025
+2335,,,,,,,,,0.002651
+2336,,,,,,,,,0.00234
+2337,,,,,,,,,0.002047
+2338,,,,,,,,,0.001795
+2339,,,,,,,,,0.001554
+2340,,,,,,,,,0.001345
+2341,,,,,,,,,0.001145
+2342,,,,,,,,,0.000974
+2343,,,,,,,,,0.000811
+2344,,,,,,,,,0.00068
+2345,,,,,,,,,0.00056
+2346,,,,,,,,,0.00044
+2347,,,,,,,,,0.00032
+2348,,,,,,,,,0.000217
+2349,,,,,,,,,0.000119
+2350,,,,,,,,,0.000028
+2351,,,,,,,,,-0.000062
+2352,,,,,,,,,-0.000134
+2353,,,,,,,,,-0.000202
+2354,,,,,,,,,-0.000263
+2355,,,,,,,,,-0.000322
More information about the grass-commit
mailing list