public override void doTransform(PointD p0, PointD p1) { JamaMatrix jamaMatrix = new JamaMatrix(1, 2); jamaMatrix.SetElement(0, 0, p0.x); jamaMatrix.SetElement(0, 1, p0.y); JamaMatrix jamaMatrix2 = IPolyPointTransformer.Polynomialize(jamaMatrix, this.polynomialDegree).times(this.matrix); p1.x = jamaMatrix2.GetElement(0, 0); p1.y = jamaMatrix2.GetElement(1, 0); }
private static IPolyPointTransformer getPolyPointTransformer(JamaMatrix sourcePoints, JamaMatrix destPoints, int polynomialDegree) { JamaMatrix am = IPolyPointTransformer.Polynomialize(destPoints, polynomialDegree); JamaMatrix matrix = PolynomialImageTransformer.SVDSolveApply(am, PolynomialImageTransformer.PointUnroll(sourcePoints)); switch (polynomialDegree) { case 1: return(new FastPoly1PointTransformer(matrix)); case 2: return(new FastPoly2PointTransformer(matrix)); default: return(new SlowGeneralPolyPointTransformer(polynomialDegree, matrix)); } }
public static void TestFunc() { double[][] array = new double[4][]; double[][] arg_15_0 = array; int arg_15_1 = 0; double[] array2 = new double[2]; arg_15_0[arg_15_1] = array2; array[1] = new double[] { 1.02, 0.93 }; array[2] = new double[] { 0.0, 1.0 }; double[][] arg_73_0 = array; int arg_73_1 = 3; double[] array3 = new double[2]; array3[0] = 1.0; arg_73_0[arg_73_1] = array3; double[][] a = array; JamaMatrix jamaMatrix = new JamaMatrix(a); PolynomialImageTransformer.RegularUnitGrid(4); JamaMatrix jamaMatrix2 = (JamaMatrix)jamaMatrix.Clone(); jamaMatrix2.SetElement(1, 0, 1.0); jamaMatrix2.SetElement(1, 1, 1.0); JamaMatrix jamaMatrix3 = IPolyPointTransformer.Polynomialize(jamaMatrix, 2); JamaMatrix jamaMatrix4 = PolynomialImageTransformer.SVDSolveApply(jamaMatrix3, PolynomialImageTransformer.PointUnroll(jamaMatrix2)); D.Say(0, "polyTransform:\n" + jamaMatrix4.ToString()); JamaMatrix unrolledVector = jamaMatrix3.times(jamaMatrix4); D.Say(0, "testSolution:\n" + PolynomialImageTransformer.PointRoll(unrolledVector, 2).ToString()); PolynomialImageTransformer.getPolyPointTransformer(jamaMatrix, jamaMatrix2, 2); LatLon p = new LatLon(0.93, 1.02); D.Say(0, "Invert test:\n" + PolynomialImageTransformer.getApproximateInverterPolyPointTransformer(jamaMatrix, jamaMatrix2, 2).getTransformedPoint(p).ToString()); }