Martin Desruisseaux created SIS-548:
---------------------------------------
Summary: Unexpected wraparound during intersection of grid
geometries
Key: SIS-548
URL: https://issues.apache.org/jira/browse/SIS-548
Project: Spatial Information Systems
Issue Type: Improvement
Components: Coverage
Affects Versions: 1.2, 1.1
Reporter: Martin Desruisseaux
Assignee: Martin Desruisseaux
Fix For: 1.3
The following test case produces an intersection like below. This intersection
covers 360° of longitude, while a much smaller intersection was expected. This
is caused by an unexpected "wraparound" applied by
{{CoordinateOperationFinder}}.
{noformat}
Column: [-179 … 181] (361 cellules)
Row: [ 3 … 3] (1 cellules)
{noformat}
Test case:
{code:java}
public static void main(String[] args) {
final CoordinateReferenceSystem crs =
CommonCRS.WGS84.normalizedGeographic();
final AffineTransform2D gridToCrs = new AffineTransform2D(1, 0, 0, -1,
-103.5, 52.5);
final GridGeometry tilingScheme = new GridGeometry(new GridExtent(1, 1),
PixelInCell.CELL_CENTER, gridToCrs, crs);
final AffineTransform2D gridToCrs2 = new
AffineTransform2D(2.7770063871146905E-4, 0, 0, -2.7770063871146905E-4,
77.00013885031936, 49.99986114968065);
final GridGeometry tileGrid = new GridGeometry(new GridExtent(3601, 3601),
PixelInCell.CELL_CENTER, gridToCrs2, crs);
GridExtent intersection =
tilingScheme.derive().clipping(GridClippingMode.NONE).rounding(GridRoundingMode.ENCLOSING).subgrid(tileGrid).getIntersection();
System.out.println(intersection);
}
{code}
This issue has been reported by Johann.
--
This message was sent by Atlassian Jira
(v8.20.7#820007)
