Created based class FourierDescClass

Author: timatet

AlgorithmsLibrary/AlgorithmsLibrary.csproj

@@ -71,6 +71,7 @@
     <Compile Include="Features\TripletComputation.cs" />
     <Compile Include="FourierDescAlgm\FourierDescAlgm.cs" />
     <Compile Include="FourierDescAlgm\FourierDescClass.cs" />
+    <Compile Include="FourierDescAlgm\FourierDescFatherClass.cs" />
     <Compile Include="FourierDescAlgm\FourierNotClosedDescClass.cs" />
     <Compile Include="FourierDescAlgm\LineSegmentClass.cs" />
     <Compile Include="FourierDescAlgm\PolygonClass.cs" />

AlgorithmsLibrary/FourierDescAlgm/FourierDescClass.cs

@@ -6,15 +6,15 @@
 
 namespace AlgorithmsLibrary
 {
-    public class Fourier
+    public class Fourier : FourierDescFatherClass
     {
         //Набор точек
         private List<MapPoint> pointCol = null;
         //Резервный набор точек. Возмомжно не нужен....
         //private List<MapPoint> newPointCol = null;
 
         private long m_PointNum = 0;
-        private long nTerm = 0;
+        
 
         //Si
         private double[] m_accuDist = null;
@@ -27,7 +27,7 @@ public class Fourier
         private double[] Ay = null;
         private double[] By = null;
 
-        private double[] d = null;
+        
 
         private double[] Power = null;
 
@@ -108,7 +108,7 @@ public void CalculateAllValue()
             GetAllDist();
             GetFourierXparameter();
             GetFourierYparameter();
-            CalculateShapeVector();
+            CalculateShapeVector(true);
         }
 
         //public void GetAllDist()
@@ -342,28 +342,6 @@ public void GetSinglePt(long nTerm, double s, out MapPoint pt)
         //    }
         //}
 
-        public double[] CalculateShapeVector()
-        {
-            long n = nTerm + 1;
-            double[] D = new double[n];
-            d = new double[n - 1];
-            ratio = new double[n - 1];
-            double sum = 0;
-            for (int i = 1; i < n; i++)
-            {
-                double CX = Ax[i] + By[i];
-                double CY = Bx[i] - Ay[i];
-                D[i] = Math.Abs(Math.Sqrt(CX * CX + CY * CY));
-                sum = sum + D[i];
-            }
-            for (int i = 1; i < n; i++)
-            {
-                ratio[i - 1] = D[i] / sum;
-                d[i - 1] = D[i] / D[1];
-            }
-            return d;
-        }
-
         public double[] CalculatePower()
         {
             long n = nTerm + 1;
@@ -525,28 +503,6 @@ public double[] CalculateEntropy()
             return Entropy;
         }
 
-        public double[] CalculateEntropy2(long Number)
-        {
-            double[] Proportion = new double[Number];
-
-            double d_sum = default(double);
-
-            for (int i = 0; i < Proportion.Length; i++)
-            {
-                d_sum += d[i + 1];
-            }
-
-            for (int i = 0; i < Proportion.Length; i++)
-            {
-                Proportion[i] = d[i + 1] / d_sum;
-            }
-
-            double[] Entropy = new double[Proportion.Length];
-            for (int i = 0; i < Entropy.Length; i++)
-            {
-                Entropy[i] = -Proportion[i] * Math.Log(Proportion[i], 2);
-            }
-            return Entropy;
-        }
+        
     }
 }

AlgorithmsLibrary/FourierDescAlgm/FourierDescFatherClass.cs

@@ -0,0 +1,63 @@
+using System;
+using System.Collections.Generic;
+using System.Linq;
+using System.Text;
+using System.Threading.Tasks;
+
+namespace AlgorithmsLibrary
+{
+    public abstract class FourierDescFatherClass
+    {
+        public double[] d = null;
+        public long nTerm = 0;
+
+        public double[] CalculateShapeVector(bool closed)
+        {
+            long n = nTerm + 1;
+            double[] D = new double[n];
+            d = new double[n - 1];
+            ratio = new double[n - 1];
+            double sum = 0;
+            for (int i = 1; i < n; i++)
+            {
+                double CX = Ax[i] + By[i];
+                double CY = Bx[i] - Ay[i];
+                dobule t = Math.Sqrt(CX * CX + CY * CY);
+                D[i] = closed ? Math.Abs(t) : t;
+                sum = sum + D[i];
+            }
+            for (int i = 1; i < n; i++)
+            {
+                ratio[i - 1] = D[i] / sum;
+                d[i - 1] = D[i] / D[1];
+            }
+            return d;
+        }
+        public double[] CalculateEntropy2(long Number)
+        {
+            double[] Proportion = new double[Number];
+
+            double d_sum = default(double);
+
+            for (int i = 0; i < Proportion.Length; i++)
+            {
+                d_sum += d[i + 1];
+            }
+
+            for (int i = 0; i < Proportion.Length; i++)
+            {
+                Proportion[i] = d[i + 1] / d_sum;
+            }
+
+            double[] Entropy = new double[Proportion.Length];
+            for (int i = 0; i < Entropy.Length; i++)
+            {
+                Entropy[i] = -Proportion[i] * Math.Log(Proportion[i], 2);
+            }
+            return Entropy;
+        }
+
+        
+
+    }
+}

AlgorithmsLibrary/FourierDescAlgm/FourierNotClosedDescClass.cs

@@ -6,34 +6,28 @@
 
 namespace AlgorithmsLibrary
 {
-    public class Fourier_NotClosed
+    public class Fourier_NotClosed : FourierDescFatherClass
     {
         #region After
-        //点集
         private List<MapPoint> pointCol = null;
-        //还原点集？(截至2017-11-7未使用)
         //private IPointCollection newPointCol = null;
-        //点的数量
+
         private long m_PointNum = 0;
-        //算子数量
-        private long nTerm = 0;
+
         //Si
         private double[] m_accuDist = null;
-        //点i-1到点Si的距离 
         private double[] m_dis_betPoint = null;
-        //总长度
+        
         private double m_totalS = 0.0;
-        //Ax,Bx,Ay,By系数数组
+      
         private double[] Ax = null;
         private double[] Bx = null;
         private double[] Ay = null;
         private double[] By = null;
-        //傅里叶算子
-        private double[] d = null;
-        //每个算子所占比例
+        
         private double[] ratio = null;
 
-        //清除距离过近的点
+        
         public Fourier_NotClosed(List<MapPoint> pointcollection, long nTerm)
         {
             this.nTerm = nTerm;
@@ -105,7 +99,7 @@ public void CalculateAllValue()
             GetAllDist();
             GetFourierXparameter();
             GetFourierYparameter();
-            CalculateShapeVector();
+            CalculateShapeVector(false);
         }
 
         //返回总长度
@@ -337,28 +331,6 @@ public void GetSinglePt(long nTerm, double s, out MapPoint pt)
         //    }
         //}
 
-        public double[] CalculateShapeVector()
-        {
-            long n = nTerm + 1;
-            double[] D = new double[n];
-            d = new double[n - 1];
-            ratio = new double[n - 1];
-            double sum = 0;
-            for (int i = 1; i < n; i++)
-            {
-                double CX = Ax[i] + By[i];
-                double CY = Bx[i] - Ay[i];
-                D[i] = Math.Sqrt(CX * CX + CY * CY);
-                sum = sum + D[i];
-            }
-            for (int i = 1; i < n; i++)
-            {
-                ratio[i - 1] = D[i] / sum;
-                d[i - 1] = D[i] / D[1];
-            }
-            return d;
-        }
-
         public double[,] GetRecoveryPoints(long FittingPointNumber, long Fouriers)
         {
             double x = Ax[0];
@@ -488,29 +460,7 @@ public double[] CalculateEntropy()
             return Entropy;
         }
 
-        public double[] CalculateEntropy2(long Number)
-        {
-            double[] Proportion = new double[Number];
-
-            double d_sum = default(double);
-
-            for (int i = 0; i < Proportion.Length; i++)
-            {
-                d_sum += d[i + 1];
-            }
-
-            for (int i = 0; i < Proportion.Length; i++)
-            {
-                Proportion[i] = d[i + 1] / d_sum;
-            }
-
-            double[] Entropy = new double[Proportion.Length];
-            for (int i = 0; i < Entropy.Length; i++)
-            {
-                Entropy[i] = -Proportion[i] * Math.Log(Proportion[i], 2);
-            }
-            return Entropy;
-        }
+        
         #endregion
     }
 }