Update documentation

Author: spatialthoughts

_images/1e11784ad3b12a6d90fa016c327bfa1dd66386e2615f36cb2b47f5b09d0b29e9.png

@@ -501,6 +501,35 @@ <h2>Setup and Data Download<a class="headerlink" href="#setup-and-data-download"
 </pre></div>
 </div>
 </div>
+<div class="cell_output docutils container">
+<div class="output stderr highlight-myst-ansi notranslate"><div class="highlight"><pre><span></span>/usr/local/lib/python3.11/dist-packages/cudf/utils/_ptxcompiler.py:64: UserWarning: Error getting driver and runtime versions:
+
+stdout:
+
+
+
+stderr:
+
+Traceback (most recent call last):
+  File &quot;&lt;string&gt;&quot;, line 4, in &lt;module&gt;
+  File &quot;/usr/local/lib/python3.11/dist-packages/numba_cuda/numba/cuda/cudadrv/driver.py&quot;, line 314, in __getattr__
+    raise CudaSupportError(&quot;Error at driver init: \n%s:&quot; %
+numba.cuda.cudadrv.error.CudaSupportError: Error at driver init: 
+
+CUDA driver library cannot be found.
+If you are sure that a CUDA driver is installed,
+try setting environment variable NUMBA_CUDA_DRIVER
+with the file path of the CUDA driver shared library.
+:
+
+
+Not patching Numba
+  warnings.warn(msg, UserWarning)
+/usr/local/lib/python3.11/dist-packages/cudf/utils/gpu_utils.py:62: UserWarning: Failed to dlopen libcuda.so.1
+  warnings.warn(str(e))
+</pre></div>
+</div>
+</div>
 </div>
 <div class="cell docutils container">
 <div class="cell_input docutils container">
@@ -613,7 +642,7 @@ <h2>Procedure<a class="headerlink" href="#procedure" title="Link to this heading
 <li><p><strong>Growth</strong>: Small Positive Change between 100 to 1000.</p></li>
 <li><p><strong>High Growth</strong>: Large Positive Change greater than 1000.</p></li>
 </ul>
-<p>We use XArray Spatial’s <a class="reference external" href="https://xarray-spatial.org/reference/_autosummary/xrspatial.classify.reclassify.html"><code class="docutils literal notranslate"><span class="pre">reclassify()</span></code></a> method to perform the reclassification from continuous values to 4 discrete classes.</p>
+<p>We use XArray Spatial’s <a class="reference external" href="https://xarray-spatial.readthedocs.io/en/stable/reference/_autosummary/xrspatial.classify.reclassify.html"><code class="docutils literal notranslate"><span class="pre">reclassify()</span></code></a> method to perform the reclassification from continuous values to 4 discrete classes.</p>
 <div class="cell docutils container">
 <div class="cell_input docutils container">
 <div class="highlight-ipython3 notranslate"><div class="highlight"><pre><span></span><span class="n">class_bins</span> <span class="o">=</span> <span class="p">[</span><span class="o">-</span><span class="mi">100</span><span class="p">,</span> <span class="mi">100</span><span class="p">,</span> <span class="mi">1000</span><span class="p">,</span> <span class="n">np</span><span class="o">.</span><span class="n">inf</span><span class="p">]</span>
@@ -624,20 +653,31 @@ <h2>Procedure<a class="headerlink" href="#procedure" title="Link to this heading
 </div>
 </div>
 </div>
-<p>The result is a 2D array of reclassified values. We can visualize it the same way as before. The <code class="docutils literal notranslate"><span class="pre">imshow</span></code> method supports only a colorbar legend which is not appropriate for a discrete classified raster such as ours. We use Matplotlib’s <a class="reference external" href="https://matplotlib.org/stable/api/_as_gen/matplotlib.patches.Patch.html"><code class="docutils literal notranslate"><span class="pre">Patch()</span></code></a> method to create a patch with appropriate labels and colors as described in Matplotlib’s <a class="reference external" href="https://matplotlib.org/stable/tutorials/intermediate/legend_guide.html">Legend guide</a>.</p>
+<p>The result is a 2D array of reclassified values. We can visualize it the same way as before.</p>
+<p>Since we have discrete pixel values, we can assign a specific color to each class using the <code class="docutils literal notranslate"><span class="pre">levels</span></code> parameter. The <code class="docutils literal notranslate"><span class="pre">levels</span></code> list defines the boundary of each interval and the <code class="docutils literal notranslate"><span class="pre">colors</span></code> list defines the colors assigned to each interval.</p>
+<ul class="simple">
+<li><p>Interval 1: Values between 1 and 2 (Decline)</p></li>
+<li><p>Interval 2: Values between 2 and 3 (Neutral)</p></li>
+<li><p>Interval 3: Values between 3 and 4 (Growth)</p></li>
+<li><p>Interval 4: Values between 4 and 5 (High Growth)</p></li>
+</ul>
+<p>The <code class="docutils literal notranslate"><span class="pre">imshow</span></code> method supports only a colorbar legend which is not appropriate for a discrete classified raster such as ours. We use Matplotlib’s <a class="reference external" href="https://matplotlib.org/stable/api/_as_gen/matplotlib.patches.Patch.html"><code class="docutils literal notranslate"><span class="pre">Patch()</span></code></a> method to create a patch with appropriate labels and colors as described in Matplotlib’s <a class="reference external" href="https://matplotlib.org/stable/tutorials/intermediate/legend_guide.html">Legend guide</a>.</p>
 <div class="cell docutils container">
 <div class="cell_input docutils container">
 <div class="highlight-ipython3 notranslate"><div class="highlight"><pre><span></span><span class="n">fig</span><span class="p">,</span> <span class="n">ax</span> <span class="o">=</span> <span class="n">plt</span><span class="o">.</span><span class="n">subplots</span><span class="p">(</span><span class="mi">1</span><span class="p">,</span><span class="mi">1</span><span class="p">)</span>
 <span class="n">fig</span><span class="o">.</span><span class="n">set_size_inches</span><span class="p">(</span><span class="mi">15</span><span class="p">,</span><span class="mi">7</span><span class="p">)</span>
 
-<span class="n">levels</span> <span class="o">=</span> <span class="n">class_values</span>
+<span class="n">levels</span> <span class="o">=</span> <span class="n">np</span><span class="o">.</span><span class="n">linspace</span><span class="p">(</span>
+    <span class="nb">min</span><span class="p">(</span><span class="n">class_values</span><span class="p">),</span> <span class="nb">max</span><span class="p">(</span><span class="n">class_values</span><span class="p">)</span> <span class="o">+</span> <span class="mi">1</span><span class="p">,</span> <span class="nb">len</span><span class="p">(</span><span class="n">class_values</span><span class="p">)</span> <span class="o">+</span> <span class="mi">1</span><span class="p">)</span>
+
 <span class="n">colors</span> <span class="o">=</span> <span class="p">[</span><span class="s1">&#39;#3288bd&#39;</span><span class="p">,</span> <span class="s1">&#39;#e0e0e0&#39;</span><span class="p">,</span> <span class="s1">&#39;#fdae61&#39;</span><span class="p">,</span> <span class="s1">&#39;#d7191c&#39;</span><span class="p">]</span>
 <span class="n">labels</span> <span class="o">=</span> <span class="p">[</span><span class="s1">&#39;Decline&#39;</span><span class="p">,</span> <span class="s1">&#39;Neutral&#39;</span><span class="p">,</span> <span class="s1">&#39;Growth&#39;</span><span class="p">,</span> <span class="s1">&#39;High Growth&#39;</span><span class="p">]</span>
 
 <span class="n">change_class</span><span class="o">.</span><span class="n">plot</span><span class="o">.</span><span class="n">imshow</span><span class="p">(</span>
     <span class="n">ax</span><span class="o">=</span><span class="n">ax</span><span class="p">,</span>
     <span class="n">add_colorbar</span><span class="o">=</span><span class="kc">False</span><span class="p">,</span>
-    <span class="n">levels</span><span class="o">=</span><span class="p">[</span><span class="mi">1</span><span class="p">,</span> <span class="mi">2</span><span class="p">,</span> <span class="mi">3</span><span class="p">,</span> <span class="mi">4</span><span class="p">,</span> <span class="mi">5</span><span class="p">],</span>
+    <span class="n">levels</span><span class="o">=</span><span class="n">levels</span><span class="p">,</span>
+    <span class="n">vmin</span><span class="o">=</span><span class="mi">1</span><span class="p">,</span>
     <span class="n">colors</span><span class="o">=</span><span class="n">colors</span><span class="p">)</span>
 
 
@@ -654,7 +694,7 @@ <h2>Procedure<a class="headerlink" href="#procedure" title="Link to this heading
 </div>
 </div>
 <div class="cell_output docutils container">
-<img alt="../_images/43dff7fffad3067d916eb2b8fbb072cf93a0e6d5a630d371b7ef84259669e674.png" src="../_images/43dff7fffad3067d916eb2b8fbb072cf93a0e6d5a630d371b7ef84259669e674.png" />
+<img alt="../_images/1e11784ad3b12a6d90fa016c327bfa1dd66386e2615f36cb2b47f5b09d0b29e9.png" src="../_images/1e11784ad3b12a6d90fa016c327bfa1dd66386e2615f36cb2b47f5b09d0b29e9.png" />
 </div>
 </div>
 <p>Last step is to save the results to disk as a GeoTiff file.</p>