浏览代码

i.atcorr: Sentinel-2A support added using the Sentinel-2A Spectral Response Functions (S2A-SRF) from https://sentinel.esa.int/web/sentinel/user-guides/sentinel-2-msi/document-library

git-svn-id: https://svn.osgeo.org/grass/grass/trunk@70161 15284696-431f-4ddb-bdfa-cd5b030d7da7
Markus Neteler 8 年之前
父节点
当前提交
d2c96f9f5d

+ 3 - 1
imagery/i.atcorr/geomcond.cpp

@@ -378,6 +378,7 @@ void GeomCond::parse()
     case 22: /* pleiades1a       * enter month,day,hh.ddd,long.,lat. */
     case 23: /* pleiades1b       * enter month,day,hh.ddd,long.,lat. */
     case 24: /* worldview3       * enter month,day,hh.ddd,long.,lat. */
+    case 25: /* sentinel2a       * enter month,day,hh.ddd,long.,lat. */
    {
 	cin >> month;
 	cin >> jday;
@@ -448,7 +449,8 @@ void GeomCond::print()
 	string(" spot7 observation           "),
 	string(" pleiades1a observation      "),
 	string(" pleiades1b observation      "),
-	string(" worldview3 observation      ")
+	string(" worldview3 observation      "),
+	string(" sentinel2a observation      ")
 	};
 
     static const string head(" geometrical conditions identity  ");

+ 2 - 0
imagery/i.atcorr/geomcond.h

@@ -77,6 +77,8 @@
 /*                                                                        */
 /*         24      worldview3       * enter month,day,hh.ddd,long.,lat.   */
 /*                                                                        */
+/*         25      sentinel2a       * enter month,day,hh.ddd,long.,lat.   */
+/*                                                                        */
 /*     note:       for hrv and tm experiments long. and lat. are the      */
 /*                 coordinates of the scene center.                       */
 /*                 lat. must be > 0 for north lat., < 0 for south lat.    */

+ 25 - 1
imagery/i.atcorr/i.atcorr.html

@@ -219,6 +219,12 @@ the ascendant node at equator</td>
 <td>enter month,day,hh.ddd,long.,lat. *</td>
 </tr>
 
+<tr>
+<td>25</td>
+<td><b>Sentinel-2A</b></td>
+<td>enter month,day,hh.ddd,long.,lat. *</td>
+</tr>
+
 </table>
 
 <p>
@@ -692,6 +698,21 @@ micrometer.</td>
 <tr><td>163</td><td>Worldview3 SWIR7 band (2235nm - 2285 nm)</td></tr>
 <tr><td>165</td><td>Worldview3 SWIR8 band (2295nm - 2365 nm)</td></tr>
 
+<!-- derived from sensors_csv/sentinel_2A_msi.csv -->
+<tr><td>166</td><td><b>Sentinel2A</b> Coastal blue band B1 ( 430nm - 457nm)</td></tr>
+<tr><td>167</td><td>Sentinel2A Blue band B2 ( 440nm - 535nm)</td></tr>
+<tr><td>168</td><td>Sentinel2A Green band B3 ( 537nm - 582nm)</td></tr>
+<tr><td>169</td><td>Sentinel2A Red band B4 ( 646nm - 684nm)</td></tr>
+<tr><td>170</td><td>Sentinel2A Red edge band B5 ( 694nm - 713nm)</td></tr>
+<tr><td>171</td><td>Sentinel2A Red edge band B6 ( 731nm - 749nm)</td></tr>
+<tr><td>172</td><td>Sentinel2A Red edge band B7 ( 769nm - 797nm)</td></tr>
+<tr><td>173</td><td>Sentinel2A NIR band B8 ( 773nm -  908nm)</td></tr>
+<tr><td>174</td><td>Sentinel2A Red edge band B8A ( 848nm - 881nm)</td></tr>
+<tr><td>175</td><td>Sentinel2A Water vapour band B9 ( 932nm - 958nm)</td></tr>
+<tr><td>176</td><td>Sentinel2A SWIR Cirrus band B10 ( 1337nm - 1412nm)</td></tr>
+<tr><td>177</td><td>Sentinel2A SWIR band B11 ( 1539nm - 1682nm)</td></tr>
+<tr><td>178</td><td>Sentinel2A SWIR band B12 ( 2078nm - 2320nm)</td></tr>
+
 </table>
 
 <h2>EXAMPLES</h2>
@@ -889,7 +910,10 @@ overview., IEEE Trans. Geosc. and Remote Sens. 35(3):675-686.
 <br>Marco Vizzari
 
 <p><em>Worldview3 addition 8/2016:</em>
-<br>Markus Neteler, mundialis, Germany
+<br>Markus Neteler, mundialis.de, Germany
+
+<p><em>Sentinel-2A addition 12/2016:</em>
+<br>Markus Neteler, mundialis.de, Germany
 
 <p>
 <i>Last changed: $Date$</i>

文件差异内容过多而无法显示
+ 1618 - 8
imagery/i.atcorr/iwave.cpp


+ 1 - 0
imagery/i.atcorr/iwave.h

@@ -83,6 +83,7 @@ private:
 	void pleiades1a(int iwa);
 	void pleiades1b(int iwa);
 	void worldview3(int iwa);
+	void sentinel2a(int iwa);
 
 public:
 	/* To compute the equivalent wavelength needed for the calculation of the

文件差异内容过多而无法显示
+ 2302 - 0
imagery/i.atcorr/sensors_csv/sentinel_2A_msi.csv