Skip to content

Limit Spatial Bounds #104

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

Open
AM1729 wants to merge 4 commits into
base: master
Choose a base branch
from
Open

Limit Spatial Bounds #104

Changes from all commits
Commits
Show all changes
4 commits
Select commit Hold shift + click to select a range
7d3438d
remove ws
AM1729 Apr 13, 2026
a4418f2
fixes for GeographicCRS, checks AxisDir
AM1729 Apr 13, 2026
3b8160d
use getDimension other then assuming 2
AM1729 Apr 13, 2026
b1c2e66
use equals instead of ==
AM1729 Apr 13, 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
17 changes: 14 additions & 3 deletions gears/src/main/java/org/hortonmachine/gears/utils/coverage/CoverageUtilities.java
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -55,6 +55,7 @@
import org.eclipse.imagen.iterator.RandomIterFactory;
import org.eclipse.imagen.iterator.WritableRandomIter;
import org.eclipse.imagen.media.range.NoDataContainer;
import org.geotools.api.referencing.cs.AxisDirection;
import org.geotools.coverage.CoverageFactoryFinder;
import org.geotools.coverage.GridSampleDimension;
import org.geotools.coverage.grid.GridCoordinates2D;
Expand Down Expand Up @@ -481,7 +482,6 @@ public static RegionMap getRegionParamsFromGridCoverage( GridCoverage2D gridCove
double[] westSouth = lowerCorner.getCoordinate();
Position upperCorner = envelope.getUpperCorner();
double[] eastNorth = upperCorner.getCoordinate();

GridGeometry2D gridGeometry = gridCoverage.getGridGeometry();
GridEnvelope2D gridRange = gridGeometry.getGridRange2D();
int height = gridRange.height;
Expand All @@ -494,8 +494,19 @@ public static RegionMap getRegionParamsFromGridCoverage( GridCoverage2D gridCove
double yRes = XAffineTransform.getScaleY0(gridToCRS);

CoordinateReferenceSystem crs = gridCoverage.getCoordinateReferenceSystem();
// check if crs is geographic
if (crs != null && crs instanceof GeographicCRS) {
// check if crs is geographic, angular
if (crs instanceof GeographicCRS) {
// limit to bounds check (lat, lon) or (lon, lat)
for (int i = 0; i < upperCorner.getDimension(); i++) {
AxisDirection dir = crs.getCoordinateSystem().getAxis(i).getDirection();
if (AxisDirection.NORTH.equals(dir) || AxisDirection.SOUTH.equals(dir)) {
westSouth[i] = Math.max(-90, Math.min(90, westSouth[i]));
eastNorth[i] = Math.max(-90, Math.min(90, eastNorth[i]));
} else if (AxisDirection.EAST.equals(dir) || AxisDirection.WEST.equals(dir)) {
westSouth[i] = Math.max(-180, Math.min(180, westSouth[i]));
eastNorth[i] = Math.max(-180, Math.min(180, eastNorth[i]));
}
}
// convert to metric
GeodeticCalculator gc = new GeodeticCalculator(crs);
gc.setStartingGeographicPoint(westSouth[0], westSouth[1]);
Expand Down