*Rhino*

## Python and C# Code for Rhino ##

benches.py

from rhino3dm import *

def move(input):
    vectH = Vector3d(2,0,0)
    vectV = Vector3d(0,1,0)
    xform = Transform.Translation(vectH)
    yform = Transform.Translation(vectV)
    b = input.Transform(xform)
    c = b.Transform(yform)
    print (b)
    print(c)
    return b,c





a = Point3d(0,0,0)

print(f'start {a}')

i = 0

while i<4:
    print (i)
    result = move(a)
    i+=1
    a = result[1]

CreateTrenches.py

import rhinoscriptsyntax as rs
import scriptcontext as sc
import Rhino.Geometry as rg
import Rhino
import math


#2015, Konrad Sobon, archi-lab package
def groupCurves(Brep_List): 
    ignore_distance = 0.1 # Assume points this close or closer to each other are touching 
    Grouped_Lines = [] 
    Queue = set() 
    while Brep_List: 
        Shape = [] 
        Queue.add(Brep_List.pop()) # Move a line from the Brep_List to our queue 
        while Queue: 
            Current_Brep = Queue.pop() 
            Current_Brep_Pts = GetExtrusionEndPoints(Current_Brep)
            Shape.append(Current_Brep) 
            for Potential_Match in Brep_List: 
                #Points = (Potential_Match.PointAtStart, Potential_Match.PointAtEnd)
                Points = GetExtrusionEndPoints(Potential_Match)
                for P1 in Points: 
                    #for P2 in (Current_Brep.PointAtStart, Current_Brep.PointAtEnd): 
                    for P2 in (Current_Brep_Pts): 
                        distance = rs.Distance(P1,P2) 
                        if distance <= ignore_distance: 
                            Queue.add(Potential_Match) 
            Brep_List = [item for item in Brep_List if item not in Queue]
        Grouped_Lines.append(Shape) 
    return Grouped_Lines

def GetExtrusionEndPoints(objId):
    obj = rs.coercegeometry(objId)
    objType = str(obj.GetType())
    Points = ()
    if objType == "Rhino.Geometry.Brep":
        line = find_average_line_between_vertices(objId)
        Points = (line.PointAt(0), line.PointAt(1))
    elif objType == "Rhino.Geometry.Mesh":
        brep = rg.Brep.CreateFromMesh(obj, True)
        line = find_average_line_between_vertices(brep)
        Points = (line.PointAt(0), line.PointAt(1))
    elif objType == "Rhino.Geometry.InstanceReferenceGeometry":
        rhinoObj = Rhino.RhinoDoc.ActiveDoc.Objects.FindId(objId)
        rhinoBrep = None
        for obj in rhinoObj.GetSubObjects():
            rhinoBrep = Rhino.Geometry.Brep.CreateFromMesh(obj.MeshGeometry, True)
        line = find_average_line_between_vertices(rhinoBrep)
        Points = (line.PointAt(0), line.PointAt(1))
        """
        atts = rhinoObj.Attributes.Duplicate()
        for obj in rhinoObj.GetSubObjects():
            geom = Rhino.Geometry.Brep.CreateFromMesh(obj.MeshGeometry, True)
            if geom:
                sc.doc.Objects.AddBrep(geom, atts)
        """
    else:
        stPt = obj.PathStart
        endPt = obj.PathEnd
        Points = (stPt, endPt)
    return Points

def find_average_line_between_vertices(brep_id):
    obj = rs.coercegeometry(brep_id)
    vertices = obj.Vertices         
    pts = [v.Location for v in vertices]    
    average_line = rg.Line.TryFitLineToPoints(pts)[1]
    #sc.doc.Objects.AddLine(average_line)
    return average_line

def ConvertToBrep(objId):
    obj = rs.coercegeometry(objId)
    if str(obj.GetType()) == "Rhino.Geometry.Brep":
        return obj
    else:
        return obj.ToBrep()
    
def OffsetCurve(id, cover):
    obj = rs.coercegeometry(id)
    stPt = obj.PathStart
    endPt = obj.PathEnd
    line = rg.LineCurve(stPt, endPt)
    dia = float(rs.GetUserText(id, 'Dimensions|Outside Diameter'))/1000
    offsetDistance = dia/2+cover
    zMoveVect = rg.Vector3d(0,0,-offsetDistance)
    line.Transform(rg.Transform.Translation(zMoveVect))
    return [line, offsetDistance]

def ExtrudeElbow(items, extrusionH, tolerance, cover):
    offset1 = OffsetCurve(items[0], cover)
    curve1 = offset1[0]
    offsetDistance1 = offset1[1]
    
    curve2 = OffsetCurve(items[1], cover)[0]
    st = curve1.PointAtStart
    end = curve1.PointAtEnd
    pts = [curve2.PointAtStart, curve2.PointAtEnd]
    
    closest0 = pts[rs.PointArrayClosestPoint(pts, st)]
    closest1 = pts[rs.PointArrayClosestPoint(pts, end)]
    closest = None
    arc = None
    
    if abs(rs.Distance(st, closest0)) < abs(rs.Distance(end, closest1)):
        closest = closest0   
        p = rs.CurveClosestPoint(curve1, st)
        arc = rs.AddArcPtTanPt(st, -rs.CurveTangent(curve1, p), closest)
    else:
        closest = closest1
        p = rs.CurveClosestPoint(curve1, end)
        arc = rs.AddArcPtTanPt(end, rs.CurveTangent(curve1, p), closest)
    arcObj = rs.coercecurve(arc)
    extrusion = rg.Extrusion.CreateExtrusion(arcObj, rg.Vector3d(0,0,extrusionH))
    brep = rg.Brep.CreateFromOffsetFace(extrusion.ToBrep().Faces[0], offsetDistance1, tolerance, True, True)
    #sc.doc.Objects.AddBrep(brep)
    return brep

def CreateExtrusions(ids, extrusionH, cover, tolerance ):
    #separate extrusions from elbows
    extrusionsIds = []
    elbows = []
    extrudedBreps = []

    for id in ids:
        obj = rs.coercegeometry(id)
        objType = str(obj.GetType())
        if objType == "Rhino.Geometry.Brep":
            elbows.append(obj)
        elif objType == "Rhino.Geometry.Mesh":
            brep = rg.Brep.CreateFromMesh(obj, True)
            elbows.append(brep)
        elif objType == "Rhino.Geometry.InstanceReferenceGeometry":
            rhinoObj = Rhino.RhinoDoc.ActiveDoc.Objects.FindId(id)
            brep = None
            for obj in rhinoObj.GetSubObjects():
                brep = Rhino.Geometry.Brep.CreateFromMesh(obj.MeshGeometry, True)
            elbows.append(brep)
        else:
            stPt = obj.PathStart
            endPt = obj.PathEnd
            extrusionsIds.append(id)
            #line_id = rs.AddLine(stPt, endPt)
            line = rg.LineCurve(stPt, endPt)
            #print line
            dia = float(rs.GetUserText(id, 'Dimensions|Outside Diameter'))/1000
            offsetDistance = dia/2+cover
            zMoveVect = rg.Vector3d(0,0,-offsetDistance)
            line.Transform(rg.Transform.Translation(zMoveVect))
            #extrusion = rg.Extrusion.Create(line, 10, True)
            extrusion = rg.Extrusion.CreateExtrusion(line, rg.Vector3d(0,0,extrusionH))
            brep = rg.Brep.CreateFromOffsetFace(extrusion.ToBrep().Faces[0], offsetDistance, tolerance, True, True)
            extrudedBreps.append(brep)
            #sc.doc.Objects.AddBrep(brep)
    counter = 0
    for elbow in elbows:
        #center = elbow.GetBoundingBox(rg.Plane.WorldXY).Center
        connectedExtrusions = []
        for extrusionId in extrusionsIds:
            extrusion = rs.coercegeometry(extrusionId)
            intersCurves = rg.Intersect.Intersection.BrepBrep(elbow, extrusion.ToBrep(), tolerance)[1]
            if len(intersCurves)>0:
                #excenter = extrusion.ToBrep().GetBoundingBox(rg.Plane.WorldXY).Center
                connectedExtrusions.append(extrusionId)
        extrudedBreps.append(ExtrudeElbow(connectedExtrusions, extrusionH, tolerance, cover))
        #add boolean union
        counter += 1
    return extrudedBreps


ids = rs.GetObjects("select objects", rs.filter.extrusion+rs.filter.polysurface+rs.filter.mesh+rs.filter.instance)

rs.EnableRedraw(False)

groupedCurves = groupCurves(ids)

index = 0
t = sc.doc.ModelAbsoluteTolerance
c = 0.05    #conduits cover
H = 3       #extrusion H

for g in groupedCurves:
    print (str(index))
    breps = CreateExtrusions(g, H, c, t)
    print (breps)
    #rhinoObj = Rhino.RhinoDoc.ActiveDoc.Objects.FindId(breps[0])    
    union = rg.Brep.CreateBooleanUnion(breps, 0.1)
    if union:
        sc.doc.Objects.AddBrep(union[0])
    else:
        for b in breps:
            sc.doc.Objects.AddBrep(b)
    #union = rg.Brep.JoinBreps(breps, 0.1)
    #print (union)
    
    index +=1  

rs.EnableRedraw(True)

CurvesToLayers.py

import rhinoscriptsyntax as rs
import Rhino
import scriptcontext as sc
import math
from System.Collections.Generic import List
import System
import System.Collections.Generic.IEnumerable as IEnumerable
import hashlib

objects = rs.SelectedObjects()

existing_layers = rs.LayerNames() or []




for obj in objects:
    layer_name = rs.ObjectName(obj).split(' ')[2]
    if layer_name in existing_layers:
        print (layer_name)
        rs.ObjectLayer(obj, layer_name)
    else:
        color = rs.CreateColor([0,0,0])
        if 'N24' in layer_name:
            color = rs.CreateColor([226,66,68])
        elif 'N28' in layer_name:
            color = rs.CreateColor([213,7,174])
        elif 'N16' in layer_name:
            color = rs.CreateColor([79,163,28])
        elif 'N20' in layer_name:
            color = rs.CreateColor([15,54,24])
        elif 'N12' in layer_name:
            color = rs.CreateColor([245,95,163])
        rs.AddLayer(layer_name)
        rs.LayerColor(layer_name, color)
        existing_layers.append(layer_name)
        rs.ObjectLayer(obj, layer_name)

explodeBlocksWithAttr.cs

//https://github.com/mcneel/rhino-developer-samples/blob/7/rhinocommon/cs/SampleCsCommands/SampleCsExplodeBlock.cs
    if (x){
      // Select a block instance
      var go = new GetObject();
      go.SetCommandPrompt("Select a block instance to explode");
      go.GeometryFilter = ObjectType.InstanceReference;
      go.Get();
        

        
      var iref = go.Object(0).Object();
      var xform = Transform.Identity;

      if (!(iref is InstanceObject)){
        A = "Error";
      }

      var atts = iref.Attributes.Duplicate();
      atts.ObjectId = Guid.Empty;


      var objects = iref.GetSubObjects();

      foreach (var obj in objects){
        var geom = obj.DuplicateGeometry();
        doc.Objects.Add(geom, atts);
      }

      doc.Objects.Delete(go.Object(0), false);
      doc.Views.Redraw();

      A = objects;
    }
    else {
      A = "waiting...";
    }

explodeBlocksWithAttr.py

import rhinoscriptsyntax as rs
import scriptcontext as sc
import Rhino
import math

def explode_block_instance(goes):
    print len(goes.Objects())
    for go in goes.Objects():
        iref = go.Object()
        xform = Rhino.Geometry.Transform.Identity
        
        if not isinstance(iref, Rhino.DocObjects.InstanceObject):
            return "Error"
    
        atts = iref.Attributes.Duplicate()
        objects = iref.GetSubObjects()
    
        for obj in objects:
            geom = Rhino.Geometry.Brep.CreateFromMesh(obj.MeshGeometry, True)
            if geom:
                sc.doc.Objects.AddBrep(geom, atts)
    
        #sc.doc.Objects.Delete(go, False)
        rs.DeleteObject(go)
        #sc.doc.Views.Redraw()
    return "done"

selection = Rhino.Input.Custom.GetObject()
selection.SetCommandPrompt("Select a block instance to explode")
selection.GeometryFilter = Rhino.DocObjects.ObjectType.InstanceReference
selection.GetMultiple(1,0)

#Rhino.DocObjects.Tables.ObjectTable.Delete(go.Objects()[0].Object()., False)

explode_block_instance(selection)

#equivalent to
ids = rs.GetObject("select extrusion", rs.filter.instance)
rhinoObj = Rhino.RhinoDoc.ActiveDoc.Objects.FindId(ids);
atts = rhinoObj.Attributes.Duplicate()
for obj in rhinoObj.GetSubObjects():
    geom = Rhino.Geometry.Brep.CreateFromMesh(obj.MeshGeometry, True)
    if geom:
        sc.doc.Objects.AddBrep(geom, atts)
        
rs.Redraw()

FindConnectedLines.py

import rhinoscriptsyntax as rs
import scriptcontext as sc
import Rhino.Geometry as rg

selectedId = rs.GetObject('select start')
selectedObj = rs.coercegeometry(selectedId)
print (selectedObj)

rs.EnableRedraw(False)

#all_object_ids = rs.VisibleObjects()
all_object_ids = rs.AllObjects()

lista = []

for objId in all_object_ids:
    line = rs.coercegeometry(objId)
    objType = str(line.GetType())
    if objType == "Rhino.Geometry.LineCurve":
        st = line.PointAtStart
        end = line.PointAtEnd
        length = rg.LineCurve.GetLength(line, 0.001)
        midPt = rg.LineCurve.PointAtLength(line, length*0.5)
        #rs.AddTextDot(i, midPt)
        lista.append((objId, st, end))
    
# print('lista')
# for l in lista:
#     print (l)



connected = []
queue = set()

for s in lista:
    if s[0] == selectedId:
        p = lista.pop(lista.index(s))
        queue.add(p)
        connected.append(p)
        break

print ('connected', connected)

while queue:
    current = queue.pop()
    currentPoints = (current[1], current[2])
    for s in lista:
        points = (s[1],s[2])
        for p1 in currentPoints:
            for p2 in points:
                distance = rg.Point3d.DistanceTo(p1,p2)
                if distance < 0.02 and s not in connected:
                    queue.add(s)
                    connected.append(s)




# print ('result', connected)
print (len(connected))


for c in connected:
    rs.SelectObject(c[0])

FindConnectedObjects.py

import rhinoscriptsyntax as rs
import scriptcontext as sc
import Rhino.Geometry as rg
import Rhino

def GetExtrusionEndPoints(objId):
    obj = rs.coercegeometry(objId)
    objType = str(obj.GetType())
    Points = ()
    if objType == "Rhino.Geometry.Brep":
        line = find_average_line_between_vertices(objId)
        Points = (line.PointAt(0), line.PointAt(1))
    elif objType == "Rhino.Geometry.Mesh":
        brep = rg.Brep.CreateFromMesh(obj, True)
        line = find_average_line_between_vertices(brep)
        Points = (line.PointAt(0), line.PointAt(1))
    elif objType == "Rhino.Geometry.InstanceReferenceGeometry":
        rhinoObj = Rhino.RhinoDoc.ActiveDoc.Objects.FindId(objId)
        rhinoBrep = None
        for obj in rhinoObj.GetSubObjects():
            rhinoBrep = Rhino.Geometry.Brep.CreateFromMesh(obj.MeshGeometry, True)
        line = find_average_line_between_vertices(rhinoBrep)
        Points = (line.PointAt(0), line.PointAt(1))
    elif objType == "Rhino.Geometry.Extrusion":
        stPt = obj.PathStart
        endPt = obj.PathEnd
        Points = (stPt, endPt)
    return Points

def find_average_line_between_vertices(brep_id):
    obj = rs.coercegeometry(brep_id)
    vertices = obj.Vertices         
    pts = [v.Location for v in vertices]    
    average_line = rg.Line.TryFitLineToPoints(pts)[1]
    #sc.doc.Objects.AddLine(average_line)
    return average_line

selectedId = rs.GetObject('select start')
selectedObj = rs.coercegeometry(selectedId)
print (selectedObj)

rs.EnableRedraw(False)

all_object_ids = rs.VisibleObjects()
#all_object_ids = rs.AllObjects()

lista = []

filterTypes= ["Rhino.Geometry.Extrusion", "Rhino.Geometry.Brep", "Rhino.Geometry.InstanceReferenceGeometry"]

for objId in all_object_ids:
    obj = rs.coercegeometry(objId)
    objType = str(obj.GetType())
    points = GetExtrusionEndPoints(objId)
    if points:
        lista.append((objId, points[0], points[1]))



# for l in lista[0:10]:
#     print (l)


connected = []
queue = set()

for s in lista:
    if s[0] == selectedId:
        p = lista.pop(lista.index(s))
        queue.add(p)
        connected.append(p)
        break

print ('connected', connected)

while queue:
    current = queue.pop()
    currentPoints = (current[1], current[2])
    for s in lista:
        points = (s[1],s[2])
        for p1 in currentPoints:
            for p2 in points:
                distance = rg.Point3d.DistanceTo(p1,p2)
                if distance < 0.08 and s not in connected:
                    queue.add(s)
                    connected.append(s)




# print ('result', connected)
print (len(connected))


for c in connected:
    rs.SelectObject(c[0])

flattenMesh.py

"""Provides a scripting component.
    Inputs:
        x: The x script variable
        y: The y script variable
    Output:
        a: The a output variable"""

__author__ = "G"
__version__ = "2023.06.06"

import rhinoscriptsyntax as rs
import Rhino
from math import pow
import ghpythonlib.treehelpers as th

newVertexes = []
lines = []
newCurves = [crv[0]]
currentCrv = crv[0]

for i in range(0,len(crv)-1):
    print '##################curve ' + str(i)
    subVert = []
    subLines = []
    divisionPoints = None
    if rs.CurveLength(crv[i]) > 1.1:
        divisionPoints = rs.DivideCurveEquidistant(crv[i+1],1)
    else:
        divisionPoints = rs.DivideCurve(crv[i+1],3)
    #for pt in crv[i+1].ToPolyline().ToArray():
    for pt in divisionPoints:
        #point on previous curve
        pt0 = currentCrv.PointAt(rs.CurveClosestPoint(currentCrv, pt))        
        run = pow( pow(pt0.X-pt.X,2) + pow(pt0.Y-pt.Y,2) , 0.5) 
        rise = abs( pt0.Z - pt.Z )
        print 'initial rise ' + str(round(run,2))
        slope = rise/run*100
        print 'initial slope ' + str(round (slope,2))+ '%'
        newPt = pt
        counter = 0
        while counter < 100:
            print '****************************************counter' + str(counter)
            print '*****************new slope ' + str(slope)
            if (slope > 3):
                print 'slope higher than 3%!!!'
                #newPt = rs.coerce3dpoint(rs.AddPoint(pt.X, pt.Y, pt.Z + run * 0.03))
                newPt = rs.coerce3dpoint(rs.AddPoint(newPt.X, newPt.Y, newPt.Z + 0.01))
                print 'point Z ' +  str(newPt.Z)
                newRise = abs( pt0.Z - newPt.Z )
                print 'new rise' + str(newRise)
                slope = newRise/run*100
                print 'new slope ' + str(round(slope,2)) + '%'
            counter +=1
        subVert.append(newPt)
        subLines.append(rs.AddLine(pt0, newPt))
    #print subVert
    currentCrv = rs.coercecurve( rs.AddPolyline(subVert) )
    newVertexes.append(subVert)
    lines.append(subLines)
    newCurves.append(currentCrv)

newVertexes = th.list_to_tree(newVertexes, source=[0,0])
lines = th.list_to_tree(lines, source=[0,0])


newPt = newVertexes
slope = lines
curves = newCurves

hideRhinoSkin.cs

            Rhino.ApplicationSettings.AppearanceSettings.MenuVisible = false;
            Rhino.ApplicationSettings.AppearanceSettings.ShowSideBar = false;
            Rhino.ApplicationSettings.AppearanceSettings.ShowOsnapBar = false;
            Rhino.ApplicationSettings.AppearanceSettings.ShowStatusBar = false;
            Rhino.ApplicationSettings.AppearanceSettings.ShowFullPathInTitleBar = false;

            for (int i = 0; i < RhinoApp.ToolbarFiles.Count; i++)
            {
                RhinoApp.ToolbarFiles[i].Close(false);
            }

            foreach (var panel in Rhino.UI.Panels.GetOpenPanelIds())
            {
                Rhino.UI.Panels.ClosePanel(panel);
            }


            var view = Rhino.RhinoDoc.ActiveDoc.Views.ActiveView;

            int index = Rhino.RhinoDoc.ActiveDoc.NamedViews.FindByName("Perspective");

            Rhino.RhinoDoc.ActiveDoc.NamedViews.Restore(index, view.ActiveViewport);

            view.Maximized = true;
            view.Redraw();

            RhinoApp.RunScript("-_CommandPrompt _Show _No _Enter", false);
            //Rhino.ApplicationSettings.AppearanceSettings.CommandPromptPosition = Rhino.ApplicationSettings.CommandPromptPosition.Hidden;

ImportInLayer.py

#credits: https://discourse.mcneel.com/t/import-dxf-files-in-rhino-and-reference-the-geometry/78832

import rhinoscriptsyntax as rs
import scriptcontext as sc
import Rhino,os

selectedFiles = rs.OpenFileNames("Import files", None,r"C:\OneDrive");

for file in selectedFiles:
    
    filename = file.split("\\")[-1]
    
    fullpath=file
    
    obj_count=0
    file_count=0
    raise_copy_count=False
    
    #get existing layer set
    exist_layers=set(rs.LayerNames())
    #import file
    print "Importing " + filename
    rs.Command('_-Import "{}" _Enter'.format(fullpath), False)
    if rs.LastCommandResult()==0: file_count+=1
    lco = rs.LastCreatedObjects(False)
    if lco: obj_count+=len(lco)
    #get all newly imported layers
    imp_layers=set(rs.LayerNames()).difference(exist_layers)
    #Some file formats (STEP) can be "layerless"! imp_layers will have len==0
    
    #create new parent layer with filename
    if not rs.IsLayer(filename):
        p_layer=rs.AddLayer(filename)
    else:
        #add suffix to make parent layer unique
        p_layer=rs.AddLayer(AddNumberSuffix(filename,copy_count))
        raise_copy_count=True
        #try to transfer layers to new parent
    if len(imp_layers)==0:
        #need to transfer imported objects to newly created parent layer
        rs.ObjectLayer(lco,p_layer)
    else:
        for imp_layer in imp_layers:
            try:
                if not rs.ParentLayer(imp_layer):
                    rs.ParentLayer(imp_layer,p_layer)
            except: pass
    rs.ExpandLayer(filename, False)
    rs.ZoomExtents()

joinMeshesByCentreline.py

import rhinoscriptsyntax as rs
import Rhino
import scriptcontext as sc
import math
from System.Collections.Generic import List
import System
import System.Collections.Generic.IEnumerable as IEnumerable
from itertools import groupby
from operator import attrgetter
from collections import defaultdict

class ReoCL:
	def __init__(self, curve, name):
		self.crv = curve
		self.n = name

class MeshCL:
	def __init__(self, mesh, name):
		self.mesh = mesh
		self.n = name

rs.EnableRedraw(False)

view = rs.CurrentView()

all_curves = rs.ObjectsByType(4, select=False)
visible_centrelines = [c for c in all_curves if not rs.IsObjectHidden(c)]

joined_curves = rs.JoinCurves(visible_centrelines, True)

reoCentrelines = []

for i in range(0,len(joined_curves)):
    reoCentrelines.append(ReoCL(joined_curves[i], i))

#print ('cls', reoCentrelines)

meshes = rs.ObjectsByType(32, select=False)
visible_meshes = [c for c in meshes if not rs.IsObjectHidden(c)]

namedMeshes = []

for mesh in visible_meshes:
    #print ('mesh:',mesh)
    centroid = rs.MeshVolumeCentroid(mesh)
    rs.AddPoint(centroid)
    shortestDistance = 9999
    curveName = 'ZZZ'
    for reoCL in reoCentrelines:
        param = rs.CurveClosestPoint(reoCL.crv, centroid)
        closestPt = rs.EvaluateCurve(reoCL.crv, param)
        distance = rs.Distance(centroid, closestPt)
        if distance < shortestDistance:
            shortestDistance = distance
            curveName = reoCL.n
    namedMeshes.append(MeshCL(mesh, curveName))
    #print ('shortest distance', curveName)

grouped_meshes = defaultdict(list)
for mesh in namedMeshes:
    grouped_meshes[mesh.n].append(mesh)


for group in grouped_meshes:
    subMesh = []
    for mesh in grouped_meshes[group]:
        subMesh.append(mesh.mesh)
    rs.JoinMeshes(subMesh)
    subMesh = []


rs.EnableRedraw(True)

NavisworksViewpoints.cs

//------------------------------------------------------------------
// NavisWorks Sample code
//------------------------------------------------------------------

// (C) Copyright 2010 by Autodesk Inc.

// Permission to use, copy, modify, and distribute this software in
// object code form for any purpose and without fee is hereby granted,
// provided that the above copyright notice appears in all copies and
// that both that copyright notice and the limited warranty and
// restricted rights notice below appear in all supporting
// documentation.

// AUTODESK PROVIDES THIS PROGRAM "AS IS" AND WITH ALL FAULTS.
// AUTODESK SPECIFICALLY DISCLAIMS ANY IMPLIED WARRANTY OF
// MERCHANTABILITY OR FITNESS FOR A PARTICULAR USE.  AUTODESK
// DOES NOT WARRANT THAT THE OPERATION OF THE PROGRAM WILL BE
// UNINTERRUPTED OR ERROR FREE.
//------------------------------------------------------------------
//
// This sample compares the time it takes to find ModelItems in a 
// model using various search methods.
//
//------------------------------------------------------------------
using System;
using System.Collections.Generic;
using System.Linq;
using System.Text;
using System.Windows.Forms;
using Autodesk.Navisworks.Api;
using Autodesk.Navisworks.Api.Plugins;
using ComApi = Autodesk.Navisworks.Api.Interop.ComApi;
using ComApiBridge = Autodesk.Navisworks.Api.ComApi;

namespace SearchComparisonPlugIn
{
   [PluginAttribute("SearchComparisonPlugIn.SearchComparisonPlugIn", "ADSK",
      DisplayName="Search Comparison",
      ToolTip="Compares the time taken for different search techniques on a Navisworks Document")]
   public class SearchComparisonPlugIn : AddInPlugin
   {

        public static double[] Elements;
        public override int Execute(params string[] parameters)
      {

            //ComApi.InwOpState10 state;
            //state = ComApiBridge.ComApiBridge.State;

            //foreach (ComApi.InwOaPath path in state.CurrentSelection.Paths())
            //{
            //    ComApi.InwOaNode node;
            //    node = path.Nodes().Last() as ComApi.InwOaNode;
            //    MessageBox.Show("UserName=" + node.UserName);
            //}

            ModelItemCollection oModelColl = Autodesk.Navisworks.Api.Application.ActiveDocument.CurrentSelection.SelectedItems;
            //convert to COM selection

            ComApi.InwOpState oState = ComApiBridge.ComApiBridge.State;

            ComApi.InwOpSelection oSel = ComApiBridge.ComApiBridge.ToInwOpSelection(oModelColl);

            CallbackGeomListener callbkListener = new CallbackGeomListener();

            



            foreach (ComApi.InwOaPath3 path in oSel.Paths())
            {
                foreach (ComApi.InwOaFragment3 frag in path.Fragments())
                {
                    // https://forums.autodesk.com/t5/navisworks-api/geometry-in-wrong-location-when-accessing-in-api/m-p/5896205/highlight/true
                    // https://github.com/jotpunktbee/SpeckleNavisworks/blob/08204aed7ee9f99d8b5ba50c5df102a4c95907a2/SpeckleNavisworks/Models/SearchComparisonPlugIn.cs 
                    ComApi.InwLTransform3f3 localToWorld = (ComApi.InwLTransform3f3)(object)frag.GetLocalToWorldMatrix();

                    Array localToWorldMatrix = (Array)(object)localToWorld.Matrix;
                    Elements = ToArray<double>(localToWorldMatrix);

                    callbkListener.LCS2WCS = localToWorld;


                    // generate the primitives
                    frag.GenerateSimplePrimitives(ComApi.nwEVertexProperty.eNORMAL, callbkListener);
                }
            }

            var result = callbkListener.coords;

            Document doc = Autodesk.Navisworks.Api.Application.ActiveDocument;
            Selection storedSelection = new Selection();


            storedSelection = new Selection();
            storedSelection.CopyFrom(doc.CurrentSelection.ToSelection());

            ModelGeometry geom = storedSelection.ExplicitSelection.First.Geometry;

            //int count = geom.PrimitiveCount;

            Point3D bb3dcenter = storedSelection.GetSelectedItems().BoundingBox().Center;



            #region
            StringBuilder output = new StringBuilder(1000);
         output.Append("Searching for Visible ModelItems which contain geometry, returning a ModelItemCollection\n");
         output.Append("Time format is (mm:ss.fff)\n");

         /////////////////////////////////////////////////////////////////////////////
         //iterative method
         DateTime start = DateTime.Now;

         ModelItemCollection searchResults = new ModelItemCollection();
         foreach (ModelItem modelItem in Autodesk.Navisworks.Api.Application.ActiveDocument.Models.CreateCollectionFromRootItems().DescendantsAndSelf)
         {
            //find the model items with the Geometry, not hidden
            if (modelItem.HasGeometry && !modelItem.IsHidden)
               searchResults.Add(modelItem);
         }

         DateTime end = DateTime.Now;
         DateTime taken = new DateTime(end.Ticks - start.Ticks);
         output.Append("Iterative method\t");
         output.Append(taken.ToString("mm:ss.fff"));
         output.Append(" to find ");
         output.Append(searchResults.Count);
         output.Append(" ModelItems\n");

         /////////////////////////////////////////////////////////////////////////////
         //using LINQ style query

         start = DateTime.Now;

         ModelItemCollection searchResults1 = new ModelItemCollection();
         searchResults1.CopyFrom(
            Autodesk.Navisworks.Api.Application.ActiveDocument.Models.RootItemDescendantsAndSelf.
            Where(x =>
               //find the ModelItems with the Geometry, not hidden
               (x.HasGeometry && !x.IsHidden)));

         end = DateTime.Now;
         taken = new DateTime(end.Ticks - start.Ticks);
         output.Append("LINQ query\t");
         output.Append(taken.ToString("mm:ss.fff"));
         output.Append(" to find ");
         output.Append(searchResults1.Count);
         output.Append(" ModelItems\n");

         /////////////////////////////////////////////////////////////////////////////
         //Searching purely with Search class

         start = DateTime.Now;

         //Create a new search
         Search s2 = new Search();

         //Add a search condition for ModelItems with Geometry only...
         s2.SearchConditions.Add(SearchCondition.HasCategoryByName(PropertyCategoryNames.Geometry));

         //...and not hidden
         s2.SearchConditions.Add(SearchCondition.HasPropertyByName(PropertyCategoryNames.Item, DataPropertyNames.ItemHidden).EqualValue(VariantData.FromBoolean(false)));

         //set the selection to everything
         s2.Selection.SelectAll();
         s2.Locations = SearchLocations.DescendantsAndSelf;

         //get the resulting collection by applying this search
         ModelItemCollection searchResults2 = s2.FindAll(Autodesk.Navisworks.Api.Application.ActiveDocument, false);

         end = DateTime.Now;
         taken = new DateTime(end.Ticks - start.Ticks);
         output.Append("Search query\t");
         output.Append(taken.ToString("mm:ss.fff"));
         output.Append(" to find ");
         output.Append(searchResults2.Count);
         output.Append(" ModelItems\n");


         /////////////////////////////////////////////////////////////////////////////
         //Mixture search
         start = DateTime.Now;

         //Create a new search
         Search s3 = new Search();

         //Add a search condition for those ModelItems with Geometry
         s3.SearchConditions.Add(SearchCondition.HasCategoryByName(PropertyCategoryNames.Geometry));

         //set the selection to everything
         s3.Selection.SelectAll();
         s3.Locations = SearchLocations.DescendantsAndSelf;

         //get the resulting collection by applying this search
         ModelItemCollection searchResults3 = new ModelItemCollection();
         searchResults3.CopyFrom(
            s3.FindIncremental(Autodesk.Navisworks.Api.Application.ActiveDocument, false)
            .Where(x =>
               //find the model items, not hidden
               (!x.IsHidden)));

         end = DateTime.Now;
         taken = new DateTime(end.Ticks - start.Ticks);
         output.Append("Mixed query\t");
         output.Append(taken.ToString("mm:ss.fff"));
         output.Append(" to find ");
         output.Append(searchResults3.Count);
         output.Append(" ModelItems\n");

            #endregion

            //MessageBox.Show(output.ToString());
            MessageBox.Show("Ciao da visual studoi");
         return 0;
      }
        public T[] ToArray<T>(Array arr)
        {
            T[] result = new T[arr.Length];
            Array.Copy(arr, result, result.Length);
            return result;
        }

    }

    class CallbackGeomListener : ComApi.InwSimplePrimitivesCB
    {
        public ComApi.InwLTransform3f3 LCS2WCS;
        public List<string> coords = new List<string>();

        
        public void Line(ComApi.InwSimpleVertex v1, ComApi.InwSimpleVertex v2)

        {
            // do your work

            //Array array_v1 = (Array)(object)v1.coord;

            //coords.Add($"{array_v1.GetValue(1)},{array_v1.GetValue(2)},{array_v1.GetValue(3)}");

            //Array array_v2 = (Array)(object)v2.coord;

            //2
            Array array_v1 = (Array)(object)v1.coord;
            Array array_v2 = (Array)(object)v2.coord;


            double w = (SearchComparisonPlugIn.Elements[3] * Convert.ToDouble(array_v1.GetValue(1))) + 
                        (SearchComparisonPlugIn.Elements[7] * Convert.ToDouble(array_v1.GetValue(2))) + 
                        (SearchComparisonPlugIn.Elements[11] * Convert.ToDouble(array_v1.GetValue(3))) + 
                        SearchComparisonPlugIn.Elements[15];

            double newPointX = ((SearchComparisonPlugIn.Elements[0] * Convert.ToDouble(array_v1.GetValue(1))) + (SearchComparisonPlugIn.Elements[4] * Convert.ToDouble(array_v1.GetValue(2))) + (SearchComparisonPlugIn.Elements[8] * Convert.ToDouble(array_v1.GetValue(3))) + SearchComparisonPlugIn.Elements[12]) / w;
            double newPointY = ((SearchComparisonPlugIn.Elements[1] * Convert.ToDouble(array_v1.GetValue(1))) + (SearchComparisonPlugIn.Elements[5] * Convert.ToDouble(array_v1.GetValue(2))) + (SearchComparisonPlugIn.Elements[9] * Convert.ToDouble(array_v1.GetValue(3))) + SearchComparisonPlugIn.Elements[13]) / w;
            double newPointZ = ((SearchComparisonPlugIn.Elements[2] * Convert.ToDouble(array_v1.GetValue(1))) + (SearchComparisonPlugIn.Elements[6] * Convert.ToDouble(array_v1.GetValue(2))) + (SearchComparisonPlugIn.Elements[10] * Convert.ToDouble(array_v1.GetValue(3))) + SearchComparisonPlugIn.Elements[14]) / w;

            double w2 = (SearchComparisonPlugIn.Elements[3] * Convert.ToDouble(array_v2.GetValue(1))) +
            (SearchComparisonPlugIn.Elements[7] * Convert.ToDouble(array_v2.GetValue(2))) +
            (SearchComparisonPlugIn.Elements[11] * Convert.ToDouble(array_v2.GetValue(3))) +
            SearchComparisonPlugIn.Elements[15];

            double newPoint2X = ((SearchComparisonPlugIn.Elements[0] * Convert.ToDouble(array_v2.GetValue(1))) + (SearchComparisonPlugIn.Elements[4] * Convert.ToDouble(array_v2.GetValue(2))) + (SearchComparisonPlugIn.Elements[8] * Convert.ToDouble(array_v2.GetValue(3))) + SearchComparisonPlugIn.Elements[12]) / w2;
            double newPoint2Y = ((SearchComparisonPlugIn.Elements[1] * Convert.ToDouble(array_v2.GetValue(1))) + (SearchComparisonPlugIn.Elements[5] * Convert.ToDouble(array_v2.GetValue(2))) + (SearchComparisonPlugIn.Elements[9] * Convert.ToDouble(array_v2.GetValue(3))) + SearchComparisonPlugIn.Elements[13]) / w2;
            double newPoint2Z = ((SearchComparisonPlugIn.Elements[2] * Convert.ToDouble(array_v2.GetValue(1))) + (SearchComparisonPlugIn.Elements[6] * Convert.ToDouble(array_v2.GetValue(2))) + (SearchComparisonPlugIn.Elements[10] * Convert.ToDouble(array_v2.GetValue(3))) + SearchComparisonPlugIn.Elements[14]) / w2;


            coords.Add($"{newPointX},{newPointY},{newPointZ}");
            coords.Add($"{newPoint2X},{newPoint2Y},{newPoint2Z}");
        
        }



        public void Point(ComApi.InwSimpleVertex v1)
        {

            // do your work

        }



        public void SnapPoint(ComApi.InwSimpleVertex v1)
        {

            // do your work

        }



        public void Triangle(ComApi.InwSimpleVertex v1, ComApi.InwSimpleVertex v2, ComApi.InwSimpleVertex v3)
        {

            // do your work

        }

        


    }
}

objectGuid.py

__author__ = "GBrogiolo"
__version__ = "2023.08.09"

import rhinoscriptsyntax as rs
import Rhino

for i in x:
    rhinoObj = Rhino.RhinoDoc.ActiveDoc.Objects.FindId(i)
    print rhinoObj.GetType()
    if "Extrusion" in str(rhinoObj.GetType()):
        print 'found'
        dia = rs.GetUserText(rhinoObj, "Dimensions|Outside Diameter", False)
    else:
        curve = i

ObjectsToLayers.py

import rhinoscriptsyntax as rs
import Rhino
import scriptcontext as sc
import math
from System.Collections.Generic import List
import System
import System.Collections.Generic.IEnumerable as IEnumerable
import hashlib

objects = rs.SelectedObjects()

existing_layers = rs.LayerNames() or []




for obj in objects:
    layer_name = rs.ObjectName(obj).split(' ')[2]
    if layer_name in existing_layers:
        print (layer_name)
        rs.ObjectLayer(obj, layer_name)
    else:
        color = rs.CreateColor([0,0,0])
        if 'N24' in layer_name:
            color = rs.CreateColor([226,66,68])
        elif 'N28' in layer_name:
            color = rs.CreateColor([213,7,174])
        elif 'N16' in layer_name:
            color = rs.CreateColor([79,163,28])
        elif 'N20' in layer_name:
            color = rs.CreateColor([15,54,24])
        elif 'N12' in layer_name:
            color = rs.CreateColor([245,95,163])
        rs.AddLayer(layer_name)
        rs.LayerColor(layer_name, color)
        existing_layers.append(layer_name)
        rs.ObjectLayer(obj, layer_name)

PipeTrenches.py

import rhinoscriptsyntax as rs
import scriptcontext as sc
import Rhino
import math

def OffsetCurve(id):
    obj = rs.coercegeometry(id)
    stPt = obj.PathStart
    endPt = obj.PathEnd
    line = Rhino.Geometry.LineCurve(stPt, endPt)
    dia = float(rs.GetUserText(id, 'Dimensions|Outside Diameter'))/1000
    offsetDistance = dia/2+cover
    zMoveVect = Rhino.Geometry.Vector3d(0,0,-offsetDistance)
    line.Transform(Rhino.Geometry.Transform.Translation(zMoveVect))
    return [line, offsetDistance]

def ExtrudeElbow(items, extrusionH, tolerance):
    offset1 = OffsetCurve(items[0])
    curve1 = offset1[0]
    offsetDistance1 = offset1[1]
    
    curve2 = OffsetCurve(items[1])[0]
    st = curve1.PointAtStart
    end = curve1.PointAtEnd
    pts = [curve2.PointAtStart, curve2.PointAtEnd]
    
    closest0 = pts[rs.PointArrayClosestPoint(pts, st)]
    closest1 = pts[rs.PointArrayClosestPoint(pts, end)]
    closest = None
    arc = None
    
    if abs(rs.Distance(st, closest0)) < abs(rs.Distance(end, closest1)):
        closest = closest0   
        p = rs.CurveClosestPoint(curve1, st)
        arc = rs.AddArcPtTanPt(st, -rs.CurveTangent(curve1, p), closest)
    else:
        closest = closest1
        p = rs.CurveClosestPoint(curve1, end)
        arc = rs.AddArcPtTanPt(end, rs.CurveTangent(curve1, p), closest)
    arcObj = rs.coercecurve(arc)
    extrusion = Rhino.Geometry.Extrusion.CreateExtrusion(arcObj, Rhino.Geometry.Vector3d(0,0,extrusionH))
    brep = Rhino.Geometry.Brep.CreateFromOffsetFace(extrusion.ToBrep().Faces[0], offsetDistance1, tolerance, True, True)
    sc.doc.Objects.AddBrep(brep)


tolerance = sc.doc.ModelAbsoluteTolerance
cover = 0.05
ids = rs.GetObjects("select extrusion", rs.filter.extrusion+rs.filter.polysurface)
extrusionH = 3
rs.EnableRedraw(enable=False) 

#separate extrusions from elbows
extrusionsIds = []
elbows = []

for id in ids:
    obj = rs.coercegeometry(id)
    if str(obj.GetType()) == "Rhino.Geometry.Brep":
        elbows.append(obj)
    else:
        stPt = obj.PathStart
        endPt = obj.PathEnd
        extrusionsIds.append(id)
        #line_id = rs.AddLine(stPt, endPt)
        line = Rhino.Geometry.LineCurve(stPt, endPt)
        #print line
        dia = float(rs.GetUserText(id, 'Dimensions|Outside Diameter'))/1000
        offsetDistance = dia/2+cover
        zMoveVect = Rhino.Geometry.Vector3d(0,0,-offsetDistance)
        line.Transform(Rhino.Geometry.Transform.Translation(zMoveVect))
        # SET EXTRUSION HEIGHT
        #extrusion = Rhino.Geometry.Extrusion.Create(line, 10, True)
        extrusion = Rhino.Geometry.Extrusion.CreateExtrusion(line, Rhino.Geometry.Vector3d(0,0,extrusionH))
        brep = Rhino.Geometry.Brep.CreateFromOffsetFace(extrusion.ToBrep().Faces[0], offsetDistance, tolerance, True, True)
        sc.doc.Objects.AddBrep(brep)

counter = 0
for elbow in elbows:
    center = elbow.GetBoundingBox(Rhino.Geometry.Plane.WorldXY).Center
    #rs.AddTextDot('elbow: '+str(counter), center)
    connectedExtrusions = []
    for extrusionId in extrusionsIds:
        extrusion = rs.coercegeometry(extrusionId)
        intersCurves = Rhino.Geometry.Intersect.Intersection.BrepBrep(elbow, extrusion.ToBrep(), tolerance)[1]
        if len(intersCurves)>0:
            excenter = extrusion.ToBrep().GetBoundingBox(Rhino.Geometry.Plane.WorldXY).Center
            connectedExtrusions.append(extrusionId)
            #rs.AddTextDot(counter, excenter)
    ExtrudeElbow(connectedExtrusions,extrusionH,tolerance)
    counter += 1

RhinoToSAP.py

import rhinoscriptsyntax as rs
import csv

#destination csv file
filename = 'C:\\temp\\ElementiNuovi_Excel.csv'

#pick the lines you want to export
selectedLines = rs.GetObjects('Select lines', rs.filter.curve)

#array to store the points ids and coordinates
points = [['Joint', 'CoordSys', 'CoordType', 'XorR', 'Y', 'T', 'Z']] 

for l in selectedLines:
    attributes = rs.ObjectName(l).split(',')
    jI = attributes[0]
    jJ = attributes[1]
    frame = attributes[2]
    
    stDot = rs.AddTextDot(jI, rs.CurveStartPoint(l))
    rs.ObjectColor(stDot, (0,255,0))
    points.append([jI, 'GLOBAL', 'Cartesian', rs.CurveStartPoint(l).X,rs.CurveStartPoint(l).Y,'',rs.CurveStartPoint(l).Z])
    
    endDot = rs.AddTextDot(jJ, rs.CurveEndPoint(l))
    rs.ObjectColor(endDot, (0,0,255))
    points.append([jJ, 'GLOBAL', 'Cartesian', rs.CurveStartPoint(l).X,rs.CurveStartPoint(l).Y,'',rs.CurveStartPoint(l).Z])
    
    frameDot = rs.AddTextDot(frame, rs.CurveMidPoint(l))
    rs.ObjectColor(frameDot, (255,0,0))

with open(filename, "wb") as csvfile:
    csvwriter = csv.writer(csvfile,  delimiter=',')
    for point in points:
        csvwriter.writerow(point)
    print "Points sucessfully written to file"

SAPtoRhino.py

import rhinoscriptsyntax as rs
import csv

#update the view once all the geometry is computed (faster)
rs.EnableRedraw(False) 

#get all the points from SAP
points = []

#https://developer.rhino3d.com/guides/rhinopython/python-csv-file/
with open('C:\\temp\\Nodi_Excel.csv') as csvfile:
    reader = csv.reader(csvfile)
    reader.next() #skip first row TABLE:
    reader = csv.DictReader(csvfile)
    for row in reader:
        if (row['Joint']):
            points.append({'Joint': int(row['Joint']),'X':float(row['XorR']),'Y':float(row['Y']),'Z':float(row['Z'])})

"""
print points
for p in points:
    rs.AddPoint(p['X'],p['Y'],p['Z'])
"""


#get all the elementfrom SAP
elements = []

#https://developer.rhino3d.com/guides/rhinopython/python-csv-file/
with open('C:\\temp\\Elementi_Excel.csv') as csvfile:
    reader = csv.reader(csvfile)
    reader.next() #skip first row TABLE:
    reader = csv.DictReader(csvfile)
    for row in reader:
        if row['Frame']:
            elements.append({'Frame': int(row['Frame']),'JointI':int(row['JointI']),'JointJ':int(row['JointJ'])}) #assume all straight lines, no curves


for e in elements:
    #print e
    jI = e['JointI']
    jJ = e['JointJ']
    ptI = None
    ptJ = None
    # find the point that matches jI or jJ
    for p in points:
        currentJoint = p['Joint']
        if  currentJoint == jI:
            ptI = [p['X'],p['Y'],p['Z']]
        elif currentJoint == jJ:
            ptJ = [p['X'],p['Y'],p['Z']]
    #if both points exists, create a line. Write start point, end point and frame ID in the object name
    if ptI and ptJ:
        line = rs.AddLine(ptI, ptJ)
        rs.ObjectName(line, str(jI) + ',' + str(jJ) + "," + str(e['Frame']))

sightseeing.cs

using System;
using System.Collections;
using System.Collections.Generic;

using Rhino;
using Rhino.Geometry;

using Grasshopper;
using Grasshopper.Kernel;
using Grasshopper.Kernel.Data;
using Grasshopper.Kernel.Types;



/// <summary>
/// This class will be instantiated on demand by the Script component.
/// </summary>
public class Script_Instance : GH_ScriptInstance
{
#region Utility functions
  /// <summary>Print a String to the [Out] Parameter of the Script component.</summary>
  /// <param name="text">String to print.</param>
  private void Print(string text) { /* Implementation hidden. */ }
  /// <summary>Print a formatted String to the [Out] Parameter of the Script component.</summary>
  /// <param name="format">String format.</param>
  /// <param name="args">Formatting parameters.</param>
  private void Print(string format, params object[] args) { /* Implementation hidden. */ }
  /// <summary>Print useful information about an object instance to the [Out] Parameter of the Script component. </summary>
  /// <param name="obj">Object instance to parse.</param>
  private void Reflect(object obj) { /* Implementation hidden. */ }
  /// <summary>Print the signatures of all the overloads of a specific method to the [Out] Parameter of the Script component. </summary>
  /// <param name="obj">Object instance to parse.</param>
  private void Reflect(object obj, string method_name) { /* Implementation hidden. */ }
#endregion

#region Members
  /// <summary>Gets the current Rhino document.</summary>
  private readonly RhinoDoc RhinoDocument;
  /// <summary>Gets the Grasshopper document that owns this script.</summary>
  private readonly GH_Document GrasshopperDocument;
  /// <summary>Gets the Grasshopper script component that owns this script.</summary>
  private readonly IGH_Component Component;
  /// <summary>
  /// Gets the current iteration count. The first call to RunScript() is associated with Iteration==0.
  /// Any subsequent call within the same solution will increment the Iteration count.
  /// </summary>
  private readonly int Iteration;
#endregion

  /// <summary>
  /// This procedure contains the user code. Input parameters are provided as regular arguments,
  /// Output parameters as ref arguments. You don't have to assign output parameters,
  /// they will have a default value.
  /// </summary>
  private void RunScript(bool x, ref object A, ref object B, ref object C, ref object D, ref object E)
  {
    if (x){
      Rhino.Input.Custom.GetObject goMesh = new Rhino.Input.Custom.GetObject();
      goMesh.SetCommandPrompt("Select mesh");
      goMesh.GeometryFilter = Rhino.DocObjects.ObjectType.Mesh;
      goMesh.SubObjectSelect = false;
      goMesh.GroupSelect = false;
      goMesh.Get();

      Rhino.DocObjects.ObjRef objrefMesh = goMesh.Object(0);
      Rhino.Geometry.Mesh signalMesh = objrefMesh.Mesh();

      Rhino.Geometry.Point3d signalBBoxCenter = signalMesh.GetBoundingBox(true).Center;
      double radius = 0.75;
      Rhino.Geometry.Sphere sphere = new Rhino.Geometry.Sphere(signalBBoxCenter, radius);

      A = sphere;

      Rhino.Input.Custom.GetObject goCurve = new Rhino.Input.Custom.GetObject();
      goCurve.SetCommandPrompt("Select alignment");
      goCurve.GeometryFilter = Rhino.DocObjects.ObjectType.Curve;
      goCurve.SubObjectSelect = false;
      goCurve.GroupSelect = false;
      goCurve.Get();

      Rhino.DocObjects.ObjRef objrefCurve = goCurve.Object(0);
      Rhino.Geometry.Curve alignmentCurve = objrefCurve.Curve();

      B = alignmentCurve;


      double t = 0;
      alignmentCurve.ClosestPoint(signalBBoxCenter, out t);

      Rhino.Geometry.Point3d pointOnAlignment = alignmentCurve.PointAt(t);

      Rhino.Geometry.Curve splitAlignment = alignmentCurve.Split(t)[1];

      double speed = 80; //km/h
      int travelTime = 10; //seconds
      double distance = speed * 1000 / 3600; //distance per second
      double driverHeight = 2.7; //above rail

      Rhino.Geometry.Point3d[] viewPoints = splitAlignment.DivideEquidistant(distance);

      Point3d[] travelPoints = new Point3d[travelTime + 1];

      Array.Copy(viewPoints, travelPoints, travelTime + 1);

      C = travelPoints;

      Line[] sightLines = new Line[travelTime];
      Vector3d up = new Vector3d(0, 0, -1);

      Vector3d[] directions = new Vector3d[travelTime];
      string[] xmlViewpoints = new string[travelTime];

      for (int i = travelTime; i > 0; i--)
      {
        Point3d driverView = new Point3d(travelPoints[i].X, travelPoints[i].Y, travelPoints[i].Z + driverHeight);
        Point3d nextDriverView = new Point3d(travelPoints[i - 1].X, travelPoints[i - 1].Y, travelPoints[i - 1].Z + driverHeight);
        Line line = new Line(driverView, nextDriverView);
        sightLines.SetValue(line, i - 1);

        Vector3d direction = nextDriverView - driverView;
        //Vector3d direction = line.Direction;
        directions.SetValue(direction, i - 1);
        direction.Unitize();


        Vector3d xaxis = Vector3d.CrossProduct(up, direction);
        Vector3d yaxis = Vector3d.CrossProduct(direction, xaxis);

        xaxis.Unitize();
        yaxis.Unitize();


        List<double[]> m = new List<double[]>();

        double[] m0 = new double[3];
        double[] m1 = new double[3];
        double[] m2 = new double[3];

        m0.SetValue(xaxis.X, 0);
        m0.SetValue(yaxis.X, 1);
        m0.SetValue(direction.X, 2);


        m1.SetValue(xaxis.Y, 0);
        m1.SetValue(yaxis.Y, 1);
        m1.SetValue(direction.Y, 2);


        m2.SetValue(xaxis.Z, 0);
        m2.SetValue(yaxis.Z, 1);
        m2.SetValue(direction.Z, 2);

        m.Add(m0);
        m.Add(m1);
        m.Add(m2);

        double qw = Math.Sqrt(1 + m[0][0] + m[1][1] + m[2][2]) / 2;
        double qz = -(m[2][1] - m[1][2]) / (4 * qw);
        double qy = -(m[0][2] - m[2][0]) / ( 4 * qw);
        double qx = (m[1][0] - m[0][1]) / ( 4 * qw);

        string quaternion = String.Format("				<quaternion a=\"{0}\" b=\"{1}\" c=\"{2}\" d=\"{3}\"/>", qw, qx, qy, qz);

        string position = String.Format("<pos3f x=\"{0}\" y=\"{1}\" z=\"{2}\"/>", driverView.X, driverView.Y, driverView.Z);

        string view = String.Format("<view name=\"{0}\">\n" +
          "  <viewpoint tool = \"none\" render = \"preview\" lighting = \"headlight\" linear = \"22.2222891780\" angular = \"0.0771756700\"> \n" +
          "     <camera projection = \"persp\" near = \"1.0000000000\" far = \"10.0000000000\" aspect = \"1.0000000000\" height = \"0.7853980000\"> \n" +
          "          <position>\n" +
          "             {1}\n" +
          "           </position>\n" +
          "           <rotation>\n" +
          "             {2}\n" +
          "           </rotation>\n" +
          "     </camera>\n" +
          "  </viewpoint>\n" +
          "</view>\n", i * distance, position, quaternion);

        xmlViewpoints.SetValue(view, travelTime - i);

      }

      Guid g = Guid.NewGuid();

      string xmlAnimation = String.Format("<animation name=\"{0}\" guid=\"{1}\">\n" +
        "       {2}\n" +
        "</animation>", "Animation Name", g, string.Join("\n", xmlViewpoints));

      E = xmlAnimation;
      D = sightLines;
    }
    else
    {
      A = "Set true to run";
    }
  }

  // <Custom additional code> 

  // </Custom additional code> 
}

splitMeshByPlanes.py

import rhinoscriptsyntax as rs
import Rhino
import scriptcontext as sc
import math
from System.Collections.Generic import List
import System
import System.Collections.Generic.IEnumerable as IEnumerable

def splitMesh(m_id, rc):
    mesh = rs.coercemesh(m_id)
    print (m_id, mesh)
    centroid = Rhino.Geometry.AreaMassProperties.Compute(rc).Centroid
    normal = Rhino.Geometry.Surface.NormalAt(rc, 0.5,0.5)
    plane = Rhino.Geometry.Plane(centroid,normal)
    try:
        rs.DeleteObject(m_id)
    except Exception as e:
        print (e.message)
    return mesh.Split(plane)


mesh_id = rs.GetObject("Select a mesh to split", rs.filter.mesh)

srf_ids = rs.GetObjects("Select a cutter", rs.filter.surface)

splittedMeshes = []

bake = []

for id in srf_ids:
    print ('id',id)
    rhinoCutter = rs.coercesurface(id)
    try:
        for s in splitMesh(mesh_id, rhinoCutter):
            bake.append(s)
            splittedMeshes.append(s)
        
    except:
        for s in splittedMeshes:
            splittedMeshes = splitMesh(s, rhinoCutter)
            if len(splittedMeshes)>1:
                bake.remove(s)
                print ('s', s)
                for p in splittedMeshes:
                    bake.append(p)


for s in bake:
    sc.doc.Objects.AddMesh(s)
    
sc.doc.Views.Redraw()

#print (sm)

vectorDirToQuaternion.py

"""
reference:
https://forums.autodesk.com/t5/navisworks-forum/quaternion-function-for-viewpoint-export-import/td-p/2687701
"""

__author__ = "G"
__version__ = "2022.11.17"

import rhinoscriptsyntax as rs
import math


#xaxis = Vector.Normalized(Vector.Cross(up, direction))
xaxis = rs.VectorUnitize( rs.VectorCrossProduct(up,direction) )

#yaxis = Vector.Normalized(Vector.Cross(direction,xaxis))
yaxis = rs.VectorUnitize( rs.VectorCrossProduct(direction, xaxis) )


m = [[] for j in range(3)]
#row1
m[0].append(xaxis.X)
m[0].append(yaxis.X)
m[0].append(direction.X)

#row2
m[1].append(xaxis.Y)
m[1].append(yaxis.Y)
m[1].append(direction.Y)

#row3
m[2].append(xaxis.Z)
m[2].append(yaxis.Z)
m[2].append(direction.Z)


qw = math.sqrt(1 + m[0][0] + m[1][1] + m[2][2])/2
qz = -(m[2][1] - m[1][2]) / (4*qw)
qy = -(m[0][2] - m[2][0])/( 4 *qw)
qx = (m[1][0] - m[0][1])/( 4 *qw)


a = qw
b = qx
c = qy
d = qz

xml = '<quaternion a="{a}" b="{b}" c="{c}" d="{d}"/>'.format(a=qw,b=qx,c=qy,d=qz)

viewpoint.xml

<?xml version="1.0" encoding="UTF-8" ?>
    <exchange xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" xsi:noNamespaceSchemaLocation="http://download.autodesk.com/us/navisworks/schemas/nw-exchange-12.0.xsd" units="m" filename="" filepath="">
        <viewpoints>
            <animation name="_Default" guid="6d1c964c-2b71-4f7e-863a-08b9da8ac9b2">
                <view name="_Default_222">
    <viewpoint tool="none" render="preview" lighting="headlight" linear="22.2222891780" angular="0.0771756700">
      <camera projection="persp" near="1.0000000000" far="10.0000000000" aspect="1.0000000000" height="0.7853980000">
        <position>            
             <pos3f x="345028.699749" y="5812832.91015" z="130.199903282"/>
        </position>
        <rotation>
            				<quaternion a="0.698754957505" b="0.151779899758" c="0.148388872071" d="0.683143553023"/>
        </rotation>
      </camera>
    </viewpoint>
  </view>
<view name="_Default_200">
    <viewpoint tool="none" render="preview" lighting="headlight" linear="22.2222891780" angular="0.0771756700">
      <camera projection="persp" near="1.0000000000" far="10.0000000000" aspect="1.0000000000" height="0.7853980000">
        <position>            
             <pos3f x="345028.699749" y="5812832.91015" z="130.199903282"/>
        </position>
        <rotation>
            				<quaternion a="0.698754957505" b="0.151779899758" c="0.148388872071" d="0.683143553023"/>
        </rotation>
      </camera>
    </viewpoint>
  </view>
<view name="_Default_178">
    <viewpoint tool="none" render="preview" lighting="headlight" linear="22.2222891780" angular="0.0771756700">
      <camera projection="persp" near="1.0000000000" far="10.0000000000" aspect="1.0000000000" height="0.7853980000">
        <position>            
             <pos3f x="345028.699749" y="5812832.91015" z="130.199903282"/>
        </position>
        <rotation>
            				<quaternion a="0.698754957505" b="0.151779899758" c="0.148388872071" d="0.683143553023"/>
        </rotation>
      </camera>
    </viewpoint>
  </view>
<view name="_Default_156">
    <viewpoint tool="none" render="preview" lighting="headlight" linear="22.2222891780" angular="0.0771756700">
      <camera projection="persp" near="1.0000000000" far="10.0000000000" aspect="1.0000000000" height="0.7853980000">
        <position>            
             <pos3f x="345028.699749" y="5812832.91015" z="130.199903282"/>
        </position>
        <rotation>
            				<quaternion a="0.698754957505" b="0.151779899758" c="0.148388872071" d="0.683143553023"/>
        </rotation>
      </camera>
    </viewpoint>
  </view>
<view name="_Default_133">
    <viewpoint tool="none" render="preview" lighting="headlight" linear="22.2222891780" angular="0.0771756700">
      <camera projection="persp" near="1.0000000000" far="10.0000000000" aspect="1.0000000000" height="0.7853980000">
        <position>            
             <pos3f x="345028.699749" y="5812832.91015" z="130.199903282"/>
        </position>
        <rotation>
            				<quaternion a="0.698754957505" b="0.151779899758" c="0.148388872071" d="0.683143553023"/>
        </rotation>
      </camera>
    </viewpoint>
  </view>
<view name="_Default_111">
    <viewpoint tool="none" render="preview" lighting="headlight" linear="22.2222891780" angular="0.0771756700">
      <camera projection="persp" near="1.0000000000" far="10.0000000000" aspect="1.0000000000" height="0.7853980000">
        <position>            
             <pos3f x="345028.699749" y="5812832.91015" z="130.199903282"/>
        </position>
        <rotation>
            				<quaternion a="0.698754957505" b="0.151779899758" c="0.148388872071" d="0.683143553023"/>
        </rotation>
      </camera>
    </viewpoint>
  </view>
<view name="_Default_089">
    <viewpoint tool="none" render="preview" lighting="headlight" linear="22.2222891780" angular="0.0771756700">
      <camera projection="persp" near="1.0000000000" far="10.0000000000" aspect="1.0000000000" height="0.7853980000">
        <position>            
             <pos3f x="345019.483332" y="5812853.12416" z="130.701917806"/>
        </position>
        <rotation>
            				<quaternion a="0.699570652419" b="0.151976964876" c="0.148225026137" d="0.682299967791"/>
        </rotation>
      </camera>
    </viewpoint>
  </view>
<view name="_Default_067">
    <viewpoint tool="none" render="preview" lighting="headlight" linear="22.2222891780" angular="0.0771756700">
      <camera projection="persp" near="1.0000000000" far="10.0000000000" aspect="1.0000000000" height="0.7853980000">
        <position>            
             <pos3f x="345010.266098" y="5812873.33706" z="131.257299826"/>
        </position>
        <rotation>
            				<quaternion a="0.697570383484" b="0.15202891007" c="0.149102471004" d="0.684142692521"/>
        </rotation>
      </camera>
    </viewpoint>
  </view>
<view name="_Default_044">
    <viewpoint tool="none" render="preview" lighting="headlight" linear="22.2222891780" angular="0.0771756700">
      <camera projection="persp" near="1.0000000000" far="10.0000000000" aspect="1.0000000000" height="0.7853980000">
        <position>            
             <pos3f x="345001.022401" y="5812893.53669" z="131.689102470"/>
        </position>
        <rotation>
            				<quaternion a="0.692313759095" b="0.154616986866" c="0.153629123869" d="0.687890499018"/>
        </rotation>
      </camera>
    </viewpoint>
  </view>
<view name="_Default_022">
    <viewpoint tool="none" render="preview" lighting="headlight" linear="22.2222891780" angular="0.0771756700">
      <camera projection="persp" near="1.0000000000" far="10.0000000000" aspect="1.0000000000" height="0.7853980000">
        <position>            
             <pos3f x="344991.573770" y="5812913.63521" z="131.831452336"/>
        </position>
        <rotation>
            				<quaternion a="0.686402474227" b="0.159744236067" c="0.160812799951" d="0.690993969431"/>
        </rotation>
      </camera>
    </viewpoint>
  </view>

            </animation>

        </viewpoints>
    </exchange>

vscodeSetup.cs

# Follow this: https://discourse.mcneel.com/t/howto-run-c-script-in-rhino-using-vs-code/87642 
# then update the .csproj dll file paths 
# and target framework see: https://docs.microsoft.com/en-us/dotnet/standard/frameworks

<Project Sdk="Microsoft.NET.Sdk">
  <PropertyGroup>
    <TargetFramework>netstandard2.1</TargetFramework>
  </PropertyGroup>
  <PropertyGroup Condition="'$(Configuration)|$(Platform)'=='Debug|AnyCPU'">
  </PropertyGroup>
  <ItemGroup>
    <Reference Include="Grasshopper">
      <HintPath>C:\Program Files\Rhino 7\Plug-ins\Grasshopper\Grasshopper.dll</HintPath>
      <Private>False</Private>
    </Reference>
    <Reference Include="RhinoCommon">
      <HintPath>C:\Program Files\Rhino 7\System\RhinoCommon.dll</HintPath>
      <Private>False</Private>
    </Reference>
  </ItemGroup>
</Project>