############################## import rhinoscriptsyntax as rs import random, time, sys ### sys.path.append("P:/") ### import DM_lib as dm ### reload(dm) ############################## rs.UnitSystem(4) # km = 5, meters = 4, cm = 3 etc rs.ShowGrid(None, 0) # grid > 0 = off rs.ShowGridAxes(None, 1) # y/y/z axen display > 0/1 = off/on rs.ViewDisplayMode(rs.CurrentView(), "wireframe") rs.Command("cplane w t enter", 0) # cPlane World Top dm.PointRadius(displayModeX=0, rad=3, styl=3) dm.printDisplay(0) # nomen est omen rs.EnableRedraw(0) ### MAC >> uncomment ! ###_________________________________________# ### basic settings for grid to fit UN_slab # ### # ! no need 2 change ! floors = H = dm.H = 40 # default=40 / incl roof slabs = L = dm.L = 11 # default=11 depth = D = dm.D = 4 # default= 4 / division in building_depth floorHeight = fH = dm.fH = 4.0 # default= 4.0 / 4.0*(H-1) = 156 meters # ############# # get gridCoords L*D*H = 1760 UnoGridCoords = dm.getUnoGridCoords() # get from DM_lib as dm #################___________________________# dm.newEmptyLayer("UNO::setUp::points", [150,150,200]) dm.newEmptyLayer("UNO::setUp::geo", [150, 20,150]) dm.newEmptyLayer("UNO::setUp::crv", [200, 100, 100]) dm.newEmptyLayer("UNO::setUp", [120,120,140]) ################## depthVec = dVec = rs.VectorUnitize(rs.VectorSubtract( dm.getUnoCoord(0, 0, 0), dm.getUnoCoord(0, 2, 0) )) lengthVec = lVec = rs.VectorUnitize(rs.VectorSubtract( dm.getUnoCoord(0, 0, 0), dm.getUnoCoord(1, 0, 0) )) ################## if 1: ### SETUP >> dont' exec @ homework ! if 1: rs.AddLine( dm.getUnoCoord(0, 0, 0), dm.getUnoCoord(1, 0, 0) ) rs.AddLine( dm.getUnoCoord(0, 0, 0), dm.getUnoCoord(0, 1, 0) ) rs.AddLine( dm.getUnoCoord(0, 0, 0), dm.getUnoCoord(0, 0, 1) ) rs.ObjectColor(rs.AllObjects()[1], [0, 222, 0] ) rs.ObjectColor(rs.AllObjects()[2], [222, 0, 0] ) rs.ObjectColor(rs.AllObjects()[0], [0, 0, 222] ) rs.ObjectPrintWidth( rs.AllObjects()[0:3], 1.0 ) p0 = dm.getUnoCoord( 0, 0, 0) p1 = dm.getUnoCoord( 10, 0, 0) p2 = dm.getUnoCoord( 10, 3, 0) p3 = dm.getUnoCoord( 0, 3, 0) unoBaseCoords = [p0, p1, p2, p3, p0] unoBaseCrv = rs.AddCurve( unoBaseCoords, 1) unoTopCoords = [dm.getUnoCoord( 0, 0, 39), dm.getUnoCoord( 10, 0, 39), dm.getUnoCoord( 10, 3, 39), dm.getUnoCoord( 0, 3, 39), dm.getUnoCoord( 0, 0, 39), ] unoTopCrv = rs.AddCurve( unoTopCoords, 1) rs.ZoomExtents() if 0: rs.CurrentLayer("UNO::setUp::points") #rs.AddPoints( UnoGridCoords ) for i in range(L): rs.ObjectColor( rs.AddTextDot( str(i), dm.getUnoCoord(i, 0, 0) ), [0, 128, 128] ) for i in range(D): rs.ObjectColor( rs.AddTextDot( str(i), dm.getUnoCoord(0, i, 0) ), [0, 128, 128] ) for i in range(H): rs.ObjectColor( rs.AddTextDot( str(i), dm.getUnoCoord(0, 0, i) ), [0, 128, 128] ) #### if 0: rs.CurrentLayer("UNO::setUp::crv") dm.UN_slab(showAll=1) for i in range(9): p0 = dm.getUnoCoord( random.randint(0,L+1), random.randint(0,D+1), random.randint(0,H+1) ) p1 = dm.getUnoCoord( random.randint(0,L+1), random.randint(0,D+1), random.randint(0,H+1) ) rs.AddCurve( [p0, p1] ) #dm.textDots( [p0, p1] ) rs.ZoomExtents() if 0: coords = UnoGridCoords[:] random.shuffle( coords ) rs.AddCurve( coords[0:], random.randint(1, 2) ) coords = [] dirL = 1 Lval = 0 for L in range(40): #print L d=0 if L>0 and L%10 == 0: dirL *= -1 print "11er", L, dirL d = L%1 #print L, d cor = dm.getUnoCoord(Lval, d, L*1) coords.append( cor ) Lval += dirL crv = rs.AddCurve( coords, 2) dom = rs.CurveDomain( crv )[1] print "dom", dom paraX = dom*0.5 print "paraX", paraX cor = rs.EvaluateCurve( crv, paraX ) #rs.AddPoint( cor ) rad = 4.44 coords = rs.DivideCurve( crv, 45, create_points=1) allCoords = [] anzCirc = 33 for cor in coords: pass paraX = rs.CurveClosestPoint(crv, cor) normEbene = rs.CurvePerpFrame(crv, paraX) circ = rs.AddCircle( normEbene, rad) coordsX = rs.DivideCurve( circ, anzCirc, create_points=1) allCoords.append(coordsX) i = 10 for i in range(anzCirc): coords = [] for list in allCoords: cor = list[i] #rs.AddPoint( cor ) coords.append( cor ) rs.AddCurve( coords, 3 ) rs.AddInterpCurve( coords, 3) def gC (l=1, d=2, h=3): cor = dm.getUnoCoord(l, d, h) return cor #second 'slide' rs.CurrentLayer("UNO::setUp::crv") p0 = gC (l=10, d=0, h=39) p00 = gC (10, 3, 39) rs.AddPoint( p0 ) p1 = rs.VectorAdd(p0, rs.VectorScale( dVec, 100.0 ) ) p2 = gC (0, 3, 39) p3 = gC (0, 0, 39) p4 = rs.VectorAdd(p3, rs.VectorScale( dVec, -100.0 ) ) p5 = gC (10, 0, 35) p6 = gC (10, 3, 35) p7 = rs.VectorAdd(p6, rs.VectorScale( dVec, 100.0 ) ) p8 = gC (0, 3, 35) p9 = gC (0, 0, 35) p10 = rs.VectorAdd(p9, rs.VectorScale( dVec, -100.0 ) ) p11 = gC (5, 0, 30) p12 = gC (5, 3, 30) p13 = rs.VectorAdd(p12, rs.VectorScale( dVec, 100.0 ) ) p14 = gC (0, 3,30) p15 = gC (0, 0, 30) p16 = rs.VectorAdd(p15, rs.VectorScale( dVec, -100.0 ) ) p17 = gC (10, 0, 25) p18 = gC (10, 3, 25) p19 = rs.VectorAdd(p18, rs.VectorScale( dVec, 100.0 ) ) p20 = gC (0, 3, 25) p21 = gC (0, 0, 25) p22 = rs.VectorAdd(p21, rs.VectorScale( dVec, -100.0 ) ) p23 = gC (5, 0, 20) p24 = gC (5, 3, 20) p25 = rs.VectorAdd(p24, rs.VectorScale( dVec, 100.0 ) ) p26 = gC (0, 3, 20) p27 = gC (0, 0, 20) p28 = rs.VectorAdd(p27, rs.VectorScale( dVec, -100.0 ) ) p29 = gC (10, 0, 15) p30 = gC (10, 3, 15) p31 = rs.VectorAdd(p30, rs.VectorScale( dVec, 100.0 ) ) p32 = gC (0, 3, 15) p33 = gC (0, 0, 15) p34 = rs.VectorAdd(p33, rs.VectorScale( dVec, -100.0 ) ) p35 = gC (5, 0, 10) p36 = gC (5, 3, 10) p37 = rs.VectorAdd(p36, rs.VectorScale( dVec, 100.0 ) ) p38 = gC (0, 3, 10) p39 = gC (0, 0, 10) p40 = rs.VectorAdd(p39, rs.VectorScale( dVec, -100.0 ) ) p41 = gC (10, 0, 5) p42 = gC (10, 3, 5) p43 = rs.VectorAdd(p42, rs.VectorScale( dVec, 100.0 ) ) p44 = gC (0, 3, 5) p45 = gC (0, 0, 5) p46 = rs.VectorAdd(p45, rs.VectorScale( dVec, -100.0 ) ) p47 = gC (5, 0, 0) p48 = gC (5, 3, 0) p49 = rs.VectorAdd(p48, rs.VectorScale( dVec, 100.0 ) ) p50 = gC (0, 3, 0) p51 = gC (0, 0, 0) coords = [p0, p00, p1, p2, p3, p4, p5, p6, p7, p8, p9, p10, p11, p12, p13, p14, p15, p16, p17, p18, p19, p20, p21, p22, p23, p24, p25, p26, p27, p28, p29, p30, p31, p32, p33, p34, p35, p36, p37, p38, p39, p40, p41, p42, p43, p44, p45, p46, p47, p48, p49, p50, p51] crvDiv = 75 crv = rs.AddCurve( coords, 2 ) coords = rs.DivideCurve( crv, crvDiv, create_points=0 ) #rs.AddPoints( coords ) for cor in coords: pass dom = rs.CurveDomain( crv ) print "dom =", dom dom = rs.CurveDomain( crv )[1] print "dom =", dom circDiv = 90 allCoordsLists = [] for cor in coords: paraX = rs.CurveClosestPoint( crv, cor ) planeX = rs.CurvePerpFrame( crv, paraX ) circ = rs.AddCircle( planeX, 3) circCoords = rs.DivideCurve( circ, circDiv, 0) rs.DeleteObject( circ ) allCoordsLists.append( circCoords ) i = 0 for i in range(circDiv): coords = [] for list in allCoordsLists: pX = list[i] coords.append( pX ) rs.ObjectColor( rs.AddCurve( coords, 2 ), [0, 0, 128] ) # Anchor selected points of curves to the base of the building base_height = 0 # Define the z-value for the base of the building (default is 0) step = 5 # Define step size to reduce the number of anchor points for i in range(0, circDiv, step): # Iterate with a step to skip some points for coord_list in allCoordsLists: pX = coord_list[i] # Get the point on the curve # Add a line from each selected point to the building base base_point = rs.AddPoint(pX[0], pX[1], base_height) # Base point at z=base_height rs.AddLine(pX, base_point) # Line from curve point to building base #____________________________:here you end ################################### ################################### #dm.zA( proz=0.95 ) rs.EnableRedraw(1) ### 4 the MACs #dm.printDisplay(1) dm.newEmptyLayer("Default")