<aclass="sourceLine"id="cb1-7"data-line-number="7"><spanclass="co">#> Geospatial Data Abstraction Library extensions to R successfully loaded</span></a>
<aclass="sourceLine"id="cb1-8"data-line-number="8"><spanclass="co">#> Loaded GDAL runtime: GDAL 2.2.2, released 2017/09/15</span></a>
<aclass="sourceLine"id="cb1-9"data-line-number="9"><spanclass="co">#> Path to GDAL shared files: /usr/share/gdal/2.2</span></a>
<aclass="sourceLine"id="cb1-10"data-line-number="10"><spanclass="co">#> GDAL binary built with GEOS: TRUE </span></a>
<aclass="sourceLine"id="cb1-16"data-line-number="16"><spanclass="kw"><ahref="https://www.rdocumentation.org/packages/base/topics/cat">cat</a></span>(<spanclass="st">"We suggest installing the pacakge mapview for interactive visualisation"</span>,</a>
<aclass="sourceLine"id="cb1-17"data-line-number="17"><spanclass="st">"of cartography from within R"</span>)</a></code></pre></div>
<spanclass="kw"><ahref="https://www.rdocumentation.org/packages/base/topics/cat">cat</a></span>(<spanclass="st">"We suggest installing the pacakge mapview for interactive visualisation"</span>,
<spanclass="st">"of cartography from within R"</span>)</code></pre></div>
<p>Download cartography from the Global Administrative Borders Database (GADM, <ahref="https://gadm.org/"class="uri">https://gadm.org/</a>) directly from within R.</p>
<divclass="sourceCode"id="cb5"><preclass="sourceCode r"><codeclass="sourceCode r"><aclass="sourceLine"id="cb5-1"data-line-number="1"><spanclass="co"># Not using this for the moment</span></a>
<ahref="#save-pre-processed-cartography-for-use-within-the-package"class="anchor"></a>Save pre-processed cartography for use within the package</h2>
<p>Prefer standard and modern Open Geospatial Consortium (<ahref="http://www.opengeospatial.org/">OGC</a>) formats: GeoPackage for vector maps and GeoTiff for raster images.</p>
<p>L’objet <code>cmr</code> est une simple liste de cartes : objets de type <code>Spatial*</code> pour les cartes vectoriels et de type <code>RasterLayer</code> pour les cartes raster.</p>
<p>Pour utiliser d’autre cartographie, utiliser les fonctions <code><ahref="https://www.rdocumentation.org/packages/rgdal/topics/readOGR">rgdal::readOGR()</a></code> et <code><ahref="https://www.rdocumentation.org/packages/raster/topics/raster">raster::raster()</a></code> respectivement pour des cartes vectoriels ou raster.</p>
<p>On peut aussi bénéficier de la fonction <code><ahref="../reference/layer_type.html">layer_type()</a></code> pour charger tout type des cartes automatiquement, comme démontré dans le code commenté.</p>
<imgsrc="mapMCDA_files/figure-html/unite-epidemiologique-1.png"alt="Unités épidémiologiques d'exemple pour le Cameroun."width="700"><pclass="caption">
Unités épidémiologiques d’exemple pour le Cameroun.
...
...
@@ -142,17 +142,17 @@ Unités épidémiologiques d’exemple pour le Cameroun.
<p>Chaque facteur de risque varie dans une échelle qui lui est propre. La densité animale, par exemple, varie entre 0 et presque 5500 têtes par <spanclass="math inline">\(km^2\)</span>.</p>
<p>Pour les cartes <em>vectorielles</em>, qui representent la localisation des entités spatiales telles que lacs ou forêts, on considère la <strong>distance</strong> à dites entités.</p>
<p>Ce package utilise pour l’instant une fonction linèaire pour la mise en échelle. Cependant, la rélation peut être directe ou inverse.</p>
<aclass="sourceLine"id="cb14-9"data-line-number="9"><spanclass="dt">ylab =</span><spanclass="kw"><ahref="https://www.rdocumentation.org/packages/base/topics/c">c</a></span>(<spanclass="st">"Échelle de risque"</span>)</a>
<spanclass="dt">ylab =</span><spanclass="kw"><ahref="https://www.rdocumentation.org/packages/base/topics/c">c</a></span>(<spanclass="st">"Échelle de risque"</span>)
<imgsrc="mapMCDA_files/figure-html/mise-echelle-1.png"alt="Mis en échelle directe ou inverse."width="384"><pclass="caption">
Mis en échelle directe ou inverse.
...
...
@@ -160,27 +160,27 @@ Mis en échelle directe ou inverse.
</div>
<p>La fonction <code><ahref="../reference/risk_layer.html">risk_layer()</a></code> calcule la carte de risque associé à un facteur concret dans une échelle donnée. Par défault c’est entre 0 et 100. Pour inverser la rélation il suffit de passer les limits dans l’ordre inverse, e.g. <code><ahref="https://www.rdocumentation.org/packages/base/topics/c">c(100, 0)</a></code>.</p>
<p>D’ailleurs, elle utilise la carte d’unités épidémiologiques pour établir les limits de calcul de risque.</p>
<p>On voudrait examiner les cartes de risque ainsi calculées. Mais pour cela, il faut qu’elles soient <em>alignées</em>. C’est à dire, qu’elles aient les mêmes <em>extents</em>, <em>résolutions</em> et <em>projections</em>.</p>
<p>On peut se servir de la fonction <code><ahref="../reference/align_layers.html">align_layers()</a></code> qui arrange tout ça pour nous. Noter que ce pas n’est pas nécessaire pour continuer, car cette fonction est automatiquement executée si besoin.</p>
<imgsrc="mapMCDA_files/figure-html/align-layers-1.png"alt="Niveaux de risque associé à chaque facteur."width="576"><pclass="caption">
Niveaux de risque associé à chaque facteur.
...
...
@@ -192,35 +192,35 @@ Niveaux de risque associé à chaque facteur.
<ahref="#ponderation-des-facteurs-de-risque"class="anchor"></a>3. Pondération des facteurs de risque</h1>
<p>Il y a 3 facteurs à considérer. Il faut comparer 2-à-2 leurs importances relatives en une échelle de 1 à 9 et representer ces rélations en une matrice qui doit avoir des 1 dans son diagonale.</p>
<p>Noter que les élements symmetriques doivent être réciproques.</p>
@@ -229,7 +229,7 @@ Niveaux de risque associé à chaque facteur.
</table>
<p>Dans cette exemple, on considère que la densité animale est 6 fois plus importante que la distance aux points d’eau, et 4 fois plus importante que la distance aux parcs. Au même temps, que la distance aux points d’eau est 3 fois plus importante que celle aux parcs.</p>
<p>Le système calcule les coefficients de pondération les plus consistents avec ces valorations par paires, avec la fonction <code><ahref="../reference/compute_weights.html">compute_weights()</a></code>.</p>
<imgsrc="mapMCDA_files/figure-html/compute-weights-1.png"alt="Pondération des facteurs de risque."width="384"><pclass="caption">
Pondération des facteurs de risque.
...
...
@@ -240,15 +240,15 @@ Pondération des facteurs de risque.
<h1class="hasAnchor">
<ahref="#calcul-de-la-carte-de-risque"class="anchor"></a>4. Calcul de la carte de risque</h1>
<p>La fonction <code><ahref="../reference/wlc.html">wlc()</a></code> (pour weighted linear combination) combine tous les facteurs de risque en utilisant les poids calculés précedament, et produise une carte de risque qui couvre toute la région.</p>