i.atcorr: filter functions updated for Rapideye (source: https://resa.planet.com/files/2014-06/Spectral_Response_Curves.pdf; contributed by craab) (trac https://trac.osgeo.org/grass/ticket/2545)

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

imagery/i.atcorr/iwave.cpp

@@ -2178,72 +2178,113 @@ void IWave::rapideye(int iwa)
 {
 
     /* "Blue" of rapideye */
-    static const float sr1[30] = {
-		.0200, .5050, .7300, .7450, .7600, .7750, .7800,
-		.8000, .8000, .8000, .8100, .8000, .8100, .8200,
-		.8300, .8400, .8500, .8650, .8800, .8900, .9200,
-		.9500, .9700, .9700, 1.0000, 1.0000
+    static const float sr1[128] = {
+		.0020, .0023, .0026, .0029, .0018, .0014, .0010, .0095,
+		.3210, .7250, .7400, .7590, .7700, .7810, .7840, .7935,
+		.7960, .8005, .8060, .8040, .8070, .8170, .8200, .8275,
+		.8400, .8470, .8620, .8765, .8860, .9105, .9280, .9415,
+		.9690, .9685, 1.0000, .9875, .4370, .0190, .0010, .0010,
+		.0010, .0010, .0010, .0010, .0010, .0010, .0010, .0010,
+		.0010, .0010, .0010, .0010, .0010, .0010, .0010, .0010,
+		.0010, .0010, .0010, .0010, .0010, .0010, .0010, .0010,
+		.0010, .0010, .0010, .0010, .0010, .0010, .0010, .0010,
+		.0010, .0010, .0010, .0010, .0010, .0010, .0010, .0010,
+		.0010, .0010, .0010, .0010, .0010, .0010, .0010, .0010,
+		.0010, .0010, .0010, .0010, .0010, .0010, .0010, .0010,
+		.0010, .0010, .0010, .0010, .0010, .0010, .0010, .0010,
+		.0010, .0010, .0010, .0010, .0010, .0010, .0010, .0010,
+		.0010, .0010, .0010, .0010, .0010, .0010, .0010, .0010,
+		.0010, .0010, .0010, .0010, .0010, .0010, .0010, .0010,
+		
     };
-
+	
     /* "Green" of rapideye */
-    static const float sr2[53] = {
-		.0100, .0183, .0267, .0105, .0132, .0158, .0184,
-		.0237, .0263, .0289, .0288, .0268, .0248, .0228,
-		.0188, .0168, .0148, .0128, .0108, .0200, .1200,
-		.8700, .8700, .8800, .8700, .8800, .8800, .8800,
-		.8900, .9000, .9000, .9050, .9200, .9250, .9300,
-		.9500, .9550, .9700, .9750, .9800, .9800, .9900,
-		.9900, .9850, 1.0000, .2750
+    static const float sr2[181] = {
+		.0010, .0010, .0010, .0010, .0010, .0010, .0010, .0010,
+		.0010, .0010, .0010, .0010, .0010, .0010, .0010, .0015,
+		.0030, .0015, .0340, .0010, .0010, .0010, .0010, .0010,
+		.0010, .0010, .0270, .0019, .0018, .0017, .0015, .0014,
+		.0013, .0011, .0010, .0020, .0070, .0400, .2720, .8260,
+		.8640, .8750, .8730, .8730, .8810, .8825, .8800, .8845,
+		.8940, .8985, .8980, .9075, .9200, .9275, .9330, .9440,
+		.9520, .9575, .9700, .9765, .9740, .9855, .9930, .9830,
+		1.0000, .9730, .8900, .0825, .0060, .0010, .0010, .0010,
+		.0010, .0010, .0010, .0010, .0010, .0010, .0010, .0010,
+		.0010, .0010, .0010, .0010, .0010, .0010, .0010, .0010,
+		.0010, .0010, .0010, .0010, .0010, .0010, .0010, .0010,
+		.0010, .0010, .0010, .0010, .0010, .0010, .0010, .0010,
+		.0010, .0010, .0010, .0010, .0010, .0010, .0010, .0010,
+		.0010, .0010, .0010, .0010, .0010, .0010, .0010, .0010,
+		.0010, .0010, .0010, .0010, .0010, .0010, .0010, .0010,
+		.0010, .0010, .0010, .0010, .0010, .0010, .0010, .0010,
+		.0010, .0010, .0010, .0010, .0010, .0010, .0010, .0010,
+		.0010, .0010, .0010, .0010, .0010, .0010, .0010, .0010,
+		.0010, .0010, .0010, .0010, .0010, .0010, .0010, .0010,
+		.0010, .0010, .0010, .0010, .0010, .0010, .0010, .0010,
+		.0010, .0010, .0010, .0010, .0010, .0010, .0010, .0010,
+		.0010, .0010, .0010, .0010, .0010
     };
-
+	
     /* "Red" of rapideye */
-    static const float sr3[27] = {
-		.0100, .0300, .2000, .7900, .8500, .8750, .8800,
-		.9000, .9100, .9100, .9200, .9300, .9450, .9500,
-		.9700, .9750, .9900, .9900, .9900, .9900, 1.0000,
-		.8600, .1150
+    static const float sr3[94] = {
+		.0010, .0011, .0012, .0013, .0014, .0015, .0016, .0017,
+		.0018, .0019, .0020, .0019, .0017, .0016, .0014, .0012,
+		.0011, .0010, .0010, .0010, .0010, .0010, .0010, .0010,
+		.0010, .0010, .0010, .0010, .0010, .0010, .0010, .0010,
+		.0010, .0010, .0010, .0010, .0010, .0010, .0010, .0010,
+		.0010, .0010, .0010, .0010, .0010, .0010, .0010, .0010,
+		.0010, .0010, .0010, .0010, .0020, .0019, .0018, .0018,
+		.0017, .0016, .0015, .0014, .0013, .0013, .0012, .0011,
+		.0010, .0025, .0090, .0620, .4370, .8595, .8540, .8785,
+		.8810, .8875, .9020, .9105, .9150, .9255, .9390, .9490,
+		.9580, .9645, .9700, .9785, .9890, .9900, .9850, .9970,
+		.9890, .9940, .4630, .0430, .0050, .0010
     };
-
+	
     /* "RedEdge " of rapideye */
-    static const float sr4[95] = {
-		.0300, .0285, .0271, .0256, .0241, .0226, .0212,
-		.0182, .0168, .0153, .0138, .0124, .0109, .0101,
-		.0108, .0111, .0114, .0118, .0121, .0124, .0128,
-		.0134, .0138, .0141, .0144, .0147, .0151, .0154,
-		.0161, .0164, .0167, .0170, .0174, .0177, .0180,
-		.0187, .0190, .0193, .0197, .0200, .0200, .0200,
-		.0200, .0200, .0200, .0200, .0200, .0200, .0200,
-		.0200, .0200, .0197, .0191, .0185, .0178, .0172,
-		.0159, .0153, .0146, .0140, .0133, .0127, .0121,
-		.0108, .0101, .0200, .1050, .4900, .9450, 1.0000,
-		1.0000, .9900, .9900, .9900, .9800, .9750, .9700,
-		.9600, .9450, .9400, .9000, .4300, .0950
+    static const float sr4[98] = {
+		.0270, .0247, .0224, .0200, .0177, .0154, .0131, .0107,
+		.0084, .0061, .0038, .0015, .0010, .0010, .0010, .0010,
+		.0010, .0010, .0010, .0010, .0017, .0025, .0034, .0042,
+		.0051, .0059, .0068, .0076, .0085, .0093, .0101, .0110,
+		.0118, .0127, .0135, .0144, .0152, .0161, .0169, .0178,
+		.0186, .0195, .0203, .0212, .0220, .0195, .0170, .0145,
+		.0120, .0095, .0070, .0045, .0020, .0016, .0047, .0079,
+		.0110, .0141, .0010, .0010, .0010, .0010, .0010, .0010,
+		.0010, .0010, .0010, .0010, .0010, .0010, .0010, .0010,
+		.0020, .0055, .0210, .1085, .4910, .9490, .9980, .9995,
+		.9980, .9915, .9870, .9855, .9820, .9740, .9660, .9640,
+		.9610, .9455, .9390, .8980, .4250, .0940, .0200, .0055,
+		.0020, .0010
     };
-
+	
     /* "NIR " of rapideye */
-    static const float sr5[137] = {
-		.0100, .0102, .0103, .0105, .0106, .0108, .0110,
-		.0113, .0115, .0116, .0118, .0119, .0121, .0123,
-		.0126, .0128, .0129, .0131, .0132, .0134, .0136,
-		.0139, .0141, .0142, .0144, .0145, .0147, .0149,
-		.0152, .0154, .0155, .0157, .0158, .0160, .0162,
-		.0165, .0167, .0168, .0170, .0171, .0173, .0175,
-		.0178, .0180, .0181, .0183, .0184, .0186, .0188,
-		.0191, .0193, .0194, .0196, .0197, .0199, .0199,
-		.0192, .0188, .0185, .0182, .0178, .0175, .0171,
-		.0164, .0161, .0158, .0154, .0151, .0147, .0144,
-		.0137, .0134, .0130, .0127, .0123, .0120, .0116,
-		.0110, .0106, .0103, .0100, .0100, .0250, .0700,
-		.5000, .8850, 1.0000, .9750, .9700, .9750, .9800,
-		.9700, .9600, .9600, .9600, .9600, .9600, .9600,
-		.9500, .9400, .9400, .9300, .9300, .9300, .9300,
-		.9200, .9200, .9100, .9000, .8900, .8850, .8800,
-		.8600, .8400, .8400, .8150, .5800, .2450, .0800,
-
-    };
-
-    static const float wli[5] = {0.438, 0.463, 0.624, 0.500, 0.520};
-    static const float wls[5] = {0.513, 0.594, 0.690, 0.737, 0.862};
+    static const float sr5[160] = {
+		.0010, .0010, .0010, .0010, .0010, .0010, .0010, .0010,
+		.0010, .0010, .0010, .0010, .0010, .0010, .0010, .0012,
+		.0014, .0015, .0017, .0019, .0110, .0010, .0010, .0010,
+		.0010, .0010, .0010, .0010, .0013, .0016, .0019, .0020,
+		.0019, .0018, .0017, .0016, .0015, .0014, .0013, .0012,
+		.0011, .0010, .0010, .0011, .0012, .0013, .0014, .0015,
+		.0016, .0016, .0017, .0018, .0019, .0020, .0010, .0010,
+		.0010, .0010, .0010, .0010, .0010, .0010, .0030, .0010,
+		.0010, .0010, .0010, .0010, .0010, .0010, .0020, .0010,
+		.0010, .0010, .0010, .0010, .0010, .0010, .0010, .0010,
+		.0010, .0015, .0030, .0010, .0010, .0010, .0010, .0010,
+		.0010, .0010, .0010, .0010, .0010, .0010, .0010, .0010,
+		.0010, .0010, .0010, .0010, .0010, .0010, .0010, .0010,
+		.0010, .0010, .0010, .0010, .0010, .0015, .0030, .0055,
+		.0120, .0275, .0680, .1920, .4970, .8825, 1.0000, .9775,
+		.9670, .9760, .9830, .9795, .9700, .9620, .9580, .9575,
+		.9580, .9585, .9560, .9535, .9480, .9425, .9370, .9320,
+		.9280, .9270, .9260, .9255, .9230, .9165, .9060, .8965,
+		.8890, .8845, .8820, .8735, .8570, .8400, .8400, .8135,
+		.5820, .2440, .0800, .0285, .0110, .0050, .0020, .0010,
+		
+    };
+	
+    static const float wli[5] = {0.420, 0.424, 0.460, 0.500, 0.470};
+    static const float wls[5] = {0.739, 0.876, 0.693, 0.743, 0.870};
 
     ffu.wlinf = (float)wli[iwa-1];
     ffu.wlsup = (float)wls[iwa-1];
@@ -2253,15 +2294,15 @@ void IWave::rapideye(int iwa)
 
     switch(iwa)
     {
-    case 1: for(i = 0; i < 30; i++)  ffu.s[75+i] = sr1[i];
+    case 1: for(i = 0; i < 128; i++)  ffu.s[68+i] = sr1[i];
         break;
-    case 2: for(i = 0; i < 53; i++)  ffu.s[85+i] = sr2[i];
+    case 2: for(i = 0; i < 181; i++)  ffu.s[69+i] = sr2[i];
         break;
-    case 3: for(i = 0; i < 27; i++)  ffu.s[149+i] = sr3[i];
+    case 3: for(i = 0; i < 94; i++)  ffu.s[84+i] = sr3[i];
         break;
-    case 4: for(i = 0; i < 95; i++)  ffu.s[100+i] = sr4[i];
+    case 4: for(i = 0; i < 98; i++)  ffu.s[100+i] = sr4[i];
         break;
-    case 5: for(i = 0; i < 137; i++)  ffu.s[108+i] = sr5[i];
+    case 5: for(i = 0; i < 160; i++)  ffu.s[88+i] = sr5[i];
         break;
     }
 }

imagery/i.atcorr/sensors_csv/rapideye.csv

@@ -1,482 +1,482 @@
 "Wavelength (nm)","Blue","Green","Red","RedEdge ","NIR "
-420,0,0,0,0,0
-421,0,0,0,0,0
-422,0,0,0,0,0
-423,0,0,0,0,0
-424,0,0,0,0,0
-425,0,0,0,0,0
-426,0,0,0,0,0
-427,0,0,0,0,0
-428,0,0,0,0,0
-429,0,0,0,0,0
-430,0,0,0,0,0
-431,0,0,0,0,0
-432,0,0,0,0,0
-433,0,0,0,0,0
-434,0,0,0,0,0
-435,0,0,0,0,0
-436,0,0,0,0,0
-437,0,0,0,0,0
-438,0.02,0,0,0,0
-439,0.07,0,0,0,0
-440,0.32,0,0,0,0
-441,0.69,0,0,0,0
-442,0.72,0,0,0,0
-443,0.73,0,0,0,0
-444,0.74,0,0,0,0
-445,0.74,0,0,0,0
-446,0.75,0,0,0,0
-447,0.76,0,0,0,0
-448,0.76,0,0,0,0
-449,0.77,0,0,0,0
-450,0.77,0,0,0,0
-451,0.78,0,0,0,0
-452,0.78,0,0,0,0
-453,0.78,0,0,0,0
-454,0.78,0,0,0,0
-455,0.78,0,0,0,0
-456,0.79,0,0,0,0
-457,0.79,0,0,0,0
-458,0.8,0,0,0,0
-459,0.8,0,0,0,0
-460,0.8,0,0,0,0
-461,0.8,0,0,0,0
-462,0.8,0,0,0,0
-463,0.8,0.01,0,0,0
-464,0.81,0,0,0,0
-465,0.81,0,0,0,0
-466,0.81,0,0,0,0
-467,0.8,0,0,0,0
-468,0.8,0,0,0,0
-469,0.81,0.03,0,0,0
-470,0.81,0.01,0,0,0
-471,0.81,0,0,0,0
-472,0.82,0,0,0,0
-473,0.82,0,0,0,0
-474,0.82,0,0,0,0
-475,0.82,0,0,0,0
-476,0.82,0,0,0,0
-477,0.83,0,0,0,0
-478,0.83,0,0,0,0
-479,0.84,0,0,0,0
-480,0.84,0,0,0,0
-481,0.84,0,0,0,0
-482,0.85,0,0,0,0
-483,0.85,0,0,0,0
-484,0.86,0,0,0,0
-485,0.86,0,0,0,0
-486,0.87,0,0,0,0
-487,0.88,0,0,0,0
-488,0.88,0,0,0,0
-489,0.88,0.03,0,0,0
-490,0.89,0,0,0,0
-491,0.89,0,0,0,0
-492,0.91,0,0,0,0
-493,0.92,0,0,0,0
-494,0.92,0,0,0,0
-495,0.93,0,0,0,0
-496,0.93,0,0,0,0
-497,0.94,0,0,0,0
-498,0.95,0,0,0,0
-499,0.96,0,0,0,0
-500,0.97,0,0,0.03,0
-501,0.97,0,0,0,0
-502,0.97,0,0,0,0
-503,0.97,0,0,0,0
-504,0.98,0,0,0,0
-505,1,0,0,0,0
-506,1,0,0,0,0
-507,0.98,0,0,0,0
-508,1,0,0,0,0
-509,0.96,0,0,0,0
-510,0.44,0,0,0,0
-511,0.11,0,0,0,0
-512,0.03,0,0,0,0
-513,0.01,0,0,0,0
-514,0,0.01,0,0,0
-515,0,0.01,0,0,0
-516,0,0.03,0,0,0
-517,0,0.05,0,0,0
-518,0,0.12,0,0,0
-519,0,0.27,0,0,0
-520,0,0.54,0,0,0.01
-521,0,0.78,0,0,0
-522,0,0.87,0,0,0
-523,0,0.87,0,0,0
-524,0,0.86,0,0,0
-525,0,0.87,0,0,0
-526,0,0.87,0,0,0
-527,0,0.88,0,0,0
-528,0,0.88,0,0,0
-529,0,0.87,0,0,0
-530,0,0.87,0,0,0
-531,0,0.87,0,0,0
-532,0,0.87,0,0,0
-533,0,0.88,0,0,0
-534,0,0.88,0,0.01,0
-535,0,0.88,0,0,0
-536,0,0.88,0,0,0
-537,0,0.88,0,0,0
-538,0,0.88,0,0,0
-539,0,0.88,0,0,0
-540,0,0.88,0,0,0
-541,0,0.88,0,0,0
-542,0,0.89,0,0,0
-543,0,0.89,0,0,0
-544,0,0.89,0,0,0
-545,0,0.9,0,0,0
-546,0,0.9,0,0,0
-547,0,0.9,0,0,0
-548,0,0.9,0,0,0
-549,0,0.9,0,0,0
-550,0,0.9,0,0,0
-551,0,0.91,0,0,0
-552,0,0.91,0,0,0
-553,0,0.92,0,0,0
-554,0,0.92,0,0,0
-555,0,0.92,0,0,0
-556,0,0.93,0,0,0
-557,0,0.93,0,0,0
-558,0,0.93,0,0,0
-559,0,0.93,0,0,0
-560,0,0.94,0,0,0
-561,0,0.94,0,0,0
-562,0,0.95,0,0,0
-563,0,0.95,0,0,0
-564,0,0.95,0,0,0
-565,0,0.95,0,0,0
-566,0,0.96,0,0,0
-567,0,0.96,0,0,0
-568,0,0.97,0,0,0
-569,0,0.97,0,0,0
-570,0,0.97,0,0,0
-571,0,0.98,0,0,0
-572,0,0.98,0,0,0
-573,0,0.98,0,0,0
-574,0,0.97,0,0,0
-575,0,0.98,0,0,0
-576,0,0.98,0,0,0
-577,0,0.99,0,0,0
-578,0,0.99,0,0,0
-579,0,0.99,0,0,0
-580,0,0.99,0,0,0
-581,0,0.98,0,0,0
-582,0,0.98,0,0,0
-583,0,0.99,0,0,0
-584,0,1,0,0,0
-585,0,0.99,0,0,0
-586,0,0.98,0,0,0
-587,0,0.97,0,0,0
-588,0,1,0,0,0
-589,0,0.89,0,0,0
-590,0,0.42,0,0,0
-591,0,0.13,0,0,0
-592,0,0.04,0,0,0
-593,0,0.01,0,0,0
-594,0,0.01,0,0,0
-595,0,0,0,0,0
-596,0,0,0,0,0
-597,0,0,0,0,0
-598,0,0,0,0,0
-599,0,0,0,0,0
-600,0,0,0,0,0
-601,0,0,0,0,0
-602,0,0,0,0,0
-603,0,0,0,0,0
-604,0,0,0,0,0
-605,0,0,0,0,0
-606,0,0,0,0,0
-607,0,0,0,0,0
-608,0,0,0,0,0
-609,0,0,0,0,0
-610,0,0,0,0.02,0
-611,0,0,0,0,0
-612,0,0,0,0,0
-613,0,0,0,0,0
-614,0,0,0,0,0
-615,0,0,0,0,0
-616,0,0,0,0,0
-617,0,0,0,0,0
-618,0,0,0,0,0
-619,0,0,0,0,0
-620,0,0,0,0,0
-621,0,0,0,0,0
-622,0,0,0,0,0
-623,0,0,0,0,0
-624,0,0,0.01,0,0
-625,0,0,0.01,0,0
-626,0,0,0.02,0,0
-627,0,0,0.04,0,0
-628,0,0,0.09,0,0
-629,0,0,0.2,0,0
-630,0,0,0.44,0,0
-631,0,0,0.72,0,0
-632,0,0,0.86,0,0
-633,0,0,0.86,0,0
-634,0,0,0.85,0,0
-635,0,0,0.85,0,0
-636,0,0,0.87,0,0
-637,0,0,0.88,0,0
-638,0,0,0.88,0,0
-639,0,0,0.88,0,0
-640,0,0,0.88,0,0
-641,0,0,0.88,0,0
-642,0,0,0.89,0,0
-643,0,0,0.89,0,0
-644,0,0,0.9,0.02,0
-645,0,0,0.9,0,0
-646,0,0,0.91,0,0
-647,0,0,0.91,0,0
-648,0,0,0.91,0,0
-649,0,0,0.91,0,0
-650,0,0,0.92,0,0
-651,0,0,0.92,0,0
-652,0,0,0.92,0,0
-653,0,0,0.93,0,0
-654,0,0,0.93,0,0
-655,0,0,0.94,0,0
-656,0,0,0.94,0,0
-657,0,0,0.95,0,0
-658,0,0,0.95,0,0
-659,0,0,0.95,0,0
-660,0,0,0.96,0,0
-661,0,0,0.96,0,0
-662,0,0,0.96,0,0
-663,0,0,0.97,0,0
-664,0,0,0.97,0,0
-665,0,0,0.97,0,0
-666,0,0,0.97,0,0
-667,0,0,0.98,0,0
-668,0,0,0.98,0,0
-669,0,0,0.99,0,0
-670,0,0,0.99,0,0
-671,0,0,0.99,0,0
-672,0,0,0.99,0,0
-673,0,0,0.99,0,0
-674,0,0,0.99,0,0.02
-675,0,0,0.99,0,0
-676,0,0,0.99,0,0
-677,0,0,0.99,0,0
-678,0,0,1,0,0
-679,0,0,1,0,0
-680,0,0,0.99,0,0
-681,0,0,0.98,0,0
-682,0,0,0.99,0,0
-683,0,0,1,0.01,0
-684,0,0,0.86,0.01,0
-685,0,0,0.46,0.02,0
-686,0,0,0.17,0.04,0
-687,0,0,0.06,0.07,0
-688,0,0,0.02,0.14,0
-689,0,0,0.01,0.28,0
-690,0,0,0.01,0.49,0
-691,0,0,0,0.74,0
-692,0,0,0,0.91,0
-693,0,0,0,0.98,0
-694,0,0,0,1,0
-695,0,0,0,1,0
-696,0,0,0,1,0
-697,0,0,0,1,0
-698,0,0,0,1,0
-699,0,0,0,1,0
-700,0,0,0,1,0
-701,0,0,0,1,0
-702,0,0,0,0.99,0
-703,0,0,0,0.99,0
-704,0,0,0,0.99,0
-705,0,0,0,0.99,0
-706,0,0,0,0.99,0
-707,0,0,0,0.99,0
-708,0,0,0,0.99,0
-709,0,0,0,0.98,0
-710,0,0,0,0.98,0
-711,0,0,0,0.98,0
-712,0,0,0,0.98,0
-713,0,0,0,0.97,0
-714,0,0,0,0.97,0
-715,0,0,0,0.97,0
-716,0,0,0,0.97,0
-717,0,0,0,0.96,0
-718,0,0,0,0.96,0
-719,0,0,0,0.96,0
-720,0,0,0,0.96,0
-721,0,0,0,0.96,0
-722,0,0,0,0.95,0
-723,0,0,0,0.94,0
-724,0,0,0,0.94,0
-725,0,0,0,0.94,0
-726,0,0,0,0.95,0
-727,0,0,0,0.94,0
-728,0,0,0,0.86,0
-729,0,0,0,0.67,0
-730,0,0,0,0.43,0
-731,0,0,0,0.23,0
-732,0,0,0,0.12,0
-733,0,0,0,0.07,0
-734,0,0,0,0.04,0
-735,0,0,0,0.02,0
-736,0,0,0,0.01,0
-737,0,0,0,0.01,0
-738,0,0,0,0,0
-739,0,0,0,0,0
-740,0,0,0,0,0
-741,0,0,0,0,0
-742,0,0,0,0,0
-743,0,0,0,0,0
-744,0,0,0,0,0
-745,0,0,0,0,0
-746,0,0,0,0,0
-747,0,0,0,0,0.01
-748,0,0,0,0,0.01
-749,0,0,0,0,0.01
-750,0,0,0,0,0.01
-751,0,0,0,0,0.02
-752,0,0,0,0,0.02
-753,0,0,0,0,0.03
-754,0,0,0,0,0.05
-755,0,0,0,0,0.07
-756,0,0,0,0,0.1
-757,0,0,0,0,0.15
-758,0,0,0,0,0.23
-759,0,0,0,0,0.35
-760,0,0,0,0,0.5
-761,0,0,0,0,0.67
-762,0,0,0,0,0.83
-763,0,0,0,0,0.94
-764,0,0,0,0,0.99
-765,0,0,0,0,1
-766,0,0,0,0,0.99
-767,0,0,0,0,0.98
-768,0,0,0,0,0.97
-769,0,0,0,0,0.97
-770,0,0,0,0,0.97
-771,0,0,0,0,0.97
-772,0,0,0,0,0.97
-773,0,0,0,0,0.98
-774,0,0,0,0,0.98
-775,0,0,0,0,0.98
-776,0,0,0,0,0.98
-777,0,0,0,0,0.98
-778,0,0,0,0,0.98
-779,0,0,0,0,0.97
-780,0,0,0,0,0.97
-781,0,0,0,0,0.97
-782,0,0,0,0,0.96
-783,0,0,0,0,0.96
-784,0,0,0,0,0.96
-785,0,0,0,0,0.96
-786,0,0,0,0,0.96
-787,0,0,0,0,0.96
-788,0,0,0,0,0.96
-789,0,0,0,0,0.96
-790,0,0,0,0,0.96
-791,0,0,0,0,0.96
-792,0,0,0,0,0.96
-793,0,0,0,0,0.96
-794,0,0,0,0,0.96
-795,0,0,0,0,0.96
-796,0,0,0,0,0.96
-797,0,0,0,0,0.95
-798,0,0,0,0,0.95
-799,0,0,0,0,0.95
-800,0,0,0,0,0.95
-801,0,0,0,0,0.95
-802,0,0,0,0,0.94
-803,0,0,0,0,0.94
-804,0,0,0,0,0.94
-805,0,0,0,0,0.94
-806,0,0,0,0,0.94
-807,0,0,0,0,0.93
-808,0,0,0,0,0.93
-809,0,0,0,0,0.93
-810,0,0,0,0,0.93
-811,0,0,0,0,0.93
-812,0,0,0,0,0.93
-813,0,0,0,0,0.93
-814,0,0,0,0,0.93
-815,0,0,0,0,0.93
-816,0,0,0,0,0.93
-817,0,0,0,0,0.93
-818,0,0,0,0,0.93
-819,0,0,0,0,0.93
-820,0,0,0,0,0.92
-821,0,0,0,0,0.92
-822,0,0,0,0,0.92
-823,0,0,0,0,0.92
-824,0,0,0,0,0.91
-825,0,0,0,0,0.91
-826,0,0,0,0,0.9
-827,0,0,0,0,0.9
-828,0,0,0,0,0.9
-829,0,0,0,0,0.89
-830,0,0,0,0,0.89
-831,0,0,0,0,0.89
-832,0,0,0,0,0.89
-833,0,0,0,0,0.88
-834,0,0,0,0,0.88
-835,0,0,0,0,0.88
-836,0,0,0,0,0.88
-837,0,0,0,0,0.88
-838,0,0,0,0,0.87
-839,0,0,0,0,0.87
-840,0,0,0,0,0.86
-841,0,0,0,0,0.85
-842,0,0,0,0,0.84
-843,0,0,0,0,0.84
-844,0,0,0,0,0.84
-845,0,0,0,0,0.84
-846,0,0,0,0,0.84
-847,0,0,0,0,0.83
-848,0,0,0,0,0.8
-849,0,0,0,0,0.71
-850,0,0,0,0,0.58
-851,0,0,0,0,0.43
-852,0,0,0,0,0.3
-853,0,0,0,0,0.19
-854,0,0,0,0,0.12
-855,0,0,0,0,0.08
-856,0,0,0,0,0.05
-857,0,0,0,0,0.03
-858,0,0,0,0,0.02
-859,0,0,0,0,0.02
-860,0,0,0,0,0.01
-861,0,0,0,0,0.01
-862,0,0,0,0,0.01
-863,0,0,0,0,0
-864,0,0,0,0,0
-865,0,0,0,0,0
-866,0,0,0,0,0
-867,0,0,0,0,0
-868,0,0,0,0,0
-869,0,0,0,0,0
-870,0,0,0,0,0
-871,0,0,0,0,0
-872,0,0,0,0,0
-873,0,0,0,0,0
-874,0,0,0,0,0
-875,0,0,0,0,0
-876,0,0,0,0,0
-877,0,0,0,0,0
-878,0,0,0,0,0
-879,0,0,0,0,0
-880,0,0,0,0,0
-881,0,0,0,0,0
-882,0,0,0,0,0
-883,0,0,0,0,0
-884,0,0,0,0,0
-885,0,0,0,0,0
-886,0,0,0,0,0
-887,0,0,0,0,0
-888,0,0,0,0,0
-889,0,0,0,0,0
-890,0,0,0,0,0
-891,0,0,0,0,0
-892,0,0,0,0,0
-893,0,0,0,0,0
-894,0,0,0,0,0
-895,0,0,0,0,0
-896,0,0,0,0,0
-897,0,0,0,0,0
-898,0,0,0,0,0
-899,0,0,0,0,0
-900,0,0,0,0,0
+420,0.002,,,,
+421,,,,,
+422,,,,,
+423,,,,,
+424,,0.001,,,
+425,,,,,
+426,,,,,
+427,,,,,
+428,0.003,,,,
+429,0.002,,,,
+430,,,,,
+431,,,,,
+432,,,,,
+433,,,,,
+434,,,,,
+435,0.001,,,,
+436,0.002,,,,
+437,0.004,,,,
+438,0.015,,,,
+439,0.065,,,,
+440,0.321,,,,
+441,0.687,,,,
+442,0.719,,,,
+443,0.731,,,,
+444,0.74,,,,
+445,0.74,,,,
+446,0.746,,,,
+447,0.756,,,,
+448,0.762,,,,
+449,0.765,,,,
+450,0.77,,,,
+451,0.775,,,,
+452,0.78,,,,
+453,0.782,,,,
+454,0.783,,,,
+455,0.784,,,,
+456,0.788,,,,
+457,0.792,,,,
+458,0.795,,,,
+459,0.796,,,,
+460,0.796,,0.001,,
+461,0.797,0.001,,,
+462,0.799,0.002,,,
+463,0.802,0.005,,,
+464,0.805,0.003,,,
+465,0.806,0.002,,,
+466,0.806,0.001,,,
+467,0.804,0.002,,,
+468,0.804,0.004,,,
+469,0.805,0.034,,,
+470,0.807,0.008,,,0.001
+471,0.812,0.001,,,
+472,0.816,,,,
+473,0.818,,,,
+474,0.82,,,,
+475,0.82,,,,
+476,0.822,,,,
+477,0.825,,,,
+478,0.83,,,,
+479,0.835,,,,
+480,0.84,,,,
+481,0.843,,,,
+482,0.845,,,,
+483,0.849,,,,
+484,0.855,,,,
+485,0.862,,,,
+486,0.869,,0.002,,
+487,0.875,,,,
+488,0.878,0.001,,,
+489,0.881,0.027,,,
+490,0.886,0.002,,,
+491,0.894,,,,
+492,0.905,,,,
+493,0.916,,,,
+494,0.924,,,,
+495,0.928,,,,
+496,0.931,,,,
+497,0.936,,,,
+498,0.947,,,,
+499,0.96,,,,
+500,0.969,,,0.027,
+501,0.97,,0.001,,
+502,0.967,,0.001,,
+503,0.97,,0.001,,
+504,0.984,,,,0.001
+505,1,,,,0.001
+506,0.995,,,,
+507,0.976,,,,
+508,0.999,,,,
+509,0.959,0.001,,,
+510,0.437,0.001,,,
+511,0.109,0.002,,,
+512,0.029,0.002,,,
+513,0.009,0.004,,,
+514,0.003,0.007,,,
+515,0.001,0.013,,,
+516,0.001,0.026,,,
+517,,0.054,,,
+518,,0.121,,,
+519,,0.272,,,0.002
+520,,0.539,,,0.011
+521,,0.784,,,0.001
+522,,0.868,,,
+523,,0.87,,,
+524,,0.864,,,
+525,,0.868,,,
+526,,0.873,,,
+527,,0.877,,,
+528,,0.876,,0.001,
+529,,0.873,,,
+530,,0.871,,,
+531,,0.872,,,
+532,,0.874,,,
+533,,0.877,,0.001,
+534,,0.881,,0.007,
+535,,0.882,,0.001,
+536,,0.883,,,0.001
+537,,0.882,,,0.001
+538,,0.881,,,0.001
+539,,0.88,,,
+540,,0.881,,,
+541,,0.883,,,
+542,,0.886,,,
+543,,0.891,0.001,,
+544,0.001,0.894,,0.001,
+545,,0.897,,,
+546,,0.898,,,0.002
+547,,0.899,,,0.002
+548,,0.898,,0.001,
+549,,0.898,,,
+550,,0.901,,,
+551,,0.905,,,
+552,,0.91,,,
+553,,0.915,,,
+554,,0.92,,,
+555,,0.924,,,
+556,,0.927,,,
+557,,0.928,,,
+558,,0.93,,,
+559,,0.933,,,
+560,,0.936,,,
+561,,0.942,,,
+562,,0.946,,,
+563,,0.949,,,
+564,,0.952,,,
+565,,0.953,,,
+566,,0.955,,,
+567,,0.96,,,
+568,,0.965,,,
+569,,0.97,,,
+570,,0.974,,,
+571,,0.977,,,
+572,,0.976,,,
+573,,0.975,,,0.001
+574,,0.974,,,0.001
+575,,0.976,,,
+576,,0.982,,,
+577,,0.989,0.001,,
+578,,0.993,,,
+579,,0.993,,,
+580,,0.988,,,
+581,,0.982,,,
+582,,0.984,,,
+583,,0.992,,,
+584,,1,,,
+585,,0.994,,,
+586,,0.976,,,
+587,,0.97,,,
+588,,0.999,,,
+589,,0.89,0.001,,
+590,,0.417,0.002,,
+591,,0.126,,,
+592,,0.039,,,
+593,,0.014,,,
+594,,0.006,,,
+595,,0.002,,,
+596,,0.001,,,
+597,,0.001,,,
+598,,,,,
+599,,,,,
+600,,,,,
+601,,,,,
+602,,,,,
+603,,,,,0.002
+604,,,,,0.004
+605,,,,,0.001
+606,,,,,
+607,,,,,
+608,,,,,
+609,,,,,
+610,,,,0.022,
+611,,,,,
+612,,,,,
+613,,,,,
+614,,,,,
+615,,,,,
+616,,,,,
+617,,,,,
+618,,,,,
+619,,,,,
+620,,,0.001,,
+621,,,0.001,,
+622,,,0.002,,
+623,,,0.003,,0.001
+624,,,0.005,,0.001
+625,,,0.009,,0.003
+626,,,0.018,,0.003
+627,,,0.038,,0.001
+628,,,0.086,,0.001
+629,,,0.2,,
+630,,,0.437,,
+631,,,0.723,0.001,
+632,,,0.856,0.001,
+633,,,0.863,,
+634,,,0.852,,
+635,,,0.854,,
+636,,,0.865,,
+637,,,0.876,,
+638,,,0.881,,
+639,,,0.882,,
+640,,,0.881,,
+641,,,0.882,,0.001
+642,,,0.885,,0.001
+643,,,0.89,,0.001
+644,,,0.896,0.016,0.002
+645,,,0.902,0.001,0.002
+646,,,0.907,,0.001
+647,,,0.909,,
+648,,,0.912,,
+649,,,0.913,,
+650,,,0.915,,
+651,,,0.918,,
+652,,,0.923,,
+653,,,0.928,,
+654,,,0.934,,
+655,,,0.939,,
+656,,,0.944,,
+657,,,0.947,,
+658,,,0.951,,
+659,,,0.954,,
+660,,,0.958,,
+661,,,0.961,,
+662,,,0.963,,
+663,,,0.966,,
+664,,,0.968,,
+665,,,0.97,,
+666,,,0.973,,
+667,,,0.976,,
+668,,,0.981,,
+669,,,0.985,,
+670,,,0.989,,
+671,,,0.991,,
+672,,,0.991,,0.001
+673,,,0.989,,0.002
+674,0.001,,0.985,,0.022
+675,,,0.985,,0.003
+676,,,0.989,,0.001
+677,,,0.994,0.001,
+678,,,1,0.001,
+679,,,0.998,0.001,
+680,,,0.989,0.002,
+681,,,0.981,0.003,
+682,,,0.989,0.004,
+683,,,0.999,0.007,
+684,,,0.858,0.012,
+685,,,0.463,0.021,
+686,,,0.174,0.039,
+687,,,0.062,0.074,
+688,,,0.024,0.143,
+689,,,0.01,0.275,
+690,,,0.005,0.491,
+691,,,0.002,0.74,
+692,,,0.001,0.914,
+693,,,0.001,0.984,
+694,,,,0.999,
+695,,,,0.998,
+696,,,,0.997,
+697,,,,0.999,
+698,,,,1,
+699,,,,1,
+700,,,,0.998,
+701,,,,0.995,
+702,,,,0.993,
+703,,,,0.99,
+704,,,,0.988,
+705,,,,0.987,
+706,,,,0.987,
+707,,,,0.986,
+708,,,,0.985,
+709,,,,0.984,
+710,,,,0.982,
+711,,,,0.979,
+712,,,,0.976,
+713,,,,0.972,
+714,,,,0.969,
+715,,,,0.966,
+716,,,,0.965,
+717,,,,0.964,
+718,,,,0.964,
+719,,,,0.963,
+720,,,,0.961,
+721,,,,0.956,
+722,,,,0.949,
+723,,,,0.942,
+724,,,,0.938,
+725,,,,0.939,
+726,,,,0.946,
+727,,,,0.936,
+728,,,,0.86,
+729,,,,0.671,
+730,,,,0.425,
+731,,,,0.234,
+732,,,,0.123,
+733,,,,0.065,
+734,,,,0.035,
+735,,,,0.02,
+736,,,,0.012,
+737,,,,0.007,
+738,0.001,,,0.004,0.001
+739,0.001,,,0.003,0.001
+740,,,,0.002,0.001
+741,,,,0.001,0.001
+742,,,,0.001,0.001
+743,,,,0.001,0.002
+744,,,,,0.002
+745,,,,,0.003
+746,,,,,0.004
+747,,,,,0.005
+748,,,,,0.006
+749,,,,,0.009
+750,,,,,0.012
+751,,,,,0.016
+752,,,,,0.023
+753,,,,,0.032
+754,,,,,0.046
+755,,,,,0.068
+756,,,,,0.102
+757,,,,,0.153
+758,,,,,0.231
+759,,,,,0.345
+760,,,,,0.497
+761,,,,,0.67
+762,,,,,0.828
+763,,,,,0.937
+764,,0.001,,,0.988
+765,,0.001,,,1
+766,,0.001,,,0.994
+767,,0.001,,,0.982
+768,,0.001,,,0.973
+769,,0.001,,,0.969
+770,,0.001,,,0.967
+771,,0.001,,,0.969
+772,,0.001,,,0.974
+773,,,,,0.978
+774,,,,,0.981
+775,,,,,0.983
+776,,,,,0.983
+777,,,,,0.981
+778,,,,,0.978
+779,,,,,0.974
+780,,,,,0.97
+781,,,,,0.967
+782,,,,,0.963
+783,,,,,0.961
+784,,,,,0.96
+785,,,,,0.958
+786,,,,,0.958
+787,,,,,0.957
+788,,,,,0.958
+789,,,,,0.959
+790,,,,,0.958
+791,,,,,0.959
+792,,,,,0.959
+793,,,,,0.958
+794,,,,,0.957
+795,,,,,0.956
+796,,,,,0.955
+797,,,,,0.954
+798,,,,,0.953
+799,,,,,0.95
+800,,,,,0.948
+801,,,,,0.946
+802,,,,,0.944
+803,,,,,0.941
+804,,,,,0.939
+805,,,,,0.937
+806,,,,,0.935
+807,,,,,0.933
+808,,,,,0.931
+809,,,,,0.929
+810,,,,,0.928
+811,,,,,0.927
+812,,,,,0.927
+813,,,,,0.927
+814,,,,,0.926
+815,,,,,0.926
+816,,,,,0.926
+817,,,,,0.926
+818,,,,,0.925
+819,,,,,0.925
+820,,,,,0.923
+821,,,,,0.921
+822,,,,,0.918
+823,,,,,0.915
+824,,,,,0.911
+825,,,,,0.906
+826,,,,,0.902
+827,,,,,0.898
+828,,,,,0.895
+829,,,,,0.891
+830,,,,,0.889
+831,,,,,0.886
+832,,,,,0.885
+833,,,,,0.884
+834,,,,,0.883
+835,,,,,0.882
+836,,,,,0.88
+837,,,,,0.876
+838,,,,,0.871
+839,,,,,0.865
+840,,,,,0.857
+841,,,,,0.849
+842,,,,,0.842
+843,,,,,0.838
+844,,,,,0.838
+845,,,,,0.84
+846,,,,,0.842
+847,,,,,0.832
+848,,,,,0.795
+849,,,,,0.712
+850,,,,,0.582
+851,,,,,0.431
+852,,,,,0.295
+853,,,,,0.193
+854,,,,,0.124
+855,,,,,0.08
+856,,,,,0.052
+857,,,,,0.034
+858,,,,,0.023
+859,,,,,0.016
+860,,,,,0.011
+861,,,,,0.008
+862,,,,,0.006
+863,,,,,0.004
+864,,,,,0.003
+865,,,,,0.002
+866,,,,,0.002
+867,,,,,0.001
+868,,,,,0.001
+869,,,,,0.001
+870,,,,,0.001
+871,,,,,
+872,,,,,
+873,,,,,
+874,,,,,
+875,,0.001,,,
+876,,0.001,,,
+877,,,,,
+878,,,,,
+879,,,,,
+880,,,,,
+881,,,,,
+882,,,,,
+883,,,,,
+884,,,,,
+885,,,,,
+886,,,,,
+887,,,,,
+888,,,,,
+889,,,,,
+890,,,,,
+891,,,,,
+892,,,,,
+893,,,,,
+894,,,,,
+895,,,,,
+896,,,,,
+897,,,,,
+898,,,,,
+899,,,,,
+900,,,,,