public int Compare(object x, object y)
{ // Typ-Casting
IGeoObject o1 = (IGeoObject)x;
IGeoObject o2 = (IGeoObject)y;
// Sortierung nach Mittelpunkten
BoundingRect re = IGeoObjectImpl.GetExtent(o1, projection, false);
GeoPoint2D p1 = re.GetCenter();
re = IGeoObjectImpl.GetExtent(o2, projection, false);
GeoPoint2D p2 = re.GetCenter();
if (dirx)
{
if (p1.x < p2.x)
{
return(-1);
}
if (p1.x > p2.x)
{
return(1);
}
}
else
{
if (p1.y < p2.y)
{
return(-1);
}
if (p1.y > p2.y)
{
return(1);
}
}
return(0);
}
static public void SameHeight(GeoObjectList gl, Projection projection, Project pr)
{
using (pr.Undo.UndoFrame)
{
BoundingRect re1 = IGeoObjectImpl.GetExtent(gl[gl.Count - 1], projection, false);
for (int i = 0; i < gl.Count - 1; ++i)
{
BoundingRect re = IGeoObjectImpl.GetExtent(gl[i], projection, false);
double hig = 1;
if (re.Height != 0)
{
hig = re1.Height / re.Height;
}
GeoPoint refPkt = projection.DrawingPlane.ToGlobal(re.GetCenter());
ModOp m = ModOp.Scale(refPkt, new GeoVector(0.0, 1.0, 0.0), hig);
gl[i].Modify(m);
}
}
}
private int debugCount; // to identify instance when debugging
#endif
public ProjectedCurve(ICurve curve3D, ISurface surface, bool forward, BoundingRect domain, double precision = 0.0)
{
#if DEBUG
debugCount = debugCounter++;
#endif
this.curve3D = curve3D; // keep in mind, the curve is not cloned, curve3D should not be modified after this
this.surface = surface;
List <GeoPoint> lpoles = new List <GeoPoint>();
List <GeoPoint2D> lpoles2d = new List <GeoPoint2D>();
GeoPoint2D cnt2d = domain.GetCenter();
GeoPoint sp = curve3D.StartPoint;
GeoPoint ep = curve3D.EndPoint;
double[] us = surface.GetUSingularities();
double prec = precision;
if (prec == 0.0)
{
prec = curve3D.Length * 1e-3; // changed to 1e-3, it is used to snap endpoints to poles
}
startPoint2d = surface.PositionOf(curve3D.StartPoint);
endPoint2d = surface.PositionOf(curve3D.EndPoint);
bool distinctStartEndPoint = false;
if ((surface.IsUPeriodic && Math.Abs(startPoint2d.x - endPoint2d.x) < surface.UPeriod * 1e-3) ||
(surface.IsVPeriodic && Math.Abs(startPoint2d.y - endPoint2d.y) < surface.VPeriod * 1e-3))
{ // adjust start and endpoint according to its neighbors
GeoPoint2D p2d = surface.PositionOf(curve3D.PointAt(0.1));
SurfaceHelper.AdjustPeriodic(surface, periodicDomain, ref p2d);
BoundingRect ext = new BoundingRect(p2d);
SurfaceHelper.AdjustPeriodic(surface, ext, ref startPoint2d);
p2d = surface.PositionOf(curve3D.PointAt(0.9));
SurfaceHelper.AdjustPeriodic(surface, periodicDomain, ref p2d);
ext = new BoundingRect(p2d);
SurfaceHelper.AdjustPeriodic(surface, ext, ref endPoint2d);
distinctStartEndPoint = true;
}
periodicDomain = domain;
if (periodicDomain.IsEmpty() && (surface.IsUPeriodic || surface.IsVPeriodic))
{
// make a few points and assure that they don't jump over the periodic seam
// if the curve3d doesn't jump around wildly, this should work. Maybe use curve3D.GetSavePositions?
GeoPoint2D[] point2Ds = new GeoPoint2D[11];
for (int i = 0; i < 11; i++)
{
point2Ds[i] = surface.PositionOf(curve3D.PointAt(i / 10.0));
}
for (int i = 0; i < 10; i++)
{
GeoVector2D offset = GeoVector2D.NullVector;
if (surface.IsUPeriodic && Math.Abs(point2Ds[i + 1].x - point2Ds[i].x) > surface.UPeriod / 2.0)
{
if ((point2Ds[i + 1].x - point2Ds[i].x) < 0)
{
offset.x = surface.UPeriod;
}
else
{
offset.x = -surface.UPeriod;
}
}
if (surface.IsVPeriodic && Math.Abs(point2Ds[i + 1].y - point2Ds[i].y) > surface.VPeriod / 2.0)
{
if ((point2Ds[i + 1].y - point2Ds[i].y) < 0)
{
offset.y = surface.VPeriod;
}
else
{
offset.y = -surface.VPeriod;
}
}
point2Ds[i + 1] += offset;
}
for (int i = 0; i < 11; i++)
{
periodicDomain.MinMax(point2Ds[i]);
}
startPoint2d = point2Ds[0];
endPoint2d = point2Ds[10];
}
if (!periodicDomain.IsEmpty() && (!surface.IsUPeriodic || periodicDomain.Width < surface.UPeriod * (1 - 1e-6)) && (!surface.IsVPeriodic || periodicDomain.Height < surface.VPeriod * (1 - 1e-6)))
{
SurfaceHelper.AdjustPeriodic(surface, periodicDomain, ref startPoint2d);
SurfaceHelper.AdjustPeriodic(surface, periodicDomain, ref endPoint2d);
}
startPointIsPole = endPointIsPole = false;
for (int i = 0; i < us.Length; i++)
{
GeoPoint pl = surface.PointAt(new GeoPoint2D(us[i], cnt2d.y));
if ((pl | sp) < prec)
{
GeoPoint2D tmp = surface.PositionOf(curve3D.PointAt(0.1));
startPoint2d = new GeoPoint2D(us[i], tmp.y);
startPointIsPole = true;
spu = true;
}
if ((pl | ep) < prec)
{
GeoPoint2D tmp = surface.PositionOf(curve3D.PointAt(0.9));
endPoint2d = new GeoPoint2D(us[i], tmp.y);
endPointIsPole = true;
epu = true;
}
}
double[] vs = surface.GetVSingularities();
for (int i = 0; i < vs.Length; i++)
{
GeoPoint pl = surface.PointAt(new GeoPoint2D(cnt2d.x, vs[i]));
if ((pl | sp) < prec)
{
GeoPoint2D tmp = surface.PositionOf(curve3D.PointAt(0.1));
startPoint2d = new GeoPoint2D(tmp.x, vs[i]);
startPointIsPole = true;
spu = false;
}
if ((pl | ep) < prec)
{
GeoPoint2D tmp = surface.PositionOf(curve3D.PointAt(0.9));
endPoint2d = new GeoPoint2D(tmp.x, vs[i]);
endPointIsPole = true;
epu = false;
}
}
if (forward)
{
startParam = 0.0;
endParam = 1.0;
}
else
{
startParam = 1.0;
endParam = 0.0;
}
#if DEBUG
this.MakeTriangulation();
#endif
}
internal void CenterPatch(LayoutPatch patch, double scale, HorizontalCenter hor, VerticalCenter ver)
{
// die Projektion ist ja zwei Anteile, unscaledProjection hält immer den Nullpunkt
// fest und skaliert nicht, und Placement platziert
BoundingRect areaext;
if (patch.Area != null)
{
areaext = patch.Area.Extent;
}
else
{
areaext = new BoundingRect(0.0, 0.0, paperWidth, paperHeight);
}
BoundingRect modelext = patch.Model.GetExtent(patch.Projection);
if (modelext.IsEmpty())
{
return;
}
GeoPoint2D modelcnt = modelext.GetCenter();
GeoPoint2D areacnt = areaext.GetCenter();
double factor, dx, dy;
patch.Projection.GetPlacement(out factor, out dx, out dy);
if (scale != 0.0)
{
factor = scale;
}
switch (hor)
{
case HorizontalCenter.left:
dx = areaext.Left - (modelext.Left * factor);
break;
case HorizontalCenter.center:
dx = areacnt.x - (modelcnt.x * factor);
break;
case HorizontalCenter.right:
dx = areaext.Right - (modelext.Right * factor);
break;
default:
break;
}
switch (ver)
{
case VerticalCenter.bottom:
dy = areaext.Bottom - (modelext.Bottom * factor);
break;
case VerticalCenter.center:
dy = areacnt.y - (modelcnt.y * factor);
break;
case VerticalCenter.top:
dy = areaext.Top - (modelext.Top * factor);
break;
default:
break;
}
patch.Projection.SetPlacement(factor, dx, dy);
}