Added VGT1 and VGT2 atmospheric correction. Thanks to Alfredo Alessandrini for the extraction from 6S sources.

git-svn-id: https://svn.osgeo.org/grass/grass/trunk@46980 15284696-431f-4ddb-bdfa-cd5b030d7da7

imagery/i.atcorr/GeomCond.cpp

@@ -367,6 +367,8 @@ void GeomCond::parse()
     case 11: /* avnir            * enter month,day,hh.ddd,long.,lat. */
     case 12: /* ikonos           * enter month,day,hh.ddd,long.,lat. */
     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. */
     {
 	cin >> month;
 	cin >> jday;
@@ -412,7 +414,7 @@ void GeomCond::parse()
 /* ---- print geometrical conditions ---- */
 void GeomCond::print()
 {
-    static const string etiq1[14] = {
+    static const string etiq1[16] = {
 	string(" user defined conditions     "),
 	string(" meteosat observation        "),
 	string(" goes east observation       "),
@@ -426,7 +428,9 @@ void GeomCond::print()
 	string(" aster    observation        "),
 	string(" avnir    observation        "),
 	string(" ikonos   observation        "),
-	string(" rapideye observation        ")
+	string(" rapideye observation        "),
+	string(" vgt1_spot4 observation      "),
+	string(" vgt2_spot5 observation      ")
     };
 
     static const string head(" geometrical conditions identity  ");

imagery/i.atcorr/GeomCond.h

@@ -55,6 +55,10 @@
 /*                                                                      c */
 /*         13      rapideye          * enter month,day,hh.ddd,long.,lat.c */
 /*                                                                      c */
+/*         14      vgt1_spot4        * enter month,day,hh.ddd,long.,lat.c */
+/*                                                                      c */
+/*         15      vgt2_spot5        * enter month,day,hh.ddd,long.,lat.c */
+/*                                                                      c */
 /*                                                                      c */
 /*     note:       for hrv and tm experiments long. and lat. are the    c */
 /*                 coordinates of the scene center.                     c */

imagery/i.atcorr/Iwave.cpp

@@ -2266,6 +2266,197 @@ void IWave::rapideye(int iwa)
     }
 }
 
+/* Following filter function created using create_iwave.py */
+
+void IWave::vgt1_spot4(int iwa)
+{
+
+    /* B0 of VGT1-spot4 */
+    static const float sr1[41] = {
+		.0001, .0003, .0005, .0007, .0010, .0013, .0028, .0050,
+		.0182, .0377, .0828, .1423, .2234, .3168, .4108, .5052,
+		.5925, .6759, .7509, .8212, .8815, .9361, .9699, .9920,
+		.9907, .9761, .9528, .9248, .8811, .8287, .7223, .5857,
+		.4394, .2876, .1820, .1024, .0579, .0331, .0196, .0125,
+		.0000
+    };
+
+    /* B2 of VGT1-spot4 */
+    static const float sr2[82] = {
+		.0017, .0017, .0016, .0016, .0016, .0024, .0031, .0039,
+		.0046, .0439, .0832, .1225, .1618, .2737, .3856, .4975,
+		.6094, .6615, .7136, .7657, .8178, .8278, .8377, .8477,
+		.8577, .8673, .8769, .8865, .8961, .9088, .9214, .9341,
+		.9468, .9595, .9722, .9849, .9976, .9895, .9815, .9734,
+		.9653, .9169, .8685, .8201, .7717, .6989, .6260, .5532,
+		.4804, .4243, .3681, .3120, .2559, .2229, .1899, .1569,
+		.1239, .1063, .0887, .0712, .0536, .0454, .0372, .0289,
+		.0207, .0175, .0143, .0111, .0079, .0069, .0060, .0050,
+		.0040, .0034, .0027, .0021, .0015, .0015, .0015, .0015,
+		.0015, .0000
+    };
+
+    /* B3 of VGT1-spot4 */
+    static const float sr3[133] = {
+		.0007, .0008, .0008, .0009, .0010, .0011, .0013, .0014,
+		.0015, .0019, .0022, .0026, .0030, .0051, .0072, .0092,
+		.0113, .0210, .0308, .0406, .0503, .0806, .1109, .1411,
+		.1714, .2271, .2828, .3385, .3942, .4584, .5226, .5869,
+		.6511, .6966, .7421, .7876, .8331, .8607, .8883, .9159,
+		.9435, .9552, .9668, .9785, .9902, .9927, .9951, .9976,
+		1.0000, .9999, .9997, .9996, .9994, .9984, .9974, .9964,
+		.9954, .9927, .9900, .9873, .9846, .9766, .9686, .9606,
+		.9526, .9445, .9365, .9284, .9203, .9008, .8813, .8618,
+		.8423, .7897, .7370, .6844, .6318, .5584, .4849, .4114,
+		.3380, .2893, .2407, .1920, .1434, .1222, .1010, .0798,
+		.0586, .0505, .0425, .0345, .0264, .0229, .0194, .0159,
+		.0124, .0111, .0098, .0085, .0072, .0065, .0057, .0050,
+		.0042, .0038, .0034, .0031, .0027, .0024, .0021, .0017,
+		.0014, .0014, .0014, .0013, .0013, .0013, .0013, .0012,
+		.0012, .0012, .0011, .0011, .0011, .0010, .0009, .0009,
+		.0008, .0007, .0006, .0006, .0005
+    };
+
+    /* MIR of VGT1-spot4 */
+    static const float sr4[141] = {
+		.0010, .0010, .0010, .0010, .0010, .0011, .0011, .0012,
+		.0013, .0014, .0014, .0015, .0016, .0016, .0017, .0018,
+		.0018, .0019, .0021, .0022, .0024, .0026, .0028, .0030,
+		.0032, .0035, .0038, .0041, .0044, .0050, .0055, .0060,
+		.0066, .0075, .0083, .0092, .0101, .0118, .0136, .0154,
+		.0171, .0203, .0235, .0267, .0299, .0366, .0433, .0501,
+		.0568, .0708, .0848, .0989, .1129, .1406, .1684, .1961,
+		.2239, .2691, .3142, .3594, .4045, .4677, .5308, .5939,
+		.6571, .7099, .7628, .8156, .8685, .8934, .9183, .9433,
+		.9682, .9761, .9841, .9920, 1.0000, .9971, .9941, .9912,
+		.9883, .9840, .9797, .9753, .9710, .9517, .9325, .9132,
+		.8939, .8713, .8488, .8262, .8036, .7457, .6878, .6299,
+		.5720, .5109, .4498, .3887, .3276, .2926, .2576, .2226,
+		.1876, .1704, .1531, .1358, .1186, .1078, .0970, .0863,
+		.0755, .0685, .0614, .0544, .0474, .0430, .0385, .0340,
+		.0296, .0269, .0242, .0214, .0187, .0172, .0157, .0142,
+		.0127, .0118, .0109, .0101, .0092, .0086, .0080, .0075,
+		.0069, .0064, .0060, .0055, .0050
+    };
+
+    static const float wli[4] = {0.400, 0.580, 0.700, 1.450};
+    static const float wls[4] = {0.500, 0.782, 1.030, 1.800};
+
+    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 < 41; i++)  ffu.s[60+i] = sr1[i];
+        break;
+    case 2: for(i = 0; i < 82; i++)  ffu.s[132+i] = sr2[i];
+        break;
+    case 3: for(i = 0; i < 133; i++)  ffu.s[180+i] = sr3[i];
+        break;
+    case 4: for(i = 0; i < 141; i++)  ffu.s[480+i] = sr4[i];
+        break;
+    }
+}
+
+/* Following filter function created using create_iwave.py */
+
+void IWave::vgt2_spot5(int iwa)
+{
+
+    /* B0 of VGT2_spot5 */
+    static const float sr1[61] = {
+		.0007, .0008, .0010, .0010, .0011, .0012, .0013, .0029,
+		.0045, .0209, .0374, .0964, .1553, .2551, .3549, .4536,
+		.5524, .6318, .7111, .7602, .8093, .8304, .8514, .8698,
+		.8882, .9441, 1.0000, .9787, .9574, .7493, .5412, .3553,
+		.1694, .1083, .0472, .0318, .0163, .0116, .0069, .0053,
+		.0036, .0028, .0020, .0017, .0014, .0011, .0009, .0008,
+		.0008, .0008, .0007, .0006, .0006, .0006, .0005, .0004,
+		.0004, .0004, .0004, .0003, .0003
+    };
+
+    /* B2 of VGT2_spot5 */
+    static const float sr2[81] = {
+		.0006, .0006, .0006, .0007, .0007, .0033, .0059, .0084,
+		.0110, .0688, .1265, .1842, .2420, .3532, .4643, .5755,
+		.6866, .7306, .7745, .8185, .8625, .8766, .8906, .9046,
+		.9187, .9307, .9426, .9546, .9666, .9749, .9833, .9917,
+		1.0000, .9957, .9915, .9872, .9829, .9387, .8945, .8503,
+		.8061, .7316, .6570, .5825, .5079, .4468, .3856, .3245,
+		.2634, .2293, .1952, .1611, .1270, .1095, .0921, .0746,
+		.0572, .0490, .0408, .0326, .0244, .0208, .0171, .0135,
+		.0099, .0085, .0071, .0057, .0043, .0037, .0031, .0025,
+		.0019, .0017, .0014, .0012, .0010, .0009, .0008, .0007,
+		.0006
+    };
+
+    /* B3 of VGT2_spot5 */
+    static const float sr3[121] = {
+		.0003, .0003, .0003, .0003, .0003, .0003, .0003, .0004,
+		.0004, .0005, .0006, .0008, .0009, .0017, .0024, .0032,
+		.0040, .0082, .0124, .0166, .0208, .0361, .0514, .0667,
+		.0820, .1175, .1530, .1885, .2240, .2774, .3309, .3844,
+		.4378, .4906, .5434, .5962, .6490, .6899, .7308, .7716,
+		.8125, .8369, .8612, .8856, .9100, .9248, .9396, .9544,
+		.9692, .9761, .9829, .9898, .9967, .9975, .9983, .9992,
+		1.0000, .9939, .9878, .9817, .9756, .9618, .9480, .9342,
+		.9204, .9030, .8856, .8682, .8508, .8236, .7965, .7693,
+		.7422, .6945, .6467, .5989, .5512, .4910, .4308, .3706,
+		.3104, .2681, .2258, .1836, .1413, .1213, .1013, .0813,
+		.0613, .0530, .0447, .0364, .0281, .0245, .0209, .0172,
+		.0136, .0120, .0104, .0089, .0073, .0065, .0056, .0047,
+		.0039, .0036, .0032, .0029, .0025, .0023, .0021, .0018,
+		.0016, .0015, .0014, .0012, .0011, .0010, .0008, .0007,
+		.0006
+    };
+
+    /* MIR of VGT2_spot5 */
+    static const float sr4[141] = {
+		.0005, .0006, .0006, .0007, .0008, .0008, .0008, .0009,
+		.0009, .0009, .0009, .0009, .0009, .0011, .0014, .0016,
+		.0019, .0022, .0024, .0027, .0029, .0034, .0040, .0045,
+		.0050, .0060, .0070, .0080, .0090, .0108, .0126, .0143,
+		.0161, .0196, .0231, .0266, .0301, .0372, .0444, .0515,
+		.0586, .0723, .0860, .0997, .1134, .1401, .1668, .1936,
+		.2203, .2668, .3134, .3599, .4065, .4664, .5263, .5863,
+		.6462, .6929, .7396, .7862, .8329, .8630, .8931, .9233,
+		.9534, .9644, .9754, .9863, .9973, .9980, .9987, .9993,
+		1.0000, .9987, .9974, .9962, .9949, .9914, .9879, .9845,
+		.9810, .9687, .9564, .9440, .9317, .9055, .8793, .8532,
+		.8270, .7811, .7351, .6892, .6433, .5746, .5059, .4371,
+		.3684, .3244, .2803, .2363, .1923, .1713, .1502, .1292,
+		.1082, .0970, .0859, .0747, .0636, .0575, .0514, .0453,
+		.0392, .0360, .0327, .0295, .0262, .0235, .0209, .0182,
+		.0155, .0139, .0123, .0106, .0090, .0082, .0073, .0065,
+		.0057, .0053, .0050, .0046, .0043, .0040, .0037, .0033,
+		.0030, .0026, .0022, .0019, .0015
+    };
+
+    static const float wli[4] = {0.400, 0.580, 0.700, 1.450};
+    static const float wls[4] = {0.550, 0.780, 1.000, 1.800};
+
+    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 < 61; i++)  ffu.s[60+i] = sr1[i];
+        break;
+    case 2: for(i = 0; i < 81; i++)  ffu.s[132+i] = sr2[i];
+        break;
+    case 3: for(i = 0; i < 121; i++)  ffu.s[180+i] = sr3[i];
+        break;
+    case 4: for(i = 0; i < 141; i++)  ffu.s[480+i] = sr4[i];
+        break;
+    }
+}
+
 float IWave::equivwl() const
 {
     float seb = 0;
@@ -2336,6 +2527,8 @@ void IWave::parse()
 	else if(iwave <= 84)    avnir(iwave - 80);
 	else if(iwave <= 87)    ikonos(iwave - 84);
 	else if(iwave <= 92)    rapideye(iwave - 87);
+	else if(iwave <= 96)    vgt1_spot4(iwave - 92);
+	else if(iwave <= 100)   vgt2_spot5(iwave - 96);
 	else G_warning(_("Unsupported iwave value: %d"), iwave);
     }
 
@@ -2352,7 +2545,7 @@ void IWave::parse()
 /* --- spectral condition ---- */
 void IWave::print()
 {
-    static const string nsat[94] = {
+    static const string nsat[102] = {
 	string(" constant        "), string(" user s          "),
 	string(" meteosat        "), string(" goes east       "), string(" goes west       "),
 	string(" avhrr 1 (noaa6) "), string(" avhrr 2 (noaa6) "),
@@ -2392,6 +2585,10 @@ void IWave::print()
 	string(" ikonos green    "), string(" ikonos red      "),	string(" ikonos NIR      "),
 	string(" rapideye blue   "), string(" rapideye green  "),	string(" rapideye red    "),
 	string(" rapideye rededge"), string(" rapideye NIR    "),
+	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  "),
     };
 
 

imagery/i.atcorr/Iwave.h

@@ -109,11 +109,19 @@ c        84  4th      "               ( 0.700-0.925 )                  c
 c        85  Green band of IKONOS     ( 0.350-1.035 )                  c
 c        86  Red      "               ( 0.350-1.035 )                  c
 c        87  NIR      "               ( 0.350-1.035 )                  c
-c        88  Blue band of RapidEye  ( 0.438-0.513 )                  c
-c        89  Green      "               ( 0.463-0.594 )                  c
+c        88  Blue band of RapidEye    ( 0.438-0.513 )                  c
+c        89  Green    "               ( 0.463-0.594 )                  c
 c        90  Red      "               ( 0.624-0.690 )                  c
-c        91  RedEdge      "               ( 0.500-0.737 )                  c
+c        91  RedEdge  "               ( 0.500-0.737 )                  c
 c        92  NIR      "               ( 0.520-0.862 )                  c
+c        93  B0 band of VGT1 (SPOT4)  ( 0.400-0.500 )                  c
+c        94  B2       "               ( 0.580-0.782 )                  c
+c        95  B3       "               ( 0.700-1.030 )                  c
+c        96  MIR      "               ( 1.450-1.800 )                  c
+c        97  B0 band of VGT2 (SPOT5)  ( 0.400-0.550 )                  c
+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  note: wl has to be in micrometer                                    c
 c**********************************************************************/
 
@@ -153,6 +161,8 @@ private:
 	void avnir(int iwa);
 	void ikonos(int iwa);
 	void rapideye(int iwa);
+	void vgt1_spot4(int iwa);
+	void vgt2_spot5(int iwa);
 
 
 public:

imagery/i.atcorr/README

@@ -12,8 +12,14 @@ Adding new band filters (see also Iwave.h):
    - then the ffu.s array is filled with each of the spectral band data
      one after the other.
      Example Lsat TM: (435-250) / 2.5 = 74
- o add iwave values, also in Iwave.h
+ 
+ o add filter function to Iwave.cpp, before IWave::equivwl()
+ o add else-if in IWave::parse() in Iwave.cpp
+ o add print strings in Iwave.cpp
+ 
+ o add iwave values and signature in Iwave.h
  o add to GeomCond.cpp and GeomCond.h
+ o add to i.atcorr.html
 
 
 Preparation of the Iwave.cpp file code:

imagery/i.atcorr/i.atcorr.html

@@ -155,6 +155,17 @@ the ascendant node at equator</td>
 <td>enter month,day,hh.ddd,long.,lat. *</td>
 </tr>
 
+<tr>
+<td>14</td>
+<td><b>VGT1 (SPOT4)</b></td>
+<td>enter month,day,hh.ddd,long.,lat. *</td>
+</tr>
+
+<tr>
+<td>15</td>
+<td><b>VGT2 (SPOT5)</b></td>
+<td>enter month,day,hh.ddd,long.,lat. *</td>
+</tr>
 </table>
 
 <blockquote>
@@ -546,6 +557,15 @@ Pre-defined satellite bands:
 <tr><td>91</td><td>RapidEye RedEdge band (0.500-0.737)</td></tr>
 <tr><td>92</td><td>RapidEye NIR band (0.520-0.862)</td></tr>
 
+<tr><td>93</td><td><b>VGT1 (SPOT4)</b> band 0 (0.400-0.500)</td></tr>
+<tr><td>94</td><td>VGT1 (SPOT4) band 2 (0.580-0.782)</td></tr>
+<tr><td>95</td><td>VGT1 (SPOT4) band 3 (0.700-1.030)</td></tr>
+<tr><td>96</td><td>VGT1 (SPOT4) MIR band (1.450-1.800)</td></tr>
+
+<tr><td>97</td><td><b>VGT2 (SPOT5)</b> band 0 (0.400-0.550)</td></tr>
+<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>
 </table>
 
 <h2>EXAMPLES</h2>
@@ -715,6 +735,9 @@ overview., IEEE Trans. Geosc. and Remote Sens. 35(3):675-686.
 <p><em>RapidEye sensors addition 11/2010:</em>
 <br>Peter L&ouml;we, Anne Ghisla
 
+<p><em>VGT1 and VGT2 sensors addition from 6S sources 07/2011:</em>
+<br>Alfredo Alessandrini, Anne Ghisla
+
 <p>
 <i>Last changed: $Date$</i>
 

imagery/i.atcorr/sensors_csv/VGT1_spot4.csv

@@ -0,0 +1,562 @@
+WL(nm),B0,B2,B3,MIR
+400,0.000128,,,
+402.5,0.000328,,,
+405,0.000528,,,
+407.5,0.000728,,,
+410,0.000992,,,
+412.5,0.001292,,,
+415,0.002808,,,
+417.5,0.005008,,,
+420,0.018248,,,
+422.5,0.037698,,,
+425,0.08278,,,
+427.5,0.14228,,,
+430,0.223444,,,
+432.5,0.316794,,,
+435,0.410816,,,
+437.5,0.505216,,,
+440,0.592544,,,
+442.5,0.675894,,,
+445,0.750892,,,
+447.5,0.821192,,,
+450,0.881476,,,
+452.5,0.936126,,,
+455,0.969944,,,
+457.5,0.992044,,,
+460,0.990656,,,
+462.5,0.976056,,,
+465,0.952848,,,
+467.5,0.924798,,,
+470,0.881132,,,
+472.5,0.828682,,,
+475,0.722344,,,
+477.5,0.585694,,,
+480,0.43938,,,
+482.5,0.28763,,,
+485,0.182024,,,
+487.5,0.102374,,,
+490,0.05786,,,
+492.5,0.03311,,,
+495,0.019624,,,
+497.5,0.012474,,,
+500,8.8633101E-017,,,
+502.5,,,,
+505,,,,
+507.5,,,,
+510,,,,
+512.5,,,,
+515,,,,
+517.5,,,,
+520,,,,
+522.5,,,,
+525,,,,
+527.5,,,,
+530,,,,
+532.5,,,,
+535,,,,
+537.5,,,,
+540,,,,
+542.5,,,,
+545,,,,
+547.5,,,,
+550,,,,
+552.5,,,,
+555,,,,
+557.5,,,,
+560,,,,
+562.5,,,,
+565,,,,
+567.5,,,,
+570,,,,
+572.5,,,,
+575,,,,
+577.5,,,,
+580,,0.0017,,
+582.5,,0.001675,,
+585,,0.00165,,
+587.5,,0.001625,,
+590,,0.0016,,
+592.5,,0.00235,,
+595,,0.0031,,
+597.5,,0.00385,,
+600,,0.0046,,
+602.5,,0.0439,,
+605,,0.0832,,
+607.5,,0.1225,,
+610,,0.1618,,
+612.5,,0.2737,,
+615,,0.3856,,
+617.5,,0.4975,,
+620,,0.6094,,
+622.5,,0.6615,,
+625,,0.7136,,
+627.5,,0.7657,,
+630,,0.8178,,
+632.5,,0.827775,,
+635,,0.83775,,
+637.5,,0.847725,,
+640,,0.8577,,
+642.5,,0.8673,,
+645,,0.8769,,
+647.5,,0.8865,,
+650,,0.8961,,
+652.5,,0.908775,,
+655,,0.92145,,
+657.5,,0.934125,,
+660,,0.9468,,
+662.5,,0.9595,,
+665,,0.9722,,
+667.5,,0.9849,,
+670,,0.9976,,
+672.5,,0.989525,,
+675,,0.98145,,
+677.5,,0.973375,,
+680,,0.9653,,
+682.5,,0.9169,,
+685,,0.8685,,
+687.5,,0.8201,,
+690,,0.7717,,
+692.5,,0.698875,,
+695,,0.62605,,
+697.5,,0.553225,,
+700,,0.4804,0.0007,
+702.5,,0.424275,0.000775,
+705,,0.36815,0.00085,
+707.5,,0.312025,0.000925,
+710,,0.2559,0.001,
+712.5,,0.2229,0.001125,
+715,,0.1899,0.00125,
+717.5,,0.1569,0.001375,
+720,,0.1239,0.0015,
+722.5,,0.106325,0.001875,
+725,,0.08875,0.00225,
+727.5,,0.071175,0.002625,
+730,,0.0536,0.003,
+732.5,,0.045375,0.005075,
+735,,0.03715,0.00715,
+737.5,,0.028925,0.009225,
+740,,0.0207,0.0113,
+742.5,,0.0175,0.02105,
+745,,0.0143,0.0308,
+747.5,,0.0111,0.04055,
+750,,0.0079,0.0503,
+752.5,,0.006925,0.080575,
+755,,0.00595,0.11085,
+757.5,,0.004975,0.141125,
+760,,0.004,0.1714,
+762.5,,0.003375,0.2271,
+765,,0.00275,0.2828,
+767.5,,0.002125,0.3385,
+770,,0.0015,0.3942,
+772.5,,0.0015,0.458425,
+775,,0.0015,0.52265,
+777.5,,0.0015,0.586875,
+780,,0.0015,0.6511,
+782.5,,4.2632564E-017,0.6966,
+785,,,0.7421,
+787.5,,,0.7876,
+790,,,0.8331,
+792.5,,,0.8607,
+795,,,0.8883,
+797.5,,,0.9159,
+800,,,0.9435,
+802.5,,,0.955175,
+805,,,0.96685,
+807.5,,,0.978525,
+810,,,0.9902,
+812.5,,,0.99265,
+815,,,0.9951,
+817.5,,,0.99755,
+820,,,1,
+822.5,,,0.99985,
+825,,,0.9997,
+827.5,,,0.99955,
+830,,,0.9994,
+832.5,,,0.9984,
+835,,,0.9974,
+837.5,,,0.9964,
+840,,,0.9954,
+842.5,,,0.9927,
+845,,,0.99,
+847.5,,,0.9873,
+850,,,0.9846,
+852.5,,,0.9766,
+855,,,0.9686,
+857.5,,,0.9606,
+860,,,0.9526,
+862.5,,,0.944525,
+865,,,0.93645,
+867.5,,,0.928375,
+870,,,0.9203,
+872.5,,,0.9008,
+875,,,0.8813,
+877.5,,,0.8618,
+880,,,0.8423,
+882.5,,,0.789675,
+885,,,0.73705,
+887.5,,,0.684425,
+890,,,0.6318,
+892.5,,,0.55835,
+895,,,0.4849,
+897.5,,,0.41145,
+900,,,0.338,
+902.5,,,0.28935,
+905,,,0.2407,
+907.5,,,0.19205,
+910,,,0.1434,
+912.5,,,0.1222,
+915,,,0.101,
+917.5,,,0.0798,
+920,,,0.0586,
+922.5,,,0.05055,
+925,,,0.0425,
+927.5,,,0.03445,
+930,,,0.0264,
+932.5,,,0.0229,
+935,,,0.0194,
+937.5,,,0.0159,
+940,,,0.0124,
+942.5,,,0.0111,
+945,,,0.0098,
+947.5,,,0.0085,
+950,,,0.0072,
+952.5,,,0.00645,
+955,,,0.0057,
+957.5,,,0.00495,
+960,,,0.0042,
+962.5,,,0.003825,
+965,,,0.00345,
+967.5,,,0.003075,
+970,,,0.0027,
+972.5,,,0.002375,
+975,,,0.00205,
+977.5,,,0.001725,
+980,,,0.0014,
+982.5,,,0.001375,
+985,,,0.00135,
+987.5,,,0.001325,
+990,,,0.0013,
+992.5,,,0.001275,
+995,,,0.00125,
+997.5,,,0.001225,
+1000,,,0.0012,
+1002.5,,,0.001175,
+1005,,,0.00115,
+1007.5,,,0.001125,
+1010,,,0.0011,
+1012.5,,,0.001025,
+1015,,,0.00095,
+1017.5,,,0.000875,
+1020,,,0.0008,
+1022.5,,,0.000725,
+1025,,,0.00065,
+1027.5,,,0.000575,
+1030,,,0.0005,
+1032.5,,,,
+1035,,,,
+1037.5,,,,
+1040,,,,
+1042.5,,,,
+1045,,,,
+1047.5,,,,
+1050,,,,
+1052.5,,,,
+1055,,,,
+1057.5,,,,
+1060,,,,
+1062.5,,,,
+1065,,,,
+1067.5,,,,
+1070,,,,
+1072.5,,,,
+1075,,,,
+1077.5,,,,
+1080,,,,
+1082.5,,,,
+1085,,,,
+1087.5,,,,
+1090,,,,
+1092.5,,,,
+1095,,,,
+1097.5,,,,
+1100,,,,
+1102.5,,,,
+1105,,,,
+1107.5,,,,
+1110,,,,
+1112.5,,,,
+1115,,,,
+1117.5,,,,
+1120,,,,
+1122.5,,,,
+1125,,,,
+1127.5,,,,
+1130,,,,
+1132.5,,,,
+1135,,,,
+1137.5,,,,
+1140,,,,
+1142.5,,,,
+1145,,,,
+1147.5,,,,
+1150,,,,
+1152.5,,,,
+1155,,,,
+1157.5,,,,
+1160,,,,
+1162.5,,,,
+1165,,,,
+1167.5,,,,
+1170,,,,
+1172.5,,,,
+1175,,,,
+1177.5,,,,
+1180,,,,
+1182.5,,,,
+1185,,,,
+1187.5,,,,
+1190,,,,
+1192.5,,,,
+1195,,,,
+1197.5,,,,
+1200,,,,
+1202.5,,,,
+1205,,,,
+1207.5,,,,
+1210,,,,
+1212.5,,,,
+1215,,,,
+1217.5,,,,
+1220,,,,
+1222.5,,,,
+1225,,,,
+1227.5,,,,
+1230,,,,
+1232.5,,,,
+1235,,,,
+1237.5,,,,
+1240,,,,
+1242.5,,,,
+1245,,,,
+1247.5,,,,
+1250,,,,
+1252.5,,,,
+1255,,,,
+1257.5,,,,
+1260,,,,
+1262.5,,,,
+1265,,,,
+1267.5,,,,
+1270,,,,
+1272.5,,,,
+1275,,,,
+1277.5,,,,
+1280,,,,
+1282.5,,,,
+1285,,,,
+1287.5,,,,
+1290,,,,
+1292.5,,,,
+1295,,,,
+1297.5,,,,
+1300,,,,
+1302.5,,,,
+1305,,,,
+1307.5,,,,
+1310,,,,
+1312.5,,,,
+1315,,,,
+1317.5,,,,
+1320,,,,
+1322.5,,,,
+1325,,,,
+1327.5,,,,
+1330,,,,
+1332.5,,,,
+1335,,,,
+1337.5,,,,
+1340,,,,
+1342.5,,,,
+1345,,,,
+1347.5,,,,
+1350,,,,
+1352.5,,,,
+1355,,,,
+1357.5,,,,
+1360,,,,
+1362.5,,,,
+1365,,,,
+1367.5,,,,
+1370,,,,
+1372.5,,,,
+1375,,,,
+1377.5,,,,
+1380,,,,
+1382.5,,,,
+1385,,,,
+1387.5,,,,
+1390,,,,
+1392.5,,,,
+1395,,,,
+1397.5,,,,
+1400,,,,
+1402.5,,,,
+1405,,,,
+1407.5,,,,
+1410,,,,
+1412.5,,,,
+1415,,,,
+1417.5,,,,
+1420,,,,
+1422.5,,,,
+1425,,,,
+1427.5,,,,
+1430,,,,
+1432.5,,,,
+1435,,,,
+1437.5,,,,
+1440,,,,
+1442.5,,,,
+1445,,,,
+1447.5,,,,
+1450,,,,0.001
+1452.5,,,,0.001
+1455,,,,0.001
+1457.5,,,,0.001
+1460,,,,0.001
+1462.5,,,,0.001075
+1465,,,,0.00115
+1467.5,,,,0.001225
+1470,,,,0.0013
+1472.5,,,,0.001375
+1475,,,,0.00145
+1477.5,,,,0.001525
+1480,,,,0.0016
+1482.5,,,,0.00165
+1485,,,,0.0017
+1487.5,,,,0.00175
+1490,,,,0.0018
+1492.5,,,,0.00195
+1495,,,,0.0021
+1497.5,,,,0.00225
+1500,,,,0.0024
+1502.5,,,,0.0026
+1505,,,,0.0028
+1507.5,,,,0.003
+1510,,,,0.0032
+1512.5,,,,0.0035
+1515,,,,0.0038
+1517.5,,,,0.0041
+1520,,,,0.0044
+1522.5,,,,0.00495
+1525,,,,0.0055
+1527.5,,,,0.00605
+1530,,,,0.0066
+1532.5,,,,0.007475
+1535,,,,0.00835
+1537.5,,,,0.009225
+1540,,,,0.0101
+1542.5,,,,0.01185
+1545,,,,0.0136
+1547.5,,,,0.01535
+1550,,,,0.0171
+1552.5,,,,0.0203
+1555,,,,0.0235
+1557.5,,,,0.0267
+1560,,,,0.0299
+1562.5,,,,0.036625
+1565,,,,0.04335
+1567.5,,,,0.050075
+1570,,,,0.0568
+1572.5,,,,0.070825
+1575,,,,0.08485
+1577.5,,,,0.098875
+1580,,,,0.1129
+1582.5,,,,0.14065
+1585,,,,0.1684
+1587.5,,,,0.19615
+1590,,,,0.2239
+1592.5,,,,0.26905
+1595,,,,0.3142
+1597.5,,,,0.35935
+1600,,,,0.4045
+1602.5,,,,0.46765
+1605,,,,0.5308
+1607.5,,,,0.59395
+1610,,,,0.6571
+1612.5,,,,0.70995
+1615,,,,0.7628
+1617.5,,,,0.81565
+1620,,,,0.8685
+1622.5,,,,0.893425
+1625,,,,0.91835
+1627.5,,,,0.943275
+1630,,,,0.9682
+1632.5,,,,0.97615
+1635,,,,0.9841
+1637.5,,,,0.99205
+1640,,,,1
+1642.5,,,,0.997075
+1645,,,,0.99415
+1647.5,,,,0.991225
+1650,,,,0.9883
+1652.5,,,,0.983975
+1655,,,,0.97965
+1657.5,,,,0.975325
+1660,,,,0.971
+1662.5,,,,0.951725
+1665,,,,0.93245
+1667.5,,,,0.913175
+1670,,,,0.8939
+1672.5,,,,0.871325
+1675,,,,0.84875
+1677.5,,,,0.826175
+1680,,,,0.8036
+1682.5,,,,0.7457
+1685,,,,0.6878
+1687.5,,,,0.6299
+1690,,,,0.572
+1692.5,,,,0.5109
+1695,,,,0.4498
+1697.5,,,,0.3887
+1700,,,,0.3276
+1702.5,,,,0.2926
+1705,,,,0.2576
+1707.5,,,,0.2226
+1710,,,,0.1876
+1712.5,,,,0.17035
+1715,,,,0.1531
+1717.5,,,,0.13585
+1720,,,,0.1186
+1722.5,,,,0.107825
+1725,,,,0.09705
+1727.5,,,,0.086275
+1730,,,,0.0755
+1732.5,,,,0.068475
+1735,,,,0.06145
+1737.5,,,,0.054425
+1740,,,,0.0474
+1742.5,,,,0.04295
+1745,,,,0.0385
+1747.5,,,,0.03405
+1750,,,,0.0296
+1752.5,,,,0.026875
+1755,,,,0.02415
+1757.5,,,,0.021425
+1760,,,,0.0187
+1762.5,,,,0.0172
+1765,,,,0.0157
+1767.5,,,,0.0142
+1770,,,,0.0127
+1772.5,,,,0.011825
+1775,,,,0.01095
+1777.5,,,,0.010075
+1780,,,,0.0092
+1782.5,,,,0.008625
+1785,,,,0.00805
+1787.5,,,,0.007475
+1790,,,,0.0069
+1792.5,,,,0.006425
+1795,,,,0.00595
+1797.5,,,,0.005475
+1800,,,,0.005

imagery/i.atcorr/sensors_csv/VGT2_spot5.csv

@@ -0,0 +1,562 @@
+WL(nm),B0,B2,B3,MIR
+400,0.0007,,,
+402.5,0.00085,,,
+405,0.001,,,
+407.5,0.00105,,,
+410,0.0011,,,
+412.5,0.0012,,,
+415,0.0013,,,
+417.5,0.0029,,,
+420,0.0045,,,
+422.5,0.02095,,,
+425,0.0374,,,
+427.5,0.09635,,,
+430,0.1553,,,
+432.5,0.2551,,,
+435,0.3549,,,
+437.5,0.45365,,,
+440,0.5524,,,
+442.5,0.63175,,,
+445,0.7111,,,
+447.5,0.7602,,,
+450,0.8093,,,
+452.5,0.83035,,,
+455,0.8514,,,
+457.5,0.8698,,,
+460,0.8882,,,
+462.5,0.9441,,,
+465,1,,,
+467.5,0.9787,,,
+470,0.9574,,,
+472.5,0.7493,,,
+475,0.5412,,,
+477.5,0.3553,,,
+480,0.1694,,,
+482.5,0.1083,,,
+485,0.0472,,,
+487.5,0.03175,,,
+490,0.0163,,,
+492.5,0.0116,,,
+495,0.0069,,,
+497.5,0.00525,,,
+500,0.0036,,,
+502.5,0.0028,,,
+505,0.002,,,
+507.5,0.0017,,,
+510,0.0014,,,
+512.5,0.00115,,,
+515,0.0009,,,
+517.5,0.00085,,,
+520,0.0008,,,
+522.5,0.00075,,,
+525,0.0007,,,
+527.5,0.00065,,,
+530,0.0006,,,
+532.5,0.00055,,,
+535,0.0005,,,
+537.5,0.00045,,,
+540,0.0004,,,
+542.5,0.0004,,,
+545,0.0004,,,
+547.5,0.00035,,,
+550,0.0003,,,
+552.5,,,,
+555,,,,
+557.5,,,,
+560,,,,
+562.5,,,,
+565,,,,
+567.5,,,,
+570,,,,
+572.5,,,,
+575,,,,
+577.5,,,,
+580,,0.0006,,
+582.5,,0.000625,,
+585,,0.00065,,
+587.5,,0.000675,,
+590,,0.0007,,
+592.5,,0.003275,,
+595,,0.00585,,
+597.5,,0.008425,,
+600,,0.011,,
+602.5,,0.06875,,
+605,,0.1265,,
+607.5,,0.18425,,
+610,,0.242,,
+612.5,,0.35315,,
+615,,0.4643,,
+617.5,,0.57545,,
+620,,0.6866,,
+622.5,,0.730575,,
+625,,0.77455,,
+627.5,,0.818525,,
+630,,0.8625,,
+632.5,,0.87655,,
+635,,0.8906,,
+637.5,,0.90465,,
+640,,0.9187,,
+642.5,,0.930675,,
+645,,0.94265,,
+647.5,,0.954625,,
+650,,0.9666,,
+652.5,,0.97495,,
+655,,0.9833,,
+657.5,,0.99165,,
+660,,1,,
+662.5,,0.995725,,
+665,,0.99145,,
+667.5,,0.987175,,
+670,,0.9829,,
+672.5,,0.9387,,
+675,,0.8945,,
+677.5,,0.8503,,
+680,,0.8061,,
+682.5,,0.73155,,
+685,,0.657,,
+687.5,,0.58245,,
+690,,0.5079,,
+692.5,,0.446775,,
+695,,0.38565,,
+697.5,,0.324525,,
+700,,0.2634,0.0003,
+702.5,,0.2293,0.0003,
+705,,0.1952,0.0003,
+707.5,,0.1611,0.0003,
+710,,0.127,0.0003,
+712.5,,0.10955,0.000325,
+715,,0.0921,0.00035,
+717.5,,0.07465,0.000375,
+720,,0.0572,0.0004,
+722.5,,0.049,0.000525,
+725,,0.0408,0.00065,
+727.5,,0.0326,0.000775,
+730,,0.0244,0.0009,
+732.5,,0.020775,0.001675,
+735,,0.01715,0.00245,
+737.5,,0.013525,0.003225,
+740,,0.0099,0.004,
+742.5,,0.0085,0.0082,
+745,,0.0071,0.0124,
+747.5,,0.0057,0.0166,
+750,,0.0043,0.0208,
+752.5,,0.0037,0.0361,
+755,,0.0031,0.0514,
+757.5,,0.0025,0.0667,
+760,,0.0019,0.082,
+762.5,,0.001675,0.1175,
+765,,0.00145,0.153,
+767.5,,0.001225,0.1885,
+770,,0.001,0.224,
+772.5,,0.0009,0.27745,
+775,,0.0008,0.3309,
+777.5,,0.0007,0.38435,
+780,,0.0006,0.4378,
+782.5,,,0.4906,
+785,,,0.5434,
+787.5,,,0.5962,
+790,,,0.649,
+792.5,,,0.689875,
+795,,,0.73075,
+797.5,,,0.771625,
+800,,,0.8125,
+802.5,,,0.836875,
+805,,,0.86125,
+807.5,,,0.885625,
+810,,,0.91,
+812.5,,,0.9248,
+815,,,0.9396,
+817.5,,,0.9544,
+820,,,0.9692,
+822.5,,,0.976075,
+825,,,0.98295,
+827.5,,,0.989825,
+830,,,0.9967,
+832.5,,,0.997525,
+835,,,0.99835,
+837.5,,,0.999175,
+840,,,1,
+842.5,,,0.9939,
+845,,,0.9878,
+847.5,,,0.9817,
+850,,,0.9756,
+852.5,,,0.9618,
+855,,,0.948,
+857.5,,,0.9342,
+860,,,0.9204,
+862.5,,,0.903,
+865,,,0.8856,
+867.5,,,0.8682,
+870,,,0.8508,
+872.5,,,0.82365,
+875,,,0.7965,
+877.5,,,0.76935,
+880,,,0.7422,
+882.5,,,0.69445,
+885,,,0.6467,
+887.5,,,0.59895,
+890,,,0.5512,
+892.5,,,0.491,
+895,,,0.4308,
+897.5,,,0.3706,
+900,,,0.3104,
+902.5,,,0.268125,
+905,,,0.22585,
+907.5,,,0.183575,
+910,,,0.1413,
+912.5,,,0.1213,
+915,,,0.1013,
+917.5,,,0.0813,
+920,,,0.0613,
+922.5,,,0.053,
+925,,,0.0447,
+927.5,,,0.0364,
+930,,,0.0281,
+932.5,,,0.024475,
+935,,,0.02085,
+937.5,,,0.017225,
+940,,,0.0136,
+942.5,,,0.012025,
+945,,,0.01045,
+947.5,,,0.008875,
+950,,,0.0073,
+952.5,,,0.00645,
+955,,,0.0056,
+957.5,,,0.00475,
+960,,,0.0039,
+962.5,,,0.00355,
+965,,,0.0032,
+967.5,,,0.00285,
+970,,,0.0025,
+972.5,,,0.002275,
+975,,,0.00205,
+977.5,,,0.001825,
+980,,,0.0016,
+982.5,,,0.001475,
+985,,,0.00135,
+987.5,,,0.001225,
+990,,,0.0011,
+992.5,,,0.000975,
+995,,,0.00085,
+997.5,,,0.000725,
+1000,,,0.0006,
+1002.5,,,,
+1005,,,,
+1007.5,,,,
+1010,,,,
+1012.5,,,,
+1015,,,,
+1017.5,,,,
+1020,,,,
+1022.5,,,,
+1025,,,,
+1027.5,,,,
+1030,,,,
+1032.5,,,,
+1035,,,,
+1037.5,,,,
+1040,,,,
+1042.5,,,,
+1045,,,,
+1047.5,,,,
+1050,,,,
+1052.5,,,,
+1055,,,,
+1057.5,,,,
+1060,,,,
+1062.5,,,,
+1065,,,,
+1067.5,,,,
+1070,,,,
+1072.5,,,,
+1075,,,,
+1077.5,,,,
+1080,,,,
+1082.5,,,,
+1085,,,,
+1087.5,,,,
+1090,,,,
+1092.5,,,,
+1095,,,,
+1097.5,,,,
+1100,,,,
+1102.5,,,,
+1105,,,,
+1107.5,,,,
+1110,,,,
+1112.5,,,,
+1115,,,,
+1117.5,,,,
+1120,,,,
+1122.5,,,,
+1125,,,,
+1127.5,,,,
+1130,,,,
+1132.5,,,,
+1135,,,,
+1137.5,,,,
+1140,,,,
+1142.5,,,,
+1145,,,,
+1147.5,,,,
+1150,,,,
+1152.5,,,,
+1155,,,,
+1157.5,,,,
+1160,,,,
+1162.5,,,,
+1165,,,,
+1167.5,,,,
+1170,,,,
+1172.5,,,,
+1175,,,,
+1177.5,,,,
+1180,,,,
+1182.5,,,,
+1185,,,,
+1187.5,,,,
+1190,,,,
+1192.5,,,,
+1195,,,,
+1197.5,,,,
+1200,,,,
+1202.5,,,,
+1205,,,,
+1207.5,,,,
+1210,,,,
+1212.5,,,,
+1215,,,,
+1217.5,,,,
+1220,,,,
+1222.5,,,,
+1225,,,,
+1227.5,,,,
+1230,,,,
+1232.5,,,,
+1235,,,,
+1237.5,,,,
+1240,,,,
+1242.5,,,,
+1245,,,,
+1247.5,,,,
+1250,,,,
+1252.5,,,,
+1255,,,,
+1257.5,,,,
+1260,,,,
+1262.5,,,,
+1265,,,,
+1267.5,,,,
+1270,,,,
+1272.5,,,,
+1275,,,,
+1277.5,,,,
+1280,,,,
+1282.5,,,,
+1285,,,,
+1287.5,,,,
+1290,,,,
+1292.5,,,,
+1295,,,,
+1297.5,,,,
+1300,,,,
+1302.5,,,,
+1305,,,,
+1307.5,,,,
+1310,,,,
+1312.5,,,,
+1315,,,,
+1317.5,,,,
+1320,,,,
+1322.5,,,,
+1325,,,,
+1327.5,,,,
+1330,,,,
+1332.5,,,,
+1335,,,,
+1337.5,,,,
+1340,,,,
+1342.5,,,,
+1345,,,,
+1347.5,,,,
+1350,,,,
+1352.5,,,,
+1355,,,,
+1357.5,,,,
+1360,,,,
+1362.5,,,,
+1365,,,,
+1367.5,,,,
+1370,,,,
+1372.5,,,,
+1375,,,,
+1377.5,,,,
+1380,,,,
+1382.5,,,,
+1385,,,,
+1387.5,,,,
+1390,,,,
+1392.5,,,,
+1395,,,,
+1397.5,,,,
+1400,,,,
+1402.5,,,,
+1405,,,,
+1407.5,,,,
+1410,,,,
+1412.5,,,,
+1415,,,,
+1417.5,,,,
+1420,,,,
+1422.5,,,,
+1425,,,,
+1427.5,,,,
+1430,,,,
+1432.5,,,,
+1435,,,,
+1437.5,,,,
+1440,,,,
+1442.5,,,,
+1445,,,,
+1447.5,,,,
+1450,,,,0.0005
+1452.5,,,,0.000575
+1455,,,,0.00065
+1457.5,,,,0.000725
+1460,,,,0.0008
+1462.5,,,,0.000825
+1465,,,,0.00085
+1467.5,,,,0.000875
+1470,,,,0.0009
+1472.5,,,,0.0009
+1475,,,,0.0009
+1477.5,,,,0.0009
+1480,,,,0.0009
+1482.5,,,,0.00115
+1485,,,,0.0014
+1487.5,,,,0.00165
+1490,,,,0.0019
+1492.5,,,,0.00215
+1495,,,,0.0024
+1497.5,,,,0.00265
+1500,,,,0.0029
+1502.5,,,,0.003425
+1505,,,,0.00395
+1507.5,,,,0.004475
+1510,,,,0.005
+1512.5,,,,0.006
+1515,,,,0.007
+1517.5,,,,0.008
+1520,,,,0.009
+1522.5,,,,0.010775
+1525,,,,0.01255
+1527.5,,,,0.014325
+1530,,,,0.0161
+1532.5,,,,0.0196
+1535,,,,0.0231
+1537.5,,,,0.0266
+1540,,,,0.0301
+1542.5,,,,0.037225
+1545,,,,0.04435
+1547.5,,,,0.051475
+1550,,,,0.0586
+1552.5,,,,0.0723
+1555,,,,0.086
+1557.5,,,,0.0997
+1560,,,,0.1134
+1562.5,,,,0.140125
+1565,,,,0.16685
+1567.5,,,,0.193575
+1570,,,,0.2203
+1572.5,,,,0.26685
+1575,,,,0.3134
+1577.5,,,,0.35995
+1580,,,,0.4065
+1582.5,,,,0.466425
+1585,,,,0.52635
+1587.5,,,,0.586275
+1590,,,,0.6462
+1592.5,,,,0.692875
+1595,,,,0.73955
+1597.5,,,,0.786225
+1600,,,,0.8329
+1602.5,,,,0.863025
+1605,,,,0.89315
+1607.5,,,,0.923275
+1610,,,,0.9534
+1612.5,,,,0.964375
+1615,,,,0.97535
+1617.5,,,,0.986325
+1620,,,,0.9973
+1622.5,,,,0.997975
+1625,,,,0.99865
+1627.5,,,,0.999325
+1630,,,,1
+1632.5,,,,0.998725
+1635,,,,0.99745
+1637.5,,,,0.996175
+1640,,,,0.9949
+1642.5,,,,0.991425
+1645,,,,0.98795
+1647.5,,,,0.984475
+1650,,,,0.981
+1652.5,,,,0.968675
+1655,,,,0.95635
+1657.5,,,,0.944025
+1660,,,,0.9317
+1662.5,,,,0.905525
+1665,,,,0.87935
+1667.5,,,,0.853175
+1670,,,,0.827
+1672.5,,,,0.781075
+1675,,,,0.73515
+1677.5,,,,0.689225
+1680,,,,0.6433
+1682.5,,,,0.574575
+1685,,,,0.50585
+1687.5,,,,0.437125
+1690,,,,0.3684
+1692.5,,,,0.324375
+1695,,,,0.28035
+1697.5,,,,0.236325
+1700,,,,0.1923
+1702.5,,,,0.171275
+1705,,,,0.15025
+1707.5,,,,0.129225
+1710,,,,0.1082
+1712.5,,,,0.09705
+1715,,,,0.0859
+1717.5,,,,0.07475
+1720,,,,0.0636
+1722.5,,,,0.0575
+1725,,,,0.0514
+1727.5,,,,0.0453
+1730,,,,0.0392
+1732.5,,,,0.03595
+1735,,,,0.0327
+1737.5,,,,0.02945
+1740,,,,0.0262
+1742.5,,,,0.023525
+1745,,,,0.02085
+1747.5,,,,0.018175
+1750,,,,0.0155
+1752.5,,,,0.013875
+1755,,,,0.01225
+1757.5,,,,0.010625
+1760,,,,0.009
+1762.5,,,,0.008175
+1765,,,,0.00735
+1767.5,,,,0.006525
+1770,,,,0.0057
+1772.5,,,,0.00535
+1775,,,,0.005
+1777.5,,,,0.00465
+1780,,,,0.0043
+1782.5,,,,0.003975
+1785,,,,0.00365
+1787.5,,,,0.003325
+1790,,,,0.003
+1792.5,,,,0.002625
+1795,,,,0.00225
+1797.5,,,,0.001875
+1800,,,,0.0015