From 85eec6108e844347aaabc6b14625ec5d200a299e Mon Sep 17 00:00:00 2001
From: lleplat
Date: Mon, 12 Feb 2024 12:41:07 +0000
Subject: [PATCH 001/866] Make AbstractTileableImageServer.readRegion return an
empty tile (instead of null) when AbstractTileableImageServer.readTile
returns an empty tile
---
.../images/servers/AbstractTileableImageServer.java | 12 ++++++------
1 file changed, 6 insertions(+), 6 deletions(-)
diff --git a/qupath-core/src/main/java/qupath/lib/images/servers/AbstractTileableImageServer.java b/qupath-core/src/main/java/qupath/lib/images/servers/AbstractTileableImageServer.java
index 6ac763628..a8c547049 100644
--- a/qupath-core/src/main/java/qupath/lib/images/servers/AbstractTileableImageServer.java
+++ b/qupath-core/src/main/java/qupath/lib/images/servers/AbstractTileableImageServer.java
@@ -397,7 +397,12 @@ public BufferedImage readRegion(final RegionRequest request) throws IOException
}
}
}
- // Maybe we don't have anything at all (which is not an error if the image is sparse!)
+
+ // If we have an empty region, try to use an empty tile
+ if (isEmptyRegion) {
+ return getEmptyTile(request.getWidth(), request.getHeight());
+ }
+
if (raster == null)
return null;
@@ -428,11 +433,6 @@ else if (xEnd - xStart <= 0 || yEnd - yStart <= 0)
raster = raster2;
}
- // If we have an empty region, try to use an empty tile
- if (isEmptyRegion) {
- return getEmptyTile(raster.getWidth(), raster.getHeight());
- }
-
// Return the image, resizing if necessary
BufferedImage imgResult = new BufferedImage(colorModel, raster, alphaPremultiplied, null);
int currentWidth = imgResult.getWidth();
From 2fa9b6cd3deac0bfcd51778427668ecbd283ea36 Mon Sep 17 00:00:00 2001
From: lleplat
Date: Fri, 23 Feb 2024 14:53:51 +0000
Subject: [PATCH 002/866] Added possibility to open Zarr images
---
gradle/libs.versions.toml | 2 ++
qupath-extension-bioformats/build.gradle | 3 +++
.../servers/bioformats/BioFormatsImageServer.java | 12 +++++++++---
3 files changed, 14 insertions(+), 3 deletions(-)
diff --git a/gradle/libs.versions.toml b/gradle/libs.versions.toml
index 0317ca9c5..fb87dec50 100644
--- a/gradle/libs.versions.toml
+++ b/gradle/libs.versions.toml
@@ -2,6 +2,8 @@
bioformats = "7.0.1"
bioimageIoSpec = "0.1.0"
+omeZarrReader = "0.4.1-SNAPSHOT"
+blosc = "1.21.5-SNAPSHOT"
commonmark = "0.21.0"
commonsMath3 = "3.6.1"
diff --git a/qupath-extension-bioformats/build.gradle b/qupath-extension-bioformats/build.gradle
index 62fce96c9..0f9a032d4 100644
--- a/qupath-extension-bioformats/build.gradle
+++ b/qupath-extension-bioformats/build.gradle
@@ -34,6 +34,7 @@ dependencies {
implementation libs.qupath.fxtras
implementation libs.controlsfx
implementation libs.picocli
+ implementation libs.jna // needed for OMEZarrReader (see https://github.com/bcdev/jzarr/issues/31)
implementation "ome:formats-gpl:${bioformatsVersion}", {
exclude group: 'xalan', module: 'serializer'
@@ -45,6 +46,8 @@ dependencies {
exclude group: 'com.google.code.findbugs', module: 'jsr305'
exclude group: 'com.google.code.findbugs', module: 'annotations'
}
+ implementation group: 'ome', name: 'OMEZarrReader', version: libs.versions.omeZarrReader.get()
+ implementation "io.github.qupath:blosc:${libs.versions.blosc.get()}:${nativesClassifier.replace("natives-", "")}"
// testImplementation("ome:bioformats_package:${bioformatsVersion}")
testImplementation "ome:bio-formats_plugins:${bioformatsVersion}"
diff --git a/qupath-extension-bioformats/src/main/java/qupath/lib/images/servers/bioformats/BioFormatsImageServer.java b/qupath-extension-bioformats/src/main/java/qupath/lib/images/servers/bioformats/BioFormatsImageServer.java
index a9b216ebf..86bd8bd83 100644
--- a/qupath-extension-bioformats/src/main/java/qupath/lib/images/servers/bioformats/BioFormatsImageServer.java
+++ b/qupath-extension-bioformats/src/main/java/qupath/lib/images/servers/bioformats/BioFormatsImageServer.java
@@ -37,6 +37,7 @@
import loci.formats.gui.AWTImageTools;
import loci.formats.in.DynamicMetadataOptions;
import loci.formats.in.MetadataOptions;
+import loci.formats.in.ZarrReader;
import loci.formats.meta.DummyMetadata;
import loci.formats.meta.MetadataStore;
import loci.formats.ome.OMEPyramidStore;
@@ -1185,10 +1186,15 @@ private IFormatReader createReader(final BioFormatsServerOptions options, final
}
IFormatReader imageReader;
- if (classList != null) {
- imageReader = new ImageReader(classList);
+ if (new File(id).isDirectory()) {
+ // Using new ImageReader() on a directory won't work
+ imageReader = new ZarrReader();
} else {
- imageReader = new ImageReader(getDefaultClassList());
+ if (classList != null) {
+ imageReader = new ImageReader(classList);
+ } else {
+ imageReader = new ImageReader(getDefaultClassList());
+ }
}
imageReader.setFlattenedResolutions(false);
From a7b184f2bf0134432c52080e394ab51fea02dfe3 Mon Sep 17 00:00:00 2001
From: Leo Leplat <60394504+Rylern@users.noreply.github.com>
Date: Thu, 29 Feb 2024 14:22:31 +0000
Subject: [PATCH 003/866] Fix accessing z-stack pyramidal images
---
.../images/servers/bioformats/BioFormatsImageServer.java | 6 ++++++
1 file changed, 6 insertions(+)
diff --git a/qupath-extension-bioformats/src/main/java/qupath/lib/images/servers/bioformats/BioFormatsImageServer.java b/qupath-extension-bioformats/src/main/java/qupath/lib/images/servers/bioformats/BioFormatsImageServer.java
index 86bd8bd83..bbff2702a 100644
--- a/qupath-extension-bioformats/src/main/java/qupath/lib/images/servers/bioformats/BioFormatsImageServer.java
+++ b/qupath-extension-bioformats/src/main/java/qupath/lib/images/servers/bioformats/BioFormatsImageServer.java
@@ -1376,6 +1376,12 @@ BufferedImage openImage(TileRequest tileRequest, int series, int nChannels, bool
if (tileWidth <= 0 || tileHeight <= 0) {
throw new IOException("Unable to request pixels for region with downsampled size " + tileWidth + " x " + tileHeight);
}
+
+ // OME Zarr files provide z scaling (the number of z stacks decreases when the resolution becomes
+ // lower, like the width and height)
+ if (ipReader instanceof ZarrReader) {
+ z /= tileRequest.getDownsample();
+ }
synchronized(ipReader) {
ipReader.setSeries(series);
From d475f92c2f01f361dc3dcb0067afbe5eb2f6bf3b Mon Sep 17 00:00:00 2001
From: =?UTF-8?q?L=C3=A9o=20Leplat?=
<60394504+Rylern@users.noreply.github.com>
Date: Fri, 1 Mar 2024 10:36:57 +0000
Subject: [PATCH 004/866] Update gradle/libs.versions.toml
Co-authored-by: Alan O'Callaghan
---
gradle/libs.versions.toml | 2 +-
1 file changed, 1 insertion(+), 1 deletion(-)
diff --git a/gradle/libs.versions.toml b/gradle/libs.versions.toml
index fb87dec50..9f0925798 100644
--- a/gradle/libs.versions.toml
+++ b/gradle/libs.versions.toml
@@ -3,7 +3,7 @@
bioformats = "7.0.1"
bioimageIoSpec = "0.1.0"
omeZarrReader = "0.4.1-SNAPSHOT"
-blosc = "1.21.5-SNAPSHOT"
+blosc = "1.21.5"
commonmark = "0.21.0"
commonsMath3 = "3.6.1"
From 21b43e0a54b070c19e6c556b00d1fdef66be1489 Mon Sep 17 00:00:00 2001
From: lleplat
Date: Fri, 1 Mar 2024 11:21:07 +0000
Subject: [PATCH 005/866] Improved fix accessing z-stack pyramidal images
---
.../bioformats/BioFormatsImageServer.java | 20 +++++++++++++------
1 file changed, 14 insertions(+), 6 deletions(-)
diff --git a/qupath-extension-bioformats/src/main/java/qupath/lib/images/servers/bioformats/BioFormatsImageServer.java b/qupath-extension-bioformats/src/main/java/qupath/lib/images/servers/bioformats/BioFormatsImageServer.java
index bbff2702a..5ddfeb606 100644
--- a/qupath-extension-bioformats/src/main/java/qupath/lib/images/servers/bioformats/BioFormatsImageServer.java
+++ b/qupath-extension-bioformats/src/main/java/qupath/lib/images/servers/bioformats/BioFormatsImageServer.java
@@ -1376,16 +1376,24 @@ BufferedImage openImage(TileRequest tileRequest, int series, int nChannels, bool
if (tileWidth <= 0 || tileHeight <= 0) {
throw new IOException("Unable to request pixels for region with downsampled size " + tileWidth + " x " + tileHeight);
}
-
- // OME Zarr files provide z scaling (the number of z stacks decreases when the resolution becomes
- // lower, like the width and height)
- if (ipReader instanceof ZarrReader) {
- z /= tileRequest.getDownsample();
- }
synchronized(ipReader) {
ipReader.setSeries(series);
ipReader.setResolution(level);
+
+ // Some files provide z scaling (the number of z stacks decreases when the resolution becomes
+ // lower, like the width and height), so z needs to be updated for levels > 0
+ if (level > 0) {
+ ipReader.setResolution(0);
+ int zStacksFullResolution = ipReader.getSizeZ();
+ ipReader.setResolution(level);
+ int zStacksCurrentResolution = ipReader.getSizeZ();
+
+ if (zStacksFullResolution != zStacksCurrentResolution) {
+ z = (int) Math.round(z / tileRequest.getDownsample());
+ }
+ }
+
order = ipReader.isLittleEndian() ? ByteOrder.LITTLE_ENDIAN : ByteOrder.BIG_ENDIAN;
interleaved = ipReader.isInterleaved();
pixelType = ipReader.getPixelType();
From d7fc43a2b7728a03f97824195ef480519c186e40 Mon Sep 17 00:00:00 2001
From: lleplat
Date: Fri, 1 Mar 2024 11:32:29 +0000
Subject: [PATCH 006/866] Improved fix accessing z-stack pyramidal images
---
.../lib/images/servers/bioformats/BioFormatsImageServer.java | 2 +-
1 file changed, 1 insertion(+), 1 deletion(-)
diff --git a/qupath-extension-bioformats/src/main/java/qupath/lib/images/servers/bioformats/BioFormatsImageServer.java b/qupath-extension-bioformats/src/main/java/qupath/lib/images/servers/bioformats/BioFormatsImageServer.java
index 5ddfeb606..cc9d98869 100644
--- a/qupath-extension-bioformats/src/main/java/qupath/lib/images/servers/bioformats/BioFormatsImageServer.java
+++ b/qupath-extension-bioformats/src/main/java/qupath/lib/images/servers/bioformats/BioFormatsImageServer.java
@@ -1390,7 +1390,7 @@ BufferedImage openImage(TileRequest tileRequest, int series, int nChannels, bool
int zStacksCurrentResolution = ipReader.getSizeZ();
if (zStacksFullResolution != zStacksCurrentResolution) {
- z = (int) Math.round(z / tileRequest.getDownsample());
+ z = (int) Math.floor(z / tileRequest.getDownsample());
}
}
From 00d05ffaefa7e69ab7fc15a7e80d1ed180387540 Mon Sep 17 00:00:00 2001
From: lleplat
Date: Fri, 1 Mar 2024 13:53:03 +0000
Subject: [PATCH 007/866] Addressed comments
---
.../images/servers/bioformats/BioFormatsImageServer.java | 9 ++++++---
1 file changed, 6 insertions(+), 3 deletions(-)
diff --git a/qupath-extension-bioformats/src/main/java/qupath/lib/images/servers/bioformats/BioFormatsImageServer.java b/qupath-extension-bioformats/src/main/java/qupath/lib/images/servers/bioformats/BioFormatsImageServer.java
index cc9d98869..e3b7705f4 100644
--- a/qupath-extension-bioformats/src/main/java/qupath/lib/images/servers/bioformats/BioFormatsImageServer.java
+++ b/qupath-extension-bioformats/src/main/java/qupath/lib/images/servers/bioformats/BioFormatsImageServer.java
@@ -1379,19 +1379,22 @@ BufferedImage openImage(TileRequest tileRequest, int series, int nChannels, bool
synchronized(ipReader) {
ipReader.setSeries(series);
- ipReader.setResolution(level);
// Some files provide z scaling (the number of z stacks decreases when the resolution becomes
// lower, like the width and height), so z needs to be updated for levels > 0
- if (level > 0) {
+ if (level > 0 && z > 0) {
ipReader.setResolution(0);
int zStacksFullResolution = ipReader.getSizeZ();
ipReader.setResolution(level);
int zStacksCurrentResolution = ipReader.getSizeZ();
if (zStacksFullResolution != zStacksCurrentResolution) {
- z = (int) Math.floor(z / tileRequest.getDownsample());
+ z = (int) (z * zStacksCurrentResolution / (float) zStacksFullResolution);
}
+
+
+ } else {
+ ipReader.setResolution(level);
}
order = ipReader.isLittleEndian() ? ByteOrder.LITTLE_ENDIAN : ByteOrder.BIG_ENDIAN;
From 59ac700c3785df3cbe59e3b886bf0e739f64c466 Mon Sep 17 00:00:00 2001
From: lleplat
Date: Mon, 4 Mar 2024 09:50:23 +0000
Subject: [PATCH 008/866] Started zarr writer
---
.../writers/ome/zarr/OMEZarrAttributes.java | 133 ++++++++
.../writers/ome/zarr/OMEZarrWriter.java | 322 ++++++++++++++++++
2 files changed, 455 insertions(+)
create mode 100644 qupath-extension-bioformats/src/main/java/qupath/lib/images/writers/ome/zarr/OMEZarrAttributes.java
create mode 100644 qupath-extension-bioformats/src/main/java/qupath/lib/images/writers/ome/zarr/OMEZarrWriter.java
diff --git a/qupath-extension-bioformats/src/main/java/qupath/lib/images/writers/ome/zarr/OMEZarrAttributes.java b/qupath-extension-bioformats/src/main/java/qupath/lib/images/writers/ome/zarr/OMEZarrAttributes.java
new file mode 100644
index 000000000..087e02ac8
--- /dev/null
+++ b/qupath-extension-bioformats/src/main/java/qupath/lib/images/writers/ome/zarr/OMEZarrAttributes.java
@@ -0,0 +1,133 @@
+package qupath.lib.images.writers.ome.zarr;
+
+import qupath.lib.images.servers.ImageServer;
+import qupath.lib.images.servers.ImageServerMetadata;
+import qupath.lib.images.servers.PixelCalibration;
+
+import java.awt.image.BufferedImage;
+import java.util.ArrayList;
+import java.util.HashMap;
+import java.util.List;
+import java.util.Map;
+import java.util.stream.IntStream;
+
+class OMEZarrAttributes {
+
+ private enum Dimension {
+ X,
+ Y,
+ Z,
+ C,
+ T
+ }
+
+ public static Map getGroupAttributes(ImageServer server) {
+ return Map.of(
+ "multiscales", List.of(Map.of(
+ "axes", getAxes(server),
+ "datasets", getDatasets(server),
+ "name", server.getMetadata().getName(),
+ "version", "0.4"
+ ))
+ );
+ }
+
+ public static Map getLevelAttributes(ImageServer server) {
+ List arrayDimensions = new ArrayList<>();
+ if (server.nTimepoints() > 1) {
+ arrayDimensions.add("t");
+ }
+ if (server.nChannels() > 1) {
+ arrayDimensions.add("c");
+ }
+ if (server.nZSlices() > 1) {
+ arrayDimensions.add("z");
+ }
+ arrayDimensions.add("y");
+ arrayDimensions.add("x");
+
+ return Map.of("_ARRAY_DIMENSIONS", arrayDimensions);
+ }
+
+ private static List
+ *
+ * Note that the image server used internally by this writer may not be the one given in
+ * {@link Builder#Builder(ImageServer, String)}. Therefore, the {@link ImageServer#getTileRequestManager() TileRequestManager}
+ * of the internal image server may be different from the one of the provided image server,
+ * so functions like {@link TileRequestManager#getAllTileRequests()} may not return the expected tiles.
+ *
*
* @param tileRequest the tile to write
*/
@@ -117,11 +131,7 @@ public void writeTile(TileRequest tileRequest) {
try {
writeTile(
tileRequest,
- getData(
- server.readRegion(tileRequest.getRegionRequest()),
- tileRequest.getTileHeight(),
- tileRequest.getTileWidth()
- )
+ getData(server.readRegion(tileRequest.getRegionRequest()))
);
} catch (Exception e) {
logger.error("Error when writing tile", e);
@@ -139,6 +149,8 @@ public static class Builder {
private final String path;
private Compressor compressor = CompressorFactory.createDefaultCompressor();
private int numberOfThreads = 12;
+ private double[] downsamples = new double[0];
+ private int maxChunkSize = 12;
/**
* Create the builder.
@@ -180,6 +192,37 @@ public Builder setNumberOfThreads(int numberOfThreads) {
return this;
}
+ /**
+ *
+ * Enable the creation of a pyramidal image with the provided downsamples. The levels corresponding
+ * to the provided downsamples will be automatically generated.
+ *
+ *
+ * If this function is not called (or if it is called with no parameters), the downsamples of
+ * the provided image server will be used instead.
+ *
+ *
+ * @param downsamples the downsamples of the pyramid to generate
+ * @return this builder
+ */
+ public Builder setDownsamples(double... downsamples) {
+ this.downsamples = downsamples;
+ return this;
+ }
+
+ /**
+ * In Zarr files, data is stored in chunks. This parameter
+ * defines the maximum number of chunks on the x,y, and z dimensions.
+ * By default, this value is set to 12.
+ *
+ * @param maxChunkSize the maximum number of chunks on the x,y, and z dimensions
+ * @return this builder
+ */
+ public Builder setMaxChunkSize(int maxChunkSize) {
+ this.maxChunkSize = maxChunkSize;
+ return this;
+ }
+
/**
* Create a new instance of {@link OMEZarrWriter}. This will also
* create an empty image on the provided path.
@@ -193,7 +236,7 @@ public OMEZarrWriter build() throws IOException {
}
}
- private Map createLevelArrays(ZarrGroup root, Map levelAttributes) throws IOException {
+ private Map createLevelArrays(ZarrGroup root, Map levelAttributes, int maxChunkSize) throws IOException {
Map levelArrays = new HashMap<>();
for (int level=0; level createLevelArrays(ZarrGroup root, Map DataType.u1;
@@ -254,7 +297,10 @@ private int[] getDimensionsOfImage(int level) {
return dimensionArray;
}
- private int[] getChunksOfImage() {
+ private int[] getChunksOfImage(int maxChunkSize) {
+ int chunkWidth = Math.max(server.getMetadata().getPreferredTileWidth(), server.getWidth() / maxChunkSize);
+ int chunkHeight = Math.max(server.getMetadata().getPreferredTileHeight(), server.getHeight() / maxChunkSize);
+
List chunks = new ArrayList<>();
if (server.nTimepoints() > 1) {
chunks.add(1);
@@ -263,10 +309,10 @@ private int[] getChunksOfImage() {
chunks.add(1);
}
if (server.nZSlices() > 1) {
- chunks.add(Math.max(server.getMetadata().getPreferredTileWidth(), server.getMetadata().getPreferredTileHeight()));
+ chunks.add(Math.max(chunkWidth, chunkHeight));
}
- chunks.add(server.getMetadata().getPreferredTileHeight());
- chunks.add(server.getMetadata().getPreferredTileWidth());
+ chunks.add(chunkHeight);
+ chunks.add(chunkWidth);
int[] chunksArray = new int[chunks.size()];
for (int i = 0; i < chunks.size(); i++) {
@@ -317,18 +363,18 @@ private int[] getOffsetsOfTile(TileRequest tileRequest) {
return offsetArray;
}
- private Object getData(BufferedImage image, int height, int width) {
+ private Object getData(BufferedImage image) {
Object pixels = AWTImageTools.getPixels(image);
if (server.isRGB()) {
int[][] data = (int[][]) pixels;
- int[] output = new int[server.nChannels() * width * height];
+ int[] output = new int[server.nChannels() * image.getWidth() * image.getHeight()];
int i = 0;
for (int c=0; c {
byte[][] data = (byte[][]) pixels;
- byte[] output = new byte[server.nChannels() * width * height];
+ byte[] output = new byte[server.nChannels() * image.getWidth() * image.getHeight()];
int i = 0;
for (int c=0; c {
short[][] data = (short[][]) pixels;
- short[] output = new short[server.nChannels() * width * height];
+ short[] output = new short[server.nChannels() * image.getWidth() * image.getHeight()];
int i = 0;
for (int c=0; c {
int[][] data = (int[][]) pixels;
- int[] output = new int[server.nChannels() * width * height];
+ int[] output = new int[server.nChannels() * image.getWidth() * image.getHeight()];
int i = 0;
for (int c=0; c {
float[][] data = (float[][]) pixels;
- float[] output = new float[server.nChannels() * width * height];
+ float[] output = new float[server.nChannels() * image.getWidth() * image.getHeight()];
int i = 0;
for (int c=0; c {
double[][] data = (double[][]) pixels;
- double[] output = new double[server.nChannels() * width * height];
+ double[] output = new double[server.nChannels() * image.getWidth() * image.getHeight()];
int i = 0;
for (int c=0; c
Date: Fri, 8 Mar 2024 16:23:07 +0000
Subject: [PATCH 019/866] Improved max number of chunks
---
.../writers/ome/zarr/OMEZarrWriter.java | 55 +++++++++----------
1 file changed, 26 insertions(+), 29 deletions(-)
diff --git a/qupath-extension-bioformats/src/main/java/qupath/lib/images/writers/ome/zarr/OMEZarrWriter.java b/qupath-extension-bioformats/src/main/java/qupath/lib/images/writers/ome/zarr/OMEZarrWriter.java
index 781b3032b..c7c155232 100644
--- a/qupath-extension-bioformats/src/main/java/qupath/lib/images/writers/ome/zarr/OMEZarrWriter.java
+++ b/qupath-extension-bioformats/src/main/java/qupath/lib/images/writers/ome/zarr/OMEZarrWriter.java
@@ -15,7 +15,6 @@
import qupath.lib.images.servers.PixelCalibration;
import qupath.lib.images.servers.TileRequest;
import qupath.lib.images.servers.TileRequestManager;
-import ucar.ma2.InvalidRangeException;
import java.awt.image.BufferedImage;
import java.io.IOException;
@@ -77,7 +76,7 @@ private OMEZarrWriter(Builder builder) throws IOException {
attributes.getGroupAttributes()
),
attributes.getLevelAttributes(),
- builder.maxChunkSize
+ builder.maxNumberOfChunks
);
this.executorService = Executors.newFixedThreadPool(builder.numberOfThreads);
}
@@ -129,9 +128,10 @@ public void writeImage() {
public void writeTile(TileRequest tileRequest) {
executorService.execute(() -> {
try {
- writeTile(
- tileRequest,
- getData(server.readRegion(tileRequest.getRegionRequest()))
+ levelArrays.get(tileRequest.getLevel()).write(
+ getData(server.readRegion(tileRequest.getRegionRequest())),
+ getDimensionsOfTile(tileRequest),
+ getOffsetsOfTile(tileRequest)
);
} catch (Exception e) {
logger.error("Error when writing tile", e);
@@ -150,7 +150,7 @@ public static class Builder {
private Compressor compressor = CompressorFactory.createDefaultCompressor();
private int numberOfThreads = 12;
private double[] downsamples = new double[0];
- private int maxChunkSize = 12;
+ private int maxNumberOfChunks = 12;
/**
* Create the builder.
@@ -211,15 +211,20 @@ public Builder setDownsamples(double... downsamples) {
}
/**
- * In Zarr files, data is stored in chunks. This parameter
- * defines the maximum number of chunks on the x,y, and z dimensions.
- * By default, this value is set to 12.
*
- * @param maxChunkSize the maximum number of chunks on the x,y, and z dimensions
+ *
+ * In Zarr files, data is stored in chunks. This parameter defines the maximum number
+ * of chunks on the x,y, and z dimensions. By default, this value is set to 12.
+ *
+ *
+ * Use a negative value to not define any maximum number of chunks.
+ *
+ *
+ * @param maxNumberOfChunks the maximum number of chunks on the x,y, and z dimensions
* @return this builder
*/
- public Builder setMaxChunkSize(int maxChunkSize) {
- this.maxChunkSize = maxChunkSize;
+ public Builder setMaxNumberOfChunks(int maxNumberOfChunks) {
+ this.maxNumberOfChunks = maxNumberOfChunks;
return this;
}
@@ -236,7 +241,7 @@ public OMEZarrWriter build() throws IOException {
}
}
- private Map createLevelArrays(ZarrGroup root, Map levelAttributes, int maxChunkSize) throws IOException {
+ private Map createLevelArrays(ZarrGroup root, Map levelAttributes, int maxNumberOfChunks) throws IOException {
Map levelArrays = new HashMap<>();
for (int level=0; level createLevelArrays(ZarrGroup root, Map DataType.u1;
@@ -264,18 +269,6 @@ private Map createLevelArrays(ZarrGroup root, Map dimensions = new ArrayList<>();
if (server.nTimepoints() > 1) {
@@ -297,9 +290,13 @@ private int[] getDimensionsOfImage(int level) {
return dimensionArray;
}
- private int[] getChunksOfImage(int maxChunkSize) {
- int chunkWidth = Math.max(server.getMetadata().getPreferredTileWidth(), server.getWidth() / maxChunkSize);
- int chunkHeight = Math.max(server.getMetadata().getPreferredTileHeight(), server.getHeight() / maxChunkSize);
+ private int[] getChunksOfImage(int maxNumberOfChunks) {
+ int chunkWidth = maxNumberOfChunks > 0 ?
+ Math.max(server.getMetadata().getPreferredTileWidth(), server.getWidth() / maxNumberOfChunks) :
+ server.getMetadata().getPreferredTileWidth();
+ int chunkHeight = maxNumberOfChunks > 0 ?
+ Math.max(server.getMetadata().getPreferredTileHeight(), server.getHeight() / maxNumberOfChunks) :
+ server.getMetadata().getPreferredTileHeight();
List chunks = new ArrayList<>();
if (server.nTimepoints() > 1) {
From 768c05cc4143ad7dadc76e1674bba0b6baf699b2 Mon Sep 17 00:00:00 2001
From: lleplat
Date: Mon, 11 Mar 2024 14:09:56 +0000
Subject: [PATCH 020/866] Enforce max number of chunks
---
.../writers/ome/zarr/OMEZarrWriter.java | 77 ++++++++++++-------
1 file changed, 49 insertions(+), 28 deletions(-)
diff --git a/qupath-extension-bioformats/src/main/java/qupath/lib/images/writers/ome/zarr/OMEZarrWriter.java b/qupath-extension-bioformats/src/main/java/qupath/lib/images/writers/ome/zarr/OMEZarrWriter.java
index c7c155232..35e8745e7 100644
--- a/qupath-extension-bioformats/src/main/java/qupath/lib/images/writers/ome/zarr/OMEZarrWriter.java
+++ b/qupath-extension-bioformats/src/main/java/qupath/lib/images/writers/ome/zarr/OMEZarrWriter.java
@@ -45,18 +45,20 @@ public class OMEZarrWriter implements AutoCloseable {
private static final Logger logger = LoggerFactory.getLogger(OMEZarrWriter.class);
private final ImageServer server;
- private final Compressor compressor;
private final Map levelArrays;
private final ExecutorService executorService;
private OMEZarrWriter(Builder builder) throws IOException {
- this.server = builder.downsamples.length == 0 ?
- builder.server :
- ImageServers.pyramidalize(
- builder.server,
- builder.downsamples
- );
- this.compressor = builder.compressor;
+ server = ImageServers.pyramidalizeTiled(
+ builder.server,
+ builder.maxNumberOfChunks > 0 ?
+ Math.max(builder.server.getMetadata().getPreferredTileWidth(), builder.server.getWidth() / builder.maxNumberOfChunks) :
+ builder.server.getMetadata().getPreferredTileWidth(),
+ builder.maxNumberOfChunks > 0 ?
+ Math.max(builder.server.getMetadata().getPreferredTileHeight(), builder.server.getHeight() / builder.maxNumberOfChunks) :
+ builder.server.getMetadata().getPreferredTileHeight(),
+ builder.downsamples.length == 0 ? builder.server.getPreferredDownsamples() : builder.downsamples
+ );
OMEZarrAttributesCreator attributes = new OMEZarrAttributesCreator(
server.getMetadata().getName(),
@@ -70,15 +72,17 @@ private OMEZarrWriter(Builder builder) throws IOException {
server.isRGB(),
server.getPixelType()
);
- this.levelArrays = createLevelArrays(
+ levelArrays = createLevelArrays(
+ server,
ZarrGroup.create(
builder.path,
attributes.getGroupAttributes()
),
attributes.getLevelAttributes(),
- builder.maxNumberOfChunks
+ builder.compressor
);
- this.executorService = Executors.newFixedThreadPool(builder.numberOfThreads);
+
+ executorService = Executors.newFixedThreadPool(builder.numberOfThreads);
}
/**
@@ -121,6 +125,8 @@ public void writeImage() {
* {@link Builder#Builder(ImageServer, String)}. Therefore, the {@link ImageServer#getTileRequestManager() TileRequestManager}
* of the internal image server may be different from the one of the provided image server,
* so functions like {@link TileRequestManager#getAllTileRequests()} may not return the expected tiles.
+ * Use the {@link ImageServer#getTileRequestManager() TileRequestManager} of {@link #getReaderServer()}
+ * to get accurate tiles.
*
*
* @param tileRequest the tile to write
@@ -139,6 +145,23 @@ public void writeTile(TileRequest tileRequest) {
});
}
+ /**
+ *
+ *
+ * Get the image server used internally by this writer to read the tiles. It can be
+ * different from the one given in {@link Builder#Builder(ImageServer, String)}.
+ *
+ *
+ * This function can be useful to get information like the tiles used by this server
+ * (for example when using the {@link #writeTile(TileRequest)} function).
+ *
+ *
+ * @return the image server used internally by this writer to read the tiles
+ */
+ public ImageServer getReaderServer() {
+ return server;
+ }
+
/**
* Builder to create an instance of a {@link OMEZarrWriter}.
*/
@@ -150,7 +173,7 @@ public static class Builder {
private Compressor compressor = CompressorFactory.createDefaultCompressor();
private int numberOfThreads = 12;
private double[] downsamples = new double[0];
- private int maxNumberOfChunks = 12;
+ private int maxNumberOfChunks = 30;
/**
* Create the builder.
@@ -214,7 +237,7 @@ public Builder setDownsamples(double... downsamples) {
*
*
* In Zarr files, data is stored in chunks. This parameter defines the maximum number
- * of chunks on the x,y, and z dimensions. By default, this value is set to 12.
+ * of chunks on the x,y, and z dimensions. By default, this value is set to 30.
*
*
* Use a negative value to not define any maximum number of chunks.
@@ -241,15 +264,20 @@ public OMEZarrWriter build() throws IOException {
}
}
- private Map createLevelArrays(ZarrGroup root, Map levelAttributes, int maxNumberOfChunks) throws IOException {
+ private static Map createLevelArrays(
+ ImageServer server,
+ ZarrGroup root,
+ Map levelAttributes,
+ Compressor compressor
+ ) throws IOException {
Map levelArrays = new HashMap<>();
for (int level=0; level DataType.u1;
@@ -269,7 +297,7 @@ private Map createLevelArrays(ZarrGroup root, Map server, int level) {
List dimensions = new ArrayList<>();
if (server.nTimepoints() > 1) {
dimensions.add(server.nTimepoints());
@@ -290,14 +318,7 @@ private int[] getDimensionsOfImage(int level) {
return dimensionArray;
}
- private int[] getChunksOfImage(int maxNumberOfChunks) {
- int chunkWidth = maxNumberOfChunks > 0 ?
- Math.max(server.getMetadata().getPreferredTileWidth(), server.getWidth() / maxNumberOfChunks) :
- server.getMetadata().getPreferredTileWidth();
- int chunkHeight = maxNumberOfChunks > 0 ?
- Math.max(server.getMetadata().getPreferredTileHeight(), server.getHeight() / maxNumberOfChunks) :
- server.getMetadata().getPreferredTileHeight();
-
+ private static int[] getChunksOfImage(ImageServer server) {
List chunks = new ArrayList<>();
if (server.nTimepoints() > 1) {
chunks.add(1);
@@ -306,10 +327,10 @@ private int[] getChunksOfImage(int maxNumberOfChunks) {
chunks.add(1);
}
if (server.nZSlices() > 1) {
- chunks.add(Math.max(chunkWidth, chunkHeight));
+ chunks.add(Math.max(server.getMetadata().getPreferredTileWidth(), server.getMetadata().getPreferredTileHeight()));
}
- chunks.add(chunkHeight);
- chunks.add(chunkWidth);
+ chunks.add(server.getMetadata().getPreferredTileHeight());
+ chunks.add(server.getMetadata().getPreferredTileWidth());
int[] chunksArray = new int[chunks.size()];
for (int i = 0; i < chunks.size(); i++) {
From abe5c4e620e4e27a133d0ea985ed8b69d88a46ee Mon Sep 17 00:00:00 2001
From: lleplat
Date: Mon, 11 Mar 2024 14:44:47 +0000
Subject: [PATCH 021/866] Added tile size
---
.../writers/ome/zarr/OMEZarrWriter.java | 73 ++++++++++++++++---
1 file changed, 63 insertions(+), 10 deletions(-)
diff --git a/qupath-extension-bioformats/src/main/java/qupath/lib/images/writers/ome/zarr/OMEZarrWriter.java b/qupath-extension-bioformats/src/main/java/qupath/lib/images/writers/ome/zarr/OMEZarrWriter.java
index 35e8745e7..7c696a59f 100644
--- a/qupath-extension-bioformats/src/main/java/qupath/lib/images/writers/ome/zarr/OMEZarrWriter.java
+++ b/qupath-extension-bioformats/src/main/java/qupath/lib/images/writers/ome/zarr/OMEZarrWriter.java
@@ -51,12 +51,16 @@ public class OMEZarrWriter implements AutoCloseable {
private OMEZarrWriter(Builder builder) throws IOException {
server = ImageServers.pyramidalizeTiled(
builder.server,
- builder.maxNumberOfChunks > 0 ?
- Math.max(builder.server.getMetadata().getPreferredTileWidth(), builder.server.getWidth() / builder.maxNumberOfChunks) :
- builder.server.getMetadata().getPreferredTileWidth(),
- builder.maxNumberOfChunks > 0 ?
- Math.max(builder.server.getMetadata().getPreferredTileHeight(), builder.server.getHeight() / builder.maxNumberOfChunks) :
- builder.server.getMetadata().getPreferredTileHeight(),
+ getChunkSize(
+ builder.tileWidth > 0 ? builder.tileWidth : builder.server.getMetadata().getPreferredTileWidth(),
+ builder.maxNumberOfChunks,
+ builder.server.getWidth()
+ ),
+ getChunkSize(
+ builder.tileHeight > 0 ? builder.tileHeight : builder.server.getMetadata().getPreferredTileHeight(),
+ builder.maxNumberOfChunks,
+ builder.server.getHeight()
+ ),
builder.downsamples.length == 0 ? builder.server.getPreferredDownsamples() : builder.downsamples
);
@@ -173,7 +177,9 @@ public static class Builder {
private Compressor compressor = CompressorFactory.createDefaultCompressor();
private int numberOfThreads = 12;
private double[] downsamples = new double[0];
- private int maxNumberOfChunks = 30;
+ private int maxNumberOfChunks = 50;
+ private int tileWidth = 512;
+ private int tileHeight = 512;
/**
* Create the builder.
@@ -234,10 +240,9 @@ public Builder setDownsamples(double... downsamples) {
}
/**
- *
*
* In Zarr files, data is stored in chunks. This parameter defines the maximum number
- * of chunks on the x,y, and z dimensions. By default, this value is set to 30.
+ * of chunks on the x,y, and z dimensions. By default, this value is set to 50.
*
*
* Use a negative value to not define any maximum number of chunks.
@@ -246,11 +251,53 @@ public Builder setDownsamples(double... downsamples) {
* @param maxNumberOfChunks the maximum number of chunks on the x,y, and z dimensions
* @return this builder
*/
- public Builder setMaxNumberOfChunks(int maxNumberOfChunks) {
+ public Builder setMaxNumberOfChunksOnEachSpatialDimension(int maxNumberOfChunks) {
this.maxNumberOfChunks = maxNumberOfChunks;
return this;
}
+ /**
+ *
+ * In Zarr files, data is stored in chunks. This parameter defines the size
+ * of chunks on the x dimension. By default, this value is set to 512.
+ *
+ *
+ * Use a negative value to use the tile width of the provided image server.
+ *
+ *
+ * The provided tile width may not be used if this implies creating more chunks
+ * than the value given in {@link #setMaxNumberOfChunksOnEachSpatialDimension(int)}.
+ *
+ *
+ * @param tileWidth the width each chunk should have
+ * @return this builder
+ */
+ public Builder setTileWidth(int tileWidth) {
+ this.tileWidth = tileWidth;
+ return this;
+ }
+
+ /**
+ *
+ * In Zarr files, data is stored in chunks. This parameter defines the size
+ * of chunks on the x dimension. By default, this value is set to 512.
+ *
+ *
+ * Use a negative value to use the tile width of the provided image server.
+ *
+ *
+ * The provided tile width may not be used if this implies creating more chunks
+ * than the value given in {@link #setMaxNumberOfChunksOnEachSpatialDimension(int)}.
+ *
+ *
+ * @param tileHeight the height each chunk should have
+ * @return this builder
+ */
+ public Builder setTileHeight(int tileHeight) {
+ this.tileHeight = tileHeight;
+ return this;
+ }
+
/**
* Create a new instance of {@link OMEZarrWriter}. This will also
* create an empty image on the provided path.
@@ -264,6 +311,12 @@ public OMEZarrWriter build() throws IOException {
}
}
+ private static int getChunkSize(int tileSize, int maxNumberOfChunks, int imageSize) {
+ return maxNumberOfChunks > 0 ?
+ Math.max(tileSize, imageSize / maxNumberOfChunks) :
+ tileSize;
+ }
+
private static Map createLevelArrays(
ImageServer server,
ZarrGroup root,
From e8f96669038246bbee90f63d9b194844f7bdf3f2 Mon Sep 17 00:00:00 2001
From: Leo Leplat <60394504+Rylern@users.noreply.github.com>
Date: Tue, 12 Mar 2024 11:16:04 +0000
Subject: [PATCH 022/866] Typo
---
.../writers/ome/zarr/OMEZarrWriter.java | 90 +++++++++----------
1 file changed, 45 insertions(+), 45 deletions(-)
diff --git a/qupath-extension-bioformats/src/main/java/qupath/lib/images/writers/ome/zarr/OMEZarrWriter.java b/qupath-extension-bioformats/src/main/java/qupath/lib/images/writers/ome/zarr/OMEZarrWriter.java
index 7c696a59f..122a5aef4 100644
--- a/qupath-extension-bioformats/src/main/java/qupath/lib/images/writers/ome/zarr/OMEZarrWriter.java
+++ b/qupath-extension-bioformats/src/main/java/qupath/lib/images/writers/ome/zarr/OMEZarrWriter.java
@@ -280,13 +280,13 @@ public Builder setTileWidth(int tileWidth) {
/**
*
* In Zarr files, data is stored in chunks. This parameter defines the size
- * of chunks on the x dimension. By default, this value is set to 512.
+ * of chunks on the y dimension. By default, this value is set to 512.
*
*
- * Use a negative value to use the tile width of the provided image server.
+ * Use a negative value to use the tile height of the provided image server.
*
*
- * The provided tile width may not be used if this implies creating more chunks
+ * The provided tile height may not be used if this implies creating more chunks
* than the value given in {@link #setMaxNumberOfChunksOnEachSpatialDimension(int)}.
*
*
@@ -392,48 +392,6 @@ private static int[] getChunksOfImage(ImageServer server) {
return chunksArray;
}
- private int[] getDimensionsOfTile(TileRequest tileRequest) {
- List dimensions = new ArrayList<>();
- if (server.nTimepoints() > 1) {
- dimensions.add(1);
- }
- if (server.nChannels() > 1) {
- dimensions.add(server.nChannels());
- }
- if (server.nZSlices() > 1) {
- dimensions.add(1);
- }
- dimensions.add(tileRequest.getTileHeight());
- dimensions.add(tileRequest.getTileWidth());
-
- int[] dimensionArray = new int[dimensions.size()];
- for (int i = 0; i < dimensions.size(); i++) {
- dimensionArray[i] = dimensions.get(i);
- }
- return dimensionArray;
- }
-
- private int[] getOffsetsOfTile(TileRequest tileRequest) {
- List offset = new ArrayList<>();
- if (server.nTimepoints() > 1) {
- offset.add(tileRequest.getT());
- }
- if (server.nChannels() > 1) {
- offset.add(0);
- }
- if (server.nZSlices() > 1) {
- offset.add(tileRequest.getZ());
- }
- offset.add(tileRequest.getTileY());
- offset.add(tileRequest.getTileX());
-
- int[] offsetArray = new int[offset.size()];
- for (int i = 0; i < offset.size(); i++) {
- offsetArray[i] = offset.get(i);
- }
- return offsetArray;
- }
-
private Object getData(BufferedImage image) {
Object pixels = AWTImageTools.getPixels(image);
@@ -531,4 +489,46 @@ private Object getData(BufferedImage image) {
};
}
}
+
+ private int[] getDimensionsOfTile(TileRequest tileRequest) {
+ List dimensions = new ArrayList<>();
+ if (server.nTimepoints() > 1) {
+ dimensions.add(1);
+ }
+ if (server.nChannels() > 1) {
+ dimensions.add(server.nChannels());
+ }
+ if (server.nZSlices() > 1) {
+ dimensions.add(1);
+ }
+ dimensions.add(tileRequest.getTileHeight());
+ dimensions.add(tileRequest.getTileWidth());
+
+ int[] dimensionArray = new int[dimensions.size()];
+ for (int i = 0; i < dimensions.size(); i++) {
+ dimensionArray[i] = dimensions.get(i);
+ }
+ return dimensionArray;
+ }
+
+ private int[] getOffsetsOfTile(TileRequest tileRequest) {
+ List offset = new ArrayList<>();
+ if (server.nTimepoints() > 1) {
+ offset.add(tileRequest.getT());
+ }
+ if (server.nChannels() > 1) {
+ offset.add(0);
+ }
+ if (server.nZSlices() > 1) {
+ offset.add(tileRequest.getZ());
+ }
+ offset.add(tileRequest.getTileY());
+ offset.add(tileRequest.getTileX());
+
+ int[] offsetArray = new int[offset.size()];
+ for (int i = 0; i < offset.size(); i++) {
+ offsetArray[i] = offset.get(i);
+ }
+ return offsetArray;
+ }
}
From 09659c68e9115d0f9af1faa00fd468ced21ef1e2 Mon Sep 17 00:00:00 2001
From: Pete
Date: Thu, 28 Mar 2024 12:08:17 +0000
Subject: [PATCH 023/866] Lazy-load ImageServer from ImageData
Rough draft at an alternative approach to address the issue raised in https://github.com/qupath/qupath/pull/1488
---
.../java/qupath/lib/images/ImageData.java | 72 +++++++++++++++----
.../src/main/java/qupath/lib/io/PathIO.java | 69 +++++++++++++-----
.../qupath/lib/projects/DefaultProject.java | 38 +++++-----
.../gui/scripting/DefaultScriptEditor.java | 4 +-
4 files changed, 131 insertions(+), 52 deletions(-)
diff --git a/qupath-core/src/main/java/qupath/lib/images/ImageData.java b/qupath-core/src/main/java/qupath/lib/images/ImageData.java
index 824dceee4..c612906b3 100644
--- a/qupath-core/src/main/java/qupath/lib/images/ImageData.java
+++ b/qupath-core/src/main/java/qupath/lib/images/ImageData.java
@@ -4,7 +4,7 @@
* %%
* Copyright (C) 2014 - 2016 The Queen's University of Belfast, Northern Ireland
* Contact: IP Management (ipmanagement@qub.ac.uk)
- * Copyright (C) 2018 - 2020 QuPath developers, The University of Edinburgh
+ * Copyright (C) 2018 - 2024 QuPath developers, The University of Edinburgh
* %%
* QuPath is free software: you can redistribute it and/or modify
* it under the terms of the GNU General Public License as
@@ -29,6 +29,7 @@
import java.util.HashMap;
import java.util.Map;
import java.util.Objects;
+import java.util.function.Supplier;
import org.slf4j.Logger;
import org.slf4j.LoggerFactory;
@@ -60,8 +61,8 @@
* @param
*
*/
-public class ImageData implements WorkflowListener, PathObjectHierarchyListener {
-
+public class ImageData implements WorkflowListener, PathObjectHierarchyListener, AutoCloseable {
+
/**
* Enum representing possible image types.
*
@@ -109,7 +110,10 @@ public String toString() {
private static final Logger logger = LoggerFactory.getLogger(ImageData.class);
private transient PropertyChangeSupport pcs;
-
+
+ private transient Supplier> serverSupplier;
+ private transient ImageServerMetadata lazyMetadata;
+
private transient ImageServer server;
private String lastSavedPath = null;
@@ -138,8 +142,9 @@ public String toString() {
* @param hierarchy
* @param type
*/
- public ImageData(ImageServer server, PathObjectHierarchy hierarchy, ImageType type) {
- pcs = new PropertyChangeSupport(this);
+ public ImageData(Supplier> supplier, ImageServer server, PathObjectHierarchy hierarchy, ImageType type) {
+ this.pcs = new PropertyChangeSupport(this);
+ this.serverSupplier = supplier;
this.server = server;
this.hierarchy = hierarchy == null ? new PathObjectHierarchy() : hierarchy;
this.serverPath = server == null ? null : server.getPath(); // TODO: Deal with sub image servers
@@ -147,14 +152,18 @@ public ImageData(ImageServer server, PathObjectHierarchy hierarchy, ImageType
if (type == null)
type = ImageType.UNSET;
setImageType(type);
-
+
// Add listeners for changes
this.hierarchy.addListener(this);
workflow.addWorkflowListener(this);
-
+
// Discard any changes during construction
changes = false;
}
+
+ public ImageData(ImageServer server, PathObjectHierarchy hierarchy, ImageType type) {
+ this(null, server, hierarchy, type);
+ }
/**
* Create a new ImageData with a specified type and creating a new object hierarchy.
@@ -231,6 +240,15 @@ public void setColorDeconvolutionStains(ColorDeconvolutionStains stains) {
*/
public void updateServerMetadata(ImageServerMetadata newMetadata) {
Objects.requireNonNull(newMetadata);
+ if (server == null) {
+ if (serverSupplier == null)
+ throw new IllegalStateException("Cannot update server metadata without a server or server supplier");
+ else {
+ logger.debug("Setting serve metadata lazily (no change will be fired)");
+ lazyMetadata = newMetadata;
+ return;
+ }
+ }
logger.trace("Updating server metadata");
var oldMetadata = server.getMetadata();
server.setMetadata(newMetadata);
@@ -345,6 +363,13 @@ private static void addColorDeconvolutionStainsToWorkflow(ImageData> imageData
* @return
*/
public ImageServer getServer() {
+ if (server == null && serverSupplier != null) {
+ logger.debug("Lazily requesting image server");
+ server = serverSupplier.get();
+ if (lazyMetadata != null && !lazyMetadata.equals(server.getMetadata())) {
+ updateServerMetadata(lazyMetadata);
+ }
+ }
return server;
}
@@ -506,14 +531,33 @@ public void hierarchyChanged(PathObjectHierarchyEvent event) {
public void workflowUpdated(Workflow workflow) {
changes = true;
}
-
-
+
+ /**
+ * Close the server if it has been loaded.
+ * @throws Exception
+ */
+ @Override
+ public void close() throws Exception {
+ if (server != null)
+ server.close();
+ }
+
+
@Override
public String toString() {
- if (getServer() == null)
- return "ImageData: " + getImageType() + ", no server";
- else
- return "ImageData: " + getImageType() + ", " + ServerTools.getDisplayableImageName(getServer());
+ String serverName;
+ if (server == null) {
+ if (serverSupplier == null) {
+ serverName = "no server";
+ } else if (lazyMetadata != null){
+ serverName = lazyMetadata.getName() + " (not yet loaded)";
+ } else {
+ serverName = "lazy-loaded server";
+ }
+ } else {
+ serverName = ServerTools.getDisplayableImageName(server);
+ }
+ return "ImageData: " + getImageType() + ", " + serverName;
}
diff --git a/qupath-core/src/main/java/qupath/lib/io/PathIO.java b/qupath-core/src/main/java/qupath/lib/io/PathIO.java
index 36c42c297..29f194935 100644
--- a/qupath-core/src/main/java/qupath/lib/io/PathIO.java
+++ b/qupath-core/src/main/java/qupath/lib/io/PathIO.java
@@ -4,7 +4,7 @@
* %%
* Copyright (C) 2014 - 2016 The Queen's University of Belfast, Northern Ireland
* Contact: IP Management (ipmanagement@qub.ac.uk)
- * Copyright (C) 2018 - 2020 QuPath developers, The University of Edinburgh
+ * Copyright (C) 2018 - 2024 QuPath developers, The University of Edinburgh
* %%
* QuPath is free software: you can redistribute it and/or modify
* it under the terms of the GNU General Public License as
@@ -78,6 +78,7 @@
import java.util.Map.Entry;
import java.util.Objects;
import java.util.Set;
+import java.util.function.Supplier;
import java.util.stream.Collectors;
import java.util.zip.GZIPInputStream;
import java.util.zip.GZIPOutputStream;
@@ -287,15 +288,29 @@ static ServerBuilderWrapper create(ServerBuilder builder, String id)
}
}
-
-
-
- private static ImageData readImageDataSerialized(final Path path, ImageData imageData, ImageServer server, Class cls) throws FileNotFoundException, IOException {
+
+
+ private static ImageData readImageDataSerialized(final Path path, ImageData imageData,
+ ImageServer server, Class cls) throws FileNotFoundException, IOException {
+ imageData = readImageDataSerialized(path, imageData, () -> server, cls);
+ imageData.getServer(); // Ensure the server is loaded
+ return imageData;
+ }
+
+ private static ImageData readImageDataSerialized(final InputStream stream, ImageData imageData,
+ ImageServer server, Class cls) throws FileNotFoundException, IOException {
+ imageData = readImageDataSerialized(stream, imageData, () -> server, cls);
+ imageData.getServer(); // Ensure the server is loaded
+ return imageData;
+ }
+
+ private static ImageData readImageDataSerialized(final Path path, ImageData imageData,
+ Supplier> serverSupplier, Class cls) throws FileNotFoundException, IOException {
if (path == null)
return null;
logger.info("Reading data from {}...", path.getFileName().toString());
try (InputStream stream = Files.newInputStream(path)) {
- imageData = readImageDataSerialized(stream, imageData, server, cls);
+ imageData = readImageDataSerialized(stream, imageData, serverSupplier, cls);
// Set the last saved path (actually the path from which this was opened)
if (imageData != null)
imageData.setLastSavedPath(path.toAbsolutePath().toString(), true);
@@ -307,7 +322,7 @@ private static ImageData readImageDataSerialized(final Path path, ImageDa
}
@SuppressWarnings("unchecked")
- private static ImageData readImageDataSerialized(final InputStream stream, ImageData imageData, ImageServer server, Class cls) throws IOException {
+ private static ImageData readImageDataSerialized(final InputStream stream, ImageData imageData, Supplier> serverSupplier, Class cls) throws IOException {
long startTime = System.currentTimeMillis();
Locale locale = Locale.getDefault(Category.FORMAT);
@@ -338,7 +353,7 @@ private static ImageData readImageDataSerialized(final InputStream stream
String serverString = (String)inStream.readObject();
// Don't log warnings if we are provided with a server
- serverBuilder = extractServerBuilder(serverString, server == null);
+ serverBuilder = extractServerBuilder(serverString, serverSupplier == null);
while (true) {
// logger.debug("Starting read: " + inStream.available());
@@ -395,23 +410,27 @@ else if (input == null) {
var existingBuilder = imageData == null || imageData.getServer() == null ? null : imageData.getServer().getBuilder();
if (imageData == null || !Objects.equals(serverBuilder, existingBuilder)) {
// Create a new server if we need to
- if (server == null) {
- try {
- server = serverBuilder.build();
- } catch (Exception e) {
- logger.error(e.getLocalizedMessage());
+ if (serverSupplier == null) {
+ // Load the server lazily
+ var builder = serverBuilder;
+ serverSupplier = () -> {
+ try {
+ return builder.build();
+ } catch (Exception e) {
+ logger.error("Warning: Unable to build server with " + builder);
+ if (e instanceof RuntimeException runtimeException)
+ throw runtimeException;
+ else
+ throw new RuntimeException(e);
+ }
};
- if (server == null) {
- logger.error("Warning: Unable to build server with " + serverBuilder);
- // throw new RuntimeException("Warning: Unable to create server for path " + serverPath);
- }
}
// TODO: Make this less clumsy... but for now we need to ensure we have a fully-initialized hierarchy (which deserialization alone doesn't achieve)
PathObjectHierarchy hierarchy2 = new PathObjectHierarchy();
hierarchy2.setHierarchy(hierarchy);
hierarchy = hierarchy2;
- imageData = new ImageData<>(server, hierarchy, imageType);
+ imageData = new ImageData<>(serverSupplier, null, hierarchy, imageType);
} else {
if (imageType != null)
imageData.setImageType(imageType);
@@ -487,6 +506,20 @@ public static ImageData readImageData(final InputStream stream, ImageData
return readImageDataSerialized(stream, imageData, server, cls);
}
+
+ /**
+ * Read an ImageData with lazy image loading.
+ * @param stream
+ * @param serverSupplier
+ * @param cls
+ * @return
+ * @param
+ * @throws IOException
+ */
+ public static ImageData readLazyImageData(final InputStream stream, Supplier> serverSupplier, Class cls) throws IOException {
+ return readImageDataSerialized(stream, null, serverSupplier, cls);
+ }
+
/**
* Read ImageData from a File into an existing ImageData object, or create a new one if required.
diff --git a/qupath-core/src/main/java/qupath/lib/projects/DefaultProject.java b/qupath-core/src/main/java/qupath/lib/projects/DefaultProject.java
index 0523f07aa..04889a001 100644
--- a/qupath-core/src/main/java/qupath/lib/projects/DefaultProject.java
+++ b/qupath-core/src/main/java/qupath/lib/projects/DefaultProject.java
@@ -4,7 +4,7 @@
* %%
* Copyright (C) 2014 - 2016 The Queen's University of Belfast, Northern Ireland
* Contact: IP Management (ipmanagement@qub.ac.uk)
- * Copyright (C) 2018 - 2020 QuPath developers, The University of Edinburgh
+ * Copyright (C) 2018 - 2024 QuPath developers, The University of Edinburgh
* %%
* QuPath is free software: you can redistribute it and/or modify
* it under the terms of the GNU General Public License as
@@ -47,6 +47,7 @@
import java.util.Map.Entry;
import java.util.Objects;
import java.util.concurrent.atomic.AtomicLong;
+import java.util.function.Supplier;
import java.util.stream.Collectors;
import javax.imageio.ImageIO;
@@ -688,20 +689,25 @@ private Path getThumbnailPath() {
@Override
public synchronized ImageData readImageData() throws IOException {
Path path = getImageDataPath();
- ImageServer server;
- try {
- server = getServerBuilder().build();
- } catch (IOException e) {
- throw e;
- } catch (Exception e) {
- throw new IOException(e);
- }
- if (server == null)
- return null;
+ Supplier> serverSupplier = () -> {
+ try {
+ var server = getServerBuilder().build();
+ // Ensure the names match
+ var name = getOriginalImageName();
+ if (name != null)
+ ServerTools.setImageName(server, name);
+ return server;
+ } catch (Exception e) {
+ if (e instanceof RuntimeException exception)
+ throw exception;
+ else
+ throw new RuntimeException(e);
+ }
+ };
ImageData imageData = null;
if (Files.exists(path)) {
try (var stream = Files.newInputStream(path)) {
- imageData = PathIO.readImageData(stream, null, server, BufferedImage.class);
+ imageData = PathIO.readLazyImageData(stream, serverSupplier, BufferedImage.class);
imageData.setLastSavedPath(path.toString(), true);
} catch (Exception e) {
logger.error("Error reading image data from " + path, e);
@@ -713,7 +719,7 @@ public synchronized ImageData readImageData() throws IOException
var pathBackup = getBackupImageDataPath();
if (Files.exists(pathBackup)) {
try (var stream = Files.newInputStream(pathBackup)) {
- imageData = PathIO.readImageData(stream, null, server, BufferedImage.class);
+ imageData = PathIO.readLazyImageData(stream, serverSupplier, BufferedImage.class);
imageData.setLastSavedPath(pathBackup.toString(), true);
logger.warn("Restored previous ImageData from {}", pathBackup);
} catch (IOException e) {
@@ -723,11 +729,7 @@ public synchronized ImageData readImageData() throws IOException
}
if (imageData == null)
- imageData = new ImageData<>(server);
- // Ensure the names match
- var name = getOriginalImageName();
- if (name != null)
- ServerTools.setImageName(server, name);
+ imageData = new ImageData<>(serverSupplier, null, new PathObjectHierarchy(), ImageType.UNSET);
imageData.setProperty(IMAGE_ID, getFullProjectEntryID()); // Required to be able to test for the ID later
imageData.setChanged(false);
return imageData;
diff --git a/qupath-gui-fx/src/main/java/qupath/lib/gui/scripting/DefaultScriptEditor.java b/qupath-gui-fx/src/main/java/qupath/lib/gui/scripting/DefaultScriptEditor.java
index e7b7be578..4735744ac 100644
--- a/qupath-gui-fx/src/main/java/qupath/lib/gui/scripting/DefaultScriptEditor.java
+++ b/qupath-gui-fx/src/main/java/qupath/lib/gui/scripting/DefaultScriptEditor.java
@@ -4,7 +4,7 @@
* %%
* Copyright (C) 2014 - 2016 The Queen's University of Belfast, Northern Ireland
* Contact: IP Management (ipmanagement@qub.ac.uk)
- * Copyright (C) 2018 - 2022 QuPath developers, The University of Edinburgh
+ * Copyright (C) 2018 - 2024 QuPath developers, The University of Edinburgh
* %%
* QuPath is free software: you can redistribute it and/or modify
* it under the terms of the GNU General Public License as
@@ -1730,7 +1730,7 @@ public Void call() {
executeScript(tab, tab.getEditorControl().getText(), project, imageData, batchIndex, batchSize, doSave, useCompiled);
if (doSave)
entry.saveImageData(imageData);
- imageData.getServer().close();
+ imageData.close();
if (clearCache.get()) {
try {
From aae4fdf8453d40340b01cbd3daf74a47ed8bc736 Mon Sep 17 00:00:00 2001
From: Pete
Date: Thu, 28 Mar 2024 12:19:28 +0000
Subject: [PATCH 024/866] Synchronize lazy loading in ImageData
---
.../src/main/java/qupath/lib/images/ImageData.java | 12 ++++++++----
1 file changed, 8 insertions(+), 4 deletions(-)
diff --git a/qupath-core/src/main/java/qupath/lib/images/ImageData.java b/qupath-core/src/main/java/qupath/lib/images/ImageData.java
index c612906b3..29f447f44 100644
--- a/qupath-core/src/main/java/qupath/lib/images/ImageData.java
+++ b/qupath-core/src/main/java/qupath/lib/images/ImageData.java
@@ -364,10 +364,14 @@ private static void addColorDeconvolutionStainsToWorkflow(ImageData> imageData
*/
public ImageServer getServer() {
if (server == null && serverSupplier != null) {
- logger.debug("Lazily requesting image server");
- server = serverSupplier.get();
- if (lazyMetadata != null && !lazyMetadata.equals(server.getMetadata())) {
- updateServerMetadata(lazyMetadata);
+ synchronized (this) {
+ if (server == null) {
+ logger.debug("Lazily requesting image server");
+ server = serverSupplier.get();
+ if (lazyMetadata != null && !lazyMetadata.equals(server.getMetadata())) {
+ updateServerMetadata(lazyMetadata);
+ }
+ }
}
}
return server;
From 1a559eb93565294a9f7df100a7e0d3ffd29b4f47 Mon Sep 17 00:00:00 2001
From: Pete
Date: Fri, 12 Apr 2024 12:31:28 +0100
Subject: [PATCH 025/866] Bump v0.6.0-SNAPSHOT
---
qupath-app/src/main/resources/VERSION | 2 +-
qupath-core-processing/src/main/resources/VERSION | 2 +-
qupath-core/src/main/resources/VERSION | 2 +-
qupath-extension-bioformats/src/main/resources/VERSION | 2 +-
qupath-extension-openslide/src/main/resources/VERSION | 2 +-
qupath-extension-processing/src/main/resources/VERSION | 2 +-
qupath-extension-script-editor/src/main/resources/VERSION | 2 +-
qupath-extension-svg/src/main/resources/VERSION | 2 +-
qupath-gui-fx/src/main/resources/VERSION | 2 +-
settings.gradle | 3 +--
10 files changed, 10 insertions(+), 11 deletions(-)
diff --git a/qupath-app/src/main/resources/VERSION b/qupath-app/src/main/resources/VERSION
index 5d4294b91..e75e24346 100644
--- a/qupath-app/src/main/resources/VERSION
+++ b/qupath-app/src/main/resources/VERSION
@@ -1 +1 @@
-0.5.1
\ No newline at end of file
+0.6.0-SNAPSHOT
\ No newline at end of file
diff --git a/qupath-core-processing/src/main/resources/VERSION b/qupath-core-processing/src/main/resources/VERSION
index 5d4294b91..e75e24346 100644
--- a/qupath-core-processing/src/main/resources/VERSION
+++ b/qupath-core-processing/src/main/resources/VERSION
@@ -1 +1 @@
-0.5.1
\ No newline at end of file
+0.6.0-SNAPSHOT
\ No newline at end of file
diff --git a/qupath-core/src/main/resources/VERSION b/qupath-core/src/main/resources/VERSION
index 5d4294b91..e75e24346 100644
--- a/qupath-core/src/main/resources/VERSION
+++ b/qupath-core/src/main/resources/VERSION
@@ -1 +1 @@
-0.5.1
\ No newline at end of file
+0.6.0-SNAPSHOT
\ No newline at end of file
diff --git a/qupath-extension-bioformats/src/main/resources/VERSION b/qupath-extension-bioformats/src/main/resources/VERSION
index 5d4294b91..e75e24346 100644
--- a/qupath-extension-bioformats/src/main/resources/VERSION
+++ b/qupath-extension-bioformats/src/main/resources/VERSION
@@ -1 +1 @@
-0.5.1
\ No newline at end of file
+0.6.0-SNAPSHOT
\ No newline at end of file
diff --git a/qupath-extension-openslide/src/main/resources/VERSION b/qupath-extension-openslide/src/main/resources/VERSION
index 5d4294b91..e75e24346 100644
--- a/qupath-extension-openslide/src/main/resources/VERSION
+++ b/qupath-extension-openslide/src/main/resources/VERSION
@@ -1 +1 @@
-0.5.1
\ No newline at end of file
+0.6.0-SNAPSHOT
\ No newline at end of file
diff --git a/qupath-extension-processing/src/main/resources/VERSION b/qupath-extension-processing/src/main/resources/VERSION
index 5d4294b91..e75e24346 100644
--- a/qupath-extension-processing/src/main/resources/VERSION
+++ b/qupath-extension-processing/src/main/resources/VERSION
@@ -1 +1 @@
-0.5.1
\ No newline at end of file
+0.6.0-SNAPSHOT
\ No newline at end of file
diff --git a/qupath-extension-script-editor/src/main/resources/VERSION b/qupath-extension-script-editor/src/main/resources/VERSION
index 5d4294b91..e75e24346 100644
--- a/qupath-extension-script-editor/src/main/resources/VERSION
+++ b/qupath-extension-script-editor/src/main/resources/VERSION
@@ -1 +1 @@
-0.5.1
\ No newline at end of file
+0.6.0-SNAPSHOT
\ No newline at end of file
diff --git a/qupath-extension-svg/src/main/resources/VERSION b/qupath-extension-svg/src/main/resources/VERSION
index 5d4294b91..e75e24346 100644
--- a/qupath-extension-svg/src/main/resources/VERSION
+++ b/qupath-extension-svg/src/main/resources/VERSION
@@ -1 +1 @@
-0.5.1
\ No newline at end of file
+0.6.0-SNAPSHOT
\ No newline at end of file
diff --git a/qupath-gui-fx/src/main/resources/VERSION b/qupath-gui-fx/src/main/resources/VERSION
index 5d4294b91..e75e24346 100644
--- a/qupath-gui-fx/src/main/resources/VERSION
+++ b/qupath-gui-fx/src/main/resources/VERSION
@@ -1 +1 @@
-0.5.1
\ No newline at end of file
+0.6.0-SNAPSHOT
\ No newline at end of file
diff --git a/settings.gradle b/settings.gradle
index de2f1d85e..5d43ef8f6 100644
--- a/settings.gradle
+++ b/settings.gradle
@@ -2,7 +2,7 @@ plugins {
id 'org.gradle.toolchains.foojay-resolver-convention' version '0.7.0' // to download if needed
}
-gradle.ext.qupathVersion = "0.5.1"
+gradle.ext.qupathVersion = "0.6.0-SNAPSHOT"
rootProject.name = 'qupath'
@@ -36,4 +36,3 @@ dependencyResolutionManagement {
}
}
}
-
From 2e4713c634bb454390d2d8c444c3ebb71dfaa107 Mon Sep 17 00:00:00 2001
From: Pete
Date: Fri, 12 Apr 2024 12:38:45 +0100
Subject: [PATCH 026/866] Fix WrappedBufferedImageServer bug
Fixes https://github.com/qupath/qupath/issues/1494
---
CHANGELOG.md | 7 +++++++
.../lib/images/servers/WrappedBufferedImageServer.java | 6 ++++--
2 files changed, 11 insertions(+), 2 deletions(-)
diff --git a/CHANGELOG.md b/CHANGELOG.md
index 56fd9deac..0341c5a87 100644
--- a/CHANGELOG.md
+++ b/CHANGELOG.md
@@ -1,3 +1,10 @@
+## Version 0.6.0-SNAPSHOT
+
+This is a *work in progress* for the next major release.
+
+### Bugs fixed
+* Tile export to .ome.tif can convert to 8-bit unnecessarily (https://github.com/qupath/qupath/issues/1494)
+
## Version 0.5.1
This is a *minor release* that aims to be fully compatible with v0.5.0, while fixing several bugs.
diff --git a/qupath-core/src/main/java/qupath/lib/images/servers/WrappedBufferedImageServer.java b/qupath-core/src/main/java/qupath/lib/images/servers/WrappedBufferedImageServer.java
index b015eabdf..23e2c03e8 100644
--- a/qupath-core/src/main/java/qupath/lib/images/servers/WrappedBufferedImageServer.java
+++ b/qupath-core/src/main/java/qupath/lib/images/servers/WrappedBufferedImageServer.java
@@ -101,8 +101,10 @@ public WrappedBufferedImageServer(final String imageName, final BufferedImage im
int nChannels = img.getSampleModel().getNumBands();
boolean isRGB = false;
for (int type : rgbTypes) {
- isRGB = isRGB | type == img.getType();
- pixelType = PixelType.UINT8;
+ if (type == img.getType()) {
+ isRGB = true;
+ pixelType = PixelType.UINT8;
+ }
}
// Warning! This method of obtaining channels risks resulting in different colors from the original image
if (channels == null) {
From 3c437522f4a5fc07b0e46686b8c067e5626f7192 Mon Sep 17 00:00:00 2001
From: Pete
Date: Fri, 12 Apr 2024 13:18:51 +0100
Subject: [PATCH 027/866] Update dependencies
These will likely cause some problems:
* Any JavaFX from 20.0.1 onwards breaks (javascript) search links within javadocs - at least when accessing javadocs from within a .jar
* Deep Java Library has probably broken the approach previously used to test if an `Engine` is available before attempting to download it
Nevertheless, we don't want to be stuck on old dependencies - so merging early in the path towards v0.6.0 will hopefully give us time to find fixes/workarounds.
This commit also updates the pref name to have the correct QuPath version.
---
CHANGELOG.md | 14 +++++
gradle/libs.versions.toml | 48 ++++++++----------
gradle/wrapper/gradle-wrapper.jar | Bin 63721 -> 43462 bytes
gradle/wrapper/gradle-wrapper.properties | 2 +-
.../java/qupath/opencv/dnn/OpenCVDnn.java | 2 +-
.../java/qupath/lib/gui/prefs/PathPrefs.java | 4 +-
6 files changed, 38 insertions(+), 32 deletions(-)
diff --git a/CHANGELOG.md b/CHANGELOG.md
index 0341c5a87..ba80f6ebe 100644
--- a/CHANGELOG.md
+++ b/CHANGELOG.md
@@ -5,6 +5,20 @@ This is a *work in progress* for the next major release.
### Bugs fixed
* Tile export to .ome.tif can convert to 8-bit unnecessarily (https://github.com/qupath/qupath/issues/1494)
+### Dependency updates
+* Bio-Formats 7.2.0
+* Commonmark 0.22.0
+* DeepJavaLibrary 0.27.0
+* Groovy 4.0.21
+* Guava 33.1.0-jre
+* ImageJ 1.54i
+* JavaFX 22
+* JNA 5.14.0
+* OpenCV 4.9.0
+* OpenJDK 21
+* slf4j 2.0.12
+
+
## Version 0.5.1
This is a *minor release* that aims to be fully compatible with v0.5.0, while fixing several bugs.
diff --git a/gradle/libs.versions.toml b/gradle/libs.versions.toml
index b3ef2b473..e6e613976 100644
--- a/gradle/libs.versions.toml
+++ b/gradle/libs.versions.toml
@@ -1,42 +1,34 @@
[versions]
-bioformats = "7.0.1"
+bioformats = "7.2.0"
bioimageIoSpec = "0.1.0"
-commonmark = "0.21.0"
+commonmark = "0.22.0"
commonsMath3 = "3.6.1"
commonsText = "1.10.0"
controlsFX = "11.1.2"
-deepJavaLibrary = "0.24.0"
+deepJavaLibrary = "0.27.0"
-groovy = "4.0.15"
+groovy = "4.0.21"
gson = "2.10.1"
-guava = "32.1.3-jre"
+guava = "33.1.0-jre"
ikonli = "12.3.1"
-imagej = "1.54f"
-
-# Compatibility with Java 17 with QuPath v0.5.x (may update to 21 in the future)
-jdk = "17"
-
-# Would ideally use JavaCPP 1.5.9 for gradle fixes: https://github.com/bytedeco/gradle-javacpp/issues/28
-# But need JavaCPP 1.5.8 for OpenCV DNN: https://github.com/qupath/qupath/issues/1406
-# When JavaCPP 1.5.10 is released we should be able to just use that for both.
-javacppgradle = "1.5.9"
-javacpp = "1.5.8"
-# Stuck on 4.6.0 because of problem with OpenCVN DNN running x86_64 build on Apple Silicon
-# https://github.com/qupath/qupath/issues/1406
-# This is fixed in 4.8.0, but that requires JavaCPP 1.5.10 (not released yet)
-opencv = "4.6.0-1.5.8"
-cuda = "11.8-8.6-1.5.8"
-# opencv = "4.7.0-1.5.9"
-# cuda = "12.1-8.9-1.5.9"
-
-
-# JavaFX 20.0.1 and later seem to break search links in Javadocs
-javafx = "20"
-jna = "5.13.0"
+imagej = "1.54i"
+
+# Compatibility with Java 21 with QuPath v0.6.x
+jdk = "21"
+
+javacppgradle = "1.5.10"
+javacpp = "1.5.10"
+
+opencv = "4.9.0-1.5.10"
+cuda = "12.3-8.9-1.5.10"
+
+# Warning! JavaFX 20.0.1 and later seem to break search links in Javadocs
+javafx = "22"
+jna = "5.14.0"
jfreeSvg = "5.0.5"
jfxtras = "17-r1"
jts = "1.19.0"
@@ -52,7 +44,7 @@ qupath-fxtras = "0.1.4"
richtextfx = "0.11.2"
-slf4j = "2.0.9"
+slf4j = "2.0.12"
snakeyaml = "2.2"
diff --git a/gradle/wrapper/gradle-wrapper.jar b/gradle/wrapper/gradle-wrapper.jar
index 7f93135c49b765f8051ef9d0a6055ff8e46073d8..d64cd4917707c1f8861d8cb53dd15194d4248596 100644
GIT binary patch
literal 43462
zcma&NWl&^owk(X(xVyW%ySuwf;qI=D6|RlDJ2cR^yEKh!@I-
zp9QeisK*rlxC>+~7Dk4IxIRsKBHqdR9b3+fyL=ynHmIDe&|>O*VlvO+%z5;9Z$|DJ
zb4dO}-R=MKr^6EKJiOrJdLnCJn>np?~vU-1sSFgPu;pthGwf}bG
z(1db%xwr#x)r+`4AGu$j7~u2MpVs3VpLp|mx&;>`0p0vH6kF+D2CY0fVdQOZ@h;A`
z{infNyvmFUiu*XG}RNMNwXrbec_*a3N=2zJ|Wh5z*
z5rAX$JJR{#zP>KY**>xHTuw?|-Rg|o24V)74HcfVT;WtQHXlE+_4iPE8QE#DUm%x0
zEKr75ur~W%w#-My3Tj`hH6EuEW+8K-^5P62$7Sc5OK+22qj&Pd1;)1#4tKihi=~8C
zHiQSst0cpri6%OeaR`PY>HH_;CPaRNty%WTm4{wDK8V6gCZlG@U3$~JQZ;HPvDJcT1V{
z?>H@13MJcCNe#5z+MecYNi@VT5|&UiN1D4ATT+%M+h4c$t;C#UAs3O_q=GxK0}8%8
z8J(_M9bayxN}69ex4dzM_P3oh@ZGREjVvn%%r7=xjkqxJP4kj}5tlf;QosR=%4L5y
zWhgejO=vao5oX%mOHbhJ8V+SG&K5dABn6!WiKl{|oPkq(9z8l&Mm%(=qGcFzI=eLu
zWc_oCLyf;hVlB@dnwY98?75B20=n$>u3b|NB28H0u-6Rpl((%KWEBOfElVWJx+5yg
z#SGqwza7f}$z;n~g%4HDU{;V{gXIhft*q2=4zSezGK~nBgu9-Q*rZ#2f=Q}i2|qOp
z!!y4p)4o=LVUNhlkp#JL{tfkhXNbB=Ox>M=n6soptJw-IDI|_$is2w}(XY>a=H52d
z3zE$tjPUhWWS+5h=KVH&uqQS=$v3nRs&p$%11b%5qtF}S2#Pc`IiyBIF4%A!;AVoI
zXU8-Rpv!DQNcF~(qQnyyMy=-AN~U>#&X1j5BLDP{?K!%h!;hfJI>$mdLSvktEr*89
zdJHvby^$xEX0^l9g$xW-d?J;L0#(`UT~zpL&*cEh$L|HPAu=P8`OQZV!-}l`noSp_
zQ-1$q$R-gDL)?6YaM!=8H=QGW$NT2SeZlb8PKJdc=F-cT@j7Xags+Pr*jPtlHFnf-
zh?q<6;)27IdPc^Wdy-mX%2s84C1xZq9Xms+==F4);O`VUASmu3(RlgE#0+#giLh-&
zcxm3_e}n4{%|X
zJp{G_j+%`j_q5}k{eW&TlP}J2wtZ2^<^E(O)4OQX8FDp6RJq!F{(6eHWSD3=f~(h}
zJXCf7=r<16X{pHkm%yzYI_=VDP&9bmI1*)YXZeB}F?
z(%QsB5fo*FUZxK$oX~X^69;x~j7ms8xlzpt-T15e9}$4T-pC
z6PFg@;B-j|Ywajpe4~bk#S6(fO^|mm1hKOPfA%8-_iGCfICE|=P_~e;Wz6my&)h_~
zkv&_xSAw7AZ%ThYF(4jADW4vg=oEdJGVOs>FqamoL3Np8>?!W#!R-0%2Bg4h?kz5I
zKV-rKN2n(vUL%D<4oj@|`eJ>0i#TmYBtYmfla;c!ATW%;xGQ0*TW@PTlGG><@dxUI
zg>+3SiGdZ%?5N=8uoLA|$4isK$aJ%i{hECP$bK{J#0W2gQ3YEa
zZQ50Stn6hqdfxJ*9#NuSLwKFCUGk@c=(igyVL;;2^wi4o30YXSIb2g_ud$
zgpCr@H0qWtk2hK8Q|&wx)}4+hTYlf;$a4#oUM=V@Cw#!$(nOFFpZ;0lc!qd=c$S}Z
zGGI-0jg~S~cgVT=4Vo)b)|4phjStD49*EqC)IPwyeKBLcN;Wu@Aeph;emROAwJ-0<
z_#>wVm$)ygH|qyxZaet&(Vf%pVdnvKWJn9`%DAxj3ot;v>S$I}jJ$FLBF*~iZ!ZXE
zkvui&p}fI0Y=IDX)mm0@tAd|fEHl~J&K}ZX(Mm3cm1UAuwJ42+AO5@HwYfDH7ipIc
zmI;1J;J@+aCNG1M`Btf>YT>~c&3j~Qi@Py5JT6;zjx$cvOQW@3oQ>|}GH?TW-E
z1R;q^QFjm5W~7f}c3Ww|awg1BAJ^slEV~Pk`Kd`PS$7;SqJZNj->it4DW2l15}xP6
zoCl$kyEF%yJni0(L!Z&14m!1urXh6Btj_5JYt1{#+H8w?5QI%%
zo-$KYWNMJVH?Hh@1n7OSu~QhSswL8x0=$<8QG_zepi_`y_79=nK=_ZP_`Em2UI*tyQoB+r{1QYZCpb?2OrgUw#oRH$?^Tj!Req>XiE#~B|~
z+%HB;=ic+R@px4Ld8mwpY;W^A%8%l8$@B@1m5n`TlKI6bz2mp*^^^1mK$COW$HOfp
zUGTz-cN9?BGEp}5A!mDFjaiWa2_J2Iq8qj0mXzk;
z66JBKRP{p%wN7XobR0YjhAuW9T1Gw3FDvR5dWJ8ElNYF94eF3ebu+QwKjtvVu4L
zI9ip#mQ@4uqVdkl-TUQMb^XBJVLW(-$s;Nq;@5gr4`UfLgF$adIhd?rHOa%D);whv
z=;krPp~@I+-Z|r#s3yCH+c1US?dnm+C*)r{m+86sTJusLdNu^sqLrfWed^ndHXH`m
zd3#cOe3>w-ga(Dus_^ppG9AC>Iq{y%%CK+Cro_sqLCs{VLuK=dev>OL1dis4(PQ5R
zcz)>DjEkfV+MO;~>VUlYF00SgfUo~@(&9$Iy2|G0T9BSP?&T22>K46D
zL*~j#yJ?)^*%J3!16f)@Y2Z^kS*BzwfAQ7K96rFRIh>#$*$_Io;z>ux@}G98!fWR@
zGTFxv4r~v)Gsd|pF91*-eaZ3Qw1MH$K^7JhWIdX%o$2kCbvGDXy)a?@8T&1dY4`;L
z4Kn+f%SSFWE_rpEpL9bnlmYq`D!6F%di<&Hh=+!VI~j)2mfil03T#jJ_s?}VV0_hp
z7T9bWxc>Jm2Z0WMU?`Z$xE74Gu~%s{mW!d4uvKCx@WD+gPUQ
zV0vQS(Ig++z=EHN)BR44*EDSWIyT~R4$FcF*VEY*8@l=218Q05D2$|fXKFhRgBIEE
zdDFB}1dKkoO^7}{5crKX!p?dZWNz$m>1icsXG2N+((x0OIST9Zo^DW_tytvlwXGpn
zs8?pJXjEG;T@qrZi%#h93?FP$!&P4JA(&H61tqQi=opRzNpm
zkrG}$^t9&XduK*Qa1?355wd8G2CI6QEh@Ua>AsD;7oRUNLPb76m4HG3K?)wF~IyS3`fXuNM>${?wmB
zpVz;?6_(Fiadfd{vUCBM*_kt$+F3J+IojI;9L(gc9n3{sEZyzR9o!_mOwFC#tQ{Q~
zP3-`#uK#tP3Q7~Q;4H|wjZHO8h7e4IuBxl&vz2w~D8)w=Wtg31zpZhz%+kzSzL*dV
zwp@{WU4i;hJ7c2f1O;7Mz6qRKeASoIv0_bV=i@NMG*l<#+;INk-^`5w@}Dj~;k=|}qM1vq_P
z|GpBGe_IKq|LNy9SJhKOQ$c=5L{Dv|Q_lZl=-ky*BFBJLW9&y_C|!vyM~rQx=!vun
z?rZJQB5t}Dctmui5i31C_;_}CEn}_W%>oSXtt>@kE1=JW*4*v4tPp;O6
zmAk{)m!)}34pTWg8{i>($%NQ(Tl;QC@J@FfBoc%Gr&m560^kgSfodAFrIjF}aIw)X
zoXZ`@IsMkc8_=w%-7`D6Y4e*CG8k%Ud=GXhsTR50jUnm+R*0A(O3UKFg0`K;qp1bl
z7``HN=?39ic_kR|^R^~w-*pa?Vj#7|e9F1iRx{GN2?wK!xR1GW!qa=~pjJb-#u1K8
zeR?Y2i-pt}yJq;SCiVHODIvQJX|ZJaT8nO+(?HXbLefulKKgM^B(UIO1r+S=7;kLJ
zcH}1J=Px2jsh3Tec&v8Jcbng8;V-`#*UHt?hB(pmOipKwf3Lz8rG$heEB30Sg*2rx
zV<|KN86$soN(I!BwO`1n^^uF2*x&vJ$2d$>+`(romzHP|)K_KkO6Hc>_dwMW-M(#S
zK(~SiXT1@fvc#U+?|?PniDRm01)f^#55;nhM|wi?oG>yBsa?~?^xTU|fX-R(sTA+5
zaq}-8Tx7zrOy#3*JLIIVsBmHYLdD}!0NP!+ITW+Thn0)8SS!$@)HXwB3tY!fMxc#1
zMp3H?q3eD?u&Njx4;KQ5G>32+GRp1Ee5qMO0lZjaRRu&{W<&~DoJNGkcYF<5(Ab+J
zgO>VhBl{okDPn78<%&e2mR{jwVCz5Og;*Z;;3%VvoGo_;HaGLWYF7q#jDX=Z#Ml`H
z858YVV$%J|e<1n`%6Vsvq7GmnAV0wW4$5qQ3uR@1i>tW{xrl|ExywIc?fNgYlA?C5
zh$ezAFb5{rQu6i7BSS5*J-|9DQ{6^BVQ{b*lq`xS@RyrsJN?-t=MTMPY;WYeKBCNg
z^2|pN!Q^WPJuuO4!|P@jzt&tY1Y8d%FNK5xK(!@`jO2aEA*4
zkO6b|UVBipci?){-Ke=+1;mGlND8)6+P;8sq}UXw2hn;fc7nM>g}GSMWu&v&fqhiViYT=fZ(|3Ox^$aWPp4a8h24tD<|8-!aK0lHgL$N7Efw}J
zVIB!7=T$U`ao1?upi5V4Et*-lTG0XvExbf!ya{cua==$WJyVG(CmA6Of*8E@DSE%L
z`V^$qz&RU$7G5mg;8;=#`@rRG`-uS18$0WPN@!v2d{H2sOqP|!(cQ@
zUHo!d>>yFArLPf1q`uBvY32miqShLT1B@gDL4XoVTK&@owOoD)OIHXrYK-a1d$B{v
zF^}8D3Y^g%^cnvScOSJR5QNH+BI%d|;J;wWM3~l>${fb8DNPg)wrf|GBP8p%LNGN#
z3EaIiItgwtGgT&iYCFy9-LG}bMI|4LdmmJt@V@%
zb6B)1kc=T)(|L@0;wr<>=?r04N;E&ef+7C^`wPWtyQe(*pD1pI_&XHy|0gIGHMekd
zF_*M4yi6J&Z4LQj65)S
zXwdM{SwUo%3SbPwFsHgqF@V|6afT|R6?&S;lw=8%
z3}@9B=#JI3@B*#4s!O))~z
zc>2_4Q_#&+5V`GFd?88^;c1i7;Vv_I*qt!_Yx*n=;rj!82rrR2rQ8u5(Ejlo{15P%
zs~!{%XJ>FmJ})H^I9bn^Re&38H{xA!0l3^89k(oU;bZWXM@kn$#aoS&Y4l^-WEn-fH39Jb9lA%s*WsKJQl?n9B7_~P
z-XM&WL7Z!PcoF6_D>V@$CvUIEy=+Z&0kt{szMk=f1|M+r*a43^$$B^MidrT0J;RI`
z(?f!O<8UZkm$_Ny$Hth1J#^4ni+im8M9mr&k|3cIgwvjAgjH
z8`N&h25xV#v*d$qBX5jkI|xOhQn!>IYZK7l5#^P4M&twe9&Ey@@GxYMxBZq2e7?`q
z$~Szs0!g{2fGcp9PZEt|rdQ6bhAgpcLHPz?f-vB?$dc*!9OL?Q8mn7->bFD2Si60*
z!O%y)fCdMSV|lkF9w%x~J*A&srMyYY3{=&$}H
zGQ4VG_?$2X(0|vT0{=;W$~icCI{b6W{B!Q8xdGhF|D{25G_5_+%s(46lhvNLkik~R
z>nr(&C#5wwOzJZQo9m|U<;&Wk!_#q|V>fsmj1g<6%hB{jGoNUPjgJslld>xmODzGjYc?7JSuA?A_QzjDw5AsRgi@Y|Z0{F{!1=!NES-#*f^s4l0Hu
zz468))2IY5dmD9pa*(yT5{EyP^G>@ZWumealS-*WeRcZ}B%gxq{MiJ|RyX-^C1V=0
z@iKdrGi1jTe8Ya^x7yyH$kBNvM4R~`fbPq$BzHum-3Zo8C6=KW@||>zsA8-Y9uV5V
z#oq-f5L5}V<&wF4@X@<3^C%ptp6+Ce)~hGl`kwj)bsAjmo_GU^r940Z-|`<)oGnh7
zFF0Tde3>ui?8Yj{sF-Z@)yQd~CGZ*w-6p2U<8}JO-sRsVI5dBji`01W8A&3$?}lxBaC&vn0E$c5tW*
zX>5(zzZ=qn&!J~KdsPl;P@bmA-Pr8T*)eh_+Dv5=Ma|XSle6t(k8qcgNyar{*ReQ8
zTXwi=8vr>!3Ywr+BhggHDw8ke==NTQVMCK`$69fhzEFB*4+H9LIvdt-#IbhZvpS}}
zO3lz;P?zr0*0$%-Rq_y^k(?I{Mk}h@w}cZpMUp|ucs55bcloL2)($u%mXQw({Wzc~
z;6nu5MkjP)0C(@%6Q_I_vsWrfhl7Zpoxw#WoE~r&GOSCz;_ro6i(^hM>I$8y>`!wW
z*U^@?B!MMmb89I}2(hcE4zN2G^kwyWCZp5JG>$Ez7zP~D=J^LMjSM)27_0B_X^C(M
z`fFT+%DcKlu?^)FCK>QzSnV%IsXVcUFhFdBP!6~se&xxrIxsvySAWu++IrH;FbcY$
z2DWTvSBRfLwdhr0nMx+URA$j3i7_*6BWv#DXfym?ZRDcX9C?cY9sD3q)uBDR3uWg=
z(lUIzB)G$Hr!){>E{s4Dew+tb9kvToZp-1&c?y2wn@Z~(VBhqz`cB;{E4(P3N2*nJ
z_>~g@;UF2iG{Kt(<1PyePTKahF8<)pozZ*xH~U-kfoAayCwJViIrnqwqO}7{0pHw$
zs2Kx?s#vQr7XZ264>5RNKSL8|Ty^=PsIx^}QqOOcfpGUU4tRkUc|kc7-!Ae6!+B{o~7nFpm3|G5^=0#Bnm6`V}oSQlrX(u%OWnC
zoLPy&Q;1Jui&7ST0~#+}I^&?vcE*t47~Xq#YwvA^6^}
z`WkC)$AkNub|t@S!$8CBlwbV~?yp&@9h{D|3z-vJXgzRC5^nYm+PyPcgRzAnEi6Q^gslXYRv4nycsy-SJu?lMps-?
zV`U*#WnFsdPLL)Q$AmD|0`UaC4ND07+&UmOu!eHruzV|OUox<+Jl|Mr@6~C`T@P%s
zW7sgXLF2SSe9Fl^O(I*{9wsFSYb2l%-;&Pi^dpv!{)C3d0AlNY6!4fgmSgj_wQ*7Am7&$z;Jg&wgR-Ih;lUvWS|KTSg!&s_E9_bXBkZvGiC6bFKDWZxsD$*NZ#_8bl
zG1P-#@?OQzED7@jlMJTH@V!6k;W>auvft)}g
zhoV{7$q=*;=l{O>Q4a@
ziMjf_u*o^PsO)#BjC%0^h>Xp@;5$p{JSYDt)zbb}s{Kbt!T*I@Pk@X0zds6wsefuU
zW$XY%yyRGC94=6mf?x+bbA5CDQ2AgW1T-jVAJbm7K(gp+;v6E0WI#kuACgV$r}6L?
zd|Tj?^%^*N&b>Dd{Wr$FS2qI#Ucs1yd4N+RBUQiSZGujH`#I)mG&VKoDh=KKFl4=G
z&MagXl6*<)$6P}*Tiebpz5L=oMaPrN+caUXRJ`D?=K9!e0f{@D&cZLKN?iNP@X0aF
zE(^pl+;*T5qt?1jRC=5PMgV!XNITRLS_=9{CJExaQj;lt!&pdzpK?8p>%Mb+D
z?yO*uSung=-`QQ@yX@Hyd4@CI^r{2oiu`%^bNkz+Nkk!IunjwNC|WcqvX~k=><-I3
zDQdbdb|!v+Iz01$w@aMl!R)koD77Xp;eZwzSl-AT
zr@Vu{=xvgfq9akRrrM)}=!=xcs+U1JO}{t(avgz`6RqiiX<|hGG1pmop8k6Q+G_mv
zJv|RfDheUp2L3=^C=4aCBMBn0aRCU(DQwX-W(RkRwmLeuJYF<0urcaf(=7)JPg<3P
zQs!~G)9CT18o!J4{zX{_e}4eS)U-E)0FAt}wEI(c0%HkxgggW;(1E=>J17_hsH^sP
z%lT0LGgbUXHx-K*CI-MCrP66UP0PvGqM$MkeLyqHdbgP|_Cm!7te~b8p+e6sQ_3k|
zVcwTh6d83ltdnR>D^)BYQpDKlLk3g0Hdcgz2}%qUs9~~Rie)A-BV1mS&naYai#xcZ
z(d{8=-LVpTp}2*y)|gR~;qc7fp26}lPcLZ#=JpYcn3AT9(UIdOyg+d(P5T7D&*P}#
zQCYplZO5|7+r19%9e`v^vfSS1sbX1c%=w1;oyruXB%Kl$ACgKQ6=qNWLsc=28xJjg
zwvsI5-%SGU|3p>&zXVl^vVtQT3o-#$UT9LI@Npz~6=4!>mc431VRNN8od&Ul^+G_kHC`G=6WVWM
z%9eWNyy(FTO|A+@x}Ou3CH)oi;t#7rAxdIXfNFwOj_@Y&TGz6P_sqiB`Q6Lxy|Q{`|fgmRG(k+!#b*M+Z9zFce)f-7;?Km5O=LHV9f9_87;
zF7%R2B+$?@sH&&-$@tzaPYkw0;=i|;vWdI|Wl3q_Zu>l;XdIw2FjV=;Mq5t1Q0|f<
zs08j54Bp`3RzqE=2enlkZxmX6OF+@|2<)A^RNQpBd6o@OXl+i)zO%D4iGiQNuXd+zIR{_lb96{lc~bxsBveIw6umhShTX+3@ZJ=YHh@
zWY3(d0azg;7oHn>H<>?4@*RQbi>SmM=JrHvIG(~BrvI)#W(EAeO6fS+}mxxcc+X~W6&YVl86W9WFSS}Vz-f9vS?XUDBk)3TcF
z8V?$4Q)`uKFq>xT=)Y9mMFVTUk*NIA!0$?RP6Ig0TBmUFrq*Q-Agq~DzxjStQyJ({
zBeZ;o5qUUKg=4Hypm|}>>L=XKsZ!F$yNTDO)jt4H0gdQ5$f|d&bnVCMMXhNh)~mN
z@_UV6D7MVlsWz+zM+inZZp&P4fj=tm6fX)SG5H>OsQf_I8c~uGCig$GzuwViK54bcgL;VN|FnyQl>Ed7(@>=8$a_UKIz|V6CeVSd2(P
z0Uu>A8A+muM%HLFJQ9UZ5c)BSAv_zH#1f02x?h9C}@pN@6{>UiAp>({Fn(T9Q8B
z^`zB;kJ5b`>%dLm+Ol}ty!3;8f1XDSVX0AUe5P#@I+FQ-`$(a;zNgz)4x5hz$Hfbg
z!Q(z26wHLXko(1`;(BAOg_wShpX0ixfWq3ponndY+u%1gyX)_h=v1zR#V}#q{au6;
z!3K=7fQwnRfg6FXtNQmP>`<;!N137paFS%y?;lb1@BEdbvQHYC{976l`cLqn;b8lp
zIDY>~m{gDj(wfnK!lpW6pli)HyLEiUrNc%eXTil|F2s(AY+LW5hkKb>TQ3|Q4S9rr
zpDs4uK_co6XPsn_z$LeS{K4jFF`2>U`tbgKdyDne`xmR<@6AA+_hPNKCOR-Zqv;xk
zu5!HsBUb^!4uJ7v0RuH-7?l?}b=w5lzzXJ~gZcxRKOovSk@|#V+MuX%Y+=;14i*%{)_gSW9(#4%)AV#3__kac1|qUy!uyP{>?U#5wYNq}y$S9pCc
zFc~4mgSC*G~j0u#qqp9
z${>3HV~@->GqEhr_Xwoxq?Hjn#=s2;i~g^&Hn|aDKpA>Oc%HlW(KA1?BXqpxB;Ydx)w;2z^MpjJ(Qi(X!$5RC
z*P{~%JGDQqojV>2JbEeCE*OEu!$XJ>bWA9Oa_Hd;y)F%MhBRi*LPcdqR8X`NQ&1L#
z5#9L*@qxrx8n}LfeB^J{%-?SU{FCwiWyHp682F+|pa+CQa3ZLzBqN1{)h4d6+vBbV
zC#NEbQLC;}me3eeYnOG*nXOJZEU$xLZ1<1Y=7r0(-U0P6-AqwMAM`a(Ed#7vJkn6plb4eI4?2y3yOTGmmDQ!z9`wzbf
z_OY#0@5=bnep;MV0X_;;SJJWEf^E6Bd^tVJ9znWx&Ks8t*B>AM@?;D4oWUGc
z!H*`6d7Cxo6VuyS4Eye&L1ZRhrRmN6Lr`{NL(wDbif|y&z)JN>Fl5#Wi&mMIr5i;x
zBx}3YfF>>8EC(fYnmpu~)CYHuHCyr5*`ECap%t@y=jD>!_%3iiE|LN$mK9>-
zHdtpy8fGZtkZF?%TW~29JIAfi2jZT8>OA7=h;8T{{k?c2`nCEx9$r
zS+*&vt~2o^^J+}RDG@+9&M^K*z4p{5#IEVbz`1%`m5c2};aGt=V?~vIM}ZdPECDI)47|CWBCfDWUbxBCnmYivQ*0Nu_xb*C>~C9(VjHM
zxe<*D<#dQ8TlpMX2c@M<9$w!RP$hpG4cs%AI){jp*Sj|*`m)5(Bw*A0$*i-(CA5#%>a)$+jI2C9r6|(>J8InryENI
z$NohnxDUB;wAYDwrb*!N3noBTKPpPN}~09SEL18tkG
zxgz(RYU_;DPT{l?Q$+eaZaxnsWCA^ds^0PVRkIM%bOd|G2IEBBiz{&^JtNsODs;5z
zICt_Zj8wo^KT$7Bg4H+y!Df#3mbl%%?|EXe!&(Vmac1DJ*y~3+kRKAD=Ovde4^^%~
zw<9av18HLyrf*_>Slp;^i`Uy~`mvBjZ|?Ad63yQa#YK`4+c6;pW4?XIY9G1(Xh9WO8{F-Aju+nS9Vmv=$Ac0ienZ+p9*O%NG
zMZKy5?%Z6TAJTE?o5vEr0r>f>hb#2w2U3DL64*au_@P!J!TL`oH2r*{>ffu6|A7tv
zL4juf$DZ1MW5ZPsG!5)`k8d8c$J$o;%EIL0va9&GzWvkS%ZsGb#S(?{!UFOZ9<$a|
zY|a+5kmD5N&{vRqkgY>aHsBT&`rg|&kezoD)gP0fsNYHsO#TRc_$n6Lf1Z{?+DLziXlHrq4sf(!>O{?Tj;Eh@%)+nRE_2VxbN&&%%caU#JDU%vL3}Cb
zsb4AazPI{>8H&d=jUaZDS$-0^AxE@utGs;-Ez_F(qC9T=UZX=>ok2k2
ziTn{K?y~a5reD2A)P${NoI^>JXn>`IeArow(41c-Wm~)wiryEP(OS{YXWi7;%dG9v
zI?mwu1MxD{yp_rrk!j^cKM)dc4@p4Ezyo%lRN|XyD}}>v=Xoib0gOcdXrQ^*61HNj
z=NP|pd>@yfvr-=m{8$3A8TQGMTE7g=z!%yt`8`Bk-0MMwW~h^++;qyUP!J~ykh1GO
z(FZ59xuFR$(WE;F@UUyE@Sp>`aVNjyj=Ty>_Vo}xf`e7`F;j-IgL5`1~-#70$9_=uBMq!2&1l
zomRgpD58@)YYfvLtPW}{C5B35R;ZVvB<<#)x%srmc_S=A7F@DW8>QOEGwD6suhwCg
z>Pa+YyULhmw%BA*4yjDp|2{!T98~<6Yfd(wo1mQ!KWwq0eg+6)o1>W~f~kL<-S+P@$wx*zeI|1t7z#Sxr5
zt6w+;YblPQNplq4Z#T$GLX#j6yldXAqj>4gAnnWtBICUnA&-dtnlh=t0Ho_vEKwV`
z)DlJi#!@nkYV#$!)@>udAU*hF?V`2$Hf=V&6PP_|r#Iv*J$9)pF@X3`k;5})9^o4y
z&)~?EjX5yX12O(BsFy-l6}nYeuKkiq`u9145&3Ssg^y{5G3Pse
z9w(YVa0)N-fLaBq1`P!_#>SS(8fh_5!f{UrgZ~uEdeMJIz7DzI5!NHHqQtm~#CPij
z?=N|J>nPR6_sL7!f4hD_|KH`vf8(Wpnj-(gPWH+ZvID}%?~68SwhPTC3u1_cB`otq
z)U?6qo!ZLi5b>*KnYHWW=3F!p%h1;h{L&(Q&{qY6)_qxNfbP6E3yYpW!EO+IW3?@J
z);4>g4gnl^8klu7uA>eGF6rIGSynacogr)KUwE_R4E5Xzi*Qir@b-jy55-JPC8c~(
zo!W8y9OGZ&`xmc8;=4-U9=h{vCqfCNzYirONmGbRQlR`WWlgnY+1wCXbMz&NT~9*|
z6@FrzP!LX&{no2!Ln_3|I==_4`@}V?4a;YZKTdw;vT<+K+z=uWbW(&bXEaWJ^W8Td
z-3&1bY^Z*oM<=M}LVt>_j+p=2Iu7pZmbXrhQ_k)ysE9yXKygFNw$5hwDn(M>H+e1&9BM5!|81vd%r%vEm
zqxY3?F@fb6O#5UunwgAHR9jp_W2zZ}NGp2%mTW@(hz7$^+a`A?mb8|_G*GNMJ)
zjqegXQio=i@AINre&%ofexAr95aop5C+0MZ0m-l=MeO8m3epm7U%vZB8+I+C*iNFM
z#T3l`gknX;D$-`2XT^Cg*vrv=RH+P;_dfF++cP?B_msQI4j+lt&rX2)3GaJx%W*Nn
zkML%D{z5tpHH=dksQ*gzc|}gzW;lwAbxoR07VNgS*-c3d&8J|;@3t^
zVUz*J*&r7DFRuFVDCJDK8V9NN5hvpgGjwx+5n)qa;YCKe8TKtdnh{I7NU9BCN!0dq
zczrBk8pE{{@vJa9ywR@mq*J=v+PG;?fwqlJVhijG!3VmIKs>9T6r7MJpC)m!Tc#>g
zMtVsU>wbwFJEfwZ{vB|ZlttNe83)$iz`~#8UJ^r)lJ@HA&G#}W&ZH*;k{=TavpjWE
z7hdyLZPf*X%Gm}i`Y{OGeeu^~nB8=`{r#TUrM-`;1cBvEd#d!kPqIgYySYhN-*1;L
z^byj%Yi}Gx)Wnkosi337BKs}+5H5dth1JA{Ir-JKN$7zC)*}hqeoD(WfaUDPT>0`-
z(6sa0AoIqASwF`>hP}^|)a_j2s^PQn*qVC{Q}htR
z5-)duBFXT_V56-+UohKXlq~^6uf!6sA#ttk1o~*QEy_Y-S$gAvq47J9Vtk$5oA$Ct
zYhYJ@8{hsC^98${!#Ho?4y5MCa7iGnfz}b9jE~h%EAAv~Qxu)_rAV;^cygV~5r_~?l=B`zObj7S=H=~$W
zPtI_m%g$`kL_fVUk9J@>EiBH
zOO&jtn~&`hIFMS5S`g8w94R4H40mdNUH4W@@XQk1sr17b{@y|JB*G9z1|CrQjd+GX
z6+KyURG3;!*BQrentw{B2R&@2&`2}n(z-2&X7#r!{yg@Soy}cRD~j
zj9@UBW+N|4HW4AWapy4wfUI-
zZ`gSL6DUlgj*f1hSOGXG0IVH8HxK?o2|3HZ;KW{K+yPAlxtb)NV_2AwJm|E)FRs&&
z=c^e7bvUsztY|+f^k7NXs$o1EUq>cR7C0$UKi6IooHWlK_#?IWDkvywnzg&ThWo^?
z2O_N{5X39#?eV9l)xI(>@!vSB{DLt*oY!K1R8}_?%+0^C{d9a%N4
zoxHVT1&Lm|uDX%$QrBun5e-F`HJ^T$
zmzv)p@4ZHd_w9!%Hf9UYNvGCw2TTTbrj9pl+T9%-_-}L(tES>Or-}Z4F*{##n3~L~TuxjirGuIY#H7{%$E${?p{Q01
zi6T`n;rbK1yIB9jmQNycD~yZq&mbIsFWHo|ZAChSFPQa<(%d8mGw*V3fh|yFoxOOiWJd(qvVb!Z$b88cg->N=qO*4k~6;R==|9ihg&riu#P~s4Oap9O7f%crSr^rljeIfXDEg>wi)&v*a%7zpz<9w
z*r!3q9J|390x`Zk;g$&OeN&ctp)VKRpDSV@kU2Q>jtok($Y-*x8_$2piTxun81@vt
z!Vj?COa0fg2RPXMSIo26T=~0d`{oGP*eV+$!0I<(4azk&Vj3SiG=Q!6mX0p$z7I};
z9BJUFgT-K9MQQ-0@Z=^7R<{bn2Fm48endsSs`V7_@%8?Bxkqv>BDoVcj?K#dV#uUP
zL1ND~?D-|VGKe3Rw_7-Idpht>H6XRLh*U7epS6byiGvJpr%d}XwfusjH9g;Z98H`x
zyde%%5mhGOiL4wljCaWCk-&uE4_OOccb9c!ZaWt4B(wYl!?vyzl%7n~QepN&eFUrw
zFIOl9c({``6~QD+43*_tzP{f2x41h(?b43^y6=iwyB)2os5hBE!@YUS5?N_tXd=h(
z)WE286Fbd>R4M^P{!G)f;h<3Q>Fipuy+d2q-)!RyTgt;wr$(?9ox3;q+{E*ZQHhOn;lM`cjnu9
zXa48ks-v(~b*;MAI<>YZH(^NV8vjb34beE<_cwKlJoR;k6lJNSP6v}uiyRD?|0w+X@o1ONrH8a$fCxXpf?
z?$DL0)7|X}Oc%h^zrMKWc-NS9I0Utu@>*j}b@tJ=ixQSJ={4@854wzW@E>VSL+Y{i
z#0b=WpbCZS>kUCO_iQz)LoE>P5LIG-hv9E+oG}DtlIDF>$tJ1aw9^LuhLEHt?BCj&
z(O4I8v1s#HUi5A>nIS-JK{v!7dJx)^Yg%XjNmlkWAq2*cv#tHgz`Y(bETc6CuO1VkN^L-L3j_x<4NqYb5rzrLC-7uOv
z!5e`GZt%B782C5-fGnn*GhDF$%(qP<74Z}3xx+{$4cYKy2ikxI7B2N+2r07DN;|-T->nU&!=Cm#rZt%O_5c&1Z%nlWq3TKAW0w
zQqemZw_ue--2uKQsx+niCUou?HjD`xhEjjQd3%rrBi82crq*~#uA4+>vR<_S{~5ce
z-2EIl?~s
z1=GVL{NxP1N3%=AOaC}j_Fv=ur&THz
zyO!d9kHq|c73kpq`$+t+8Bw7MgeR5~`d7ChYyGCBWSteTB>8WAU(NPYt2Dk`@#+}=
zI4SvLlyk#pBgVigEe`?NG*vl7V6m+<}%FwPV=~PvvA)=#ths==DRTDEYh4V5}Cf$z@#;<
zyWfLY_5sP$gc3LLl2x+Ii)#b2nhNXJ{R~vk`s5U7Nyu^3yFg&D%Txwj6QezMX`V(x
z=C`{76*mNb!qHHs)#GgGZ_7|vkt9izl_&PBrsu@}L`X{95-2jf99K)0=*N)VxBX2q
z((vkpP2RneSIiIUEnGb?VqbMb=Zia+rF~+iqslydE34cSLJ&BJW^3knX@M;t*b=EA
zNvGzv41Ld_T+WT#XjDB840vovUU^FtN_)G}7v)1lPetgpEK9YS^OWFkPoE{ovj^=@
zO9N$S=G$1ecndT_=5ehth2Lmd1II-PuT~C9`XVePw$y8J#dpZ?Tss<6wtVglm(Ok7
z3?^oi@pPio6l&!z8JY(pJvG=*pI?GIOu}e^EB6QYk$#FJQ%^AIK$I4epJ+9t?KjqA+bkj&PQ*|vLttme+`9G=L%
ziadyMw_7-M)hS(3E$QGNCu|o23|%O+VN7;Qggp?PB3K-iSeBa2b}V4_wY`G1Jsfz4
z9|SdB^;|I8E8gWqHKx!vj_@SMY^hLEIbSMCuE?WKq=c2mJK
z8LoG-pnY!uhqFv&L?yEuxo{dpMTsmCn)95xanqBrNPTgXP((H$9N${Ow~Is-FBg%h
z53;|Y5$MUN)9W2HBe2TD`ct^LHI<(xWrw}$qSoei?}s)&w$;&!14w6B6>Yr6Y8b)S
z0r71`WmAvJJ`1h&poLftLUS6Ir
zC$bG9!Im_4Zjse)#K=oJM9mHW1{%l8sz$1o?ltdKlLTxWWPB>Vk22czVt|1%^wnN@*!l)}?EgtvhC>vlHm^t+ogpgHI1_$1ox9e;>0!+b(tBrmXRB`PY1vp-R**8N7
zGP|QqI$m(Rdu#=(?!(N}G9QhQ%o!aXE=aN{&wtGP8|_qh+7a_j_sU5|J^)vxq;#
zjvzLn%_QPHZZIWu1&mRAj;Sa_97p_lLq_{~j!M9N^1yp3U_SxRqK&JnR%6VI#^E12
z>CdOVI^_9aPK2eZ4h&^{pQs}xsijXgFYRIxJ~N7&BB9jUR1fm!(xl)mvy|3e6-B3j
zJn#ajL;bFTYJ2+Q)tDjx=3IklO@Q+FFM}6UJr6km7hj7th9n_&JR7fnqC!hTZoM~T
zBeaVFp%)0cbPhejX<8pf5HyRUj2>aXnXBqDJe73~J%P(2C?-RT{c3NjE`)om!
zl$uewSgWkE66$Kb34+QZZvRn`fob~Cl9=cRk@Es}KQm=?E~CE%spXaMO6YmrMl%9Q
zlA3Q$3|L1QJ4?->UjT&CBd!~ru{Ih^in&JXO=|<6J!&qp
zRe*OZ*cj5bHYlz!!~iEKcuE|;U4vN1rk$xq6>bUWD*u(V@8sG^7>kVuo(QL@Ki;yL
zWC!FT(q{E8#on>%1iAS0HMZDJg{Z{^!De(vSIq&;1$+b)oRMwA3nc3mdTSG#3uYO_
z>+x;7p4I;uHz?ZB>dA-BKl+t-3IB!jBRgdvAbW!aJ(Q{aT>+iz?91`C-xbe)IBoND
z9_Xth{6?(y3rddwY$GD65IT#f3<(0o#`di{sh2gm{dw*#-Vnc3r=4==&PU^hCv$qd
zjw;>i&?L*Wq#TxG$mFIUf>eK+170KG;~+o&1;Tom9}}mKo23KwdEM6UonXgc
z!6N(@k8q@HPw{O8O!lAyi{rZv|DpgfU{py+j(X_cwpKqcalcqKIr0kM^%Br3SdeD>
zHSKV94Yxw;pjzDHo!Q?8^0bb%L|wC;4U^9I#pd5O&eexX+Im{
z?jKnCcsE|H?{uGMqVie_C~w7GX)kYGWAg%-?8|N_1#W-|4F)3YTDC+QSq1s!DnOML3@d`mG%o2YbYd#jww|jD$gotpa)kntakp#K;+yo-_ZF9qrNZw<%#C
zuPE@#3RocLgPyiBZ+R_-FJ_$xP!RzWm|aN)S+{$LY9vvN+IW~Kf3TsEIvP+B9Mtm!
zpfNNxObWQpLoaO&cJh5>%slZnHl_Q~(-Tfh!DMz(dTWld@LG1VRF`9`DYKhyNv
z2pU|UZ$#_yUx_B_|MxUq^glT}O5Xt(Vm4Mr02><%C)@v;vPb@pT$*yzJ4aPc_FZ3z
z3}PLoMBIM>q_9U2rl^sGhk1VUJ89=*?7|v`{!Z{6bqFMq(mYiA?%KbsI~JwuqVA9$H5vDE+VocjX+G^%bieqx->s;XWlKcuv(s%y%D5Xbc9+
zc(_2nYS1&^yL*ey664&4`IoOeDIig}y-E~_GS?m;D!xv5-xwz+G`5l6V+}CpeJDi^
z%4ed$qowm88=iYG+(`ld5Uh&>Dgs4uPHSJ^TngXP_V6fPyl~>2bhi20QB%lSd#yYn
zO05?KT1z@?^-bqO8Cg`;ft>ilejsw@2%RR7;`$Vs;FmO(Yr3Fp`pHGr@P2hC%QcA|X&N2Dn
zYf`MqXdHi%cGR@%y7Rg7?d3?an){s$zA{!H;Ie5exE#c~@NhQUFG8V=SQh%UxUeiV
zd7#UcYqD=lk-}sEwlpu&H^T_V0{#G?lZMxL7ih_&{(g)MWBnCZxtXg
znr#}>U^6!jA%e}@Gj49LWG@*&t0V>Cxc3?oO7LSG%~)Y5}f7vqUUnQ;STjdDU}P9IF9d9<$;=QaXc
zL1^X7>fa^jHBu_}9}J~#-oz3Oq^JmGR#?GO7b9a(=R@fw@}Q{{@`Wy1vIQ#Bw?>@X
z-_RGG@wt|%u`XUc%W{J
z>iSeiz8C3H7@St3mOr_mU+&bL#Uif;+Xw-aZdNYUpdf>Rvu0i0t6k*}vwU`XNO2he
z%miH|1tQ8~ZK!zmL&wa3E;l?!!XzgV#%PMVU!0xrDsNNZUWKlbiOjzH-1Uoxm8E#r`#2Sz;-o&qcqB
zC-O_R{QGuynW14@)7&@yw1U}uP(1cov)twxeLus0s|7ayrtT8c#`&2~Fiu2=R;1_4bCaD=*E@cYI>7YSnt)nQc
zohw5CsK%m?8Ack)qNx`W0_v$5S}nO|(V|RZKBD+btO?JXe|~^Qqur%@eO~<8-L^9d
z=GA3-V14ng9L29~XJ>a5k~xT2152zLhM*@zlp2P5Eu}bywkcqR;ISbasT#;HZSf
z2m69qTV(V@EkY(1Dk3`}j)JMo%ZVJ*5eB
zYOjIisi+igK0#yW*gBGj?@I{~mUOvRFQR^pJbEbzFxTubnrw(Muk%}jI+vXmJ;{Q6
zrSobKD>T%}jV4Ub?L1+MGOD~0Ir%-`iTnWZN^~YPrcP5y3VMAzQ+&en^VzKEb$K!Q
z<7Dbg&DNXuow*eD5yMr+#08nF!;%4vGrJI++5HdCFcGLfMW!KS*Oi@=7hFwDG!h2<
zPunUEAF+HncQkbfFj&pbzp|MU*~60Z(|Ik%Tn{BXMN!hZOosNIseT?R;A`W?=d?5X
zK(FB=9mZusYahp|K-wyb={rOpdn=@;4YI2W0EcbMKyo~-#^?h`BA9~o285%oY
zfifCh5Lk$SY@|2A@a!T2V+{^!psQkx4?x0HSV`(w9{l75QxMk!)U52Lbhn{8ol?S)
zCKo*7R(z!uk<6*qO=wh!Pul{(qq6g6xW;X68GI_CXp`XwO
zxuSgPRAtM8K7}5E#-GM!*ydOOG_{A{)hkCII<|2=ma*71ci_-}VPARm3crFQjLYV!
z9zbz82$|l01mv`$WahE2$=fAGWkd^X2kY(J7iz}WGS
z@%MyBEO=A?HB9=^?nX`@nh;7;laAjs+fbo!|K^mE!tOB>$2a_O0y-*uaIn8k^6Y
zSbuv;5~##*4Y~+y7Z5O*3w4qgI5V^17u*ZeupVGH^nM&$qmAk|anf*>r
zWc5CV;-JY-Z@Uq1Irpb^O`L_7AGiqd*YpGUShb==os$uN3yYvb`wm6d=?T*it&pDk
zo`vhw)RZX|91^^Wa_ti2zBFyWy4cJu#g)_S6~jT}CC{DJ_kKpT`$oAL%b^!2M;JgT
zM3ZNbUB?}kP(*YYvXDIH8^7LUxz5oE%kMhF!rnPqv!GiY0o}NR$OD=ITDo9r%4E>E0Y^R(rS^~XjWyVI6
zMOR5rPXhTp*G*M&X#NTL`Hu*R+u*QNoiOKg4CtNPrjgH>c?Hi4MUG#I917fx**+pJfOo!zFM&*da&G_x)L(`k&TPI*t3e^{crd
zX<4I$5nBQ8Ax_lmNRa~E*zS-R0sxkz`|>7q_?*e%7bxqNm3_eRG#1ae3gtV9!fQpY
z+!^a38o4ZGy9!J5sylDxZTx$JmG!wg7;>&5H1)>f4dXj;B+@6tMlL=)cLl={jLMxY
zbbf1ax3S4>bwB9-$;SN2?+GULu;UA-35;VY*^9Blx)Jwyb$=U!D>HhB&=jSsd^6yw
zL)?a|>GxU!W}ocTC(?-%z3!IUhw^uzc`Vz_g>-tv)(XA#JK^)ZnC|l1`@CdX1@|!|
z_9gQ)7uOf?cR@KDp97*>6X|;t@Y`k_N@)aH7gY27)COv^P3ya9I{4z~vUjLR9~z1Z
z5=G{mVtKH*&$*t0@}-i_v|3B$AHHYale7>E+jP`ClqG%L{u;*ff_h@)al?RuL7tOO
z->;I}>%WI{;vbLP3VIQ^iA$4wl6@0sDj|~112Y4OFjMs`13!$JGkp%b&E8QzJw_L5
zOnw9joc0^;O%OpF$Qp)W1HI!$4BaXX84`%@#^dk^hFp^pQ@rx4g(8Xjy#!X%+X5Jd@fs3amGT`}mhq#L97R>OwT5-m|h#yT_-v@(k$q7P*9X~T*3)LTdzP!*B}
z+SldbVWrrwQo9wX*%FyK+sRXTa@O?WM^FGWOE?S`R(0P{<6p#f?0NJvnBia?k^fX2
zNQs7K-?EijgHJY}&zsr;qJ<*PCZUd*x|dD=IQPUK_nn)@X4KWtqoJNHkT?ZWL_hF?
zS8lp2(q>;RXR|F;1O}EE#}gCrY~#n^O`_I&?&z5~7N;zL0)3Tup`%)oHMK-^r$NT%
zbFg|o?b9w(q@)6w5V%si<$!U<#}s#x@0aX-hP>zwS#9*75VXA4K*%gUc>+yzupTDBOKH8WR4V0pM(HrfbQ&eJ79>HdCvE=F
z|J>s;;iDLB^3(9}?biKbxf1$lI!*Z%*0&8UUq}wMyPs_hclyQQi4;NUY+x2qy|0J;
zhn8;5)4ED1oHwg+VZF|80<4MrL97tGGXc5Sw$wAI#|2*cvQ=jB5+{AjMiDHmhUC*a
zlmiZ`LAuAn_}hftXh;`Kq0zblDk8?O-`tnilIh|;3lZp@F_osJUV9`*R29M?7H{Fy
z`nfVEIDIWXmU&YW;NjU8)EJpXhxe5t+scf|VXM!^bBlwNh)~7|3?fWwo_~ZFk(22%
zTMesYw+LNx3J-_|DM~`v93yXe=jPD{q;li;5PD?Dyk+b?
zo21|XpT@)$BM$%F=P9J19Vi&1#{jM3!^Y&fr&_`toi`XB1!n>sbL%U9I5<7!@?t)~
z;&H%z>bAaQ4f$wIzkjH70;<8tpUoxzKrPhn#IQfS%9l5=Iu))^XC<58D!-O
z{B+o5R^Z21H0T9JQ5gNJnqh#qH^na|z92=hONIM~@_iuOi|F>jBh-?aA20}Qx~EpDGElELNn~|7WRXRFnw+Wdo`|#
zBpU=Cz3z%cUJ0mx_1($X<40XEIYz(`noWeO+x#yb_pwj6)R(__%@_Cf>txOQ74wSJ
z0#F3(zWWaR-jMEY$7C*3HJrohc79>MCUu26mfYN)f4M~4gD`}EX4e}A!U}QV8!S47
z6y-U-%+h`1n`*pQuKE%Av0@)+wBZr9mH}@vH@i{v(m-6QK7Ncf17x_D=)32`FOjjo
zg|^VPf5c6-!FxN{25dvVh#fog=NNpXz
zfB$o+0jbRkHH{!TKhE709f+jI^$3#v1Nmf80w`@7-5$1Iv_`)W^px8P-({xwb;D0y
z7LKDAHgX<84?l!I*Dvi2#D@oAE^J|g$3!)x1Ua;_;<@#l1fD}lqU2_tS^6Ht$1Wl}
zBESo7o^)9-Tjuz$8YQSGhfs{BQV6zW7dA?0b(Dbt=UnQs&4zHfe_sj{RJ4uS-vQpC
zX;Bbsuju4%!o8?&m4UZU@~Z