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());
}