internal static double SmCalculator(EdgeWithScores eachEdge)
{
var edgesOfFace1 = new List <Edge>(eachEdge.Edge.OwnedFace.Edges);
var edgesOfFace2 = new List <Edge>(eachEdge.Edge.OtherFace.Edges);
edgesOfFace1 = edgesOfFace1.OrderBy(a => a.Length).ToList();
edgesOfFace2 = edgesOfFace2.OrderBy(a => a.Length).ToList();
double smallArea, largeArea;
if (eachEdge.Edge.OwnedFace.Area >= eachEdge.Edge.OtherFace.Area)
{
largeArea = eachEdge.Edge.OwnedFace.Area;
smallArea = eachEdge.Edge.OtherFace.Area;
}
else
{
largeArea = eachEdge.Edge.OtherFace.Area;
smallArea = eachEdge.Edge.OwnedFace.Area;
}
var r11 = edgesOfFace1[0].Length / edgesOfFace1[1].Length;
var r12 = edgesOfFace1[0].Length / edgesOfFace1[2].Length;
var r13 = edgesOfFace1[1].Length / edgesOfFace1[2].Length;
var r21 = edgesOfFace2[0].Length / edgesOfFace2[1].Length;
var r22 = edgesOfFace2[0].Length / edgesOfFace2[2].Length;
var r23 = edgesOfFace2[1].Length / edgesOfFace2[2].Length;
var similarity = Math.Abs(r11 - r21) + Math.Abs(r12 - r22) + Math.Abs(r13 - r23); // cannot exceed 3
var areaSimilarity = 3 * Math.Abs(1 - (smallArea / largeArea));
var SM = similarity + areaSimilarity;
return(SM);
}
internal static double AbnCalculator(EdgeWithScores eachEdge)
{
double ABN;
if (eachEdge.Edge.InternalAngle <= Math.PI)
{
ABN = (Math.PI - eachEdge.Edge.InternalAngle) * 180 / Math.PI;
}
else
{
ABN = eachEdge.Edge.InternalAngle * 180 / Math.PI;
}
if (ABN >= 180)
{
ABN -= 180;
}
if (ABN > 179.5)
{
ABN = 180 - ABN;
}
if (Double.IsNaN(ABN))
{
var eee = eachEdge.Edge.OwnedFace.Normal.dotProduct(eachEdge.Edge.OtherFace.Normal, 3);
if (eee > 1)
{
eee = 1;
}
ABN = Math.Acos(eee);
}
return(ABN);
}
private static void EdgeFuzzyClassification(EdgeWithScores e)
{
var ABN = AbnCalculator(e);
var MCM = McmCalculator(e);
var SM = SmCalculator(e);
var ABNid = CatAndProbFinder(ABN, listOfLimitsABN);
var MCMid = CatAndProbFinder(MCM, listOfLimitsMCM);
var SMid = CatAndProbFinder(SM, listOfLimitsSM);
if (ABNid.Count == 2 && ABNid[0].SequenceEqual(ABNid[1]))
{
ABNid.RemoveAt(0);
}
e.CatProb = new Dictionary <int, double>();
foreach (var ABNprobs in ABNid)
{
foreach (var MCMProbs in MCMid)
{
foreach (var SMProbs in SMid)
{
double Prob;
int group = EdgeClassifier2(ABNprobs, MCMProbs, SMProbs, edgeRules, out Prob);
if (!e.CatProb.Keys.Contains(@group))
{
e.CatProb.Add(@group, Prob);
}
else if (e.CatProb[@group] < Prob)
{
e.CatProb[@group] = Prob;
}
}
}
}
}
internal static double McmCalculator(EdgeWithScores eachEdge)
{
var cenMass1 = eachEdge.Edge.OwnedFace.Center;
var cenMass2 = eachEdge.Edge.OtherFace.Center;
var vector1 = new[] { cenMass1[0] - eachEdge.Edge.From.Position[0], cenMass1[1] - eachEdge.Edge.From.Position[1], cenMass1[2] - eachEdge.Edge.From.Position[2] };
var vector2 = new[] { cenMass2[0] - eachEdge.Edge.From.Position[0], cenMass2[1] - eachEdge.Edge.From.Position[1], cenMass2[2] - eachEdge.Edge.From.Position[2] };
var distance1 = eachEdge.Edge.Vector.normalize(3).dotProduct(vector1, 3);
var distance2 = eachEdge.Edge.Vector.normalize(3).dotProduct(vector2, 3);
//Mapped Center of Mass
var MCM = (Math.Abs(distance1 - distance2)) / eachEdge.Edge.Length;
return(MCM);
}
