Added IKONOS atmospherical correction

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

imagery/i.atcorr/GeomCond.cpp

@@ -365,6 +365,7 @@ void GeomCond::parse()
     case 9: /* liss  ( IRS 1C)   * enter month,day,hh.ddd,long.,lat. */
     case 10: /* aster            * enter month,day,hh.ddd,long.,lat. */
     case 11: /* avnir            * enter month,day,hh.ddd,long.,lat. */
+    case 12: /* ikonos           * enter month,day,hh.ddd,long.,lat. */
     {
 	cin >> month;
 	cin >> jday;
@@ -410,7 +411,7 @@ void GeomCond::parse()
 /* ---- print geometrical conditions ---- */
 void GeomCond::print()
 {
-    static const string etiq1[12] = {
+    static const string etiq1[13] = {
 	string(" user defined conditions     "),
 	string(" meteosat observation        "),
 	string(" goes east observation       "),
@@ -422,7 +423,8 @@ void GeomCond::print()
 	string(" etm+     observation        "),
 	string(" liss     observation        "),
 	string(" aster    observation        "),
-	string(" avnir    observation        ")
+	string(" avnir    observation        "),
+	string(" ikonos   observation        ")
     };
 
     static const string head(" geometrical conditions identity  ");

imagery/i.atcorr/GeomCond.h

@@ -51,6 +51,8 @@
 /*                                                                      c */
 /*         11      avnir             * enter month,day,hh.ddd,long.,lat.c */
 /*                                                                      c */
+/*         12      ikonos            * 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

@@ -2040,6 +2040,148 @@ void IWave::avnir(int iwa)
     }
 }
 
+/* Following filter function created using create_iwave.py */
+
+void IWave::ikonos(int iwa)
+{
+
+    /* "Ikonos Green" of ikonos */
+    static const float sr1[274] = {
+		.0008, .0011, .0015, .0012, .0009, .0008, .0007,
+		.0007, .0010, .0013, .0015, .0017, .0018, .0019,
+		.0028, .0040, .0052, .0064, .0077, .0082, .0088,
+		.0119, .0147, .0175, .0262, .0350, .0602, .0854,
+		.1929, .2617, .3305, .3750, .4195, .4638, .5082,
+		.6132, .6437, .6742, .6973, .7204, .7406, .7608,
+		.7601, .7742, .7883, .8225, .8566, .8730, .8894,
+		.9436, .9718, 1.0000, .9934, .9868, .9562, .9255,
+		.7613, .6519, .5424, .4282, .3139, .2590, .2041,
+		.1519, .1188, .0856, .0644, .0432, .0371, .0310,
+		.0266, .0244, .0222, .0221, .0221, .0242, .0262,
+		.0296, .0283, .0270, .0248, .0226, .0240, .0254,
+		.0305, .0289, .0274, .0243, .0213, .0215, .0217,
+		.0296, .0319, .0342, .0306, .0271, .0233, .0195,
+		.0195, .0193, .0190, .0165, .0140, .0107, .0073,
+		.0060, .0056, .0051, .0065, .0079, .0101, .0123,
+		.0200, .0210, .0221, .0214, .0208, .0199, .0190,
+		.0155, .0126, .0097, .0083, .0068, .0059, .0049,
+		.0054, .0085, .0115, .0169, .0222, .0231, .0240,
+		.0121, .0080, .0040, .0031, .0022, .0024, .0026,
+		.0023, .0029, .0035, .0030, .0025, .0024, .0023,
+		.0040, .0038, .0036, .0032, .0027, .0025, .0024,
+		.0026, .0024, .0022, .0019, .0016, .0016, .0016,
+		.0020, .0020, .0020, .0023, .0025, .0031, .0036,
+		.0056, .0061, .0066, .0071, .0075, .0078, .0080,
+		.0095, .0089, .0083, .0085, .0088, .0090, .0093,
+		.0107, .0116, .0125, .0128, .0131, .0131, .0131,
+		.0119, .0112, .0105, .0100, .0095, .0091, .0088,
+		.0084, .0083, .0082, .0081, .0080, .0082, .0084,
+		.0081, .0081, .0080, .0080, .0080, .0074, .0069,
+		.0054, .0051, .0048, .0046, .0045, .0043, .0041,
+		.0036, .0035, .0035, .0034, .0033, .0030, .0028,
+		.0028, .0026, .0024, .0024, .0023, .0025, .0026,
+		.0026, .0026, .0026, .0021, .0016, .0018, .0019,
+		.0024
+    };
+
+    /* "Ikonos Red" of ikonos */
+    static const float sr2[274] = {
+		.0006, .0007, .0008, .0008, .0007, .0008, .0008,
+		.0006, .0007, .0008, .0007, .0006, .0008, .0009,
+		.0006, .0010, .0013, .0013, .0012, .0019, .0027,
+		.0010, .0011, .0013, .0016, .0019, .0023, .0027,
+		.0034, .0037, .0040, .0049, .0057, .0074, .0091,
+		.0113, .0113, .0113, .0122, .0130, .0151, .0172,
+		.0245, .0306, .0367, .0435, .0503, .0565, .0628,
+		.1094, .1512, .1931, .2421, .2911, .3653, .4396,
+		.6507, .7286, .8066, .8391, .8717, .8801, .8884,
+		.8898, .8824, .8749, .8635, .8521, .8722, .8923,
+		.9977, .9819, .9661, .9369, .9077, .9005, .8932,
+		.9459, .9715, .9972, .9785, .9598, .9186, .8774,
+		.7497, .6259, .5020, .3926, .2832, .2215, .1599,
+		.1068, .0947, .0825, .0732, .0640, .0557, .0474,
+		.0320, .0273, .0227, .0196, .0166, .0155, .0145,
+		.0148, .0155, .0162, .0174, .0186, .0193, .0199,
+		.0188, .0170, .0152, .0138, .0125, .0122, .0119,
+		.0109, .0122, .0135, .0150, .0165, .0174, .0183,
+		.0159, .0135, .0112, .0088, .0063, .0059, .0055,
+		.0034, .0036, .0039, .0037, .0035, .0032, .0030,
+		.0031, .0027, .0023, .0029, .0034, .0035, .0036,
+		.0053, .0059, .0065, .0068, .0072, .0080, .0089,
+		.0095, .0088, .0080, .0075, .0070, .0067, .0064,
+		.0061, .0059, .0056, .0054, .0053, .0057, .0061,
+		.0074, .0070, .0067, .0067, .0067, .0081, .0094,
+		.0128, .0135, .0143, .0151, .0159, .0165, .0171,
+		.0172, .0168, .0163, .0158, .0152, .0143, .0134,
+		.0120, .0113, .0106, .0103, .0100, .0096, .0092,
+		.0098, .0093, .0089, .0087, .0084, .0081, .0078,
+		.0075, .0072, .0070, .0064, .0059, .0056, .0054,
+		.0046, .0043, .0039, .0038, .0036, .0034, .0033,
+		.0029, .0029, .0028, .0028, .0028, .0025, .0022,
+		.0025, .0024, .0022, .0024, .0025, .0023, .0021,
+		.0017, .0019, .0021, .0021, .0020, .0021, .0022,
+		.0014
+    };
+
+    /* "Ikonos NIR" of ikonos */
+    static const float sr3[274] = {
+		.0000, .0003, .0005, .0004, .0004, .0005, .0006,
+		.0003, .0006, .0008, .0007, .0007, .0005, .0004,
+		.0015, .0018, .0020, .0025, .0030, .0033, .0035,
+		.0040, .0041, .0041, .0060, .0078, .0081, .0084,
+		.0082, .0066, .0051, .0038, .0026, .0035, .0043,
+		.0055, .0055, .0055, .0042, .0029, .0023, .0017,
+		.0015, .0026, .0036, .0040, .0043, .0034, .0025,
+		.0020, .0017, .0013, .0014, .0015, .0015, .0015,
+		.0021, .0022, .0024, .0025, .0027, .0030, .0034,
+		.0029, .0029, .0029, .0035, .0041, .0050, .0060,
+		.0072, .0083, .0093, .0093, .0093, .0090, .0087,
+		.0095, .0099, .0103, .0118, .0133, .0169, .0204,
+		.0285, .0316, .0347, .0388, .0428, .0502, .0575,
+		.0855, .1052, .1249, .1506, .1764, .2174, .2584,
+		.4188, .5285, .6381, .7213, .8045, .8414, .8784,
+		.9196, .9364, .9532, .9545, .9559, .9424, .9289,
+		.9389, .9628, .9866, .9917, .9968, .9524, .9079,
+		.7698, .6808, .5919, .5067, .4215, .3634, .3053,
+		.2262, .1923, .1583, .1315, .1046, .0860, .0674,
+		.0447, .0373, .0299, .0269, .0239, .0230, .0220,
+		.0184, .0171, .0159, .0161, .0162, .0159, .0157,
+		.0165, .0168, .0171, .0167, .0163, .0142, .0121,
+		.0091, .0079, .0068, .0065, .0062, .0059, .0056,
+		.0042, .0042, .0042, .0046, .0050, .0045, .0040,
+		.0040, .0033, .0027, .0028, .0030, .0028, .0025,
+		.0033, .0029, .0026, .0026, .0026, .0027, .0028,
+		.0027, .0028, .0028, .0027, .0026, .0026, .0026,
+		.0025, .0024, .0024, .0024, .0023, .0026, .0029,
+		.0027, .0029, .0031, .0035, .0039, .0040, .0040,
+		.0040, .0037, .0035, .0031, .0026, .0024, .0022,
+		.0016, .0015, .0014, .0014, .0013, .0013, .0014,
+		.0011, .0011, .0012, .0012, .0012, .0011, .0010,
+		.0016, .0013, .0010, .0012, .0013, .0014, .0014,
+		.0029, .0023, .0018, .0020, .0023, .0018, .0013,
+		.0009
+    };
+
+    static const float wli[3] = {0.350, 0.350, 0.350};
+    static const float wls[3] = {1.035, 1.035, 1.035};
+
+    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 < 274; i++)  ffu.s[40+i] = sr1[i];
+        break;
+    case 2: for(i = 0; i < 274; i++)  ffu.s[40+i] = sr2[i];
+        break;
+    case 3: for(i = 0; i < 274; i++)  ffu.s[40+i] = sr3[i];
+        break;
+    }
+}
+
 float IWave::equivwl() const
 {
     float seb = 0;
@@ -2108,6 +2250,7 @@ void IWave::parse()
 	else if(iwave <= 71)	irs_1c_liss(iwave - 67);
 	else if(iwave <= 80)	aster(iwave - 71);
 	else if(iwave <= 84)    avnir(iwave - 80);
+	else if(iwave <= 87)    ikonos(iwave - 84);
 	else G_warning(_("Unsupported iwave value: %d"), iwave);
     }
 
@@ -2124,7 +2267,7 @@ void IWave::parse()
 /* --- spectral condition ---- */
 void IWave::print()
 {
-    static const string nsat[86] = {
+    static const string nsat[89] = {
 	string(" constant        "), string(" user s          "),
 	string(" meteosat        "), string(" goes east       "), string(" goes west       "),
 	string(" avhrr 1 (noaa6) "), string(" avhrr 2 (noaa6) "),
@@ -2160,7 +2303,8 @@ void IWave::print()
 	string(" aster 4         "), string(" aster 5         "), string(" aster 6         "),
 	string(" aster 7         "), string(" aster 8         "), string(" aster 9         "),
 	string(" avnir 1         "), string(" avnir 2         "), string(" avnir 3         "),
-	string(" avnir 4         ")
+	string(" avnir 4         "), string(" ikonos green    "), string(" ikonos red      "),
+	string(" ikonos NIR      ")
     };
 
 

imagery/i.atcorr/Iwave.h

@@ -106,6 +106,9 @@ c        81  1st band of avnir        ( 0.390-0.550 )                  c
 c        82  2nd      "               ( 0.485-0.695 )                  c
 c        83  3rd      "               ( 0.545-0.745 )                  c
 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  note: wl has to be in micrometer                                    c
 c**********************************************************************/
 
@@ -143,6 +146,7 @@ private:
 	void irs_1c_liss(int iwa);
 	void aster(int iwa);
 	void avnir(int iwa);
+	void ikonos(int iwa);
 
 
 public:

imagery/i.atcorr/i.atcorr.html

@@ -146,6 +146,12 @@ the ascendant node at equator</td>
 <td>enter month,day,hh.ddd,long.,lat. *</td>
 </tr>
 
+<tr>
+<td>11</td>
+<td><b>ikonos</b></td>
+<td>enter month,day,hh.ddd,long.,lat. *</td>
+</tr>
+
 </table>
 
 <blockquote>
@@ -518,6 +524,10 @@ micrometer.</td>
 <tr><td>83</td><td>avnir band 3 (0.545-0.745)</td></tr>
 <tr><td>84</td><td>avnir band 4 (0.700-0.925)</td></tr>
 
+<tr><td>85</td><td><b>ikonos</b> Green band (0.350-1.035)</td></tr>
+<tr><td>86</td><td>ikonos Red band (0.350-1.035)</td></tr>
+<tr><td>87</td><td>ikonos NIR band (0.350-1.035)</td></tr>
+
 </table>
 
 <h2>EXAMPLES</h2>