Hola amigos,
Tengo 30 puntos discretos en X, Y, Z que definen el plano 1, y otros 30 puntos en X, Y, Z que definen el plano 2. Necesito calcular el ángulo entre los planos 1 y 2. ChatGPT proporciona el siguiente script JSL, pero se queda atascado en la fila 33. ¿Podrían ayudarme a revisar los dos archivos adjuntos (*.jmp y *.jsl) para que este script funcione? Muchas gracias.
// Use current data table
dt = Current Data Table();
// Subset for Plane1 and Plane2
dt1 = dt << Subset( Where( :Plane == "Plane1" ), Output Table Name( "Plane1" ) );
dt2 = dt << Subset( Where( :Plane == "Plane2" ), Output Table Name( "Plane2" ) );
// Fit plane 1: z = a1*x + b1*y + c1
b1 = dt1 << Fit Model(
Y( :Z ),
Effects( :X, :Y ),
Personality( "Standard Least Squares" ),
Emphasis( "Effect Leverage" ),
Run( "Fit Model" )
);
coeff1 = b1 << Get Estimates;
a1 = coeff1["X"];
b1_ = coeff1["Y"];
// Fit plane 2: z = a2*x + b2*y + c2
b2 = dt2 << Fit Model(
Y( :Z ),
Effects( :X, :Y ),
Personality( "Standard Least Squares" ),
Emphasis( "Effect Leverage" ),
Run( "Fit Model" )
);
coeff2 = b2 << Get Estimates;
a2 = coeff2["X"];
b2_ = coeff2["Y"];
// Normal vectors: [a, b, -1]
n1 = [a1, b1_, -1];
n2 = [a2, b2_, -1];
// Compute dot product and magnitudes
dotProd = Sum( n1 * n2 );
mag1 = Sqrt( Sum( n1 ^ 2 ) );
mag2 = Sqrt( Sum( n2 ^ 2 ) );
// Compute angle in degrees
cosTheta = dotProd / (mag1 * mag2);
thetaRad = ACos( cosTheta );
thetaDeg = thetaRad * (180 / Pi());
// Display result
New Window( "Angle Between Planes",
Text Box( "Angle (degrees): " || Char( thetaDeg ) )
);
No tengo idea del resto del script, pero la línea 33 no está escrita en JSL.
n1 = [a1, b1_, -1];
n2 = [a2, b2_, -1];
Probablemente debería ser más como
n1 = a1 || b1_ || -1;
or perhaps
n2 = a2 |/ b2_ |/ -1;
El ejemplo n1 debería crear un vector de fila y el n2 un vector de columna.
Y el producto escalar probablemente también esté mal, debes buscarlo y averiguar si la multiplicación elemento por elemento de JMP es lo que necesitas; el * es una multiplicación de matrices que podría no ser el operador que necesitas.
