diff --git a/modules/library/referencing/pom.xml b/modules/library/referencing/pom.xml index 80e1e60..2ef6308 100644 --- a/modules/library/referencing/pom.xml +++ b/modules/library/referencing/pom.xml @@ -189,6 +189,11 @@ gt-metadata ${project.version} + + jgridshift + jgridshift + 1.0 + diff --git a/modules/library/referencing/src/main/java/org/geotools/referencing/ReferencingFactoryFinder.java b/modules/library/referencing/src/main/java/org/geotools/referencing/ReferencingFactoryFinder.java index 9240d28..40ed772 100644 --- a/modules/library/referencing/src/main/java/org/geotools/referencing/ReferencingFactoryFinder.java +++ b/modules/library/referencing/src/main/java/org/geotools/referencing/ReferencingFactoryFinder.java @@ -18,36 +18,37 @@ package org.geotools.referencing; import java.io.IOException; import java.io.Writer; -import java.util.Set; -import java.util.Locale; import java.util.Collections; import java.util.LinkedHashSet; +import java.util.Locale; +import java.util.Set; + import javax.imageio.spi.ServiceRegistry; +import org.geotools.factory.FactoryCreator; +import org.geotools.factory.FactoryFinder; +import org.geotools.factory.FactoryRegistry; +import org.geotools.factory.FactoryRegistryException; +import org.geotools.factory.GeoTools; +import org.geotools.factory.Hints; +import org.geotools.metadata.iso.citation.Citations; +import org.geotools.referencing.factory.gridshift.GridShiftLocator; +import org.geotools.resources.Arguments; +import org.geotools.resources.LazySet; import org.opengis.metadata.Identifier; import org.opengis.metadata.citation.Citation; -import org.opengis.referencing.Factory; import org.opengis.referencing.AuthorityFactory; -import org.opengis.referencing.crs.CRSFactory; +import org.opengis.referencing.Factory; import org.opengis.referencing.crs.CRSAuthorityFactory; -import org.opengis.referencing.cs.CSFactory; +import org.opengis.referencing.crs.CRSFactory; import org.opengis.referencing.cs.CSAuthorityFactory; -import org.opengis.referencing.datum.DatumFactory; +import org.opengis.referencing.cs.CSFactory; import org.opengis.referencing.datum.DatumAuthorityFactory; -import org.opengis.referencing.operation.CoordinateOperationFactory; +import org.opengis.referencing.datum.DatumFactory; import org.opengis.referencing.operation.CoordinateOperationAuthorityFactory; +import org.opengis.referencing.operation.CoordinateOperationFactory; import org.opengis.referencing.operation.MathTransformFactory; -import org.geotools.factory.Hints; -import org.geotools.factory.GeoTools; -import org.geotools.factory.FactoryFinder; -import org.geotools.factory.FactoryCreator; -import org.geotools.factory.FactoryRegistry; -import org.geotools.factory.FactoryRegistryException; -import org.geotools.metadata.iso.citation.Citations; -import org.geotools.resources.Arguments; -import org.geotools.resources.LazySet; - /** * Defines static methods used to access the application's default {@linkplain Factory @@ -120,7 +121,8 @@ public final class ReferencingFactoryFinder extends FactoryFinder { CRSAuthorityFactory.class, MathTransformFactory.class, CoordinateOperationFactory.class, - CoordinateOperationAuthorityFactory.class}); + CoordinateOperationAuthorityFactory.class, + GridShiftLocator.class}); } return registry; } @@ -486,7 +488,20 @@ loop: for (int i=0; ; i++) { { return getFactories(CoordinateOperationAuthorityFactory.class, hints); } - + + /** + * Returns a set of all available implementations for the + * {@link GridShiftLocator} interface. + * + * @param hints An optional map of hints, or {@code null} if none. + * @return Set of available grid shift factory implementations. + */ + public static Set getGridShiftLocators( + final Hints hints) + { + return getFactories(GridShiftLocator.class, hints); + } + /** * Returns the first implementation of {@link MathTransformFactory} matching the specified * hints. If no implementation matches, a new one is created if possible or an exception is diff --git a/modules/library/referencing/src/main/java/org/geotools/referencing/factory/gridshift/ClasspathGridShiftLocator.java b/modules/library/referencing/src/main/java/org/geotools/referencing/factory/gridshift/ClasspathGridShiftLocator.java new file mode 100644 index 0000000..6b42908 --- /dev/null +++ b/modules/library/referencing/src/main/java/org/geotools/referencing/factory/gridshift/ClasspathGridShiftLocator.java @@ -0,0 +1,42 @@ +/* + * GeoTools - The Open Source Java GIS Toolkit + * http://geotools.org + * + * (C) 2002-2012, Open Source Geospatial Foundation (OSGeo) + * + * This library is free software; you can redistribute it and/or + * modify it under the terms of the GNU Lesser General Public + * License as published by the Free Software Foundation; + * version 2.1 of the License. + * + * This library is distributed in the hope that it will be useful, + * but WITHOUT ANY WARRANTY; without even the implied warranty of + * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU + * Lesser General Public License for more details. + */ +package org.geotools.referencing.factory.gridshift; + +import java.net.URL; + +import org.geotools.metadata.iso.citation.Citations; +import org.opengis.metadata.citation.Citation; + +/** + * Default grid shift file locator, looks up grids in the classpath + * + * @author Andrea Aime - GeoSolutions + * + */ +public class ClasspathGridShiftLocator implements GridShiftLocator { + + @Override + public Citation getVendor() { + return Citations.GEOTOOLS; + } + + @Override + public URL locateGrid(String grid) { + return getClass().getResource(grid); + } + +} diff --git a/modules/library/referencing/src/main/java/org/geotools/referencing/factory/gridshift/GridShiftLocator.java b/modules/library/referencing/src/main/java/org/geotools/referencing/factory/gridshift/GridShiftLocator.java new file mode 100644 index 0000000..282308a --- /dev/null +++ b/modules/library/referencing/src/main/java/org/geotools/referencing/factory/gridshift/GridShiftLocator.java @@ -0,0 +1,38 @@ +/* + * GeoTools - The Open Source Java GIS Toolkit + * http://geotools.org + * + * (C) 2002-2012, Open Source Geospatial Foundation (OSGeo) + * + * This library is free software; you can redistribute it and/or + * modify it under the terms of the GNU Lesser General Public + * License as published by the Free Software Foundation; + * version 2.1 of the License. + * + * This library is distributed in the hope that it will be useful, + * but WITHOUT ANY WARRANTY; without even the implied warranty of + * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU + * Lesser General Public License for more details. + */ +package org.geotools.referencing.factory.gridshift; + +import java.net.URL; + +import org.opengis.referencing.Factory; + + +/** + * Provides a hook to locate grid shift files, such as NTv1, NTv2 and NADCON ones + * + * Andrea Aime - Geosolutions + */ +public interface GridShiftLocator extends Factory { + + /** + * Locate the specified resource. + * + * @param grid the grid name/location + * @return the fully resolved URL of the grid or null, if the resource cannot be located. + */ + URL locateGrid(String grid); +} diff --git a/modules/library/referencing/src/main/java/org/geotools/referencing/factory/gridshift/NADCONGridShiftFactory.java b/modules/library/referencing/src/main/java/org/geotools/referencing/factory/gridshift/NADCONGridShiftFactory.java new file mode 100644 index 0000000..fafda54 --- /dev/null +++ b/modules/library/referencing/src/main/java/org/geotools/referencing/factory/gridshift/NADCONGridShiftFactory.java @@ -0,0 +1,440 @@ +/* + * GeoTools - The Open Source Java GIS Toolkit + * http://geotools.org + * + * (C) 2002-2012, Open Source Geospatial Foundation (OSGeo) + * + * This library is free software; you can redistribute it and/or + * modify it under the terms of the GNU Lesser General Public + * License as published by the Free Software Foundation; + * version 2.1 of the License. + * + * This library is distributed in the hope that it will be useful, + * but WITHOUT ANY WARRANTY; without even the implied warranty of + * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU + * Lesser General Public License for more details. + */ +package org.geotools.referencing.factory.gridshift; + +import java.io.BufferedReader; +import java.io.EOFException; +import java.io.File; +import java.io.FileInputStream; +import java.io.IOException; +import java.io.InputStream; +import java.io.InputStreamReader; +import java.net.URL; +import java.nio.ByteBuffer; +import java.nio.ByteOrder; +import java.nio.channels.Channels; +import java.nio.channels.ReadableByteChannel; +import java.util.StringTokenizer; +import java.util.logging.Logger; + +import org.geotools.factory.BufferedFactory; +import org.geotools.referencing.factory.ReferencingFactory; +import org.geotools.resources.i18n.ErrorKeys; +import org.geotools.resources.i18n.Errors; +import org.geotools.util.SoftValueHashMap; +import org.geotools.util.logging.Logging; +import org.opengis.referencing.FactoryException; + +/** + * Loads and caches NADCON grid shifts + * + * @author Andrea Aime - GeoSolutions + * + */ +public class NADCONGridShiftFactory extends ReferencingFactory implements BufferedFactory { + + static final class NADCONKey { + String latFile; + + String longFile; + + public NADCONKey(String latFile, String longFile) { + super(); + this.latFile = latFile; + this.longFile = longFile; + } + + @Override + public int hashCode() { + final int prime = 31; + int result = 1; + result = prime * result + ((latFile == null) ? 0 : latFile.hashCode()); + result = prime * result + ((longFile == null) ? 0 : longFile.hashCode()); + return result; + } + + @Override + public boolean equals(Object obj) { + if (this == obj) + return true; + if (obj == null) + return false; + if (getClass() != obj.getClass()) + return false; + NADCONKey other = (NADCONKey) obj; + if (latFile == null) { + if (other.latFile != null) + return false; + } else if (!latFile.equals(other.latFile)) + return false; + if (longFile == null) { + if (other.longFile != null) + return false; + } else if (!longFile.equals(other.longFile)) + return false; + return true; + } + + } + + /** + * The number of hard references to hold internally. + */ + private static final int GRID_CACHE_HARD_REFERENCES = 10; + + /** + * Logger. + */ + protected static final Logger LOGGER = Logging.getLogger("org.geotools.referencing"); + + /** + * The soft cache that holds loaded grids. + */ + private SoftValueHashMap gridCache; + + /** + * Constructs a factory with the default priority. + */ + public NADCONGridShiftFactory() { + gridCache = new SoftValueHashMap(GRID_CACHE_HARD_REFERENCES); + } + + public NADConGridShift loadGridShift(URL latGridURL, URL longGridURL) throws FactoryException { + NADCONKey key = new NADCONKey(latGridURL.toExternalForm(), longGridURL.toExternalForm()); + synchronized (gridCache) { // Prevent simultaneous threads trying to load same grid + NADConGridShift grid = gridCache.get(key); + if (grid != null) { // Cached: + return grid; // - Return + } else { // Not cached: + grid = loadGridShiftInternal(latGridURL, longGridURL); // - Load + if (grid != null) { + gridCache.put(key, grid); // - Cache + return grid; // - Return + } + } + throw new FactoryException("NTv2 Grid " + latGridURL + ", " + longGridURL + + " could not be created."); + } + } + + private NADConGridShift loadGridShiftInternal(URL latGridURL, URL longGridURL) + throws FactoryException { + // decide if text or binary grid will be used + String latGridName = latGridURL.getFile(); + String longGridName = longGridURL.getFile(); + try { + if ((latGridName.endsWith(".las") && longGridName.endsWith(".los")) + || (latGridName.endsWith(".LAS") && longGridName.endsWith(".LOS"))) { + return loadBinaryGrid(latGridURL, longGridURL); + } else if ((latGridName.endsWith(".laa") && longGridName.endsWith(".loa")) + || (latGridName.endsWith(".LAA") && longGridName.endsWith(".LOA"))) { + return loadTextGrid(latGridURL, longGridURL); + } else { + throw new FactoryException(Errors.format(ErrorKeys.UNSUPPORTED_FILE_TYPE_$2, + latGridName.substring(latGridName.lastIndexOf('.') + 1), + longGridName.substring(longGridName.lastIndexOf('.') + 1))); + // Note: the +1 above hide the dot, but also make sure that the code is + // valid even if the path do not contains '.' at all (-1 + 1 == 0). + } + } catch (IOException exception) { + final Throwable cause = exception.getCause(); + if (cause instanceof FactoryException) { + throw (FactoryException) cause; + } + throw new FactoryException(exception.getLocalizedMessage(), exception); + } + } + + /** + * Reads latitude and longitude binary grid shift file data into {@link grid}. The file is + * organized into records, with the first record containing the header information, followed by + * the shift data. The header values are: text describing grid (64 bytes), num. columns (int), + * num. rows (int), num. z (int), min x (float), delta x (float), min y (float), delta y (float) + * and angle (float). Each record is num. columns 4 bytes + 4 byte separator long and the file + * contains num. rows + 1 (for the header) records. The data records (with the grid shift + * values) are all floats and have a 4 byte separator (0's) before the data. Row records are + * organized from low y (latitude) to high and columns are orderd from low longitude to high. + * Everything is written in low byte order. + * + * @param latGridUrl URL to the binary latitude shift file (.las extention). + * @param longGridUrl URL to the binary longitude shift file (.los extention). + * @throws IOException if the data files cannot be read. + * @throws FactoryException if there is an inconsistency in the data + */ + private NADConGridShift loadBinaryGrid(final URL latGridUrl, final URL longGridUrl) + throws IOException, FactoryException { + final int HEADER_BYTES = 96; + final int SEPARATOR_BYTES = 4; + final int DESCRIPTION_LENGTH = 64; + ReadableByteChannel latChannel; + ReadableByteChannel longChannel; + ByteBuffer latBuffer; + ByteBuffer longBuffer; + + // ////////////////////// + // setup + // ////////////////////// + latChannel = getReadChannel(latGridUrl); + latBuffer = fillBuffer(latChannel, HEADER_BYTES); + + longChannel = getReadChannel(longGridUrl); + longBuffer = fillBuffer(longChannel, HEADER_BYTES); + + // ////////////////////// + // read header info + // ////////////////////// + // skip the header description + latBuffer.position(latBuffer.position() + DESCRIPTION_LENGTH); + + int nc = latBuffer.getInt(); + int nr = latBuffer.getInt(); + int nz = latBuffer.getInt(); + + float xmin = latBuffer.getFloat(); + float dx = latBuffer.getFloat(); + float ymin = latBuffer.getFloat(); + float dy = latBuffer.getFloat(); + + float angle = latBuffer.getFloat(); + float xmax = xmin + ((nc - 1) * dx); + float ymax = ymin + ((nr - 1) * dy); + + // skip the longitude header description + longBuffer.position(longBuffer.position() + DESCRIPTION_LENGTH); + + // check that latitude grid header is the same as for latitude grid + if ((nc != longBuffer.getInt()) || (nr != longBuffer.getInt()) + || (nz != longBuffer.getInt()) || (xmin != longBuffer.getFloat()) + || (dx != longBuffer.getFloat()) || (ymin != longBuffer.getFloat()) + || (dy != longBuffer.getFloat()) || (angle != longBuffer.getFloat())) { + throw new FactoryException(Errors.format(ErrorKeys.GRID_LOCATIONS_UNEQUAL)); + } + + // ////////////////////// + // read grid shift data into LocalizationGrid + // ////////////////////// + final int RECORD_LENGTH = (nc * 4) + SEPARATOR_BYTES; + final int NUM_BYTES_LEFT = ((nr + 1) * RECORD_LENGTH) - HEADER_BYTES; + final int START_OF_DATA = RECORD_LENGTH - HEADER_BYTES; + + latBuffer = fillBuffer(latChannel, NUM_BYTES_LEFT); + latBuffer.position(START_OF_DATA); // start of second record (data) + + longBuffer = fillBuffer(longChannel, NUM_BYTES_LEFT); + longBuffer.position(START_OF_DATA); + + NADConGridShift gridShift = new NADConGridShift(xmin, ymin, xmax, ymax, dx, dy, nc, nr); + + int i = 0; + int j = 0; + for (i = 0; i < nr; i++) { + latBuffer.position(latBuffer.position() + SEPARATOR_BYTES); // skip record separator + longBuffer.position(longBuffer.position() + SEPARATOR_BYTES); + + for (j = 0; j < nc; j++) { + gridShift.setLocalizationPoint(j, i, longBuffer.getFloat(), latBuffer.getFloat()); + } + } + + assert i == nr : i; + assert j == nc : j; + + return gridShift; + } + + /** + * Returns a new bytebuffer, of numBytes length and little endian byte order, filled from the + * channel. + * + * @param channel the channel to fill the buffer from + * @param numBytes number of bytes to read + * @return a new bytebuffer filled from the channel + * @throws IOException if there is a problem reading the channel + * @throws EOFException if the end of the channel is reached + */ + private ByteBuffer fillBuffer(ReadableByteChannel channel, int numBytes) throws IOException { + ByteBuffer buf = ByteBuffer.allocate(numBytes); + + if (fill(buf, channel) == -1) { + throw new EOFException(Errors.format(ErrorKeys.END_OF_DATA_FILE)); + } + + buf.flip(); + buf.order(ByteOrder.LITTLE_ENDIAN); + + return buf; + } + + /** + * Fills the bytebuffer from the channel. Code was lifted from ShapefileDataStore. + * + * @param buffer bytebuffer to fill from the channel + * @param channel channel to fill the buffer from + * @return number of bytes read + * @throws IOException if there is a problem reading the channel + */ + private int fill(ByteBuffer buffer, ReadableByteChannel channel) throws IOException { + int r = buffer.remaining(); + + // channel reads return -1 when EOF or other error + // because they a non-blocking reads, 0 is a valid return value!! + while ((buffer.remaining() > 0) && (r != -1)) { + r = channel.read(buffer); + } + + if (r == -1) { + buffer.limit(buffer.position()); + } + + return r; + } + + /** + * Obtain a ReadableByteChannel from the given URL. If the url protocol is file, a FileChannel + * will be returned. Otherwise a generic channel will be obtained from the urls input stream. + * Code swiped from ShapefileDataStore. + * + * @param url URL to create the channel from + * @return a new PeadableByteChannel from the input url + * @throws IOException if there is a problem creating the channel + */ + private ReadableByteChannel getReadChannel(URL url) throws IOException { + ReadableByteChannel channel = null; + + if (url.getProtocol().equals("file")) { + File file = new File(url.getFile()); + + if (!file.exists() || !file.canRead()) { + throw new IOException(Errors.format(ErrorKeys.FILE_DOES_NOT_EXIST_$1, file)); + } + + FileInputStream in = new FileInputStream(file); + channel = in.getChannel(); + } else { + InputStream in = url.openConnection().getInputStream(); + channel = Channels.newChannel(in); + } + + return channel; + } + + /** + * Reads latitude and longitude text grid shift file data into {@link grid}. The first two lines + * of the shift data file contain the header, with the first being a description of the grid. + * The second line contains 8 values separated by spaces: num. columns, num. rows, num. z, min + * x, delta x, min y, delta y and angle. Shift data values follow this and are also separated by + * spaces. Row records are organized from low y (latitude) to high and columns are orderd from + * low longitude to high. + * + * @param latGridUrl URL to the text latitude shift file (.laa extention). + * @param longGridUrl URL to the text longitude shift file (.loa extention). + * @throws IOException if the data files cannot be read. + * @throws FactoryException if there is an inconsistency in the data + */ + private NADConGridShift loadTextGrid(URL latGridUrl, URL longGridUrl) throws IOException, + FactoryException { + String latLine; + String longLine; + StringTokenizer latSt; + StringTokenizer longSt; + + // ////////////////////// + // setup + // ////////////////////// + InputStreamReader latIsr = new InputStreamReader(latGridUrl.openStream()); + BufferedReader latBr = new BufferedReader(latIsr); + + InputStreamReader longIsr = new InputStreamReader(longGridUrl.openStream()); + BufferedReader longBr = new BufferedReader(longIsr); + + // ////////////////////// + // read header info + // ////////////////////// + latLine = latBr.readLine(); // skip header description + latLine = latBr.readLine(); + latSt = new StringTokenizer(latLine, " "); + + if (latSt.countTokens() != 8) { + throw new FactoryException(Errors.format(ErrorKeys.HEADER_UNEXPECTED_LENGTH_$1, + String.valueOf(latSt.countTokens()))); + } + + int nc = Integer.parseInt(latSt.nextToken()); + int nr = Integer.parseInt(latSt.nextToken()); + int nz = Integer.parseInt(latSt.nextToken()); + + float xmin = Float.parseFloat(latSt.nextToken()); + float dx = Float.parseFloat(latSt.nextToken()); + float ymin = Float.parseFloat(latSt.nextToken()); + float dy = Float.parseFloat(latSt.nextToken()); + + float angle = Float.parseFloat(latSt.nextToken()); + float xmax = xmin + ((nc - 1) * dx); + float ymax = ymin + ((nr - 1) * dy); + + // now read long shift grid + longLine = longBr.readLine(); // skip header description + longLine = longBr.readLine(); + longSt = new StringTokenizer(longLine, " "); + + if (longSt.countTokens() != 8) { + throw new FactoryException(Errors.format(ErrorKeys.HEADER_UNEXPECTED_LENGTH_$1, + String.valueOf(longSt.countTokens()))); + } + + // check that latitude grid header is the same as for latitude grid + if ((nc != Integer.parseInt(longSt.nextToken())) + || (nr != Integer.parseInt(longSt.nextToken())) + || (nz != Integer.parseInt(longSt.nextToken())) + || (xmin != Float.parseFloat(longSt.nextToken())) + || (dx != Float.parseFloat(longSt.nextToken())) + || (ymin != Float.parseFloat(longSt.nextToken())) + || (dy != Float.parseFloat(longSt.nextToken())) + || (angle != Float.parseFloat(longSt.nextToken()))) { + throw new FactoryException(Errors.format(ErrorKeys.GRID_LOCATIONS_UNEQUAL)); + } + + // ////////////////////// + // read grid shift data into LocalizationGrid + // ////////////////////// + NADConGridShift gridShift = new NADConGridShift(xmin, ymin, xmax, ymax, dx, dy, nc, nr); + + int i = 0; + int j = 0; + for (i = 0; i < nr; i++) { + for (j = 0; j < nc;) { + latLine = latBr.readLine(); + latSt = new StringTokenizer(latLine, " "); + longLine = longBr.readLine(); + longSt = new StringTokenizer(longLine, " "); + + while (latSt.hasMoreTokens() && longSt.hasMoreTokens()) { + gridShift.setLocalizationPoint(j, i, + (double) Float.parseFloat(longSt.nextToken()), + (double) Float.parseFloat(latSt.nextToken())); + ++j; + } + } + } + + assert i == nr : i; + assert j == nc : j; + + return gridShift; + } + +} diff --git a/modules/library/referencing/src/main/java/org/geotools/referencing/factory/gridshift/NADConGridShift.java b/modules/library/referencing/src/main/java/org/geotools/referencing/factory/gridshift/NADConGridShift.java new file mode 100644 index 0000000..6520f9d --- /dev/null +++ b/modules/library/referencing/src/main/java/org/geotools/referencing/factory/gridshift/NADConGridShift.java @@ -0,0 +1,157 @@ +/* + * GeoTools - The Open Source Java GIS Toolkit + * http://geotools.org + * + * (C) 2002-2012, Open Source Geospatial Foundation (OSGeo) + * + * This library is free software; you can redistribute it and/or + * modify it under the terms of the GNU Lesser General Public + * License as published by the Free Software Foundation; + * version 2.1 of the License. + * + * This library is distributed in the hope that it will be useful, + * but WITHOUT ANY WARRANTY; without even the implied warranty of + * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU + * Lesser General Public License for more details. + */ +package org.geotools.referencing.factory.gridshift; + +import org.geotools.referencing.operation.builder.LocalizationGrid; + +/** + * A NADCON localization grid + * + * @author Andrea Aime - GeoSolutions + * @author Rueben Schulz + */ +public class NADConGridShift extends LocalizationGrid { + + /** + * The minimum longitude value covered by this grid (decimal degrees) + */ + private double minX; + + /** + * The minimum latitude value covered by this grid (decimal degrees) + */ + private double minY; + + /** + * The maximum longitude value covered by this grid (decimal degrees) + */ + private double maxX; + + /** + * The maximum latitude value covered by this grid (decimal degrees) + */ + private double maxY; + + /** + * The difference between longitude grid points (decimal degrees) + */ + private double dx; + + /** + * The difference between latitude grid points (decimal degrees) + */ + private double dy; + + public NADConGridShift(double xmin, double ymin, double xmax, double ymax, double dx, double dy, int width, int height) { + super(width, height); + this.minX = xmin; + this.maxX = xmax; + this.minY = ymin; + this.maxY = ymax; + this.dx = dx; + this.dy = dy; + } + + /** + * Returns a hash value for this transform. To make this faster it does not + * check the grid values. + * + * @return a hash value for this transform. + */ + @Override + public final int hashCode() { + final long code = Double.doubleToLongBits(minX) + + (37 * (Double.doubleToLongBits(minY) + + (37 * (Double.doubleToLongBits(maxX) + + (37 * (Double.doubleToLongBits(maxY) + + (37 * (Double.doubleToLongBits(dx) + + (37 * (Double.doubleToLongBits(dy))))))))))); + + return (int) code ^ (int) (code >>> 32); + } + + /** + * Compares the specified object with this math transform for equality. + * + * @param object the object to compare to + * @return {@code true} if the objects are equal. + */ + @Override + public final boolean equals(final Object object) { + if (object == this) { + // Slight optimization + return true; + } + + if (super.equals(object)) { + final NADConGridShift that = (NADConGridShift) object; + + return (Double.doubleToLongBits(this.minX) == Double.doubleToLongBits(that.minX)) + && (Double.doubleToLongBits(this.minY) == Double.doubleToLongBits(that.minY)) + && (Double.doubleToLongBits(this.maxX) == Double.doubleToLongBits(that.maxX)) + && (Double.doubleToLongBits(this.maxY) == Double.doubleToLongBits(that.maxY)) + && (Double.doubleToLongBits(this.dx) == Double.doubleToLongBits(that.dx)) + && (Double.doubleToLongBits(this.dy) == Double.doubleToLongBits(that.dy)); + } + + return false; + } + + /** + * The minimum longitude value covered by this grid (decimal degrees) + */ + public double getMinX() { + return minX; + } + + /** + * The minimum latitude value covered by this grid (decimal degrees) + */ + public double getMinY() { + return minY; + } + + /** + * The maximum longitude value covered by this grid (decimal degrees) + */ + public double getMaxX() { + return maxX; + } + + /** + * The maximum latitude value covered by this grid (decimal degrees) + */ + public double getMaxY() { + return maxY; + } + + /** + * The difference between longitude grid points (decimal degrees) + */ + public double getDx() { + return dx; + } + + /** + * The difference between latitude grid points (decimal degrees) + */ + public double getDy() { + return dy; + } + + +} diff --git a/modules/library/referencing/src/main/java/org/geotools/referencing/factory/gridshift/NTv2GridShiftFactory.java b/modules/library/referencing/src/main/java/org/geotools/referencing/factory/gridshift/NTv2GridShiftFactory.java new file mode 100644 index 0000000..d9bcedd --- /dev/null +++ b/modules/library/referencing/src/main/java/org/geotools/referencing/factory/gridshift/NTv2GridShiftFactory.java @@ -0,0 +1,223 @@ +/* + * GeoTools - The Open Source Java GIS Toolkit + * http://geotools.org + * + * (C) 2002-2012, Open Source Geospatial Foundation (OSGeo) + * + * This library is free software; you can redistribute it and/or + * modify it under the terms of the GNU Lesser General Public + * License as published by the Free Software Foundation; + * version 2.1 of the License. + * + * This library is distributed in the hope that it will be useful, + * but WITHOUT ANY WARRANTY; without even the implied warranty of + * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU + * Lesser General Public License for more details. + */ +package org.geotools.referencing.factory.gridshift; + +import java.io.File; +import java.io.FileNotFoundException; +import java.io.IOException; +import java.io.InputStream; +import java.io.RandomAccessFile; +import java.net.URL; +import java.util.logging.Level; +import java.util.logging.Logger; + +import org.geotools.factory.AbstractFactory; +import org.geotools.factory.BufferedFactory; +import org.geotools.referencing.factory.ReferencingFactory; +import org.geotools.resources.i18n.ErrorKeys; +import org.geotools.resources.i18n.Errors; +import org.geotools.util.SoftValueHashMap; +import org.geotools.util.logging.Logging; +import org.opengis.referencing.FactoryException; + +import au.com.objectix.jgridshift.GridShiftFile; + +/** + * Loads and caches NTv2 grid files. Thisthat incorporates a soft cache mechanism to keep grids in + * memory when first loaded. It also checks NTv2 grid file format in {@link #isNTv2Grid(String)} + * method. + * + * @author Oscar Fonts + */ +public class NTv2GridShiftFactory extends ReferencingFactory implements BufferedFactory { + + /** + * The number of hard references to hold internally. + */ + private static final int GRID_CACHE_HARD_REFERENCES = 10; + + /** + * Logger. + */ + protected static final Logger LOGGER = Logging.getLogger("org.geotools.referencing"); + + /** + * The soft cache that holds loaded grids. + */ + private SoftValueHashMap ntv2GridCache; + + /** + * Constructs a factory with the default priority. + */ + public NTv2GridShiftFactory() { + super(); + ntv2GridCache = new SoftValueHashMap(GRID_CACHE_HARD_REFERENCES); + } + + /** + * Constructs an instance using the specified priority level. + * + * @param priority The priority for this factory, as a number between + * {@link AbstractFactory#MINIMUM_PRIORITY MINIMUM_PRIORITY} and + * {@link AbstractFactory#MAXIMUM_PRIORITY MAXIMUM_PRIORITY} inclusive. + */ + public NTv2GridShiftFactory(final int priority) { + super(priority); + ntv2GridCache = new SoftValueHashMap(GRID_CACHE_HARD_REFERENCES); + } + + /** + * Performs a NTv2 grid file lookup given its name, and checks for file format correctness. + * + * @param name The NTv2 grid file name + * @return {@code true} if file exists and is valid, {@code false} otherwise + */ + public boolean isNTv2Grid(URL location) { + if (location != null) { + return isNTv2GridFileValid(location); // Check + } else { + return false; + } + } + + /** + * Creates a NTv2 Grid. + * + * @param name The NTv2 grid name + * @return the grid + * @throws FactoryException if grid cannot be created + */ + public GridShiftFile createNTv2Grid(URL gridLocation) throws FactoryException { + if(gridLocation == null) { + throw new FactoryException("The grid location must be not null"); + } + + synchronized (ntv2GridCache) { // Prevent simultaneous threads trying to load same grid + GridShiftFile grid = ntv2GridCache.get(gridLocation.toExternalForm()); + if (grid != null) { // Cached: + return grid; // - Return + } else { // Not cached: + if (gridLocation != null) { + grid = loadNTv2Grid(gridLocation); // - Load + if (grid != null) { + ntv2GridCache.put(gridLocation.toExternalForm(), grid); // - Cache + return grid; // - Return + } + } + throw new FactoryException("NTv2 Grid " + gridLocation + " could not be created."); + } + } + } + + + + /** + * Checks if a given resource is a valid NTv2 file without fully loading it. + * + * If file is not valid, the cause is logged at {@link Level#WARNING warning level}. + * + * @param location the NTv2 file absolute path + * @return true if file has NTv2 format, false otherwise + */ + protected boolean isNTv2GridFileValid(URL url) { + RandomAccessFile raf = null; + InputStream is = null; + try { + + // Loading as RandomAccessFile doesn't load the full grid + // in memory, but is a quick method to see if file format + // is NTv2. + if (url.getProtocol().equals("file")) { + File file = new File(url.getFile()); + + if (!file.exists() || !file.canRead()) { + throw new IOException(Errors.format(ErrorKeys.FILE_DOES_NOT_EXIST_$1, file)); + } + + raf = new RandomAccessFile(file, "r"); + + // will throw an exception if not a valid file + new GridShiftFile().loadGridShiftFile(raf); + } else { + InputStream in = url.openConnection().getInputStream(); + + // will throw an exception if not a valid file + new GridShiftFile().loadGridShiftFile(in, false); + } + + return true; // No exception thrown => valid file. + } catch (IllegalArgumentException e) { + // This usually means resource is not a valid NTv2 file. + // Let exception message describe the cause. + LOGGER.log(Level.WARNING, e.getLocalizedMessage(), e); + return false; + } catch (IOException e) { + LOGGER.log(Level.WARNING, e.getLocalizedMessage(), e); + return false; + } finally { + try { + if (raf != null) + raf.close(); + } catch (IOException e) { + } + + try { + if (is != null) + is.close(); + } catch (IOException e) { + } + } + } + + /** + * Loads the grid in memory. + * + * If file cannot be loaded, the cause is logged at {@link Level#SEVERE severe level}. + * + * @param location the NTv2 file absolute path + * @return the grid, or {@code null} on error + * @throws FactoryException + */ + private GridShiftFile loadNTv2Grid(URL location) throws FactoryException { + InputStream in = null; + try { + GridShiftFile grid = new GridShiftFile(); + in = location.openStream(); + grid.loadGridShiftFile(in, false); // Load full grid in memory + in.close(); + return grid; + } catch (FileNotFoundException e) { + LOGGER.log(Level.SEVERE, e.getLocalizedMessage(), e); + return null; + } catch (IOException e) { + LOGGER.log(Level.SEVERE, e.getLocalizedMessage(), e); + return null; + } catch (IllegalArgumentException e) { + LOGGER.log(Level.SEVERE, e.getLocalizedMessage(), e); + throw new FactoryException(e.getLocalizedMessage(), e); + } finally { + try { + if (in != null) { + in.close(); + } + } catch (IOException e) { + // never mind + } + } + } + +} diff --git a/modules/library/referencing/src/main/java/org/geotools/referencing/operation/MathTransformProvider.java b/modules/library/referencing/src/main/java/org/geotools/referencing/operation/MathTransformProvider.java index dc01761..eb01af2 100644 --- a/modules/library/referencing/src/main/java/org/geotools/referencing/operation/MathTransformProvider.java +++ b/modules/library/referencing/src/main/java/org/geotools/referencing/operation/MathTransformProvider.java @@ -238,7 +238,7 @@ public abstract class MathTransformProvider extends DefaultOperationMethod { * Put the identifiers into a properties map suitable for {@link IdentifiedObject} * constructor. */ - private static Map toMap(final ReferenceIdentifier[] identifiers) { + protected static Map toMap(final ReferenceIdentifier[] identifiers) { ensureNonNull("identifiers", identifiers); if (identifiers.length == 0) { throw new IllegalArgumentException(Errors.format(ErrorKeys.EMPTY_ARRAY)); diff --git a/modules/library/referencing/src/main/java/org/geotools/referencing/operation/transform/NADCONTransform.java b/modules/library/referencing/src/main/java/org/geotools/referencing/operation/transform/NADCONTransform.java index 84b31ad..dbee749 100644 --- a/modules/library/referencing/src/main/java/org/geotools/referencing/operation/transform/NADCONTransform.java +++ b/modules/library/referencing/src/main/java/org/geotools/referencing/operation/transform/NADCONTransform.java @@ -16,49 +16,42 @@ */ package org.geotools.referencing.operation.transform; -import java.io.BufferedReader; -import java.io.EOFException; -import java.io.File; -import java.io.FileInputStream; import java.io.IOException; -import java.io.InputStream; -import java.io.InputStreamReader; import java.io.ObjectInputStream; import java.io.PrintWriter; import java.io.Serializable; +import java.net.URI; import java.net.URL; -import java.net.MalformedURLException; -import java.nio.ByteBuffer; -import java.nio.ByteOrder; -import java.nio.channels.Channels; -import java.nio.channels.ReadableByteChannel; -import java.util.StringTokenizer; import java.util.prefs.Preferences; -import org.opengis.parameter.GeneralParameterValue; -import org.opengis.parameter.ParameterDescriptor; -import org.opengis.parameter.ParameterDescriptorGroup; -import org.opengis.parameter.ParameterNotFoundException; -import org.opengis.parameter.ParameterValue; -import org.opengis.parameter.ParameterValueGroup; -import org.opengis.referencing.FactoryException; -import org.opengis.referencing.operation.MathTransform; -import org.opengis.referencing.operation.MathTransform2D; -import org.opengis.referencing.operation.Transformation; -import org.opengis.referencing.operation.TransformException; - import org.geotools.metadata.iso.citation.Citations; import org.geotools.parameter.DefaultParameterDescriptor; import org.geotools.parameter.Parameter; import org.geotools.parameter.ParameterGroup; import org.geotools.referencing.NamedIdentifier; +import org.geotools.referencing.ReferencingFactoryFinder; +import org.geotools.referencing.factory.gridshift.GridShiftLocator; +import org.geotools.referencing.factory.gridshift.NADCONGridShiftFactory; +import org.geotools.referencing.factory.gridshift.NADConGridShift; import org.geotools.referencing.operation.MathTransformProvider; import org.geotools.referencing.operation.builder.LocalizationGrid; import org.geotools.resources.Arguments; -import org.geotools.resources.i18n.Errors; import org.geotools.resources.i18n.ErrorKeys; +import org.geotools.resources.i18n.Errors; import org.geotools.resources.i18n.Vocabulary; import org.geotools.resources.i18n.VocabularyKeys; +import org.opengis.parameter.GeneralParameterValue; +import org.opengis.parameter.ParameterDescriptor; +import org.opengis.parameter.ParameterDescriptorGroup; +import org.opengis.parameter.ParameterNotFoundException; +import org.opengis.parameter.ParameterValue; +import org.opengis.parameter.ParameterValueGroup; +import org.opengis.referencing.FactoryException; +import org.opengis.referencing.NoSuchIdentifierException; +import org.opengis.referencing.operation.MathTransform; +import org.opengis.referencing.operation.MathTransform2D; +import org.opengis.referencing.operation.TransformException; +import org.opengis.referencing.operation.Transformation; /** @@ -140,6 +133,11 @@ public class NADCONTransform extends AbstractMathTransform implements MathTransf * Serial number for interoperability with different versions. */ private static final long serialVersionUID = -4707304160205218546L; + + /** + * The factory that loads the NADCON grids + */ + private static NADCONGridShiftFactory FACTORY = new NADCONGridShiftFactory(); /** * Preference node for the grid shift file location. @@ -170,43 +168,14 @@ public class NADCONTransform extends AbstractMathTransform implements MathTransf /** * Latitude grid shift file names. Output in WKT. */ - private final String latGridName; + private final URI latGridName; /** * Longitude grid shift file names. Output in WKT. */ - private final String longGridName; - - /** - * The minimum longitude value covered by this grid (decimal degrees) - */ - private double xmin; - - /** - * The minimum latitude value covered by this grid (decimal degrees) - */ - private double ymin; - - /** - * The maximum longitude value covered by this grid (decimal degrees) - */ - private double xmax; - - /** - * The maximum latitude value covered by this grid (decimal degrees) - */ - private double ymax; - - /** - * The difference between longitude grid points (decimal degrees) - */ - private double dx; - - /** - * The difference between latitude grid points (decimal degrees) - */ - private double dy; + private final URI longGridName; + /** * Longitude and latitude grid shift values. Values are organized from low * to high longitude (low x index to high) and low to high latitude (low y @@ -226,6 +195,12 @@ public class NADCONTransform extends AbstractMathTransform implements MathTransf private transient MathTransform2D inverse; /** + * The grid driving this transform + */ + NADConGridShift grid; + + + /** * Constructs a {@code NADCONTransform} from the specified grid shift files. * * @param latGridName path and name (or just name if {@link #GRID_LOCATION} @@ -240,42 +215,40 @@ public class NADCONTransform extends AbstractMathTransform implements MathTransf * (ie file extentions are unknown or there is an error reading the * grid files) */ - public NADCONTransform(final String latGridName, final String longGridName) + public NADCONTransform(final URI latGridName, final URI longGridName) throws ParameterNotFoundException, FactoryException { + if(latGridName == null) { + throw new NoSuchIdentifierException("Latitud grid shift file name is null", null); + } + + if(longGridName == null) { + throw new NoSuchIdentifierException("Latitud grid shift file name is null", null); + } + this.latGridName = latGridName; this.longGridName = longGridName; - - //decide if text or binary grid will be used - try { - final URL latGridURL = makeURL(latGridName); - final URL longGridURL = makeURL(longGridName); - - if ((latGridName.endsWith(".las") && longGridName.endsWith(".los")) - || (latGridName.endsWith(".LAS") && longGridName.endsWith(".LOS"))) { - loadBinaryGrid(latGridURL, longGridURL); - } else if ((latGridName.endsWith(".laa") && longGridName.endsWith(".loa")) - || (latGridName.endsWith(".LAA") && longGridName.endsWith(".LOA"))) { - loadTextGrid(latGridURL, longGridURL); - } else { - throw new FactoryException(Errors.format(ErrorKeys.UNSUPPORTED_FILE_TYPE_$2, - latGridName.substring(latGridName.lastIndexOf('.') + 1), - longGridName.substring(longGridName.lastIndexOf('.') + 1))); - // Note: the +1 above hide the dot, but also make sure that the code is - // valid even if the path do not contains '.' at all (-1 + 1 == 0). - } - - gridShiftTransform = gridShift.getMathTransform(); - } catch (IOException exception) { - final Throwable cause = exception.getCause(); - if (cause instanceof FactoryException) { - throw (FactoryException) cause; + + URL latGridURL = locateGrid(latGridName); + URL longGridURL = locateGrid(longGridName); + + this.grid = FACTORY.loadGridShift(latGridURL, longGridURL); + this.gridShiftTransform = grid.getMathTransform(); + } + + URL locateGrid(URI uri ) throws FactoryException { + String grid = uri.toString(); + for (GridShiftLocator locator : ReferencingFactoryFinder.getGridShiftLocators(null)) { + URL result = locator.locateGrid(grid); + if(result != null) { + return result; } - throw new FactoryException(exception.getLocalizedMessage(), - exception); - } + }; + + throw new FactoryException("Could not locate grid file " + grid); } + /** * Returns the parameter descriptors for this math transform. */ @@ -320,332 +293,6 @@ public class NADCONTransform extends AbstractMathTransform implements MathTransf } /** - * Returns a URL from the string representation. If the string has no - * path, the default path preferece is added. - * - * @param str a string representation of a URL - * @return a URL created from the string representation - * @throws MalformedURLException if the URL cannot be created - */ - private URL makeURL(final String str) throws MalformedURLException { - //has '/' or '\' or ':', so probably full path to file - if ((str.indexOf('\\') >= 0) || (str.indexOf('/') >= 0) || (str.indexOf(':') >= 0)) { - return makeURLfromString(str); - } else { - // just a file name , prepend base location - final Preferences prefs = Preferences.userNodeForPackage(NADCONTransform.class); - final String baseLocation = prefs.get(GRID_LOCATION, DEFAULT_GRID_LOCATION); - return makeURLfromString(baseLocation + "/" + str); - } - } - - /** - * Returns a URL based on a string representation. If no protocol is given, - * it is assumed to be a local file. - * - * @param str a string representation of a URL - * @return a URL created from the string representation - * @throws MalformedURLException if the URL cannot be created - */ - private URL makeURLfromString(final String str) throws MalformedURLException { - try { - return new URL(str); - } catch (MalformedURLException e) { - //try making this with a file protocal - return new URL("file", "", str); - } - } - - /** - * Reads latitude and longitude binary grid shift file data into {@link - * grid}. The file is organized into records, with the first record - * containing the header information, followed by the shift data. The - * header values are: text describing grid (64 bytes), num. columns (int), - * num. rows (int), num. z (int), min x (float), delta x (float), min y - * (float), delta y (float) and angle (float). Each record is num. - * columns 4 bytes + 4 byte separator long and the file contains num. - * rows + 1 (for the header) records. The data records (with the grid - * shift values) are all floats and have a 4 byte separator (0's) before - * the data. Row records are organized from low y (latitude) to high and - * columns are orderd from low longitude to high. Everything is written - * in low byte order. - * - * @param latGridUrl URL to the binary latitude shift file (.las extention). - * @param longGridUrl URL to the binary longitude shift file (.los extention). - * @throws IOException if the data files cannot be read. - * @throws FactoryException if there is an inconsistency in the data - */ - private void loadBinaryGrid(final URL latGridUrl, final URL longGridUrl) - throws IOException, FactoryException - { - final int HEADER_BYTES = 96; - final int SEPARATOR_BYTES = 4; - final int DESCRIPTION_LENGTH = 64; - ReadableByteChannel latChannel; - ReadableByteChannel longChannel; - ByteBuffer latBuffer; - ByteBuffer longBuffer; - - //////////////////////// - //setup - //////////////////////// - latChannel = getReadChannel(latGridUrl); - latBuffer = fillBuffer(latChannel, HEADER_BYTES); - - longChannel = getReadChannel(longGridUrl); - longBuffer = fillBuffer(longChannel, HEADER_BYTES); - - //////////////////////// - //read header info - //////////////////////// - //skip the header description - latBuffer.position(latBuffer.position() + DESCRIPTION_LENGTH); - - int nc = latBuffer.getInt(); - int nr = latBuffer.getInt(); - int nz = latBuffer.getInt(); - - xmin = latBuffer.getFloat(); - dx = latBuffer.getFloat(); - ymin = latBuffer.getFloat(); - dy = latBuffer.getFloat(); - - float angle = latBuffer.getFloat(); - xmax = xmin + ((nc - 1) * dx); - ymax = ymin + ((nr - 1) * dy); - - //skip the longitude header description - longBuffer.position(longBuffer.position() + DESCRIPTION_LENGTH); - - //check that latitude grid header is the same as for latitude grid - if ( (nc != longBuffer.getInt()) - || (nr != longBuffer.getInt()) - || (nz != longBuffer.getInt()) - || (xmin != longBuffer.getFloat()) - || (dx != longBuffer.getFloat()) - || (ymin != longBuffer.getFloat()) - || (dy != longBuffer.getFloat()) - || (angle != longBuffer.getFloat())) { - throw new FactoryException(Errors.format(ErrorKeys.GRID_LOCATIONS_UNEQUAL)); - } - - //////////////////////// - //read grid shift data into LocalizationGrid - //////////////////////// - final int RECORD_LENGTH = (nc * 4) + SEPARATOR_BYTES; - final int NUM_BYTES_LEFT = ((nr + 1) * RECORD_LENGTH) - HEADER_BYTES; - final int START_OF_DATA = RECORD_LENGTH - HEADER_BYTES; - - latBuffer = fillBuffer(latChannel, NUM_BYTES_LEFT); - latBuffer.position(START_OF_DATA); //start of second record (data) - - longBuffer = fillBuffer(longChannel, NUM_BYTES_LEFT); - longBuffer.position(START_OF_DATA); - - gridShift = new LocalizationGrid(nc, nr); - - int i = 0; - int j = 0; - for (i = 0; i < nr; i++) { - latBuffer.position(latBuffer.position() + SEPARATOR_BYTES); //skip record separator - longBuffer.position(longBuffer.position() + SEPARATOR_BYTES); - - for (j = 0; j < nc; j++) { - gridShift.setLocalizationPoint(j, i, longBuffer.getFloat(), latBuffer.getFloat()); - } - } - - assert i == nr : i; - assert j == nc : j; - } - - /** - * Returns a new bytebuffer, of numBytes length and little endian byte - * order, filled from the channel. - * - * @param channel the channel to fill the buffer from - * @param numBytes number of bytes to read - * @return a new bytebuffer filled from the channel - * @throws IOException if there is a problem reading the channel - * @throws EOFException if the end of the channel is reached - */ - private ByteBuffer fillBuffer(ReadableByteChannel channel, int numBytes) - throws IOException { - ByteBuffer buf = ByteBuffer.allocate(numBytes); - - if (fill(buf, channel) == -1) { - throw new EOFException(Errors.format(ErrorKeys.END_OF_DATA_FILE)); - } - - buf.flip(); - buf.order(ByteOrder.LITTLE_ENDIAN); - - return buf; - } - - /** - * Fills the bytebuffer from the channel. Code was lifted from - * ShapefileDataStore. - * - * @param buffer bytebuffer to fill from the channel - * @param channel channel to fill the buffer from - * @return number of bytes read - * @throws IOException if there is a problem reading the channel - */ - private int fill(ByteBuffer buffer, ReadableByteChannel channel) - throws IOException { - int r = buffer.remaining(); - - // channel reads return -1 when EOF or other error - // because they a non-blocking reads, 0 is a valid return value!! - while ((buffer.remaining() > 0) && (r != -1)) { - r = channel.read(buffer); - } - - if (r == -1) { - buffer.limit(buffer.position()); - } - - return r; - } - - /** - * Obtain a ReadableByteChannel from the given URL. If the url protocol is - * file, a FileChannel will be returned. Otherwise a generic channel will - * be obtained from the urls input stream. Code swiped from - * ShapefileDataStore. - * - * @param url URL to create the channel from - * @return a new PeadableByteChannel from the input url - * @throws IOException if there is a problem creating the channel - */ - private ReadableByteChannel getReadChannel(URL url) - throws IOException { - ReadableByteChannel channel = null; - - if (url.getProtocol().equals("file")) { - File file = new File(url.getFile()); - - if (!file.exists() || !file.canRead()) { - throw new IOException(Errors.format(ErrorKeys.FILE_DOES_NOT_EXIST_$1, file)); - } - - FileInputStream in = new FileInputStream(file); - channel = in.getChannel(); - } else { - InputStream in = url.openConnection().getInputStream(); - channel = Channels.newChannel(in); - } - - return channel; - } - - /** - * Reads latitude and longitude text grid shift file data into {@link - * grid}. The first two lines of the shift data file contain the - * header, with the first being a description of the grid. The second line - * contains 8 values separated by spaces: num. columns, num. rows, num. z, - * min x, delta x, min y, delta y and angle. Shift data values follow this - * and are also separated by spaces. Row records are organized from low y - * (latitude) to high and columns are orderd from low longitude to high. - * - * @param latGridUrl URL to the text latitude shift file (.laa extention). - * @param longGridUrl URL to the text longitude shift file (.loa - * extention). - * @throws IOException if the data files cannot be read. - * @throws FactoryException if there is an inconsistency in the data - */ - private void loadTextGrid(URL latGridUrl, URL longGridUrl) - throws IOException, FactoryException { - String latLine; - String longLine; - StringTokenizer latSt; - StringTokenizer longSt; - - //////////////////////// - //setup - //////////////////////// - InputStreamReader latIsr = new InputStreamReader(latGridUrl.openStream()); - BufferedReader latBr = new BufferedReader(latIsr); - - InputStreamReader longIsr = new InputStreamReader(longGridUrl.openStream()); - BufferedReader longBr = new BufferedReader(longIsr); - - //////////////////////// - //read header info - //////////////////////// - latLine = latBr.readLine(); //skip header description - latLine = latBr.readLine(); - latSt = new StringTokenizer(latLine, " "); - - if (latSt.countTokens() != 8) { - throw new FactoryException(Errors.format(ErrorKeys.HEADER_UNEXPECTED_LENGTH_$1, - String.valueOf(latSt.countTokens()))); - } - - int nc = Integer.parseInt(latSt.nextToken()); - int nr = Integer.parseInt(latSt.nextToken()); - int nz = Integer.parseInt(latSt.nextToken()); - - xmin = Float.parseFloat(latSt.nextToken()); - dx = Float.parseFloat(latSt.nextToken()); - ymin = Float.parseFloat(latSt.nextToken()); - dy = Float.parseFloat(latSt.nextToken()); - - float angle = Float.parseFloat(latSt.nextToken()); - xmax = xmin + ((nc - 1) * dx); - ymax = ymin + ((nr - 1) * dy); - - //now read long shift grid - longLine = longBr.readLine(); //skip header description - longLine = longBr.readLine(); - longSt = new StringTokenizer(longLine, " "); - - if (longSt.countTokens() != 8) { - throw new FactoryException(Errors.format(ErrorKeys.HEADER_UNEXPECTED_LENGTH_$1, - String.valueOf(longSt.countTokens()))); - } - - //check that latitude grid header is the same as for latitude grid - if ( (nc != Integer.parseInt(longSt.nextToken())) - || (nr != Integer.parseInt(longSt.nextToken())) - || (nz != Integer.parseInt(longSt.nextToken())) - || (xmin != Float.parseFloat(longSt.nextToken())) - || (dx != Float.parseFloat(longSt.nextToken())) - || (ymin != Float.parseFloat(longSt.nextToken())) - || (dy != Float.parseFloat(longSt.nextToken())) - || (angle != Float.parseFloat(longSt.nextToken()))) { - throw new FactoryException(Errors.format(ErrorKeys.GRID_LOCATIONS_UNEQUAL)); - } - - //////////////////////// - //read grid shift data into LocalizationGrid - //////////////////////// - gridShift = new LocalizationGrid(nc, nr); - - int i = 0; - int j = 0; - for (i = 0; i < nr; i++) { - for (j = 0; j < nc;) { - latLine = latBr.readLine(); - latSt = new StringTokenizer(latLine, " "); - longLine = longBr.readLine(); - longSt = new StringTokenizer(longLine, " "); - - while (latSt.hasMoreTokens() && longSt.hasMoreTokens()) { - gridShift.setLocalizationPoint(j, i, - (double) Float.parseFloat(longSt.nextToken()), - (double) Float.parseFloat(latSt.nextToken())); - ++j; - } - } - } - - assert i == nr : i; - assert j == nc : j; - } - - /** * Transforms a list of coordinate point ordinal values. This method is * provided for efficiently transforming many points. The supplied array * of ordinal values will contain packed ordinal values. For example, if @@ -685,15 +332,14 @@ public class NADCONTransform extends AbstractMathTransform implements MathTransf double y = srcPts[srcOff++]; //check bounding box -//issue of bbox crossing +- 180 degrees (ie input point of -188 longitude); -//abs(x - xmin) > 0 , rollLongitude() ??? - if (((x < xmin) || (x > xmax)) || ((y < ymin) || (y > ymax))) { - throw new TransformException(Errors.format(ErrorKeys.POINT_OUTSIDE_GRID)); + if (((x < grid.getMinX()) || (x > grid.getMaxX())) || ((y < grid.getMinY()) || (y > grid.getMaxY()))) { + throw new TransformException("Point (" + x + " " + y + ") is not outside of ((" + grid.getMinX() + + " " + grid.getMinY() + ")(" + grid.getMaxX() + " " + grid.getMaxY() + "))"); } //find the grid the point is in (index is 0 based) - final double xgrid = (x - xmin) / dx; - final double ygrid = (y - ymin) / dy; + final double xgrid = (x - grid.getMinX()) / grid.getDx(); + final double ygrid = (y - grid.getMinY()) / grid.getDy(); double[] array = new double[] { xgrid, ygrid }; //use the LocalizationGridTransform2D transform method (bilineal interpolation) @@ -784,32 +430,13 @@ public class NADCONTransform extends AbstractMathTransform implements MathTransf return inverse; } - /** - * Returns a hash value for this transform. To make this faster it does not - * check the grid values. - * - * @return a hash value for this transform. - */ @Override - public final int hashCode() { - final long code = Double.doubleToLongBits(xmin) - + (37 * (Double.doubleToLongBits(ymin) - + (37 * (Double.doubleToLongBits(xmax) - + (37 * (Double.doubleToLongBits(ymax) - + (37 * (Double.doubleToLongBits(dx) - + (37 * (Double.doubleToLongBits(dy))))))))))); - - return (int) code ^ (int) (code >>> 32); + public int hashCode() { + return grid.hashCode(); } - - /** - * Compares the specified object with this math transform for equality. - * - * @param object the object to compare to - * @return {@code true} if the objects are equal. - */ + @Override - public final boolean equals(final Object object) { + public boolean equals(Object object) { if (object == this) { // Slight optimization return true; @@ -817,17 +444,11 @@ public class NADCONTransform extends AbstractMathTransform implements MathTransf if (super.equals(object)) { final NADCONTransform that = (NADCONTransform) object; - - return (Double.doubleToLongBits(this.xmin) == Double.doubleToLongBits(that.xmin)) - && (Double.doubleToLongBits(this.ymin) == Double.doubleToLongBits(that.ymin)) - && (Double.doubleToLongBits(this.xmax) == Double.doubleToLongBits(that.xmax)) - && (Double.doubleToLongBits(this.ymax) == Double.doubleToLongBits(that.ymax)) - && (Double.doubleToLongBits(this.dx) == Double.doubleToLongBits(that.dx)) - && (Double.doubleToLongBits(this.dy) == Double.doubleToLongBits(that.dy)) - && (this.gridShiftTransform).equals(that.gridShiftTransform); + + return this.grid.equals(that.grid); + } else { + return false; } - - return false; } /** @@ -964,14 +585,14 @@ public class NADCONTransform extends AbstractMathTransform implements MathTransf * parameter value. The default value is "conus.las". */ public static final ParameterDescriptor LAT_DIFF_FILE = new DefaultParameterDescriptor( - "Latitude_difference_file", String.class, null, "conus.las"); + "Latitude difference file", URI.class, null, null); /** * The operation parameter descriptor for the "Longitude_difference_file" * parameter value. The default value is "conus.los". */ public static final ParameterDescriptor LONG_DIFF_FILE = new DefaultParameterDescriptor( - "Longitude_difference_file", String.class, null, "conus.los"); + "Longitude difference file", URI.class, null, null); /** * The parameters group. @@ -1017,8 +638,8 @@ public class NADCONTransform extends AbstractMathTransform implements MathTransf throws ParameterNotFoundException, FactoryException { return new NADCONTransform( - stringValue(LAT_DIFF_FILE, values), - stringValue(LONG_DIFF_FILE, values)); + (URI) getParameter(LAT_DIFF_FILE, values).getValue(), + (URI) getParameter(LONG_DIFF_FILE, values).getValue()); } } } diff --git a/modules/library/referencing/src/main/java/org/geotools/referencing/operation/transform/NTv2Transform.java b/modules/library/referencing/src/main/java/org/geotools/referencing/operation/transform/NTv2Transform.java new file mode 100644 index 0000000..5548917 --- /dev/null +++ b/modules/library/referencing/src/main/java/org/geotools/referencing/operation/transform/NTv2Transform.java @@ -0,0 +1,409 @@ +/* + * GeoTools - The Open Source Java GIS Toolkit + * http://geotools.org + * + * (C) 2002-2012, Open Source Geospatial Foundation (OSGeo) + * + * This library is free software; you can redistribute it and/or + * modify it under the terms of the GNU Lesser General Public + * License as published by the Free Software Foundation; + * version 2.1 of the License. + * + * This library is distributed in the hope that it will be useful, + * but WITHOUT ANY WARRANTY; without even the implied warranty of + * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU + * Lesser General Public License for more details. + */ +package org.geotools.referencing.operation.transform; + + +import java.io.IOException; +import java.io.ObjectInputStream; +import java.io.Serializable; +import java.net.URI; +import java.net.URL; +import java.util.logging.Level; +import java.util.logging.Logger; + +import org.geotools.metadata.iso.citation.Citations; +import org.geotools.parameter.DefaultParameterDescriptor; +import org.geotools.parameter.Parameter; +import org.geotools.parameter.ParameterGroup; +import org.geotools.referencing.NamedIdentifier; +import org.geotools.referencing.ReferencingFactoryFinder; +import org.geotools.referencing.factory.IdentifiedObjectSet; +import org.geotools.referencing.factory.gridshift.GridShiftLocator; +import org.geotools.referencing.factory.gridshift.NTv2GridShiftFactory; +import org.geotools.referencing.operation.MathTransformProvider; +import org.geotools.util.logging.Logging; +import org.opengis.parameter.ParameterDescriptor; +import org.opengis.parameter.ParameterDescriptorGroup; +import org.opengis.parameter.ParameterNotFoundException; +import org.opengis.parameter.ParameterValue; +import org.opengis.parameter.ParameterValueGroup; +import org.opengis.referencing.FactoryException; +import org.opengis.referencing.NoSuchIdentifierException; +import org.opengis.referencing.operation.MathTransform; +import org.opengis.referencing.operation.MathTransform2D; +import org.opengis.referencing.operation.TransformException; +import org.opengis.referencing.operation.Transformation; + +import au.com.objectix.jgridshift.GridShift; +import au.com.objectix.jgridshift.GridShiftFile; + +/** + * The "NTv2" coordinate transformation method (EPSG:9615). + *

+ * This transformation depends on an external resource (the NTv2 grid file). If the file + * is not available, a {@link NoSuchIdentifierException recoverable NoSuchIdentifierException} + * will be thrown on instantiation. + * + * @see {@link IdentifiedObjectSet IdentifiedObjectSet exception handling}. + * @source $URL$ + * @version $Id$ + * @author Oscar Fonts + */ +public class NTv2Transform extends AbstractMathTransform implements MathTransform2D, Serializable { + + /** Serial number for interoperability with different versions. */ + private static final long serialVersionUID = -3082112044314062512L; + + /** Logger */ + protected static final Logger LOGGER = Logging.getLogger("org.geotools.referencing"); + + /** + * The original grid name + */ + private URI grid = null; + + /** The grid file name as set in the constructor. */ + private URL gridLocation = null; + + /** + * The grid shift to be used + */ + private GridShiftFile gridShift; + + /** + * The factory that loads the grid shift files + */ + private static NTv2GridShiftFactory FACTORY = new NTv2GridShiftFactory(); + + /** + * The inverse of this transform. Will be created only when needed. + */ + private transient MathTransform2D inverse; + + /** + * Constructs a {@code NTv2Transform} from the specified grid shift file. + * + * This constructor checks for grid shift file availability, but + * doesn't actually load the full grid into memory to preserve resources. + * + * @param file NTv2 grid file name + * @throws NoSuchIdentifierException if the grid is not available. + */ + public NTv2Transform(URI file) throws NoSuchIdentifierException { + if (file == null) { + throw new NoSuchIdentifierException("No NTv2 Grid File specified.", null); + } + + this.grid = file; + + gridLocation = locateGrid(grid.toString()); + if(gridLocation == null) { + throw new NoSuchIdentifierException("Could not locate NTv2 Grid File " + file, null); + } + + // Search for grid file + if (!FACTORY.isNTv2Grid(gridLocation)) { + throw new NoSuchIdentifierException("NTv2 Grid File not available.", + file.toString()); + } + } + + URL locateGrid(String grid) { + for (GridShiftLocator locator : ReferencingFactoryFinder.getGridShiftLocators(null)) { + URL result = locator.locateGrid(grid); + if(result != null) { + return result; + } + }; + + return null; + } + + /** + * Returns the inverse of this transform. + * + * @return the inverse of this transform + */ + @Override + public synchronized MathTransform2D inverse() { + if (inverse == null) { + inverse = new Inverse(); + } + return inverse; + } + + /** + * Transforms a list of coordinate point ordinal values. This method is + * provided for efficiently transforming many points. The supplied array + * of ordinal values will contain packed ordinal values. For example, if + * the source dimension is 3, then the ordinals will be packed in this + * order: + * (x0,y0,z0, + * + * x1,y1,z1 + * ...). + * + * @param srcPts the array containing the source point coordinates. + * @param srcOff the offset to the first point to be transformed in the + * source array. + * @param dstPts the array into which the transformed point coordinates are + * returned. May be the same than {@code srcPts}. + * @param dstOff the offset to the location of the first transformed point + * that is stored in the destination array. + * @param numPts the number of point objects to be transformed. + * + * @throws TransformException if an IO error occurs reading the grid file. + */ + @Override + public void transform(double[] srcPts, int srcOff, double[] dstPts, + int dstOff, int numPts) throws TransformException { + bidirectionalTransform(srcPts,srcOff, dstPts, dstOff, numPts, true); + } + + /** + * Inverse transform. See {@link #transform(double[], int, double[], + * int, int)} + * + * @param srcPts the array containing the source point coordinates. + * @param srcOff the offset to the first point to be transformed in the + * source array. + * @param dstPts the array into which the transformed point coordinates are + * returned. May be the same than {@code srcPts}. + * @param dstOff the offset to the location of the first transformed point + * that is stored in the destination array. + * @param numPts the number of point objects to be transformed. + * + * @throws TransformException if an IO error occurs reading the grid file. + */ + public void inverseTransform(double[] srcPts, int srcOff, double[] dstPts, + int dstOff, int numPts) throws TransformException { + bidirectionalTransform(srcPts,srcOff, dstPts, dstOff, numPts, false); + } + + /** + * Performs the actual transformation. + * + * @param srcPts the array containing the source point coordinates. + * @param srcOff the offset to the first point to be transformed in the + * source array. + * @param dstPts the array into which the transformed point coordinates are + * returned. May be the same than {@code srcPts}. + * @param dstOff the offset to the location of the first transformed point + * that is stored in the destination array. + * @param numPts the number of point objects to be transformed. + * @param forward {@code true} for direct transform, {@code false} for inverse transform. + * + * @throws TransformException if an IO error occurs reading the grid file. + */ + private void bidirectionalTransform(double[] srcPts, int srcOff, double[] dstPts, + int dstOff, int numPts, boolean forward) throws TransformException { + + boolean shifted; + + if (gridShift == null) { // Create grid when first needed. + try { + gridShift = FACTORY.createNTv2Grid(gridLocation); + } catch (FactoryException e) { + throw new TransformException("NTv2 Grid " + gridLocation + + " Could not be created", e); + } + } + + try { + GridShift shift = new GridShift(); + for (int i=0; i file = new Parameter(Provider.FILE); + file.setValue(grid); + + return new ParameterGroup(Provider.PARAMETERS, + new ParameterValue[] { file } + ); + } + + /** + * Inverse of a {@link NTv2Transform}. + * + * @version $Id$ + * @author Oscar Fonts + */ + private final class Inverse extends AbstractMathTransform.Inverse + implements MathTransform2D, Serializable + { + /** Serial number for interoperability with different versions. */ + private static final long serialVersionUID = -4707304160205218546L; + + /** + * Default constructor. + */ + public Inverse() { + NTv2Transform.this.super(); + } + + /** + * Returns the parameter values for this math transform. + * + * @return A copy of the parameter values for this math transform. + */ + @Override + public ParameterValueGroup getParameterValues() { + return null; + } + + /** + * Inverse transform an array of points. + * + * @param source + * @param srcOffset + * @param dest + * @param dstOffset + * @param length + * + * @throws TransformException if the input point is outside the area + * covered by this grid. + */ + public void transform(final double[] source, final int srcOffset, + final double[] dest, final int dstOffset, final int length) + throws TransformException { + NTv2Transform.this.inverseTransform(source, srcOffset, dest, + dstOffset, length); + } + + /** + * Returns the original transform. + */ + @Override + public MathTransform2D inverse() { + return (MathTransform2D) super.inverse(); + } + + /** + * Restore reference to this object after deserialization. + * + * @param in DOCUMENT ME! + * @throws IOException DOCUMENT ME! + * @throws ClassNotFoundException DOCUMENT ME! + */ + private void readObject(ObjectInputStream in) + throws IOException, ClassNotFoundException { + in.defaultReadObject(); + NTv2Transform.this.inverse = this; + } + } + + /** + * The {@link NTv2Transform} provider. + * + * @author Oscar Fonts + */ + public static class Provider extends MathTransformProvider { + + private static final long serialVersionUID = -3710592152744574801L; + + /** + * The operation parameter descriptor for the "Latitude and longitude difference file" + * parameter value. The default value is "". + */ + public static final DefaultParameterDescriptor FILE = new DefaultParameterDescriptor( + toMap(new NamedIdentifier[] { + new NamedIdentifier(Citations.EPSG, "Latitude and longitude difference file"), + new NamedIdentifier(Citations.EPSG, "8656") + }), + URI.class, null, null, null, null, null, true); + + /** + * The parameters group. + */ + static final ParameterDescriptorGroup PARAMETERS = createDescriptorGroup(new NamedIdentifier[] { + new NamedIdentifier(Citations.EPSG, "NTv2"), + new NamedIdentifier(Citations.EPSG, "9615") + }, new ParameterDescriptor[] { + FILE + }); + + /** + * Constructs a provider. + */ + public Provider() { + super(2, 2, PARAMETERS); + } + + /** + * Returns the operation type. + */ + @Override + public Class getOperationType() { + return Transformation.class; + } + + /** + * Creates a math transform from the specified group of parameter + * values. + * + * @param values The group of parameter values. + * @return The created math transform. + * @throws ParameterNotFoundException if a required parameter was not + * found. + * @throws FactoryException if there is a problem creating this + * math transform. + */ + protected MathTransform createMathTransform(final ParameterValueGroup values) + throws ParameterNotFoundException, FactoryException + { + return new NTv2Transform(value(FILE, values)); + } + } +} diff --git a/modules/library/referencing/src/main/resources/META-INF/services/org.geotools.referencing.factory.gridshift.GridShiftLocator b/modules/library/referencing/src/main/resources/META-INF/services/org.geotools.referencing.factory.gridshift.GridShiftLocator new file mode 100644 index 0000000..b0a85b9 --- /dev/null +++ b/modules/library/referencing/src/main/resources/META-INF/services/org.geotools.referencing.factory.gridshift.GridShiftLocator @@ -0,0 +1 @@ +org.geotools.referencing.factory.gridshift.ClasspathGridShiftLocator \ No newline at end of file diff --git a/modules/library/referencing/src/main/resources/META-INF/services/org.geotools.referencing.operation.MathTransformProvider b/modules/library/referencing/src/main/resources/META-INF/services/org.geotools.referencing.operation.MathTransformProvider index b190d22..eef231d 100644 --- a/modules/library/referencing/src/main/resources/META-INF/services/org.geotools.referencing.operation.MathTransformProvider +++ b/modules/library/referencing/src/main/resources/META-INF/services/org.geotools.referencing.operation.MathTransformProvider @@ -10,6 +10,7 @@ org.geotools.referencing.operation.transform.GeocentricTransform$ProviderInverse org.geotools.referencing.operation.transform.MolodenskiTransform$Provider org.geotools.referencing.operation.transform.MolodenskiTransform$ProviderAbridged org.geotools.referencing.operation.transform.NADCONTransform$Provider +org.geotools.referencing.operation.transform.NTv2Transform$Provider org.geotools.referencing.operation.transform.SimilarityTransformProvider org.geotools.referencing.operation.transform.WarpTransform2D$Provider org.geotools.referencing.operation.projection.EquidistantCylindrical$Provider diff --git a/modules/library/referencing/src/test/java/org/geotools/referencing/factory/gridshift/NTv2GridShiftFactoryTest.java b/modules/library/referencing/src/test/java/org/geotools/referencing/factory/gridshift/NTv2GridShiftFactoryTest.java new file mode 100644 index 0000000..7e26600 --- /dev/null +++ b/modules/library/referencing/src/test/java/org/geotools/referencing/factory/gridshift/NTv2GridShiftFactoryTest.java @@ -0,0 +1,111 @@ +/* + * GeoTools - The Open Source Java GIS Toolkit + * http://geotools.org + * + * (C) 2002-2012, Open Source Geospatial Foundation (OSGeo) + * + * This library is free software; you can redistribute it and/or + * modify it under the terms of the GNU Lesser General Public + * License as published by the Free Software Foundation; + * version 2.1 of the License. + * + * This library is distributed in the hope that it will be useful, + * but WITHOUT ANY WARRANTY; without even the implied warranty of + * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU + * Lesser General Public License for more details. + */ +package org.geotools.referencing.factory.gridshift; + +import static org.junit.Assert.*; + +import java.net.URL; + +import org.junit.Before; +import org.junit.Test; +import org.opengis.referencing.FactoryException; + +import au.com.objectix.jgridshift.GridShiftFile; + + +/** + * Unit tests for {@link NTv2GridShiftFactory} public methods + * + * @author Oscar Fonts + */ +public class NTv2GridShiftFactoryTest { + + private static final URL TEST_GRID = NTv2GridShiftFactoryTest.class.getResource("BALR2009.gsb"); + private static final URL INEXISTENT_GRID = NTv2GridShiftFactoryTest.class.getResource("this_NTv2_grid_does_not_exist"); + private static final URL MALFORMED_GRID = NTv2GridShiftFactoryTest.class.getResource("malformedNTv2grid.gsb"); + + private NTv2GridShiftFactory factory; + + /** + * Instantiates the test {@link NTv2GridShiftFactory} + * + * @throws java.lang.Exception + */ + @Before + public void setUp() throws Exception { + factory = new NTv2GridShiftFactory(); + } + + + /** + * Test method for {@link org.geotools.referencing.factory.gridshift.NTv2GridShiftFactory#isNTv2Grid(java.lang.String)}. + */ + @Test + public void testIsNTv2GridAvailable() { + assertFalse(factory.isNTv2Grid(null)); + assertFalse(factory.isNTv2Grid(INEXISTENT_GRID)); // Will log a FINEST message + assertFalse(factory.isNTv2Grid(MALFORMED_GRID)); // Will log a WARNING message + assertTrue(factory.isNTv2Grid(TEST_GRID)); + } + + /** + * Test method for {@link org.geotools.referencing.factory.gridshift.NTv2GridShiftFactory#createNTv2Grid(java.lang.String)}. + */ + @Test + public void testCreateNTv2Grid() { + + // Try to create a null grid + boolean factoryExceptionThrown = false; + try { + factory.createNTv2Grid(null); + } catch (FactoryException e) { + factoryExceptionThrown = true; + } + assertTrue(factoryExceptionThrown); + + // Try to create a grid from an inexistent resource + factoryExceptionThrown = false; + try { + factory.createNTv2Grid(INEXISTENT_GRID); // Will log a FINEST message + } catch (FactoryException e) { + factoryExceptionThrown = true; + } + assertTrue(factoryExceptionThrown); + + // Try to create a grid from a corrupt NTv2 grid file + factoryExceptionThrown = false; + try { + factory.createNTv2Grid(MALFORMED_GRID); // Will log a SEVERE message + } catch (FactoryException e) { + factoryExceptionThrown = true; + } + assertTrue(factoryExceptionThrown); + + // Create a grid from the test file + factoryExceptionThrown = false; + GridShiftFile grid = null; + try { + grid = factory.createNTv2Grid(TEST_GRID); + } catch (FactoryException e) { + factoryExceptionThrown = true; + } + assertFalse(factoryExceptionThrown); + assertNotNull(grid); + assertTrue(grid.isLoaded()); + } + +} diff --git a/modules/library/referencing/src/test/java/org/geotools/referencing/factory/gridshift/TestGridShiftLocator.java b/modules/library/referencing/src/test/java/org/geotools/referencing/factory/gridshift/TestGridShiftLocator.java new file mode 100644 index 0000000..75e4d2c --- /dev/null +++ b/modules/library/referencing/src/test/java/org/geotools/referencing/factory/gridshift/TestGridShiftLocator.java @@ -0,0 +1,68 @@ +package org.geotools.referencing.factory.gridshift; + +import java.io.File; +import java.io.FileOutputStream; +import java.io.IOException; +import java.net.URL; +import java.util.logging.Level; +import java.util.logging.Logger; +import java.util.zip.GZIPInputStream; + +import org.geotools.metadata.iso.citation.Citations; +import org.geotools.util.logging.Logging; +import org.opengis.metadata.citation.Citation; + +/** + * Test helper that can load compressed grids + */ +public class TestGridShiftLocator implements GridShiftLocator { + + static final Logger LOGGER = Logging.getLogger(TestGridShiftLocator.class); + + @Override + public Citation getVendor() { + return Citations.GEOTOOLS; + } + + @Override + public URL locateGrid(String grid) { + GZIPInputStream is = null; + FileOutputStream fos = null; + try { + URL compressed = getClass().getResource(grid + ".gz"); + if (compressed != null) { + is = new GZIPInputStream(compressed.openStream()); + File out = new File("./target/" + grid); + fos = new FileOutputStream(out); + byte[] buf = new byte[1024]; + int read; + while ((read = is.read(buf)) > 0) { + fos.write(buf, 0, read); + } + + return out.toURI().toURL(); + } + } catch (IOException e) { + LOGGER.log(Level.SEVERE, "Failed to unpack the grid", e); + } finally { + try { + if (is != null) { + is.close(); + } + } catch (IOException e) { + // + } + + try { + if (fos != null) { + fos.close(); + } + } catch (IOException e) { + // + } + } + + return null; + } + +} diff --git a/modules/library/referencing/src/test/java/org/geotools/referencing/operation/transform/NADCONTransformTest.java b/modules/library/referencing/src/test/java/org/geotools/referencing/operation/transform/NADCONTransformTest.java new file mode 100644 index 0000000..4d056ae --- /dev/null +++ b/modules/library/referencing/src/test/java/org/geotools/referencing/operation/transform/NADCONTransformTest.java @@ -0,0 +1,148 @@ +/* + * GeoTools - The Open Source Java GIS Toolkit + * http://geotools.org + * + * (C) 2002-2012, Open Source Geospatial Foundation (OSGeo) + * + * This library is free software; you can redistribute it and/or + * modify it under the terms of the GNU Lesser General Public + * License as published by the Free Software Foundation; + * version 2.1 of the License. + * + * This library is distributed in the hope that it will be useful, + * but WITHOUT ANY WARRANTY; without even the implied warranty of + * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU + * Lesser General Public License for more details. + */ +package org.geotools.referencing.operation.transform; + +import static org.junit.Assert.*; + +import java.net.URI; + +import org.junit.Before; +import org.junit.Test; +import org.opengis.parameter.ParameterValueGroup; +import org.opengis.referencing.NoSuchIdentifierException; +import org.opengis.referencing.operation.TransformException; + +/** + * Unit tests for {@link NADCONTransform} public methods + * + * @author Andrea Aime - GeoSolutions + */ +public class NADCONTransformTest { + + private static final String STPAUL_LAS = "stpaul.las"; + private static final String STPAUL_LOS = "stpaul.los"; + private static final double TOLERANCE = 18E-7; + + private static final double[] TEST_POINT_SRC = {-170.25, 57.125}; + private static final double[] TEST_POINT_DST = {-170.252283, 57.125655}; + + + private NADCONTransform transform; + + /** + * Instantiates the test {@link NADCONTransform} + * + * @throws java.lang.Exception + */ + @Before + public void setUp() throws Exception { + transform = new NADCONTransform(new URI(STPAUL_LAS), new URI(STPAUL_LOS)); + } + + /** + * Test method for {@link org.geotools.referencing.operation.transform.NADCONTransform#getSourceDimensions()}. + */ + @Test + public void testGetSourceDimensions() { + assertEquals(transform.getSourceDimensions(), 2); + } + + /** + * Test method for {@link org.geotools.referencing.operation.transform.NADCONTransform#getTargetDimensions()}. + */ + @Test + public void testGetTargetDimensions() { + assertEquals(transform.getTargetDimensions(), 2); + } + + /** + * Test method for {@link org.geotools.referencing.operation.transform.NADCONTransform#getParameterValues()}. + */ + @Test + public void testGetParameterValues() { + ParameterValueGroup pvg = transform.getParameterValues(); + + // Descriptor is the same as provider's + assertEquals(pvg.getDescriptor(), NADCONTransform.Provider.PARAMETERS); + + // One single value + assertEquals(pvg.values().size(), 2); + + Object value = pvg.parameter("Latitude difference file").getValue(); + assertTrue(value instanceof URI); + assertEquals(value.toString(), STPAUL_LAS); + + value = pvg.parameter("Longitude difference file").getValue(); + assertTrue(value instanceof URI); + assertEquals(value.toString(), STPAUL_LOS); + } + + /** + * Test method for {@link org.geotools.referencing.operation.transform.NADCONTransform#NADCONTransform(java.net.URI)}. + */ + @Test + public void testNADCONTransform() throws Exception { + + try { + new NADCONTransform(null, null); + } catch (NoSuchIdentifierException e) { + assert true; + } + } + + /** + * Test method for {@link org.geotools.referencing.operation.transform.NADCONTransform#inverse()}. + * @throws TransformException + */ + @Test + public void testInverse() throws TransformException { + assertSame(transform, transform.inverse().inverse()); + + double[] p1 = new double[2]; + double[] p2 = new double[2]; + + transform.inverse().transform(TEST_POINT_DST, 0, p1, 0, 1); + transform.inverseTransform(TEST_POINT_DST, 0, p2, 0, 1); + assertEquals(p1[0], p2[0], TOLERANCE); + assertEquals(p1[1], p2[1], TOLERANCE); + } + + /** + * Test method for {@link org.geotools.referencing.operation.transform.NADCONTransform#transform(double[], int, double[], int, int)}. + * @throws TransformException + */ + @Test + public void testTransform() throws TransformException { + double[] p = new double[2]; + transform.transform(TEST_POINT_SRC, 0, p, 0, 1); + assertEquals(p[0], TEST_POINT_DST[0], TOLERANCE); + assertEquals(p[1], TEST_POINT_DST[1], TOLERANCE); + } + + /** + * Test method for {@link org.geotools.referencing.operation.transform.NADCONTransform#inverseTransform(double[], int, double[], int, int)}. + * @throws TransformException + */ + @Test + public void testInverseTransform() throws TransformException { + double[] p = new double[2]; + transform.inverseTransform(TEST_POINT_DST, 0, p, 0, 1); + assertEquals(p[0], TEST_POINT_SRC[0], TOLERANCE); + assertEquals(p[1], TEST_POINT_SRC[1], TOLERANCE); + } + +} diff --git a/modules/library/referencing/src/test/java/org/geotools/referencing/operation/transform/NTv2TransformTest.java b/modules/library/referencing/src/test/java/org/geotools/referencing/operation/transform/NTv2TransformTest.java new file mode 100644 index 0000000..1f53afb --- /dev/null +++ b/modules/library/referencing/src/test/java/org/geotools/referencing/operation/transform/NTv2TransformTest.java @@ -0,0 +1,153 @@ +/* + * GeoTools - The Open Source Java GIS Toolkit + * http://geotools.org + * + * (C) 2002-2012, Open Source Geospatial Foundation (OSGeo) + * + * This library is free software; you can redistribute it and/or + * modify it under the terms of the GNU Lesser General Public + * License as published by the Free Software Foundation; + * version 2.1 of the License. + * + * This library is distributed in the hope that it will be useful, + * but WITHOUT ANY WARRANTY; without even the implied warranty of + * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU + * Lesser General Public License for more details. + */ +package org.geotools.referencing.operation.transform; + +import static org.junit.Assert.*; + +import java.net.URI; + +import org.junit.Before; +import org.junit.Test; +import org.opengis.parameter.ParameterValueGroup; +import org.opengis.referencing.NoSuchIdentifierException; +import org.opengis.referencing.operation.TransformException; + +/** + * Unit tests for {@link NTv2Transform} public methods + * + * @author Oscar Fonts + */ +public class NTv2TransformTest { + + private static final String TEST_GRID = "BALR2009.gsb"; + private static final String INEXISTENT_GRID = "this_NTv2_grid_does_not_exist"; + private static final double[] TEST_POINT_SRC = {3.084896111, 39.592654167}; + private static final double[] TEST_POINT_DST = {3.083801819, 39.5914804}; + private static final double TOLERANCE = 18E-7; + + private NTv2Transform transform; + + /** + * Instantiates the test {@link NTv2Transform} + * + * @throws java.lang.Exception + */ + @Before + public void setUp() throws Exception { + transform = new NTv2Transform(new URI(TEST_GRID)); + } + + /** + * Test method for {@link org.geotools.referencing.operation.transform.NTv2Transform#getSourceDimensions()}. + */ + @Test + public void testGetSourceDimensions() { + assertEquals(transform.getSourceDimensions(), 2); + } + + /** + * Test method for {@link org.geotools.referencing.operation.transform.NTv2Transform#getTargetDimensions()}. + */ + @Test + public void testGetTargetDimensions() { + assertEquals(transform.getTargetDimensions(), 2); + } + + /** + * Test method for {@link org.geotools.referencing.operation.transform.NTv2Transform#getParameterValues()}. + */ + @Test + public void testGetParameterValues() { + ParameterValueGroup pvg = transform.getParameterValues(); + + // Descriptor is the same as provider's + assertEquals(pvg.getDescriptor(), NTv2Transform.Provider.PARAMETERS); + + // One single value + assertEquals(pvg.values().size(), 1); + + // Value accessible through its identifiers + Object value = pvg.parameter("8656").getValue(); + assertTrue(value instanceof URI); + assertEquals(value.toString(), TEST_GRID); + + value = pvg.parameter("Latitude and longitude difference file").getValue(); + assertTrue(value instanceof URI); + assertEquals(value.toString(), TEST_GRID); + } + + /** + * Test method for {@link org.geotools.referencing.operation.transform.NTv2Transform#NTv2Transform(java.net.URI)}. + */ + @Test + public void testNTv2Transform() throws Exception { + + try { + new NTv2Transform(null); + } catch (NoSuchIdentifierException e) { + assert true; + } + + try { + new NTv2Transform(new URI(INEXISTENT_GRID)); + } catch (NoSuchIdentifierException e) { + return; + } + } + + /** + * Test method for {@link org.geotools.referencing.operation.transform.NTv2Transform#inverse()}. + * @throws TransformException + */ + @Test + public void testInverse() throws TransformException { + assertSame(transform, transform.inverse().inverse()); + + double[] p1 = new double[2]; + double[] p2 = new double[2]; + + transform.inverse().transform(TEST_POINT_SRC, 0, p1, 0, 1); + transform.inverseTransform(TEST_POINT_SRC, 0, p2, 0, 1); + assertEquals(p1[0], p2[0], TOLERANCE); + assertEquals(p1[1], p2[1], TOLERANCE); + } + + /** + * Test method for {@link org.geotools.referencing.operation.transform.NTv2Transform#transform(double[], int, double[], int, int)}. + * @throws TransformException + */ + @Test + public void testTransform() throws TransformException { + double[] p = new double[2]; + transform.transform(TEST_POINT_SRC, 0, p, 0, 1); + assertEquals(p[0], TEST_POINT_DST[0], TOLERANCE); + assertEquals(p[1], TEST_POINT_DST[1], TOLERANCE); + } + + /** + * Test method for {@link org.geotools.referencing.operation.transform.NTv2Transform#inverseTransform(double[], int, double[], int, int)}. + * @throws TransformException + */ + @Test + public void testInverseTransform() throws TransformException { + double[] p = new double[2]; + transform.inverseTransform(TEST_POINT_DST, 0, p, 0, 1); + assertEquals(p[0], TEST_POINT_SRC[0], TOLERANCE); + assertEquals(p[1], TEST_POINT_SRC[1], TOLERANCE); + } + +} diff --git a/modules/library/referencing/src/test/resources/org/geotools/referencing/factory/gridshift/BALR2009.gsb b/modules/library/referencing/src/test/resources/org/geotools/referencing/factory/gridshift/BALR2009.gsb new file mode 100755 index 0000000000000000000000000000000000000000..5aa0b70a91ce3e93416ccbef24b48aa5f273906c GIT binary patch literal 101552 zcmZ6!1yq$?xUP-eVvB{?fCY975;mZK*qvCYs9^V6>s{}=7Am%aEf%(-qCOTT7}(wU z+1>x`KKppj^Y1bC8ROb>!RO~U&h^ZBKhK;iAhd0d!1jJVxg7q_{{(=8|9|Y?vW;Jc z|6Ok$+#{q@ke^GipHEO zmqYwQ+6ULG@8Z(dyH#NO9{w(Wf8FDdj^1rs2K*0aeyaa`*%bUg>mh&Fz1IF)4{>=| zr+A@bSPu?u)+4~XtzR?mHh$ji{epw={{)1%xC8_S`2GFo`Ly@*4)ODK@$~TU#Gj{T zXpry!@V`G_k2c=;_}|~};IPBNJ3!%luQGW1^=K0q(ENYb2cDNOl_S;(e!54?0H6Q! z^U~Yv@2CI%AO7d}f${g_6Brs0veW5*{x1&Q-dV*XZJ?Kf!<0jR|F_T&tLPO6{o~V_ zb-6N$JTb2>*{_nTIqBk7J^bJ2Zdc?oaYhQFZu$(+2qalHx$d){Y>~1;o-{FpW)@8`!9w#mTp zBBA6Zt!y%28CmOSlO>PH$^C88wet9XKVNmUO^%Kr zU!QK1n|sIy=GbKCFLK%nn{@J-@bBldZ?nlGGr5P@q+1&KgH zRcgkOpFf80Ka+RAvC5eM`g+ngt9)Ba_RgeB)<5J7dEn15`ucDYT^=Qpr@QE~D=Yon z#vMNFPA*wXm(07#UmNQ(B@g|4uemM{`;xD<)1_Q8`2p&IcR~93H|jynaB{T>T{0ae zj~u5(ccYdG8`!1{bHF+pgB7c?9`GqAstFlQ-_wr9%n&dBibY ztP$kq=fKCwZZ~uZC{8~Yd92IDVdOk-@$Y+-?D#{MDn;pMr>us23n#ngH)PU5vR6q% zcIKy_yHzq|RbTS-8iq93L*Crjkfv6lMTr>ioETlA>R*@ z3tcwEHz)nP{=Ok!LdYXu81X#M@NqQxZaJGQK1{ygX_MdC z>EoYGZSp#h+&#c1Gv|{(cD2d=TjV{1ZL+;A{k&?NP5KNZZ?@Rv?PhY?LYs8?K>oQ7 z^`IX8Ty2j{icTVTPO*t`lsxvbO&a8&pT|G4NxeXFuMaj^G>=?1qg|}G$VvI^;#!`5 zb}4O_*dgSuRqgU&7kOYq_|}2GZQar?c`A}K!T;7aUuchtD@7-`a_KFq9l}#4ZgNk)NNl%h==O6}RnD;yrodOS=>* zP5(~$W|z=5c&iM^O+WWtZcQ{Z1fLQXiS-J$iAdjDo{drUR zy2?^rM$99>SgXs57vx>rbeY?fzP2Ufyb((tenOYc&&b^_>TRu$Ttz(FDIr=$o zGegXy$XZ)N#vCO(cQWL10s8q=UqiO^BF~P%pJxj>3g>|}zsS>Y9yrp9zTOgRi2nj| z_7#Ryen_@&Mt!PAKl|=8q-r$T{e&T{kCOc^8FC{x{k#DFzp3oUUEUbdA&%>Ju>TKv zLcW;OCU>gR&p(}P@+p!WT*)RY_LHyI#(D<&c!sY{zBVTZ1>0nkL9W@`CS}se6-U@) zuoL}Uc8X2R-N+^EHkr4QoIB1YE1!~IY_iFbD)jT#{Wft@_TwE-+horkuHU?7le53b z^`F`#+Lyi_{t4#?D|tdDyVSo#?wQ{%7xU7$rQ!dQwaELc*riBkay{(-k4_=SHMfgr z0{LwS_V?$>b$Z*S?KkrH;dc36mi|qdV3**wWETyek0Q@rV3#Ya$j)o*^6(rvZM$81 z{3NdxyVP-`kJqKyrCSGb`YpShnoh0?|A%iQ?|}ck?~(^#|L>iLKJPrnD*G$b*VCq0 zW!?<(O3f;5Zj*P!TIEG;`g%{CRmR%L+Y+qu@(Fq6cB=$5rLS$tR#~xtJSN2|kKT~G zp0`Tj7W8#9)B~^OeokzuOaI>F(jmH}C6SkR$N9vOey%!5m)CvC zCr9ZLw~rh#33*OF`niGy`Bp#jlexO&-ACRTr_18J^s{ZfF86wqNAE;_wwpZokS>uq z>F2?x@%Ph>JpQsSOSh62-P6T66a9SRCF(&aIV(ykw3cKbMi?;|-~(?8nnLArDga;|<{dyW{BV3-JH&qvWO7|LeKv z=S2?-g7jT+z)xr%TTqXj>F28l zY|^MN`Q{n?^EZ%_uG=Kx9eFhT|GFf7>j?j^@FrWa|Ia>v{L;}bRSa_NQg(6PNgnEM zm%BH~%j()?SZ4Y*)z>a3tCBOfho8g9ZF|@yVH$b!5W8ecBv&10m$%o+7fp5v%u4^p z&$COun&ic+kl*wqAK7M?X$ILD{x7+oycPa``;0sZ`~N>h=<}eqR;kgNzFr-Md}THH zNiVD9|3xk`#44$w^mSd-gXlzZbJT-e+34p0)Ps3F$SqM1oRY}RP!Fv6=;u192cHL$ z-B1rY9VWY=9_%VcUl&6?_%cFa)C0HEWEa!}@3Qo@C+dOMIC5*$gOV4?BTx@6yV2Jx zP!IY{CEr9nIB<TJ&|YmoAMovR7+es@@}? z57p&bZTfn6Pvk>ZveRJr^d9;AXkFgdqOa#o(WQrl9B##a`!+dXzAlA4>1*#5$hT&Y zTW->&%~f)5^b0Py(a*Lc$hXFm&!u92pGtPUuFFbi`gzGCU4n*_DLuM=acnSF5JB+?w0RInOPWD}HNU`T+=R`wVdD742k^f|mCWo9h zq@9pAApiHsMjw}cYRHndnL()IPxWBKW>e(Nu?iLUp@=mlzyHYYZGlcx%Wz&#GE7- z*lLp++34Ff@PE8J*$Mtn3?Prd{y$R5S9UzMNzO%Fzw;jd{)fol|JWq$2{||Ne+NhU zw^VWDA@#`hD%s^-FYaPMPG6s?W|dz@$)5GC^1UQ|J<-Q1 zmqwG5P!Hy%k?)`$1XQB0KcgP}nM}@vda(W)IS1-Ns3(1$9rfUwg`D+oJ$OLQjCxSB zA$^?*_2BSaa#qxX#{ZEWQ4cow($~(Y2S1mRYoZ?1{zMK!J@9EuUynsSsJ4c@4)x&C zPx4*VgZ?4(bs6*zHf$shM*m=2M*8^(`Um4Xk;|Z;{A&x@D2sk_R{A+pRb2{qCC{v@ z%ld8PLd~$h&q_aUZv)?UA@>SHzQ2|1(p#6And#?uLv-=&NWL@{`OZf2xoNt@{w6;# zbjjI{zRrn$L6KGDzx{$WAIY1y=#u0|U%T$t#cw|ONDBG~Psrnu4_&WKKZo8yUNnOo z{6d$u>Ew~0b*Wv3e%^!gK+G_5IY;<-KY4cvLpCb=@sY0Zud*NSRNIhcarF5R3Z9)P@m_eyfsXq)&yAzR`9_J!%&_j7HM zxgOaE`~Qwz$kB;589kZ25c$vG)#OztZ4!K%yyUV?>bxUcAK0XDG5UAx8=E}vCJ+8; zli0!J(K+mre?ECx5xX=xOny?%E?HiYyL;NDbus!kYg4=Y_9NeIW0wr0$Zt9$&sk2c zH^?qmQ_081+GWiT@;dl`j~jj7>4zp4H2OP3R;y%rL^kEKO5Vow_1ofB`8|*9RlzFh zugT-8TV-Q&`g$em!T9CmWYmMuZ{$?e1K0NS^%c~Eyc@|kQ4jvhL_gmiY?ZTJ$@fQF zW%Dlbqsdk=<)xorSgg`~0Qu8g?B@@W|JDP;iGFrMJ#ZL7_CP&oaf%#-deF5DeLWiW zpg}Zw1?s`$^W@8@2Q6Ib>wM@Z_n%A-MnAdn74mY_gZ)+L>o2GWFQ$>ZV7_S2HS!tM zgWA>T>$c^Pr%xx}MLqa_l{}*!@}(;Db5kGGiz#Hs09_`hlYex?d{9ODIa6QMk0^3Q z%oDvlO&*B;@!*p5^VylmpN5iy&@Z@kfczQ#f}?rp=OYQaRO?3GiF_!7vL8P!x`ZkF z@$9LnC(3?&2+jkYl>PWO^o#R6pwG9w*QKG7kH`E$e_F}M7w5))`Z)c45qZ!_NBX%& zCG=-Ilckm+M^}-f;eV%Bhmh|t zFr>#4^60gO^t(^azsr!@73t?qM-4eQnB3r^Aue0VyY3nC@B_Ip^8P$^>F1DNhOCb! z56xwhQ~Sx?i`itNl8?I~@Bi6^eon82`JgG}F3oK6@i6%e{Qo@z{aOzG4=zXUkNtlw zKXU9?n;cT|mG#JfN?EzSVYW@ew~+0~|BqcF_t}Ww_k&z13H@ayFMoW@CV{Q#;|&*V zawC#Fg`&5&_ZsZZ~HSyj}zWh^@s*d!t zb2ju-2a-ecTP53J^0*RK$x)oXo{f4?cocah>OrluLMovOK@KE}l`%w?7G@`GMpdMtIM?Qsmko<~#1@)kbKYjgZ8RmWB z$RF2PWv()Rm=pD&TL67s5%u8S8gdKNgIvGJBT)}-2h-Q9Q4hj4kRRfFFx-KD_CP(T z*nw=I9t_$@{)u|f%z=KMj(V^ygj}W=^7aJsP1FO&pX43YbjcA&UvH?7{rw8^F<)Iq zeIREJg0Fq)>w%ajIy8^`sXy`^Wk0@tq%Jp<{kR$Zf&w$?^B7H+qoi_g}e3 zE{wc?Y#I7_({Do(dyz}!w#l*;mQK+SKLQVinPh|JLH~IvHw=`^7jVz=T+$6sKv-{LdbvC;`0;9OLw3@ zkwET*`Oog@WWQ9KxcwpbzhRS6mFeSy=m(VUNDf3lpwLWm;Y@aEzMWhM{osVVpa(e;|>XXx`CJ2g32p1vM&NR#~&$yrZpvh50a+$BwxSEsMf-_c}| znVj#rCTs4K8-CCvtpR=A>6a!sW68r%4+g&?k3~JW=TBdcM?L7VoIDQo;N3U!Xw-w5 z!Sweq)Pv%i$^B3d;L_L^3m|Pb1pye@gA=HDD zCF$pEs0SHFlYh)XKlv>A<5HYo%G1~H*CMYOPyV{qDp427x%MOPaiy;-A4mQ*iQN61 zRe~;&S6xT{!Hs^-gnIB~B6%wEqE{ElB~cHCSD>G7pdJiX_Tz_853Vcw@l@1=ggFW#1rR14|tg`MSx%DWktZhSI7oKR9MyttBXTh&O$VX$WQmH+Cow(R4y%WeA zR$Ha$Z?Yu9@4@u-`y{LUSV#6bY84k{KfW*3Dv8Q|yz4cqtXWCFyFRo^=1=4bs0U~L z>FZ$BgP#k?dot@X?-@BTpDuI%6V=81|{pr6ywKTaJ<9@!T6 z8-(1Vvo1$+)6c<}C(YQ2Y(!$dbv60xR9))5CC{}Zk7`IicU^}2WHZSj8_};mPaX^Z z7cWjfAOAZq+K1c#`~P>F$X6celKzPti}}!2P3Y(FjE3}_MGno6`%|gpR+tC%yT`Aor4z0=eI%;CiLqGrOqe*ao^37qIWIIGQM{D@|ps#bz&}4oDd58^s zioAWHCaGoV>w7CT@tr{aje79nGCB7iP0m%PuX7*P=N91=n9~e#O>vWtC z@+~4KKGS602l6E3MHK?*Yk%ZLn-a*GvSOYm6MenM5qU{Ba<5`mnYf3Xv7A*penQ4EyF5>lzcp70qRX9`Z<1uRj!UB7u$sMhBB{mVmIcKTi?tr|=iGDuPSC=k*$mfs;{oG8>1^@efB}c;l z{>|v;T$mTlpp#Ey|3CUNc^mq{G0yb!9?XY!89=^vU6<`!$T^!rN6lYd4Dyfzxe=p|J-Z*{9fe!tDNZXMneqo z3nLGWHe^aHxhwoX?L65T{+BHDX&UzbPs);KY%=7tvL6rIYlwR{uD4FX{pNAx_UWiM z%KSv+~MS^uDCzDkh~aq zP}os&L)`yb`kL&D^T4Fy^lyKh2hRGF-}OPBG>W_*=YgC`|8v_!+$TQE^&^^1hJGVI zTVTU|D*C(p8uXjnlTUhR;!vHwe%(-$ekSsHe@%|uB~NIliC;bX`cD^4{=|^G_SfXw zbMn#%O&a;o*T=_ca&am77V5#i&*TTF2lthE@~fx^T{hC+hfxnQXQiLxQ4b#WBoD*+ z;Msn%C(Z}ii_+ImPH9qWB-xg(NmwemBKim6uJrZRCz^DgPIiB*iR&G5>^Due)T6II zWVFiDIb`45Ryp?{Subprw*K_>IcKZ%j3+y~TIKZ@a;PWzw}JHass>gWv5xHEYn2dX zKOTwvC|=3OvmigJzkxna?uq%-Kjfw1m^TfkzqdwOCC@tYuSr(<`i(pb`H^cI`no;l zjg~DV_l~#9K;=H?{tZ^qn$Xu>Q4ccO$jy&f<-#rUxHRNnmFee?=pSs4B5!|+dD7$L zT_3G7ssR1`DT6NSx|2sa>N0!-*}W9@+n>lKD&u_xFZ$W1zAh)MB*{j_Dv9HkOTnv4k4fn<7Df5Q79Vk!B(9gxs z>#}hWdBI)u7nJ_ym)CfY;|teak@vT6L_Zfo-ZN?nc}oe*#~dX;cf+L1FM|M{To$Fszt zUu@xe^2joB=l1lDV2=kg>$u~;d#IHR4yQi{ErUj8-)wRjH3FJY(HW{{_oUNTru3si^ zM}OEYBYoWt=Yf{h$PaKH@a;_goL!S;LG-m#Ax$PHkUuzU66rucFLKjlM@Mr0TADQ2 zM()={lZrX%=fy2F8Qq7xzr7}I2gt{}YSOAGeSHA+AVrxcUxRwE>@3$Op&q1EqOW~X z4=PP3|1vZ=afiGLdC?5zerFxzMcwAp-@7+xQujT%>`vtQ0rd6Q1E@Fa$-7gK_hg}; zAD+|1)|34G8qPcW$r&EtJWzyw{{BLfpa}9)lPR@dU?09lwUgRO^ z)o-xZ%IF&f&T~0C6~nhKk7DlDdvL)DEEJS2cjP^g8t5j{C}-7Z}`E4_b)Qh$B!3a zo~aG_>N?Ep%p)i5L4W%eITCq~xg7oc5qVF$LF5WgbUDA7{PrXAr?=$hnQ(u%27O!- zdC$}ljUzL8%g8hF)3vzH*LkjgFJ0SnrFp=yw z2J`j{$@{0{JqhK0P;dCZ;{~qA!v9s@l0TyV+a*7J8@1bzvPxd=aTxb~yKp_=4F0^+ z$u*D%ecnW_f7cM(Rq~qWxIdPW{$2DD_ZO;>3*kOkYG?BITsFyWB3E)kKVSzr0{wxu zcgY()(BI8PUne)h{o6X^*)5R2_aaB!w@Bzn`g`hs7O|fqxB6m{oMq_i+ZnKbjw1Wy zMt*&PY%Zcn26y^;vkUw=ot*5hiN|g79@GP^E`1$`deCPMc_iw=@>k^As0XcD($^1h zK4`LrJh7i9(=*V|8HQ`}r5kxj6nwjnoIF*NhDGV?cbX>mMv-&HYBE@vKg_mFlk_U| z_475F{4$fbBqA?*MDC9KsIC`%{SE!)%1g+zPHM9F6FKJvO$I6ZaoctLd6j&;DC$G5 z4D|I5)Q5e_{P+OWhc#Qc-VF8ONEZ6J4eCQq<$V=9>chw_TrcW`d_5!m{Jae6M<}^K z6{}2JN1lT`Y3x^Wq&NBlE$HjrZO{*1NbcGZ_YaD^e1Ooi~AGX$i;VBCEqXdIGhJwH>aOVT(F954*Abr*uil@1B*uPQkqB!XWaC#_;oEvS$G9r#&LacE|isHTpTr zaNJK&`kP}WBhNX^_4CMkx+w3d{6@d`nR0*g=@!%{B~KY~2ziw9o{A|Azt4p}&U+j0 zNhtFv-!Sk0eL2^?k@tMQO@5mh@3G~kKTE^^ja|v#(eEwlLl#%Of1~6n$B_SY8^iT5 zZ$oC;$$J9uzK(KV{~`Qe`zY7DV*kHKnO`Y|{HK|cuZ*0D`(GvLUwe!prM<|b7h}J# zn=cu1p&!{D`B3s~a-shW znSX$s`zzi@dQQHZ$tJc!^mPwMn@nj+jwp`%UxUaE_F3fnT>9JfxJANWkkijuB*>e- zZh6BZ2N#kTJhDj2Tk_G@7Kv_2Utj)ekrOM)X&E&+{gb={^`KiAeLV~Hz`C96gL>eV zm%e_E^FgHYepXaXO$wjj`m09R@0O>poBJalno6D#sL9?tx;hY*=+sJuvUP#DGKX=wO$vT*gExx)1;I6o|`PG5gd*5ss#+~Bw-W$%*Lq-rv)Hht}TMU!_Lxjg#M!|sxUpK0>E z2K{{PohGAZlH~{Tv@7HSS+W1EL_c51hkn9Xa;B2l-=~n*SHe6{A^Q0U&I`+XkpsPO zA1#sWkNjxpckEB80MAn!>lPrtT)sY`1w^3l(_gedQ;6hXhY zv@&nF4f#(|Wqv$QalF^Cl0Ux@^PbC;`NMs6aKAR4>pT4M{^u+5tWew^%R&Ec?uB_! zSMuNdzh@iry-4IgqsYxBz_&}t6L20FrM$1z-D=2#H(YNOYsjtQ^lwX?4=VeUE3LXI__MflGZ%6vsBTkd;{mG$g(GN%_ zAIEv(Wiqn{_j|FfdVcUNJa=LmTU_W!Mv`#;@}Tjf(2eLd!iRnk_H zk0AfQ{fgY^n^hd@)6Z9O;(g2+qly# zJ~+|0&Al+sue^6#pcUqOI&PP#;?SCQm{>wAGdV zZ88de4JPNBg!^8T$YU&q^jc3YH4ne<5_$D<)QjKb-iq8|y+td^@(JQWS&Rvo??;TZ^?IXKG@ZkzHWBiBJLZ>Gap!_eh&J2^-GIfA3$FG(IQV$ z$i049B&;HR{VfanrOJCNUGr*k`7zfuoEJ9x(brSaPY74;1Nopnyv$5LpRbDZL0@vU z+VJraa$lSuup#jJzwcW{`@rAI{P@^b*pJ`fdN}Gui8}OkFPtY<#E=K}!1@dFoc`!% zH>0m#4@Es%M7EAb|Mv}f;zUhe`Ow!FW}-eR^M(siFWM;YW8FbNqGSX5dCLm)Gc4pH zI8U^{MlOW@+>!sZokCx~T#x$-XULCme;{X3`Z<3Z@|)h|KKHD0Ad!6TJ^Cl#$zEA> zS>R7U#}~!CIv59<(m~JFtZz2bB5pA;HLttX$8~19|>#^56G?HrylU8-sbEob>a} zDR_^l7CBAB`J@|p%UrzQHj~`c--7qdxj&BxvPh;@^!1p|7FiNU?%CHOT|bg*4MRO> zO<#XPJs7i!ycYFfhmx zYmwS+^!3vt7CB-juRdduyU)mOmn<@>C4D{imPLAPAkV~k!IX>s4n{vghLG=ku*kDi za+Mzzxmk_AZiGB(sGVFOhbH6RkdquWxg9`X*F(LC+)Q3n5`N1?Kff-k$*}=sA2;Mr zN6G7PzKAYGKYP}}{(lVl4eG}`W!~@;>PPPK^mnJ$nj}S$-GVV6dWJl%3-Y5<^mF&V zsDH!AcTg{K9Uy0q(&TJT?ynzC)1-eGITH2Z=SuRf#hM&=MNU|UJghGL+-w*6-O9XS z^ifURPjlV(q9#uY(BHB5QNOy7dm=9?wwk;I=YezoksB1S$`0i|er@>wxpE)hR15FL zout3r{H+pGfPOB7en2VZ{!a_c3%Mn5{Sf-WJCygv$J#Nk)0BShi@c|jPM)?4_q(o< zd7rj5^8K@ixZVT#PE93GDO(ut-zf8i11sVFlk$E_X7qF4D*epk z=;tOn(U;>oVSnF%{1y35E5+Ys@&4P-v0Q&KMHkn(pStj^JlVw{HKwUr(8nblf{Mp^~8Sus4`DJ8~gd~-MJot_k8xvAdkt8 z`%g;0b0qo!D=%`re-^yAqs*_|&x!qbY5G?yh`hWxxvVqhB_qf?-LU^&OrBpG^YqGm zc_{h?2VQY~Q5)2YqV#joFwB>Fk>~V8{TN8TJP`eh7;+XDi=-B&pRc=Fq^NRVe|}91 z-fQ5xS7Yq2mG=w#A+Qfxmt5%nUaBz>KX^F-1_@{(s3Nmt%G z3`5@Zun~Rjfb+%C#pHx<7U}qf-2IP5Cbg%(>*KspRr&siw#cLEW~QIF+g2``DKUS_H=c|u-(Rbwc%{1BFgno{1tBJoxE{;4Z=@Pkb zKb!}g>E|&ck!LCW&DK-UZ{Nc8bku{4U&xt|C#`BiKcCr*`Jh?kQU`I~N+rLN*#FN`?)y}mf&0A~=+mNxE?boSxAS7u8>PQlWh>qn_2B2nUDBnxGJp8< zJ?@8%<$4VAok7dVX6*knoh19<{_l$yR}`L`pzBv z`_CnR!Fi;Ja({Ck@+yah^z+&G$iI~D|4I6edUccQ;hAv0aHXF=7M2vHy=QLO*wTjC#|9eB>+o`5VX+aNl6*C-QU` z?9aXF>#NxBZ?uskaGxw(nKumVj`tYc>E~X^_cMs$mrjwIiM@82r^&R9p6{7T;4$_M^e z_S**LcdjY@%7bBe?=2I3`2+JjKMIl8VxB)%c^@l!5bonD?~!jAhW8qjJbcMmyw~i@ zpFcep-%mM^?7K#n8FR_~*Wx~~kpILZPkKfki2DE;ve2(9&@Vn!njEqd^Qb=LG02O$ z4kGVO!#rvXdGSqMF6|-TLp|xP8d`9~EM+J)vRqh8ptcLyiA+E>Qh5uaW>vCQeIXjgc-ohelACddw zykKiVU+3+J{q1IQ_a4Y&3eeYW`=kCuk#FMsFzPNj+bH;7y$=*+k)@mH?}n3*7Zsqd z1EwRti6R#@TcnS&A782CeBe!A|3*F8w~E{t{fU#A>E}_4EaKLeymOgFwj3sxz?P`M04yl%=14??yg0njC|EcQs``G5@37_bG10yy-mpd&CONmp&y=+JpPR%J)O$PQ~x{ReXIP_X!4*6F(tunoX{Z z`?}Y6k#iKod;eF-pDO9{K*>iEYGHp|fWCa$1oQApzP%g!{o@U}?u&f?MN4v(AiTd5 zLQW3G_mlM?=j@F9eJJ?{=KC92$`Df5Q=@%aIH=-Y4je3;_zA(M6V z$GJWh^Pzqr6 zcHJz!a?sc556sf9H~GO!votwCe)!QW=ZewSr+=9xTNK%z75nw8>&e~GPlzPD_QrhYHF9DfL#lfHJCi1SWQa<-k=2=aHFM|$ichoU~@|3*H2(<03q z(a)_=586#8???Zl?-6pAJeY@7<_!xWKT1^Y|D5y0`N2kCJNs!eH=X>pBl^R|=GZ`Uh~d7>00AFqLa@H6E;{++{`w6)UjkFILsaGmV+T9ZeW=;u$F@m|DO zvIBgdN7;{eL_ar23HteEE8LG7N?wCJC--4;pK*BaAt!zNeYREJ)+9GVzcxUb2Y+%1 z^ZK*6zUwmP{Wp_4y|9Xlvi~0U3-`gka(zc$F$oIFB^CJH*afw_5>mJJeAFpQceNFmx zEzS?Ib;&LI!2fl~qmdWQ>q~Z-hWlV9a$6hns!im{i_qUqCqGz)`v*VCQ_(Npq}<;R z-Tn9Phkgz|i1$dxk^7*3(0(O(;2G4vljQa$v$##8&&S%$a^M{K>H@P|FHK($UT&6- zBgt;-&Ej%`oO!!hGCR}PulJed`xx?_qh`5tfqeL^8TYB^YwKmRbTY_Q@0jJxEArwe zX899He;;{mmK!_Bi$0sBRWbUy&~LNUnnHHVVv)#aYuf48H*zvZ_`gJkk1)Qiwc z^z}WQC(`GTmz1%{l;7kBIA6ROM1LQ2!~Ruy@33iAi`-~JUq^UYWRdbd)_|H8IiHt) zu3rc9n~~%n^(@lrD)}h#tL*jY>-9}6QeFAJh!k&}m)>#R!_OkWTGHRyTUw-FJb5VU z(-9>PpNaZZ(wqK{!uiH#kcZ%W)A%acG644vlsxdTrw^z+K`7Ad!myu)OX5z6;P z6vlbvxjX%wvPPv z?B^qu`pYp{oldlh3F5=wv(sB_v?0( zYfr?VU&&Wy&cuADGQSdz{y-z;zUTQ^e7}ow-*fE}%$Hi}>vZHn<93ppqJKQ&CV6L1 zvkXzb&vV#7vp77W&nHHhrE^XC`eKw>{AQB7Pc_RDB`+_mnZ?(GzJ4^vEcq?ueM`)e z>k)b8YO^?c)7N!yUMRJkd~g?hoPoa1j`Ko;zT_M!W~qLPyge0veh>Qk<0Z4iE+ikj zX_jQ=`+}Tse&{}o{(kn4}qb?$7)OE!~B=dehXg7kA8)T39U$^G+KWaVXY z3eGcbHR$VM1u&nayE`N;!%@+4lRZHqukdVR~G&9RIYDz#h<4z{p^PN z)VDJ^M}3@sl=;FTKA4BP!}X9rxFMxq`&e$*6_C$8QYutZ@;YYu`;hR4f98h-*CNn zP22~oM?Ysq-hWE@K2D!rnv_@O$L~a;KV6Exz8i!6yz;(byY=||L9RDGg72HjM?bH^ zJm;Cd1WrTc&}fX2Y-wG{6{5E={67h zf2Ds}X^U01wE(Anz@T`%=pLv+r@g>YDO> zKYs9g#|!+tLw(GrEAxXn8zHYy?)Ri(|6W&_4{s6(-xZ)=hk!e}k?%CeyrwVt4Dz4( z?a1Yj|0H)I-^YBv+eq>x)QdR_$ex)0v>Yb?g#RbJBBvwY2`oLJ<9ZW`tVm%yK!6ydCxIQhEA$6#6Ac$B=J5!+-B% zk3&OIAAQL=2BM!ajeH9JZ-1CPItKHv8R+9@YcQ|WjO>tXkzC4teCb8} z^G|Vo$#a~y9O>^j4w(P%NH(01&#WX*@X#dYH94h)CR>&Fv7+GnRXY8>9`6(R-68jh z)nuk8eZ4ag{oq;TWv6hz?iRVoW4v#wynj^y^O+{)`+cHHVV->ieO|Qz@|iet=D+W^ zohF|_Klbw{a+?^;Pb%|)bvI#trWJi#9e!^*lDzpI_W#QFvpxM{m4w|~pOXvwB|n-zh5avp5FrgSN#1P{;uG~^;F3=p$Ya`* z6J5>tzDV*953{UO-bb8OAO3AaU)ODBmP%{LUs{>v=O6OsV6!ytMt?W&YL=JDpdA^utmZ3tQx3ZY!LN)rjDDtIt@npYPvz*LN zU#Blb-ZPVYXes`Fl=sMUFE`7_2>ScPO4N_1S1U4+AR_13FZB(=q)(!DE-Vl+fX0E z>1&4_`0uTJzt2pZkIE|ddv2lr&1y$KcR>AHI*;swJni0P@_(pz*_7{xn27UI;Q;#i z=UubhH<4SRUNksNZuHYEYrc~==EnO%wG=-)TcpKs@)S?_eFJ$r@};Fu$*<7ASX_aA z&OQS9Pe1Z?}&-oBYU;n`TW`#B6)TMauS-Jlcx)=G68~tlsL0%I|?(-4- z+-c;4n7^;Eg**ZM*yT6LQ9+o;$wq%p8g7*vHOL`0t8D8^-VDD_R_42lpT_&2>$(2z zu@&!Wlk@z+dp`HcH}WHIQTmllUGV)N3ZHSs_mls}&;R{Ct@TPC-laO;((wik^`H`zS+(4#DRW3o6jt99mrErAC4xHTa-1+@Er7Ymde=w z^(U{XX_o#e?AMuC$4(!~8 z{(e3M`Or3US=7JJ8R_Fw)6McVfLwVN=A9IOuSb0xsO-NJZ8$F#qJLM+!{3+U@1t>M z*{6IzL>|viGFr`O3( z;L{{!UVIetmU7DdpJLd*mio!h@56jd`O5Td3GAmQw7c z>{c52Nh*0K&KL5HoTmuB-^7*v%~%NYEgi{S3WAmW_%ra!JzS40jQ7|clP@Bl8dZS) zep(#&sTz|r9yUq2RQi1N36pFpN?(6TGl^$E^5rWg3ExF_xoeU}+3D-OPfaqT8+i)q zL#tinZr@C@+>yTa$cQ{cx&Qejr&)TQ<$9lj$lt5e*D)o`@*tMnv7A{R{URs0BmWMk zzXy7n_p6cnfqN_O(@q7Ol=sEszyloV=PTer6UbHCnWge0 zay)pbKYd-kJ^Iti{md=k#BW@0je4lHqQ5JIm}P|G?>V8UKR39Zvm^HF#pvULT`@1E zPJ!;`ncU3)WcBnu9aYoJb5QRFXW`ZdDvI- zfLmtV|EHhFyf=%(U~(e*5f_yG_#yPS)%OlNG{F6e3iP>WF#4+l$>rdCXQkiiG86sm z&s=X9hyFq{`g_41%wNTj`(D6%70UNDV3LAF1Rkqwt=?eI;L+jQ2r1DEZ1z+_ze%|g-)0KQ> zdJD|sjNtlS^jlNrk(a~o$CP=&#I5+g^Rrwpe*y2|-y~OgW0kCT$WeID#$Q=4SO)Le z-{ATb+UuKd^CF$$RYfRECoLqXdNxb%vO}k98Avb;Peb6LP zy~x#3ACi+4K5r6x0s7kQx=CDylQTX*JxL?4d|{Gap7izW_a>ROfSmT-BrP-0*Huv; zQiqYh;`~tl9=VvK`Tw~uG_?@+^J(;XC7dS?x1!H8mcr+i_sI*om}RDtm-j4dmc{4k z^Dn3`2SVxZ*X7Of;Rd-^1>_|?=*9Y|$bG?G z-Rb9R;MdALN{5PODf*J@N#FzC^z$q5UFE&P5|xn0+~oQO@J6M7S5pKK?uo@N$*B~LMM z9(YxfzOIP<|HfeQ05QvED>)bX6A{PBYu@5~mVrJ#2;UFzBL74`I6&!lR>1zg`w6aB zN1m@c(%1EdS!A`6uZ%*T(?5ahWmn+-@mKQH{pjzuVxE$Md7XE0qzve_BufzIT0`7nR3lNVR zmFeUy&#ltNf&OiaenGD4=Q2IRmBmO?FlMj9|NxPo(_YIsk^1dKf z|6!80qv`MK;5!cV^*8WjD>(^#LYY@-3GTCy{=SKNmDP#9jtB3WP2LUOk(qvW$Is7? zAio4xx=h~p(9`F^%Dxy@3iEPdR}3Hel8@=Da7{?TN|S~yRvBX{&MOYZCBb2v}rR`Qe&{mfFu zi+-+ve9K`pd2kH+!^(Zm0SRWA_<-wU52HU{l79Yn!z}NV_iKlJG|SG#T#vy#QO3LE zPRMf(RivN4V19JX2y${5-V0UUTlp4Yk!gZZL?If4zjrW|D`IVM4HR+X;K5Vug?`h_0li#9@-yrW9)|m^3k_99gt7Bk`3&i<64q+^kndJN+@@2b8oKBKIE`-nA>F@n3Op+~zoM!{_ zo$uuH+u+AR^!FRohce3d!v!3|c_NJdzIn_f!>*8xQzmKJnf~sOYLcT5$hC34=sleN z_C!7L{X_1QZW52#^!Z8fgL3ruO7Lc--&qX&+?W3DgWq@c9Qp7ClY9!Hzazo7FOt`R zPqv`1qw(|2c9SpQ_qQ%VfA_g;lCn|cwaCwYoh9eGYLXnq>Er2GU)O`|1nwJ8z5_0Q zot%uH|E%2S(QcaLh8z8x_6YqMW&eE`d0d6DT=)8I68GigXSwhm#aZ%2oS)YHB=2)I z%j0VF@wECl4=MM59=AdNC6?1^Rjj`n5a9k~5FC$h|YJZXomP;3T#fn4<;pyH zd*n4!dU3rC<|(6<{-@Vd)QeSIKb{MD`ayF0>X>I&=EJ?);rp_b`NMM&$kX5R^Z5_+Fe7l33(ZwpGbI{ip`&gynKytLdRk|D| z=LxZjPZ|2!c_Q}T2Kn&}?4RG0Ge=tGq#u1gYA*bGom>?0VM%BDJAH*!3SA)|T7&&< zPx^a3&J(AUc`WlbVLv{R{*K&im5Ivyl_gtnp4dpAKLEFKqt7kiBKOEi`1;ag>F+~` zPYtrs*U#|tIw<>b9ek)DeVq>WizZja_e)jYfAB+oY?L2;U1hsf&h8^0O0>%GJoNFO zomQ#To$LvATTB*w{_QpLQsmpJ+S9M!!S5QAgOjatL&<;pA%5;w{2hS!IC=#?pL)+K z*OmF=uU>%Pa6S67Rl*$TUkBvZq!amKe$>Y;Qn@1nRkGPLdsWhf4`N`nhxp@-^P%1&`3*xtv__2kvvcBnK8i|4n21 z`l?fee2*sYMt^3RbaEX(%mXb+|9VZqe782_S4%N}Xf*lLz6klag52aT^7&`Ti?N^I z_>r7i7WZFE(Vr(<<2^-XK1$VLk@8lle;$|>DP}Z3Ump3I1DnZT(O*;M6nXpiNSUhS z$9I-MzFF~iYutxAXG4Fk>>4G_?aA+kW4=K-^6W6|uj`ZlM&bSR_T+Oh@MQq`pKljL zkc(nHU>o^jbd*G=k~R3Uq%HkBVKUy^twY{88sEP^`QR|jN1H<~F&O6&A^*a>hQb-- zpx$_IUFlC5-V^z&F7$I6=680nkdv^UIgxw=U*9{OeBy~IbIx*qe(Q}XKa0}WJHD9G z+@D-L6Y>E`1z^%Qv@&JTJQ z`uyH>tMq(K{;XLgNO?af3(glui__mltl01FBYT1?^`pO!npJ#*$unOgAGeBp`WO1Ml=&)+?Qvf~*^eJV{d|qm zf4IpV{aAzO>v{cfA9*>sGWwTF+#x@jiTX=v`nked^t1Yt_n=>=-461&^WoC_FS)`) z^yeu3h@HNpUuiY{-3Ir+DtsY3)WH2QC0~Br3-fZ6e0hHK-+!J$KVL?^K4K&J(nj20 zyGp)&7WZ%ekauD~U$HX%Td@G~o`xSnE zl|;&ySL9&SUp;N;&r82C|GxE9f+qvVGhc^>|L?Z=Xbfk&+)8#v#LzDB-){?g&LhUAiGRsH&rO-tiT*D2(G;82s;TW(UdNLvdbtNv<;z`L!(cZO~NY(>%$Y!x4X`knb|wj zcabx29{F>d9Fk;};z~Zn1N-|?PW1DfLsm&q-v5~e|8I`qdM5aP$^r5c4^r-XemIT(r#)Q1gM7|E8mlkLA#Z{Ms4w60GEw z=)YXElYD+G=I7iecUl-BUODO4@VyZdT9@1u_h&OH`SOa`&xcRr`qwh>^%`<7%tP}z zMt+HT`OZo`vkvm}6O{hJgteGgUzq+phxrG#HORkGBPF{#dDmOaua&zvmE`~p?H+|EBWvTh4AwxalJ@>%uC)# zjz)f|&`olD_9$7Ljs8uxg^wGM@7hGkB!6;6@a%=;IpBnoF=`+=$BXS`)8?y zdR-CvyIBqRVUy`6;~{BmOA$%~s&m-RSEt`240L1dV1%6kvay{zK=i0c;+&%P_~_Y~-F6-#aU_im6?Li&+CkPjOb zO)iP^NUp8qw{vj*xj~kt$RF9#xB1pu#ikbdNSsynDEX8T$XC@E&Gib%|2&E%x5j=y z`+4%IBUX9ymHh2G{`^Y+VZ~R-kM*LT4die7&LvkU6)xSB{*^g((chVc{@#!Md_fQL z$gY?_r}QIEMZP}eF4tR(36~tz= z(a)sh%hPcmcbt+hk1qp1ALi#5yJCK{(qBGw0Olts{pFV|nCH-fK3=>I_1>Z6u_w`Q zsmxm`@jOCuD*fdyIsgBDJj4m_4=VYTZnz&iL+v*i5-G*1(Wl*}|LeCTPhE`rw0`8U z1l+%xME-+*%dmyy)8`}Q#4hqRv29<>q!5ON{^I2p5(z;&pn1* z|9+&jizREAzaZzyHxPdc|0K`D&(Bho{{4L|QhN3zyMn)jlShD$>?P~h!OzJv&~KEt zIQ`xGHu{^q$XmOa@=&=S`nta<&=IE8?Ln?Q!IX<}c$dp|777 zvr6k_%q0YCJ2@Ro-Q&wWX zA5LC@e9jj#J=PD-Q8C4%t}9xIc}B54ao;Cq28n1@2vI6Dn*t3_}M>JDUim` zH!6gAK_%$t$Jo!W?o0ksGh7^%e#Fy_F`x7k*SCA1U%nN6ea;v6ndg$52ckavkemzm zzhabp`SwNV=k%bTV-v%r)fjTGik^Ge8rz-^C{lkj<`4075!y!T$3{&fX6=s+$6wwXdUz~wiSH{tum-X>c% z;Qd`Y`Z|Anr0i=&GKaDQY3_o??nu3OOkP{w+aYUs!0$u3_YRYq5TqymLFge{uS|=s{E7 zOd?-7jq}1I@-CbQ&UT~E7u_>u@ny3AQ&aqU(C4`j530W+PyA%c@kshS{yXAFL;8H_ zZ&PYzklSZM|K)D_{WRjo2aP^=utEO9mp&hc_0%f#`IEn_@@q;yyihRq^KSHS2kh@l4kSOCf&QY& zWXDMC-&c|s!vA*H$>mm9WxOr@>%0l^umZWmF7)U6l6PT$AEV4$xp@Hfgoj*zosRt3 zFS6}d)MM(=&y6x;e!&Ft`<&>9KT7_Ae0^|U`ne+RGo0#8UWWet@f*pr;eVg+WGn8U z)vZh)M-9dOSReA`8R3$5I{D%XtSj?v3-3pN=xwfdxrOsZHu~1(+rNAfISTok^L@xU z8e$%_K_20Y_aPI=qtSme z61;lzn&Zxf&Ke$ za?VAjc=n>dAFM@v_b~a~HdC_Kq`%|$nR0zSIpDY{9@XgY&o~d20eF=Tu>&rtI7-#9NG zCU41u`>i?X-*08ZWllYEYxHA|>`&g+1^qh8d&_%+a36O+*WX2;zeedt%)TL9<`$uE z_aDc+qE_TH55wg}5cy-K2syNf{HrwP9|$=X_o<3JBd_R*e)`;&|NZ_G;CpLrvK9T7 zK0f61P56GH&*z6w1?dNCHnJI$sRv(zxN0E6!JZ{9O&D7g(Ia(D{{PJ zq&yr-9#J_`^2{a&*NBv5d&s%)Udr%va$J4P*U3%a&TaH>z9zXM;>)(dCIjH$I=9N@wJjT{4w}z1yi@-X<`y(M{YC)7KX`6|77{OkAOdQI&2yDIYzJ&@nos@&%s7KV5dOkZyR zw^i=v4@Ught<0~j(ir`@vGj5FXvC*L@~0jD@-yV~2l4maLY{pb@$o%*!&&5Sn$y=C zk>5$2MXri`e#z_PYLBfVh3I1w=YfGO$+7T#gOTJT6*2!v=|>E15iZ4){7SoC*nj`v z=WQ`x;iJ;O(sLgAJ-gD!Hptgcn@&D)E?n%C{=tdMZ(SS>s+;0Ek^HwQ`eROzTXn$xv>1Ipua7D7#*+6AGbPI{^2+hp z-+IvB?KM-H93h{b3BP&K=M@&Ap74k~IM$TWbLscF>rIhv^m#Jk!9Zm`?&}?R}> zB~whL-@DIaoDYOu{d&7xR+*#B55ES#|FofBmsUhQ zxIB3n>dW1lk=r02U%VT6I`;2YCBKpkpC3{9CgQ`={rtSq4E;pPef*c$-}h-uUx)hR zK9my>RUq!mcA}4cm!e+Vmb?=Cd(Ua)`fG9DSGn)G zay{;6WTKz{#9L*DH~G~zId;ewdAaF5&;^A~pIrk_W(!91hJWS0S$Po%uJIcQ3_Oq|d4Ey&kzI845sh{5mH2XdLD#Y$31NhJ4XAvJ?7y9{(cGJB<0AmFV9h zs0UB(O#X5X@nkZ&|7EYYoh6C09>0q>_?9 z&ws>_q_On%X6*m#Z726l$NAtTx%x*#idUwuw`Mh^!%%XqVx~MwApdm2{`?ntZevsG zb)mn-3;E{lu)fNK?9Jr@!w^G$qdxvWtfNMN9fTAqx9_Tl(B9 z#+0OVa;2rF1S{{y4OoT#9p$~|LhDTl|H=J&|IMb9ElQvF+m8HBRrO6;x+(qMj_~_2Wj=VdyQX|n-e0@;1pQ9Rdz!u8nUbO8!}op1 zc|oaX*2{|fNoDBUlgQU>Z%FP5zi;S5PKDpomHw3B4Un%;`V9*@p}#PJpD+Ixe_uQL z))x7n2TDEjEb{l?{kiVZ3;S>7{?8Qbzq@bY`iddAA9j^&LH^#Nyx%GCeWH?Yi8Rpv z(v$vfZAE={A-N**JK7bpb0qRD%KOaW*xw%*HqG_Fd{g;3$gc>w>jK>8vZH_V#ey4< zOQOD8TA7a+bqM`D^SB;;1^q}z$&;R2WnKn3^)Jp3W$545#V~KD6L|#Y>%5yn&e9I= zy=)?%9vm){uahld;ZiCS{rVBUU+F}Cyggj9bs<+cgZjZl@>KL+euyO>{et<>C&;IB zVSe!|aW^Sshu*lR872_Lu~i2PC>NBXx5;z^@-wBC@%%Dly;dkq=$it8WG8PcjO{p^=! z$d_K^luw3yUPi8!!<3tk$r;5>S+DfJd~!0S$V~b?pphwyUXo*5o6=sne{{tg`HG|T zyJKJMS3T+TxkF8{e@1pizUSv!`n~-m&)RY+;>GM4&ObJuwpZq*$%GZTl zFMSo~i*@AqRQ!2PlRMo$fQVDOoGwzLByYuhJ1dn?^rp@<#o61G!@te7*8MQ1Nc) zSO3O!EBxQAB7K}N5cLD4|1iwoDvOo(Jck6}zVRl0{^l^`r(TgOjKJToA$@%_2=@;n z$!Bn0n4sMEJQ<37kn%or4djFJx22DpAU+JAKn`1i`*It|JvSphDDN@59Y%ar>YWF# z;eJ{J`gqVQoCgMz@1Z~Q=RAc=;JpGNm#v3+%gX$e;qbk+D1B=|Kjzw&gMBc-0%(gX8$l=-vAI^+JWKYeVXJ}`6+IRX8i+xL-s1V+f* zC*-&gypK|V{@pM>LW(yd2jP5iCXoC`$NjNIRj}bej;~Yg!4>!`g~S) zLqfmM=c`H?(zy+NJ+PJ`P3Drjx*9U=Ci!P4Lp+Pq$C>;LiS;E{4>m-bO%4k)B*$^` zPvm=cD));D!S~0?)6WME7&26;k9uD<#BCnen>OO)8}`5k#BEEpVt`xe-@z6-G?Lpq09ql z6^i^)b^1JTnki|$>GPiv@No?DFHIMivgQQ$>;9{;pMOZd_uFd9=IiwPoqhOu)9Bj| zM-dMa$ZbxUa()}R;yF{s?;*bee@!K)oHxby1KH(@Dc6+xX2o0h^Ax6k?>)fZU&*gz zMLx$t=|4R7&6FiG`FY;0*$30M=&78vZC)tpX`9%JA-^^N)KJJ74{2L{|5>gBE z7?gg)kcQ~Le9q7BYhjh@Md;tZ9;jC~Cnt75{@#zA0QQO@=LWYtPi_N#th@)e9_-SO z{@#Ol^KJ*ZD&oxrTlzQ&@#VBKPpd><+Hd_Vstc{67dc=l`YX%Rw{{75pWchS;1K3vhLC?Rr<{%oRCkc7v%ciet&%vYafLh!`TKIE=<8Q+3<(P$ z_x@$bscqz6IZ^La-ot2H6n<9TtE*YTlzeOH^Fg)IKj}!H+cz~O`WShF2YkPfem~g- z``rcfyHh{-JBj=CyMf4myrSQWjx%MX6Mb%@p+4H5K3_ZszE|#l4qa)=r7`q*+0CXr zu1NoG-GzRAjV$|+Z&LE%a>$e~H@RNvAkHW6$es72-!vo7oTb z)2#Gw?~#}%RA1p4@VT-d-!vcbGluIQHsU_TQSwXFckjI;dtXKUq7;4m4EdaKp5%#t zF#lo<`3Ca$#tQP4D&cbVB)N4n+_!p1w#WROK}G4?;>hnjZ%%d}7cRy9$(j}MWhVJN z@;wE2kxQa~r{W#**d)BSqRbCZJBs&Noao=6i+CT)n_TNoxHM7b&$>KC{%9@N*S-&z zE=qsl^dG1{EB7^@X2Jfy4E_5U{r)YL{uHAy=H-mz`ela*NmxTZQ!zq%TqFk|UhMuw z-sFsYP*wU{ZxSJ=dy$0dmeVu$!m)072^&?%*o+mH)tV;`f`q)03AqhU@xrGhM zJ(qmA0`fIy$;X@xsqlw<1-@^p+$S0{z>r}b>F35%4RIYs9=#ZT-$b_AhIn$9d>Q$A zk4p6QxXXq-2qB-jZ%AZ3xd7_5_3x37e>bFZA^JNgJNnssk(~>nf94Q5vWzKvI@0I+ ztD`@%IDP)O5#}Kjq0j5IhM!&N^UL12Z>97zet_@uCUU>NXavp^FX;FAlT6v0gZuM7 z;iw1yA*auW50(2tx8eJJP3h--ap)(SLyp;G%BHjA-rxZl4syV0Qzl#{cfXAJ5Xydh&J9y)DEsj<;N2ed^DS^e zCErrv5%OD#&kKDtWm{hQIK&3|_0Hsvd9dG}LHNChKIL{k)K9C@&rW%f z9~n(XVm0Hd@=qDc{m#PGt zC6nU}*Z=lYe{kMduFPABKtIrzsazj~{Lq6f7whr9;v%j`CL$gw z{e?@BA9}Cc*Q|UR?~~cnw=Iw_ny2)qWFS6N7|L~5Er5cupUMpI>->04de~e47s05emLI{ z>t}MN4agt4(a&i~sOQcmpBF>s>>*D}K|S{ddFdVGYrc_dzC=8zM1R-+ihjSNAM!P` zC%FFi|F2h&DbC9L^3fCF>n_}{M{B4jz98F2q2E)vA9QOr`YDxs`@Ol?zdxkUXV1fW zUi!8{j4A%ge*84zM}DPW<>XRR^hx}@>ss7DUQW)r8TtCd#5$s@oa%KWUF`1)_1>Er!4U%Z=2?)Dk$o5{t0 zqrR@(2g;Hi^B;23zaJ0}+PaXpm%)5HrC;SQ@=qs~{_zE_RxwWU^L4u5e(opo{K4>X zIr?|Tc;w$Zkp=nuMw7`ER^oiJo*cFVpTA6wKaTv+53+F+_3jGv?Uxr;dEK6zj((jg zq2va6FyDR!`H2JOCm$z2uYvx}*W`xS|L;-eUEM^!zikux`DD*Bh~ z^#3dUg%1!9&fVeq#3?w%elJ)+5L$LZ3@9O3$(ADn0E5_p+hV~sAMN`Cxv zk}lr;=;sfob$K02o_|M|HuuTH-|AAfIDOmR#*jaq$ft{;UvdUHH}>z|Tgdl1V84Eg zoHoP||E%=woe72vs7{^`Wyq3Ha-rphB=04UTy4nS3*^RI4LO-kjzPXY;VU^E_1St# zKg`=y^k*vlj&q(O-Ylih=c7Nzu@ik>Hj^oLTXDahnA?;m%6k}L#c^M20ev1*3H^(g zLjU*gkzUV~w>I>7rIx1rDovk%?|}S63GUBN^fYDc8}j!-xF5EVK6e>``%=Z}>&c_x zKV|-5uJNdUp5yxcN$6kCLH~XY!+lhxfBdy(N@ahpKaViwo{}G5g!r&4k?YqNF>wE zx`h8GchHe9$k1i&B>H(-F67e{ zfA_3n$ZBO?%Y~N6-`LQv>rj84*_<5K%aGyY$*Yi$KfRTleWD?MualiokF`D_&q6)6 z-*58b6^58a>F2nOhOB8uezez++Jnd^&KUA@3pwn%As3YSjaAYOS(k@Cw?{sIsU!F2 zuE_5Mx1-OuXEmi*7=7Nx9=_Ye{LAIim`|p>FZa7D=Ib1w-v>9qeqDL*Gu{hMuV@+8BQ9UA#GIQv2JL=*WhrC-G<0za=2{d{yL`iUl! z%fy({OX(L+UxfUU4gKu63}4@xe0DkRk172rfh+L)yTbJi;5&Kg+e(N}D;kjfV^MEV z?sJ}AZHikI*O$cM^Lxm4n{l3cOssE4NORHqLO&;un{(4)@^*?{- zLw~5$Brh6=dHOxb53HzfDEX9|i!o1C$)^Nw!u;ne{QLssd-DGvd&2+sl=0O<^u02$jJEO^;r|9xR$)_}p(WU!0u50nSyl+5Xd&BP|Mw6SO z{yHk2Jo%?C*&mQU6vzFe0`zMT_Uq+6$SYk8>0*+XH#OwR7sP<8VEzF!4r)MJ%?Mwh0TkNlqN!`*RSQ09-%&3 z>Ye)#A7(21@ddvB`t9lSlf5u+!Ja;T*4LE3t;m1wmr-brtc&jhyol?l*^!>wQH1W;yvR z`jG~oCT}c+eEmoAri!>vP@2Boj{61eJjrn$nCCW%9NYu(T!aTG_^z)7KxR0gulOJsb*`S(<^wGsMPlzu}`7gI9%($9Mv!@rBk?+}kRDf4RAv_O8p9Q{0? z73MVuk}tcPGF|B}ywL{p9DZ^=$`kX2%G1X#ZIORc`o$X{zTHvY2fEx5>&yB1E?wZy zQ{?x)s9%31U+sfrSUx!1`hXW_na7K`hD zJ{q^ol-ae(TQ?y-_9d^_i+X(|+2OP)J$93GrJ|oVo$U79l&pE^+kf+c<}@Ui$%%f; z0pz*GvEPm&KXpRAS?NzH&=~dRyIjwNd{2X{^zo*ic(1GmIa?t1``yTSC*Zxxsbn#6 ze{CJP_yWvdR{B#)uEjjgPh3B>9rqJT(YFqVP|tNImpx~dVL{{@s0Y+sNN(~N=Z(YU zuJ3W4c}_0-3-=G~>E9LEkq>H0K96`%b0~QO&IeBO$TiU~6m^i?zb^7OFUS*ZbQ$VJ zKWjyF@tQ-9ucpiD>*Or0bZM{T!?*R+C0B3y`RZt0jx8mRLq2}EQor=xq>Do#`gqm} zT?Q!g%ze^y`Kb8&{1;uu9^vQ9XE&tRcd`^gKHMq9^}qjKzf?tiQ`ui9wKQZyN3Iw5 zLjHRoc@lhoRwIAyVaSM0WN$y52QHCw_#1Nf4>|t`eLQ~EqE&XAAw>GLVu@aK1=&)X#%vg$MW^HH2n zcG2&H&*T1X-%0=beJwzK$mJbbL;p{mAnwO2KSMlG-e2yF_~579|6K9KkW^)U?dD&G z%)dl`M`t!=RRQ{VKvwu$$)_~RfqtDxuGg|d{Z)C-Z6e~!30wL(B)2Kc+{qL3Air#q zpXD{B=UMXA{OFfa`oqf?gb%yZ&&vv<{uxWoSq$;u1=+U*=FOF$Z(~YfzwAo>P}-DM z0px_T$e%0yDSsHu`<~l20zj z`}0b_`2J0}PqvlofjAGGR_%mWPpEryw8CGe$ ziTv-qt%}#kk%$Ka{*p&M)?`E&eckh`Cb7zWpUb&)*{tNlFO=1#O9%S8or^AAXOmlW z(Ph&`^28C?pXa52Bdxj=YDb=ndg~Y^e^U1V>aW7}j4Qg#SN7lIpX&0eI(_Q=U6)Qh z$<^{0lBUets$2&B^O0O{1mBlZ`pMVSLH^}9*N-(sJ@ySbUo-UIl%k)@wKAltGGF&; z8{B_N;QGUMhSaP~e_IhR4l3_qeCQ58&!*2yB3@hyrq7=vUJUL{pEn9`S0gHZdV*KWR?qkeei@KdFGLCBA-+&ot%dAL{p`{`Rp3zJB*;eZBbtypv=eWciWKW z%KJYR@1lNDn|_YHhyK+u9mJ@Dh=`O7x}_N1Ml0=&*gD`(Uku1t%z5RZ*bk` zCGNxJr*Cs2-<7K5!|x!T-R;Zu7a#ENYmf&ZzL{Ie!+xM2HkIuB7x(9E=-YH#oL6d) zBM`q1^d!3?{`@t_IZEOFZ6f(zMbrZxkayQcyvR-8PHl?#*@#@i3-y5^cl1LU$7!B5Ds z_fhZ3OW*!MzF%)ZuKye7t$yStxiK%_B-_FN5nIWtobbN!4f2nMm}mWmT%--=rB$Mj zKO_HNxFdP-0QB=H{o+4@trESG>t&~)eyr3t^F-i2!#l3;$9dpgG5WUF8k{Fuk%w%> zeYD}^e*3I)egWCxB<9Z@Ca=AW`wuV4@9ttgX<@7DfB*k(Jg!OaUG#bX8=6%9LH_bW zlUyz6=azppIcy zpanU1VMEf#k%tsFSsf7AyP5S#j@L_yO?$3i{k>80M|H)8ChH zo|tr-+yVKeO8x2cNgBTY4RVbr*_d9tWoakcUp||!bh%Kk#Dk5^5gecV;)2l zeLfWNB+nJHJL=DEiqgkx;tgrphrDMC@^P!lg7}l?Ir;f^)VpiZ&+T>^Qg#A4BMIk~ zBjl95_LhLqhF{7Ig3B)5#7l-Lrqz()Hmzs`1zZ;9u$Lm z?KSe2HK>35Ca3K{e{UuF_OIZ4(~+F%68f>mk?+C(TUL_i!2fejlI?$?zVntm$PWFI zMd(|-81}o($sv`jQpTTL*%{}7S>*VZR++BUH-C1-e2V*AFWv|5g=MF2Z31zhxi)!k zDCSQn_kWt>JdjJ{`kYz+=5>-2mttP@4RZW?e7{We_oPJh163vOU!%$2m-KhLotiAI zONpO2kNj&G#P-Th=Q)L$DY{w@-z zORc)}?euBb&j*rw!}lkbkxy*U@g5v`X|gUC-;?j6pR#RX`q=N9E(e;BFWuGU?NG8K z^7&PklM6oArP3X8>({!JaH6kEyvLtEiv0PbF7vX|*VVu2lDLHY`WNzBt?BR1h%YO% z(BEBiV7^TX`9oghXI9YX1&d(5%?SE@J>p5u`t*5JMbvLUk~h~dBqNSKe_J2t126i! z9_kG?%Kp3$&Koj>{&w$(dcp(p(Qep3x2LZY`WW(kEBWC-+#kzHKX(f-WMUt3jvzxG z#gp@lG9=3v^7OHWG-yFTn>dd|TFGrE;XHJa{Cf)a_r>UEnFikvAn&l?zQAVk4a1Ok z%KOUok?)#TmcG3bZHT=mdD9&1ua)_m0tq>19g2)EW z4|Yq)gR@|s(Q)#r{OGTJM~*Fx`w}JS+l%o3=+@+RP0)`zlKjgH=Z7Wat36C{RO*{D z6#bw|zj)F()PI%vDRBn!Kg#@+$obf>`_sq0)|%307CCRCDeZQXR|)!e?vwLj|DQh_ zed}`%{X)t-!;n|_`*q{`!0*U+DfP{2*>V3Wj_aoiAwPeC>{1TC{X!0_h5j2y`gd_t zyr<$pj_|^Koe*+oU&MnYrPIP%KD$k*Q{_d-2(v>kohG+Y<6A$iL}eE&h@s>ts&n@j$U zdhOujF@4{CzE~xSe>|8# z&e9O`c+Qb$H8EtQ@_yT@7KT(&`crzip+8l*ubGN?Gw3&c-nuRB^LoN0O8M4cVATemflfPIt(&5uZG5>D%3i zN6(aeN^ltBR}ZeIAigZHkRQx6#7H2=EyVskm0Wcd{(WreTmN|UJJlll;{4&&oBRgx z;DSkZNI`r`B&XiOeU5wN+{i!m%|Rcp|Ag~TJ@O(O%meC6PJ#a?S;-9?@t(^L@~Il= zkGV|_ZfeTLEc9);cBr3MC*SIcdWTZqv<6^)kW$~QGTszR96!Iyz&z*_auxW0{AY65 z8q{+g=-b`fkw0opK6?oD1*N`u_&oZT7I1ymZQKt$NN)NZ?=w6hAOHOCy=D5gY8Ldz zDE)@P`H@di?*GhoK>w$epPyG1^GLUoLy-@rRVo&)vi`8M=*yJx!09Yb#LR+sB{$Vs1Y zK5I%}C;!l8;39Hs#FO@K$bGUJGEdo$7qWwoX4BUeoKJ?_BOkLzK1I22{{iQfPUGq8 zswGihK28pE#OHI+w|~oHzgmmDs|xnlzT|)BnVrgg&1ZGdzZA#y&_;&rOCb+zj(MbC z$jR=8L@D`{$ac72=0zVr>Vp2;G30;!Rx?(RdklhaPm^biK)m@#_8EtHbfxLvW7Bbd z^CDkHK5Y9~^4LX&+>Rx$z1^#p!D&$2mVSPU^Fr+~a@TIUTzW+;KyJ`Y)?OSQG+o5*K1aG-!B8Sx};5V=xiL!xJrm(@Z&dJlOC;zf&xKW1G zRtF4uo(&^-NSv2W(nm&J)Uy~jA=;w=$nxy)Y>(xI3 ze$MQnniLJ<{<`2~?7AkomHqbH zN1FUn`b(0L-yh;nKkxpjp86~!g;0O0J3=s_4#n}*A(Qlwv)FYzPQ~YZ$|#>MrQi< z5At7!s*)34qTkksd<^kp;#9K7PxSk4AiLTkKXsYBJrDf+le`k~!Lbtko2wl1U!BPr z$ZvI-Ozwzy@F|X53g-pm64|>m^7TK+&rtvHsYoAR4#a-H6WMk=_S+N5<8?#IDE)?x z^Dy7^Jl8W~v0wj8ev|;;m!@yuq29b(c^|?K{_n254^adDKfRcr_l5rpEB7}i!T%o0 z{mtbzm|v>g-%QMp`lk#1d!RJtcMc#YBmW-~L0;b&_YZcEbx%`f-6D7LMSVc&H?$j! z`?OB<@9t4}->eI{%QREgDffTwMAv1AQWUCgDDU1O0s9C_dks z-1HLiPr+ow0!dp)-gv_z7yi=cL60m_*OR_plVOpmbI4^gYcevGyt{xVPaNsvGUYW{ zs4*GcLS;Q|Nvey+&oZ`rtZfR2e4><+-qBLdy{TcD2);aon z*KgwOKCg$`A0C$KYhtl zksm9xo&0CAE^(RY+i2v+K35^vi$K24oBVsWE)6G=zbr)jTTQNjc#?F6JS$EY=MQAp z1mxdJ(6=oRKVG{SiTJdT+~FeX;X*!1`MQBXl=@=0CJ2E^6hiU)dyf+MKbw7Fy4oJNcNb5daNCNd?Nz;@p|Opi|}4% zA9BGs+?UYFHxn_>ZVTCdqeZGWq@Ty`vPkN5a^9ox`C0OTOBQ)tkiPDdW|5nH$Vu-l z(psrsM#KMhO8#WEJ^CM=>EkW%f9FZ$0B21eA124RYm#A0AAj`GQ#A z!&XhECy;gY3vG9zuYLAvQa^#*?4TxZm45PG$Un6>!Tq?P{@|zV$0wy|GWt1vKJ%(3 zdneN0-w=PIa?sa%)A09K@+-%mX>w1Q$7T0UlbEISw}$uXXa-=gHlPu9gePv!o7Hq^(nM9|-E zEp_p@OEx@pXlWl@-F2BBKt9z^mk)Eull|fIWby#SkGJXMX*iE8vZHVB zhvB?dk37Q)|MnpV&C%t!K~7kL`p!0T%e9zamrAaQ_;AXGz8${@_wlNeOCQzczAw2P z&JWpVkS&M@wKtPp9_!NUI@#?#&WC@=%YNyYuTTFr%4tXkA97LDBXUe7S1%7gt|yPF zjr-LX$&Z>LzyF;)2K9of<>=dg_2z8t$oKqlz8FK^Iu`k*Em)+aDHe;{+*2bCPT>APN9EvCVB03^xr0thd;u6D5c-9dIstf+34HHe+&t) zNsh^l`@g>AW+hBUW4~RNKJL;F^OU^EzlS4! zj3&!O%x$Rtw9JHmcTf|!AKqvCA1dEJN<_q>b zXpw!&{8+~ni)c0J=c4y4a%ns{+Z&6l67tqx7Ri>Ae!i7U6EAOa#}b+(EBz(um688{ z$o1if2c9+P=W8vnK9O9gttOMskZU16RCb`BTlCiCWGLAW@xt{Q`Cx!14HbV|gEh$$ zM}NN@qscp^Kc&D#^#9DIzn#J~sZxNx?x1P%XaPAk9R613C)Pwhs`@7SJLi1#|2xy) z4VP-N>M;4hD)`uwK3}*2`~M5%eTX+}I?>-1aDJI|o_rkfrI|Z@JsbJ0%e%?{;)z!U z`Z^T(u2;%?%MJ!vgyy@>rIIsLoAdmWt`gwNxdC+g1kNn785I>ghBD*1e zw6~>iTOfYyt3mcc{_CVKIiw=;P1DGUH4qOrl1txul>SaRz@$j2%5&AdUnbUwxPa^rP*_ntgshAw+b(6_asQEzc4*I1-WR4{qU z8eR0okEVm5e{{YjS0r5Be6NZx>uaeZ-aA>;dvef#gy+4-AMQSNo3mbbvf4 z8|K?7_kRu+M1N;K`nFzK%u{MeK2j6?ul>mbo1#A{f}FoC@-aKer@P~R!)>xY1oKX9 z=wJ6SxGz?foMi^?Gy9N_&Vt`3lcy~=WbImVgLvF0IY)lA2j>ap{kFv?@cv~f`qu3# z?iaX|pQT~mY7jZaarAdKAB*@Y{V4POuzznsf8QTrk@r#L3zIA| zK9$@u!Xig1(9gLRTV%hI|Ngb!B7KzmIFFMoBH8KVImazx=R-bs#Uddq$wP4-aC=F1 zMLfvakbbWB-6CEFxu&fq3D?P$@@n$43jJKJs3w`Glbum77=4F4rjjP#mHYQuYv8=G zhW^f75BZ%u^!46Gn(UrNuH0P1yiKy^hI;ZC`ul~aCI>!}$9B}D&UpH}epmGW|028g z(!|?Be=qK@No02Vy3J7hc_Yal!J5S7qOX^PYH})^TyHAkNoM*wFU~JB$C2?Cj3hoK z`!2x!vo7@Y))n~uCX*f4-c^c7<4?umuN4fv=r~uv@R__0VqFx+yil5J2 z1NowNrDkW9WDhJGI9{kEWJ+@G_f zZ|#<1eo|fX+6|bG)SKKa3Hcq3eEp~)w#xmVpi7vas@%uVejoixzxnyNSGfOCk^XJ> z9rHaplD}ueyyWrZ{e?_f9ZQ~D4)4dDCWqI?`;PC)mckZ^x=WvDDQl5Z)#zu38u0Bz z@_85d_%!)%8;e{oLVqvoY7yrl0e zUy*|l53(rr$|Z>wDHlaw|H1j-#$9sMX^X_wrk|T2J~XzH?eAEmVmi4v>Ic7@)7OD7 zE#kVGTn6=o7CGqWK3^?TWh~kCH{#u6@;e(%JbKaBkvWimQTh#Q=S97?GkyK3AbfR^ z{JuEq37zTh+GWt+bBTNa@us9V{T+*VlYEhU9p{xh?db1;O*Q$T)H|CXpB1dkGw*|b zrhHrI?|q0T@x|$D@4m>tDD!Lc`)gA11=l|!e)zScuQSG@pGlc-Xfp%(sciIj#Yj!^ zD(^$|jL~G?9Jp4-p^APa=CFALX}}oR$OUh4bVssHfNdqHt+ljylr651mk7 z@FI7tkNo%;@<+sj8_UTT+T(r5ljO7>crWZdx$F@1n--^U=a0htBR6tgE9IGhe%NF|8#P|{5)T-U#n_JNEo>{>cN-Slhaxo((ocVzz6e| zzLGQg8&bRseLEo-_Y*wHHdBzlA4QIiz_EeU07Z9b)US&7W4J9hr9k)kE{7;h6F14aR0Y6 zyU` zRqB^TH&~?ID*AiaR*U4erJtQ~p13ibyl}rojw$ub4M*_j@u0sWPFci%Ke=y;MY`9e zuU)QMO^1nsEB@|J!ET5O-hxdzaKY5zo{~BJfJ1w ziLxKJYX_g!r>{>UUX+_fUepiugD2#$0K_{7`gYuCO(HzV+7#r6Mv>QN5ZK>^8Z%#hG5BDpElINYkd223N$N8b-KJx!Hb*EuHrfnF& z8B5thmJp3?^pdeeS!#@}NS2Z{l#Uve7+JDY_w)RdWScDU4k=@Xls#gQD_(0DTI@=e zvggfGmNKVL=i}S&_&eG@ocD8G*L~m5w!1psk3>HzME|G`^VX{b^@bz;8vP~D$@GyL z*nbpFr=Xt{e~8X(WfJH6bYtv~x?Y`m8-;m;M;E#U`bFvE=^L0&*VgGhs0aF{&#UzWeAAA=s~Ti_V&ldZ04%)_Yl`^zA^uT8H_B7ag65 z`*g?NM5h(+&&_^$!KlPqY#D1W;Hq7&&{Ysh-VP20qt|WX5z2Tgair2JV z7C!ICe74L}Qe!>c`?-?%LOSMy64zeL=i=|sPgzO3*44x!pPp&0$y|Hpv!^ZcXc>LU zUK9Vj^r@dTiRj9FZZlYuZqf83PfcGJOACbPd!SOPevT&0;ZcFuX6WKKfDPt23|r#%uD|QhMhO z^utQ&vB{bQhcds%9z^|P*sqd$Op|*)%2X-O(SiKTW$0 zFyVcv^cc(w$_@98hfXrdacAc5`l%+F@(X<${h%+a=u#8<{f6`SU1M>-li@sm=r){Z zHe=o%Ou_olnf{1*K)^V<-+7Zfh@e;HU_PHj-!H&?;yT^A80W=5(I(7;?b_b|Gbs1j`L+*&e37m4^*p|elZmNBrE1^us8C$Fa7^Ooc9W*O<{Qd zWh~uR$9i*|&RL85wI9nJ7da_MuiO5zOXc`)md z=1ZB^9g~!#X3_gomG~O=e~%F*ofk8&%`Pi(yiTXw#Mf`feC~jHVa8(mTfI;k0NvcWLFypHhX25NFRj`kg*$;Rr;??oP(+&1*X*Ns7c=LJ7cpQy~v>-N1#FX(b2>KViP z-Ns@4@Hf0)xd-YKcfFN_D6hx7ZG`!Q-kbK>jQ&^vosRk-)o}i&=RsYD zrSbF2>AKv0K+n$9#j7Us*8Y|*X@>LoRi0oTG>M-VVt=gNGWu2p*1LW5&06S>+@kX^ zFDR?TJTAsOd2f5#4Cf<;kEC0pA9Q~q-ESE7*Ze^Tc$*|Ula30+Jlb&I_})2~2OI7i zzk&I1N5gsiiq-i3hV%Fxw&1?jjm&57B%D7;r(=#_UR^|=IfMB?J?3r2H4~n{=)QTl zU&)V-dVzk@3Oct8>&+qB4)frt_vq*b$op!{+me>JPr4KBk9n})SUS%M^WX@2id!Vs zBl`0g?8nTaw*^Fs%V*ki7Us8Yncr3sxIbwqZ9YWF_%!CX+bAU`>M);=_$X-~Ko1XA zlJqy-V2%=(`poC=!j+T<(lu8od6Z6DZ%|@wc;9`S1SM(Hnb+pYO75Pe4;@x=)`t09 z7wg0J`E=DwN{(mKQQ1nybYVU(zlWc9HC^S260>6Z(@Q1)^<-X$y;HLJ2p#@eiG>aG z+79)F<#Kwpxh4zV(Pu2s|L|vC=Qq?OIh(G8`s0`j^V+$UCV5BcB-A4@9hujYI^p>^ zfev!Oyr2Q|I@ei~QWL!&>r0FG^v9v77yOy$X`|6Ex=I)NV82lhK7WS?Y4UC>tx%r~ zsLXuMUyS*-AH5dqMU3J7Z84}na+@)q|F>0>m$T{aNtl0R(pjmRIMrj`#tQmLz3I@) znq>RafwwgIw343nSd*f|^qM!Cq~y^@zn~s6W8S{O{KByd-4pdfgYk59D_xps^h?Ym zoRjH^y>*Gqq3^k3ziB!Bs~74ITjsCDWb~&zX_=-=@A7hZ$hV%G?H{rhTADF*Ef0|^?aQZjQ0~Q2xaQ0j=U_d&~ozgg%COz{qE` z2j&63mdxw7m^@6Dl*GyJ3Au^rrU(DH)PVTg+Cn!La|^8TG@?VCHp8)DL^k(n)La z^S5R`pWm!xo1$lIS5jI)|GryE_kqmoAqSP1cF|GCl$2R8pC^1i;Xo_N=uJ{hA)PcM2H>W7Al>5?6q z?A}2qAJD}53LTn`-}gP;KNIy!3+8QetOu)H>91H1JZI2PG4FQWMBlE^WO@djh58`l zIi1r)mp%=dw-0_mJ@6C#)*jFAfpopT==c9dcN>QJ?%%YJw=QcQ)2p#x$=ZT>n}K>j z??IO-^!q2%KjJ*a<7nD%3+k7H^tvQn67SL{v0v4tD)YABj4uBe?#C*}e4v-%eyj%Q z2gMrh$NK&S&T|{?$7)@MdBIhFJ)3Gc|NDV%WNE_l8S~ftd+ZMyM3>p&{K5>ns3-2@ z-bCMUG091xFO0YoLGMxWiML%6<68|G~?L@qP zEuVJ6{?IIQ=Iy-WCTVO>pSpl~?*zKx4V-r~(Q_Z+yx$(W@D1vto3vkrNqj0XpP$#l beRS>V@TRz5Xe8aj_FMm*b{RgfrSadtG0%oM literal 0 HcmV?d00001 diff --git a/modules/library/referencing/src/test/resources/org/geotools/referencing/factory/gridshift/BALR2009_copyright_and_disclaimer.txt b/modules/library/referencing/src/test/resources/org/geotools/referencing/factory/gridshift/BALR2009_copyright_and_disclaimer.txt new file mode 100755 index 0000000..19a3798 --- /dev/null +++ b/modules/library/referencing/src/test/resources/org/geotools/referencing/factory/gridshift/BALR2009_copyright_and_disclaimer.txt @@ -0,0 +1,25 @@ +BALR2009.gsb - Copyright and disclaimer +======================================= + + +Thi grid file is free ; you can redistribute or include +in commercial/personal applications. + +Creative Commons License CC-By FjgMatesanz + +In the event that you wish to alter the grid file which +appears at this ZIP file , for a purpose other than +personal, in-house or non-commercial use, you can apply to +the Instituto Geográfico Nacional at the following address +for formal permission. + +Instituto Geográfico Nacional +General Ibañez Ibero 3 +28003 +Madrid, Spain +Subdirección General de Producción Cartográfica + +Whilst all effort is made to ensure the accuracy and +currency of this grid file, IGN and CNIG does not accept +responsibility for any actions arising from the use of +this data. diff --git a/modules/library/referencing/src/test/resources/org/geotools/referencing/factory/gridshift/malformedNTv2grid.gsb b/modules/library/referencing/src/test/resources/org/geotools/referencing/factory/gridshift/malformedNTv2grid.gsb new file mode 100644 index 0000000..0665830 --- /dev/null +++ b/modules/library/referencing/src/test/resources/org/geotools/referencing/factory/gridshift/malformedNTv2grid.gsb @@ -0,0 +1 @@ +This is *not* a NTv2 grid file. \ No newline at end of file diff --git a/modules/library/referencing/src/test/resources/org/geotools/referencing/factory/gridshift/stpaul.las b/modules/library/referencing/src/test/resources/org/geotools/referencing/factory/gridshift/stpaul.las new file mode 100644 index 0000000000000000000000000000000000000000..f55268d7cf89b4218b6bdbe3ffb4acc75b97fbb8 GIT binary patch literal 3696 zcmb`}_jlA~+Q;#LMyAhXdas$@6k$ODeIr3YG!#KPN=K0*1Qet<38A+bg#eKzgg~e( zCDahIf*=v3M9PXR%OW6%6v?vpKJy1WzwDgzojGTInDd(Zx~|W4XE)7i)+XECymN>4 zO`COSp5<=eyhSVV#sBxpZkpAieU{rDPD@KOh^Ll#($b=bjvsroNm^Rp*HWMV_vQcW z?QB&YR^&?9)=jKqa#>RG6w%j4D@Cqt9Uf5-s)O zdb9@(XtPZeo-(tvt(D?(8zbTlo}cRE{81PEp19fU_VTX=KI&&>GNXl`P0a!%UJmk3 zEJPz)n8S4+g09w;c@sO0`Y6)#2AZ23b?*?A4)|Dq*#tB$nB26DC?!vb`AFUZGlh4`{zn7-Z!?d4H+-HUSnLJaTWIFpMLC>AFnVxcPl1iBHSi-TJQu43J=yG08y>p`8mz9*< zQ_;|%p`t|wug=o)@FyMFeghBEjQns)IIo$w_KbyrCoLR4X=P3=JF%7ydd+t7QL&4j zd)yp6=;6WvFAaD5_;XVx&#m;6m>aYd-?lOA6s^3(rkmD z_MZosF(pWk!6BNx6K2jU5vB#AbXCP@{6`Ef#u;-YA?7K`hWwOspkKn3fRsTo8DkpA z+0jbD<-SVtW~w-`Tg}Rc8uq`U#W_RA$fJ4=-!#zns*&`=CcLZ7)J(IG+tbR~H*6HX zXJ>e>gBztz794X?ea=m-pFEs6>m_``$EX9DY}o8)*0KQK%m}i4c!+J?!ZgW>Fd`OZ zo<2tDlNh(J#u<1#fo?~V-wIRCe)T03HZmipY4;wJ7(l(l`9cr7|aKp^7 z6BgdtVP(}?8!xW2Q)ioljfb4<6wcsTHzUq@h@JNG+A$xi_hsS{=h?F`z|EOK>W&IA z`@Jyj-iWaA#VBX&B0JOL=x)b(>uiFJ2a+T?#k(zw^@eh>sb&GD$1;^Tw9}E{VE_#d(euPFw2;uhx&^aK*SHiL>}_oOaa- z5=TUSZcaJ750|k2Qz>(EWtf)9slQV}tFubZKUGl@(~#3WgP-%X{JKLY&PLDkcMQ}# zGScIri8Z&({CwHMOI220uC(#D!*+U9xks z;Uxh+o))B1^ru_Yy+oX6sc?37$2ca76TC0x=u(1VFLb_U%?X<4#F z$5*9#zAHEIuE@?Gel)S{M>G4XEcE@(N(XVC#pQOE9&(U(&`I6_7hmsl6WimVXor{E zB|fC9GYJ*=8Il(ub8L{$1_@`EFil%TNUs;AwQx33#Kpdjv+(DH$cQBScBP!ogQ;CK z>9~xtvvMX~S8(Tvl0^nJLtfC}ZIeNn*p)p*F2z>pIlaX|={_URR+!jaX+~3NVZ{+K zM~7_O*>C4UnFHHiC$ZfwhW^z}pKTu2mUzin>!U#Iqc#h~98C*wYgCYW146tn@^h=m zrL4M9I*9(PEu2Fh#@T!|!Htti@mncp%S#dlUy+iyCga^(az1^mV2MmcU#A*leGQ+z zl|jS)T8hT$=)6Er)>j73Z88%2tBKWn%ur_Gx4l*h_SmT2Wk>dn1NlxTh1*>8Ep_vt z*n?H<&6Uf2t7PQa4Jcy zeJSU|ClXrLNU3}(qi(u_UzJL6|51?}Rx`bkh6dt1hxXC(=O`VeGxY4tH!$TZ(VxX8 z#+91M*(P#+yOloMZ4BFH=hhYn*ETzOxx~e`^=>lOcyO)ovh#Bvr{-qz!8AW#j1Dm6 z!yxs0hPWj5=323jcEzLA6}i+z%u$V~`yAm^S0$OUKjqwRl`_U6Lur*GaVR+MRWdcI zLfSyh*(?n$I%F`Yua@pPI_gc)Q+J+$AHFbhzsN-CIx}sGEz~WsvcJSe>k>PPjSez4 zILTV);LQ^csRY#%j8);ew~tu zevF^Ap#ds-2btb6M0SfX%^F6?4MoW^$H=V{XP%h%_k}a`Lz3?frkt;Kl(M&-jQ(we zGh4xmHcCp{t0?cJrl_ZetNq2^9Hteu(s6&Xo`rb^%nOZV7nraWn)!I8Mf8D{E~{+# zR@u2z;gOt!I#a?Lrha3*)RG;g&y2 z17i&Rb8)ta{PYOtC&Ib+P|A7cBT@IkGHTwJQ!+q7RevP~15~7cpyv0%8WJNi=rmf( z;|V&-rt68!HgJ1^k)!!03cfV+PN9WI%dN~FTlMnL57HXvY9aHbt820MTv@g$dYH{RQw*NQq21n;aqw+<$QLI zRO~PrU8f4?L@UpJTPpU^ zGAmt|*;u^P&aYoM81T81yNg`RoaZKz=b>V%mpkLX#MDk5 zi0_n`5ZNKIRgzrzU!PjEstih$tx$l8A_?@!ihra9W6_^Fs7l8KPuPh*7yAGOtTK;PF{4kEBG8n@@Yp8SGIo-XnUZNB1m` zFm+9zU&+Ak=4j!TT72>ojh{WJYt4=L{0Lj znBmcWiASj~Judy=F;>o-_6KnW1 ztK;A0{5n-puY)4HjlxehLrO;Cg0=uBB;nJ#jn9fsK8<_$eA~|_`URiN zAwF)T&%)6@Q^)w6dd27G@xnXNXZ&Pe&Mf3?x>)hk2a3Ffied{C=du)A-&O3&Q0$!{ zj1v^y2Pp=pD9#J}SNQ=a<_7#sz$EcRxBWiFj`@r}=d)7u-X!xmyL`rqS2DB3ClY$A z&!Qhhw;etucl$Ki>*F2p;gHWihka@v^|>e!FVA+8Xd+QTqKHIUi92;)a!GON55<|w ziXK-*(*i~64aI_+ic>P4xu&T9yCN!A(R!Pr#d5{HvlXklC^}n(J`(W6zJP**0r^`3 z#w`lSm@K*v3Ao%RpwUxOiw*%T+ll9s#5=75?n?|v7XQzcT#qype>Dj(j|W_D6mY6R z!0!412V(*v>II}l%P5}%jWvy$Xkemw9m2F<&J>0x^EB)J(vo+l{Vax^#?9!GQZQb%h0rtI=t9WlXjmbUpzJXpdxOQ_*^{GcAjF?>x$vy z6uFYa%l#EMdn%fUkIE)1rnXeP*i8Dmp<+lq#h{w%F6Ut3+4PX1PBp{%7{f&jWm62b z(hQ#rH58p>$au>zV1Z%m$A;dk<=$t8wB@q?Jwvw{h6bYy?{za|#TmRJhVi2L`PG`{ zQm07q&9Eml2Wn|PEv1?0X?93YR5_)H-ltf)RT$SRHm{OCUaGifzVyaBg|&#Bs`yRv zxufU6E{4Uz6fxK^Z;WAjx?%Y|!#3dVuN(3& z7&;#m*7b&;=E{1>t)Qi$eR-+HRZZ+(P4Jm!Nrl=Y9uGx8|DkXwklsre#UZ zQ%=+MmZIVn#k32ex8$);{PX&rLeBQWq?QPSXtZ&&)Kj>63D<@{49kRNMoG)ss+Py% zES2Le_1jwdwznKfw)Adgd8o1FlHAX_$C8+Dcw)QcxX5sNtjxO_io}X#(|mQ5QiUAM}7yXBh1 zuip#PdP})aEFzv zL)%h@wZ>5Cy5`y6G(AsgewBC1pW8+6&4ru?dRw9fT5628+_wZxGD4e%a+KqqWeBef5~I5@w1u#bk|qeiWt; zL=W#vOOb5Lo;{X^CoT83@xQs zazyj|@rL(@8#W7nZ1+OWv9j++!qn!P)XsC%E#~M`)^WC~V_1x1c{4|38%Jtq$BrJt z+1HWyoZ}e@ub<;WcSofTj{WhDUiBPt6&#l(w|bW>OAbj5zOnRMVo86~k~vEJpC-N6 zPW&QX?psaRONsBa<&b!K=y~bEqlR@m4b{Ij*w5vCxvY@$L?zk3rlWD3V{G$6rt)1J zZPOh62RfogIie>zil;mIP3L&$-?IK)M;8fX{`*W>H`TG{CC8uvj?7eHPH@btBm0$c z6iDtPM2C(CEeAJSs(m8db1e&}SdL2%_o_MKkCGbJ5P$jN{}06i5?z)%CaiIsUoRSDJ09KX*nU7XJ}Gh5k@1`3 z+GWRs1&&W{I(puc9{R_Tazk>s>Nt?^xc?W&BS#$Z!oNZG$^F<-`d!E8qVzt4BaW-751H<$$LV)GXGua%{K};tDJF6zaV@4 z;TR_z#Y171+A!sdhe<9Src{M6k5vxyYE+n2wZe>x4zs3im?{z<)eiGW^)PQg5+>rI zFq7q(PC=MO!oKdD)Zmb#U5=ycI`NLwrg(;<@?<%0nDF(HT8XdI;vMfxU-p)MnDd~c zQ4#U@Ez9&jER}w-oIPB~SwYU&A`JcIeA6Jzdrjo*gfL@Thk31im}*_Z#B>WYqF0y= zeZm~;7iRu5Ve(|0E&spXGfZ}Bm{uuasy-2BqC7uAa_CVxOuJGtHYGtzYaTE riEnbga!e7u-x6%vc^!7%p?&-J3k$9POq$&xxf(d;6Ab~0hGk@$hv#eDjffTy`sf6ki_G&bmh^M8e2tsAwLHX(^zihV2X%Z36Aq^^wo@t*3sbpR%1e2# zm*Q8wtRLv3{5l^sYx;SV?oY0zYbLkbW$|)m7GD<1W@S<~>QmW-_sC(~(i|$4&E>v7m*O{aS=S^;X?2jC zNkOU|4D##0L8f)cBi)jRYknSqGkMG^mCxQj`4sl&6Y^s|OK#@#p?U#^umTPaFQC)* z0v0|ipiGNlcy+_LHhmcVj}2pak>Ql+J{+}cI6r+qocDhWC#2#C{(IDB%v7Q;s-lle zMWa7d6d9=I{(3d%>T1|AM#E1pHOw|@$v&3me@SWSj z!9Of~7?4EN`Xn~iNrutMTz`?w9K99aF)PX4Y-E0CV}5x%C$sE)e$P&HjDwYX9K31c zWWXFJ*(F>oOmVUDs*6d3-9&D3bG3nouH!vKy!O!Al)~=cQfSyMl}<}ic~r(rl+R0w zVBKA?-m=EWlWKl0=lU7*m!HC6Y4qNdhNez{h+zTRJP1%IGMz(P(pg?FgT*5<*#9Vl zk};X&ZqKAg<1Ffq%i^0CS!~c`lW`!MX>D^TIxC06g>$)O%a!w#%cS~2YAJ#o85?BA zo*)IUf{baC#|l#(S7ztY>|`FfCGxr3E1!had>*gNC*yiPUsNk#Wk>-GE8y}k1^Ay9 zP`3FnR>lvbhhTm5=rFqFg!GeEH3CRn$*Jp6w+XSP=oI!{T9k}5v;s(5lu#rYv> z_H9#BGRRdsyG(=q;$j*oGA4Ey!`)Xu=8IR+}1 zFw)jzq}yeq{QOLGSZ|_1odk-HO5p9Y1fFTlJpR?p`;LiJT$o7j(iXCP7EWKc5Hd80 z-!~_b*C?5R6O(E8Z!*m+R(hYYlGMk>cR$#8Sj{fq)y{$^cK+5l=y%w`$j(mIFLtu8 ztcx}2EZ8z)on^hjmh(p08@?xlfXt-Ru;_cuN+ z{ph1!Wj`%4{XD+mCpajL`>WIFRwaNtGr*!70S*mJr~K-4(kf^0MP>%qZe(zKP$pll z&17%2ESBVEap|utnRnU5Z_Q?7!yNqMa+vfohiZmgJ|4|wQLP}Rh#(OogT(I$GU9oV zV=eM%sn27{j67QXmdD|u`S^O|)83QM%jNlOy_(O+Dg}fMEucn00e7|*upp_@Rlzz-u&%u96V3^TM{SSXHiE3Cl(lZR?0{rpOLcH zjg%W=;`U?x&Gb{}G)5430Bo6IJVnd5$ zwoFgvbYUx{oK_+(Sy?u~#%CLBc_IM{#E!OdPyuB~vgva*Z5K^F@i zxY(g^Gii^T(ycvwJ;OttLMcp2PGQOE6cT!-vRAO)CFeRJ)yq$pyv*tAqvr}Algs(} z(d%d2B|ig5uqc+1Fn_5>JrEd>q~KF}d znc;5*YZ=>xB8_phI37pI9!eg5ujEVx72C2^EWM{z0v+p(cv_Vxo4F z1bR$LAm!f#HYJ)VdfJSqcOoxVBr>bA1x3(8lZO_n#3fO8e-i!MClj2P%*B#cdibpT zc*9DwP?_Jq*eKJ?j(w_~Js<5nu{kJm!NHAwPA05zQmmGXP9t1Ydg@}S+Rf<$Zkab8 z>d*Gjq(}-ottk{alforA*C8^)qvXDR?efz8oR>PieC%B8qseD}ly*P5Q+~p`q|ttE z8f6Ow_>d5w*pUF;+R6N$mQKpsbjIs5D0E29b(>7uOwZ)WyG&Y}v-tD3ENXPk=IEkq zik8V?i!X<3*KQ`!$cbAM^ON zOFl0g`851KANA#YCR8e5%isd8<`qy(u&ygucNVO>3fApphS!t(`n_O%U9kQS=brz@ zQq2;_%hPdO{7T7zRZ3P?RWUAKMc{80u_`rf52|_7LBp>38oZ^nGUv4HxvrHN8!z)+ z&UPam`4e=Ud84C}SSy z%}j~&W`6iO5#Q=WqH9`YMp%e{W?^)E62BcuqGQ)&ep;GL^YT_UWC`B)tZa?3(Q&Wf z-PX>%Id;l_=HT;G2Y+03U>odY^UuQ1>br;@>tg;tE{5yey!zEm;V(TbnB(EFU_DZ> zzIZByjdEXC6mIr;aWB)XUXGpc@^dF2mf1f3dhe6j69Ch&YB15?1;!<@Rkhv)XU`Ph)iz$ok@Od7AtpWVQZN!-#MH7kJ+5C<}mwQ z4ky3P#l9w&PUV8!>K|lHW{?HzgZy+a$ocwt)QHO?d14;B4(8GDL!R7C`8=}=&ff`7 z`y-zQf_0`~eOj=tAy{h#>+yp14}$e}!Fs=7eW>UsoSRIGrC6aj?%LvHCdaX%pOU$2 zl;qS>VH~NV+jA9_b!sjhQ8TlPhW<-5+$*akFI~&$ceE6Q3oqLi&)3a#jG3b2_ncW$2911tV5Cp1iD`RHJZ+PJZFU0x7BdSEGh@7L=G#Gu4BnVX zQhf`H##wms+JefQBv>S2|0&Jytgd3-D|9A?Egoodmoyv|6 zsq{5@@g4HwYUQKlL?8Q}`RJk$oVWTZTPKb2!8H2Zmh(L*fbYiu^(v$@#GlTSKhkOZ zwao7o8S+eKQaO;xqU)LL`!^2l4j} zGA|J1&ow~?{uN|j-8_aU^0+ZRkE{cE)PA4GwodtI1?v|}<*pX2s|(ii1?y@u!$%3$ zH~yR9!p*e8!^h5$8UEWaGVq?3AEVmK%`OPm{S1s) zZD3CgBc+EMQ9UuTTWzB2ArnWwOdxH6Xs*&`8u-mLylIBeM27v6$c?5JVy9ZT_0dAE zJ&9HqleqFtG85J()1#i1KgL?IzP559!AA8{Hk$OdQ*5Q3d6ga9%XM(^zJuf#C)0L0 ziEQCw#1t3&>tdM6&A_8>M#{PFATvDrqlbG5DZCQRHK}bX+oz;5?qw<^RbG1S@Y15b zkEi)Q=HBtqX^@{wtNbi0o5qjsG~S*`$Q=MS>>ry#aF zLE6{Jqfks97slnWZoi!AcX{Xp>o)%f>n(z{Q?RZf_w{D zPdL9mpx}=VvFw{4%Z5^MqEF-S--tsQqNLMSB}JPEFPo%d;yV@JB&m6JM$Lpi8fyHY zVOCWwCGxenA89!d7f*%#@r1M!tmo+1U0lyQx1Nrd^|%Ha*u36Ag}O#kMjN^N!bnhW zqRlZAFS;fO4@=)y$|txF?iZyIHq2Ix5^K+(qm{L$&m{w19e4KipoI)mv?GuW-lWX1kW`Pm2;pOeL$ z64|s(6@7g@n>Qgj)Z3oJl9stlpPh@VSdcG!1*zl}o&95wFK!2^Q#+52(RmCOtaXC5 zTd?*E)^5Q%QD%6A+}ACHn>`e)7YYyWvaLX7Pyyex9L74sC!BA@D>!gO!OvY{S-d2c zVddgT%7~-ipK;WURPuC(lFcnt*r%(gQ&`PGyBgyKHTU{y(67;Oyr!0RBeVpcYB{Tk zr^c`GM0C_KZh?+trSw$x>Cs--vwnzy+M5kbX=tS31S3b@7@2K0;XP#%9+tq=m4bH_ zGu`sd?E2eGtU8fmzb10LlZBraTNqz132$Z+L+&O~B_^4Jdyb#L$W43U;UgxAzJB3h zVQdPnUs5@G!|D1 z@GLDr;VS`de4WnRmFY6iGe}R*;OCne{4z9?(VJzKH_9SzVivs1VsCOb7tUv68<2zj z=Nv9J&J|vlOGJ^N=&2y9Q-b`nD#$mtg5vq)@i;1vzJm35!TN|`{Z_DUEZpoX;bwz` zhj$UpRbBM;CDGSo#7k-aqJTrKKH+@yRgBDJ1q)6n7}GNr|MFOrmEvfh8^>o4;y9;J zGHv}xSFmWHNRX|(|VAGDH}9Aud9VIS|KUA_=f~T68ipI#w9SC7nUMxljpAV$d1EzE>AXqidxH2B3G}RO=KLtL{99&B#zY33 zNTf$E3td-Q=vXz0dId?me4HfMB{Sr3GCMk3`Dd}r@3J<$>2hD+wlOQrF5KSE%SH~& z;~j*&lye>LWZ*$3A%gXqsV*YqzWzGi&BVRpr!@7DJ6d@7Lk}H8#Yv3Xt^Q-zE7HU~I+BtT4ubT-dQ z=kvlDD6JWc7u_|cPbQilgo}&D?lnA%*3Yw8Wyog3$!vO%LxZ(B!cB5{J}H;$g+y0( z4|2~F6pkb3TCo03utscie1##?p97ncV$(M(e{PCrV)CDSv zmR6JRQ}g1w8s$(88#ZhBtdW-J38K5+XnB+nPv?{I;gm4aTT zUh|^w<0EsSkCKJ_v^My8vDZ(aNgCHirP1b58ve)t`?dtAAzVCjR64&uO{a+_1HI_3 z+>V)~F3O~Nxh%$KWr@d{#pmj5>Kx6MbC!eV#~eI$a)sCD(onE|+AYWd!Fsk}?Gmhe z*3RRdV7)@P*-+u(H-v|qy5@6C^mXkO`3x4{OeY?xI=p~EqY7xU`xDONTf{JYdJH** z6<8e#zP+fR!#A;1UKh))+Hrh0DvtOUaa7VP*>Y6Lpsp&;EKw0zPR;IgH4W~l@rP?T zu}wquW?Jl1wCs4VrDAeCw%_B~-$%!nKj_$5RgWTHPno|3Yo$T-tbuDC#6wwN#8leE zMZXFEEfdYc6L`NpftM}Jd^S_~S+PX=c@nW+6@E6v!sIO$hBpy^bV?HMJ|tnVB{S@N zvdjZ3r&e3JRl`Q^FdLKqwoxa}PSrhj(QyuHOmf7Zhvol?ZBkU9P=h1cP!c-Exywyc+xcCUC&UM{rsF>#EKKkxe} zG1SkipZwe_pT;F`8kaAoQIr5fmj_t?MLKQL(^-2mFDu-je)ZgUgNo5~7!q4vf?qg6FKf7l6F$h0v zZAhc-p){;*16-UT+Do{2D_c5m&dJ&Cm%)^E8N%a5TTjen>c>ply0eJBk;Rk9Z2sCW z{#n-?Ca%n(Y3*Dt3D$=M>soSO+fsw97jE`IczDHn;*pBJZX#YvN%74N7ZKmAXFiX_ z4-fk>pP9Gwg=Z9SEBX`8QR-;=9g2pIF?3rPL*p_EN~S5ebyLB*&{)!biKRo+I8IHD z!|`5bd6H7*yOO{AsO0%lvA3$4O8KI>9;w+Lr=jeA4Tesfhvc31lBl;BYrHHJ6!5 ztdz*&U?RUhO5~!-!q!6;teuiLvnYvYpC|J}Aemyflf`$l61vq&?M5=Y$Jw~>!X|gQ zogISp2RYXsnc=lwI4Bh7q><=ra~&7ka)gIplUe?io5l0pT>IC}Ql*D4HhXAPIfY`L z6t+wLVQQI5xsj`W;b)x|iih^z%Vw>QNjrT^sqg1{fgj&pKQW=Avwu#bWZeKq zM#>z27T`{NI(3evOOBO+d1VICH8OcII#WD+xwGwAw7QnX)(G)Z1na4S^$t1L4`hb7 zoSaK*;bz-=2I=PuvT?QGd{@qOgFHrwZ#F|b(t=}ogcZ-{O>eoY{rTKnozItl$+@og z3Fr1NqG+R!7L5}vIZlk6tr#ACq2Sj{1vCCqppA;9;;vZMw~S-Jj5zU{m4rK#9J-*S zLw}WcDk>h=6mN8dn%Pg)k}GHk|5Y@0M=h-vYFSe{o`L>&Uf+mkZ-|anTXk%1qUX=a zdis3O^NZDh;=FPR}IfEl(n%QZlo1lCj@Q=I1CYi94)#o7y-x!A9y!8((YfbQP?1%^mC*=b-Ns zhnycLKWucOt>U7N&qa~nT{QaA&BlptPCj&#_N|AXzZV`Ze#!uCisV!ghDput!EEtRa+vZ{4p9wqDLOru5Yb$pi@x5J z7G%)6Ap7qJso6M>2yGrt$v$)^@@Q8wpYwD~97MVk8StkQ!9*)k6i(V`G`MFP3KQ;>1UaqeBTL-+7c&yrLv^kcwY7sAy16 z?&&dV&itdMvq8hsV;Y)v(@G|!Me#+vpl3gE&wX$xvl}F8O=qK5DDp>z_uII^py{3snIFp0&4;*wF>SV}| zPHbgd{E{qp_8}KvHFgUpb5n3xeDSV=vt(V>A9=706dt}<{L*(R+>A}7@y1kAD|mS< zzW9G$>6`X)rle6x+8aQ?eV6ptrI@#K9Jx09o}a5kFFePiS~ zijldZpyF@^$DSzgYGSGIYb+Bx#!+rz9C@XcJoYK&UQzPXP|2;gs8BajvtxpqhHvF; zn>ExwrD0<)E#?)Xu`9<@D;Q7xhw((k>BQU7(W--HaL4@@gwh_gTqmY2&J3{olFnqPEL(X6Ir92crrcG`;N*uh&V$5+_Xyx#*#FiNEP0 zq>7uaZZ|&hQ))MsoV)N>wowo4&nbP9!Tr;yhtmDX}+SACmGg3?R9U%W)u_R%`{ z$v;~u{7m<=s$tt%$v*mO3HwRQ z&>HbrhsU$(shsb49Se`>2=A(=-FJFARgn3eV_?n$1OLVu$v7zdtdoiPi%oR;B0;i{ z1cuxfjU8*|-F~y=5s8w2CK6i4f;-JZ+goy9hbGZ_t8j72u8)jM77Z@COKX*0gOysX zY|Is`+y5gPTX=Z;ZFU^B9n8scNXJ1kkRDF*W;*5Wc5*(<#oASpa}{&5TkEE?^c8xS z^Uy8P!?>NIwX3BtQ10wH$5N=@Je9=ZsZ_b0N@Ra8RhLP&S<1&br;p2LeXRM)PtTQp z604+%_D`d)WFSk70jm8TVAz0kHVM|_+GOx>vCQzQnH(RRDVjTrD1R2qMRPsVXUojW zrrD4j@s`CS?Ie0@WiG4A1z9XPmn~m%&Fw*^{1fD6+dN7nAMqa;?=&nW?5uQ4L3&8YW%RkTpQEkM-h(*Oj|^Ogw3? z;;ELPlYdJ*(%yRUjrCNhZjk4|z=dZ9EP5lK9XGPEr-`_gCJI+eV77Epsy$10ee$Teh;GX#5QAw2EnM9yjGW{k=FXeTz{Qp+g95#&zN0 zuSH+q5#Q`!6$cx=4mO__jor@4>#Y?m% z4{asK)~ixj^K**$1F7_sJNwuPIosl?44Edqg*RSIdLR13K00>tvvZN3x8>4=L!>eK zaT*H@0hS2XHaXXKx67x4NcCDVgG%-MX4h-DuHVhjTdF zJD0)hbNQ-LkX9i<{vIKnl-)AJ-vp`OQTizMJW73^$FVEokyibLv+YR)*_ud3{2EC? zrzkQPMWOyYnht?z-rb63XIKobZ83b&Ou-LR6m_inU!OEEmz6~m}j3M$W3Fttc5g`KguF2-`9e;i+}jpN5!N@|W&GX0s7 z>haQ3IIN<7XElWvtGV#GhEo9zk8f+}8lh$J4lV6l#!H_pp0h=D{N>hB@v08%5IuLc z=m|6xy*16i^TI}UxQvXvVwC=b_@i4*ENhxT9m%fa3Y!Uc%KZMrOz6Nwx^GC7Sz%%9 zXbUHvTS(F-(dS?iA#Ib%lTOM*=^yC~R#wQlUMAcuWsHpu4}_nIAO5rCTuG8&-!Vzo zY^Q^X)t!>vIYqlUxg@%4Osb1754(hGx*6k>&oP@VXQd zhonlUNVIk(FWDJh%HH;}F3d;!9X{e)`zbMBxY-wJG|Eqdmub|r1(&5(1FNzLV%t6Akee@TN?P@0)d$Hsoa$mn5oyO3QX>{=gm?9m+?#6T~UX(skcm`JwX0Wt( zCZQWM+1xUV>x;7F{+5huoM`L9;*I)qB!9@Ew=P$-r=06XK~897ZqEqP|3r{pCG)7> zH;?n_pK!kXOE}M(M(}D%g!qUNys$=c{ahs5zK&wb>L|i%MpJl1v~Z_r(U38m`!xn- zCxzst3VN4`Wurfq8aHDZ7ZOK_t&(dtR?=XSlCA$Lv07Bp-&4_|k6OG>=^$0pu&+QP z+EpVuTe8fkP3PxPHMg~7JiWh8R z?XM;pb`if!^33Mt&C)S3bM~&8e91G{?UJ+I(jsTx!l3sSzL%ati{Hdg=_nY_PG*l_ z{Xu4UL*e1-#-BV=|1BFIv2%2uo%1j4^a^v(ahZdL?;PYrI9ajG$?1QcybE_xVVR4X zl6SQVcQf{TH-CtyvM0*J&7VB{S~`VI(oxwUo=QQ-RG!XFWqFBD-Luy>yfl>zP~0NU+El+&*11;@ku4;YPByDefKP4*JSZ` z<7_(25pMQH4&N8#h;J{y+k}T-%cVjK(OO376wa02@bBWIl+B~ofKNF87!}UWUEv&U z6~W1w5z;Y?WUDKZNqu+Hv)T|tm)~Mo+f6~U?-We`B9_lG zW6A$BRywcJQP>ej;}%MOo2Fz$Ar%I@bi*&Gu=Q87e63pg7?N9$mY%{(jr44_l6PoH z=p8TKPdq)V>$pB#$KvNY0tP)vC-m}n2J}A|Fi1}!W~7nU{}{PvlsxmKN%&g=dsZdj ztR{I!f#~bUX7M=_g>wlPZzsBXwuN(rlbDs9g!^O?UArW6a!xXa57G}8tf%djOjCTb z_W9Byy=J3FPdm@1*(v(S&i;OqU(a?>O)~86{hc(QEnTEXPPPj_YckzM3-L-dJ>A0V zq@!|2`omv&aLtt&{;!8f$+F)aNTE}^RLPrD`S&x)yVAVu{>v-9W}n=TKA!cJ9DB2L zR777#d?(%ES^<_$3GlX5I+1zlOnRNp3Rec%e`Ziwn@Qvu(bv*T*}GHj>yF}~tjMNc zog5M-=TNzLE(w8LbPsZQ(>6%uB=J)g3;(_-`A5Z1IJ+N)F(NLUnft?;*FJ*T^CB2s zGE(MEBwepZ@@7bsDV9KPEFbU3GDZi(;vX0o0}=k=1S7#Og@z>o$;`b{)a^SzP1l4o|hU_$qe^a|G}u(qz`TVu@f z{K$Q6NMz&DM3OsOh+H7Lt3(nL9Z4)Wlf>L^$t;sz;qi~jJd+uITCi@=(8hGhxhno) zV?`%BYscA{dQ0Ybcj=^vPr5=hcE9e@B^~D^{F;*j;b-@xPpUg5m^YI9I?c_R6K=kd zEc?H%_}ix*ep9CKLV630+NV-vK`P@V1KE}1CF`k|>1H1%F8R0?F53E#pFh7!Lm@r! z+ie0^zYp+N?R2_Gzc92+2IcZIm?!%HJDizh+|HEeBujcu(nFFS>7C8eFKm^AYGDqe zD&`9J%jLwYTwZpP?x-WkqGjTnT?_t?hjkksMlblPx|56h=>1ggam^erp*bb2J` zZbvdUJW8^=D1L7iP14k8UVn&|9gY~@o{fYCO?i9FCvcT2jkfJ zrIN)91n<()N%2bu>6YlOFg5MBOE0ClhAY!0yDp-oqjX3YUD48MNWAc{c(ygwNoP*S z_9A*tODE;URXsIB4D{Y+Ahd;a!e<&OTHM5UDJEWDGf_Mwf%RJx_}E1J(aC1Hhs~U@ zBy#0cA}xDLPhqKr=4Fyt;Yp&?dFiI~OlI}MWRir2-b3EpDO95yIyqH4mVTkNk@FR zhwhI(|I=MH4y6#@IaRoADz~b788pUA%a301seN>R=wqA7FFh~0uVw!u>Zs&g{iHYi zOMo+yUmstQF8$>U={9Gutc+xtd6~q#%%raPqq*Xbj*QBdyjb!~+51Qrf3${V*P1Ef zhl>~fujC#pB=_hkxyOHe`&FGVDvSw}nH)xYV>oS&htr@(1f`ZmaJxbztFj|WxECoM zjwpWK9Ytj8Xs*l@yo<*0z!`(}QVjRLQQ%vrpirGy7LAIPeS}zw>f<$!8l^k5E z|)b!jdc(>E=?K};WOKB;Zre*GJEy_sAJ)~>asjYaU^JIRP)=Neq ze)yl#Q-~40wa5%R(E?Aqb%oD6f)U$C)x(atr+bG>idg1wY>Rq&>YU^NJ zp5)h}u^-oy9+J(;aPh?peQ@&2x6&P+Cuob%yMKJo8iBCb{O#H0js|q-P_&VbRv}cc#m( zPzH5YWiUi~q}8Wok|)|aOS0?kPqLV7%9ieEHs5@kBblds+DIQ|VJ_N|L7ouAn*Isr zU7JEV+c1m=6T*1?HjIA~!?|}loU>m=uyIucxm6>j%M;1pk0P0%jH2a%DE79G77U`f z_*o2^)EMqw6^%VuLD5YLmNyU{K3=kqH?h(;h~v&lIorLIG+n7Aw~9)3I8>r>RnlEh zlkuyXS)Dc9mJX?^yq34pE8P7+%W|dUT!%zsch*U+q9a0j3X6k!{+8^zgVsROQ3E;M zg@Z3Q(x{4w#rY;KJT}R$O9CwqC(x?1neP{wB~wV`lrIs(wM3o{w6JuIoNei(yv>r_ z`lj@d1|&s8GF^q^w#2~6o@x|%f;QnZYGO=)}xU}`avFs%6*+Hozl`5Qb-+~%B)?f{MALW zk2PM3wD2))iFAi+`8hkyPw%p6;xngVekU9}RqpG10cPmZsdYA;_JgGxzEe79(o^_P zZ>_7D#q04|+$@yM2+`K(?#L%LhY{kVJm`_j=Rf9>Bc9p00iSTr3=3uUwoo=S3uF7# zFwquath0qP|9m)E{UV536G4qyk{#?oMl@UjoFyib;n!nrsGf34())k@k*U-;Tc6{}y!$DrooaWzePX;`#U z!=P$fst%WK;d8BY{^J>ZBAzk5bv%%sLU>I*>qhCR_)3q_Y+%V5gK#48%hniiNXGSW zjES0(aXn2)kb5(Mu-;~NtdRbZL&4inDmA3B~e&#Pb&XrC^F4>ho;oRnK2+`4@ICqDV)jEv4Sz)9X3&-FN zr~l<}>J5tE{)Pya))Nn9tnjjzk+N4OJp6bRTf0Xy==*3cmyaQk6+`8_G31TPbyE{#rkXfhD1jA{ zan+KHtJ(lFqt}@+)K0`XB9VKK6U8gDNS9P*c=IG`Oh}^fvn2Wc(jVQJOu5QdG$~dp zpS1F(we0OlS7F`_8-==xz8)$2T&E@5tnMIP>7dOlhs;Yk-_@MbwQ+Jo_5$kEb+K`h zOMZTC`lY&MXU9#KJ|3LQJXEceB0H!l(r-`UN+Dq&|D7g&;8&1d?y$}nMWbCP=?a;Kq&n?gwbh!7%fVL^SLjazpjU~ zX=nt2EfKV7EWB)Dr0ms2QrZ$F-#dyfUqy3zWwhx27%Jt(u=`;Q*>MVn?pM&MeXMw& zu?#L5#~814Qm)658KRV4uJnbQsVF&3Mccw^Mmp7$zN}``w{llY*KBkX$*reqc~v-` zXm>nEu8Kb@+4bFRI?1Z^vI8W2v(F54@f-L~GOj0)vd6W@$kVna(F>Aul}g~1KY=l~ z68K&+uHM_s(v=XuO!CZ?Zxi{)WZ|o07Rq)^qQk5t*8MBHT-sz#>`W$1y5Xi=+4s0= z#na11o0&E&k8I*8+c_avUp;DPdKu{|eCuG+2nR0?IE3Rlxft(c%W?5k2D!-H?Gi7? zO_d3fgWQ*WpP?Q;Z1RW)n8N#|l7TcxW&Xlc^mV-KpY3I2WgkC}_u(n-=UJwo)Te&h zCrWO8DUCBB0bcA5i0_&%+FiOw)iWqHCWD{;&7h`3c-hrVhK6NPb$1r)+Gk6imCd&m zatO)Gq5IPuKIlK;Jg-IwM~8=S=V=HpwV@O{5-J_SFrF?BO2#$*k%6KbBVk93n7f%^nTZpX5-5?Ez{N+RyHvu# z51LurA(2D#66ySz?0ZO_Isbx%&-)|^2apcwXUPneu9^C{WFH-@WJ@>v(Q~WpT-x|+ zg^l8GZCvSRmtKJ#%^^GKB}7{fbTDYNLwrAn^hTW2p5-K^nv2gT%1*!=7coO6$6nx; zx#ecM(Ie;I!=lI(p2@CHN<=Dol6y=VA|28LUViB95nK>1EVRkAzFHqF^nA(gYjhykw%653>2jIh-L4xSQ=H5o`Up+ z)eqxntCYL?pb~9I6`K~R=vz+qt}~@uC|$F|u^Kria$k4S^2-vfJfHDA%@w`%cRb^@ z(la}zL*G+R^eR1FYZ$0H%E0ZH2F96<)RG;0@7E^gtToB^`m~>kq<|dQm$klcar&ju$*V<6>ckMBR<)Nb-V0*d}b$e zh@FPxWe5ADoebFr*kqA?t}PD9*PRRhb|f=LPuY zX#j;mboHrpI`qwu{*8E}^=0pRQYKraXI8;@c$b^nk5w5nNZn}3Zu&pVN|UaF1$0GC4YybQHka{7{R6vk;KoB0eo7Hf*r}Pw7X=y0gb>UI*JbWpbn{|{qtCQZ6 z?DcHWlikq3+{p&k7BaFo5|Ue2wQ-Mv~qz8mM`!8;%8?0yb;!NI*5jO!r#J&Q8=AYI{A+2W5r%;LFp6mkw{Q>fb~oY&qO!XIHn<^2;w zB{L79!?X~|RzoGj45iwIQ0bq9adk}?D{6&H-yxj#&%>q97QxIT5%lj8$=$_~OfMV7 z*XdCdx+9t^JeqOaqp>xQK|f9KE~H?rUBQ_P!p)>Fyj%Lh33Y{=jTNk=L#i?>S$10b z3J5p*NhRM^EqgC&?!HpfD^dLL-!<}H1g&^pvZvV~p0<#3r*+Jz-^VOui8n~A2LWM%pbE9P+N6fU&MdlaNg+CloGd9vdp?^HM- z?^JlS&A}DPGdE6^ezSaX5mz~bkN<>JCgTk?C@aikV4__!~?IA zO1%;CUIpP~Ma*7?o$)e5_OFMm@o~DgpZ;U~(m|5D+MLFQ(`mAY9gu!hfUUyI#*R+M z`6?ZYIRle$a9v;70a%kMy~Zpu#%3}0O%}5)+0s)JtmT(3YX}qW4&mqMp`6)0l-sT4 zuAUXb&0?YAMTN5Sawzi$h2h;8Mvwa86dEUf_^WVSrU>evh~RV&;bqH2V^@sgS#}g# z?nN;oCR+BUqcOILk)25CAQe{#w^Yzo`ofzA$D-OCON+*FG>{HyzYlRtvMDKdQOULe zDrRg{Nv}c8o=LK^_fgGv4vp+vYiKXob@N}evJ)52oLTX*UoQKZl5sWoOUEUJo<#@s z_`4X0S!RH$qN_(3Dg4sN3EAx#f7XOa{L$_|OE0{!nXjgaZ&oBxG6UJUzLH49V2hkd z3(p!SNp~`d&xbPkQ0I+$0qDa9r|^XX~VUBpuSmnJIjCD@C5uRAz5Z z<;#ZBGaKg>-%`3}Mj!XZ1MewX>)>~OMoBL`kd?*{chkhL3dp-U#2;-dKKNY8F-vAJ zLOA%2o5Ih+Go}BS$#u!M?6b0XRa`n|qOJeKxm1-Qbjllo>Cq6fl|z|%U?}rCgh;sE%8ppdNQboX%sA}Dlw6fw;mB*U-xH!D zVw;L!3$^%eYLZH5q>HJM&aQ?EQChJl`yi=b!RAjs|t?dLp_@FFOFI^x}O< zmbuO#p0$y=Q)IuVsN`FzCLFg-^ofz3>w^i5?P3HiN)q|)<5*_pNQ#c2y4 zrBiriW0LF~B=h^MWH!it{W-5yG?M7+E;f!!e>CjAjsBr_=I^k}evJdg0tdS!$1ER^ zTw{fXU+bcobj)T-ULCPh^i?tO$Aouxn(q;v>yhU}_}ID>ip!b4k)FznYqFO+ zz$?3@Uba^CN!Bc#!n?AU8!6{`yPp|N(+HcKMzVM)=cQvd?2Noeqfa_xen{t!YSJ$o zmLWT_8Pb!=B>hMx{;twP`A*KXXzTxQ{_&^5T&ppJ;v z5W1HOrFKRrukM6$BqB_7op|9b!l^zzoJ)lwm}ZaQ+Y1qV)<06Z_>oMl6@_)A_@mFG z7@&(L>}WI@U1RupX$;lME0~<6pv67$#A9UNV{feF*|N(uU-IkHO1=t+w*FH|p=gz4 zAu2|-S4-zm?(4GBGt1I&@qvaWYOVCj1aEmK&XXm%M5R6WsC6dNiVI`iM6qBtBvvv?Xm|Yc}Ee4Jfq@)Uw5!L#3>!Z zPyPS?V_iIwU5^l%+so!jc6rAwof41iD9R4NF%Rht#ShJv_bNzUeY1yPEPcZ0LSFF@ zy*Q4@xt8we#hE_Fe)MtO;%D04B_;cA<};x!Uy9}J{}(`UU&#kmW6P+ zVyJZ4Lz!_u6tf~sFc96&MGd?h|g(Gm%)2_#7OBZ2x_GyfhJ4&FCW_5l-l+f?TFT=|F(ZXKS) zlQ&81w99+0E{k6_$jZjG^3IHEHp!iA%#nM#hV%*bzuJZCIAkZ>L4k1Y>p^+<(F3Pw zdKWgyI38AXQzFqV?{IQ6^}Wn->4d+X<)O?&kLb}9R?SPH-J2AyN>4a+i}13VULv!- z@$e*<1?HfbC)zVd{Dco$NLiYLUCH00Y>HS>EWP`j*;iO8u zcF|ftsd-*oL)xuHbUj;~p;4#d3jLXQIGOr5jC=LkEPAi9_0&8(5AIs;jmlAnTRCcj-?5zed(`WB zdT4&0+zs515KwAONZ*!rpUq<-)r^qb{W=<9tIn%29&){^z8(Qtr93>ttZHsF8zd0M$yyMH(9 zrqe>++M*uMTjiW+ReTj2{RlQC9aw-QTTO&$!jlO>2P4ZTY1%Rj8~8Ui*A-h z?s`8SS6`o2Eca(XuKbKX%DwvL=4}1hmA;I5Ib<<%wQMjQDVOL-iKp}Edok1Wpr&ZfA1-W zvF@k(v-@eY5{o9TSxtlnfn+qJ6N zhgQ|BVbkM9HnQ6`oio|BeU+VixkLR&J7mAfy(`_R6MLLm(vv&ae3yDPaO>f>-DC*a z-(HV4Z1HGg$4DKW7O5v6iNf28(iV7Ovn^Uz)uEQ~Kj6Pjj z?bF|F{B(Z#wFD+5BN!|HF&I*N99c@T*DrG4njWvZ^rS4p-n|_E~J>>0lE}^}1ul8CET?#kaeZIsAkT zZQ2h1X4gyY&_~BR^z$8uW@kGU+|Qj$E>&OR(!R!=(?i|lP2GAw&ZDv4dC+DfxqC*c zac$Uec;Tz3*xT-CMQ@D8@9?VcR5Y_1F**iUSav*yjuIbSgHOTsew~==*9VWpssc`7 z$!WO4X!Ox<;`CQnI&tRFPuVa*XWvawubT-vLKgSHa}Too297w zIJD6kseE(rt^Z2ZA?{lnkED_RUCTIrimiW-oGEpTL9H3@9}>pds4>YSJc% zPfj1@9=_r5Co9<3zV8dkbr!*n^Bv#GS=|;d>uqv9%w<)&)0;7w4y1od^piKF8EZmt zMWym&%&El{o?~<(gYD$>^g$Fo$=K)SkJKbU8%P83`YLK1wIO z_*+|}@w2_UGQ%qaT;VtU9=URjiRscp)`F` zo~FK+(9Fo|`F~B1K#*ZA-P{l1_;1G#Ip$w52$hWi$L@S0FiVd$YPo~Gjp zKP#SnJs(SshBJsBHt{X+^>mg`nm!ml6qG6+>>0$F~>L*^f z@fdws$!zv#jF!jxbZ0x~HB3szr+)M!7_(uqnsJRaPmj}o`{T4Sf?l&PV1}QEQ}__h ztTq|kO0>}{i7G)KP5p&CHk#SmE#x!1<8je_G`Vdmop0p0>f>3DOj9ym=I!!y?Z2GP z9VSDyj%SdO&m><#kAfkKvx1JRPB3Q7Wxq9sEhE2Gd^1PCX6CByetb-k-&&N%`JJzB z1L++-nJsJUC|DQsG($H{Nq~#1Y)JPLt6LFX28~q88QB%FxSx{th;QuB$3zS?{a(nEo>9>oT2hN8Wl&xvpF* zCm)8t^(~&q%M}VuMIWvDV4Wuv8}a>&x@$3^bDPk)%_CMYFwTQFBfZ4UBGn~T54uu|d(xJs%c9OkqXEs~w(%yD%?HcdaZ+G4B z1|D_&*`t2&!iSf^3pb0>U&EvH#kD9hDA8KJ8xEmF z?VIzucU(ZHR|ORNcn}sM$hj2M&<2HiIj>Otw{qv=y#Ao92)?pNvtKIKXQRk>oGr$) zFVT_tB|1E^d(+Vf($UYJ;}M(8McX2{Ii-#wWO*?M*M|xuwO~je&El-yAFYvag*BHmn>CN&-!(=DZ^E0=S9tw^PZz`e>ISFq zO`BLPo)D|6)#$DZ@k4jU9dEpLk-L5Z-(xeGy_2ozrJRUHiay%*tt9em_}5wRC@_T! zqEa*y_UJ%gG?aCzy2gF$doo-H=pX%??z;N*GhpWENO>reSu0an@Cg^wvlRVHmNMy? zeUAKAZbveiGjhm!;%k+`jULa{JaZl%59<(D5|Ggr0eGRNaJG|=;n8oi! zsVtxFqhF%+w9%_gtGt@|3^{;tF{*Vh2A$id;lKLmr-D=X%Fp+UzVMG?>0pi1qJezF zFUB!*$Ln9RxO2i2^wpR6nH>^gUK4esL6ZI)k)*NzC22*7%=Lw2wM|V?>@O)g=t@=P zb}|4;(;I8jWPT-GeHP$vwaU=PQ!=A>zHZRH#O9&7y)d7hO6BX-2l4E(^Pya^O8WRI77(%TaV6i%TU{+$^GCIPI(k?N0MzMliNN@pH7L=it5oiTokP? zzeQ`K&8tDMW`mx`ph2;S>@X_ zt?89c-i3T-2Xyc`8Jg836UHc0A63oL4f>;NpM^0at9Kb6b85F7eZgA4Q9D6snP zH{8-UPYIvGj@E#y7k#Yh7Y1f#)bu<9Z=gghCo|C?q3%!ooSF`nYBpx=a~o zV;T7cn4@{+s)|1PLD_?K-nu_Rw{g5WE;Y#C)Tp=LF>1~gqt+*w$R?R|q??&MkePli zi+-rEXv1lXmb$H)xWP(q16uebn_hj$4#Q$+=C*5ZxI@cVI5=;d@KtaM*PQB??9w|w zxX2m0wQsJQzo$oq@WLn0F^hZA%(h1A%F9tY4numcRy1rKbNH!fWq7!9dgLiMoviiUJS|Qj z%d;L9yeZwIWd(e5@T^~Yz$NtgJ)rpBL4Expos`#u+U+jXnq};7bhFX9blB~H1!oRl zII>s=_=c}zts|Gi9n}x%{j!jrI}%cyO7+aAr5ae9Ug455`lVro`||x>RIcmwAFQ*h zAVMR4iO`o32A%xOpeikmdS;YSz5X%cmz$Kd+oa^z&CI)ICDunDt+c53d5fO)TJ`uA zt1iA|)7}qlT2%vmw1_$Us9hejgPGf*Caqz^$2ehp;lneS!+&%s`3=15g>Ic~=#hJf zM@KI)hsQ;#=k`cF)h$X-&w(MW8%<`4^*&GL!{^0g@+!J(4E!&%Sv{W)4f5eP`1E`- zT;X27&f{JG@>#5U)2C2vQk-5Vi+ibpF6oPOnJ z$DPYkd`z|uZbb*@yl#uXb>|-6apv$AJLyR2m`4UXPy7GMBag**d>I@HY}wH40{R>Z z^g%0hlz{=gaws5mr1Tl zJ~Elsez4$2LwYQ%RACEC^Pf7ogpE+%`S`s9iP7 z?7DEmu03|VkG1UWc22Dr@6_76PHoG94L{&g^S*B0Q*dkZQyzMzJlcQNL*Hklth*u= z_BvkX{3!ZSqREbuv%efopO06Ecfy;&A6;93@6p)DpY!3l@RRA6=b&GWP4p?Ojm2Av z;n$vc!Pi{$>6@aD*z)-&U!e! zC~pDm4=aS3Dpbd+MRLHF)!T~4g_iR6yT$l@#pwF%Yqtb-<{g_&U% z%-o4A{C8#sziLt4B&%+J50mn$O&jNMR@bqU54EEy!wY*HGH-OyC+d`CvQxgQXea@f zW*>Ixgu$&&|IZ6QJKCcyH$B`DBb5b1n%^f%GnYiE3S8m1;n8$v@RLS2<&Sup^cCJ- z7NcijQWlJ%bC#~FSp|M|Ip$XfXRL;A!oTj!IXxF1g?y%KR6H5yz7}e{>k^|J2U; zsyUwb9xml;gslJ_29A(B_D*!tx3|zE-HINC%Ah_w8q~dBg=9Dj)dioccOuNu*F~z0 zpSce$7EGtE$C{xu+`s+@Y{ck`2^lG^7&I#9x zbt4oxAVL*?L})^!K})|jXlne?wRfeN^U9*n-m_@x4GaCZ zR;BN;%G%SWm*?AHmF;ktcAY;Dm*REkw=E7l7$<$`PThW(Tu+Hh{@+|$Xhjc)7k;a) zN43Y3)4WSgGdEI$U`Qv0qnCY#fBkf{ZjO%D!GEK5i*Cx;gZP>TJdf4z!Z0akCiwXK zp@oOyJuh8X-^Frn$LYpmdU0CD>-_k5dPx(slFsl>XyFg1k?T2>sHL#NDO>r5 z_n?zY!M&aCxY?CbMz5D%f@E%cP4k8 z`=IwIlKXZ0h1pv4M2-%WF`pgDfx*bt-PyT%zgnJNOovJNhI@9Kd^H=APcAxNJ9@(% zO)ns8TY%;qP}?;DJq0tGUBA{&&sKCuaS!0Xz#rh)jxnjCbieYIU%;B3@ zYj-1hG6uq=oDAuTnZ6_5@2S<4_j?9CSZB+gFje#n*R1*B+Swq2ewqlp(FnTc4C?cZ zLH1XS@=iA@uBJ)xA(Ntxo9NLu>&-Q0GW`~{7-!*Zv*4*%>5;JVJ8#p1#Wsy?Y}eaE znZqyJRT}4@x7?w!E>7kkCz&vpj+MLA;*^WGLfrawy&HDHquY}_T2nPr6M~U=zLEOc z9HquzqNB7W-}3?eM^*4M176iS>gBDL7@b=eBQH!!gK3=Aa0+w!lh3?>cO6f^Fd6&% z6UzO!dp5O;R|DRC*j_thxE&L8#_-j!~`e+9lX&-uzRwb)Jd){V3ADxIk zdcQ;Z_q(aw0bo&B>l-oYy39AccFzp*-0=TT!T-}|=DnV&$5OIXu``Rl_H29(-e`h7 z8ke0zM?sFBeuXz&KFHP1tGV}1^2Fju>ZL*DZX)VqnF8eIOrB$MjuWc z-U0{&@yCPw`9gg=z7WnBKE+$ac~}GoR;*!zi@7fp%g&wa25bFz;}D(+nVwT2-L;l# z{>oCgj|c0#H7iUP_l2ps!s%}fS3;8r4Hy=o5C4iFuVm2X?FN0{)u?^5xOdeyX;(iJ zZ&#VL(Ph@M^=5tiyhZO%w5YU-m8`QBw!})$qfPBr*wnDOoqi3w9==8vAlU&&;ZVKq z=;3qGN1t$^!L#0HUGO|^t%euw)e$dqnnx=iiPT3hqzg{ab>+a{T2F_@3(RIyqP4z; zmv^V&k4}0u+#REJn_|@K6}qlw`QQWm+zshUx#Cx|v{>!k9}9~^-eWC2c6eNqX2z>) z13IpTC&0BP=uR+Et4<}78>2IPcM@6CWR0LVBlD#c)t{510dQu;XyMmb>p$4nefWm| z%z3?p*=%Yz{H?henpi(mgWe{qcOer-CQAul)1&YlZ!?kgn0GCk|8KmH4LK@mm8%bj z@E+HxT>TxMr;K@dDy@;Phzz*VZ}K&*ZGlpUvEHW&)Gs`svvUKQUn8jRGK2K)1^In0 z)b7!ADO@hpvv6kjxqn-;V$~f~EFT(aGuFEIXJ}`v^|{KBoM%FGJkUwG8l9AXdR7Tj zW`3B)9SYNz;o7A1YMN|y}paR-NPu&V-BxpQjfPxyj^9|v(aWU zEM~gUEHDxl*g>my7h2I&;1oGi@c85Oa_4yZ`Ct^9Bxp6+p0!mH zbq3#~?>Vy9aN%$KL|>sHS&x2|3_G2or|43c%36QMzV3M^RWAq9R7?){CUf|_uhTW4 zBi+%{Gi0chsdlAw*OAv-WX^(3rPr)!wq6;Ot(m9ky0YZx=cPGX_XOuPxtA<-jau5!q)!K%H0`2EANb80{;gT%oh-_l zVL?;joG!Jh&k6WKyG^gGwc%sh>8RxAt{v{dp*IgWl>MetQNAEipT!Y*Bz!{%$*p===ZZ zRLIHHr~5KBq!;?=Trym+!ckq{@*F)1WVjB%AXSC` zZ@&ay`0)bm4iqSVZ-GQkFAhCr&p!eBMuXakCZ61?kbPgM^=AuJZ0Ajv&x_QdNwHpJ ztv_O|b6M*-pVLhVe{^neh^%1ffBX6APJI?w!ur%eEmW+GZ@nH zp^=*NSEMc_M5z?6uw##CJu@Hw8rE#*P_Ksl%@6$1{a&*Q@gBemxRh4tDC3#Ws^VRj z@CF;+b8e;x_YvLjZkUwBs~LJLf!>U5%D40>e9gzQRkB%}1 zzVNSrIz|Q6Y7Ojv^Fq}fT&UV-3Kg{S=EG-j;f;%R8m{m#YyBK^`1+M4{J!&k&(H_! z{O*o7wX|Db?VR10Oh}j>>mR0eXTs>D4(Co3&bbhw?8ytcvJGQ|{=bKR# zEluh_%7or%QYaNJWw%-BuUp{uEQ+pgm7@{{>AY1@UK@R0Fi0KkbW}2zJ?wymaPU5$ z15E&4_)GNgHt@n@UAlJLrPOS<&h5wlc+(^6GLO2zkoFouw&!}JK2MEOkG<$-Z$vXs zN2@)o+4K?UC^x)vXU1sGff#j*@bPwnPmXqe-oNl``NQZa<*`aW7pqQjaWd?p$F2{3 z;h)Eo$>80!DG77~CF;eu;KKh()Lzy)$XcH>CadncWc}J91(Lkf7TtA!@E(AR_cO_R^lhCbjmXl#OIfNHl??};EmsTll)*WA zmu}%1Mm&&3bQM0HrzNoAoyga3>6)*?DZIb-Z$6wN`K?U_THGce|8Vv;eBmr7{iFCF o?eR79DhpNf6uP*DUdk0@GFj`k?CVjS*WK+UdTDhDtV&4#2V_-SLjV8( literal 0 HcmV?d00001 diff --git a/modules/plugin/epsg-hsql/src/test/resources/org/geotools/referencing/factory/gridshift/stgeorge.los b/modules/plugin/epsg-hsql/src/test/resources/org/geotools/referencing/factory/gridshift/stgeorge.los new file mode 100644 index 0000000000000000000000000000000000000000..4e3337fae2072deab703e779dfe2f87e2321c049 GIT binary patch literal 30256 zcmeFZRg}|b8}%KRe>|DN-L1G6$;OKnFYfLR#odd0ad(H}Uc9&!EAFL0ahZ94efZA5 z)A!(=wbGen($y~ZeQmq$P+6r){#ZFx>Do1_sbr1PrBYQZUG@iF`2YK=tWssGl}eTB zXOqcPl!uvz$&~fbhzX;Lm`sgJeEa^t7ysWs{xA6x$p52LeZO~VcSff^M4ZZPa;n{P zhn8G*==NcUl5TXU)?A0Cj^=p}hq5#&sZ&$_Vc5T0ESC-%H zTDH-ywDas*JjSjJz3f`w%&sC;>^fG&t{R!`x)WV@sPhRkJCkxJ_NN+H^LdP3dhm^?Yg7u^U#UIBr$XtyZ00XjO*sR*mgr z)yI}rwXJT|BoYSNr5CLKF$(zT5y-I!<6k+CMt?QK%y7A6I&nY6i>NhPwH zv@fAaS!^Z+?!-se?TC*SnH(S8*fu^IDH$Jalr%m%=HIyJhLds8GmGP*fA)=w-mD%M zJ)1Qyy4M;PU48v)bn51>(LQ6pMr$zC+% zM{VX}r>;+MDp^0LYPEK1at)`BM4gJy=2X$dPW^0mYV#|HUf*=6=n02LZFA`IB8Rd~ zaAC82t&cjCT4Gn?Np>yhZ&#+a zc5SF-SMHK_ZOUO+_GES~aoLsRtxaR^*!1QPo9geh>9=LzJjJGV18v&d&ZbapUN2?S z+MG7UC9|oT+osiTt#aM9s?%w!Ztb+H#&WB!PqC`wAge+htlC${sv%{pYM#fcW~sp6 zXVt7v7X9(aqR1tS+8(s%(t3+(%(dwD7>oM!wkShOi>_9;Xu2$FlEb1r$t`kwEV}c- ztSt}C8g$XD@(0awZ!l~3T(g>tHS29}v%0l3>t1!UO317++4)^kvo5>Ly8hOrLw8M@ ze#WHQyG=5$G-=i}lTr;aX-Y?vUe-0KXjzj6=P_wdN|T;>ObWe>kA8nFK3ZWxe6&*U z_-KhL@zE@q;-djmeDv;>xahX^anX?@#YIB};-csManVr^zeWr0{Th8a?Q3*( zhp*8NrN2h=Cja_BoZr27s?dF>dY^S_({87pu5c>#G^c70c4}$|r_R)IDpbaeo zlftRThEsXpJ2c`R^L@slV!IuhzrrEUG>3W*cIa^jhnm-Q=tdcbYUgq2T1tl+dL4T3 z!LC*h?0SCAuFiYydbP@~_A~5yIMlAXo$WeZ->$+x*tIC1UA8oKHTT=K^OH@%M__xw zrWN~`_tiE25xo%BQhux8J6MpRGFl$f`yctuhYq^IEH> z&9-XvNULUcw`xxltKL_(s!S29mSu()5?ZyyW>ryU{s}X`><^2&?6RozN{iCYu*g5$ zB2zbu-Zr-AMP<%H5sOY_v1m;qi-tKYs{V$#ziZZpGiDXtW7f`9W~H5JR?p#PZSHE; z*+yoauV~i3LSUc4ta1^v{*E`P-*b~bUp1-G?K111zr;sJT#S#FUK1bvcW8WcPQCbOzryk+o^RU zoO;{UsWOe2`3g?mFW^-1bWY6+Iu-ZDq1I0vx^l^(G6x;nx6Yx0zdE#Kq(iy8JG8m6 zL-{K?w5^asg);Cx4FALRE1ogKaQ(=`;J(3*KG?Nwv|ZtzcC~KG&y~Uad%G%Uv}=}m1)@S{yDzqhG!CiEcIrXyyX7QL`) z#dWKWAG6A|&8o_at=c%zsvHBXI@r#tcFcUrvR3__*Q%3gtlAv3YH1u;KDX$%YZje4 zX3^bk7JXS_kq^F)9Sr`Rz`ub-56WA#sh~w2GVr??i)NZFGX62E^$oL@9yjacHn3i7 z)-L#dcz^i4ty$M0@+CSAF0(m$t6`m&ufu*9Sn6HU6>&!pd4 zn>4kCNv)zL<;iB!C%FEg6|Vm)K00q(d^FSe_~_CW@zK0R;ry8RXo2T((LINl^2vACqM(NPEG&Eq2LXN#vXIXy%oM+=#XPPeBQ^Q z4=o+)SRKwU?oh>S4sA%pd+iQ&eP!2)n|7rp6H2TwmC* zY4&@o5ttDcOu>QOJN-n6jFTf?ecC9G+(VP3?qh)5qM|X6LkLE5DAKjKb{(m?p zf9q1+J1$NB!=+O@T(T^4sqkc%x({$^LtB@g)^aIBNtasYaA{pqmp(gPs_@3CCAY!& zlvDL~IJIx7Q<)|^HG6uIStOjJDcjsi53oT#CIheuxkFcr^e80V^RWViZba=ew-&u7t zsa0RyRs}y;YLA`w5d%x<2A|jDn8l;jW_Yt>j$^C=5y;+YPT}^-D>#B zrCASMI(yzF_da-hwM)ZhxO8}^OPH37acMkaH!f^JpL?) zo)32@XIF4;=+KJt4qeahkSnc21p*GW`;2}+vg_IfyAtdN`!#m0oN1S9xLrNE*!8l3 zU46>im8yVUx6|6SE@;>6uQn}yYSZZ}Hl_F-pSQ`T)cE>y`1+ZBZ5r6dranK}G`=+a zpU0*TsnPl%UO&#NG=HP%@O|=g2x+j z->jBr&5G_e%d-N1KLuEk2A5oOU1~HMjC;WEOh)2V5{J7wGC zRL^-%Jsjgy&0b`D?Kd-tMuXw-X<^gH>NZV*@AKud>1Ils zM&R)>e6?!FGpnN4z!j~x{$`byTGa}F-(|2>bvvQ+4Xip;p7#{AYI=HnU>MxvnP;@V z!=Dy?g6}(SvgqJEi{hC1)XaRQX5_seEqY&=*+=Imgvov6x{puHnsW(`hwCq`F{?TL z?&DCiMt3qRQjg!2F)L$kykBzszuT-UXnY$q-h#$2_|2q3i<$TF-*7JU+O4j)+*)$N ztv|Q96&DrQJ7ON_xVj-rLagMJ|<_ z;L@hPE@f)%(##qz8FH!bcP>3k>{4xqOTWK%D(7EL%{mF+|AwDm!kkaU!}W8jY#XOW z)pY8(oHA!8!zY339dN{JhgRKk=nwMTbMl;ju|s($;NSZ?G?_ei3XM-dp6iUpUqj>T zqwz0Z*)@ke*NQw>cDr5WmfF>3vRzvT+Lf+@UAybr)w3MDpWm)h>F|DGyC#_Ix=WU? zj?Vve(xxF~xJ{lu&H5l_yFwB3C_=HQvmD-`pEj9!S*S5U&XT@AgcmFpm3U#+-aovgtTXw{YLwcn z*Iu&*elRJ~{cku=FJtJB+=i?v4dwM1YWv=;h4*)?ho z{Lah|?QPfQmUdZcGWTTpP;R?^Pi5CH0lPZHkUFZW7DAzc>0IrwF_2Vz~`MO+bt&Bl_lHl>tl%&JVD`;WKZ%r=x|q@jl042^GO=wJmy?+O^onGVc@ zhSq#>>-l51@?K)554d$^ty^hlgYgKrE_G$r8@V;Ff?LjlZgo%R*4?06KYnH2pTPB( zU7C8(rMPu2HT#vB9*LHBcd2J%myVI={DoZlF}+LU$aDM2b5F^0k;}~UL8r=-=bDq} zx{~Jxbp`K+oCP#~HyVFEjZ;=MJ`WoI^CP;s3*@x@4lP~d(D!Kk3p9Qo8o!~DL;EXo zHVQdZ2d;k-V=ehl7zMhUw%=EB#5 z{bR8I1om-Y{{-w0g8cxnPXYG7g8kFsc)#v8nVZ4^KiX8YC_X=nO@AiDvsFoAt(s%9$`98M#oPZwms4|{MZ0EM zlx#RycP7W9@jsyPTXI{JJcUK0-Q+&9e7oCpH7Ct#^&9%V$gEy?yZ>;m`H!J-*A4AH zYUssgLn#*+sy2?-dK)^{+>o=Xp|V8{{`-b*#2U(GF*M`_^L@>&+(+G-w%M&O^WADS z&aI2R-BJs;R#$Z^K{3un7PlTIaI3P_txf;9qN;LyTf(090gb!&%u)I3m_ksP#6L9@@ zGCMt9^{I4qgYkbI$#nH>>RpaAlGmoX_`J#R{eq8Hjd{R4pR>xg$Evw2tjatUejkY6 zBh!z8>tD#K@_75ji8upRy1subT11}zh{jjmXweFMy?>NNZMvcB4K2!Bo-9Xq8%%A{ z7`&YaZ#VMZH=MVX^yq00kCG(wsJzRg0dEYgyKU&C;>XmV0R|2hrTegl`^b}QE(Zq3=jY%g=G%VaoxfLle{;pKkfER=MsdQN;m zGPhp1+{#0q>qwrPPo6tTo_o2Jb3mR;MxM(|p36p_`>upbIkSU%5_BJp&+?zf-*PH? z!l^pjoa(;_yvO77`#AMeD<^u*nV_?Ee8(9`?9^p6{>{J4_#KBflIQBu-4$5nQ0bWt z^%~*O1v;GC_&aAcI{GM&9B@B;U&K$}2lh|EzTOSH&ZG4enEM^8>`FGA)Mb|o)Oj|29%Qj_zT`yF5L?9XhPL!Yw~kG}wqUw%GwJ{IopjkiPVvs9tSDMIce z+f|6zbRK_S7M)*ziJX1_EnfrX?mGvY+|YuAk8RHDG@Py`Rh6=YsPGan2hJ|Z&Met-IE$Nl|$!eWTn$hU{gZNfBwL#3HZCChphTcms1dp z?=ak|MV+m>Q_rd_baPnR}57?Y-qy1^3E^l6?=qXEby-VyKvjGx8kmzUo%C-`#4xkr^k?of_@d>mF_yEVtVYj?-E`31b?_YB1^a1pKt&4o^(!}8|jp~A?H=wWm0bS4M(#O=~d_Q{d z$*IDR=-@Aa?SB0HTBpvEReDYimy1T{2zxn!6bU4L5PG$Z?W}~+qz>e<0 zQMx*?Pl3mKG@V?>+^_BF&@U|=Y6|D4EaTAjd~`YK$$sd4A1j=|&gRQ)yPBLq_xI4{ zt){P?1>d9hO}p8Zr?Fi}E0XIA;qlYkRRF$^`(o4l$LRY-n?~)k>Gmp{GSTNWhU+J^ zC(HlDYvejtcAK&$wy6dhKky&&I+^a$?^ZcBT2*Wg86IzUlq~OSXjPLR=HHCkT-u8-pPzBAM= zv7se)LziE{PZB)%Qa599Mb{Yqaymf1(| z=eKp}G#>v@IXpoDhrWP)0kBVt*E{l$U8!%fr$z5yY$MApw(GYEb{*`C*Msj%vLwYs~>yX9VR#euD4#H$2n+I?zPPKOq-g}<8?#h$JDWDX(_k>U;iM9O@@Q)hqvo; z!>Ub3t@?Kp*=`=YxY7Jfe^(Pvzo!D4Uw}@JuY%LOYC71f$sN4fTgR)XrM*g)+pB6R zyc%bCb?O~h-eacE!0Wp_`mh|lr+Tzx5cA)`qt3ewg zvxdg(#=ot=zfCjLZ7|y3!O*=rhWzBYY~;BTDGXIKc)UZ)?~>h4yXDylzb|!b?$pES8 zcSBB%GjkT$(f*9yXF7+k+waiV_256>p?~mt!}~dO73}}4&3=}-f0W0er+B@^K8He| z>?-)suFU7{y3fpaz~f(<${uHsUFF-s0qA^xq5Ik3goJ#z;Q8_Pzo79G4%swqolVPT zv8x@%F1{0ePF=d5GB%Z=#~GE(rejWcQA=I8T4#)&5K1 zd%&xhbzYU8?N#p)Uajlu)#HXjGYFO6!#==v9p`9xZ?Dk>!#{jSqNqXst(? zW_vVo1ZSYDNA(+dw7-H!=?Z!@ES*QUf*xi2YN*!}JpCm@uMZeXvCiPXW2o*3L#?`! z-x?ZfT;5QH{D!ioF=X-^y7&=K_rR?t=iExS$E_#i_$}o4e&qOy?cK^;8&6+~2i`6V z-mVVbZrK}V`VPJ9AI$hpmwsAKPd^o(H`t|4o$&8)eWMEezK~1lGr4pJUw_|@CwNO= zPnNF__LUCP%WuKoG4~0`^J9iPWy0TgXzA4NwVZli2A{{S?qLR}W+dP}4yPK@-8a5Z zXLEr&7zenM0rnNSo3VThIj)yOrJ6f5ql!Zdi#XIJ6FXdHzBBu}ZBOZNFR`oJZ`VI? z{iEr2{W;jKRkvBfn2} zXvzSGEX;hVnhy4iboN>Cd)zW6y=zr*JuGw@8PoH2Td-^#xtr|&&?`G3a4Q<--1KmwNd_743c)NOdyCm^$-G54+yX@AUgKnK( z=a%VLw`z>!=WghEW4EeQa;rYxZb?RZ9yH!fj!*K6&gL&LX6EnH-6dH=PfLeWcdSeG zz=kx5P!fcH;! z%qsPJHrI17Szw00kV7{Jblw7%)Oo0U)q%U27c|RO*hDLh3B!u z8T}3Cy$^l5f8Hm@KA*C#_Nn#^pT-UK>0l?HKGyRoZ#kcO=JjcBDxds5pBjJkYUcy5 z5}x;}*IuvwT;(Nac{OjSSB}nbe0{Hum-8w^KF&aD^xX&cpFFDe(4!IO!FR7mXIFXj za5}kth)1_OdUUWZ+Fr&(Mu6kd_{ z&+Ko;=O4c8R8#gg!!|pW9_%ZT?LYTH?_05}sqWOCVot?maVl3VeXog5@0mkKFFDj7 zozJ|Qjvm}6gL^q}cZ2&e=6xc#w*>b};Ql?h7XkN@2ie7t=Q_>A(+$PTceLwk9bPBL zPswT5-%0JN{E>mK24nC(}ht! zrRd>P`zAi^ujEsL!f<>BpU#Il2k~Cbc;?l|E6n*}w0?tEspfk1^JuUB=;2k;ro3L+ ztJ#IUI+4Mv*I}=`@g606MmKZCBj+J-T?hWN$?#}=Ry6)v1AhJi|3{9WoeHh@lKt>@ zIq`P8P8%w})6fI@n~C%{<>_w{w=ra^X((2PO5*8%#nY#9;OSp8?|-@V_LN&4@OBx| z_|((kd2)PKWL+&bAo1f5$#ADSp81Qmyy!JDUIGEc-h8`$OEv$}tr!7zqAt zok~@cKA(Ko7u*wp`}%l?syxNpGw%zS_gc(*0Nl?Hq|b%t$JJ!U!M!KA4*~b7;JzK) zZ-aZ_5O;3q?Wd6C9uEWW&fLMQ$G*0#U0vwx|Ksf;6Z~4%*RSKP{CZp6uVhg_cRc(W zkjSr1HosoH1k0N~^*Zj;?^}IJy2z()<9#~U$EWPAd>UWfr`N@Os+P^C?TLI!Z1<`6 zE3Z!91n(1GwciHcFY@Z)c&`%lW%gTnRk1pLzBpW;)vFu{;Q^~x$I$pe*O=!c9-Sh` z*C)r{9pzEe?%+?` zQ)xp#Ws!8h3kJN#}C*-U%wa~pW;>tbiM=FKWxU_*KlilX}4Sj z+?tWet)G%`CcJJXjC1KB{p}#U-X;2+vxm^~E%fzE*wy3r60@UUNVcEf$fX38nEk?Z zei>Y99p?T<96cU-UlOlp>~m_{N<2Nd-x$DoXyeo~=Dq)a^Nz-Ej=&3X=rx)CCf;t- zUWcZw1na3hz_K@i+Ht{o$gk_|{R-Cs^U{7z$mQ4h`4rp1r|xxpI$PSObh+XC6!5*_lj*%z`R{q@5xg3-%d0WVy&5$I-XG{y z?RIedPhS0n#!txZRkkEv?IFjPBFFE&0q(~~)cf0K>=W^P4~lH=*u zq(u_~9)&-{-w)yM^M(@bGt_=HdQNZm+b~1hyTJQ+`neSiO(<;WRwg`PLPG5a?C)%4KeLRProX#Lk8^_k?3WJgX3+a3)!cNC_<{Uxl>+~oiO~bMTYrCa zslyZW`zoBy?B4s|Cc$K}QSDVT4*~szJb9(g-Z&%Uj)%@2U z?iPDgn*L_i7I2?Wt{Y2d*OTu}$?cVy_d@V{dc0l8qnBj4#!n5|uh7#Sg4;J5T85_| zJI>JIKKMGgK07o2Q(5-1}r2tSDP5W*K0XhZ!-K3-%kYh zk5yf&T7-^|dH)*Z4##IOe!zYXo)6>g_b+v7_(bky^kq-m(y6pn(R^^v2kzCueHgg! z0rz-tFAdMng6H4B^K}M0bfP_4PKIAq(xFsz_W#ZEws8TS>mAU`<^d(B8c>m<0kzH? z(Dc}V&X@!8zF?-W`N=DdSO}T5cw*&Bx1+h2MMmNM$SmYr?zN(7>!?k86H2#_q9Bq z^nN7WT`xnq=u3{SANLrmH8{x5np>+(e(Z$gOsh;Qv9~(;(y3Zs1m`3T~Ay z=+?n>aDR}w|LoF^hjex4TpF>D=nVT?^58SQ#bFC z<>>BL?chP4-^xAvi}<=vE%5f>-Ur-|fqN==zCS#F2cEBZiLQ1(IgRda_YCgh4gH35 ziaP=2|0AFpI|AyrG@u2O0y@?|pm%Ko%2_L*_9X&Zo1Jrzgfrj_XxzVkJ^0J7=qbMz zZ1g&-X`q~FIJStqmqZ{nzP9`VI;YJU`x0@xgR>-O2B*c)gyX^_2~c#OtlfL%*L6-N)x=pu_F{ zhEC@me10CDXZELWAj{9=9E{`6-OpsX7W6n(@cG}nRW-d^P4W5l;rgn~dqeWwkl(m( zwa}%!_`BUb$@)#`amaP?`1+@*+245Kd3;?`a4!Sy1Ht_yxF=)Y`!Me}nD>(K{8qZ# zoap<;7&rhwzwH_K>#uynd1$_%W~B~lhcBob9|MYi5KzYR0oC3Mm#+%w$n=05Lj(G; zQ-Cfapr>U6DxNoh4-V*yH$cwvYyShk(wz5e*j~T>TIE;48Gelz3deUsp^8P3Y4Ai%-A) z?bQZwKLYNT!ToP=H^4n7xHkj$_23>uUpp1tzju4(f9p};jz?w5@CV88vzO7`O`)?N z%z^~W>!oByD<*@2&50-sNCt0%bU!Ry_s z=+-23zFjIFhFi1Wz~#5dc*ohtqVHekv6mZ-zwb_-=Uzrydb~Q|J_FqEgL`qZ+(K{< zoMAt=lgtj!_aMW6M&HM_b}C^F?pWjJ^Jk^c`|x-*6MFQW z$D>&GbcbFW`uzdk4bAVx{J;9m(3EBDZNWZ%C|-cszeK0AmyGuc>|cQWaC+Q*e)c>c z-Q+Ykovd4N+uSO-fSwkgR~(O@5T4Ic)~&(0==SMw4>{bLPoJ~k27G?hrN^7V8Qj}| z`(bcT2kz6r-CMw=k!f6V!t;~i`7{ssdRP=AaE-#$U5XcbhA>OqYu9@L4fK?M>8Ro@oWnwJ4tZU$8MctF2x4JhHFfCh{Y z=&wG^c`N3;Ivig-AY}{C%LcUH#`Bkc#oT1(kNcH&t6vEg`t@#{UpsrF>u9_Sjh{%4 zkIUd!-7uQ})u+pNyLT79^*63nKK?9ydN>fBXWq9k?~|GLep%6Wc)l?_-wd7~N`~Ku zzK6km0=VaA-V=3U7uSHU?gu)({NSG!uOH;~IP&{*x_Gi&BDg;DPP&^_bar#-YT^3{ zgUImR$aXC~+5`T7;`6S+`}xUs%`$lupU9&xPJG~7b~X2yX=Xp$A%oBM849m3v}cx~ zMWgBOnEQ{-*yX_an+nqVq~>*l`Tv(3kIy@c$M3fu&o_%b?NGNG;O$b?q2I;V=gfiT zlj&C2@%%68=E40OxR>AT(gkp@I|g27S2qxzPfU*ctsvT-&Ls!FZWF%lr$Dv=RX z&vyY`zYE4^0-Cccpij#Ksyr3E2Z4S2fD+dUs7&dATILF*zkquK<~i=gjSV|^N8qQk-0l>_&8;C>z41HUrYBk}Rw>1&&K)wQx$KNj)oMi#Hqk>w(W zS4%$8#Xj?>5WYSe`?>k-X}00*hpu4W!GGu^^4xHbdf@fyG5^-u9=$1r&gbHOMoN!X zHX0Ct8o8BJl_L6Uru!Y)MD;p{9lJ36AbnagccI#eyx5~1!*P57B ze^43L2GxC5P)CNt^Ie1FEM6}U-xmmKUfLl0te_fx325VEwEH3)egIt81~hM0KzoM= z^hX!;zCl2%eh8?4zJLm*=8T~6EkBUg$?@Nx!K?2CbG+S!iEuo)p9S}I;64D{i-G(5 zke|DWKGlQgo%FSv(D&Bpd)CS1^!{{tt$cb}&8I*yv^_I;$N1C-+}ktn)8YBIbawML z($CNHs@*s*eUDd9+K}7ou)ksEo73AB$?8>Ge0^%aSKCZp)q4YXKVh!%_XRJ4Gnu{& z^KV*7pEnKfH;~NNo}GOok1EyjC|4!+HOzkbJYb(19-!Ar^}*0J?qjWn?|;~U&g1i! z)9EZF-+fn|j<=|x1sM(138U{{-Re)DbM6B7?bz8IUWJbX_b=eSge+f;dH0qEb2_{M z@O&|sTOa7}_R!l8KIu|ve4VkFUT!@6j-USz=PB1iT6r|2Lz_doF+Zd)V?#>QJEY3Z zLmF5G%!`EdIa5f*V?&%f-t%`**{=mP>PS!zHw9I4K72nmC|hr|y?Idks{|D+8syIy ze-|6nIdf27{thV7wSZC{A-iot+vkG!Xt3@PP!cqLens*dIlggP{JuY+qjNXS2ALZ3|y}hd58ZQU-x5|4} ztQb0g*3V(~%cOxH5|IPEaKIOjrafZj&!Xw<=4-F_sK|1U{+LgWo9@v>=q}Me=Dpn$-8QDX6o+Lzn5z?xEgYw)C zl1YO)xt%#*8dRrALGA1x)TcH<<){_ZPbGrtm>up2KU?**GN->kUcEB z`0C7j!OdQ6Sp@b|;P(;uyMA7!r>{Rro;`i#j z;{Me`kGfn&?|=7b@=fNiw7gp> z^MiXD_IP}?9Trxf&S6ceAJ)2ZVV%z#)`wJKrSgSU?IV~!2)U`nAOX|F=HvRf%{bQ zdop}oqS|EoGIaQPed>^g=Rv&QSFetu@ma64i}?eL@%4?^&&}M8uixp_+^t^KS;?88 z$J^N-PG~`oU!6Nx%>RQDWPp5L{ln}Z4RMF#t4BrQ{p)Cc@onHd*Q2(~`nyit(ZJ(f zL*pOjfa{ZbwAleK;OkSM?H`Ym-NF41dG0pa{*ZZp+yXCOm5vvlf18>2#KHmSyY+=z z@qh9;IegunjrhB{+{eMsukG;-=RdE8_4sgDaT~%)IycOHtFW5(2y0A}FuVS+-WCcg zZ-%gXhr>D?7gEf#klJ1eY44$s5^V^n_neS!jtZ%8kC0|I3F$?pkO~%t^E2?CFwf(H ziaZM{%M~7nf{Ixm)Qj0ctsjB*qw$X$;Ol+}@^?Q3^)_Wt$vx---tNEt<_)-S0r$EK z!GA10T~9i?#sOum5KtfVeRoQHAHFW{8$X{N0@Gu7b~?M-bIEZd{VLJbue;>;#by2K zm)ozt$^F{k;(WaEsSvoI0r!_X$aLttV;a5f5T6EiqMNVp)3kCv)yYSmqpLj~q1SQq zeh2Si_QxCS<$0bM-9HKcU%EiH{|!B7{-@*fmrnKSS2A9Fd*=Now4o$DZ%(hCCjfJN ze)q@RvmxL8v4@=;8b5ZbNBQvg@$Edi^OHxHN^&N^{Ux{?E;2v!{uA2XA5XV%JJ}AN zzctYyD|6oxea~By`?z9X{~f%Ln1`J^x3A#xn>>zx!@1md5!Ff@QG0ttqh5uz{AO6E zkB9YsYgm~Uh1FzySo8XXb*p7qS*nNCuXtGJvWAs0QCLH4VcmX-kGmPtyyIy5){v?! z3TeUkkRJ6x=UXxN)j12rLmHDcq|pgO>TC_E)IVhQ>p`tR;|rqkYv!Wiqk_uY9qb$7 z2G3F!2@{urp)_;zXGaxBA_)}=<4PLlm&gS)(yPzb#v)%_TuO7Bn`;!ptpPJ zSMO_R_+h^at@rETEWa8J^DAj5zx;Lms!`f6BbQ(0!9CvXml>XS(A!+3x5>7LeauR{ z{8XP-4)SR(yP1>VKC%i}^O#uOC-+ah&SOzk`nys--OB*~JIMfey*>}{cihXkc9EX` zAbZ&Pba_MR@zHx1n7=QD=cmtknS%F{>+iltx9{WS(fB`hbH5hctK#iu4fg0s2Y$!A zk1XTSncN;F#M9Rz%TIl8=*&HOS$1`K+0`{B$BmfEXJ_zrr`mG{YQqDixMQ7@OrM;) z8hmXZ9MQ235#6X0(VNl{I;DsTr--Pz5z*9lVRlntdCr9SJX2UxmxpzKN|<{S@O=9) zcN%!TbXco$g=I?-Rx=~a=Z`~j-3zJGnUF^A3TfB!5Px3?dOi@nZx_;B5c>cHte+TYo(f7@#*~hVyd%oHy+d?>fo=<K;?VmY3%sg7ft?LJpPX4ghtEq#pOYgm ze9xVXJwC7Mexk2M+gn}0!|n5^X|m89f62UU&8wFaac<(hGjYsR{gbM?VJ@>%!sf$b`9%9!?03R2&;X8u(qV-90bBj z`6Z+xk3%YdF{E<)L(0Aed}oHVeP~Ga(fIpy@%Cjxnw=}8lgUE5gSWf#I;c(H-V@vt zf%~t_yMuYJ56{o3Pp(7XFXtrF;p?u`-JE?vFLyPdr3V9=v?ico(*qhzejm}69A5)( zUkqK(98j%DKo`FHHTSV!zhCgH#Xh=uaGyX=yN{msF!MgNo?l78{S-W(3)~Cf>oUbL z+y8<+`qmomPwv_`rO&0q34hPtHY2;82>W}wyRY1_PkR+UKl}~n6~mh7>^p(*!f{y9*frZ4y$%^epe-|rA5F!Q&<~f-~o+ZDyzf5-f9%7Wg*{UyD;PG_I_A$M=U z{l_EhXtpx*3)s($@+l#j-xi*4U4nCxgKi((_c80HPwd_vGv|2xT==|B8`#ay6?fU!*?~h>a!F?K@-ow1N%^A|_By@K+ zy1N%a{d1LBK8Sx?Lzc(SR~>}MYm3*f9@ONbWVwv=Iw5lXXL{HB__)&necc{VuEhZ@ z`6Zx6y~*#*0(w@N9W1zapr5ZuKVKxyuhvh=arE;i(DyW3{5nQ{uQ;3g7n9(4eEl=d zz_I>hIQsgvt>|)VljS%EtMT|TWcnU4WP0|qRhj=USA4q5y!Tq=(`@EFJ6>-KxG!zs z)A|ZNEn(gpXY}cx7_y#)e(xWzR*~yB9>>pZ<1X!DxO^h_YT^2`I7h<$Q zM#0W8%3D81&C11SM&1~mO%5~y&MtD@udCXQeq@5{U|JqWLF;5mAKZ7d$|XR;hU-Y|In92#$B=2IoZ6T0bi-m#av z$DP_U;7+D%2G=K;PL>kYTP42&Wh-8rHJkpis)tfh@OWcbk$)kd=gfTOJTj< zA6Cybc=j21`=Mb~?TF^rfe%WB<)^>d1MZc4>|7YXZ#%bJslnWd3u|1i}3K{g38mA`EHEA$IA~O!qy$jmQPKt=Q8^|J=#A;+J@<4y5$F5>a*ul)M= zF5Vu^|8&r=F(={BHwlQ z>Ib^}t(BO0`rEGbH@~Ct|KU8WU#w=fj@7CfvDzPv)zxgV`jRMC>Fu$q_9{jrZpP@~ z@fbO_#i;zE7)>7^qq}`#l)qIB83S%F9>cycM!zPC(Q6z0|B8HeGoq8nBeHFaDEA^f z-uQ^>^&!8tjHp;OIG|WWN6`4@X#7)CMAe>$HRB3+9t`Ug-fr7WG8?$(2KU3vdtT;! z)OTSWN*LB%GkKp3fBjNOTlc~BD?`eSpTFG?f8R2sj8#H9TqvZbY4QE^HoM+|`)%g^ zL{JsB1XTwwKMPN1bPMWm13X9$+W8&Zkv)`n72Mk{)!rP4V+T`V=jK zm(NW0izUxn$npQs-O=ah824|s(B&9iB-Clu_`hRYzN1x zdxuzh{aBqU9jo|UvE-~+^)+I3@LdcX0)L-j-gm{ww1S_f#%RN!7`Z#dsA(O3E*+x} zxnlTtn0q5e3*O<~?nZRtG}!LM*DZ@^(Bz2n55Uu-@%3x+dK7P$Euw-6BPwNyDA(U% zS+3&Uz`fLexHIqGAz_t+=iAf@t26rEI7?WC@O37-o9&P2X3vG>CdVzIt1U);zuJp4 zKyQ<&0$HBE?r;iv8yB4pUOpLKe$3IJijdt4;N{1Id+F}vH|9NirJ(v3X67>mRWBjE zturVm*x$k1%_Ga%&jz&PP(V$7<2-N=qs82Sx=je^3jTgZmw@&)A=g!9ClB_@&h9oT z?}z*Qyk!r=%qKg;{kq+FfR)^}A>U<(>woP+o^QmCwxXXqZg74^vOOCAv)QjXX#9rj z^z_GkTD_GHZxMGbC-_vdAA1};{mGi#!zcmI!}Z&ea$m#c(~&oH@^|p@e|Qz7yZdjR zd!RrQFeP1C~Lkh&|Oxjqb z3C61Dml&OT9HVrXV)$H9j4rR`=h-o8KO#nJyE69;W8|+8qbvntR3I%J5P2IMKliFqu_oXPyYk+zW6fy4$s?GgE{*CX+T() zTI1uYh1C#0AE2|_M2@RYSNqT1ZykO1Z6W>1Zf4%-kcM@mlWh>v`v1?%cg4%6!^@`w z_nzRM?>IPb4Jr!Hub2?TXVKrZ52^(@?rTN%GR1=GN|tl5m%Gh=W{)eVR(QOPf0N

f-t7axUa0_opWB2iVc!^Y=aRD+7LS{9*QR z8`-_j_iO7odYeA<`1rfQHRx>R*B{y0-zCKtxZwJ?}`DbdcNNMoItjg36!W>0%b3jKozqjQ0D{*G{cfWhhM+Ed~B>r^^Vny=COKSg`bPY>i?8<*3nfRS`#l;D3W{c@pEx^D7+R6 zPzFkoqAA6SyS2F6$D_ChDTF|2ap#rb4nbQqcyT8b4KC&Tv0yD0izNTenb|Y5_dSWp zZ|kYed{1sG`T=wOyrqo$Ek&)h_Ppt;q(lA{xDPLz`Yx|--GWp3Yr>=&aZh-AIp6IBf(HKzOD`3O&VQla4vSZYK5;UEBfw{UfA=SbF$TA;sb4yEh4`2e=<#2R{zJpJE5U zoE>~F=6kJ#pbFvV_x+Dv4(ubQlH}J0nh3@wX=vhnduxi8gviuBwZw#|ODJ8z(2^sHlgo203 z^_$`O;s_nb-#r@!-^uo76EweCgwmJ559DBH8{j76Ew`3;;qz(ku6DD3U5&>_>xZHB zZ+gM~mLGBMd&JStn;p$q6ty)rqpkSVwjLx|GOobmzbs|ii0&`2bY+sI z*nyV5X2ze0;@|R?rscQvr)B9n-tI!8smI%bY?vBy$`)z3vygFLka&ebT~7;O=?3Wo`kgcLP(YOhZF?& z@ErE;qeH6LE2N1nLh3?4f4xjdsq%+(96w(f&PRd$quW8HIR~ei^MS2k%`U#qujKd1 zLERt5uC537H?7I|tKToSwhy=Ux~ok#vw5#I z-~Ynajn!|x#$HyRc=S{jJAYes*Q%jSFF-Vodm zoy5oQAkU%mU8b7qL56SLo<658&q5^6MlPO*2vhZ5f;(O9zCR6(r=xF1Z&QAZp~^kU zcTIT)su*fioGusK?}7VgNqBhpUhz;!Gtl{-^F!o!}$fFQ)& zm-nrM`zNld{^O|MSx2km9VM=Dly44rk8u>&+tI%*9o4Gk=+`oilJYpbZ`e`nds}z! zk6NB0`~&&eaRPO`E2}LdP6Ii@1NZc@%e7}O^#at z=S#7hDG|+`1^d`M%=w$GL)ug?r2kgt>m@?kn}<0M_p8GFVsQVjTX?V)YY|k zeRO{!ypI`zzhmZCMDw3l;Rc6YO-bhco*?%caDFG8FAnE_;~wMvTJ|)|`Phlv+zjLn zqazuwK7C%LfEsRln1_hYzfm&(0@$k7##@1;WfAKU73(bnLD@OquCp1;|uH^El8zfCW}^T3S1 z7YXKhZA~}u^zYDeyxraeOLMpL41oLFsTQBzwREEcUhX?foy%Do%Y2_1;yJ+A4WPTJ zf7DdrO{UV!H+6NqsU3ajbXu62QQg!WdYeD9n9BMY&p{HI9NhCBB)7#G>JQ)3($hXd z=f7@k$fT!jS`K|L$exb;J}Ql&^DpS-z&_>-etkb$&JD))CFJ#4AqJ-E zh#Xffq{z}C`jL=^gMGD6(f_A#8{Strj<4VG@&4)+-R)X@ek11nlh#2EtBbGyihtq! zOE~`noZmvWPY36tj?&$14bY+E~`q$Hqz%c z#NWN#!+qtNk2sIc=4obnPb*%!+H=F@xp0+io2xuaT}4fIHF%h-l|Q<=(a2SfZ(Ox0 z?rKpMR}a(jT)c2J;ku*CM;#U0>S*W^w0@eS48t7N`4J!A$kFD?j;Fu78-LI0*598&S?{gQ?&rdZpc^H1KtD#}+W@gr8R|EDZ3o++28M+1UgYVeK z-sk@2B7W{L`?qcE<5%M6@p+-~+*u7o2f+SSLp)v;?(Mz|DSNJvzBJ%G9iB;ddxhO> z3_ic&hM>+bz`sWab$2Kp|3~hv;Cz{?=m4BQmm{dq2hP7^h9`pQ8D@Gs{SDrJ$Xqhs zxPVIc3n*pVfa1w?{mL`<3o_sFckS?ZA9Op|+f%NVo=Vm9ct@b8wz)n19P+g6ovVbq zu0A>K>hs;M+O9&|V_aPw=_+?mSDl-?+W4)jl%-tN$>nNZ(AA^2j!NBiH29RGt-Bn( zUJ3VOnE4|ex&=qQnxXsE9n~x0C}$2wmzeP*nepyTTjP$|y0yht&c(K>M%$`AnA}c( z6Y-s`J>cF3+;7A8DEPhvoxgFMu5PQPqKhn*oouNB-A!a0`r2BSikGtVWe)nA^p@H@ zXTN%dKIRa5zs}UB;C>F=uk|+Nwlp=O2D!Zyy*=}NOwiOf@9^r0hCU^~XJR)q8_o~F z+x;)vP%$$6?e6q+E$Q^?@+^q^i(-bp$<7Q2^8Zf^y`!_ObQ#=_z?6&q}U)`6oA58}aT7@%PaIJs!g6Dms(l@pm6^o_WyI%5|RN zfAe&9f~Uv*J*98wsc4j^hULM$fTuXe)1Bm=az1p`@^3KR55L#CDm>3s|M9L4^mFC3 zb=9h_t7YZjdVZb_$K`jO@Ussb?LUu>?}N*KFz4qwIyu(SvOZwm%26g}d|YWqMal8g zf{y-vV=Mh_^7#o{&NlcC?mNM~Gq}Hj?|q_dT_|g-KweuNLUeh|_YF7c;OTCzZL*X! zpYD#%E^|Lim0DTqRm0+4k(P32C)201)DhfkUNzMozOUF|swg`Dh`#O#^Zi?VU9x(n z!ot%PVPBKUBx{)(`U-u&V`vkdU6K8Uwy&qBhxeIg;P*xtYS$Z2N2ZIaXDEj5E*_t^ z2kd_Z``YhAIt}m3;PL;9*IW5P^B2+OO+)jU_dj(DsZcX|ditE6=zXfZWIziYz~^=3 z9%CPSItP#66wc3r^YL`|*O~K4or8McFsN6RgSw5k-<6dct@LEPS9rQx^s*=5^$xmw zGTjH9b7b*($GopP&poxd=4sFoPg6I8?IKTCqCKS@3fDV(>e9f|f=Ya?sHZHMJhe*W zX~|PpkFU7;>aa^c>+1YMv^^SaAA+`b;#sH<{*}=8qRjn_;GP;CNOJfbl%p^Hbo9vv zhj-?~?}?7iljEzlbu_n*qqC8Yp5pD^(BIsBYikF%_XT$k+!um-MsOcGh)$;+nXay_ zO3e3exxqZZGx5sO{p;lSBjk5-T+Hv5UX8cZm;CO5`%Q2^Pj8bvCtV&oe}~;n9P@n} zyBTM@sr}1LCH!Kl;t06klbH|a7uV*wAj@?uY-%ppKgQoRe`Dwzysyhnb_=uKLGv3g zF*KC!Zv04gGhiPH_Lt~!TEP3A^z~`7;0e(DYOh18a)-Y5G~H}`h;D_>o<1l01okxp z(EAQN5A^mKE79q(r#q54L>Itb2hNX$^UF?fXUClXa|O@AEOvGygPMxB??j*XA2NMX zX?`XTGv4G`c>fXSH9dUoY3A!}bzhH4`1C!#G6j5y zt84#pwL7b;sp(wRV8%bcMutCvhu`F=$wEi-Cp$WZx4YGj*(BR>QcA@!qn$zVm@BfW7l$hV( z_ZpG)QyD5pzPoh=ywUrVTSFSO6n{4Z56_-%7M;#jIPa0~YQgz#`1=Vacel*>8TZNR z=eecbA5_gaI-TFR!I%)__dEqvzCB$3jy{J!W%m23w$@jRUw!o%=WB9bUrSp1ymQva z2m1=<^Hss})%_D+3;zM%vz{{UMbrQ1ap&Rb#u!gU`grQ!%G35*9-j#W@4TKSo1Tum zca{7fH2f?xJRY8}an&u>Rj1MXzZV|Ag{wqn{H&6$isW?lM}(`auO0Qc>1gROM@P3f zy10n$W{RUlgYfludN&H42X}${7WiHm+$WLY?=s)Te4j>lb7vCW+yGn4@$)&reL1+- z2lp)SJ@A^`hR(Mkzn>w$_b0#qfS>PAjyvAjQlEzSc{=*l#VwU%ADa%&=YLBtlV~dA zIdVET7`yRzyTSe}*xv*DWbi(1EBH^AzX$fa(EP!fO=VBdPX0BX?Jl^Z`E&Qu<&o_U z|HfVIBrqOqXliHlzLB8}<04PQ#hX=&R2)? zt=sc;vRzl^eA^=I@-mX=(~t+Bq3u`M-5lZWa!XLB7JtO~K*q36rw%JI$ybUizWl#@ zeZJ9Gtp&b%PWCl@5Pa|8E2*BZJo44Fu&)^3*U41AtS9jM5_*2f)BW|HzMPNlPx7>L zpr>RVJbtg0r+yVYEiB~ekmu=IN>5K7;nyy>x^lqPjyP9inDG^u@z==l-O2HZwcvOe zGJS41kGH$<2ETU8;kFn}-|DEvVn@Mfx;tjtKA!(m{63sd zel?_=aK0ito2VrreLo%Cne!#Pa*Nv}qK-$G=3lOg zZ*;Y8zN=i5T=f~?Y75@(dL4Qj`kSYDTwM&gihGN1zwIc`33v|foxuI|6#6-EPXq1~ z!QBA&+2CFr+)sl0X!!mEI$sN&Z^nF|#(ckzuj|{94)1&By=;{)Y-^TptMI3Ew@J)* zy4pTyy#eQMljZKsKnF&!oA1ePhB^O1@0-)#WGQ0lVkS#-@OZUeb5BEuGxLn8ocl~I zjDzR+yzi%&nljYnci*v_`+<8avfb{Id_Ol`kI8*rGE)y97@B^8e*Tc5wHxSf7t-J3 z@soy<>wl!jX^hW@^NDai1DwwV=Zn8(mf!t2=QGjiyj@LJk0s}m>Ff0gDJ^}@mfGCq z@aIl^Sg+TFm2OU0S;mA_s&`m*T87oBW>}+3hqWwsSjR$Py?o~@U!t!jr+rP`?d#Yo zU+HJV`%%84d-}X9##go);9bhsqFiWx&{u(X9vwUL{FJ8yyFA@n=_zH5r-+g0dk>z8 zrk;*6<71fdQOx)k#=;!)CaTxAUfu<>+rRws-Yu6dW&$&&$gkFz9&R zk>S9-F}S|~_o?7s1il{w_d(!ZBFfQwaKDw;(LV-S@UN|!?B(K4vV+-aE7=O0_n4sX zKij;&+*azQ?BlB0I#8reF=XRgm_P3?p2l4hB$Z_PkczU}=!!14PZfQ&l zOWA6(m!-ogUWgl=u%$Gg@%2vB-`WdU#%z4&;N*XKR=^JrOK$8kHhT6!-_u$m)C{$^tZ4gCWcjHKv+@j!WtL_ z*UN`>sQ{_N4XbvFu*N;~b>MGbX%FDny_p`kG(Y*OhX> zU44tEFErTIzu>+b+=r6oO60@mf%|N5uXD#yYI^zu+Z^6s?r7jtM}67POz!AtC!Jkt zI=j9F;l1nVTymZX_HxhJ%hla)>kNCjsq<`2ptIYG)@LHauY&U(+0T4e%vS5Hd_RM& zOXRuI=>6|!;rxE?G1k%B;qkiA;hg$Eo{O%`{HB&ReQV(hEv3Zgd5)!UO7!3fH+7dy z9Y1Vp`W9~MmYQ^x==#sxTH*6nG$-5DG&Q^oGa;|3{5E|Ld--_}4ef{Xcky`1;C#9T hA8+s8Gw1JhVwcmv(7dk=bt2Oj!rLdNWuHrz_y5K?=Wzf4 literal 0 HcmV?d00001