Newer
Older
<!DOCTYPE html>
<html class="writer-html5" lang="en" >
<head>
<meta charset="utf-8" /><meta name="generator" content="Docutils 0.18.1: http://docutils.sourceforge.net/" />
<meta name="viewport" content="width=device-width, initial-scale=1.0" />
<title>Preparing the land cover raster — modspa-pixel 0.1 documentation</title>
<link rel="stylesheet" href="../_static/pygments.css" type="text/css" />
<link rel="stylesheet" href="../_static/css/theme.css" type="text/css" />
<link rel="stylesheet" href="../_static/graphviz.css" type="text/css" />
<link rel="stylesheet" href="../_static/custom.css" type="text/css" />
<!--[if lt IE 9]>
<script src="../_static/js/html5shiv.min.js"></script>
<![endif]-->
<script src="../_static/jquery.js?v=5d32c60e"></script>
<script src="../_static/_sphinx_javascript_frameworks_compat.js?v=2cd50e6c"></script>
<script src="../_static/documentation_options.js?v=2709fde1"></script>
<script src="../_static/doctools.js?v=888ff710"></script>
<script src="../_static/sphinx_highlight.js?v=dc90522c"></script>
<script src="../_static/js/theme.js"></script>
<link rel="index" title="Index" href="../genindex.html" />
<link rel="search" title="Search" href="../search.html" />
<link rel="next" title="Preparing the soil data raster" href="soil_data.html" />
<link rel="prev" title="Preparing the weather data cube" href="weather_cube.html" />
<link href="../.static/custom.css" rel="stylesheet" type="text/css">
</head>
<body class="wy-body-for-nav">
<div class="wy-grid-for-nav">
<nav data-toggle="wy-nav-shift" class="wy-nav-side">
<div class="wy-side-scroll">
<div class="wy-side-nav-search" >
<a href="../index.html" class="icon icon-home">
modspa-pixel
</a>
<div role="search">
<form id="rtd-search-form" class="wy-form" action="../search.html" method="get">
<input type="text" name="q" placeholder="Search docs" aria-label="Search docs" />
<input type="hidden" name="check_keywords" value="yes" />
<input type="hidden" name="area" value="default" />
</form>
</div>
</div><div class="wy-menu wy-menu-vertical" data-spy="affix" role="navigation" aria-label="Navigation menu">
<p class="caption" role="heading"><span class="caption-text">Contents:</span></p>
<ul class="current">
<li class="toctree-l1"><a class="reference internal" href="../setup/setup.html">Setting up the correct environnment and parameters</a></li>
<li class="toctree-l1 current"><a class="reference internal" href="input_preparation.html">Preparing the inputs for Modspa-Pixel</a><ul class="current">
<li class="toctree-l2"><a class="reference internal" href="ndvi_cube.html">Preparing the NDVI data cube</a></li>
<li class="toctree-l2"><a class="reference internal" href="weather_cube.html">Preparing the weather data cube</a></li>
<li class="toctree-l2 current"><a class="current reference internal" href="#">Preparing the land cover raster</a><ul>
<li class="toctree-l3"><a class="reference internal" href="#parcel-mode">Parcel mode</a></li>
<li class="toctree-l3"><a class="reference internal" href="#pixel-mode">Pixel mode</a></li>
<li class="toctree-l3"><a class="reference internal" href="#land-cover-preprocessing">Land Cover preprocessing</a></li>
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
</ul>
</li>
<li class="toctree-l2"><a class="reference internal" href="soil_data.html">Preparing the soil data raster</a></li>
<li class="toctree-l2"><a class="reference internal" href="samir_params.html">Preparing the SAMIR parameter csv file</a></li>
</ul>
</li>
<li class="toctree-l1"><a class="reference internal" href="../samir/run_samir.html">Run the SAMIR model</a></li>
<li class="toctree-l1"><a class="reference internal" href="../file_structure.html">Project file structure</a></li>
</ul>
</div>
</div>
</nav>
<section data-toggle="wy-nav-shift" class="wy-nav-content-wrap"><nav class="wy-nav-top" aria-label="Mobile navigation menu" >
<i data-toggle="wy-nav-top" class="fa fa-bars"></i>
<a href="../index.html">modspa-pixel</a>
</nav>
<div class="wy-nav-content">
<div class="rst-content">
<div role="navigation" aria-label="Page navigation">
<ul class="wy-breadcrumbs">
<li><a href="../index.html" class="icon icon-home" aria-label="Home"></a></li>
<li class="breadcrumb-item"><a href="input_preparation.html">Preparing the inputs for Modspa-Pixel</a></li>
<li class="breadcrumb-item active">Preparing the land cover raster</li>
<li class="wy-breadcrumbs-aside">
<a href="../_sources/preprocessing/land_cover.rst.txt" rel="nofollow"> View page source</a>
</li>
</ul>
<hr/>
</div>
<div role="main" class="document" itemscope="itemscope" itemtype="http://schema.org/Article">
<div itemprop="articleBody">
<section id="preparing-the-land-cover-raster">
<span id="land-cover"></span><h1>Preparing the land cover raster<a class="headerlink" href="#preparing-the-land-cover-raster" title="Link to this heading"></a></h1>
<p>There is currently no automatic script to generate land cover rasters. This part has to be done manually by the user.</p>
<section id="parcel-mode">
<h2>Parcel mode<a class="headerlink" href="#parcel-mode" title="Link to this heading"></a></h2>
<p>The land cover information has to be included in the parcel definition shapefile. Each parcel should have a land cover class identification number (integers starting by one). Extraction is then done automatically by the NDVI extraction scripts. This will later be used to build an artificial raster to run the models.</p>
<div class="admonition warning">
<p class="admonition-title">Warning</p>
<p>Development in progress</p>
</div>
</section>
<section id="pixel-mode">
<h2>Pixel mode<a class="headerlink" href="#pixel-mode" title="Link to this heading"></a></h2>
<p>For the pixel mode, the user has to find a land cover raster and project it on the same grid as the <a class="reference internal" href="ndvi_cube.html#ndvi"><span class="std std-ref">NDVI</span></a> and <a class="reference internal" href="weather_cube.html#weather"><span class="std std-ref">weather</span></a> datasets. Each class should be represented as an integer (starting at one) and saved as a Geotiff file (for <code class="docutils literal notranslate"><span class="pre">OTB</span></code>). Once you have a valid raster, you can reproject it and clip it using the <a class="reference external" href="https://www.orfeo-toolbox.org/CookBook/Applications/app_Superimpose.html">OTB Superimpose</a> function. Here is an example of the command you can run to do that:</p>
<div class="highlight-bash notranslate"><div class="highlight"><pre><span></span><span class="o">(</span>modspa_pixel<span class="o">)</span><span class="w"> </span>$<span class="w"> </span>otbcli_Superimpose<span class="w"> </span>-inr<span class="w"> </span>reference_satellite_image.tif<span class="w"> </span>-inm<span class="w"> </span>lc_raster_input.tif<span class="w"> </span>-out<span class="w"> </span>lc_raster_output.tif<span class="w"> </span>-interpolator<span class="w"> </span>nn
</pre></div>
</div>
<p>Where:</p>
<ul class="simple">
<li><p><code class="docutils literal notranslate"><span class="pre">-inr</span></code> is the reference image on which you want to superimpose your raster (e.g. Sentinel-2 image or extraction of Sentinel-2 image).</p></li>
<li><p><code class="docutils literal notranslate"><span class="pre">-inm</span></code> is the raster you want to superimpose.</p></li>
<li><p><code class="docutils literal notranslate"><span class="pre">-out</span></code> is the name of the output raster that will be created</p></li>
<li><p><code class="docutils literal notranslate"><span class="pre">-interpolator</span></code> is the chosen interpolator (nn stands for nearest neighbor), other interpolators will not work as they will create new float values between the integer values.</p></li>
<li><p><em>additional parameters like the data type can be found in the OTB documentation</em></p></li>
</ul>
<p>You will then have a land cover raster on the same grid as the other inputs.</p>
</section>
<section id="land-cover-preprocessing">
<h2>Land Cover preprocessing<a class="headerlink" href="#land-cover-preprocessing" title="Link to this heading"></a></h2>
<p>Additional processing can be done with the <code class="docutils literal notranslate"><span class="pre">custom_inputs.ipynb</span></code> <strong>notebook</strong> found in the <code class="docutils literal notranslate"><span class="pre">preprocessing</span></code> directory.</p>
</section>
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
</section>
</div>
</div>
<footer><div class="rst-footer-buttons" role="navigation" aria-label="Footer">
<a href="weather_cube.html" class="btn btn-neutral float-left" title="Preparing the weather data cube" accesskey="p" rel="prev"><span class="fa fa-arrow-circle-left" aria-hidden="true"></span> Previous</a>
<a href="soil_data.html" class="btn btn-neutral float-right" title="Preparing the soil data raster" accesskey="n" rel="next">Next <span class="fa fa-arrow-circle-right" aria-hidden="true"></span></a>
</div>
<hr/>
<div role="contentinfo">
<p>© Copyright 2023, auclairj.</p>
</div>
Built with <a href="https://www.sphinx-doc.org/">Sphinx</a> using a
<a href="https://github.com/readthedocs/sphinx_rtd_theme">theme</a>
provided by <a href="https://readthedocs.org">Read the Docs</a>.
</footer>
</div>
</div>
</section>
</div>
<script>
jQuery(function () {
SphinxRtdTheme.Navigation.enable(true);
});
</script>
</body>
</html>