public VoxelGrid(RhCommon_Scene ModelSurfaces, List <IOn3dPoint> Pts, int VG_Domain)
{
VoxelCT = VG_Domain;
VoxelInventory = new List <int> [VoxelCT, VoxelCT, VoxelCT];
Voxel = new List <On3dPoint> [VoxelCT, VoxelCT, VoxelCT];
IRhinoObject[] Surfaces = new IRhinoObject[ModelSurfaces.Count()];
OnBoundingBox TightOBox = new OnBoundingBox();
On3dPoint BoxMin = new On3dPoint();
On3dPoint BoxMax = new On3dPoint();
ArrayOn3dPoint Points = new ArrayOn3dPoint();
for (int i = 0; i < Pts.Count; i++)
{
Points.Append(Pts[i]);
}
RhUtil.RhinoGetTightBoundingBox(Surfaces, ref TightOBox, false, null);
TightOBox.Set(Points, true);
OverallBBox = new OnBoundingBox(new On3dPoint(TightOBox.Min().x - 1, TightOBox.Min().y - 1, TightOBox.Min().z - 1), new On3dPoint(TightOBox.Max().x + 1, TightOBox.Max().y + 1, TightOBox.Max().z + 1));
this.X_Incr = (OverallBBox.Max().x - OverallBBox.Min().x) / VoxelCT;
this.Y_Incr = (OverallBBox.Max().y - OverallBBox.Min().y) / VoxelCT;
this.Z_Incr = (OverallBBox.Max().z - OverallBBox.Min().z) / VoxelCT;
//For((int XBox = 0; XBox < VoxelCT; XBox++)
Parallel.For(0, VoxelCT, XBox =>
{
RhUtil.RhinoApp().SetStatusBarMessagePane(string.Format("Voxelizing: {0}%", Math.Round((double)XBox / VoxelCT - 1, 2) * 100));
for (int YBox = 0; YBox < VoxelCT; YBox++)
{
for (int ZBox = 0; ZBox < VoxelCT; ZBox++)
{
BoxMin = new On3dPoint((OverallBBox.Min().x + this.X_Incr * XBox) - X_Incr / 10, (OverallBBox.Min().y + this.Y_Incr * YBox) - Y_Incr / 10, (OverallBBox.Min().z + this.Z_Incr * ZBox) - Z_Incr / 10);
BoxMax = new On3dPoint((OverallBBox.Min().x + this.X_Incr * (XBox + 1)) + X_Incr / 10, (OverallBBox.Min().y + this.Y_Incr * (YBox + 1)) + Y_Incr / 10, (OverallBBox.Min().z + this.Z_Incr * (ZBox + 1)) + X_Incr / 10);
this.Voxel[XBox, YBox, ZBox] = new List <On3dPoint>();
this.Voxel[XBox, YBox, ZBox].Add(BoxMin);
this.Voxel[XBox, YBox, ZBox].Add(BoxMax);
this.VoxelInventory[XBox, YBox, ZBox] = new List <int>();
for (int Index = 0; Index < ModelSurfaces.Count(); Index++)
{
OnBoundingBox TestBox = new OnBoundingBox();
TestBox = new OnBoundingBox(this.Voxel[XBox, YBox, ZBox][0], this.Voxel[XBox, YBox, ZBox][1]);
if (BoxIntersection(ModelSurfaces, TestBox, Index))
{
this.VoxelInventory[XBox, YBox, ZBox].Add(Index);
}
}
}
}
});
}
public VoxelGrid(RhCommon_Scene ModelSurfaces, List<IOn3dPoint> Pts, int VG_Domain)
{
VoxelCT = VG_Domain;
VoxelInventory = new List<int>[VoxelCT, VoxelCT, VoxelCT];
Voxel = new List<On3dPoint>[VoxelCT, VoxelCT, VoxelCT];
IRhinoObject[] Surfaces = new IRhinoObject[ModelSurfaces.Count()];
OnBoundingBox TightOBox = new OnBoundingBox();
On3dPoint BoxMin = new On3dPoint();
On3dPoint BoxMax = new On3dPoint();
ArrayOn3dPoint Points = new ArrayOn3dPoint();
for (int i = 0; i < Pts.Count; i++)
{
Points.Append(Pts[i]);
}
RhUtil.RhinoGetTightBoundingBox(Surfaces, ref TightOBox, false, null);
TightOBox.Set(Points, true);
OverallBBox = new OnBoundingBox(new On3dPoint(TightOBox.Min().x - 1, TightOBox.Min().y - 1, TightOBox.Min().z - 1), new On3dPoint(TightOBox.Max().x + 1, TightOBox.Max().y + 1, TightOBox.Max().z + 1));
this.X_Incr = (OverallBBox.Max().x - OverallBBox.Min().x) / VoxelCT;
this.Y_Incr = (OverallBBox.Max().y - OverallBBox.Min().y) / VoxelCT;
this.Z_Incr = (OverallBBox.Max().z - OverallBBox.Min().z) / VoxelCT;
//For((int XBox = 0; XBox < VoxelCT; XBox++)
Parallel.For(0, VoxelCT, XBox =>
{
RhUtil.RhinoApp().SetStatusBarMessagePane(string.Format("Voxelizing: {0}%", Math.Round((double)XBox / VoxelCT-1, 2) * 100));
for (int YBox = 0; YBox < VoxelCT; YBox++)
{
for (int ZBox = 0; ZBox < VoxelCT; ZBox++)
{
BoxMin = new On3dPoint((OverallBBox.Min().x + this.X_Incr * XBox) - X_Incr / 10, (OverallBBox.Min().y + this.Y_Incr * YBox) - Y_Incr / 10, (OverallBBox.Min().z + this.Z_Incr * ZBox) - Z_Incr / 10);
BoxMax = new On3dPoint((OverallBBox.Min().x + this.X_Incr * (XBox + 1)) + X_Incr / 10, (OverallBBox.Min().y + this.Y_Incr * (YBox + 1)) + Y_Incr / 10, (OverallBBox.Min().z + this.Z_Incr * (ZBox + 1)) + X_Incr / 10);
this.Voxel[XBox, YBox, ZBox] = new List<On3dPoint>();
this.Voxel[XBox, YBox, ZBox].Add(BoxMin);
this.Voxel[XBox, YBox, ZBox].Add(BoxMax);
this.VoxelInventory[XBox, YBox, ZBox] = new List<int>();
for (int Index = 0; Index < ModelSurfaces.Count(); Index++)
{
OnBoundingBox TestBox = new OnBoundingBox();
TestBox = new OnBoundingBox(this.Voxel[XBox, YBox, ZBox][0], this.Voxel[XBox, YBox, ZBox][1]);
if (BoxIntersection(ModelSurfaces, TestBox, Index))
{
this.VoxelInventory[XBox, YBox, ZBox].Add(Index);
}
}
}
}
});
}
public void SetBrep(IOnBrep brep)
{
if (null == brep || !brep.IsValid())
{
return;
}
m_brep = brep;
int face_count = m_brep.m_F.Count();
m_points.Reserve(face_count * SURFACE_ARROW_COUNT * SURFACE_ARROW_COUNT);
m_normals.Reserve(face_count * SURFACE_ARROW_COUNT * SURFACE_ARROW_COUNT);
m_points.SetCount(0);
m_normals.SetCount(0);
for (int i = 0; i < face_count; i++)
{
IOnBrepFace face = m_brep.m_F[i];
IOnBrepLoop loop = face.OuterLoop();
if (null == loop)
{
continue;
}
OnInterval udomain = face.Domain(0);
OnInterval vdomain = face.Domain(1);
if (loop.m_pbox.IsValid())
{
OnInterval domain = new OnInterval();
domain.Set(loop.m_pbox.m_min.x, loop.m_pbox.m_max.x);
domain.Intersection(udomain);
if (domain.IsIncreasing())
{
udomain.Set(domain.Min(), domain.Max());
}
domain.Set(loop.m_pbox.m_min.y, loop.m_pbox.m_max.y);
domain.Intersection(vdomain);
if (domain.IsIncreasing())
{
vdomain.Set(domain.Min(), domain.Max());
}
}
bool bUntrimmed = m_brep.FaceIsSurface(i);
ArrayOnInterval intervals = new ArrayOnInterval();
bool bRev = face.m_bRev;
for (double u = 0.0; u < SURFACE_ARROW_COUNT; u += 1.0)
{
double d = u / (SURFACE_ARROW_COUNT - 1.0);
double s = udomain.ParameterAt(d);
intervals.SetCount(0);
if (bUntrimmed || RhUtil.RhinoGetIsoIntervals(face, 1, s, intervals) > 0)
{
for (double v = 0.0; v < SURFACE_ARROW_COUNT; v += 1.0)
{
d = v / (SURFACE_ARROW_COUNT - 1.0);
double t = vdomain.ParameterAt(d);
bool bAdd = bUntrimmed;
for (int k = 0; !bAdd && k < intervals.Count(); k++)
{
if (intervals[k].Includes(t))
{
bAdd = true;
}
}
if (bAdd)
{
On3dPoint pt = new On3dPoint();
On3dVector du = new On3dVector();
On3dVector dv = new On3dVector();
On3dVector dir = new On3dVector();
if (face.EvNormal(s, t, ref pt, ref du, ref dv, ref dir))
{
m_points.Append(pt);
if (bRev)
{
dir.Reverse();
}
m_normals.Append(dir);
}
}
}
}
}
}
}
