From 46b16c917e404befe471998f56ba0fe6dea66001 Mon Sep 17 00:00:00 2001
From: Dane Springmeyer <dane@dbsgeo.com>
Date: Tue, 3 Apr 2012 16:06:36 -0700
Subject: [PATCH] support nodata for paletted images and allow user to set
 nodata on-the-fly - closes #1160 and #1161 - refs #688 and refs #730 and refs
 #50 and refs #1018

---
 CHANGELOG.md                                  |   4 +
 plugins/input/gdal/gdal_datasource.cpp        |   9 +-
 plugins/input/gdal/gdal_datasource.hpp        |   1 +
 plugins/input/gdal/gdal_featureset.cpp        |  79 +++++++++++++++---
 plugins/input/gdal/gdal_featureset.hpp        |   5 +-
 tests/data/good_maps/tiff_colortable.xml      |  16 ++++
 tests/data/good_maps/vrt_colortable.xml       |  15 ++++
 tests/data/raster/dataraster.vrt              |  45 ++++++++++
 .../images/support/tif_colortable.png         | Bin 0 -> 26726 bytes
 .../images/support/vrt_colortable.png         | Bin 0 -> 7308 bytes
 tests/python_tests/raster_colormapped_test.py |  38 +++++++++
 11 files changed, 196 insertions(+), 16 deletions(-)
 create mode 100644 tests/data/good_maps/tiff_colortable.xml
 create mode 100644 tests/data/good_maps/vrt_colortable.xml
 create mode 100644 tests/data/raster/dataraster.vrt
 create mode 100644 tests/python_tests/images/support/tif_colortable.png
 create mode 100644 tests/python_tests/images/support/vrt_colortable.png
 create mode 100644 tests/python_tests/raster_colormapped_test.py

diff --git a/CHANGELOG.md b/CHANGELOG.md
index 4e42b0dda..e84eb42f1 100644
--- a/CHANGELOG.md
+++ b/CHANGELOG.md
@@ -9,6 +9,10 @@ For a complete change history, see the SVN log.
 
 ## Mapnik 2.1.0
 
+- GDAL: allow setting nodata value on the fly (will override value if nodata is set in data) (#1161)
+ 
+- GDAL: respect nodata for paletted/colormapped images (#1160)
+
 - PostGIS: the primary key, for tables containing one, is now auto-detected allowing for globally unique feature id values (#804)
 
 - Fix Markers rendering so that ellipse height/width units are pixels (previously were unintentially radii)
diff --git a/plugins/input/gdal/gdal_datasource.cpp b/plugins/input/gdal/gdal_datasource.cpp
index 1bdf85ce1..29303bd5e 100644
--- a/plugins/input/gdal/gdal_datasource.cpp
+++ b/plugins/input/gdal/gdal_datasource.cpp
@@ -78,7 +78,8 @@ inline GDALDataset* gdal_datasource::open_dataset() const
 gdal_datasource::gdal_datasource(parameters const& params, bool bind)
     : datasource(params),
       desc_(*params.get<std::string>("type"), "utf-8"),
-      filter_factor_(*params_.get<double>("filter_factor", 0.0))
+      filter_factor_(*params_.get<double>("filter_factor", 0.0)),
+      nodata_value_(params_.get<double>("nodata"))
 {
 #ifdef MAPNIK_DEBUG
     std::clog << "GDAL Plugin: Initializing..." << std::endl;
@@ -242,7 +243,8 @@ featureset_ptr gdal_datasource::features(query const& q) const
                                               nbands_,
                                               dx_,
                                               dy_,
-                                              filter_factor_));
+                                              filter_factor_,
+                                              nodata_value_));
 }
 
 featureset_ptr gdal_datasource::features_at_point(coord2d const& pt) const
@@ -261,5 +263,6 @@ featureset_ptr gdal_datasource::features_at_point(coord2d const& pt) const
                                               nbands_,
                                               dx_,
                                               dy_,
-                                              filter_factor_));
+                                              filter_factor_,
+                                              nodata_value_));
 }
diff --git a/plugins/input/gdal/gdal_datasource.hpp b/plugins/input/gdal/gdal_datasource.hpp
index 4aa841371..4b2d2ea9d 100644
--- a/plugins/input/gdal/gdal_datasource.hpp
+++ b/plugins/input/gdal/gdal_datasource.hpp
@@ -57,6 +57,7 @@ private:
     mutable int nbands_;
     mutable bool shared_dataset_;
     double filter_factor_;
+    boost::optional<double> nodata_value_;
     inline GDALDataset* open_dataset() const;
 };
 
diff --git a/plugins/input/gdal/gdal_featureset.cpp b/plugins/input/gdal/gdal_featureset.cpp
index 7546adf3f..14bf11ee5 100644
--- a/plugins/input/gdal/gdal_featureset.cpp
+++ b/plugins/input/gdal/gdal_featureset.cpp
@@ -53,7 +53,8 @@ gdal_featureset::gdal_featureset(GDALDataset& dataset,
                                  int nbands,
                                  double dx,
                                  double dy,
-                                 double filter_factor)
+                                 double filter_factor,
+                                 boost::optional<double> const& nodata)
     : dataset_(dataset),
       ctx_(boost::make_shared<mapnik::context_type>()),
       band_(band),
@@ -65,6 +66,7 @@ gdal_featureset::gdal_featureset(GDALDataset& dataset,
       dy_(dy),
       nbands_(nbands),
       filter_factor_(filter_factor),
+      nodata_value_(nodata),
       first_(true)
 {
     ctx_->push("value");
@@ -241,8 +243,17 @@ feature_ptr gdal_featureset::get_feature(mapnik::query const& q)
 
                 float* imageData = (float*)image.getBytes();
                 GDALRasterBand * band = dataset_.GetRasterBand(band_);
-                int hasNoData;
-                double nodata = band->GetNoDataValue(&hasNoData);
+                int hasNoData(0);
+                double nodata(0);
+                if (nodata_value_)
+                {
+                    hasNoData = 1;
+                    nodata = *nodata_value_;
+                }
+                else
+                {
+                    nodata = band->GetNoDataValue(&hasNoData);
+                }
                 band->RasterIO(GF_Read, x_off, y_off, width, height,
                                imageData, image.width(), image.height(),
                                GDT_Float32, 0, 0);
@@ -340,8 +351,21 @@ feature_ptr gdal_featureset::get_feature(mapnik::query const& q)
 #ifdef MAPNIK_DEBUG
                     std::clog << "GDAL Plugin: processing rgb bands..." << std::endl;
 #endif
-                    int hasNoData;
-                    float nodata = red->GetNoDataValue(&hasNoData);
+                    int hasNoData(0);
+                    double nodata(0);
+                    if (nodata_value_)
+                    {
+                        hasNoData = 1;
+                        nodata = *nodata_value_;
+                    }
+                    else
+                    {
+                        nodata = red->GetNoDataValue(&hasNoData);
+                    }
+                    if (hasNoData)
+                    {
+                        feature->put("NODATA",nodata);
+                    }
                     GDALColorTable *color_table = red->GetColorTable();
 
                     if (! alpha && hasNoData && ! color_table)
@@ -380,12 +404,25 @@ feature_ptr gdal_featureset::get_feature(mapnik::query const& q)
 #ifdef MAPNIK_DEBUG
                     std::clog << "GDAL Plugin: processing gray band..." << std::endl;
 #endif
-                    int hasNoData;
-                    float nodata = grey->GetNoDataValue(&hasNoData);
+                    int hasNoData(0);
+                    double nodata(0);
+                    if (nodata_value_)
+                    {
+                        hasNoData = 1;
+                        nodata = *nodata_value_;
+                    }
+                    else
+                    {
+                        nodata = grey->GetNoDataValue(&hasNoData);
+                    }
                     GDALColorTable* color_table = grey->GetColorTable();
 
                     if (hasNoData && ! color_table)
                     {
+                        if (hasNoData)
+                        {
+                            feature->put("NODATA",nodata);
+                        }
 #ifdef MAPNIK_DEBUG
                         std::clog << "\tno data value for layer: " << nodata << std::endl;
 #endif
@@ -422,23 +459,41 @@ feature_ptr gdal_featureset::get_feature(mapnik::query const& q)
 #ifdef MAPNIK_DEBUG
                         std::clog << "GDAL Plugin: Loading colour table..." << std::endl;
 #endif
+                        unsigned nodata_value = static_cast<unsigned>(nodata);
+                        if (hasNoData)
+                        {
+                            feature->put("NODATA",static_cast<int>(nodata_value));
+                        }
                         for (unsigned y = 0; y < image.height(); ++y)
                         {
                             unsigned int* row = image.getRow(y);
                             for (unsigned x = 0; x < image.width(); ++x)
                             {
                                 unsigned value = row[x] & 0xff;
-                                const GDALColorEntry *ce = color_table->GetColorEntry(value);
-                                if (ce )
+                                if (hasNoData && (value == nodata_value))
                                 {
-                                    // TODO - big endian support
-                                    row[x] = (ce->c4 << 24)| (ce->c3 << 16) |  (ce->c2 << 8) | (ce->c1) ;
+                                    // make no data fully alpha
+                                    row[x] = 0;
+                                }
+                                else
+                                {
+                                    const GDALColorEntry *ce = color_table->GetColorEntry(value);
+                                    if (ce)
+                                    {
+                                        // TODO - big endian support
+                                        row[x] = (ce->c4 << 24)| (ce->c3 << 16) |  (ce->c2 << 8) | (ce->c1) ;
+                                    }
+                                    else
+                                    {
+                                        // make lacking color entry fully alpha
+                                        // note - gdal_translate makes black
+                                        row[x] = 0;
+                                    }
                                 }
                             }
                         }
                     }
                 }
-
                 if (alpha)
                 {
 #ifdef MAPNIK_DEBUG
diff --git a/plugins/input/gdal/gdal_featureset.hpp b/plugins/input/gdal/gdal_featureset.hpp
index d4a2e9db8..b812da815 100644
--- a/plugins/input/gdal/gdal_featureset.hpp
+++ b/plugins/input/gdal/gdal_featureset.hpp
@@ -29,6 +29,7 @@
 
 // boost
 #include <boost/variant.hpp>
+#include <boost/optional.hpp>
 
 class GDALDataset;
 class GDALRasterBand;
@@ -47,7 +48,8 @@ public:
                     int nbands,
                     double dx,
                     double dy,
-                    double filter_factor);
+                    double filter_factor,
+                    boost::optional<double> const& nodata);
     virtual ~gdal_featureset();
     mapnik::feature_ptr next();
 private:
@@ -67,6 +69,7 @@ private:
     double dy_;
     int nbands_;
     double filter_factor_;
+    boost::optional<double> nodata_value_;
     bool first_;
 };
 
diff --git a/tests/data/good_maps/tiff_colortable.xml b/tests/data/good_maps/tiff_colortable.xml
new file mode 100644
index 000000000..c102ff363
--- /dev/null
+++ b/tests/data/good_maps/tiff_colortable.xml
@@ -0,0 +1,16 @@
+<?xml version="1.0" encoding="utf-8"?>
+<Map srs="+init=epsg:32630" background-color="rgb(255,255,255)">
+    <Style name="2011_5km_vrt_nodata_style">
+        <Rule>
+            <RasterSymbolizer/>
+        </Rule>
+    </Style>
+    <Layer name="2011_5km_vrt_nodata" srs="+init=epsg:32630">
+        <StyleName>2011_5km_vrt_nodata_style</StyleName>
+        <Datasource>
+            <Parameter name="file">../raster/dataraster.tif</Parameter>
+            <Parameter name="nodata">20</Parameter>
+            <Parameter name="type">gdal</Parameter>
+        </Datasource>
+    </Layer>
+</Map>
\ No newline at end of file
diff --git a/tests/data/good_maps/vrt_colortable.xml b/tests/data/good_maps/vrt_colortable.xml
new file mode 100644
index 000000000..d20367c72
--- /dev/null
+++ b/tests/data/good_maps/vrt_colortable.xml
@@ -0,0 +1,15 @@
+<?xml version="1.0" encoding="utf-8"?>
+<Map srs="+init=epsg:32630" background-color="rgb(255,255,255)">
+    <Style name="2011_5km_vrt_nodata_style">
+        <Rule>
+            <RasterSymbolizer/>
+        </Rule>
+    </Style>
+    <Layer name="2011_5km_vrt_nodata" srs="+init=epsg:32630">
+        <StyleName>2011_5km_vrt_nodata_style</StyleName>
+        <Datasource>
+            <Parameter name="file">../raster/dataraster.vrt</Parameter>
+            <Parameter name="type">gdal</Parameter>
+        </Datasource>
+    </Layer>
+</Map>
\ No newline at end of file
diff --git a/tests/data/raster/dataraster.vrt b/tests/data/raster/dataraster.vrt
new file mode 100644
index 000000000..26f6b5cfe
--- /dev/null
+++ b/tests/data/raster/dataraster.vrt
@@ -0,0 +1,45 @@
+<VRTDataset rasterXSize="2283" rasterYSize="1913">
+  <SRS>PROJCS[&quot;WGS 84 / UTM zone 30N&quot;,GEOGCS[&quot;WGS 84&quot;,DATUM[&quot;WGS_1984&quot;,SPHEROID[&quot;WGS 84&quot;,6378137,298.257223563,AUTHORITY[&quot;EPSG&quot;,&quot;7030&quot;]],AUTHORITY[&quot;EPSG&quot;,&quot;6326&quot;]],PRIMEM[&quot;Greenwich&quot;,0],UNIT[&quot;degree&quot;,0.0174532925199433],AUTHORITY[&quot;EPSG&quot;,&quot;4326&quot;]],PROJECTION[&quot;Transverse_Mercator&quot;],PARAMETER[&quot;latitude_of_origin&quot;,0],PARAMETER[&quot;central_meridian&quot;,-3],PARAMETER[&quot;scale_factor&quot;,0.9996],PARAMETER[&quot;false_easting&quot;,500000],PARAMETER[&quot;false_northing&quot;,0],UNIT[&quot;metre&quot;,1,AUTHORITY[&quot;EPSG&quot;,&quot;9001&quot;]],AUTHORITY[&quot;EPSG&quot;,&quot;32630&quot;]]</SRS>
+  <GeoTransform> -1.4637000000000000e+04,  5.0000000000000000e+02,  0.0000000000000000e+00,  4.8596780000000000e+06,  0.0000000000000000e+00, -5.0000000000000000e+02</GeoTransform>
+  <VRTRasterBand dataType="Int16" band="1">
+    <NoDataValue>-9.99000000000000E+02</NoDataValue>
+    <HideNoDataValue>0</HideNoDataValue>
+    <ColorInterp>Palette</ColorInterp>  
+    <ColorTable>
+        <Entry c1="0" c2="0" c3="0" c4="0"/>
+        <Entry c1="0" c2="0" c3="255"/>
+        <Entry c1="0" c2="255" c3="0"/>
+        <Entry c1="255" c2="0" c3="0"/>
+        <Entry c1="128" c2="128" c3="128"/>
+        <Entry c1="128" c2="128" c3="128"/>
+        <Entry c1="128" c2="128" c3="128"/>
+        <Entry c1="128" c2="128" c3="128"/>
+        <Entry c1="128" c2="128" c3="128"/>
+        <Entry c1="128" c2="128" c3="128"/>
+        <Entry c1="128" c2="128" c3="128"/>
+        <Entry c1="128" c2="128" c3="128"/>
+        <Entry c1="128" c2="128" c3="128"/>
+        <Entry c1="128" c2="128" c3="128"/>
+        <Entry c1="128" c2="128" c3="128"/>
+        <Entry c1="128" c2="128" c3="128"/>
+        <Entry c1="128" c2="128" c3="128"/>
+        <Entry c1="128" c2="128" c3="128"/>
+        <Entry c1="128" c2="128" c3="128"/>
+        <Entry c1="128" c2="128" c3="128"/>
+        <Entry c1="128" c2="128" c3="128"/>
+        <Entry c1="128" c2="128" c3="128"/>
+        <Entry c1="128" c2="128" c3="128"/>
+        <Entry c1="0" c2="0" c3="255"/>
+        <Entry c1="0" c2="255" c3="0"/>
+        <Entry c1="255" c2="0" c3="0"/>
+    </ColorTable>
+    <ComplexSource>
+      <SourceFilename relativeToVRT="1">dataraster.tif</SourceFilename>
+      <SourceBand>1</SourceBand>
+      <SourceProperties RasterXSize="2283" RasterYSize="1913" DataType="Int16" BlockXSize="256" BlockYSize="256" />
+      <SrcRect xOff="0" yOff="0" xSize="2283" ySize="1913" />
+      <DstRect xOff="0" yOff="0" xSize="2283" ySize="1913" />
+      <NODATA>-999</NODATA>
+    </ComplexSource>
+  </VRTRasterBand>
+</VRTDataset>
diff --git a/tests/python_tests/images/support/tif_colortable.png b/tests/python_tests/images/support/tif_colortable.png
new file mode 100644
index 0000000000000000000000000000000000000000..aa0021b38d731fbd7d0b259d07341fff6893cd8b
GIT binary patch
literal 26726
zcmeEt`9GB1ANRRtVQd*?7e?A_g+WEkD4`N6yNYf*S-P!}=Ax9lTWC?1FhzEeC0m-Q
z?%PsiPxk3nmSoG$n0b!x^DjI<KEE)p7uT8Ve9q^5&inIzzt3Nm<|ab?()<7bp_8UY
z=K#RKj~GDUfxi~Ry4L~lR6A*8XchK!I=!##*Wo|i0$+A^KPXV$7g?CGn9X?UH>~Wg
zs2_P7yN#zR&9`nmV+Xd$+4<n%oh)10aT1Y$!_3#N7Sti|Lll*?%nJbgMU0hzU-17w
z{67yugDrP+gu7W)HA&;H&;|y8ArWVozW$T@?})(}3YWbgah105<B8MTyt;evFU0Pm
zWu8Yrb2D$OKNkTnGyO*ytxgymA(})8|Epg{{e?V@J^#1gsC;l5J6)u@h095C=bK@K
z$_vA1alYpwu${H7v~FKc?}(D1Wk=&Z#!V6-nl5Y`fjuXE6N@ORP=o8gEFFwwFL7sF
zZf*T7#@nBIN~;V6q;0M~of;WYOFGO7Q-I%?XE!ZGXy#XvQBTzqi{yyWHm96u=T|h2
z!}n~*R32i!=t#`hq%!|qb)c<%a^fs!a@R5$n*)*Tt%+E{i6nQnQ224(b2a0axR-Kh
zQxwi9<!#=^Q7Zfg%V+&EOl9wNDMDgK#J%-{9~iA%7^o3O-o{g_DJRYkr~dy1t1`bH
zd9HOc=QoHNubjx^YI3*>G)|j8agL8P*T!YH(JVvLw&pSI?a!3YmWnhI`S2vp>0$6r
z?y^63`9xt`FJ-ff5+kRQzY%x`{^RR{dW2^A&v%oB&j$Wl9LFc`oW>&4HW{2L2FHD3
zw5@;E?_z$4`SAT}*Uxg;UO{;P*KE<d@(@e-KS9U!NFJcNz%b@^4<7vRmX@PVrI~QW
z5!6v+RBmaU!d(mf_||+yduy51V9@PO2p+{CkLO47Y2KVF;}0vRFz7|1r4ovPF>K7M
zM~aTRXGWpQAePDG2kmC-7PR#)iDDeWCs<n|0~evGc$uf#3XjMENd;8;+9LO{4Tf1%
zrrsV|ycNX*@U{MKxOSs<;gXxUvpDb+r`|-R7AOcPBeqAXZ6(`*B^^|J{bH?9Du;o=
zlmX|IYeC-`+d+W>35)L?*D)6ka`I3yD#@H#h9lyw7mnE)WgZtA*r-zl0ST`FT8J$r
z(FfRBKz-r5|0{ZSmyW{o9CM*d8+ud_O~0eWbv^K=ZAo(VYrYK?<pnbrxtbv0IJ0QY
z7ND;G75)v#x?U#7TDwiCZ}#YOc7(*;B2A2n(H?}tgnHHCcKuwqfdC=}s(<vgjO#2G
zRM7SSOM;ynr-vAj+14Y$@X#trVc?LSF`CEYq7gVHJbV}r`n1`@+m2vx)_`UISweeF
z7PLR}Cg$BnT3ZT5`!pRMhUQNqvV1thw(~Q84d&7AOa>C^^}VXV=*xnK@w&*utj_X9
zA`^+S2~wK=$8+e@<rmq$&rIw`UbW<$zbM<RIzc6N@8aGp)|)xe!g#-08BOtbtX}jL
z2L)vy6M!CdQWW6D$t_Gja{ZE$w0rPcY;;3hfu<}GS2>Zuc=P1@y}&i8g_YD9FTjdH
z2NJ6*mSoOz^++TaG`C5*Cc!chkDgAsa@@G1&gB>}=qYFQbTd_@g~1Q@<E*ftYCDPr
zC#49bf~B9HT0|rH#w1PccHiZU>E|BhNMR0d_9&h}ZC6~{aI(y*AdzaFXi0152@?RQ
z2?WTj=Vx&Xj;t5Fc1O^gQ7waOMVBK4!y4Wcl1Tg*sl?ih(BtoZDiH|LIg1-k7jnY5
zQyRSoMib0=!v^3JU!+P*Hym%L+7DNCvTfh}&C9RpZs)aL_~XduJ{QK@jp0UPl+_uZ
z3JcS*buPE2mCL!-U{%*L;<$OR^1kE8e_sR@lAOP=v!>XgY21@UbpdSQz3?TUlBQd(
zz}e&eE}UTjhjX)WcNQq^TvY?Hr?BXLT-arj51kyb7>ZEE2%cjlX9TR0T|^gm^&aQJ
zJa^#3of9(BIt>ilZ8kS~y>AaSyzOG|0d#x#WL)}8ff-d4SG1D<+1ctQwW406%Te2*
zM13>8z>FL5jD^M%qlfoLobK+5W9$I#WMpxL?0!;cBpOxjD#6@L8gsh*XG#SvX|e=7
zu-|W2<Ne?Np2#b?;c@*zhmCYquwv2<;LNt~AMzE)>{9fcWixOfDT{M}F!#Fd&+LKA
zoP8Ly@gK!H5(bba2=EQtWIq?gjYODKUK}fN1F8f`LsI`W6T$OjWA9ZF0P|a{Q4Eh3
zHG6#5;N*QX8KEu3j-+uihosj~bNI3a3@)*zby7JP16u!VS4Bgz=ps1H^y0>&Y71?j
zH*KYogO7?4>6MX5YoRH>nm+<G2disH2rv9zmo=@$C4eGz8XwqADCJKt+pFo4A&m=b
zWqcgm)p#F4W7j<Rdu$918=Mt=2YflA;Ap7Ht+nIr3cpi-doQdUV=6E70%r-FRdz=`
z9#v_j@FQnIX@9ylH@103hr0tXQ#abE%&$@iKgqMa;o#8Ho_g+drxI{JLo*sXT<*SI
zDX+_ntVf&CDB&k0Bj~AxOzoi*QQ8SdTQqv<Qc^yx8%J#q*_^*L{Ioez7q|*YWv!rX
z?gG1vPth_oxN%-Z(Lewk9{niUr8HK1s?f}7<5C0wQh?BZfFs;-aP=YL9J1$wD$W&U
z>DYWWt;wwWd)!@t>b>FVAT}l+6(vXvdl~gp$EKjbf<>k3zQpRDB&=I)&OdSRSnvWw
z4Vo}$32G8K{)5CMTbBxeqeNm#Lpz{<G%X+N_u}X&&}$FI?(FS3i13i|F0kd@bz5@*
zwVQ1VtM99Gfluvc#x+XBh(tVi={p!KspeYTL}*E!IuO?UTJ3jzF&gzspEoF#xTy4i
zYuBe?4XC;Ej_G|r&!)JQcDLu(*od-EnAAuM+KpL+YQ<sD5Pj===WHc?7B|W;g3m{O
ze4M=>e4(P4H;7JLbZw>;$)!96fIKFTM~8?3GAc3A7XeF~R(#{K7OghBLjYt<<zp~5
z2vy6+h7ahq(0DM}v;Xus4RMESOCDl%?FqWpZ24X3$Mz74aH#!48=z|AOl4WVgQA$D
z+{~2nW#(NSAz~Cr#lzKV89%xqWea{66j5`0GiJwC2}UIJ=C=)LJ$hCQeaFAtsGPo+
z2ZuMS_F<?oYaR&wkMCm5f_946z1N4Yanr}RAA~8_??1+~?t;dtPqltFDD_9^H=JaX
zH233z5J2gg&9pbInH)bIl~Nf{&_c`A2ouwA62jPc?5|$7Q>3W~un;g^y9z<Y?-Rp5
zF?BY!$syWLjZ3;r2eJcCwfTInJAy_00Xoaza{hk1F;lnvnws#0`>x}xPtL=99#m|M
zlE(K{2pQCRw~bKhEmnPu$_M-ilJXQR$U1vXfm&<h5%=~ZThX`gGUpN@?=#M>_bVbh
zQ(+7wdxj27-2biLJ}_yUE8zj0lhzz%ePi-X;%-d|Hkx16a`A2%<SITG9w4i|JA<|C
zV&o$GL6|0x8||%lf4y>iy4K*uRAoQAUY&Skn-JzG{hG&vUFqN5Uh@WZGkk=#X-x3w
z(I&lRqC^JRq|*L1^=X(ZUU}Y7ghY=Aj^2=u!uX4$g$Ny}e?bC(`laJ?F6&cH+;zPZ
zf>(YlS>x-<dy{f~VL`tYB?eeg=v055$-urBgK7)BJE=&YaOxTSY&6=FLm7DGE=1{k
zJinya#nY!6ztlgh3aWUSo*XR!n6}=Fa@Xz^ep>V}=}MVub>R#LPU>e4tD-!hGdN_)
zU-8N%Zmv?uI(NoALL8VM=9!S}6D~-Wza)g@rcpfF#l{W>CuXSPHm~R+=+skxz{){(
z0j5-^m=wsn1=ip2gX9>v4)XYd{kRqYT6H=^xI!Q|*)wJhi=i8?V~bu~?saLPg&Kn1
zs0=?&QVu1g>B7bYcd};ebNy1CAPCDrvSpE)d{8SI@MiC9NW%E-)EyEO&g!Uym!+lt
z9%^P|3cbqaGLPxe-t?Z=cU$g`Uq3~pmXv#g!>vXt2&=;7=IhfmP-Vd8Y_3+{y}~R~
z$k+ykmC)F0rm`=3uD*<x@VOgx6XX*EekD(40BI1n_r<@y><5QQu{Jf1|6SU5F0VQ~
z{Jjsy!N{3}K{9kv0*D)%#iOxjZKQqN_5$LZcQfr#i&OY7Rh(7aT?-ug)1^}PWZ`Hw
zqg*amHdIGIJ)|?E6|L|BN?3IPa0q8fE`-t1wO?4NaD$lz=!Voope2r@Z@e~LHNr6;
z46Nv-D+aIc1VBrPBRpz4;=)OBbG7({vr1uj;4<*QV*Dg+JPJRv>d`{`l~=3_{7x_l
zaO60^Gz8j<g23t+h!nC|VhiGmG+~)DB|2Bav>&zx+ktD`Z_SlDN_NgcL+0k{MfYt4
zRB}={HkszG()W;zh~>|D4jtzKGx2|(;J5oIZ4~dAwr%bJ?;BnBT4nSm*aj*@e0ad_
z<9|m3ULZy+{bTW9WseU>gpkJr{O>iZf?q@clpk}@jlXwqUc`cf1S>a;vh(H2Nn2Os
zx#F408tZ^ZO8nme9IW5a#~$5qs2i!mV3b?U$ld(lYtnizVCwqYgZ8)nUzE`ufR*O|
zRwG!Gu+!P9UAna%!poziL=Dpf10kT9@2+w>QRA<cBc)^BHEGzZCmasp0tPS%M3St3
z$l_1!3ez7OI}qw=ZYB<NhID=5NBO~=Bu9_{0DV_$JrxHuAuz$zwZQ0fcma!3*Tjh-
z$&@EAdvv`I2+$4D3OHx<Zk`>+<lfRGAV#5rAfWfL1embtD#xRuz`ccQML5jNqif%#
zcs&)5#V}{gn9XDUvYEZriteQt9`_t@;NE!<Bu&MFZJ5I#@$-1gg=ImUDGyB<cusc%
zs><`H2x@B6j0Jr^Xj${N$u~rC09Z!-VL{hM`qkXaNtP*mcE0TLu{PN!+W(bk-BBgB
zJWZP~9}D5xspspJ-7qZ_(2Hx&2a|8Y1@B1@oFX)t$#`b^cdwEGDLwZ9m>c#vhC%bz
zOj_7eB@E6@cSAr-ln8C?m<nJ0$3+bbnperB(r{hp7n8Or6#cyxA1Qh!7x?fQ3!ZNH
z>gQ_xR|Ee=dnF#-MCZ+;F6>41sxA@c0F?)qy+YdsjPC9aMNkh;CI;mt995)!p4+3Z
zu!A)pGoQhqT02|Xq^9w<BxcI4B=Cm*%=*fQDEb$6BO}xulko><RKsn|MUUkCT8<??
zBuuS6N$kCRnn)!hqrus3n8Vn|Eyde0S|FCnJOCqi*8O124unaTWu?p}B-3`wwBr$x
zh4G;6z;lc?shc_KfdK)>P&`48jHn1Q`WH=F<C+S@Qut$!J${5wd>qj6NbQir!x=Gk
zKjJ@zAk@4(aNzO08^XjQ<jkG{ggPLkjs;Sq2~6sbbMKlKcFRk)CGmjKwR->_at_na
zGFMiQ43Ce3ECT`HmnF~bX&HGa7JH2OOY|LIgTEUHaYYJH(w_%v{!{bvqF3y+EZ~*&
z7dKkrZ7&1e-&JoiT3gGr=pp3_77fn56HXY#>A@f@%9GeB2grx84*=6=#znsNPj7m7
zxk4*jk>o7(Hs~|i*lgf?^Y6h)t_*nDI9L_Wq4vB<bYqv3N(b1sG8|R#6oZHWdZTJF
zf3@E4SK`c5eWzCVB^9WjFm}e{+G{D(F^Yhgl4GyIl}Es22EPjIRT1@wLsZ{~@+IFK
z3sRMn5uBP2{mc_1O+*~DzYgFg0P5~?@@MXOQ2jC=YqqflN--E%IW2Ah|8*x_KVe*e
zJ$QNiHx<~5PBcS)Dgz4!A9M-=%d94dewn<Wb1j|<TDJ#GSC;7pakxC_#)0fCimWp$
z5-`R5uKtT!6HEKq9%qh{-QHDFFrb{6sF9qPZ$!rIguNe2`E@;!bcQc%>}|b>8m7h7
zNx6$7MTi9tf^USFAz=tp&))z-<ruWb60VsXi4GC=MO8wb)4p?}Uhf7sbq_%{;~xCF
zao7wIc1Dq~NL`3P+?yPo@Ji8Znfr0BhR(lk)f?~7$oN%(N?<1Pe2aJ!NmNW!EZ?bd
zb-n0J4NhMhS~qr7-O3sl#t$G&>-xtUk0uc<?T?pA15JC0{;|F*cb1O~PjU<Ig;I`@
z|KgPl69#rU1_rl7-oh;+?G-cG=F01?z!U&IM7sWsrMxaIMj988PS})<0eraXwteBl
zh1z;|+ph21m4bRc%f=(JfT+Z-Yi-i(=I^?cPY`^{55!WoO$EwyacyO4rW=ISXXQYm
zpwq^NGU4lGjyPzcqE+OVw?UhWtxS>xrpOveJxoZP*=0&m^s!<?pc@5d|2*k4a_azS
z8suf|0MwhnU)*^gup=fm{~2aL<_AwohCsk3TOo_P6Bse6g*)b=QTLpSSsc?dC=TiP
z9L9rD*5s*jZ*A6876t)I7q;;?2t`-7Mpv(+^Ql?Ar+>V%Rq<IDA%y7y&w<ksEsfEr
zUtlJJ%vw>kfsq2!WtdD4$^83O&IH^$65wb<n;UkwWU><fqV(FS(G@1s?9Q(kq;0=}
zu3M>hP+APAJ)8JQN|KBy;sUxb+Xby$Gk0Yct{kx@)nPD2hh>4L*F=gIdq(j1UV%ge
z{iw=}57fjkk-d>RpuMhUIvOpKn$_W&=_vzd*2Nt8lYtQm%>Hlp-UD7_up8yUCP^vH
zq!vp5YkMI<9N1NYj647w4r*e7m_6!wxAvHd0-Yad5m<6EDLi?%6##ntX-dg$)E*rE
zmu(S+9$0V^LN3^#ocL6l7|zIGJjLgae=AilhMom;#$up-%0(Nv@^@xnkX(TBAU%MZ
zzzt_eWOIE0RWn#v)JVaZ2bhz&J4nv3DJaQ|0={*&d_1s=O#y9-)Cy<qI^Esj=|#12
zH9m32s5%*fKpJNhrBHjlHT;avBT5>6;w8QmpyAft?2;8ZUOryQC0OkF5jj(qZol`#
zA-Cm}v;V0AyI@SJhR0o;%1{1=_{<E-%dGJbT2}rI2SLv>d;ui7TUEuG;^(Ernw<LH
z3n9V$q3QPas1#t4;OzHHe#S68jW2!poEj<0V`Gl>EGWV^xg$pbQ+CR34WbGRoLWX>
zxbByIX^b!LY_aHamES^Dnb0|Lbn4e+@a{9h?Y{sXldppwLsxV!@_`PY1M3+OV)fp!
z5?mDpO`;fg>otcY;+zqXOgO|324BTb=YjwlLs<RSB{0`AZBmYm&StMPyCaa&w62fl
z-By;rD%&uN$EfU}TcbW*SZP8WB~a<RsSOFkn?$wL7JGIEJJ@>)4~%jU8$cDu9Hn(F
zi^Zn;xH~+5zKt~SnIGi6#J9&~lDSt?=7L*go%`DmAdItIV7T5ZbvNnS*Rt6r`8G+U
z$sMD5ZKXPMNMPhA+iR(c2kGnvwq)iRUQ;PViJYi(H$5skrT=D&6$qWS0?aUd(@KT2
zfMpOW3_Q;x4-`BE%T9?AEmQFAB$5L!3lF9bzZiW!on08yFna)y_`r7DpmAG$Libg7
z<AZgrl;8bHIhpywUoWrP<oZvi_#1+lsHSez`XKWi`dF-qdQ2yf6e;YUt57l3mCDm(
zBnmrwdL~0#)>RC9Dg_n$E#Bd#z}gLx(|KZDkqVdpj5+-iT#nyPI719-_Q4&Y`z+F)
zuX&hoth27SHMOM*!~<CQs-II}@5=*wFOeMy-3q3}0C<Yplk&w$SkT#u>H?C($=Mjx
z8-(?1p^fau`xHDnA93Wy!ENf(27sm8%l}~|`p%)qYd@Q4PrkJ!(;C`6j@?B3&pL}^
zc3(Y@yChi(Gn2P3HCl?8__8l42DRbY&ky>wAb!m7Leq^i@zKAu(9mmt{QJtW3y3G$
zcY%)klqw40&^qH~9oDQ)0#ulX>;Xn&CqgA@-<PtoHkGM!sk9g0OlE~j-l0~)%onJ~
zPmC(D{Vrd*I;(k?s}H|sdw?>q>jL*7IF@GCiKmktr(elki))vbicp{yJV?M$c`!B@
zSPxQ;&}KeA+X*i8vXF)$EU@9ND|mK)(e2J!;Ox~0lMrHmnW*nHGfd-8Mxq1-UA5Bm
zla=<yy=R&!IEw+-9Hc^*==*cu1AIG&lIjrC;SohWb^r9IAMOfUDap!mLnhjWkI+7`
z5r72&=~D$f@$yEloIk79W&!Mu4kFA03soehBk`*c=omMVV+df(9r@JH^9IPnnX3}!
z;=hNc(Ei!^!VCs%*a6+YRq-W2b><7W62Onw5}+6``o#CeT0;6IQBZJ?(aH8Qd8~as
zVLot7UOVy=leRvYTUs$bVo7pCD$Gz<f{0{0o%Gx8Q@1T(wp1_$!C<d7jITjxKF{u2
zX=g1t?_-_ZnsK||SFC-9erW0*`wY<@X1C<rF;u(fzVFRmwl=5^La(j>AnP@)NqgDV
z`4iZ<0<H4pqw!Nro47@l))AjosBHu~NbkHLqQ<}Ajw=!=vw1bE?K)L+sX8{5vSK&2
zAsPjQz;hzp8*QSAJfK6BYSq7x1=(<<b-xGVya#7wl<@a^<|b`gzowgOx5^9D0ju+)
z@wGy_%43YN*6aHr$G$A=?z$VlT%O}}ne(`h`}@;DziUyj;3RE-or6Dm`iHW~SIB#j
zqB!~+!{R&b&O;h2mWo#Hd^E6wx>2Ox1#HBb0=Qv43x`9`%f3f}LOJ*_4whS2LplbN
zogWjomXJl`XLV8+IC5|=)q~R7vvQ!+1}Wk}WiS&bp6Of2Y?jC^nf$$3A;K+MGa)$_
z3C~?21}fuPs5h4nTR84vB_qIZ4+y>IrohAiUu{#^4rqQ=1)e3L+d&)=q6x@bIbK$O
z`+C{n>DuXh8lYzFl0v8Ew*#vWXzug{6Qy_WR#za`$iGk`xnvRrS3CxUXJ`@HotXl7
zl}AxGTYQQDy3zksZGSp~&If)0OaY|-5vpCg0!}nQ$_iG(U3ZUBgZ{c)!Ad;6+10_F
z&=?JUHX!4K08z{vKssY#6~+AN5}*3kebpW$j&0~6%rrb)abBK~=4kAs|Kq-^?*3;~
zJc)FKipQDuj{nkv{nOE$DeXo6uBuCc+Q8~P8u|ievwRoV;nmtthbnCMFscsP8^aCg
z_5{27MVF-c60=R)RLGku;_PnId+fpl^vQcKZPq^Mi^#A|;)^xMwj6``2*KqxOpya<
z4ru%ZHtUxM?LQ>ABV|B6n(AJ~d`KdBb1w97|8C(zwae}#R3%bHF)v5bojFg;wP?lt
zjDEHns6O-lq4Y2x(vh?q0>}5qv-d)r-|vko6G{QIZXX)7O^uLv7hCS`>n33Xm;|Ul
zgz4gpTAWU829DJD9lQY5tI`4;?&%7$9IZ|4zv-jH-lGX2@9!^d)*)yrju-62*_}R+
zHCBI)oT%{brp<m+GGe2@b>zfEJy<(UbXLQr$AH#<id&)<9qIt<Hl#DP2zfR+u>lS2
z%426**)`X=e)5C#Wu6M3-kppW)>)h0pt=GJ4s`6I`lMImFd3T?wF>CRhv~mj)EzyS
zbMF}X<)u9#7DL5h%IpCjVPdwvZAb`^?NM1GJ&MV{FiyHx7<KLb{D`+dSQNGvUn$n&
zik(L~&r!DlcM=nexwJo7Y5$VBzJ2F8)&suqL(;$95ltdw5~s_**x&+)%riF!iOy&E
zj6$v*l>)|uuY5p`kbGvtWa~I#4Djj{)+6nd?#rC*AgqN^Kx*&jb-o1>G=4?HUL7Xg
zlWYgH1E+ubnKBu<bj_Pm5f{sYh!agIvpjH3-2rPi;j1rvI4FX8vtI^mQi|2DC#=~b
z)Fm6{;3h_s=p04mA%rx0mshHORwib<OSEEKNzM>Yk(jB{&(?Q?yl7DDUrr^Fh17jG
zJGTaU2nULMfn_R>Qij<2TT%-_(8|#8d%dP~!Z^qN<Anes$qiJfqWgA~VK5bKb~m?}
zWy$vC({tbh)MmVnuwZ|=USWEc_3;Spz!I&nZiGUe7DC^yD}o}e$W~bI?0lCyN2?uq
z9C<EBG0ko}$K1Cy^<MD2olz6ub(7?KA{vi$_~<>D?#MY_6iwwf>c}M9NAJU;Jqgwj
z<~&|IsYnw7wHxf?-NRs_a1#*gD`&kIa%V5QA+;oy)q;G&@l<iBWDK-aU-R<{d}Mky
z8$$1?-c<flh}h8-wG@vWK3#+)fY`{}Y(?KHg#uW@Tg%*Y6d1YegUgW_Aq?k5NTfRB
ziA9?0rV4qF#pP*o3TPTnqHeRPIJJKts`r{7B<5h=5SeO8i%itx#J_jf(G+W7O;mTB
zJEZrsS7&uTnHEZ+(YJz{#EY64wbFRNo3Ytq>>QQapl*vYi@te`HC5JJRT3as*9F>x
z*)b!SflcF%cu-+XO+7F-9Pun%nMwY=;>V}DRV|9tueasNk7e3KVCl4`fweXs;OuKm
zN`1ZZx(5r+?S-wX>8F8;d*|Mq-errWz7}WpAG~}E<L8n+!Qm$Z`p~T=3<^n;yQwst
z`h!viW;CdT1EE@7tx5lMkQjqnAmAX*b$rGafv1Iv04YWEq6pD`TeZv<NCI@<WY;qy
zo=oCJgX;ING}z%0l`i(Bg<Q}+^n_vr7aH>uhZ1mBl8i>S2aGXRCC<Qc;hTo35ZF}-
zYgHV{SuN@)0`e+04h9$>lwxS@QHO%%7KzS-bFccnB>C2r{iTUfmvAE)z8Jn|Yt!Uu
za)&_ac8sR)Umz-&Ek~%%#U9nmv0ovyYby3Yp#&Q@odsA3;`Nc)?{(I}dlEs5)&hYg
zf><A#ZH`;K%s$;S4E8@0gM1@^)Oad)8?Zxg<-T9>;Jzwap=voG*dzdKZhsI!&RxX7
zf<!Q-JF6Y6C_~6l<7G$ycPst=1WAJy#sxUX2S0IuZ)=nt2)#ho8YD)cO=DDtd8ZI*
zX2jET;?zRE_8gw-k<$s4Ik@b6)O%qBfuy3<z4gnSVwh)-tPTMysk&aa9&qLnw5bn;
ztsSf=pns4%bYse{08mr6m;S^)|Lux73P@5cGa;rjp(&XiH+g<xc82!9IB)-G3>(kc
z=?*9qMZJsfjom_^4f?bCu$TcY$>$IZqN+GkK<^Q3AMa4Y7l$R@1pFW+4)Egpb=UnC
z+H%(O*Fsn7<#U#KsD`aZ62K}Ms!=4L<<v=!_<IqbJ>&bu-)+&{s+Or6z3@!j>4C5e
z%)L%g`^!-*F8?ug;EHK9@VQAZT5Id(bv0Bdw{MGEbNV?GL9v8UiVl~CTWf2`&fCAX
zzuUh>M(BL(2qDD%5`av@c3cZo<F2<WAzoO8ZQCXmieJR!La`EZvosvs-$TOIFvn4M
zmqVaPhx-tY694l>)cx0<L}%Y^M&$rD*o7Q5LTvX@(<RCK={sMv##Jcc^AAu5)cz+y
zu{79MF3THpO`zu?@x1jb0Z5F}VD;JSQh&qsz;K#u?Cai*(0m7urkb~njkt3Wwtb9S
zdB1)+#TcP}A{bqEzFzeBY+?=;G8|<9N<>K>vb2A`Ys@k(!pS{-p!;bD;{0R{iq~gd
z%36lVYNh<3dIRTa4-SGj-v@CS8y9R5pnqeA`*bo8e9Px!-9;J$l?>rgUyzFwHO{Eh
z;@D5%8IX8h<Y56CwWDQ3hD<E|#8z`@JD{}q`tOK?Ue>79tIx#*$z0@R@2fFxJz#0j
zwq1WS{^+Kd*TOE9hjd?^-D39PIXKNOj0k9AMlT*NS#xwGb&bx>gi^Lrb?4n=D^VV0
z#Q9F2?J;VI>!gN#iO97*^*Onbi~gP0bv&%P<2=P73er9+1ZLV=H&VjO<Q!hd=;d~_
zQPa#zXO0~<PY*+!du|(A-xXcH+L$F98p_y!juo)|k<>Gr<IK-ioGYF@_oaqrb!c94
zx>(+d^o$kha?mLXQMM>X6(p8HMXh`NO>rEG!*92{+DMAJDS@eQB$xndbkvSkCtbN#
z2<aVRIX68k2A%?>&}Et3lFV3SH%**VG5&;M4qckv@G3++?-V1<uQOtR*vxOljK@TU
z`01a=^#7)|wMLuspe;?A(Ysszw>*+=f_WHRhBxQ2k`<sX9wFyACGoX)4xU@@KTj}K
zK(T1ukO(M!X)!zQxp`*&dT|?+WrNNu@Pkf%WcJL>!-TJ1tM|_`Y=<IeKfxD^JsZ@$
zJqD+i;mHZL!Y~EZTO<vtzQqOqeUTqXtxB2{RcV}zfh=}I4>$|j4XXOutitcTQWIxg
zqtFD3cFEw~0JTu4T`;m?6`B(2N+N5vs&1+_hahM-zk5{Pg$I<CvZTbqs~-MX&P||F
zjYUDWD#{D8>kx0CzR`Em-NGr2bD9!5^GZC<G_#&}D}d)o$uX#DhjYWCV#CN#+c5@A
z^>$m@A<o9g#jSc->#_&9wm4zIlzE^CAvp@;O98VvRF?<~Jeb)mS`H|=Y>4LWC1HRF
z;jjqYb9S+C&s{1J%NfBV`%e6$^(G-9XH!4zU(<o6a0T?V9xMO{a61)HaZJ2VP0b(7
zSM@>R?*m&~lPQrL*~Rtj-k+P?iQr#RK<UHtA#F<xSU+L~KnyI#r-|)DjAI4;GIP4q
z>TGYqR}q_tvbk^i4c-e3QH>ijkIS~>ti-qdS$gd69>1({-7q0If95ZbZwQzTcDZ;S
zc2TdryAT+a4ZIy7rf-V{hInw{=ZuR&6&zGKIE6b&WL91dRro_3Xq{jBCFb7@$>=3y
z^lCsC2o+~K@Hr#QWPI|bsyPOo*PjeFB!l5Ur~;Vh!1DecC|1`w#dJc~Q=w7ea|;Y=
zPH2i>F9f7CPE||r7q$|x3{1y^&!yRWe9pz_)*nK<KX2&2m{LnyuL9^Lla49x>qj?J
z-tG%EsPUpVX;FKALMkMT1#-A&;h;XdgWmv@!F%XnT|c|<Jf|m%D*`6C9=3#l#(&g!
zzx)01!40VvLxQZq6A`;>k&eA-jjNMODNvviXJXI^9zzkroI9a*<7}SQw9snUH9g2a
z(y<uj^P87|-R0H6r~d97+S46xTeZ%K9>#*uftnv^XC!i%AbF_fTMph`Q;bZcyQOCP
zoUJvlnL4SPnWB&66}<>9Cpx=r4eJR}z}!46*uQR}J!5SM^@O*RW~6VkQ}$|`o)BFO
z_Zgr%LQF%Z9r7}Y`NC5pEBH!=PbB+SjZY${#SNJm89G$|-%F+Gr{M^-U$Uw8vZdAj
z?>FD-(h>QRg8M>%bzxy#_nI)SC6#*$wvC_)R)1uwT9w*FMS~)JCTZ(>Q?&&+RE(&8
z!(5%y3+tl%!B4ILY(S&MqT<-WVk|`#;%@Zic1a)HuEvS<pRoDl!F0|RDS_}vF8~#c
zMRN2S!vV&}^V*p;tu5uWH91t0(tFYal2r1a<_ecVtRyW+7|pY{Co)q6ZTC0z90~g=
zTcjJ^5?pezE@SE|<ldK{lSmL}X*x9{^Gp%8T?su!N}*T?y*!O0TVD(><m=^PDy)j#
z;6@(FirvE7Jed6iPa?G(BZ_m^l<%*!VDPSR$-RH#0F)~EjcnAIW_w|J-HL}jmccP{
z;6^G!k;15+e~ZArxT4qYGDigD>2P~MVlytx(D}ORA#Gyu4j7n4efm{y{An)Sc7B1P
zplA~~b1cun2}%}2!oxZ78Ux~9kcdax0&C1u$@O{7K7l>fV-&}ohaeX8#quI9I!qUH
zm%Y%YR8FKv3P0~HgM5b<(8W{l!LpD{vC3ZNxgHk+E81>Xr;<+g9smjJm3bFWMQ+Q)
zu;!Hg!-V*C#j%o$!sX{NgZr)f^%K*!xBtXeFZ?(nu^9LyKGQf@oT^LvqDl2p+zTMG
zfm1vsMX<vIx&3wsA_hn);W{o<gepB9F7)QVSrNB;l)%lC34faj*oplboyh!ziyNb0
zk<prOfW@(|{XD7N7KOTVtO0tTQeg_`dN}<NegZ||w6Nd_JjhNsFr-0QZ^lwss)%I5
z=|b1s+YK+c#qXK?gv1+i2^Z_eOb9Gv?t?Wf+#x-1{m%1%kmF2M!LU)>g1nv=AJ{ck
zVqK7skfME~<eN!Jd5ediA^KRPbh)~>#i@WVR#mF|^Y~`p`Tgym#UvcQ7)fvXoItH1
zz?C@GxRsW<De-GKLhNfr45sfcH>C%2VPCOW%~ina%WlNnv>(ls`eA!1bY<<+F;Fa#
z$fBlhW98sfNm<hiD5Bryl0%hQazvXiTBzuNpEJDAOrUX5K0oF2nkXR;7ZrE4?buB`
zaCn#LzTSQ1Xxgf)z^41yvPCrTeWYUj<}``^j{3d)fq}dLp?WsXI_PBXX;0q|cUx(-
z8`KFASZW%*Z9T5*;#!hJWA+16FYeh(=jUB{-T-3$1SBZ#x20n91#tl>QTqukcU6YJ
z8|=ux1euhVm;QizuZ|nM64~yumGs>+NBVn+MBme2*tSfpY^67Fe!(uj@9}!#VA&)8
zq%{)*mu)$!L{<*|F!m7(-O#IYO~m4+j||Ha8r$yFJ=R)y0+^2Kd<%EP;h`k2l&X1Q
z$jdfl(R#0fRL!9RW3JrRUK9gH#_WA@>wCNuFA3nXtqyyP?Ie26kxqV8HkFcl=}skw
zn-9)+>DU5lD1%qu9HxV5Yf$MvVEE;5!E?>cHt4R?aBQIXaB^k1vScX3@Q0BPJgx~`
z!e0m<BLA!Yc6MlZo}AS<4yl~OrGvuWux~OrSFTu;dyIK<Y!(59{UZ-G<6AZ|{Wy!{
zK|FGaSgNVHoe&@e1TfY0gqgOW3&$i8(9iJ(K#*z(9iob2lK~cYqGO*G#S4PC=|=@+
zies4wDEU^R$!!X5h&iN3Z7b)+sY=FDqZez=lIVIO;g5I!Nn{mI(E;;alPVx7f<6C>
z@s)korHL&S=QSUCE|by_mf2tNoce$N*bV~5W-oqn@@0$lJ`R~FX7zthAy?6zpFUi1
z6)utnlj3#TMF%rC2qB$qUpmkEoqFEMp6%VrrO~Z?@Ba9?Mn%!mZNJn)2@fdLjT8q*
zr<yIaxnj`NfbJucqJZ`#i^tAS_m(~$xVoZNO^K*=P2A{nxLCFwRF6)~2tRZA-#GEx
zA|UogB}>*IQq$Zm0F+y3OV?>jHLzK8-E$|n-dz*>pe9^EFl#J=n~}ypCw{{iKzYzc
zJgbg!D;`1#gL8#vzG}ul98cChUn)1Naow9L2a>gEg<spsU$q?wbP@uiBlRq4dmE2;
zXp$JTy&OKjxsb_apMYY8+s)BM5PHVbx5;mPG`)6WM&Mc2)_e6F0Cq3pf4#peG8^_b
z@2{}CAG4z=36nkzBN}BT05x@O8Nt9F){_(<qpm$LFjqq3F6L4qJA^PRBQG0I7Ygqt
zGy$qrh3Bg}?oGL8umEO<xB~+zjB-mH`LU;4jT@)$EcdUm5`-DI1}$n=_B!K$UpFm{
zyPqHV@~7c+J#4H)xxD(uTj}U#oWBplR6gAF)^_arIhR*=KkeQh3p08lGbSyz&jv%)
z+SMB-4t5n1CXOC`g3$kG=wh3TVv2@vH9mJ8@4;56d+svF61GgMhu!p6auFu~HIvqk
zz#l_#4UXMO`EyQQ5LQx4zFChs@j*JZ-3pFV-OANrZ#lL?r4^i?-dTYs>c`z!V7hSp
z->o9>m&~{Np|9oX8?XhO^|5~Jn6=TgllxL+wdkAUd36;)|4;facdVy0NAt_oML4t2
zk;oO^lSJwh{$Ua`>JA?CyBj^7PnaDTS}{yG4;Rv$49Nk|MpzVWh3*7Q(aJ?%PN+<1
z?g$K!na&@FfG2eS0{C8MxC8!-&F51d7JI46ghc>_4u3$O?H!smOxQXc&Xvo~ojR|7
z#^K}8L!qXOthjmGghdsz5dq{eUmJ$H7zx!cUQD$$?3V|@T<f8c>mpFZG)g=MyNPom
z&PGF#GK*V~%EUWFbh&v<KILYc_<-Bzr17=XBQrYTB_a{!d+b4zyM#MJTDVZW4Nb#T
zn{l_?8qWnVMXHdU>;dRmg9=DWt(sb6(J%ACGPPR$W!Y-)3J);LKc_5z-EroUK*Hqd
zUVBkslgt=T&4{Xx;NIzoom7SmbZ=-Xq*$P{0(|ki^>|*fRakB7SU2y}DwO~Y=A{FR
zw#l@hzlWgLSZ+aUf#H(#<sSh^Xz4!qmI|rC;5YWD3^*{_TLzv?_(gTfDw>>Vg|2ut
zOE}0WXeNwz-G%p2iF?z(T3OlM16XR;M_UdCo_~JvTPtvdgi56taE{E#=#?!wQLCn&
z&I4;o4-%J#?GQuygSRc*+4jGI`r<Wo`9E3}g{HiP2ma==qCLvrm3zMids%`Ic7B!&
znZ2AVKWP#mYzn}Bs)87>b2AXGAR6TOM>gk~F|2c=T6z-F+Ee>dXs%O3rn1amrjP%Q
zeZ^huvS8#(OKZ7}4MWPs`Va_nI;e?PH^&Z+wVi<;JLbLr5YT>kD+#@F*PyZWz=svT
z%7_PRCJ5^`nm*Qcj}6<{)(1vWEF-{r838*kW%Q~%Dgqj&qx2)70fK&fYnMOFC1;#m
z^8lnjKJS}U8e@Ez(9h>8gXCuZr#SFnZ3NI?a2u1rpgZ-P2&vvaako#JQ!=H(aBEpl
z=)S!xQoHoVnx}YZQ`Kj6I6ppb-ZqV$yj*4FwW<O|G;0-v8Qmwf(A&v%5J&xZTZVSV
zV$j#+nI4%m0RHbN#d7nu+20g-<kQO&;^(ruB%T^xGg*;aNJ!o&(#(NqB;{-f4!D;I
zTLXGzrY4~J^J~2X?XanZw~K;eA)vmRDhp?To_;c`#TCRc**;rVtxsqcPsDcU8DOPw
zeh&*Jrc#E$egUMqOa6b)Xxt47>lw2^=pSI5i0NJOXJVz#{j-|CW-GJSM&Uly<lfDP
zT92~v*IWecHcwb{gqjMn(stonCJ$`40}7ocsBraMw6H+0YrU{#1dE7tt+Z%6j!}Fj
z$_w)NpjK=nGxZl>R`zJp_+dSXBhgzoJV?}a{TBItNYx*8u=R-8l$rjeImYORC>H`n
zZ;sT%c`NAjUlDVcM^HgQPmQ$T0#LRM%(-_?EVUxg2q}s)zG_^80AB2N=(nS{!VFMM
zy7<#Xuf~E5&C}J7Xbl$(pcjV-3m!S^Cu{)`&>7k&y_xU;o-R=9wd+HYDO!ar?d}m)
zZnJVlnQnF6!wsVa#%|d=u3Ft+6-Croe@m_&WU#k?5rqDB^z;T6Tv<rKpu1b!UrZtg
zLC}hC`$8}lT}b`uKiBn0JJN#pi3Zm@F;9oq>Ri?^tIbVv_fMfMbu%Sxc)S}dDKP1}
zDfL1DFV;-R7<JnIta*pGoHjM>`Wx_d=tc7;x`&#pvs7>)?h;yEz2`VtYqq{0&-<=3
zIO~OZs<<#J_0%R7TpvjjG57C!VgoIYnF!jqNDF1~q#yKack@;8z2LOTyQEdDLsQye
z>B~tXCTep}_vhV!2nw2%I4YVDBG2V00iW*{djB>}IAPJ04q7OUNB_9PSWj)qKauGg
zy^h={{WE8wViHsscSB6n&g+KncZ&CyRkekLoJD<iiryaLna$I{xp=dr+Z*JeZr+@V
z<JpSbSHEdGo5BlsY)rAxrvjaD+16v3QD}<Dd9eV9yOK5*LM1k>nW#0C17lMt(~^yk
z7>{5->zvc+bZQGe=P3j9LMVKHS8Ew=*Nc=}mDgd9P$%Ssjq+S#h<13P`w*)5gBSN1
zZk(WWhARikT{@eSGYq=NN#)Mhei}l+ZhtZZI#stA_Uq7yNbau?JM=pjkcwBFG}fQI
z=1k92(Z~OC%s*lK8UX1>lTErTXYo^SWUD+H4L-BUQ)9OA>y;R0iWt(d#NCH5Rbk(V
z+EnclV-Fc0O+8XIm8mo9^r}PHn46hDFa6W+@0FQ{+_3B1<_(;i8Dp1AY1zX}4Xm#(
z*RN|un%4cRc-OnT7$)?V_IZI~XKRS<H#Sd80soH7UiQ^2cts|1Bz7ibune%?8Z&<S
zP2}cvu3BjR&fE2GBpDx(e%r~OJstzSi`U<4p?L6kEkTiXV45;Klu*}U5wrGTm2z6q
zX!y`yKm8u=<pZ4g6%NmhXz2D58?zW+x6t9L0@v})P{1@NOcY_Ede%SoWO%x*wxg}(
z4W)%kHu5nvKe-Zt2X?B(AAW~C-xt!I`YvL&r3&uDvBgkPZmSs<dd`bZQt+^gp`4%K
z1*>K!3s0EKuLWX&>h|eUo8a2QGb1i@HkU^v?_b-P*txE3YIj6Gu@wWIUXB8mZeJkh
zd<*z~$P@wDN4On}N+2L~J*y$(3{2#Vefj3!--f*9ayAaG&(zJ{*r=*&<%88^zCdSk
zUkWV0Vb9khNTf780AE2s8~1!y4p9I25m+JAP<i~?3n{KACu@G>r(eXeX`!$;HbKdH
zaTK4`$CC*n_8*+QV2_W1*P1+}p^a;0Neg=?srO>TsuJ*w*bL0#4*Y7oKS01#tW&(%
z4n5q+MxBZ4R})Q7{0~w1`mfxA6*a8*r>X^c(|&XohYqEneUMuR$vrQ9QpuZUKPbGw
zh1VQ`Cw`vfW<u4Kn9oV>%T`(Qw1w6~PTfPC|IMUKw?-^|cH3y?Upo+pRqctZkSVC%
zMsAQlx#^Qd%O^7oonI3zQ&;Y+WgeWYO^nOLFngV7zQ>r=@3S|;ws;axd)XmAS$`hj
zK57NN;*m@OAoniEtwjK;!+K+D-zmTv!VzL<nOXJ=3<R3QYX4a*O(zY=%-MOR4}HAf
z%imQJe+uLci$Uo_gtTe}nl6plkFK<A`b=>q->7+Tw6MgJT$d__hmTfp^RFZK#%T4Q
zg%AnJxRMj!{(R%Ut$7<<l>426Ixe4k(r5fWBmn3xBg4%dZBJlPnHxvju4X{D8#FE4
znDnD%bWm)o|3yF*hKfa!x#cI$3`sAIC~`HHC_mSs=D9L;GT>T!E0n3`4%?5VYL9<3
zUvXVn@#EL>;nb0wZ$&kYW^J+BPjN{3!e&M{W$0UrCavEvkVqQ($E2=C7_g4i-Ww$M
zKRL(M+O%-8z@j%_-iHUtxxUY#wz)%Y@H!9k00Dp^ENBy=y~*NwFEmbY?Q6I0!i)wQ
zgX-2yna#l7FPD2yMquS~&3Kq(VcPqN5uahPem;;P9rf`F(;IvZR)e%hUk9b0djyYO
z{UDxC9@m}V2ZvoRN&FF%wzQuQ?B5+AeGjS!J{);47jR9_DPht!WYB)p$#=0c;ujC3
z<bL6zl0X1I&Zu|Q!-RWbedOW#_+!pVugHYjjlHoF+BcftD5@u>*yuvfp+9r<k;bwt
zF`!`W^z4h21^$HuDs!8poe&rz68rkvd^q<Pxk8|Gu*in}M1vN`X?vvgyUh%18gp}t
zeK0T+QA`NxhS~8Y;$*{{PPRof3h&kNAxRJ_tY%HKO-#62`YOTg>fE;N+8d`9#vea)
z=sI=B9p+f6RX3ry@x-Mz=%R&waClLL6`3grjQY>=Au3d72@n+h?w2UE1h)UYG9w3q
zWSmX{H1LGbu+MUTg-8_*7kXaFt8e7lRzM$r)MUn9HrD(4Yr#WrXg{!`qKBUEmDS#C
znleG?(ga;-?l)t$@B@|2+5Gk#sCRLT>-2{<cP?_1r2g&tPrsz7?&LWzVQ?r8zSspD
zqn1CRidixbZD&T^*(#dbd(LC5)!XsF+8y3K6d}$5V9VZTeZ5N+96BdczAb3!9RWP4
z-8eikJFq^?b-0`SM(ZO($f+)A@m)UGaazdkfJN9n&CWG<Cl@@hI!@Km1uQE<m^3JQ
zBW3BzrELS*KjDueA&0yg9})nEEmHfsxFc=Nl3dUA)HV4G1&B4otsA#&AeH<;nMB|9
zX6%9+bkf`5`5v!jZUYal+P=SKvpWxA9^Rhi!&L(;egJ(9!BzQn{x>=Mo0%qE<3>W4
zG}#ks&+6LlJYVLy^3ct=BM$>+f}!8wu2_zL{hh|e`jSd%j50r87NA0v*e>+<!<X-V
zO+l5%c5NB~-g6J;Kd>{Ob-YnbxF^X7@sFJuJp}vw;<7kJ>7KEY`muN~RL>@r5}9yU
z>5e%eFWI0}lcW`bpl(oncO6NCmrX(9viFlQ?&KiXb1|!7IcUWrv6vsMvg(e#KIB}d
z+7;Vy0K^6jt7z^hZ#<s$k)2f%x(oM|03a&GK>ZIYIau^6G?0UOmjNL#{2!&2ZL)uI
z<`XH?9XL0zg)O*FFjgVfsqx)PZU2~!-v})o0<6dpvCQM*7-i(YNq7^AD+ulMDF~=#
zL(nj@lC->kVy2NEej?5(Z3QZ$FU1<pi<LuJa>kp`R55wC$%E(>@lShs0VWEHbboLl
z<!1=ep~VX|UFY5w^7-qJP<fu7Pd@ANn(GW0ju5zsrLu}KsBG3|fzGwh5sP07%nW`z
zz`M(=XIP+G&NyR<L4Q?xxIyhm@(G9^qeU{vAfP@or1$-dfq+9U{{p<ZU{qKJ$g<xK
zfIM*Pqs|G$nHT4JNd<#AEIY%Vq;&*dL!kRQ;z7_R<>JfWdXD~kUhu}Ye2k0!$ACHK
z+4Ip0^1#+<jODl1DbZaB>&B4OXxp2Y5*T=b#IBX#(7MRQ`cOB}a0GS)yC<>694Iuj
z7w{|21$o;vR)1vIti($PjvXe?emgRv*FyPS(?&(DDWpmh@?spkjKY+HrnbSxF$8*>
z-hhC~^?2B`h#}z!IDxkBX=c|AB*9kAK}<V3@2tDY_o2JSF>AJNm#TX&;K7#-*C(Rh
z)AN$Lyg2JI3k7F9X!n9QMkWkEc_Hu=LPPGnpSdCMvzIrS`vRI?l>gwdzqy$F?SM;y
zm6wU)Mz~6jQ+Dq~1XiA6444hSG4)DnL0*B50okMdcrd68ra#{r+7u7EE$7F1vNf$(
z1*_Do<*K;+h%RQ|9>>0>EuWLS8STf0Z3<WHEXb2v^MUq+rRVQjl>xm96ORF}!EFOb
zT6o&>0Oc}}l;bZI@%<HaXP<Z>Vw=jnBEE|s$p|qjiL%auxUjTLMjQ;rp#3&X7g-Rh
z;(P$>YzX_Lb=a)XFWQ5NXw=L@9)DYwsP5BL&NHh${L$spFN3f~r7?<l6HJ-;)rut0
zIgLUx;Sjdd1__(NQ~e3uL;i-AQT;Ljx5_{}f~WxJL%5dsk1EEiDhSobC1~oK9;9Gg
zySX^Pn&FwK7c1oj2ZC4ZiB`}nwDXzSzNgv^7aI-+ua*<lWyHoVUgmr;eX~osKfLi*
zrKQHW#gnEJ2s0v13<QjbB-V)qd=Fth-MVrmd^Q3bpBTE^PUF<D7ssIy+UNTYH^VO_
zxm+Kfnz<1<6F&F;;um&}_i8R+`f^6iQANlg4>brf39!YT*lKYVA5GrM1mMdye}Fo`
zOBoh>fL*4)+fOG40caO^VoZJvSeN)L@gTK!L(Q|qz`0KzQfl#MdnTNwa28UVu(vjQ
z>befq&?6>fz`Q7&(WaBO+^h5cS>5-m5r|PjrkqOs?Z%L+F!^{u?asuPfc^>i4{1ZY
zl_|qSgH?wEs{z(=oFhA$Bx$xhFkP)AgDdX7{c|&NF+BN*L?N8bIvEi4zrf4u_UKs$
zA!za(Wmt{AjDop?+p{$XvFTshJCb3|5OlrA^#XLMHC0p3`%%m*0htlGYtV(dxSpc_
z`Bi$SC{tfZ-H+gE$PBZk>WVRkhI|KA(3j#9GbN<``EW7vuiR-%(FUy9R|23wlTWVP
zK?vy_zm&e*2AJ?uu-K6sY9BOSpT%x@1+}0w3QgrG`ckyIzjw&P&=La9!{^Ad2GdOi
z@EQu13`#$?6pT%UCNp?IUz@UY|M;6O_GEZi?eUC5JBT%pz3L{-QRyvXYml<@({1!m
zm=_ywyJWF-8-)*V)$K7jwKbm0zU!N`FUhvYPznAwKOyAY(&tqu78W8-u4fqTi>Vc<
zz~b@b+b>$B&wTG;(8qZC5*Q2ut^9FBR0$n9jcJ+*G2Bu;_qR&UHRsIwjzZ%J9_iy}
z%Ywv|&U{Ty>{?*jcTXa*nR>EO`*omdQ#rg4%s=Ofsh&g*2(%9!4F33WcW{;PIT&pd
zuN+0Ff*qYLhapm+zh^`7sexS%h$j~u9lwQ>5jjG8{NM2$CD4A1c~u+6sQL=imz~w1
z3e<4dU4&4)4~lN-4GCRwd{AKM2+<b-At|&U<7f4Mwk0z>6nA=^Lf@l!p5aM8KIngv
z+(NH>JbsC#emfT%b}^egc)?+)u!TLXMd<yWyJqWmfC7>;FR)6zUO6Ap#Tg!&%@(@4
zn^=mH&OIR_z(S3hV}N3qkynN`I}cp#Enb!isgl^5OQY8pV6MS|7hG|mDV4TZp0CCJ
z^}C(`FPpIJh;%VQ_kNw%5dK^cN68ykpHPc+P;+ct3>2s8@tEo<Nib&zG4R<3ae*)9
zZ$xh$@-|HNiXo~ktAH@WoPfUuEAkUd1~-w8ZJY#?U*4^VGMD;-S1j~441J{Ob@via
z!+FGMDPFwHR57EZMdjTc0+t%^&Ybg?xzOau2DaL06%~3yJKlfzKka>2R8wic_p_4#
z5k#Ya(uty@0v4(w!~{n{!GZ-4kTTLi+JF=pNW_^@k!Ar6Dr5vi1VIG`1T;9lBPd|#
zyddorhoaI98cLFV{@b&@yYJ$=`7X{{=f>p<+1WeKe)hlpN@f^<z?m(ARH^lrT)1c%
zkAQDW1hq9ey{&j`;Ew8d3Qr-k)i)S*W?1(HNTE`La{iYs&M!+%7B9}kC1<}|xt%QI
z+-laMwYZiRb_<7&LUPNRnvPhFsBiG`<nzEYxL+rVtW^0KAf~kR?vqX$U)qTV#hBib
zu~0gYGDShtpVJ`DH0V*5k;W2GTOh}f%xY%;%9#o8>nCqQeWobqwoH0S=-JKDl}hpe
zyVUCkmT))`bl|d&nczJgd|#=f8Qn{&LZf4KjK6da;tSb>OuI_jJ$>wzmj*VZw1!pl
z{3mf5pd_hIT_+C9o>QXJ6-W+QhS2+6RQ_wPgY%rz669ctX;Q+*@o+Wz9bo}@3t$l0
z+uTvNxi&V`evgzFq%nf~!|jP_R{gZ>T7nUZJXdEc6p@I306Wm5D?sb)SElMfqc?*>
zYN&td*jrNVAD~@BLo-mIv_8;$EvRF4H8BL2Rfkf?j3gF%JbjCWQX}?T&l-sAu6a9-
zo<~78EEfs$FJ#Itq^56?yAIP*8e+iy+~J}obkF$H1JlkbjP@o?#ylaE)xF}v*0bQ?
z0OAwCAeNE*U7u!vjx&s2$7SBZMc$ad)sTdd`_%wu%#`>uEOKglad@fiJE!M(@fDkt
zGb4;BRb0y<$MpxvpbW}3V66wHpoi;Tt2L}%jUl6tO(C2%2DT2qCwlfLhP7buO^#gM
zxiHfx-cPgF(jN!~qJ3|+VZv4vWBpIhsp=#bFe~A`n(9$ABw@V?5c8=#7HbL{dbY-%
zLV7g?aj&v^@wu~GeP6%qMi|X$>OL&f8J*QrtkNJkzOtp=Pj%1*RjPLrlSI05;5-Fu
z8?jhL3Cm@N<ms4K4x4Mfi{}@J29Q@c9XcBzbkcy}_Zob8w26qj=I{lb@#3`sF)R-A
zkQtPHF{2;+==7bP=`uhtwamrG=Kmpv1tqHK?^MDlXdy0Z0g$e(lt|IJG7_=`B-;N*
zWxI~WEOa(*;kFiQgOjp9^>na-fkB=2)^a4HJE9xUZL$UD56QGdwU47ovlqEz{yC}Q
zr|+2xDDSi!qst=u_7YyVV{{xZov1_n>-m9tp1;66Wb7Ed&a1khT$k($*-b4}8CAXC
zUD6P;3e~oww%Q$TKa4(Iwvn=r3)4_gqo=qnXqmEiE{q}o)vFi(l{{wJ-LFkJzLeV<
zHRz$~-GJsWSCC#UdsfX@G4{X&S3}dqf2Wl^tkj$EBX$K9#gy1hmgS>uyE(^c`>Z%M
zq2DX>6STT`*|tblw8swJK=EwW#}J)$V0-iM>e9CM*N_@~*K#=3-L3&rkSsY4#U*_*
z0a`I+N^R8Q2qSyXT3O2xK?S$wi2W<^Qy?pRcCofT&zVdJbL;~qCr*U|$iY}Ee%8Fz
z;H4hWJk~qjnwA&32a5*;p1wUfK2qS-Lv^GD^Cri|a#CHxsFEqW+91qv_^~S&O$vGS
zgZ;QwDH!DRW*E*LSmo)v?x{TG?zTriYHhQhF+Wk$y|+!;)wZ(;rQeco*djNp3`J<2
z9;y7f((QSGw#)m0Lg_;)KMU8iJdvPmJp2+V=JIn$SkV*G2LH}j(e1d`pUX5PVIzV~
zPvm)^614?Tu`02O>K&W35h{!xg$MG$#WdN*CWg;E2Y2#V1L-(ul_b>%B}(9+fGle!
z99V$s_g>$#3Kb_-1qyF)Y(djZ*(b=5Crw4P^mkj)a^YreIz`b}^vj>gNmsm|(4Q)F
z|7oG;3cp|(NV15MZn+%leL9lhr($m{2vsJTaMC>=)|8&7!`$D$_22)AYi7mu9cyAL
zn7WN<#^!Lnnr<endlup_Z_~_%dKT3M767*Oh9WS((OdX}6Dc)OZRn6|iAZc^MpnW8
zTbuS8rgK|_Ip@x`-YimN*FSxe#r8ea?`7R#qGyDibaliyo!W(p+GQexKR)zMyr!LK
zr}U7lN_=x`z#WRXiYtB{JA$qrHpJDaV-$`&^<uI(7ug3u<OPOjE(R4JC?c8$BMh{u
zpaH8Y1QrHy>PK$c<i|rp%DNU9MKA5?Vr_HYyzVBpOMVEKFDF||qnS3Vp+`n5G^xqB
z=r%rs8y9K5Qoqn+0FyU4t3$BrdmVHe3bOdozj_=gu3R`#eEjE%dBM%+-E$R^Z)Let
zkMG<8ZGy}RwN;Z5L#)3R`({q$7ox`LdN-(fOsKX{7QFk%)T7xF#+5~DJVy=BCe=IM
zV*hjB<DLIK63J36u@iFh_R8|Z#xac`skKUT07f;PSHnEdgUK52pf+k?vAjcS8^y+S
z<1%@k1}YXjIki8>v^&&+7cnx`6SaM(7D2AzoZ!kXKY66R0B>t$m%G0t0hzXXaVu31
z`k~&)XtDOOo2CM#X0TN{v>jT9i@TBOkDD}rSie~btSk<#K+^{f$JFg1Q$qU_H~cu%
zJmk;Iomk$|A8loS1@56<v!><9jJHsr{jz`yQSAg5qHIGJUygb=7x&l5H;a3F%TN7W
z;}72GZCr7G{3?nHJkAqW5BwFhHk6FyvBTSN7}2DxJGva83U?v2f~ParW^Jw-tdWPW
zba3G$LM;~{O;22HT6=0;Sv$6YW7TQ6CsjA7r3ho)viX)*$RJRT6+%(+pNxq2+GsjS
z`*PGTM_mJIKr@y-JTY1H-4JZY?lRw_0ona6kE$6*uDJ(@0himn<sF8OX09P)T^Ijd
z^FzYPwb6}bzm=FB1R}5ZPtC`s5aR3t64;L2D)-u!V4Zc3K7Q35jj9-r>-~GPvG9;m
zF^Q8rZ!^;Ou!HaOrBVvex5P?wX>w_!mcHD^whoAbb~#GPIt!kVpMlvBnCj=F&i-xQ
zHMo()TGy$$sbu&-KU_Q1CUsuY+cO>u|GA)ZoHk@c^3viI><EL7T4YBnhsJM*t+91P
zQn(;%{<Gm366>4Tdk|bWjJ{;?Wq{~ypF27$5e?$sgj7<O;+Z@%0f{p@SF@-~bs-X2
zzTQjozvujVGZ~p%&`rixAkAk9<e`~AM@9Pxf%g4z4{PSfnyKfD2X!fZbYL#Xg6#s<
zhno)9kLjER!?UpK4^ryAZgLf)d=Hun#cKBE{o<!))~P@kZt}%UXctqZh4X2@v&R9D
z_E`g*?DQn;nf|ugGis}rMv81`<zUIS<IvF!TP4B1|GLxf61|G_L;%Cew%SEKVZlU`
zra%ef>}LJ&bt@P55my2oGXsN82h6oP`d3@EoPzb+5&nuD(7BNf-0t%)^@M|<AL${1
zo>TaV2>zr!N$A<t-gV8Av~g+c9Djb~)>id3EZAZU(Tu}Qh>8gNzLC<70QveyR>y2E
zN<TOp5it&(*Mj=MHz#jY{MU}d>TGELlh`32*cv!tPO%50XTb!Xl@4=^m}Flf;aP68
z$5v)xj<@&~>>AM--{zV(bAG_Y+h{}L?6Y!(ko$&U95xU@_!b$?K0}`0YSw;{SIwg4
z_xoYFJ8S<Q>n;bN#l<@O@InEVZ33O0E=s>s=|Lp6^r()@1}Q~^k5E+yQ-;$BG0js2
zvtC25rqFZ~#(bdE(N#8G06R8Dh9~!rQA4%%+otTdEq=T;3R|;Y+qd61D$PLY(4y=A
zM4*LbFf(=xjfb)I@Kp7p8udHyxErc?W8s)K`D<2!4DvRXEmKnn+|pyt23h4e#M5cS
zah%V=IN_V<n>pnlwbrRNpxx0-vehl3Iofz_lDR?_Le}36><jZYMlz1B4Bt<Ma~_Jw
zo!XL#+WMtct^Dg!0WsP@h~^DXef_wR%mVI!9zi_?mQ#^~{?ZK?@DK9(aX<#ak7zDU
zLjrSsi;>TdUy9cd^w`ukf0}4YFIlWiu-w#?CXH};`%77$zzfX@9Vc_jXCOzxV7ZqH
z+Dv_jk5lJ(vJfBLoG!2>78Y$erhe)ch#vUGVPvb?#Oqf%rj8*Sm;voM?TFPoFL0Fz
zf=TQ8l`oP}Fn->aR<&zzGJ%??4*i8$YvlClED$9+fv^YeL+o$!rH<t(F5w<s_U2$R
z1vN|v;%VR~a@hf7o&(wjw7Sv0ejO`&&<C~LvlRk;#DSKsC1I)|UmLuudIkaoqSq{G
zkw>bdICPCE?5g)nP~EnJeT(-R7|%Qd4!OUQUlh!@)H~r)@uvF=A6D?qK8lLR1zG&n
zsQ5_ODwOvvVj}l%wmWJGq2xi9h$9sPw;%BccKc?iw5aH_5xekCk-}p3T@^xJ=#Fke
zenjc=sNxx(Ph#nl?O}M2ex*CY>5k+Bl0sLiv02o*$M*_Ff*^r9-)8}g8J<|<16otH
zsO}+ZvBtNfWc;i(?xM)?W1>be!w%o}<$`+&)DZV!)($xA5q5_?q}u{6I5-tIhpl}#
zU|YXFMG7`&owzb3?<spLYAPX3hlt5Q;f3Q)cZU!{jUpR~zS+bjV17KiTiX4T*bRAV
z=?3459OXMhu^y^|pp5WZxVd<Nw-t81ql4-$(}c!yW~bOG_02=kZNtTVF12r_jp=kX
z4=DJ&pzGNoD#cjdvmPKeX=O&c5y6MBraayL;EFWqM1VBI`WYlm+3m+@m2akV+YT_}
zmiZx07Jh<qUq$MJp^@cQpkLu^cHw@!@xC5hEAPUN{kep1WC~TZ)(9J9#=d^stkIL8
zzVn|2NhAq{Ofn;kdB6bdJ{_Gyv66%A-so`WOB#s<8mb~3w75TPh6khbcm7#-X5dHn
z{UfU+maq5_ar{@&UFMmt5+KJ-)Dp2il0)DT`pq7fTOsY!8NPomu<QG3Wk_%<1HJ(+
zy=72n;2~6Hd(TcW@a`aBd^eL$jCcUN@271D!&6E;Ov^F@_P!W|2|bM`i3fq>O};M+
z27_CK0#aRmGQsm8#fqNhjn>tO)5s`jfLT3%7&MWD(_^>jZcSPu&$nW`Cs&Us)S*$q
zCryD<U^%jt8)`~vx+l?BWL_;1s7sTUvLYjcV8#?OCEUFoZ=IBmnz99N*<~o*5Vu2!
z`Z2W5=%FkE$2srM15I^J4IU)a7^A-ioT;1y!%nIR-jZs+L+^|iDwpQC%gvsItX1%l
zANXmor-OMwCo!Smo2<dL?rnaTQJEdaID=U<;mWa1{GWg2|21n)79vhJM8-Nd1$_*Z
zl{FZz=5QkKB|N+9lWzKg)^miRWE_C+`-m$$AeK!QbV_SS;5J;@#F2*wxb7N*%dqK}
ztfXjxD97N)J^w?-3*d+KguKeKcpF}T>c)SIcQ}U1MnW+Tcxhb9zP#@%fN5}%6zu;T
zc3rR9ci`H9PR{5ZSs22r1mNw>IjOdO&?91*m(s=q0*?bQLZBn*Z;cttm63_CTV8s9
zJ}y(@JY=E(91?q_@6Y|ljQR$KZwU91EH*HXDQMO@U*1q40hY-M>cv#=+~0j)8!r@~
zAvdW_3Lo;psKQr7o3(&QcN`!Eu68=2XH`sRF%b_sHd@^?k1lbUzH#|KY!iu&GoEA0
z-ZTO6M2(<T+Ox^F8HQh%Qug5YxRy;GvZ%R5AFLfAOVAGbI2kl*3Z~*y$76@0@M{ln
zr&O`NHQp6_9!&OD6Q^5Y{P1@)QSmTpFH9c@10od~pivzij<=Z63Fbmk6E4VN=Z%H_
z;C1&-&18an;xgEeZ`N~zqk{S0T=a+Tz7a#>aaRI5NNsSlve@VCrpmu|lfa*adF@+T
z7T}2`Q&P#4L>RNBbnu0E%~&+JygpqwoZZt`Z`@ld%h}H5%OXd0HdtuL0B^K*;sf%$
zFh4l@t9Fs!yhn$aaZBQ#GweS%EE4aE?w2T17*Kj_><&Q<Tl8-m;%__f2i=u`ofNxM
z{|%&q#}3CAUB|<AR_sp82$6!o&A(FAX)VpZBfVcfH6Jrj4N3d#eRwGMIDjC?f3m{Q
zh_nGnXI>yatmA%Nh8{J+#;Q90t;x@uk(tFef!jn+krs|S1<yw0+?te<6WsygGB_fG
z{AC|LE19in3c6$R#ba~Gbubz;22gGA7(hUitGS}lxBH_9F{`aWam(~qK`~1$PDzg{
zzR+lfz+?OynBPH;tQ}$JxbDoq%jIWE({A!9Xd9av2;FF=Kjd(*^I9$t4){uq0T$vo
zTJG?Y)&!yA3{so9=H9eB#_VN8^In8A{!ZvU%g&Rm8H0Y<^HS|SsDf|5V7pZC9nwcO
z9~Tyr(v$zJ9)p44iXCce=y9o*?^(7+2`=LLiMj1wVi&W#8Ucs58dAF5aYd9J84sEb
zs%4M3?3F;xRTkaVmD+;FkYl2iZW{M?dVHpRJsY$n9dNbTk$lrHa)@PrnOM=Rd0&k8
zMr?{W^TI_u$6E?hf&*lGD)bK~XI7fB7lc)xw^nzgSx6^$krH`QS%cOlaqR{rT72+L
z_JG)r&g-2s9%t)+W$2$xDE7*?ldKcoh>{3De~}Yj;1RP|yMHwo`XQE`#FZ;kbbo5|
zNRDZ_0l14voYC~5CHum7jDekRTOZ}mCscn(dpB1f;m1^EQV6o?iMt@3=puRBNOE93
zbDf)Mn|BN)?>nEk)WExW&?w<-tv0Q#*Y@()^4(0qL3WjfxDGb!V@Ogf!$&MEWIR)*
z2Q$CHfYa>G9WBKO-|L1z2?o4x@g{>yd-$+yq;#{StWA3Ku@1$GP|Y=Z=MR-d%mG(z
zr1V8ydlUcEhn)TAZRBPAmJR@^&O9`&eDvpfAGPPlhy3!z3=6z9LRV{1;(SIVL!%I<
zHEM=ZYguu3sw%#L%U<`j+2L=yMarHge&Abv+w{}!_X(vHX<5?EjkA2oH`ms=iWIhG
zF8d~vqE!mpQ^~adbA+E6CFzjd?ZAbea(3ifS0ibIwUvGn4xEi`Y(@EEuc6FY=X>lR
zF~59r&=zkWR%PZ+?m~%l2UxBL_h*hZ`!2C`Mq4ymg1;O=Ck#JAkhM*7a*{$D6}0!f
z--oHf&jr5=z|A3_-8Po9XjubJp_D?n_L(oKpL37Nwia@(-IIM1G+)tL>NmiDDK5H)
z49<)Gbh8!NKr83NJ|eXFsFSD-D4eBjj|yd_!Jf<5;dAF7r2KMMjWkuu_Ze&35(O4U
zninec=z()rXW~W{Dy~R_L9@e5lG4WrO0&&Jvoq?(HA@^PR34G3#j&wzOWOE?ycXk;
z(qKOeGJ;K9k)|<Ww|!#fxX4a7USP)0ocHM7BN8bk(gD7k_TO8E+-PCyVDm?g0=X*5
z@(cobN^d2ztbA*b)Uj6we3#R{%TKO6@&s-99Tx3jlgA%pzraWs<`Ac)KfY5B(D7S1
zBD@hY-k3jM5wT6X%&{okL#cnzhMIU17NmX*Di0GC?A+;0z*>Rv({a2^CSrvHnbk2@
z*kkv*PXDWJBdPxt9k={&4~_o8AE}9baF|U1yxEtmhtZ9eZ{fy{n3v@yfWVjtkjGFu
zu!E$!u)mgrktMO^bPB<&`;v+rrcAlJs2IsYG92#Os$`R68B0rcDs1zTy8JZw`M!s7
zf<P0V+5$V)RM7k9+130w8s#Z0WLgFVTS@g+kkijpo`l5Sq@T>Z>)?!v>JA&lH_z$r
zjlW`^5d5#LDGTFYx4GN(4_f^9Mb>{B&0$LqY62J1{_$N;n}eL|Qhu|S@4&9L!NIjZ
z1ROtYMSv;8ms97SWr&B>6Po9q)zgC;Ow{q1=9Zv`i#4)V_8?Y?`E3Vvx=ZGB-${Q$
zPRsW$))rQdynKdOipO`pq9X6i83)3%Art5!^s#VQo}fG*%m_U?r3W`fB~Iw%FcYW9
zbWw)ft(tOtXZD+vUmkSLIHfF1)p$g)`E;s0l+>kHugly=oI7!peL|>$As(o`>|Q#c
zd)TKI3i3@!O=lAVYTRFoEAH+JDcAY$3hX-U5<LeWY;7aC2&;3!bqwP}<=I^yfR6db
zLggnFN$YU(7ExL6!9RSu#Q7xV5|l25cm?sawEYmeNR;~4!#5D`0VlL#J_zwPFDaIV
zQAv5#R-?fw?&yuESu$d^1VN?KU3f~JM@F?Zc<hb3*aL$WE8Ul%EAfOg#tV+DTy|_i
z@N*C2_~1)9j^q}JTO!>bpthL&Dg(}+<}D?&usV^$t6j)YulcL@<=Uc#)5_kR#{Rq&
z=V*t0_<nhU{bZE9ux2b5&gJxRKMZQl>7(nYPOrCd;FZ}<8*OJrdl1kZ$t@T<a`xFi
zAHL(zC~z+KLACJx>nLwm9Ktwb66^s;GBMRrDPwPf@yOKEo)quN0>qJ{m8(#%=-~v(
zfYnk@EI=m}X1NW0n*OWMI=U)WS0i<^75b$x!For$2o%^|U0L0glRtX)S8WKC9{rDp
zNF<lY3xM%Jr#eQDD*!Wvv;%S{(6+a+uA6Wm$2%&^p%rDKM6}%|XMd~HYp-vuCit0W
zt@vPDl5|YW7DdY@@{Hj4j0G^Wm!@NnCKVp-{2Gd(Xfi!`62Dq6hFOA!qb|S*i&@+t
zzDT*Z8F9{>*A3b;<@UUu%Z#3I%PE?TgAES0=eF=*bdLk#v>D@k?0`z8ZzU^ea?s*S
zWfSh`{XyU`VAnY}up)u*(PF1TV?3)4=RQ;21T>QEijRUmdF-tJBm<}`bk7BCxY0at
z=@+8kH|4w)*>r!7YoE*5B$XddQ=SdqB>wz$2|UT&kmgFS$)KZ~>-11d_EZn)W*17<
z{k(!yWzOFDakPp3>m{!F{GvE0O@~=f{5VM612>CWU<ZEd9xR`i*b4bF!xwl?#HVyD
zK|`_6E2-&#AehuRh&QM##*U7e=Vm6yU;kF=+EEo6YcuW8o$NcOzf|MZaw{Fk9JNRL
zLGGmUNFLin38Q-+lnZNISSOs4eqQ&`ND7WKW4o5qz9-a}d~4JqkF7@yt5Bc-ig#J;
zD1%d>KzC{xJ=o57zL|T0?hhzO%)>3ieEh37_9O^r02N5KUv*z2C9IO1`QS^%+Fc=x
zUqpVg;A?O7-73U6T5i@CIoi6a(c#q2n7mA%rZn~XPB5vz3@!v08;liMDX|-r!)UB2
z#chemNrUxAoy*$oj*Y~R0j!75*z=)q$;AOB0OVWZ$bim(HfM1_!Uae9c*09!#e^(p
zO!#yGg-pdRp@wYa@vlWa{Uo&nKYn5|2woeD4>+GZf}-oR)Hn^F!iy+*>e#ZSc?W~5
znajc4e>E&@3)5th76VGrq9Hd};e(1`P?G+H(v}+4pCbcTsz8`)=zxsd7{&-PD$$x9
zQnQegWk0Wmh)Hh%AqUUus4VYe)Qdi^b<ToL`YACp`gWh<#cBZqjwU42{xmRIs*$m!
z(K5W}T1Df#Zos!oe*k6!dqA@)19Zi~_(`~AI?!s+ig3{q7&uy}48e~h%vGyLM&=m9
zp9sv01%SZ`XHHT8r&U}pQp;VmB#iiv!7{DesDPZQRXOz#9;*iVi+!*_R_(+ez{>&b
z=VoG|5>5QRkJ9No6<IgPsaY7>a-Pn`J5wixB>0>RSc>$!;n8)jD=Eo6oaFIAzRtm%
zEc~mqf0P2h6QKv40{urAywh1mFqHI^!WoDy##$C$M7Yy*-S@D=VK4+M53nQW_Isf8
z>OkY7d%SOiJy-yR2z5cVVmtcz4bX5_S_MbsgwFya8(#G%UVWgB#)ilGm`9=Lu+jmn
ziiI^Jw0*B<z7#{mOGWUK%o9h-i8Wki+xIsSuk3BhluflRRm11bs!Cwf#FxVossNnQ
zC=NZNLxvxWDBlZLBflQ2n`1Aix*ThXUkabrvI@?rw+gJT0y$Z`nsgHBV`pHs%J=m}
zJ<poY9a|1BI1n%F?-19o|LyzYbwg5av4$cip9voj1pC*&?9mJ22Qp?mF+~n>|NQxI
zcH*A2B*5pxH%{b9;8Wl&6stuz0Nz4P{04ce<y`eQ%G}p(558FJfq(Uwm^cT2c>`}5
zW4Qt9MjSE=N%8-3D0~!+G678m96+b`5_f=jsoVp;Ch^k$^ELcW{-6*E@e>_`-<RO*
zjpqv*r0>5m(+#uW@Bdj11Lgny&;P$gc!kTf<&Vmgx?Z{kqfxYf@1Z><yFC;C3pd#V
A*#H0l

literal 0
HcmV?d00001

diff --git a/tests/python_tests/images/support/vrt_colortable.png b/tests/python_tests/images/support/vrt_colortable.png
new file mode 100644
index 0000000000000000000000000000000000000000..48afae7fbf293ec94788e2d338fb4475db54c93a
GIT binary patch
literal 7308
zcmeHsX;f3mvv}VG2ufI76V{-BA}X6GqL2gyB<P6Zf(jTEM^HvUk)0&hZP-^qTnWRX
zq9QnofCxbu+yEhhqU?%505L3K%lbOb`Tx)F%lq^`y>s4|d#kIe>sD7+S69`I+qH9>
z>U^#F0032*vx6G|7zkp3vJ!j=!dry^=<K06*zXO$Ioj26=)2RWhr3Oqb!XqtxqpaV
zxW#33xs}1u<QJ1qJ3=+zc=1*nXHm!qdX}6^K|pT=6k>P`0P2?i9YDqzAm+$qO5mER
z4Iq=HZNohxCKhe@vvQ$yriVpu#MQLvU#<F$r3(HYmRfG~8f&|5I?>{;P;DwBu{qaw
zut>m?Cybx@`HdqM)vL(v_ov*Dtp~(|W!T=XI@u0>HvuG?NWjKT3xT&`mnyo_Obktv
zH+&p8c+MsP;PkN$#8%I0BOvki7es^csL32SJ^u~pwEO_!an-#kn-t{3n=RiVTTPZ0
z#?{NSDLOGNMiso)1WBbKES(DW&+aVno!baL3r8XYVu5+`3S7;0f#R+|F*R@SKvMzC
zTfZ4x$u;}K5(At*G=Qt7I{t(-p2~GRK)Dyy(NdMypV;@=PxU&|h-7M0yehgH*W2yr
zGg5rw!J)rrBBn#5z|&4G5T?fAB_y}6X|HU$U)F&cdtZPA4S>YzoVz_IH%~@5Fn{{<
zCksW&ME8$~e&)UWZ1PO32>0geL>b66Q+DHj3{P^56MwE*e@BTbIJYXPbx(2e^=#or
zYQ%l4cceuV{WCYF!0OCon^#hshpKBdJ}0{o(LeFPC49IKOi45WHzp(h*{#!461pV3
zCU`%tMzD=mac{hx&0af6!5x)o0e`dov}Y~d*D(Lg%GodGKS!F@0mDl;TfRSV;LIBJ
ztA%GH`m1bvBdQhkS6S)*2?y{Qh&azziJe@)0M73^I@{^k;km0E-T*5j|N7cfl8~Ha
zb|BK=NW!PT3Oq4Vj174pZ3FaNbsWv6eRD9nNj!)et?xN&iuBLUNAn+I-hI{qW?Zwg
zv`9kGEkwUibbIpTs2KBUZYi4o6)W9j1#VO5qNw!N$QpWO59WZ#KfQI2f#y```13NX
z_bQ>|8d{jvi#KISzPfsVYvxHa`Kam3iau2kNQ%(af|&g@ldjC*V#U%x>j)x9|M1kd
zC;4U=*QG1K92gF5Sr91`pQ!%{Hcl8`m+emfd+}f22fN!(^&QMosu8r!FmD<rb?^3N
ziN6HfOg!b@Y2hpxyGB|~(`+Xnp)gCQo1RCJTIfAO-<JV|t>L20$WB&glqEEu;i|uj
zRm4q&!ttkU*6quamrv_Tm)3OmWaU2?^%87fa;AfL1<|KjcT@(EBK#wi&2R64qObD^
z*@li`{pVS}%C<}qF<s8GI%GUBM8W*_cC5EFKbHBG=C;?!Oz0a^lUH_m*)2U#ztTaV
zvMYW0CdOBKWpa-WCZ)Uk#j$5a9;Pf6ePP{7+HPIpHQ6<KTiXkF%zIBlvbfwis%)xH
z^kG)$BN#l?QWS{%VMoiQ-%Q_oe`^n)GZ8qqUYv<8Y`fH-l62<#_$|tcG`9aKM-FD}
zC!h8~xk?wfn~GPyqP}JA8*8A14QwhAou;Iyqj{77KI7&I>BgNiukV_wqgkc6b;7z;
zwDlc209V||QR|KG(JF5r7(HwwGwWh)Zp$bVMf2HqlhW>~$$r@vP(JY}D##x+IC6}Q
z3Xgjb8mbaxh1j`Q?l>3Ki`#zI-x_%{5qwZG`LeBGu5WJYK(r}&z*k)CW!wl-^RZsK
zz^UdPBcfO9%<l?odS)_ThcZd@3SGO3^<etvRUYAiVdjsfOR`s4vL3w#_lF&qw%JdK
zCTdEL`78X<{b(h5rphAhyF~CF3cFBrSnQgN=H>Mci#D&DzZ_R{O%%oC1SXfL&_}~(
zHfn)DoeG~^wr%>}dyOhi=N60|fCdpr?79Xk2|rY@zB8vqog9flxfT3$B5!4qVP-%q
zGIm@+AlQe}XJ+0m2;y0I;L*8|Q$m?Rxp)b>p6j}pX7c5$nf8no*zjk~C-p=o{z^)1
zTEa}xU(<6lp(*#?fLP`kXSbrAs;GfM?|Y9L<=ARjR?3S-GjGMADGg9<>mo0son9sE
z@cSDvxqMma>zsGsKqrnuit;Xlq2%(G^pj;Fb$&mKdh3cyk&<Ygcy%->LDu$K_bPO?
zz{s%v&ykxq6CWO~P)j*?nZ=n}P+|?Fo;=gTVix&l4j#}VTet~hT$Rvk?7#JhSOQC_
zvle{TMp&k20|O@0I&4g^)bvo3kTlb|>iA?-4Odg#5fus^bbi`OVwvB;QKxXTf2A1E
zeA}Pp@Fi7MGw#E|t86e2JF{b7e|<mtJhIa1Lo0w`%=yEYZZGp(vs6bXY7AX*WAtY0
zfVP4ZN@_D=Ud3rAZA*9Jk17-Wz?#63j*hv34PMgE*ob4mbAd$Ra(?$c8*oSqi3#w%
zBwwALb0)HDt6<U2^hn-B8M43okS<c$qjdM|xJ>vNYXJ;TluR>2ZNoyr6LU-`t8r^V
zf~@H^$zdGXwf)zu6O)*OK?)8j7fL(I#IjOy@$SXipwxkrYFJp7EIwaT@tCtTrw%EZ
zmiaSw<{4XplX1)iG`HeRE#U4%9ktMq*d%_N26v5wb$P*FwxsHKX=lOi)$73_ZDew$
zN_@aIddbnnSvcx%T(j#2>cTH9Q(>^#mjQoC!wq|VyAPj~uWp_^6ad{iuty-f-^0dM
zSb>uX%ojS0-Zp7>C1js8`#74wS)gzCDXkdOss~Crx=@|_jOXWW>H*h{B5&L2Ezf2v
z_+1ZI2G)Lz7eA1F*N?4kgvVJMxNobuR(7X);n*97%~g0O)?dRN<J%ASzer)8;nagv
z;g1PmJiCcC;epxrv(`Su27F@vxHZ%E^0K~g9y9+QfrX{s`_K)ps-NpB{=QNNF&GH|
z8PRM$`F)G)dEUx7?maKDDkyZ>hC>&WoR+CRxR2}-GO9%fTq$E**NJtSz>QL*hFTls
zgl{(RF5-!P3|4D<KF0Raiflh~U8mRCS%5@PxWl)4{@6dt;noVd{QJeC)$_19%ebz>
zbxUOX(FD%kqB)h-$AAtI;m&ZI7@4t1b1HA8l$%^#?j<hRa7Y)ZCZRmh;2_iG74UNq
zERyXgP`||MFi7ei0LE6B?l0OHe~SKhvdbGe898l5U=?yss-bJQapr5xEkS?aQ6~qI
z#<qKmN7nQ6H^I$~>XuG!<|y?Ju_Zph%0P&I9*JWQiQb`m@ARuoyi@+X=i9eHU|au7
z-Xr(X1X=rQxzU|7iO8My@L+CW4ALCSt36WOPGHeOI&^`N6V-<xrkAlPc=S*BF}@|p
zmu|MgJni0jSnN?gHkpkyUk!6@=hfE2eJOsY95l{PN1Cs}@g@3jWf5WYOIgv(BBp3#
z1K{lvs5I0;>|+COz9WedRUNSCjwZMj&wNxdyWPfr^rdXd=5!p={D?<@mv4VV9Sg5S
zeVt1WPhY%b<`?5Q&CMah2b?t)DjL%k?7o6D|HUKxWoW0{56RHST_D*PoHhKH3>TC3
zFL{axbmN%ITt2kyPC=UE{zZmp0@Va26nU)mcVK1u?;JMLe3wUH8#?|I49)AvQ8pTi
zi-pEqA*?f|mFB4`Q4MYr90e+S(oZ@>;ZY?Aj=53c{wMW7SwnIjW*uA{w)Y6vYzNIm
zmthI?mP)o{<}c(@u)X$Nvjene+I{9g4{~PaS6b*{dl7k#ut~dbogBI4N`78i3^IoD
z2u}=k<2>=`yIgv6`je$Zu*8v5V{}4u2Oiy$OK(koLTJzgijJIeqZ6wSSc3WS%yZj6
z>@PV67B7`7&&rP|sRMWQg^FghJ^vpn?$yQSoaYg$3}5JBbJBSPPJX8x%H?6HZtr`Q
z<+xvvu^ie;b*t{h%b^QYw}-uJ<WMby1_sr1j)1)j8ARuIx<4#{NgQ*0`@-iw5G8q;
zmA~r1FX;S#Lhf=EjD?rXRt1*426qjGbr1@=i0q;XU<&6t8H`J|1SYX$ykWocPjfl6
z*f6Q<Qv}u*F57{mmi}&PZhU;7U~qgS_>hh7zq(+zxLqSx_Z-qhcE!sjEn7CkEKZQl
z5J0sPb#@S|s0lJ1NQ={dTC1X2ORMJtBTc!mL+S4X>jAlA&IDx0jm7Y$nY5h&4opSh
zvJT8oAeR_s`ahf4Cfh>2`$BHQ=JdS}Pk`AklsvfBSUGv5#|+HIiga=?a&uz0yH(`k
zK<iPtP5fwD26y3N4Ta6AR!+MTkX`nIbIS^4IUBG!IXvvvM^)9o=}2t+ZvH>S?>8jv
z#r2LoN>iG)FF^T;$VW~j;xqTSgeoM*?ZqNInvkomoLOH-U#l2)GeS)h5a0=X)?Jt|
zn`v!8wLV{Ydytzpn1KY@3zq%Wx3B6m-&4#8q57WByV7){Xt1Sw@*ODG7v?lG3(SQ4
zZtlkg-pYP@KUPV^mcTTYY%Hk$KpCM#kI!TRFXd$8aa#{^4<1#?RmY!lpQsndm`+=B
zj5HXmRU0<mdlmd1lp70sn;+ei<#I4S`a-gsXOPp-;4LI*j{wP?dTBNr_$O7s$O}Ec
z&V|CqqnKRv1(Jtx>~U?-KJUGWq!U^=j=W><ps&Q2bx<~OlwU7xk3aoHP@Wwj)dZIu
zNGlGmo_j2LOb~U^4S57`(%{vH_P$6#gdbSEkdc$v+&IFA!RRr4!WvjP$|9pYd+Z3R
z@I<Iz>5dY9V}$_E$aO$K?;y<U9+6md2|Irgvkpg1oR|H~rAr!>ok<L>Aqy|&PbKOy
zZp?O9-Re|a^gi!E<4)=Pb4cn2Vft5#4PbsO`9S1Iv*WYc1leB;dhdf(d8BK4zSta;
zSNnTL{9%DvGU`TZR&KTx1`R!u@MbLA`M%w$a^gKTG={UV*8{n}t;5vr_T#pQVxgo#
zdm*Dd(QOI|I-=mIi%~*v`*rlcp7NPK)LoJBR2Xdz0^-S@PQz))ShKP@5X3+AVtP{c
zAFkKVYPrJJR;p3u#?<?r@^a-%CkvcvuGxEo*l91uAA`kxSwkUWAv-(UK&i$NGR1aP
z%(5fckf~i`qxd+$>yB0m%HPumR}a3uhBRZaz2~8pttY(}`=_9D#EnzO74X_gqK#pX
zq!FH3#p4MC+5UfJ;TmE;KG!E(VDO+S_jCfG`|InF(bbH{2^=I9UPkLqmYT(uvb<`B
ztg+s2{94-mqepTX-w6vFVWCW2&c!v{wL5KZJN9e?IBV1XWLmg6a=d$!u~T?5QSQZ-
zdDE?RaIxV>@kEevcCdsjdF1#0Q<HeRXG0VV@~m_u9XFWyeyJpGfj6-IuqVd$)WO3)
zJj!RwrfoHuK&NUlg}jj)^W}{lW0~!G%t<_g?M=g@7nc%d?<O!tnl7KkQm=2a11{?D
zoJP7zcBF(lPzEC-v+)#gIG1#`l>UZj4u&ZYAD3M>K*L-5G$c1ihHQ=k=N>yU29F9;
z3X#<JnXikcCrkFNQi5f<eZGas2DSan9fDKZk5(#y(F$Ss$4I@WCXKM*Cns7R&Wb^i
zT~%$jHgRU#h~Npop>JP*#OIQ|*j^*yt`AqcB-8eSLCR6ipZg1nL38!G6!)aTBlRC(
zfZ|=C`5XsNTAG@qs(|bxThG4&GmuMo^m8^5B;I)!{OHIkd1nAA(9r}IYjG%B672AF
z(77EkCom42_~9vLpJI-)3M~N9{q`VD2>THSR&6ttHvJBQ+^-D6#G~e%d}`A%u&cwk
zNgbPW4!b?#%y@B6F-s<!>;dK&<KHs|FUt9h3jZjTtp&udHb;O}&;gps4WwpKB#TT0
zsyh4A-#y&`;;UIist_4l$bi%iaT_xG8|r-v89M=-^VUmRIRmVD8j-#r`d2h<1n3MN
z4dljT=1+VSyCR(P*_xmySKWlxe?#_4WnVYtYy+5G%GgCSVIx#pp50HSnjL~8P>c0u
zQVQvUl~YuMDtbvCXj%3tcvO{(i_d4C*V6>rx$28)48<x7umMJnbYUE;bEU)_Y*<6@
z<l++ZLw#<_90FD-)zG-Or2Of34U85uY;QjAf|J{g-}|bE*z%}pZ%|ledAcVN245H^
z(^?N7ET_LZ0z5Ss&CVu!hI)+5u`_Qs;?c{w>T78I*QJ+Vv}lxzJaE2WtBDtope=Gm
zJzMIS(*55kdR8Vgwc&Vl(Jx+M#cWfth&=-vQ6jxv(;bxM*xnpot&7_isi)Z2p!cUp
zZ}&h4*#z5*@-DczRWPdslu8~)ES7DN&HI{mtg`AL5#St12h!g?5ogOLzR`z00EHB&
ze=xuoE(+nsT**Hp(L%I)WM*^4xa2jpr*y6%B?iKh^u15#ilyw7b6Nx$J7u1#B?xsO
z$#+g1CiSO_ttslt$H9YGa*m;0+hu};Lej=)0u^pdTK>(cfQuceNr?F(c{jy_?^Le(
zg$Or`E>1twR=FTx>jtHoI|-cp(Ba?sY}WB19959ul|C@#2WO>0-%&Z79{KdSszKM&
zG?-aEA2aHI80KGgB%@=4CG<F$aNM-d7syY8Mrwvx>=h<(AYsxc6xW&qt&MoYB+`==
zQy`YC0rP-E8KID<oIatj))E*No@rhV%QzFIA{QHwtFA!v!Pz<6z+<egO!HB!qREl-
zX?+S+ki3G6i^*qR=+*@3XDx~L!*VxGbsk<#bo}M9InKi?iI;x)>ZYZM|HwVB#Q7wJ
z7>A6l6?&(&2AD6$dan?sLRT<1!Fn$fmZY`rG`EyPHU9xQjZF3#h=a?;6}2^ZbO~SW
z0uM_YE@->NG7$z>L~DC{E7imoZPgc=r$_oRzp}d%fnSQ;U+fL1n#C_EBv6$Mh?)Ol
zO1md~rF2&NI6$G@aE>GgH+9$%YYE{zMu6N497!$^Q!yL10q8H1yU_*O@G>kA9sxRw
z8RoR2nZYPo0*uSxlm*zoM*Tc*EA3>}N4j{>SJrHW^*%C6jU`{r|KuU5V?9PRBy+mY
zR8td?u{pl4n2Sqo_4nWidjJ|*06SK=WUDPKQ|T?K5*0dv{?LAZ&bVtJXO=$_O0`r%
z^F69Ww`7qOv#c&P=M|%HxBkfAt@U+A8j}cX@vIvx8B6*)TL~Qd9XLPh<g*5w0Ggez
z^GmLa5)iEJ`ZXo(%y<m50Jg(olC0cX9GMa0DbnX^ENQCebr`xnI~%O$G@Ef$H6#_f
z)vBqU@etYUoVP({f+Bzcqd&Yq-2RNLoMkcIr*<}D^@V8K->rRhETO@#5ARDM$cCcP
z2B=aYM*140W~}ICAA#N~Kkg_CqEBF?ngLLZd3t+5#9``wxUvlUQ1}$sJHcnGy&Q8E
z!K%U7T<iE@^0#71QLy*i%(BfeMYb-gdVkIzN@)xggs9!L{w|OI<%VABcbh<J>)x7m
za*~2A{BEm5&0{$DO#xP1sf*hMk&4;SF!Sfd!e^wdKI~@>Jwiia8&GeaUjiJYp-=If
z<$i5Tlc_+5(Kq@z*_Q`ba6RKB=$!K{*d6I5xHG^y-0y=lYDI$Qopu*<)%^?4)~hfZ
zi!oqSiwQC`7~b2B4oMq0!q?dA$HC3TjM(k&QNysddV~(Z6w;A2;8u@}u9N2+eg#Xh
z?Y{Ets@G+$;*%y|Y%FZfJfzIKfb2RT_`2E=CxB_Np-|&WE=K4L_AX{<Y+rydU<qR`
zlDNNA_boh5g@s`Z5qYQ4tPX1&CDcicmw=3!8&LC<1a%<%Dw4AsD?ADm%z@QWDjko4
z>iIdlF~Z-0f(5YJLgnC5eI+zfop%G#_@@kMjiSN~l<I^=s>w0%YT(ZVUWR#Ruy$cq
zKil<vH7SP$eKV!)Mlp!)LOcSMsOLmsb(I`uzXX8=*YX{Ns_$6uXIM}1>Mjo9V>p`%
zt298#vcdN9-6dY)`v1Uw!PSW4IpTw?#st}Sp5nu4$2s0$NPlp}eSrwu?gJFw+5)Dj
ziMsSCXYDrA01wJQVZ!*A6!&4DXJfUW<(gCH5+a_gBsXRANjYnwUUTp&eMJPW!95_1
z2+qM@6A}wyZu@qS0dG@7;kWik0Lxb^W8s#LlfR%hWeT2R<m5V_bqtixVq1>`t8aax
zzs3R;%XwH>6!^UnhbbpkgKapfmlAdB09`S*8dT0=6B~hFLo6LR(DzbZ1J?5m&RBSa
zVb~px0H7I%&9S$WhG78Qi{_le@Xv+B5@>yz03Zh%tQNK&fS?KZ9k}141myPY8{lmA
zR`KbrU8m%&oII)a#1cTQUVU}uI|x9~IF6X$Q)RFl+dI-1c@N<03cJ%{k=}3W<||f4
zZOety`%XU(8+ePDtEz!lec>dBq6Cha44e3Nh@85RAhn^0aKtPFrg_0O0DLTgk`Dl|
z05sPpu)s7CU@HYGGFv0y7Sz7-@6Z=v<G43(<Jb)o-}vrA?Dy}pu5SXzuxR~%^M|Tr
zcw|Oltav`a*UxLZm%_I;0Ny)xh5_d{0H7D~;QxUCff0Bz&}3UQ^__`Jgnt<U&2gv0
J!!6$D{|joF+Y|r*

literal 0
HcmV?d00001

diff --git a/tests/python_tests/raster_colormapped_test.py b/tests/python_tests/raster_colormapped_test.py
new file mode 100644
index 000000000..c700cd5f4
--- /dev/null
+++ b/tests/python_tests/raster_colormapped_test.py
@@ -0,0 +1,38 @@
+#coding=utf8
+import os
+import mapnik
+from utilities import execution_path
+from nose.tools import *
+
+def setup():
+    # All of the paths used are relative, if we run the tests
+    # from another directory we need to chdir()
+    os.chdir(execution_path('.'))
+
+def test_vrt_rendering():
+    m = mapnik.Map(512,512)
+    mapnik.load_map(m,'../data/good_maps/vrt_colortable.xml')
+    m.zoom_all()
+    im = mapnik.Image(512,512)
+    mapnik.render(m,im)
+    actual = '/tmp/vrt_colortable.png'
+    expected = 'images/support/vrt_colortable.png'
+    im.save(actual)
+    expected_im = mapnik.Image.open(expected)
+    eq_(im.tostring(),expected_im.tostring(), 'failed comparing actual (%s) and expected(%s)' % (actual,'tests/python_tests/'+ expected))
+
+def test_tif_rendering_nodata():
+    m = mapnik.Map(512,512)
+    mapnik.load_map(m,'../data/good_maps/tiff_colortable.xml')
+    m.zoom_all()
+    im = mapnik.Image(512,512)
+    mapnik.render(m,im)
+    actual = '/tmp/tif_colortable.png'
+    expected = 'images/support/tif_colortable.png'
+    im.save(actual)
+    expected_im = mapnik.Image.open(expected)
+    eq_(im.tostring(),expected_im.tostring(), 'failed comparing actual (%s) and expected(%s)' % (actual,'tests/python_tests/'+ expected))
+
+if __name__ == "__main__":
+    setup()
+    [eval(run)() for run in dir() if 'test_' in run]