{"id":1478,"date":"2012-12-11T14:44:04","date_gmt":"2012-12-11T14:44:04","guid":{"rendered":"http:\/\/www.new-territories.com\/blog\/ncertainties\/col12\/?page_id=1478"},"modified":"2012-12-19T01:01:06","modified_gmt":"2012-12-19T01:01:06","slug":"code-5-add-title","status":"publish","type":"page","link":"https:\/\/www.new-territories.com\/blog\/ncertainties\/col12\/code\/code-5-add-title\/","title":{"rendered":"EarthClouds"},"content":{"rendered":"

MAIN CODE\u00a0<\/strong><\/p>\n

import rhinoscriptsyntax as rs<\/p>\n

import math<\/p>\n

class waveSource():<\/p>\n

def __init__(self, GUID, YVALUES):<\/p>\n

self.id<\/a> = GUID<\/p>\n

self.waveValues = YVALUES<\/p>\n

self.alpha = 0<\/p>\n

def explicitWaveVec(self,pt):<\/p>\n

if rs.IsPoint(self.id<\/a>):<\/p>\n

self.pos = rs.PointCoordinates(self.id<\/a>)<\/p>\n

else:<\/p>\n

t = rs.CurveClosestPoint(self.id<\/a>,pt)<\/p>\n

self.pos = rs.EvaluateCurve(self.id<\/a>,t)<\/p>\n

dist = rs.Distance(self.pos, pt)<\/p>\n

sine = self.waveValues[int(dist)]<\/p>\n

vec = rs.VectorCreate(pt,self.pos)<\/p>\n

vec = rs.VectorUnitize(vec)<\/p>\n

vec = rs.VectorScale(vec, sine)<\/p>\n

return vec<\/p>\n

def updateExplicit(self):<\/p>\n

#temp = self.waveValues[0]<\/p>\n

self.waveValues.pop(0)<\/p>\n

#print len(self.waveValues)<\/p>\n

class medium():<\/p>\n

def __init__(self,GUID, FRICTION):<\/p>\n

self.id<\/a> = GUID<\/p>\n

self.friction = FRICTION<\/p>\n

self.pts = rs.PointCloudPoints(self.id<\/a>)<\/p>\n

def update(self,waveSources):<\/p>\n

newPts = []<\/p>\n

for pt in self.pts:<\/p>\n

sumVec = [0,0,0]<\/p>\n

for myWaveSource in waveSources:<\/p>\n

sumVec = rs.PointAdd(sumVec,myWaveSource.explicitWaveVec(pt))<\/p>\n

sumVec = rs.VectorScale(sumVec,1-self.friction)<\/p>\n

newPts.append(rs.PointAdd(pt,sumVec))<\/p>\n

#rs.DeleteObject(self.id<\/a>)<\/p>\n

self.id<\/a> = rs.AddPointCloud(newPts)<\/p>\n

rs.AddLayer(“MyLayer####”)<\/p>\n

rs.ObjectLayer(self.id<\/a>,”MyLayer####”)<\/p>\n

self.pts = newPts<\/p>\n

 <\/p>\n

def Main():<\/p>\n

waveCrv = rs.GetObject(“select the curve for the curve in x-y plane”, rs.filter.curve)<\/p>\n

yValues = []<\/p>\n

rs.EnableRedraw(False)<\/p>\n

for i in range(500):<\/p>\n

temp = rs.AddLine([i,-100,0],[i,100,0])<\/p>\n

ccx = rs.CurveCurveIntersection(waveCrv, temp)<\/p>\n

if ccx:<\/p>\n

yValues.append(ccx[0][1][1])<\/p>\n

else:<\/p>\n

yValues.append(0)<\/p>\n

rs.DeleteObject(temp)<\/p>\n

sources = rs.GetObjects(“select the sources”, rs.filter.point+rs.filter.curve)<\/p>\n

ptClouds = rs.GetObjects(“select the mediums”, rs.filter.pointcloud)<\/p>\n

myWaveSources = []<\/p>\n

for source in sources:<\/p>\n

myWaveSources.append(waveSource(source, yValues))<\/p>\n

myMediums = []<\/p>\n

for ptCloud in ptClouds:<\/p>\n

name = rs.ObjectName(ptCloud)<\/p>\n

if name:<\/p>\n

myMediums.append( medium(ptCloud, float(name)) )<\/p>\n

else:<\/p>\n

myMediums.append( medium(ptCloud, 0) )<\/p>\n

for i in range(20):<\/p>\n

for myMedium in myMediums:<\/p>\n

myMedium.update(myWaveSources)<\/p>\n

for myWaveSource in myWaveSources:<\/p>\n

myWaveSource.updateExplicit()<\/p>\n

 <\/p>\n

 <\/p>\n

Main()<\/p>\n

 <\/p>\n

 <\/p>\n

 <\/p>\n

COLOR CODE<\/strong><\/p>\n

 <\/p>\n

import rhinoscriptsyntax as rs<\/p>\n

import random<\/p>\n

 <\/p>\n

class robot():<\/p>\n

def __init__(self, PATH, RADIOUS, TOTALITERATIONS):<\/p>\n

self.path = PATH<\/p>\n

self.radious = RADIOUS<\/p>\n

self.steps = TOTALITERATIONS<\/p>\n

self.pathPts = rs.DivideCurve(self.path, self.steps+1)<\/p>\n

self.pos = self.pathPts[0]<\/p>\n

self.folder = rs.BrowseForFolder()<\/p>\n

self.averageDistances = []<\/p>\n

def updateCamera(self, currentI):<\/p>\n

view = rs.CurrentView()<\/p>\n

rs.ViewCameraTarget( view, self.pathPts[currentI], self.pathPts[currentI+1] )<\/p>\n

self.pos = self.pathPts[currentI]<\/p>\n

def saveFrame(self, currentI):<\/p>\n

if currentI < 100 and currentI > 10:<\/p>\n

result = rs.CreatePreviewImage(self.folder + “\\\\robotView0<\/a>” + str(currentI) + “.jpg”)<\/p>\n

elif currentI<10:<\/p>\n

result = rs.CreatePreviewImage(self.folder + “\\\\robotView00<\/a>” + str(currentI) + “.jpg”)<\/p>\n

else:<\/p>\n

result = rs.CreatePreviewImage(self.folder + “\\\\robotView<\/a>” + str(currentI) + “.jpg”)<\/p>\n

def look(self, myEarthclouds):<\/p>\n

sumDistances = 0<\/p>\n

count = 0<\/p>\n

for myEarthcloud in myEarthclouds:<\/p>\n

dist = rs.Distance(self.pos, myEarthcloud.pos)<\/p>\n

if dist <self.radious:<\/p>\n

#rs.ShowObject(myEarthcloud.id)<\/p>\n

currentDistance= myEarthcloud.colorByDistance2Robot(self)<\/p>\n

count = count +1<\/p>\n

else:<\/p>\n

rs.HideObject(myEarthcloud.id)<\/p>\n

rs.HideObjects(myEarthcloud.ptids)<\/p>\n

currentDistance= 0<\/p>\n

sumDistances = sumDistances + currentDistance<\/p>\n

if count > 0 :<\/p>\n

avDistances = sumDistances\/count<\/p>\n

else:<\/p>\n

avDistances = 0<\/p>\n

self.averageDistances.append(avDistances)<\/p>\n

 <\/p>\n

class earthcloud():<\/p>\n

def __init__(self, PTCLOUDGUID, POINTNEIGHBORHOODRADIOUS):<\/p>\n

self.id<\/a> = PTCLOUDGUID<\/p>\n

self.radious = POINTNEIGHBORHOODRADIOUS<\/p>\n

self.friction = random.random()<\/p>\n

self.pos = [0,0,0]<\/p>\n

self.pts = rs.PointCloudPoints(self.id<\/a>)<\/p>\n

for pt in self.pts:<\/p>\n

self.pos = rs.PointAdd(self.pos, pt)<\/p>\n

self.pos = [self.pos[0]\/len(self.pts),self.pos[1]\/len(self.pts),self.pos[2]\/len(self.pts)]<\/p>\n

rs.HideObject(self.id<\/a>)<\/p>\n

self.ptids = []<\/p>\n

for pt in self.pts:<\/p>\n

self.ptids.append(rs.AddPointCloud([pt]))<\/p>\n

rs.HideObject(self.ptids[-1])<\/p>\n

def colorByDistance2Robot(self, robot):<\/p>\n

min = 5000000000000<\/p>\n

max = 0<\/p>\n

allDist = []<\/p>\n

sumDistances = 0<\/p>\n

count = 0<\/p>\n

for i in range(len(self.pts)):<\/p>\n

pt = self.pts[i]<\/p>\n

dist = rs.Distance(pt, robot.pos)<\/p>\n

if dist<min : min = dist<\/p>\n

if dist>max : max = dist<\/p>\n

allDist.append(dist)<\/p>\n

if dist<robot.radious:<\/p>\n

sumDistances = sumDistances + dist<\/p>\n

count +=1<\/p>\n

if count>0:<\/p>\n

sumDistances = sumDistances\/count<\/p>\n

for i in range(len(self.pts)):<\/p>\n

#pt = self.pts[i]<\/p>\n

ptid = self.ptids[i]<\/p>\n

normalizedDistance = 1-(allDist[i]-min)\/(max-min)<\/p>\n

rs.ShowObject(self.ptids[i])<\/p>\n

rs.ObjectColor(self.ptids[i], [255*normalizedDistance,0,255-255*normalizedDistance])<\/p>\n

return rs.Distance(robot.pos, self.pos)<\/p>\n

 <\/p>\n

 <\/p>\n

def Main():<\/p>\n

crv = rs.GetObject(“select the path for the robot”, rs.filter.curve)<\/p>\n

radious = rs.GetReal(“plase type the radious of the robot”, 5)<\/p>\n

frameCount = rs.GetInteger(“please type the amount of frames”, 200)<\/p>\n

myRobot = robot(crv, radious, frameCount)<\/p>\n

ptClouds = rs.GetObjects(“select a bunch of point clouds”, rs.filter.pointcloud)<\/p>\n

myEarthclouds = []<\/p>\n

for ptCloud in ptClouds:<\/p>\n

myEarthclouds.append(earthcloud(ptCloud, 0.7))<\/p>\n

rs.HideObject(myRobot.path)<\/p>\n

for i in range(frameCount):<\/p>\n

myRobot.look(myEarthclouds)<\/p>\n

myRobot.updateCamera(i)<\/p>\n

myRobot.saveFrame(i)<\/p>\n

for myEarthcloud in myEarthclouds:<\/p>\n

rs.ShowObject(myEarthcloud.id)<\/p>\n

rs.DeleteObjects(myEarthcloud.ptids)<\/p>\n

 <\/p>\n

pts4Distances = []<\/p>\n

for i in range(len(myRobot.averageDistances)):<\/p>\n

distanceValue = myRobot.averageDistances[i]<\/p>\n

pts4Distances.append([i,distanceValue,0])<\/p>\n

p = rs.AddPolyline(pts4Distances)<\/p>\n

rs.ObjectColor(p,[255,0,0])<\/p>\n

rs.ShowObject(myRobot.path)<\/p>\n

Main()<\/p>\n","protected":false},"excerpt":{"rendered":"

MAIN CODE\u00a0 import rhinoscriptsyntax as rs import math class waveSource(): def __init__(self, GUID, YVALUES): self.id = GUID self.waveValues = YVALUES self.alpha = 0 def explicitWaveVec(self,pt): if rs.IsPoint(self.id): self.pos = rs.PointCoordinates(self.id) else: t = rs.CurveClosestPoint(self.id,pt) self.pos = rs.EvaluateCurve(self.id,t) dist = rs.Distance(self.pos, pt) sine = self.waveValues[int(dist)] vec = rs.VectorCreate(pt,self.pos) vec = rs.VectorUnitize(vec) vec = rs.VectorScale(vec, sine) return […]<\/p>\n","protected":false},"author":2,"featured_media":0,"parent":9,"menu_order":0,"comment_status":"closed","ping_status":"open","template":"","meta":{"footnotes":""},"class_list":["post-1478","page","type-page","status-publish","hentry"],"_links":{"self":[{"href":"https:\/\/www.new-territories.com\/blog\/ncertainties\/col12\/wp-json\/wp\/v2\/pages\/1478","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/www.new-territories.com\/blog\/ncertainties\/col12\/wp-json\/wp\/v2\/pages"}],"about":[{"href":"https:\/\/www.new-territories.com\/blog\/ncertainties\/col12\/wp-json\/wp\/v2\/types\/page"}],"author":[{"embeddable":true,"href":"https:\/\/www.new-territories.com\/blog\/ncertainties\/col12\/wp-json\/wp\/v2\/users\/2"}],"replies":[{"embeddable":true,"href":"https:\/\/www.new-territories.com\/blog\/ncertainties\/col12\/wp-json\/wp\/v2\/comments?post=1478"}],"version-history":[{"count":7,"href":"https:\/\/www.new-territories.com\/blog\/ncertainties\/col12\/wp-json\/wp\/v2\/pages\/1478\/revisions"}],"predecessor-version":[{"id":1734,"href":"https:\/\/www.new-territories.com\/blog\/ncertainties\/col12\/wp-json\/wp\/v2\/pages\/1478\/revisions\/1734"}],"up":[{"embeddable":true,"href":"https:\/\/www.new-territories.com\/blog\/ncertainties\/col12\/wp-json\/wp\/v2\/pages\/9"}],"wp:attachment":[{"href":"https:\/\/www.new-territories.com\/blog\/ncertainties\/col12\/wp-json\/wp\/v2\/media?parent=1478"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}