新增按格网添加计曲线注记和示坡线

加密后软件工具/CAD插件工具/源代码/HCTools/HCTools/SlopeL.cs

@@ -83,6 +83,8 @@ namespace HCTools
                 return;
             }
             drawSlpeLine();
+            MessageBox.Show("绘制完成。");
+            this.Close();
         }
 
         private void drawSlpeLine()
@@ -251,85 +253,80 @@ namespace HCTools
                 pline.GetClosestPointTo(pl.StartPoint, new Vector3d(0, 0, 1), false))
                 ).ToList();
             //判断离pl最近等高线与pl高程关系，决定示坡线方向
-            if (pllist_sorted[0].Elevation > pl.Elevation)
+            if(pllist_sorted[0].Elevation == pl.Elevation|| 
+                    ((pllist_sorted[0].StartPoint == pllist_sorted[0].EndPoint||pllist_sorted[0].Closed==true)
+                    && pllist_sorted[0].Length<pl.Length))
+            {
+                bool flag = inside;
+                int i = 0;
+                do
+                {
+                    if (pllist_sorted[i].Elevation > pl.Elevation&&
+                       ((pllist_sorted[i].StartPoint == pllist_sorted[i].EndPoint || pllist_sorted[0].Closed == true)
+                        && pllist_sorted[i].Length < pl.Length))
+                        inside = true;
+                    else
+                        inside = false;
+                    i++;
+                }
+                while (flag!=inside);
+            }
+            else if (pllist_sorted[0].Elevation > pl.Elevation)
                 inside = true;
             else
                 inside = false;
             return inside;
         }
 
-        private void markSlopeline(Polyline pl, bool inside,Database db,int blc)
+        private void markSlopeline(Polyline pl, bool inside, Database db, int blc)
         {
-            //拟合多段线上所有点获取一条多段线内近似最长连线
-            int n = pl.NumberOfVertices;
-            double sumX = 0, sumY = 0, sumXY = 0, sumX2 = 0;
-
-            for (int i = 0; i < n; i++)
-            {
-                Point3d point = pl.GetPoint3dAt(i);
-                sumX += point.X;
-                sumY += point.Y;
-                sumXY += point.X * point.Y;
-                sumX2 += point.X * point.X;
-            }
-            //斜率
-            double a = (n * sumXY - sumX * sumY) / (n * sumX2 - sumX * sumX);
-            double b = (sumY - a * sumX) / n;
-            double x_base1 = pl.Bounds.Value.MaxPoint.X;
-            double x_base2 = pl.Bounds.Value.MinPoint.X;
-            Point3d line_pt1 = new Point3d(x_base1, x_base1 * a + b, pl.Elevation);
-            Point3d line_pt2 = new Point3d(x_base2, x_base2 * a + b, pl.Elevation);
-            Line fitline = new Line(line_pt1,line_pt2);
-            Vector3d k = new Vector3d(1, a, 0);
-            //double xComponent = 1.0;
-            //double yComponent = b;
-
-            //// 归一化向量，使其成为单位向量。
-            //double length = Math.Sqrt(xComponent * xComponent + yComponent * yComponent);
-            //// 防止除以零
-            //if (Math.Abs(length) > double.Epsilon)
-            //{
-            //    xComponent /= length;
-            //    yComponent /= length;
-            //}
-            //else
-            //{
-            //    // 处理斜率可能为无穷大或未定义的情况（例如垂直线）
-            //    if (b > 0)
-            //    {
-            //        xComponent = 0;
-            //        yComponent = 1;
-            //    }
-            //    else if (b == 0)
-            //    {
-            //        xComponent = 1;
-            //        yComponent = 0;
-            //    }
-            //    else
-            //    {
-            //        xComponent = 0;
-            //        yComponent = -1;
-            //    }
-            //}
+            //绘制最小外接矩形
+            Polyline box = new Polyline();
+            Point3d ptMax = pl.Bounds.Value.MaxPoint;
+            Point3d ptMin = pl.Bounds.Value.MinPoint;
+            box.AddVertexAt(0, new Point2d(ptMin.X, ptMin.Y), 0, 0, 0);
+            box.AddVertexAt(1, new Point2d(ptMax.X, ptMin.Y), 0, 0, 0);
+            box.AddVertexAt(2, new Point2d(ptMax.X, ptMax.Y), 0, 0, 0);
+            box.AddVertexAt(3, new Point2d(ptMin.X, ptMax.Y), 0, 0, 0);
+            box.Closed = true;
 
-
-            //由于IntersectWith()函数在离原点太远的位置会有误差，所以将需要判断的闭合等高线移动到原点附近
             Vector3d acVec3d_to = pl.StartPoint.GetVectorTo(new Point3d(0, 0, 0));
             var pl_clone = (Polyline)pl.Clone();
-            var xline_clone = (Line)fitline.Clone();
+            var box_clone = (Polyline)box.Clone();
 
             pl_clone.TransformBy(Matrix3d.Displacement(acVec3d_to));
-            xline_clone.TransformBy(Matrix3d.Displacement(acVec3d_to));
+            box_clone.TransformBy(Matrix3d.Displacement(acVec3d_to));
 
-            //Vector3d acVec3d_back = pl_clone.StartPoint.GetVectorTo(pl.StartPoint);
             Point3dCollection ptcoll_inter = new Point3dCollection();
+            Point3dCollection ptcoll_moveback = new Point3dCollection();
             Point3dCollection ptcoll = new Point3dCollection();
-            pl_clone.IntersectWith(xline_clone, Intersect.OnBothOperands, ptcoll_inter, IntPtr.Zero, IntPtr.Zero);
+            var plane = new Plane(Point3d.Origin, Vector3d.ZAxis);
+            pl_clone.IntersectWith(box_clone, Intersect.OnBothOperands, plane, ptcoll_inter, IntPtr.Zero, IntPtr.Zero);
             foreach (Point3d pt in ptcoll_inter)
             {
                 Point3d pt_new = pl.GetPointAtParameter(pl_clone.GetParameterAtPoint(pt));
-                ptcoll.Add(pt_new);
+                ptcoll_moveback.Add(pt_new);
+            }
+
+            Point3d pt1 = ptcoll_moveback[0];
+            Point3d pt2 = ptcoll_moveback[1];
+            double maxDistance = pt1.DistanceTo(pt2);
+            for (int i = 0; i < ptcoll_moveback.Count - 1; i++)
+            {
+                for (int j = i + 1; j < ptcoll_moveback.Count; j++)
+                {
+                    double distance = ptcoll_moveback[i].DistanceTo(ptcoll_moveback[j]);
+                    if (distance > maxDistance)
+                    {
+                        maxDistance = distance;
+                        pt1 = ptcoll_moveback[i];
+                        pt2 = ptcoll_moveback[j];
+                    }
+                }
             }
+            ptcoll.Add(pt1);
+            ptcoll.Add(pt2);
+            var k = (ptcoll[0] - ptcoll[1]).GetNormal();
             double spxlength = 1.6;
             if (blc == 500)
                 spxlength = spxlength * ((double)500 / 2000);