public HatchNestingNode(CompositeCurve2d boundary, Extents2d extends,
Point2dCollection ptsColl, List <double> ptParams, Hatch hatch) : base(ptsColl)
{
Boundary = boundary;
Extents = extends;
//base.Point2dCollection = ptsColl;
this.ptParams = ptParams;
Hatch = hatch;
}
private void CreatePositionCurve(List <Point2d> polyPts)
{
List <Curve2d> segments = new List <Curve2d>();
for (int i = 0; i < polyPts.Count - 1; i++)
{
Point2d p0 = polyPts[i];
Point2d p1 = polyPts[i + 1];
LineSegment2d segment = new LineSegment2d(p0, p1);
segments.Add(segment);
}
PositionCurve = new CompositeCurve2d(segments.ToArray());
}
/// <summary> 构造函数 </summary>
/// <param name="docMdf"></param>
/// <param name="roadCurve"></param>
/// <param name="groundCurve"></param>
public LongitudinalSection(DocumentModifier docMdf, CompositeCurve3d roadCurve, CompositeCurve3d groundCurve)
{
_docMdf = docMdf;
RoadCurve2d = roadCurve.Get2dLinearCurve();
GroundCurve2d = groundCurve.Get2dLinearCurve();
//
StartStation = RoadCurve2d.StartPoint.X;
EndStation = RoadCurve2d.EndPoint.X;
// 求交点
Intersects = new CurveCurveIntersector2d(RoadCurve2d, GroundCurve2d);
GetIntersects();
//_road = Curve.CreateFromGeCurve(roadCurve) as Polyline;
//_ground = Curve.CreateFromGeCurve(groundCurve) as Polyline;
}
/// <summary>
/// 计算道路横断面的某一侧中，从路面中心到边坡外边缘的范围内，属于填方的区域在 AutoCAD 几何中的 X 范围
/// </summary>
/// <param name="sec"></param>
/// <param name="left"></param>
/// <param name="slp"> 某一侧边坡，其值可能为null，表示此侧没有边坡线 </param>
/// <param name="cGround"></param>
/// <param name="db"></param>
/// <param name="edgeXleft">此侧边坡的填方左边界</param>
/// <param name="edgeXright">此侧边坡的填方右边界</param>
/// <returns>如果没有填方区域，则返回 false </returns>
public bool GetFillSlopeXRange(bool left, SlopeLine slp, CompositeCurve2d cGround, Database db,
out double edgeXleft, out double edgeXright)
{
edgeXleft = 0.0;
edgeXright = 0.0;
var centerFill = IsCenterFill();
var slopeFill = slp == null || slp.XData.FillCut;
if (!centerFill && (slp == null || !slopeFill))
{
return(false);
}
// 确定进行搜索的左右边界：路基边缘（或边坡脚） 到 道路中线
double roadEdge = left ? LeftRoadEdge.X : RightRoadEdge.X;
double slopeEdge = roadEdge; // 边坡的坡脚的 X 值
if (slp != null && slp.XData.Slopes.Count > 0)
{
var data = slp.XData;
slopeEdge = data.Slopes[data.Slopes.Count - 1].OuterPoint.X;
}
if (centerFill && slopeFill)
{
// 道路中心与边坡均为填方
double edgeX1 = CenterX;
edgeXleft = Math.Min(edgeX1, slopeEdge);
edgeXright = Math.Max(edgeX1, slopeEdge);
}
else
{
// 说明 坡脚与道路中心这二者中有一个为挖方，另一个为填方
var roadSurfHandle = left ? LeftRoadSurfaceHandle : RightRoadSurfaceHandle;
var roadSurf = roadSurfHandle.GetDBObject <Polyline>(db);
var cRoad = roadSurf.Get2dLinearCurve();
var inters = new CurveCurveIntersector2d(cRoad, cGround);
double iX;
if (inters.NumberOfIntersectionPoints > 0)
{
iX = inters.GetIntersectionPoint(0).X;
}
else
{
// 这种情况极少会出现，但测试中确实会有，即自然地面线与挖方坡底边沟相交，而不与路面相交
iX = roadEdge;
}
// 自然地面与路面的交点
var roadWidth = Math.Abs((roadEdge - CenterX));
var innerRatio = Math.Abs((iX - CenterX) / roadWidth);
if ((centerFill && innerRatio > 0.5))
{
// 靠道路中心为填方，边坡为挖方
if (left)
{
edgeXleft = iX;
edgeXright = CenterX;
}
else
{
edgeXleft = CenterX;
edgeXright = iX;
}
}
else if (!centerFill && innerRatio <= 0.5)
{
// 靠道路中心为挖方，边坡为填方
if (left)
{
edgeXleft = slopeEdge;
edgeXright = iX;
}
else
{
edgeXleft = iX;
edgeXright = slopeEdge;
}
}
else
{
// 填方区域太小
return(false);
}
}
return(true);
}
public static double GetHatchArea([NotNull] this Hatch hatch)
{
double area = 0;
try
{
area = hatch.Area;
}
catch
{
var nLoop = hatch.NumberOfLoops;
for (var i = 0; i < nLoop; i++)
{
double looparea = 0;
var loopType = (int)hatch.LoopTypeAt(i);
if ((loopType & (int)HatchLoopTypes.Polyline) > 0)
{
var hatchLoop = hatch.GetLoopAt(i);
var bulgeVertex = hatchLoop.Polyline;
using var pPoly = new Polyline(bulgeVertex.Count);
for (var j = 0; j < bulgeVertex.Count; j++)
{
pPoly.AddVertexAt(j, bulgeVertex[j].Vertex, bulgeVertex[j].Bulge, 0, 0);
}
pPoly.Closed = (loopType & (int)HatchLoopTypes.NotClosed) == 0;
looparea = pPoly.Area;
if ((loopType & (int)HatchLoopTypes.External) > 0)
{
area += Math.Abs(looparea);
}
else
{
area -= Math.Abs(looparea);
}
}
else
{
var hatchLoop = hatch.GetLoopAt(i);
var cur2ds = new Curve2d[hatchLoop.Curves.Count];
hatchLoop.Curves.CopyTo(cur2ds, 0);
using var compCurve = new CompositeCurve2d(cur2ds);
var interval = compCurve.GetInterval();
double dMin = interval.GetBounds()[0], dMax = interval.GetBounds()[1];
if (Math.Abs(dMax - dMin) > 1e-6)
{
try
{
looparea = compCurve.GetArea(dMin, dMax);
if ((loopType & (int)HatchLoopTypes.External) > 0)
{
area += Math.Abs(looparea);
}
else
{
area -= Math.Abs(looparea);
}
}
catch
{
// Разбиваем кривую на 1000000 точек. Надеюсь, что такой точности
// будет достаточно.
var pts = compCurve.GetSamplePoints(1000);
var np = pts.Length;
for (var j = 0; j < np; j++)
{
looparea += 0.5 * pts[j].X * (pts[(j + 1) % np].Y - pts[(j + np - 1) % np].Y);
}
if ((loopType & (int)HatchLoopTypes.External) > 0)
{
area += Math.Abs(looparea);
}
else
{
area -= Math.Abs(looparea);
}
}
}
}
}
}
return(Math.Abs(area));
}
private Node BoundaryTraverse(Node currNode, List <Point2d> polygon, HashSet <Node> visitedThisTraverse)
{
if (visitedThisTraverse.Contains(currNode))
{
throw new Exception("Ошибка логики обхода полигона");
}
//Определить к какому контуру штриховки относится точка в которую пришли
//затем выбор направления обхода вдоль контура полигона так чтобы обходить внутри текущего диапазона
//(использовать прирощение параметра в 2 раза меньшее чем длина текущего диапазона и определять находится ли точка контура с данным прирощением в текущем диапазоне)
//обход вдоль контура до следующей точки пересечения в соответствии с коллекцией по параметру композитной кривой для текущего контура и выбранным направлением обхода
HatchNestingNode hnn = currNode.Boundary;
CompositeCurve2d boundaryCurve = currNode.Boundary.Boundary;
LinkedListNode <Node> currLLNode = currNode.BoundaryJunction;
//Соседние узлы должны браться как будто это циклический список!!!
LinkedListNode <Node> nextLLNode = Common.Utils.NextVertCycled(currLLNode); //следующий узел в направлении возрастания параметра кривой (возможен переход через 0)
LinkedListNode <Node> prevLLNode = Common.Utils.PreviousVertCycled(currLLNode); //следующий узел в направлении убывания параметра кривой (возможен переход через 0)
if (nextLLNode == null || prevLLNode == null || nextLLNode == currLLNode || prevLLNode == currLLNode)
{
throw new Exception("Ошибка логики обхода полигона");
}
double currParam = currNode.CurveParameter;
double nextParam = nextLLNode.Value.CurveParameter;
double maxParam = boundaryCurve.GetInterval().UpperBound;
double nextCheckParam = currParam < nextParam ? //в списке параметр идет по возрастанию
(currParam + nextParam) / 2
: //если следующий параметр меньше текущего, значит происходит переход через ноль
(currParam + (maxParam + nextParam - currParam) / 2) % maxParam; //TODO: ПРОВЕРИТЬ ПРАВИЛЬНОСТЬ!!!
//(как вариант - за средний параметр можно просто взять ноль (но нужно учесть, что 0 может быть одной из границ интервала))
Point2d nextCheckPt = boundaryCurve.EvaluatePoint(nextCheckParam);
Node result = null;
bool orderForward = true;//проход вдоль контура в направлении возрастания параметра
//Если nextCheckPt попадает точно на границу интервала, то считать, что она попала в интервал
//(ТОЛЬКО ЕСЛИ НЕ СОВПАДАЕТ ПО X С currNode!!!)
if ((start < /*=*/ nextCheckPt.X || Utils.LengthIsEquals(start, nextCheckPt.X)) &&
(end > /*=*/ nextCheckPt.X || Utils.LengthIsEquals(end, nextCheckPt.X)) &&
!Utils.LengthIsEquals(currNode.Location.X, nextCheckPt.X))
{
//if (start < nextCheckPt.X && end > nextCheckPt.X)
result = nextLLNode.Value;
}
else
{
result = prevLLNode.Value;
orderForward = false;//точки добавляются в полигон в порядке убывания параметра
}
//все вершины, лежащие внутри интервала должны перейти в полигон в нужном порядке
double intervalStartParam = orderForward ? currParam : result.CurveParameter;
double intervalEndParam = orderForward ? result.CurveParameter : currParam;
List <Point2d> addToPolygon = hnn.GetPointsInInterval(intervalStartParam, intervalEndParam);
if (!orderForward)
{
addToPolygon.Reverse();
}
polygon.Add(currNode.Location);
polygon.AddRange(addToPolygon);
visitedThisTraverse.Add(currNode);
return(result);
}
private void Create3dProfile(object arg)
{
try
{
if (doc != null)
{
List <ObjectId> toHighlight = new List <ObjectId>();
using (doc.LockDocument())
{
Editor ed = doc.Editor;
Database db = doc.Database;
using (Transaction tr = db.TransactionManager.StartTransaction())
{
//найти координату X начала профиля (крайнюю левую)
double minx = double.PositiveInfinity;
List <ObjectId> allSelected = new List <ObjectId>(soilHatchIds);
foreach (ObjectId id in allSelected)
{
Entity ent = null;
try
{
ent = (Entity)tr.GetObject(id, OpenMode.ForRead);
}
catch (Autodesk.AutoCAD.Runtime.Exception) { continue; }
Extents3d?ext = ent.Bounds;
if (ext != null)
{
Point3d minPt = ext.Value.MinPoint;
if (minx > minPt.X)
{
minx = minPt.X;
}
}
}
//Штриховки должны быть заранее раскиданы по слоям в соответствии с ИГЭ!
//пересчет всех точек штриховок в координаты:
//X - положение отностиельно начала профиля с учетом горизонтального масштаба профиля,
//Y - отметка расчитанная согласно базовой отметке с учетом вертикального масштаба профиля
Matrix2d transform =
new Matrix2d(new double[]
{
HorScaling, 0, 0,
0, VertScaling, 0,
0, 0, 1
})
* Matrix2d.Displacement(new Vector2d(-minx, -ElevBasePoint.Value.Y + ElevationInput / VertScaling));
C5.IntervalHeap <HatchEvent> eventQueue = new C5.IntervalHeap <HatchEvent>();
List <HatchData> allHatchData = new List <HatchData>();
List <Point2dCollection> selfintersectingLoops = new List <Point2dCollection>();
foreach (ObjectId id in soilHatchIds)
{
//получить все точки штриховок с учетом возможных дуг, сплайнов и проч
//Для каждой штриховки создается набор композитных кривых, состоящих из линейных сегментов
Hatch hatch = null;
try
{
hatch = (Hatch)tr.GetObject(id, OpenMode.ForRead);
}
catch (Autodesk.AutoCAD.Runtime.Exception) { continue; }
List <CompositeCurve2d> boundaries = new List <CompositeCurve2d>();
List <Extents2d> extends = new List <Extents2d>();
List <Point2dCollection> ptsCollList = new List <Point2dCollection>();
List <List <double> > ptParamsList = new List <List <double> >();
for (int i = 0; i < hatch.NumberOfLoops; i++)
{
HatchLoop hl = hatch.GetLoopAt(i);
if (!hl.LoopType.HasFlag(HatchLoopTypes.SelfIntersecting) &&
!hl.LoopType.HasFlag(HatchLoopTypes.Textbox) &&
!hl.LoopType.HasFlag(HatchLoopTypes.TextIsland) &&
!hl.LoopType.HasFlag(HatchLoopTypes.NotClosed))
{
List <Curve2d> curves = Utils.GetHatchLoopCurves(hl);
List <Curve2d> compositeCurveElems = new List <Curve2d>();
double _minx = double.PositiveInfinity;
double _miny = double.PositiveInfinity;
double _maxx = double.NegativeInfinity;
double _maxy = double.NegativeInfinity;
Action <Point2d> updateExtends = new Action <Point2d>(p =>
{
_minx = p.X < _minx ? p.X : _minx;
_miny = p.Y < _miny ? p.Y : _miny;
_maxx = p.X > _maxx ? p.X : _maxx;
_maxy = p.Y > _maxy ? p.Y : _maxy;
});
Point2dCollection ptsColl = new Point2dCollection();
List <double> ptParams = new List <double>();
double currParam = 0;
foreach (Curve2d c in curves)
{
if (!(c is LineSegment2d))
{
Interval interval = c.GetInterval();
PointOnCurve2d[] samplePts = c.GetSamplePoints(interval.LowerBound, interval.UpperBound, 0.02);
Point2d[] pts = samplePts.Select(p => transform * p.Point).ToArray();
for (int n = 0; n < pts.Length - 1; n++)
{
LineSegment2d lineSeg = new LineSegment2d(pts[n], pts[n + 1]);
compositeCurveElems.Add(lineSeg);
ptsColl.Add(pts[n]);
ptParams.Add(currParam);
updateExtends(pts[n]);
currParam += lineSeg.Length;
}
}
else
{
LineSegment2d lineSeg = (LineSegment2d)c;
lineSeg.TransformBy(transform);
compositeCurveElems.Add(lineSeg);
ptsColl.Add(lineSeg.StartPoint);
ptParams.Add(currParam);
updateExtends(lineSeg.StartPoint);
currParam += lineSeg.Length;
}
}
CompositeCurve2d boundary = new CompositeCurve2d(compositeCurveElems.ToArray());
Extents2d ext = new Extents2d(_minx, _miny, _maxx, _maxy);
boundaries.Add(boundary);
ptsCollList.Add(ptsColl);
ptParamsList.Add(ptParams);
extends.Add(ext);
}
}
//контуры штриховок не могут иметь самопересечений!
#region Проверка на пересечения
//проверка на самопересечения
//bool badBoundaries = false;
HashSet <int> badBoundaries = new HashSet <int>();
HashSet <int> splitBoundaries = new HashSet <int>();//Если 2 контура в одной штриховке пересекаются, то разносить их по разным штриховкам
//List<HatchData> decomposeHatchData = new List<HatchData>();//TODO: самопересекающиеся полигоны нужно разбить на отдельные по количеству самопересечний.
for (int i = 0; i < boundaries.Count; i++)
{
CompositeCurve2d b = boundaries[i];
CurveCurveIntersector2d intersector = new CurveCurveIntersector2d(b, b);
if (intersector.NumberOfIntersectionPoints > 0)
{
//если происходит только наложение???
badBoundaries.Add(i);
selfintersectingLoops.Add(ptsCollList[i]);
}
}
if (boundaries.Count > 1)
{
//проверка на взаимные пересечения.
//Исп RBush для того чтобы избежать проверки на пересечение каждого с каждым и квадратичной сложности
//(работает только если контуры разнесены)
//Не брать в расчет пересечения по касательной
RBush <Spatial> boundariesRBush = new RBush <Spatial>();
List <Spatial> spatialData = new List <Spatial>();
for (int i = 0; i < extends.Count; i++)
{
spatialData.Add(new Spatial(extends[i], i));
}
boundariesRBush.BulkLoad(spatialData);
foreach (Spatial s in spatialData)
{
IReadOnlyList <Spatial> nearestNeighbors = boundariesRBush.Search(s.Envelope);
if (nearestNeighbors.Count > 1)
{
CompositeCurve2d thisCurve = boundaries[(int)s.Obj];
foreach (Spatial n in nearestNeighbors)
{
if (!s.Equals(n))
{
CompositeCurve2d otherCurve = boundaries[(int)n.Obj];
CurveCurveIntersector2d intersector
= new CurveCurveIntersector2d(thisCurve, otherCurve);
if (intersector.NumberOfIntersectionPoints > 0 ||
intersector.OverlapCount > 0)
{
bool matches = false;
//Проверить, что кривые не накладываются друг на друга по всей длине (то есть полностью совпадают)
if (intersector.OverlapCount > 0)
{
//сумма длин всех интервалов перекрытия равна общей длине кривой
double thisCurveOverlapLength = 0;
double otherCurveOverlapLength = 0;
for (int i = 0; i < intersector.OverlapCount; i++)
{
Interval[] intervals = intersector.GetOverlapRanges(i);
Interval thisOverlapInterval = intervals[0];
thisCurveOverlapLength += thisOverlapInterval.Length;
Interval otherOverlapInterval = intervals[1];
otherCurveOverlapLength += otherOverlapInterval.Length;
}
Interval thisCurveInterval = thisCurve.GetInterval();
Interval otherCurveInterval = otherCurve.GetInterval();
if (Utils.LengthIsEquals(thisCurveOverlapLength, thisCurveInterval.Length) &&
Utils.LengthIsEquals(otherCurveOverlapLength, otherCurveInterval.Length))
{
matches = true;
}
}
if (!matches)
{
splitBoundaries.Add((int)s.Obj);
splitBoundaries.Add((int)n.Obj);
}
else
{
badBoundaries.Add((int)s.Obj);
badBoundaries.Add((int)n.Obj);
}
}
}
}
}
}
}
splitBoundaries.ExceptWith(badBoundaries);
List <HatchData> splitHatchData = new List <HatchData>();
if (badBoundaries.Count > 0 || splitBoundaries.Count > 0)
{
List <CompositeCurve2d> boundariesClear = new List <CompositeCurve2d>();
List <Extents2d> extendsClear = new List <Extents2d>();
List <Point2dCollection> ptsCollListClear = new List <Point2dCollection>();
List <List <double> > ptParamsListClear = new List <List <double> >();
for (int i = 0; i < boundaries.Count; i++)
{
if (!badBoundaries.Contains(i) && !splitBoundaries.Contains(i))
{
boundariesClear.Add(boundaries[i]);
extendsClear.Add(extends[i]);
ptsCollListClear.Add(ptsCollList[i]);
ptParamsListClear.Add(ptParamsList[i]);
}
}
foreach (int index in splitBoundaries)
{
splitHatchData.Add(new HatchData(
new HatchNestingNode(
boundaries[index],
extends[index],
ptsCollList[index],
ptParamsList[index], hatch)));
}
boundaries = boundariesClear;
extends = extendsClear;
ptsCollList = ptsCollListClear;
ptParamsList = ptParamsListClear;
}
#endregion
//определяется вложенность контуров штриховки
//ЕСЛИ ШТРИХОВКА СОСТОИТ ИЗ 2 И БОЛЕЕ КОНТУРОВ, КОТОРЫЕ НЕ ВЛОЖЕНЫ ДРУГ В ДРУГА,
//ТО ЭТИ КОНТУРЫ ДОЛЖНЫ РАССМАТРИВАТЬСЯ КАК ОТДЕЛЬНЫЕ ШТРИХОВКИ!!!
HatchNestingTree hatchNestingTree
= new HatchNestingTree(boundaries, extends, ptsCollList, ptParamsList, hatch);
List <HatchData> currHatchData = hatchNestingTree.GetHatchData();
currHatchData.AddRange(splitHatchData);//добавить контуры, полученные из взаимно пересекающихся контуров
//currHatchData.AddRange(decomposeHatchData);
allHatchData.AddRange(currHatchData);
//Каждая штриховка имеет диапазон по X от начала до конца по оси.
//В общую очередь событий сохраняются события начала и конца штриховки
foreach (HatchData hd in currHatchData)
{
hd.AddEventsToQueue(eventQueue);
}
}
//Трассу разбить на отрезки, на которых будут вставлены прямолинейные сегменты штриховок
Polyline alignmentPoly = null;
try
{
alignmentPoly = (Polyline)tr.GetObject(AlignmentPolyId, OpenMode.ForRead);
}
catch (Autodesk.AutoCAD.Runtime.Exception)
{
return;
}
int segments = alignmentPoly.NumberOfVertices - 1;
List <AlignmentSegment> alignmentSegments = new List <AlignmentSegment>();
Action <Point2d, Point2d> addAlignmentSegment = new Action <Point2d, Point2d>((p0, p1) =>
{
double start = alignmentPoly.GetDistAtPoint(new Point3d(p0.X, p0.Y, 0));
double end = alignmentPoly.GetDistAtPoint(new Point3d(p1.X, p1.Y, 0));
if (Math.Abs(start - end) > Tolerance.Global.EqualPoint)//TODO: Это спорный момент - ведь может быть большое множество очень коротких участков подряд!
{
Vector2d startDir = p1 - p0;
alignmentSegments.Add(new AlignmentSegment(start, end, p0, startDir));
}
});
for (int i = 0; i < segments; i++)
{
SegmentType segmentType = alignmentPoly.GetSegmentType(i);
Point2d startLoc = alignmentPoly.GetPoint2dAt(i);
Point2d endLoc = alignmentPoly.GetPoint2dAt(i + 1);
switch (segmentType)
{
case SegmentType.Line:
addAlignmentSegment(startLoc, endLoc);
break;
case SegmentType.Arc:
CircularArc2d arc = new CircularArc2d(startLoc, endLoc, alignmentPoly.GetBulgeAt(i), false);
Interval interval = arc.GetInterval();
PointOnCurve2d[] samplePts = arc.GetSamplePoints(interval.LowerBound, interval.UpperBound, 0.1);
for (int n = 0; n < samplePts.Length - 1; n++)
{
addAlignmentSegment(samplePts[n].Point, samplePts[n + 1].Point);
}
break;
}
}
//проход по каждому отрезку трассы (диапазон отрезка - диапазон длины полилинии)
HashSet <HatchData> currentHatchData = new HashSet <HatchData>();
foreach (AlignmentSegment alignmentSegment in alignmentSegments)
{
if (eventQueue.Count == 0)
{
break; //штриховки закончились
}
//Получить те штриховки, диапазон по X которых пересекается с диапазоном текущего отрезка
//(СКАНИРУЮЩАЯ ЛИНИЯ: события - начало, конец штриховки)
HashSet <HatchData> intervalHatchData = new HashSet <HatchData>(currentHatchData);//штриховки пришедшие из предыдущего участка остаются все
//Собрать все события до конца сегмента
//Если при проходе от начала до конца сегмента какая-то штриховка проходится полностью от начала до конца, то
//все ее контуры без изменений должны быть переданы для создания М-полигона!!! (для них обход графа не нужен!)
HashSet <HatchData> startedInsideInterval = new HashSet <HatchData>();
List <HatchData> hatchesCompletelyInsideInterval = new List <HatchData>();
while (eventQueue.Count > 0)
{
HatchEvent nextEvent = eventQueue.FindMin();
if (nextEvent.Position > alignmentSegment.End)
{
break;
}
else if (nextEvent.Start)
{
//добавить штриховку в текущий набор
HatchData hd = eventQueue.DeleteMin().HatchData;
currentHatchData.Add(hd);
//добавлять в набор текущего интервла только в том случае,
//если сканирующая линия еще не дошла до конца интервала
if (nextEvent.Position < alignmentSegment.End &&
!Utils.LengthIsEquals(nextEvent.Position, alignmentSegment.End)) //Допуск нужен
{
startedInsideInterval.Add(hd);
intervalHatchData.Add(hd);
}
}
else
{
//убрать штриховку из текущего набора
HatchData hd = eventQueue.DeleteMin().HatchData;
currentHatchData.Remove(hd);
if (startedInsideInterval.Contains(hd))
{
hatchesCompletelyInsideInterval.Add(hd);
}
}
}
foreach (HatchData hd in hatchesCompletelyInsideInterval)
{
HatchSegmentData hsd = new HatchSegmentData(hd);
alignmentSegment.HatchSegmentData.Add(hsd);
hsd.Polygons = hd.GetAllBoundaries();
}
intervalHatchData.ExceptWith(hatchesCompletelyInsideInterval);
foreach (HatchData hd in intervalHatchData)
{
HatchSegmentData hsd = new HatchSegmentData(hd);
alignmentSegment.HatchSegmentData.Add(hsd);
//для каждой штриховки выполнить построение и обход графа сегмента штриховки
HatchSegmentGraph graph = new HatchSegmentGraph(alignmentSegment.Start, alignmentSegment.End, hd, doc.Editor);
//сохранить наборы полигонов для текущего диапазона
hsd.Polygons = graph.Result;
}
}
//для каждого диапазона создать полученные полигоны в 3d
BlockTable bt = tr.GetObject(db.BlockTableId, OpenMode.ForWrite) as BlockTable;
BlockTableRecord ms
= (BlockTableRecord)tr.GetObject(bt[BlockTableRecord.ModelSpace], OpenMode.ForWrite);
BlockTableRecord btr = new BlockTableRecord();//Каждый профиль в отдельный блок
btr.Name = Guid.NewGuid().ToString();
ObjectId btrId = bt.Add(btr);
tr.AddNewlyCreatedDBObject(btr, true);
foreach (AlignmentSegment alignmentSegment in alignmentSegments)
{
List <Entity> flatObjs = new List <Entity>();
foreach (HatchSegmentData hsd in alignmentSegment.HatchSegmentData)
{
//PlaneSurface
hsd.GetNesting();
Dictionary <Point2dCollection, List <Point2dCollection> > bwhs = hsd.GetBoundariesWithHoles();
foreach (KeyValuePair <Point2dCollection, List <Point2dCollection> > bwh in bwhs)
{
//создать Region
Region region = null;
using (Polyline poly = new Polyline())
{
for (int i = 0; i < bwh.Key.Count; i++)
{
poly.AddVertexAt(i, bwh.Key[i], 0, 0, 0);
}
poly.Closed = true;
DBObjectCollection coll = new DBObjectCollection();
coll.Add(poly);
try
{
DBObjectCollection regionColl = Region.CreateFromCurves(coll);
foreach (DBObject dbo in regionColl)
{
region = (Region)dbo;
break;
}
}
catch { }
}
//из Region создать PlaneSurface
if (region != null)
{
using (PlaneSurface planeSurface = new PlaneSurface())
{
planeSurface.CreateFromRegion(region);
planeSurface.LayerId = hsd.Hatch.LayerId;
planeSurface.ColorIndex = 256;
ObjectId planeSurfaceId = ms.AppendEntity(planeSurface);
tr.AddNewlyCreatedDBObject(planeSurface, true);
//вырезать отверстия в PlaneSurface
foreach (Point2dCollection holePts2d in bwh.Value)
{
if (holePts2d.Count < 3)
{
continue;
}
using (Polyline poly = new Polyline())
{
for (int i = 0; i < holePts2d.Count; i++)
{
poly.AddVertexAt(i, holePts2d[i], 0, 0, 0);
}
poly.Closed = true;
ObjectIdCollection trimPolyColl = new ObjectIdCollection();
trimPolyColl.Add(ms.AppendEntity(poly));
tr.AddNewlyCreatedDBObject(poly, true);
List <Point2d> ptsList = new List <Point2d>(holePts2d.ToArray());
Point2d pickPt2d = Utils.GetAnyPointInsidePoligon(ptsList,
Utils.DirectionIsClockwise(ptsList));
try
{
AcadDB.Surface.TrimSurface(planeSurfaceId, new ObjectIdCollection(), trimPolyColl,
new Vector3dCollection()
{
Vector3d.ZAxis
}, new Point3d(pickPt2d.X, pickPt2d.Y, 0),
-Vector3d.ZAxis, false, false);
}
catch /*(Exception ex)*/
{
//Вывод в командную строку
Utils.ErrorToCommandLine(doc.Editor,
"Ошибка при попытке вырезания отверстия в поверхности" /*, ex*/);
}
//Удалить все объекты, добавленные в чертеж!
poly.Erase();
}
}
flatObjs.Add((Entity)planeSurface.Clone());
//Удалить все объекты, добавленные в чертеж!
planeSurface.Erase();
}
region.Dispose();
}
}
}
foreach (Entity ent in flatObjs)
{
ent.TransformBy(alignmentSegment.Transform);
/*ms*/
btr.AppendEntity(ent);
tr.AddNewlyCreatedDBObject(ent, true);
ent.Dispose();
}
}
BlockReference br = new BlockReference(Point3d.Origin, btrId);
ms.AppendEntity(br);
tr.AddNewlyCreatedDBObject(br, true);
if (selfintersectingLoops.Count > 0)
{
Utils.ErrorToCommandLine(doc.Editor,
"Отбраковано самопересекающихся контуров штриховок - " + selfintersectingLoops.Count + " (отмечены на профиле)");
ObjectId layerId = Utils.CreateLayerIfNotExists("САМОПЕРЕСЕКАЮЩИЙСЯ КОНТУР ШТРИХОВКИ", db, tr,
color: Color.FromColorIndex(ColorMethod.ByAci, 1), lineWeight: LineWeight.LineWeight200);
Matrix2d returnTransform = transform.Inverse();
foreach (Point2dCollection pts in selfintersectingLoops)
{
using (Polyline selfIntersectingPoly = new Polyline())
{
selfIntersectingPoly.LayerId = layerId;
selfIntersectingPoly.ColorIndex = 256;
for (int i = 0; i < pts.Count; i++)
{
Point2d pt = pts[i].TransformBy(returnTransform);
selfIntersectingPoly.AddVertexAt(i, pt, 0, 0, 0);
}
toHighlight.Add(ms.AppendEntity(selfIntersectingPoly));
tr.AddNewlyCreatedDBObject(selfIntersectingPoly, true);
//selfIntersectingPoly.Highlight();
}
}
}
tr.Commit();
}
}
ps.Visible = false;
if (toHighlight.Count > 0)
{
using (Transaction tr = doc.Database.TransactionManager.StartTransaction())
{
foreach (ObjectId id in toHighlight)
{
Entity ent = (Entity)tr.GetObject(id, OpenMode.ForRead);
ent.Highlight();
}
tr.Commit();
}
}
}
}
catch (System.Exception ex)
{
GeologyConvertationCommand.ClosePalette(null, null);
CommonException(ex, "Ошибка при создании 3d профиля геологии");
}
}
