public XYZ ProjectLineToIntersection(Line line, StreetIntersection intersection)
{
if (intersection == null)
{
return(null);
}
foreach (var curve in intersection.CurveLoop)
{
curve.Intersect(line, out IntersectionResultArray resultArray);
if (resultArray != null && !resultArray.IsEmpty)
{
var result = resultArray.get_Item(0);
return(result.XYZPoint);
}
}
return(null);
}
private CurveLoop CreateIntersectionCurveLoop(List <StreetSegment> intersStreetSegments, StreetIntersection intersection)
{
var arc = Arc.Create(intersection.Origin, osmStore.DefaultStreetWidth, 0, Math.PI * 1.99, XYZ.BasisX, XYZ.BasisY);
var lineOffsetDic = new Dictionary <Line, double>();
var offsetLines = new List <Line>();
foreach (var segm in intersStreetSegments)
{
var line0 = segm.Line.CreateOffset(osmStore.DefaultStreetWidth / 2, XYZ.BasisZ) as Line;
offsetLines.Add(line0);
var line1 = segm.Line.CreateOffset(osmStore.DefaultStreetWidth / 2, -XYZ.BasisZ) as Line;
offsetLines.Add(line1);
}
foreach (var line in offsetLines)
{
arc.Intersect(line, out IntersectionResultArray resultArray);
if (resultArray != null && !resultArray.IsEmpty)
{
var result = resultArray.get_Item(0);
var intPoint = result.XYZPoint;
var pResult = arc.Project(intPoint);
var param = arc.ComputeNormalizedParameter(pResult.Parameter);
lineOffsetDic.Add(line, param);
}
}
var orderedDic = lineOffsetDic.OrderBy(d => d.Value).ToList();
if (orderedDic.Count == 0)
{
return(null);
}
orderedDic.Add(orderedDic[0]);
var loop = new CurveLoop();
for (int i = 0; i < orderedDic.Count - 1; i++)
{
var start = arc.Evaluate(orderedDic[i].Value, true);
var end = arc.Evaluate(orderedDic[i + 1].Value, true);
if (start.DistanceTo(end) < 0.05)
{
return(null);
}
var line = Line.CreateBound(start, end);
loop.Append(line);
}
return(loop);
}
public List <StreetIntersection> CreateIntersections(List <StreetSegment> streetSegments)
{
var uniquePoints = new List <XYZ>();
foreach (var line in streetSegments.Select(s => s.Line))
{
var start = line.GetEndPoint(0);
var end = line.GetEndPoint(1);
if (!geometryService.ListContainsPoint(uniquePoints, start))
{
uniquePoints.Add(start);
}
if (!geometryService.ListContainsPoint(uniquePoints, start))
{
uniquePoints.Add(end);
}
}
var intersections = new List <StreetIntersection>();
foreach (var point in uniquePoints)
{
var count = 0;
var intersStreetSegments = new List <StreetSegment>();
foreach (var segment in streetSegments)
{
if (segment.Line.Distance(point) < 0.0001)
{
count++;
intersStreetSegments.Add(segment);
}
}
if (count > 1)
{
var intersection = new StreetIntersection();
intersection.Id = osmStore.MoveNextId();
intersection.StreetIds = intersStreetSegments.Select(s => s.Id).Distinct().ToList();
intersection.NumberOfStreets = intersection.StreetIds.Count;
intersection.StreetSegmentIds = intersStreetSegments.Select(s => s.SegmentId).Distinct().ToList();
intersection.NumberOfStreetSegments = intersection.StreetSegmentIds.Count();
intersection.Origin = point;
var streetNames = intersStreetSegments.Select(s => s.Name).Where(n => !string.IsNullOrWhiteSpace(n)).Distinct();
var nameCount = streetNames.Count();
if (nameCount == 1)
{
intersection.Name = streetNames.ElementAt(0);
}
else if (nameCount >= 2)
{
intersection.Name = $"{streetNames.ElementAt(0)} with {streetNames.ElementAt(1)}";
for (int i = 2; i < nameCount; i++)
{
intersection.Name += $"and {streetNames.ElementAt(i)}";
}
}
intersection.CurveLoop = CreateIntersectionCurveLoop(intersStreetSegments, intersection);
if (intersection.CurveLoop != null)
{
intersections.Add(intersection);
}
}
}
return(intersections);
}