Skip to content

use HM for Resampling and CRS related Operations #306

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Merged
AM1729 merged 3 commits into from
May 18, 2026
Merged

use HM for Resampling and CRS related Operations #306

Changes from all commits
Commits
Show all changes
3 commits
Select commit Hold shift + click to select a range
0c3e90e
use HMCRSRegistry instead of .Decode
AM1729 May 18, 2026
b487c1e
fixes
AM1729 May 18, 2026
776b39b
fixes
AM1729 May 18, 2026
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
38 changes: 19 additions & 19 deletions adapters/klab.ogc/src/main/java/org/integratedmodelling/klab/stac/STACEncoder.java
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -32,6 +32,9 @@
import org.hortonmachine.gears.utils.RegionMap;
import org.hortonmachine.gears.utils.crs.CrsUtilities;
import org.hortonmachine.gears.utils.geometry.GeometryUtilities;
import org.hortonmachine.gears.utils.crs.HMCrsRegistry;
import org.hortonmachine.gears.utils.crs.HMCrsTransformer;

import org.integratedmodelling.klab.Authentication;
import org.integratedmodelling.klab.Observables;
import org.integratedmodelling.klab.api.data.IGeometry;
Expand Down Expand Up @@ -375,20 +378,18 @@ public boolean test(HMStacAsset asset) { // Assuming for now that "eo:bands" wou
}
HMRaster outRaster = collection.readRasterBandOnRegion(regionTransformed, predicate, items, true,
MergeMode.SUBSTITUTE, lpm);
CoordinateReferenceSystem targetCRS = HMCrsRegistry.INSTANCE.getCrs("4326");
if (!HMCrsRegistry.crsEquals(outRaster.getCrs(),targetCRS)) {
var transformer = new HMCrsTransformer(outRaster.getCrs(), targetCRS);
transformer.setAcceptLenientDatumShift(true);
outRaster = transformer.transform(outRaster);
}

coverage = outRaster.buildCoverage();
}
CoordinateReferenceSystem targetCRS = CRS.decode("EPSG:4326");
if (bandIndex != null) { // Which means theat it's a Multi Band COG
coverage = (GridCoverage2D) Operations.DEFAULT.selectSampleDimension(coverage, new int[]{bandIndex});
}
if (!CRS.equalsIgnoreMetadata(
coverage.getCoordinateReferenceSystem(),
targetCRS)) {

coverage = (GridCoverage2D) Operations.DEFAULT.resample(
coverage,
targetCRS);
}
} catch (Exception e) {
throw new KlabResourceAccessException(
"Cannot build output for static collection " + collectionId + ". Reason: " + e.getLocalizedMessage());
Expand All @@ -413,7 +414,7 @@ public boolean test(HMStacAsset asset) { // Assuming for now that "eo:bands" wou
throw new KlabResourceAccessException("Cannot access to STAC collection " + collectionUrl); // Fail
// fast
}

IObservable targetSemantics = scope.getTargetArtifact() instanceof Observation
? ((Observation) scope.getTargetArtifact()).getObservable()
: null;
Expand Down Expand Up @@ -505,21 +506,20 @@ public boolean test(HMStacAsset asset) { // Assuming for now that "eo:bands" wou
scope.getMonitor().error("No STAC assets were found. Please check the spatial/temporal coverage of the resource");
throw new KlabIllegalStateException("No STAC assets were found. Please check the spatial/temporal coverage of the resource");
}
CoordinateReferenceSystem targetCRS = HMCrsRegistry.INSTANCE.getCrs("4326");
if (!HMCrsRegistry.crsEquals(outRaster.getCrs(),targetCRS)) {
var transformer = new HMCrsTransformer(outRaster.getCrs(), targetCRS);
transformer.setAcceptLenientDatumShift(true);
outRaster = transformer.transform(outRaster);
}

coverage = outRaster.buildCoverage();
if (bandIndex != null) { // Which means theat it's a Multi Band COG
coverage = (GridCoverage2D) Operations.DEFAULT.selectSampleDimension(coverage, new int[]{bandIndex});
}
CoordinateReferenceSystem targetCRS = CRS.decode("EPSG:4326");

manager.close();
encoder = new RasterEncoder();
if (!CRS.equalsIgnoreMetadata(
coverage.getCoordinateReferenceSystem(),
targetCRS)) {

coverage = (GridCoverage2D) Operations.DEFAULT.resample(
coverage,
targetCRS);
}
((RasterEncoder) encoder).encodeFromCoverage(resource, urnParameters, coverage, geometry, builder, scope);
} catch (Exception e) {
e.printStackTrace();
Expand Down
Loading