Hello All,
Inspired by Bruce's notification on the demo, I perpared a quick Vpython version for the windmills part of it.
I used PyAudio to capture mic input and used it to control head axis, and rotation speed.
My laptop has a mono built-in mic, so I can only see the blades rotation. If you have a stereo built-in mic, I believe you should see both.
I did not elaborate on using a web-cam's mic for stereo input since I could not see any support for it by PyAudio, and I did not want to go into details of
wave library and WinMM api.
The video part is more tricky of course. It needs a lot more work on image processing, and also some new capabilities on VPython side, to paint the background with a live video image, rather than a plane color, as it is today, which is well beyond my scope.
Start the application, and just talk to the screen ( the mic)...
Kadir
# ecoMILLs in VPython
from visual import *
import pyaudio
import sys
from array import array
###### Set the Scene
scene.autocenter = 0
scene.center = (0,10,0)
scene.range = 20
scene.autoscale = 0
scene.width = 1024
scene.height = 768
scene.title = "ecoMILLs"
field = box(size=(50,1,50), color=(0,1,0))
###### Prepare Audio
threshold = 200
maxValue = 0
chunk = 2048
FORMAT = pyaudio.paInt16
CHANNELS = 2
RATE = 44100
RECORD_SECONDS = 1
p = pyaudio.PyAudio()
stream = p.open(format = FORMAT,
channels = CHANNELS,
rate = RATE,
input = True,
output = True,
frames_per_buffer = chunk)
######### Setup a windmill
class WindMill():
def __init__(self, offset=(0,0,0), color1=(1,0,0), color2=(1,1,0),
color3=(0,0,1)):
self.wndml = frame(pos=offset)
self.head = frame(frame=self.wndml)
self.wings = frame(frame=self.head)
self.pole = cone(frame=self.wndml, radius=0.5, axis=(0,10,0))
self.dyno1 = cone(frame=self.head, pos=self.pole.axis-(0,0,1), radius=0.5,
axis=(0,0,-1))
self.dyno2 = cylinder(frame=self.head, pos=self.pole.axis-(0,0,1),
radius=0.5, axis=(0,0,1))
self.bld1 = pyramid(frame=self.wings, pos=self.pole.axis+(0,2.5,0),
size=(0.5, 5.0, 0.05), color=color1)
self.bld2 = pyramid(frame=self.wings, pos=self.pole.axis+(0,2.5,0),
size=(0.5, 5.0, 0.05), color=color2)
self.bld2.rotate(axis=(0,0,1), angle=pi*2/3, origin=self.pole.axis)
self.bld3 = pyramid(frame=self.wings, pos=self.pole.axis+(0,2.5,0),
size=(0.5, 5.0, 0.05), color=color3)
self.bld3.rotate(axis=(0,0,1), angle=-pi*2/3, origin=self.pole.axis)
self.shaft = cylinder(frame=self.wings, pos=self.pole.axis, radius=0.1,
axis=(0,0,0.1), color=(0.4,0.4,0.4))
axis = (1,0,0)
self.rotate_wings(pi/5+random.uniform(-pi/4,+pi/4))
def rotate_head(self, axis):
self.head.axis = axis
rotang = pi/64.0
def rotate_wings(self, rotang):
self.wings.rotate(axis=self.dyno1.axis, origin=self.pole.axis,
angle=rotang)
##### Create WindMills
mills = []
mill1 = WindMill()
mill2 = WindMill((10,0,3), color1=(0,1,0), color2=(0,1,1), color3=(1,0,1))
mill3 = WindMill((-10,0,3), color1=(1,1,1), color2=(1,1,0), color3=(1,0,0))
mills.append(mill1)
mills.append(mill2)
mills.append(mill3)
left = 1
right = 1
headax = (1,0,0)
rotang = pi/64.
while left >0 or right > 0:
for i in range(0, 44100 / chunk * RECORD_SECONDS):
data = stream.read(chunk)
asInts = array('h', data)
for j in range(0,len(asInts),2):
left += asInts[j]
right += asInts[j+1]
if left > threshold or right > threshold:
headax = ((1,0,1-(right/left)))
rotang = pi/64 + (pi/64)*(left+right)/10000
for mill in mills:
mill.rotate_head(headax)
mill.rotate_wings(rotang)
left = 1
right = 1
stream.stop_stream()
stream.close()
p.terminate()
________________________________
From: Bruce Sherwood <Bru...@nc...>
To: VPython <vis...@li...>
Sent: Tuesday, June 30, 2009 8:03:28 AM
Subject: [Visualpython-users] Very cool 3D demo
Fans of 3D, you might find this as intriguing as I did:
http://ge.ecomagination.com/smartgrid/#/augmented_reality
Print the "Solar Panel Marker"
page, then follow the rest of the instructions on the page.
You need to have a webcam on your computer; the webcam views the
special page that you print and follows the page as you move it around,
thereby giving you navigation in the 3D scene.
Bruce Sherwood
|
