static public void CreateTunnel()
{
Document doc = Application.DocumentManager.MdiActiveDocument;
Editor ed = doc.Editor;
Database db = doc.Database;
//===========
Matrix3d curUCSMatrix = doc.Editor.CurrentUserCoordinateSystem;
CoordinateSystem3d curUCS = curUCSMatrix.CoordinateSystem3d;
double TunnelDia = 6.04;
PromptEntityResult per = ed.GetEntity("Select polylines");
ObjectId oid = per.ObjectId;
if (per.Status == PromptStatus.OK)
{
using (Transaction tr = db.TransactionManager.StartTransaction())
{
BlockTable bt = (BlockTable)tr.GetObject(db.BlockTableId, OpenMode.ForRead, false);
BlockTableRecord btr = (BlockTableRecord)tr.GetObject(bt[BlockTableRecord.ModelSpace], OpenMode.ForWrite, false);
DBObject objPick = tr.GetObject(oid, OpenMode.ForRead);
Entity objEnt = objPick as Entity;
string sLayer = objEnt.Layer.ToString();
ObjectIdCollection oDBO = CADops.SelectAllPolyline(sLayer);
List <string> data = new List <string>();
foreach (ObjectId id in oDBO)
{
//get 3d polyline
Polyline3d poly3d = tr.GetObject(id, OpenMode.ForRead) as Polyline3d;
//get points on 3d polyline
List <Point3d> pts0 = CADops.GetPointsFrom3dPolyline(poly3d, doc);
List <Point3d> pts = pts0.Distinct(new PointComparer(3)).ToList();
//get vectors on points
List <Vector3d> vectorAlongPath = new List <Vector3d>();
List <Vector3d> vectors = CADops.getVectors(pts, doc, ref vectorAlongPath);
List <LoftProfile> loftProfiles = new List <LoftProfile>();
for (int i = 0; i < pts.Count(); i = i + 4)
{
Point3d pt = pts[i];
//ed.WriteMessage($"TUN_CreateTBM => {pt.X}, {pt.Y}, {pt.Z}; ");
Vector3d v = vectors[i].GetNormal() * TunnelDia / 2;
Matrix3d mat = Matrix3d.Displacement(v);
Point3d npt = pt.TransformBy(mat);
//create a 2d line in XY plane
Line ln = new Line(npt, pt);
btr.AppendEntity(ln);
tr.AddNewlyCreatedDBObject(ln, true);
double ang = ln.Angle; //vectors[i].GetAngleTo(curUCS.Xaxis);
//ed.WriteMessage($"angle {ang}\n");
Region acRegion = new Region();
try
{
using (Circle acCirc = new Circle())
{
acCirc.Center = new Point3d(pt.X, pt.Y, pt.Z);
acCirc.Radius = TunnelDia / 2;
acCirc.TransformBy(Matrix3d.Rotation(Angle.angToRad(90), curUCS.Xaxis, acCirc.Center));
acCirc.TransformBy(Matrix3d.Rotation(ang, curUCS.Zaxis, acCirc.Center));
// Add the new object to the block table record and the transaction
btr.AppendEntity(acCirc);
tr.AddNewlyCreatedDBObject(acCirc, true);
DBObjectCollection acDBObjColl = new DBObjectCollection();
acDBObjColl.Add(acCirc);
DBObjectCollection myRegionColl = new DBObjectCollection();
myRegionColl = Region.CreateFromCurves(acDBObjColl);
Region acRegion1 = myRegionColl[0] as Region;
// Add the new object to the block table record and the transaction
btr.AppendEntity(acRegion1);
tr.AddNewlyCreatedDBObject(acRegion1, true);
LoftProfile lp1 = new LoftProfile(acRegion1);
loftProfiles.Add(lp1);
}
}
catch (Autodesk.AutoCAD.Runtime.Exception ex) { ed.WriteMessage(ex.ToString()); }
}
LoftProfile[] lps = new LoftProfile[loftProfiles.Count()];
for (int i = 0; i < loftProfiles.Count(); i++)
{
lps[i] = loftProfiles[i];
}
try
{
// =========== create loft solid
Solid3d sol = new Solid3d();
sol.SetDatabaseDefaults();
LoftOptions lpOptions = new LoftOptions();
//LoftProfile lpGuide = new LoftProfile(acLine);//create loft profile
//LoftProfile[] guideToLoft = new LoftProfile[1] { lpGuide };
//sol.CreateLoftedSolid(lps, guideToLoft, null, lpOptions);//guide to loft can not be null, this parameter is not optional
sol.CreateLoftedSolid(lps, null, null, lpOptions);
sol.Layer = sLayer;//assign layer
//=============save it into database
Autodesk.AutoCAD.DatabaseServices.Entity ent = sol;//this is created for using hyperlink
btr.AppendEntity(ent);
tr.AddNewlyCreatedDBObject(ent, true);
}
catch (Autodesk.AutoCAD.Runtime.Exception ex) { ed.WriteMessage(ex.ToString()); }
}
tr.Commit();
}
}
}
static public void ListVertices()
{
Document doc = Application.DocumentManager.MdiActiveDocument;
Editor ed = doc.Editor;
Database db = doc.Database;
HostApplicationServices hs = HostApplicationServices.Current;
string path = hs.FindFile(doc.Name, doc.Database, FindFileHint.Default);
//ed.WriteMessage("\nTo List all the points in the 3D string, type command LISTVERTICES.");
long adjustX = 0;
long adjustY = 0;
long scaler = 1;
PromptEntityResult per = ed.GetEntity("Select polylines");
ObjectId oid = per.ObjectId;
if (per.Status == PromptStatus.OK)
{
Transaction tr = db.TransactionManager.StartTransaction();
DBObject objPick = tr.GetObject(oid, OpenMode.ForRead);
Entity objEnt = objPick as Entity;
string sLayer = objEnt.Layer.ToString();
path += "_" + sLayer;
ObjectIdCollection oDBO = CADops.SelectAllPolyline(sLayer);
//Handle handseed = db.Handseed;
try
{
using (tr)
{
long shtNo = 1;
List <string> data = new List <string>();
foreach (ObjectId id in oDBO)
{
//DBObject obj = tr.GetObject(per.ObjectId, OpenMode.ForRead);
DBObject obj = tr.GetObject(id, OpenMode.ForRead);
Entity ent = obj as Entity;
string layerName = ent.Layer.ToString();
//create list for storing x,y,z value of point
string StrId = "";
// If a "lightweight" (or optimized) polyline
Autodesk.AutoCAD.DatabaseServices.Polyline lwp = obj as Autodesk.AutoCAD.DatabaseServices.Polyline;
if (lwp != null)
{
//StrId = lwp.ObjectId.ToString();
StrId = lwp.Handle.ToString();
ed.WriteMessage("\n" + StrId);
// Use a for loop to get each vertex, one by one
int vn = lwp.NumberOfVertices;
for (int i = 0; i < vn; i++)
{
// Could also get the 3D point here
Point2d pt = lwp.GetPoint2dAt(i);
string temp = CADops.scaleNmove(pt.X, adjustX, scaler) + "," + CADops.scaleNmove(pt.Y, adjustY, scaler) + ", ," + StrId + "," + layerName;
data.Add(temp);
ed.WriteMessage("\n" + pt.ToString());
}
}
else
{
// If an old-style, 2D polyline
Polyline2d p2d = obj as Polyline2d;
if (p2d != null)
{
StrId = p2d.Handle.ToString();
ed.WriteMessage("\n" + StrId);
// Use foreach to get each contained vertex
foreach (ObjectId vId in p2d)
{
Vertex2d v2d = (Vertex2d)tr.GetObject(vId, OpenMode.ForRead);
string temp = CADops.scaleNmove(v2d.Position.X, adjustX, scaler) + "," + CADops.scaleNmove(v2d.Position.Y, adjustY, scaler) + "," + v2d.Position.Z * scaler + "," + StrId + "," + layerName;
data.Add(temp);
ed.WriteMessage("\n" + v2d.Position.ToString());
}
}
else
{
// If an old-style, 3D polyline
Polyline3d p3d = obj as Polyline3d;
if (p3d != null)
{
StrId = p3d.Handle.ToString();
ed.WriteMessage("\n" + StrId);
// Use foreach to get each contained vertex
foreach (ObjectId vId in p3d)
{
PolylineVertex3d v3d = (PolylineVertex3d)tr.GetObject(vId, OpenMode.ForRead);
string temp = CADops.scaleNmove(v3d.Position.X, adjustX, scaler) + "," + CADops.scaleNmove(v3d.Position.Y, adjustY, scaler) + "," + v3d.Position.Z * scaler + "," + StrId + "," + layerName;
data.Add(temp);
ed.WriteMessage("\n" + v3d.Position.ToString());
}
}
}
}
//create dataarray to populate in excel
int no = data.Count;
ed.WriteMessage("\n Number of Point:" + no);
//createCSV(data, shtNo);
shtNo += 1;
}
Excel.createCSV(data, path);
// Committing is cheaper than aborting
ed.WriteMessage("\ncsv file has been created under path " + path);
ed.WriteMessage("\nData format: Easting, Northing, Elevation, String ID, Layer Name.");
tr.Commit();
}
}
catch { }
}
}
//project to level 0 to find the intersect point
static public void CheckElevationDifference()
{
string msg = string.Empty;
double checkDistance = 0.55;
Document doc = Application.DocumentManager.MdiActiveDocument;
Editor ed = doc.Editor;
Database db = doc.Database;
HostApplicationServices hs = HostApplicationServices.Current;
string outputPath = hs.FindFile(doc.Name, doc.Database, FindFileHint.Default);
string logPath = outputPath + "_log.log";
System.IO.File.WriteAllText(logPath, msg);
//get user input
//select point input file
ed.WriteMessage("Select the Excel file contains the door step points");
string excelPath = FileOps.SelectFile();
#region get Excel Sheet Name
PromptStringOptions pStrOpts = new PromptStringOptions("\nEnter Sheet Name: ");
pStrOpts.AllowSpaces = true;
PromptResult pStrRes = doc.Editor.GetString(pStrOpts);
string shtName = pStrRes.StringResult;
#endregion
List <Point3d> doorStepPts = Excel.getAllpoint(doc, excelPath, shtName);
PromptEntityResult per = ed.GetEntity("Select polylines");
ObjectId oid = per.ObjectId;
ObjectIdCollection oDBO = new ObjectIdCollection();
List <Polyline3d> poly3ds = new List <Polyline3d>();
List <Polyline3d> poly3ds0 = new List <Polyline3d>();
if (per.Status == PromptStatus.OK)
{
Transaction tr = db.TransactionManager.StartTransaction();
BlockTable bt = (BlockTable)tr.GetObject(db.BlockTableId, OpenMode.ForRead, false);
BlockTableRecord btr = (BlockTableRecord)tr.GetObject(bt[BlockTableRecord.ModelSpace], OpenMode.ForWrite, false);
DBObject objPick = tr.GetObject(oid, OpenMode.ForRead);
Entity objEnt = objPick as Entity;
string sLayer = objEnt.Layer.ToString();
oDBO = CADops.SelectAllPolyline(sLayer);
foreach (ObjectId id in oDBO)
{
Polyline3d pl = new Polyline3d();
System.IO.File.AppendAllText(logPath, $"{id}\n");
poly3ds0.Add(CADops.CreatePolylineOnXYPlane(doc, id, ref pl));
poly3ds.Add(pl);
}
}
List <Point3d> output = new List <Point3d>();
List <string> data = new List <string>();
if (poly3ds.Count() > 0)
{
Transaction tr = db.TransactionManager.StartTransaction();
BlockTable bt = (BlockTable)tr.GetObject(db.BlockTableId, OpenMode.ForRead, false);
BlockTableRecord btr = (BlockTableRecord)tr.GetObject(bt[BlockTableRecord.ModelSpace], OpenMode.ForWrite, false);
using (tr)
{
List <Vector3d> vectorsAlongPath = new List <Vector3d>();
List <Vector3d> vectors = CADops.getVectors(doorStepPts, doc, ref vectorsAlongPath);
//foreach (Point3d pt in doorStepPts)
for (int i = 0; i < doorStepPts.Count(); i++)
{
Point3d pt = doorStepPts[i];
Point3d pt0 = new Point3d(pt.X, pt.Y, 0);
Vector3d v = vectors[i].GetNormal() * 5;
Matrix3d mat = Matrix3d.Displacement(v);
Point3d npt = pt0.TransformBy(mat);
v = vectors[i].GetNormal() * -5;
mat = Matrix3d.Displacement(v);
Point3d npt2 = pt0.TransformBy(mat);
//create a 2d line in XY plane
Line ln = new Line(npt, npt2);
btr.AppendEntity(ln);
tr.AddNewlyCreatedDBObject(ln, true);
msg = $"pt => {pt.X}, {pt.Y}, {pt.Z}\n";
#region get intersect point from point to polyline
Point3d ptNearest = Point3d.Origin;
for (int j = 0; j < poly3ds0.Count(); j++)
{
Polyline3d p3d0 = poly3ds0[j];
Polyline3d p3d = poly3ds[j];
Point3d ptTemp = new Point3d();
#region get the alignment object and find the nearest point to the nominated point
Point3dCollection pts3D = new Point3dCollection();
p3d0.IntersectWith(ln, Intersect.OnBothOperands, pts3D, IntPtr.Zero, IntPtr.Zero);
try
{
if (pts3D.Count > 0)
{
double para = p3d0.GetParameterAtPoint(pts3D[0]);
//ed.WriteMessage($"{pts3D[0]}, {para}\n");
ptTemp = p3d.GetPointAtParameter(para);
}
}
catch { }
#region get the point with lower Z
if (ptNearest == Point3d.Origin)
{
ptNearest = ptTemp;
}
else
{
if (ptNearest.Z > ptTemp.Z)
{
ptNearest = ptTemp;
}
}
#endregion
}
#endregion
#endregion
msg += $"ptNearest: {ptNearest.X}, {ptNearest.Y}, {ptNearest.Z}\n";
try
{
double diff = ptNearest.Z - pt.Z;
if (Math.Abs(diff) <= checkDistance)
{
Point3d newPt = new Point3d(pt.X, pt.Y, ptNearest.Z + checkDistance);
output.Add(newPt);
data.Add($"{newPt.X},{newPt.Y},{newPt.Z}, less, {diff}, {ptNearest.Z}");
msg += $", Z diff: {diff} => less than {checkDistance} => {newPt.X}, {newPt.Y}, {newPt.Z}\n\n";
}
else
{
Point3d newPt = pt;
output.Add(newPt);
data.Add($"{newPt.X},{newPt.Y},{newPt.Z}, more, {diff}, {ptNearest.Z}");
msg += $", Z diff: {diff} => more than {checkDistance} => {newPt.X}, {newPt.Y}, {newPt.Z}\n\n";
}
}
catch { }
System.IO.File.AppendAllText(logPath, msg);
}
tr.Commit();
}
}
Excel.createCSV(data, outputPath);
ed.WriteMessage("\ncsv file has been created under path " + outputPath);
}