|
@@ -2,9 +2,39 @@
|
|
|
|
|
|
<em>i.eb.soilheatflux</em> calculates the soil heat flux approximation (g0)
|
|
|
after Bastiaanssen (1995). The main reference for implementation is Alexandridis, 2009.
|
|
|
+
|
|
|
It takes input of Albedo, NDVI, Surface Skin temperature, Net Radiation (see
|
|
|
-<em>r.sun</em>), time of satellite overpass, and a flag for the Roerink empirical
|
|
|
-modification from the HAPEX-Sahel experiment.
|
|
|
+<em>r.sun</em>), time of satellite overpass, and a flag for the Roerink empirical modification from the HAPEX-Sahel experiment.
|
|
|
+
|
|
|
+The "time of satellite overpass" map can be obtained as follows:
|
|
|
+<ul>
|
|
|
+<li>MODIS: a related sub dataset is included in each HDF file, and simply
|
|
|
+ to be imported as a raster map;</li>
|
|
|
+<li>Landsat: to be generated as map from the overpass time stored in
|
|
|
+ the metadata file (given in Greenwich Mean Time - GMT), see below.</li>
|
|
|
+</ul>
|
|
|
+
|
|
|
+For Landsat, the overpass map can be computed by using a two-step method:
|
|
|
+
|
|
|
+<div class="code"><pre>
|
|
|
+# 1) extract the overpass time in GMT from metadata file
|
|
|
+
|
|
|
+i.landsat.toar -p input=dummy output=dummy2 \
|
|
|
+ metfile=LC81250452013338LGN00_MTL.txt lsatmet=time
|
|
|
+# ... in this example approx. 03:12am GMT
|
|
|
+
|
|
|
+# 2) create map for computational region of Landsat scene
|
|
|
+g.region rast=LC81250452013338LGN00_B4 -p
|
|
|
+r.mapcalc "overpasstime = 3.211328"
|
|
|
+</pre></div>
|
|
|
+
|
|
|
+<!--
|
|
|
+<h2>EXAMPLE</h2>
|
|
|
+
|
|
|
+<div class="code"><pre>
|
|
|
+i.eb.soilheatflux ...
|
|
|
+</pre></div>
|
|
|
+-->
|
|
|
|
|
|
<h2>SEE ALSO</h2>
|
|
|
|
|
@@ -14,6 +44,7 @@ modification from the HAPEX-Sahel experiment.
|
|
|
<a href="i.emissivity.html">i.emissivity</a>,
|
|
|
<a href="i.eb.hsebal01.html">i.eb.hsebal01</a>,
|
|
|
<a href="i.eb.evapfr.html">i.eb.evapfr</a>
|
|
|
+<a href="i.landsat.toar.html">i.landsat.toar</a>
|
|
|
</em>
|
|
|
|
|
|
<h2>REFERENCES</h2>
|