public void testInterpolate() { MCoordinate mc0NaN = MCoordinate.create2d(0.0, 0.0); MCoordinate mc0 = MCoordinate.create2dWithMeasure(0.0, 0.0, 0.0); MCoordinate mc2NaN = MCoordinate.create2d(2.0, 0.0); MCoordinate mc5NaN = MCoordinate.create2d(5.0, 0.0); MCoordinate mc10NaN = MCoordinate.create2d(10.0, 0.0); MCoordinate mc10 = MCoordinate.create2dWithMeasure(10.0, 0.0, 10.0); // Internal coordinate measures are not defined, outer measures are // 0-10, total 2d length is 10 MLineString line = mgeomFactory.createMLineString(new MCoordinate[] {mc0, mc2NaN, mc5NaN, mc10}); MLineString lineBeginNaN = mgeomFactory.createMLineString(new MCoordinate[] {mc0NaN, mc2NaN, mc5NaN, mc10}); MLineString lineEndNaN = mgeomFactory.createMLineString(new MCoordinate[] {mc0, mc2NaN, mc5NaN, mc10NaN}); assertTrue(DoubleComparator.equals(line.getLength(), 10)); assertTrue(DoubleComparator.equals(lineBeginNaN.getLength(), 10)); assertTrue(DoubleComparator.equals(lineEndNaN.getLength(), 10)); line.interpolate(mc0.m, mc10.m); lineBeginNaN.interpolate(mc0.m, mc10.m); lineEndNaN.interpolate(mc0.m, mc10.m); assertTrue(line.isMonotone(false)); assertTrue(line.isMonotone(true)); assertTrue(lineBeginNaN.isMonotone(false)); assertTrue(lineBeginNaN.isMonotone(true)); assertTrue(lineEndNaN.isMonotone(false)); assertTrue(lineEndNaN.isMonotone(true)); double[] expectedM = new double[] {mc0.m, 2.0, 5.0, mc10.m}; for (int i = 0; i < expectedM.length; i++) { double actualMLine = line.getCoordinateSequence().getOrdinate(i, CoordinateSequence.M); double actualBeginNaN = lineBeginNaN.getCoordinateSequence().getOrdinate(i, CoordinateSequence.M); double actualEndNaN = lineEndNaN.getCoordinateSequence().getOrdinate(i, CoordinateSequence.M); assertTrue(DoubleComparator.equals(expectedM[i], actualMLine)); assertTrue(DoubleComparator.equals(expectedM[i], actualBeginNaN)); assertTrue(DoubleComparator.equals(expectedM[i], actualEndNaN)); } // Test Continuous case by interpolating with begin and end measures // equal double continuousMeasure = 0.0D; line.interpolate(continuousMeasure, continuousMeasure); double[] measures = line.getMeasures(); for (int i = 0; i < measures.length; i++) assertTrue(DoubleComparator.equals(measures[i], continuousMeasure)); }
public void testGetClosestPoint() { try { if (!arbitraryLine.isMonotone(false)) { Coordinate mc = arbitraryLine.getClosestPoint(new Coordinate(1.0, 2.0), 0.1); assertTrue(false); // should never evaluate this } } catch (Exception e) { assertTrue( ((MGeometryException) e).getType() == MGeometryException.OPERATION_REQUIRES_MONOTONE); } try { // check reaction on null string MCoordinate mc = nullLine.getClosestPoint(new Coordinate(0.0, 1.0), 1.0); assertNull(mc); // must return the very same coordinate if the coordinate is a // coordinate of the line arbitraryLine.measureOnLength(false); int selp = (int) (arbitraryLine.getNumPoints() / 2); MCoordinate mcexp = (MCoordinate) arbitraryLine.getCoordinateN(selp); MCoordinate mctest = arbitraryLine.getClosestPoint(mcexp, 1); assertEquals(mcexp, mctest); // must not return a point that is beyond the tolerance mctest = controlledLine.getClosestPoint(new Coordinate(20.0, 20, 0), 1.0); assertNull(mctest); // check for cases of circular MGeometry: lowest measure should be // return. ringLine.measureOnLength(false); assertTrue(ringLine.isRing()); assertTrue(ringLine.isMonotone(false)); assertTrue(ringLine.getMeasureDirection() == MGeometry.INCREASING); MCoordinate expCo = MCoordinate.create2dWithMeasure(0.0, 0.0, 0.0); MCoordinate testCo = ringLine.getClosestPoint(expCo, 0.1); assertTrue(DoubleComparator.equals(testCo.m, expCo.m)); ringLine.reverseMeasures(); testCo = ringLine.getClosestPoint(expCo, 0.1); assertTrue(DoubleComparator.equals(testCo.m, expCo.m)); ringLine.measureOnLength(false); int n = ringLine.getNumPoints() - 1; ringLine.setMeasureAtIndex(n, 100.0); ringLine.setMeasureAtIndex(0, 0.0); testCo = ringLine.getClosestPoint(expCo, 0.001); assertTrue(DoubleComparator.equals(testCo.m, 0.0)); // get two neighbouring points along the arbitraryline arbitraryLine.measureOnLength(false); int elem1Indx = (int) (Math.random() * (arbitraryLine.getNumPoints() - 1)); int elem2Indx = 0; if (elem1Indx == arbitraryLine.getNumPoints() - 1) { elem2Indx = elem1Indx - 1; } else { elem2Indx = elem1Indx + 1; } // testsuite-suite whether a coordinate between these two returns exactly MCoordinate mco1 = (MCoordinate) arbitraryLine.getCoordinateN(elem1Indx); MCoordinate mco2 = (MCoordinate) arbitraryLine.getCoordinateN(elem2Indx); double d = mco1.distance(mco2); double offset = Math.random(); mcexp = MCoordinate.create2dWithMeasure( mco1.x + offset * (mco2.x - mco1.x), mco1.y + offset * (mco2.y - mco1.y), 0.0); mctest = arbitraryLine.getClosestPoint(mcexp, d); mcexp.m = mco1.m + offset * (mco2.m - mco1.m); assertEquals(mcexp.x, mctest.x, 0.001); assertEquals(mcexp.y, mctest.y, 0.001); assertEquals(mcexp.z, mctest.z, 0.001); double delta = Math.random(); MCoordinate mcin = MCoordinate.create2dWithMeasure( mco1.x + offset * (mco2.x - mco1.x) + delta, mco1.y + offset * (mco2.y - mco1.y) + delta, 0.0); // returned point is on the line mctest = arbitraryLine.getClosestPoint(mcin, d); assertEquals(mcin.x, mctest.x, delta * Math.sqrt(2)); assertEquals(mcin.y, mctest.y, delta * Math.sqrt(2)); } catch (Exception e) { e.printStackTrace(); assertTrue(false); // should never reach this point } }