[
https://issues.apache.org/jira/browse/SIS-548?page=com.atlassian.jira.plugin.system.issuetabpanels:all-tabpanel
]
Martin Desruisseaux updated SIS-548:
------------------------------------
Description:
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) {
CoordinateReferenceSystem crs = CommonCRS.WGS84.normalizedGeographic();
AffineTransform2D gridToCrs = new AffineTransform2D(1, 0, 0, -1, -103.5,
52.5);
GridGeometry tilingScheme = new GridGeometry(new GridExtent(1, 1),
PixelInCell.CELL_CENTER, gridToCrs, crs);
AffineTransform2D gridToCrs2 = new AffineTransform2D(
2.7770063871146905E-4, 0,
0, -2.7770063871146905E-4,
77.00013885031936,
49.99986114968065);
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.
was:
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.
> 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.1, 1.2
> Reporter: Martin Desruisseaux
> Assignee: Martin Desruisseaux
> Priority: Major
> 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) {
> CoordinateReferenceSystem crs = CommonCRS.WGS84.normalizedGeographic();
> AffineTransform2D gridToCrs = new AffineTransform2D(1, 0, 0, -1, -103.5,
> 52.5);
> GridGeometry tilingScheme = new GridGeometry(new GridExtent(1, 1),
> PixelInCell.CELL_CENTER, gridToCrs, crs);
> AffineTransform2D gridToCrs2 = new AffineTransform2D(
> 2.7770063871146905E-4, 0,
> 0, -2.7770063871146905E-4,
> 77.00013885031936,
> 49.99986114968065);
> 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)
