881781_ja-JP

liliconcom

2つの平面間の角度を計算する

皆さんこんにちは

平面1を定義するX、Y、Z軸の離散点が30個あり、平面2を定義するX、Y、Z軸の離散点が30個あります。平面1と平面2の間の角度を計算する必要があります。ChatGPTには以下のJSLスクリプトが付属していますが、行番号33で止まってしまいます。このスクリプトを実行するために、添付ファイル（*.jmpと*.jsl）を確認していただけると助かります。よろしくお願いいたします。

// 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 ) )
);

自動化とスクリプト基本的なデータ分析とモデリングRe: 2つの平面間の角度を計算する

スクリプトの残りの部分についてはわかりませんが、行 33 は JSL で書かれていません。

n1 = [a1, b1_, -1];
n2 = [a2, b2_, -1];

おそらくもっとこうすべきだろう

n1 = a1 ||  b1_ || -1;
 or perhaps
n2 = a2 |/ b2_ |/ -1;

n1 の例では行ベクトルを作成し、n2 の例では列ベクトルを作成します。

Re: 2つの平面間の角度を計算する

また、ドット積もおそらく間違っているので、JMP の要素ごとの乗算が必要かどうかを調べて判断する必要があります。 * は行列乗算であり、必要な演算子ではない可能性があります。