[GRASS-SVN] r70602 - in grass-addons/grass7/raster/r.modis: . r.modis.download r.modis.import
svn_grass at osgeo.org
svn_grass at osgeo.org
Thu Feb 16 10:00:10 PST 2017
Author: veroandreo
Date: 2017-02-16 10:00:10 -0800 (Thu, 16 Feb 2017)
New Revision: 70602
Modified:
grass-addons/grass7/raster/r.modis/r.modis.download/r.modis.download.html
grass-addons/grass7/raster/r.modis/r.modis.html
grass-addons/grass7/raster/r.modis/r.modis.import/r.modis.import.html
Log:
r.modis.* manuals: fix typos and add examples
Modified: grass-addons/grass7/raster/r.modis/r.modis.download/r.modis.download.html
===================================================================
--- grass-addons/grass7/raster/r.modis/r.modis.download/r.modis.download.html 2017-02-16 12:51:15 UTC (rev 70601)
+++ grass-addons/grass7/raster/r.modis/r.modis.download/r.modis.download.html 2017-02-16 18:00:10 UTC (rev 70602)
@@ -1,7 +1,7 @@
<h2>DESCRIPTION</h2>
<em>r.modis.download</em> downloads selected MODIS products. The module
-can download several tiles at once and also multiple observation days
+can download several tiles at once and also multiple observation dates
for each MODIS product.
<h2>NOTES</h2>
@@ -14,7 +14,7 @@
<ul>
<li> First time user: The user has to register at
<a href="https://urs.earthdata.nasa.gov/users/new">https://urs.earthdata.nasa.gov/users/new</a>;
- then to login and change to his/her profile page. Once there, under the "My applications"
+ then login and change to his/her profile page. Once there, under the "My applications"
tab, the user needs to approve the following applications: "LP DAAC Data Pool"
and "Earthdata Search".</li>
<li> If the user is already registered, he/she just needs to login and enable the
@@ -25,7 +25,7 @@
<p>
In order to download the desired MODIS product(s), the username and password must be
provided through the <em>settings</em> option. There are two ways: either pass a file
-in which the first row is the username, and the second row is the password, like as below:
+in which the first row is the username, and the second row is the password, as showed below:
<div class="code"><pre>
your_NASA_username
your_NASA_password
@@ -35,8 +35,8 @@
<b>Warning</b>: As per NASA policy no more than ten simultaneous connections
are permitted in parallel with their server. Hence, no more than ten processes
-of <em>r.modis.download</em> should run at the same time to avoid that your
-actual IP address will get blacklisted.
+of <em>r.modis.download</em> should run at the same time to avoid that user's
+actual IP address gets blacklisted.
<p>
By default, the downloaded files are stored in the path in which the settings
file was saved. The user can change this directory with the <em>folder</em>
@@ -60,7 +60,7 @@
</pre></div>
<p>
-Reading the user and password options from a file (this will download per default
+Reading the user and password options from a file (this will download by default
the "lst_terra_daily_1000" product). MODIS data will be downloaded to
the folder where the SETTING file is:
Modified: grass-addons/grass7/raster/r.modis/r.modis.html
===================================================================
--- grass-addons/grass7/raster/r.modis/r.modis.html 2017-02-16 12:51:15 UTC (rev 70601)
+++ grass-addons/grass7/raster/r.modis/r.modis.html 2017-02-16 18:00:10 UTC (rev 70602)
@@ -4,18 +4,18 @@
It uses the <a href="http://www.pymodis.org">pyModis</a> library and either the
<a href="http://gdal.org/">GDAL library</a> or the
<a href="https://lpdaac.usgs.gov/tools/modis_reprojection_tool"> MODIS Reprojection Tool</a> (MRT)
-software to convert, mosaik and process MODIS data.
+software to convert, mosaic and process MODIS data.
<p>
The <em>r.modis</em> suite requires the pyModis library and either the GDAL or MRT
software to be installed.
<p>
-The <em>r.modis</em> suite offers three modules as interface to MODIS data. Each modules
-is dedicated for a specific operation. The module <em>r.modis.download</em>
+The <em>r.modis</em> suite offers two modules as interface to MODIS data. Each module
+is dedicated to a specific operation. The module <em>r.modis.download</em>
is used to download MODIS HDF products from NASA servers. These files can then
be imported with <em>r.modis.import</em> which supports import of Level 3
-MODIS products as a single image or as a mosaik into GRASS GIS.
+MODIS products as a single image or as a mosaic into GRASS GIS.
<p>
-Subsequently, the user can create a temporal dataset using <em>t.create</em> and register
+Subsequently, the user can create a temporal dataset using <em>t.create</em> and, register
the maps with <em>t.register</em>. NOTE: In order to work with the temporal framework of
GRASS GIS the flag <em>w</em> must be set during the import with <em>r.modis.import</em>.
<p>
@@ -24,7 +24,8 @@
<!--
The suite process Level 2 data using the <em>r.modis.process</em> module.
-->
-The suite imports Level 3 (georeferenced) products either as single image or as daily mosaik.
+The suite imports Level 3 (georeferenced) products either as single images or as mosaics for
+each date.
<h2>Supported MODIS products</h2>
These products are currently supported:
@@ -38,7 +39,7 @@
(related <a href="https://lpdaac.usgs.gov/dataset_discovery/modis/modis_products_table/mod11a1_v006">MOD11A1</a> and
<a href="https://lpdaac.usgs.gov/dataset_discovery/modis/modis_products_table/myd11a1_v006">MYD11A1</a>
product pages)</li>
- <li><b>Land Surface Temperature daily 1 Km (Terra/Aqua)</b>: data are composed from the daily
+ <li><b>Land Surface Temperature eight day 1 Km (Terra/Aqua)</b>: data are composed from the daily
1-kilometer LST product (MOD11A1/MYD11A1) and stored on a 1-km Sinusoidal grid as the
average values of clear-sky LSTs during an 8-day period.<br>
MOD11A2/MYD11A2 is comprised of daytime and nighttime LSTs, quality assessment,
@@ -47,7 +48,7 @@
(related <a href="https://lpdaac.usgs.gov/dataset_discovery/modis/modis_products_table/mod11a2_v006">MOD11A2</a> and
<a href="https://lpdaac.usgs.gov/dataset_discovery/modis/modis_products_table/myd11a2_v006">MYD11A2</a>
product pages)</li>
- <li><b>Land Surface Temperature eight day 6 Km (Terra/Aqua)</b>: products provide per-pixel
+ <li><b>Land Surface Temperature eight day ~6 Km (Terra/Aqua)</b>: products provide per-pixel
temperature and emissivity values in a sequence of swath-based to grid-based
global products. The MODIS/Terra-Aqua LST/E Daily L3 Global 6 Km Grid
(Short name: MOD11B1/MYD11B1), is tile-based and gridded in the Sinusoidal projection,
Modified: grass-addons/grass7/raster/r.modis/r.modis.import/r.modis.import.html
===================================================================
--- grass-addons/grass7/raster/r.modis/r.modis.import/r.modis.import.html 2017-02-16 12:51:15 UTC (rev 70601)
+++ grass-addons/grass7/raster/r.modis/r.modis.import/r.modis.import.html 2017-02-16 18:00:10 UTC (rev 70602)
@@ -8,7 +8,7 @@
library. Please install it beforehand.
<p>
-The input file is given as list of full paths to the MODIS HDF files,
+The input file is given as a list of full paths to the MODIS HDF files,
one per line. The input file(s) have to be inside the folder where
the HDF files are stored.
<p>
@@ -17,7 +17,7 @@
<p>
The <em>mrtpath</em> option is the path to the main folder of the
MODIS Reprojection Tools (MRT) binaries, i.e. the folder which contains
-the bin/ and the data/ folder since these two folders are essential for
+the bin/ and the data/ folder, since these two folders are essential for
obtaining a successful result.
<p>
@@ -105,13 +105,82 @@
g.gui.timeline ndvi_16_5600m
</pre></div>
+<h3>Example of a complete workflow</h3>
+Download the data: MOD11A1 from 2016-12-23 to 2016-12-31, tiles h18v04,h18v05
+
+<div class="code"><pre>
+r.modis.download settings=$HOME/SETTING product=lst_terra_eight_1000 \
+ tiles=h18v04,h18v05 startday=2016-12-23 endday=2016-12-31
+</pre></div>
+
+Import mosaics of LST Day and QC Day bands for all dates in the list of files:
+
+<div class="code"><pre>
+r.modis.import -m files=$HOME/listfileMOD11A1.006.txt \
+ spectral="( 1 1 0 0 0 0 0 0 0 0 0 0 )"
+</pre></div>
+
+Extract and apply the mandatory QA band (for more details see
+<a href="i.modis.qc.html">i.modis.qc</a>):
+
+<div class="code"><pre>
+for map in `g.list type=raster pattern=*_QC_Day` ; do
+ i.modis.qc input=${map} output=${map}_mandatory_qa \
+ productname=mod11A1 qcname=mandatory_qa_11A1
+done
+
+for m in `g.list rast pat=*2016*LST_Day_1km` ; do
+ # get name of product and date from filenames
+ i=`echo $m | cut -c 1-16`
+ # apply qa flags
+ r.mapcalc --o expression="${m} = if(${i}_mosaic_QC_Day_mandatory_qa < 2, ${m}, null())"
+done
+</pre></div>
+
+Create the time series (i.e.: temporal database) and register maps:
+
+<div class="code"><pre>
+t.create type=strds temporaltype=absolute output=LST_Day_daily \
+ title="Daily LST Day 1km" \
+ description="Daily LST Day 1km MOD11A1.006, December 2016"
+t.register -i input=LST_Day_daily \
+ maps=`g.list type=raster pattern=*2016*LST_Day_1km separator=,` \
+ start="2016-12-21" increment="1 day" --o
+</pre></div>
+
+Verify list of maps and dates and visualize timeline:
+
+<div class="code"><pre>
+t.rast.list LST_Day_daily
+g.gui.timeline LST_Day_daily
+</pre></div>
+
+<p>
+It is also possible to create a time series using the list of maps with start and
+end time written by <em>r.modis.import</em> with the <em>w</em> flag.
+
+<div class="code"><pre>
+# Import mosaics of LST Day for all dates using the list of downloaded files from
+# r.modis.download (see above) and get a list of the imported files along with
+# dates to use with t.register in the temporal framework
+r.modis.import -mw files=$HOME/listfileMOD11A1.006.txt \
+ spectral="( 1 0 0 0 0 0 0 0 0 0 0 0 )"
+
+# Create time series and register maps
+t.create type=strds temporaltype=absolute output=LST_Day_daily \
+ title="Daily LST Day 1km" \
+ description="Daily LST Day 1km MOD11A1.006, December 2016"
+t.register input=LST_Day_daily file=$HOME/tmp/grass7-user-5370/tmp_rGPcg
+</pre></div>
+
<h2>SEE ALSO</h2>
<em>
<a href="r.modis.html">r.modis</a>,
<a href="r.modis.download.html">r.modis.download</a>,
<a href="i.modis.qc.html">i.modis.qc</a>
+<a href="i.modis.qc.html">i.modis.qc</a>
<!-- <a href="r.modis.process.html">r.modis.process</a>,-->
</em>
<p>
More information about the grass-commit
mailing list