Tag: Improvements

Making the world a better place, one piece at a time

Storm Door Brace: Now, With Inserts!

The mid-1950s wood doors on our house have wood storm doors with interchangeable wood-framed glass and screen panels. Twice a year, the diligent homeowner will swap the panels to match the season; during the last 60+ years, the glass panels remain undropped.

The back door has a diagonal tension brace to hold the door in shape; the door may be slightly distorted or the frame slightly out of square. In any event, the brace obstructs the panel, so the semiannual ritual includes loosening the brace and removing four screws. During the last 60+ years, the screw holes have required repair / filling several times; about five years ago, I plugged them with epoxy putty and drilled them to fit the screws.

That repair having aged out, I was about to renew the epoxy when I realized that I now have brass inserts that would work even better, if I replaced the original wood screws with 10-32 machine screws.

I cut the screws to the exact length using the brace and brass insert as a fixture:

Storm door – screw cutting

The vacuum cleaner nozzle to the lower right inhales the debris from the Dremel cutoff wheel that would otherwise fill the shop; I used up the last half of a wheel on four stainless steel screws.

Because each end of the brace has two screws, I knew that I couldn’t just drill out the four holes, plant four inserts, and be done with the job: the first insert on each end could go pretty nearly anywhere, but the second insert must match the brace hole spacing. The only way I know how to do that is to epoxy the first two inserts in place and let them cure, drill the other two holes slightly oversize, mount those inserts on the brace, butter them with epoxy, put the brace in place, tighten the first two screws, snug the brace, and hope I didn’t epoxy the brace to the door or the screws to the inserts.

Slips of waxed paper between the brace and the door prevented the first problem and oiling the screws prevented the second. It’s not the best-looking job I’ve ever done, but nobody will ever see the inserts behind the brace:

Storm door – inserts

Now, we’re ready for winter and I’m ready for spring!

Most likely, the new owners (whoever and whenever they may be) will never use these inserts, as they’ll replace all the windows & doors, plus sand & refinish the hardwood floors, before moving in …

2016-11-05
Kitchen Spatula Search

A long long time ago, we bought a kitchen spatula that’s served us well ever since:

Spatula Search – original

To give you an idea of how old that poor thing is, the back of the handle bears a Japan stamp. I’ve re-set the rivets several times, the blade has rusted as badly as you think, and we recently, very reluctantly, decided it has passed its best-used-by date.

The 3 x 4.5 inch blade is 19 mil = 0.45 mm plated carbon steel, stiff enough to remain flat and springy enough to bend a little, with a 9 inch = 230 mm steel handle ending in a plastic overmold.

These days, it’s essential to the cutting, flipping, and serving of the morning’s omelet-like substance, made of eggs, bacon, veggies, green leafy things, plus this-and-that, in the cast-iron pan. Mary chops the disk into quarters with the reasonably sharp edge, maneuvers the reasonably bendy blade under each quarter, flips them over, tops with bacon & cheese, pauses for consolidation & melting, then pops them onto plates. Yum!

Omelet in cast-iron pan

So we set out to buy a replacement.

Here’s what we’ve tried and rejected so far:

Spatula Search – overview

I’ve used this one for many years to flip pancakes on a succession of non-stick griddles, a service at which it excels. The edge isn’t sharp enough to cut the green-and-leafy and the completely inflexible blade cannot be maneuvered under the omelet quarters:

Spatula Search – heavy solid plastic

This one gets deployed for burgers and their ilk, also in the cast-iron pan. The blade, although sharp enough, is completely rigid:

Spatula Search – heavy slotted metal

On the other paw, a slightly concave 7 mil = 0.18 mm spring steel blade is much too thin and, well, springy. Although very sharp, you cannot apply enough cutting force without suddenly bending the blade and, if the omelet quarter isn’t positioned exactly right, the blade will bend underneath it and dump breakfast on the stovetop. The alert reader will notice a missing weld between the blade and the bottom wire handle:

Spatula Search – thin spring steel

This very thin plastic blade has similar problems with poor cut-ability and excessive flexibility:

Spatula Search – thin springy plastic

This one looked really promising and worked almost perfectly. Regrettably, its nylon blade bears a 400 °F rating and the bottom of the omelet reaches nearly 450 °F. You can see what happens to the reasonably sharp edge as it scrapes across the pan:

Spatula Search – heavy slotted nylon

The omelet cooks at the temperature it cooks at, which part of the specifications is not subject to further discussion.

So, we’re stumped. Having trawled the usual online and big-box stores, we’ve been unable to find a replacement. Simple steel blades aren’t available. Trendy silicone-bonded stainless steel blades combine the worst of all worlds: won’t cut and won’t flip. Pretty nearly anything you don’t see above seems obviously unsuitable for our simple needs: too big, too small, or too melty.

We’ll consider all recommendations and suggestions! Thanks …

2016-11-02
Monthly Science: CR2032 Lithium Cell Life

One of the Hobo dataloggers asked for a new battery during its most recent data dump. The old battery dates back to January 2015:

Maxell CR2032 lithium cell – 22 month life

That was when a batch of Energizer cells failed in quick succession: it wasn’t the datalogger’s fault. I’ve been handling the cells a bit more carefully, too, although that certainly doesn’t account for the much longer life.

With batteries, particularly from eBay, you definitely can’t tell what you’re going to get or how long it’ll last; that’s true of many things in life.

2016-11-01
Inside Another 9 V Battery

A long time ago, I discovered some quasi-AAAA cells inside 9 V batteries:

Inside a batteries.com 9V battery

It occurred to me that I should dismantle a defunct Rayovac Maximum 9 V alkaline battery from the most recent batch (*) to see what it looked like:

Rayovac Maximum 9V battery – interior

Surprise!

A closer look at those pancake cells:

Rayovac Maximum 9V battery – detail

They look like separate cells bonded into a stack, although there’s no easy way to probe the inter-cell contacts; the leftmost cell probably died first.

(*) Which has apparently outlived the Rayovac Maximum brand, as they don’t appear on the Rayovac site.

2016-10-31
Hand Sprayer Hose Kink Prevention

Mary’s new half-gallon sprayer arrived with a kink in the hose just below the handle, which is about what you’d expect from a non-reinforced plastic tube jammed into the smallest possible box containing both the sprayer and its wand. Fortunately, the Box o’ Springs had one that just fit the hose and jammed firmly into the handle:

Sprayer hose with kink-resisting spring

The kink slowly worked its way out after being surrounded by the spring and shouldn’t come back.

That was easy…

2016-10-26

Raspberry Pi Streaming Radio Player: Improved Pipe Handling

My Raspberry Pi-based streaming radio player generally worked fine, except sometimes the keypad / volume control knob would stop responding after switching streams. This being an erratic thing, the error had to be a timing problem in otherwise correct code and, after spending Quality Time with the Python subprocess and select doc, I decided I was abusing mplayer’s stdin and stdout pipes.

This iteration registers mplayer’s stdout pipe as Yet Another select.poll() Polling Object, so that the main loop can respond whenever a complete line arrives. Starting mplayer in quiet mode reduces the tonnage of stdout text, at the cost of losing the streaming status that I really couldn’t do anything with, and eliminates the occasional stalls when mplayer (apparently) dies in the middle of a line.

The code kills and restarts mplayer whenever it detects an EOF or stream cutoff. That works most of the time, but a persistent server or network failure can still send the code into a sulk. Manually selecting a different stream (after we eventually notice the silence) generally sets things right, mainly by whacking mplayer upside the head; it’s good enough.

It seems I inadvertently invented streaming ad suppression by muting (most of) the tracks that produced weird audio effects. Given that the “radio stations” still get paid for sending ads to me, I’m not actually cheating anybody out of their revenue: I’ve just automated our trips to the volume control knob. The audio goes silent for a few seconds (or, sheesh, a few minutes) , blatting a second or two of ad noise around the gap to remind us of what we’re missing; given the prevalence of National Forest Service PSAs, the audio ad market must be a horrific wasteland.

The Python source code as a GitHub Gist:

	from evdev import InputDevice,ecodes,KeyEvent
	import subprocess32 as subp
	import select
	import re
	import sys
	import time
	import logging

	Media = {'KEY_KP7' : ['Classical',False,['mplayer','–quiet','-playlist','http://stream2137.init7.net/listen.pls'%5D%5D,
	'KEY_KP8' : ['Jazz',False,['mplayer','–quiet','-playlist','http://stream2138.init7.net/listen.pls'%5D%5D,
	'KEY_KP9' : ['WMHT',False,['mplayer','–quiet','http://live.str3am.com:2070/wmht1'%5D%5D,
	'KEY_KP4' : ['Classic 1000',True,['mplayer','–quiet','-playlist','http://listen.radionomy.com/1000classicalhits.m3u'%5D%5D,
	'KEY_KP5' : ['DCNY 911',False,['mplayer','–quiet','-playlist','http://www.broadcastify.com/scripts/playlists/1/12305/-5857889408.m3u'%5D%5D,
	'KEY_KP6' : ['WAMC',False,['mplayer','–quiet','http://pubint.ic.llnwd.net/stream/pubint_wamc'%5D%5D,
	'KEY_KP1' : ['60s',True,['mplayer','–quiet','-playlist','http://listen.radionomy.com/all60sallthetime-keepfreemusiccom.m3u'%5D%5D,
	'KEY_KP2' : ['50-70s',True,['mplayer','–quiet','-playlist','http://listen.radionomy.com/golden-50-70s-hits.m3u'%5D%5D,
	'KEY_KP3' : ['Soft Rock',True,['mplayer','–quiet','-playlist','http://listen.radionomy.com/softrockradio.m3u'%5D%5D,
	'KEY_KP0' : ['Zen',True,['mplayer','–quiet','-playlist','http://listen.radionomy.com/zen-for-you.m3u'%5D%5D
	}

	CurrentKC = 'KEY_KP3'

	Controls = {'KEY_KPSLASH' : '//////',
	'KEY_KPASTERISK' : '******',
	'KEY_KPENTER' : ' ',
	'KEY_KPMINUS' : '<',
	'KEY_KPPLUS' : '>',
	'KEY_VOLUMEUP' : '*',
	'KEY_VOLUMEDOWN' : '/'
	}

	MuteStrings = ["TargetSpot","[Unknown]","Advert:","+++","—","SRR","Srr","ZEN FOR"]

	MuteDelay = 8.0 # delay before non-music track; varies with buffering
	UnMuteDelay = 7.5 # delay after non-music track
	Muted = False # keep track of muted state

	MixerChannel = 'PCM' # which amixer thing to use

	logging.basicConfig(format='%(asctime)s %(levelname)s:%(message)s',filename='/tmp/Streamer.log',level=logging.INFO)
	logger = logging.getLogger()

	# set up event inputs and polling objects
	# This requires some udev magic to create the symlinks

	k = InputDevice('/dev/input/keypad')
	k.grab()

	kp = select.poll()
	kp.register(k.fileno(),select.POLLIN + select.POLLPRI + select.POLLERR)

	v = InputDevice('/dev/input/volume')
	v.grab()

	vp = select.poll()
	vp.register(v.fileno(),select.POLLIN + select.POLLPRI + select.POLLERR)

	# set up file for output tracing

	lw = open('/tmp/mp.log','w') # mplayer piped output

	# set the mixer output low enough that the initial stream won't wake the dead

	subp.call(['amixer','sset',MixerChannel,'10'])

	# Start the player with the default stream, set up for polling

	print 'Starting mplayer on',Media[CurrentKC][0],' -> ',Media[CurrentKC][-1][-1]
	logging.info('Starting mplayer on %s -> %s',Media[CurrentKC][0],Media[CurrentKC][-1][-1])
	p = subp.Popen(Media[CurrentKC][-1],
	stdin=subp.PIPE,stdout=subp.PIPE,stderr=subp.STDOUT)

	pp = select.poll() # this may be valid for other invocations, but is not pretty
	pp.register(p.stdout.fileno(),select.POLLIN + select.POLLPRI + select.POLLERR)

	print ' … running'

	#——————–
	#— Play the streams

	while True:

	# pluck next line from mplayer and decode it

	if [] != pp.poll(10):
	text = p.stdout.readline()

	if 'ICY Info: ' in text: # these lines may contain track names
	lw.write(text)
	lw.flush()

	trkinfo = text.split(';') # also splits at semicolon embedded in track name
	# logging.info('Raw split line: %s', trkinfo)
	for ln in trkinfo:
	if 'StreamTitle' in ln: # this part contains a track name
	NeedMute = False # assume a listenable track
	if 2 == ln.count("'"): # exactly two single quotes = probably none embedded in track name
	trkhit = re.search(r"StreamTitle='(.*)'",ln)
	if trkhit: # true for valid search results
	TrackName = trkhit.group(1) # the string between the two quotes
	print 'Track name: ', TrackName
	logging.info('Track name: [%s]', TrackName)
	if Media[CurrentKC][1] and ( (len(TrackName) == 0) or any(m in TrackName for m in MuteStrings) ) :
	NeedMute = True
	else:
	print ' … semicolon in track name: ', ln
	logging.info('Semicolon in track name: [' + ln + ']')
	else:
	print ' … quotes in track name: ', ln
	logging.info('Quotes in track name: [' + ln + ']')
	if NeedMute:
	print ' … muting:',
	if Media[CurrentKC][1] and not Muted:
	time.sleep(MuteDelay) # brute-force assumption about buffer leadtime
	subp.call(['amixer','-q','sset',MixerChannel,'mute'])
	Muted = True
	print 'done'
	logging.info('Track muted')
	else:
	print ' … unmuting:',
	if Muted:
	if Media[CurrentKC][1]:
	time.sleep(UnMuteDelay) # another brute-force timing assumption
	Muted = False
	subp.call(['amixer','-q','sset',MixerChannel,'unmute'])
	print 'done'
	logging.info('Track unmuted')

	elif 'Exiting.' in text: # mplayer just imploded
	lw.write(text)
	lw.flush()
	print 'Got EOF / stream cutoff'
	logging.info('EOF or stream cutoff')
	print ' … killing dead mplayer'
	pp.unregister(p.stdout.fileno())
	p.terminate() # p.kill()
	p.wait()
	# print ' … flushing pipes'
	# lw.truncate(0)
	print ' … discarding keys'
	while [] != kp.poll(0):
	kev = k.read
	print ' … restarting mplayer: ',Media[CurrentKC][0]
	logging.info('Restarting mplayer')
	p = subp.Popen(Media[CurrentKC][-1],
	stdin=subp.PIPE,stdout=subp.PIPE,stderr=subp.STDOUT)
	pp.register(p.stdout.fileno(),select.POLLIN + select.POLLPRI + select.POLLERR)
	print ' … running'
	logging.info(' … running')


	# accept pending events from volume control knob

	if [] != vp.poll(10):
	vev = v.read()
	lw.write('Volume')
	lw.flush()
	for e in vev:
	if e.type == ecodes.EV_KEY:
	kc = KeyEvent(e).keycode

	# print 'Volume kc: ',kc

	if kc in Controls:
	print 'Vol Control: ',kc
	try:
	p.stdin.write(Controls[kc])
	except Exception as e:
	print "Can't send control: ",e
	print ' … restarting player: ',Media[CurrentKC][0]
	logging.info('Error sending volume, restarting player')
	pp.unregister(p.stdout.fileno())
	p = subp.Popen(Media[CurrentKC][-1],
	stdin=subp.PIPE,stdout=subp.PIPE,stderr=subp.STDOUT)
	pp.register(p.stdout.fileno(),select.POLLIN + select.POLLPRI + select.POLLERR)
	print ' … running'
	logging.info(' … running')


	# accept pending events from keypad

	if [] != kp.poll(10):
	kev = k.read()
	lw.write("Keypad")
	lw.flush()
	for e in kev:
	if e.type == ecodes.EV_KEY:
	kc = KeyEvent(e).keycode

	if kc == 'KEY_NUMLOCK': # discard these, as we don't care
	continue

	# print 'Got: ',kc

	if (kc == 'KEY_BACKSPACE') and (KeyEvent(e).keystate == KeyEvent.key_hold):
	if True:
	print 'Backspace = shutdown!'
	p.kill()
	logging.shutdown()
	q = subp.call(['sudo','shutdown','-P','now'])
	q.wait()
	time.sleep(5)
	print "Oddly, we did not die…"
	else:
	print 'BS = bail from main!'
	logging.shutdown()
	sys.exit(0)
	break

	if KeyEvent(e).keystate != KeyEvent.key_down: # discard key up & other rubbish
	continue

	if kc in Controls:
	print 'Control:', kc
	try:
	p.stdin.write(Controls[kc])
	except Exception as e:
	print "Can't send control: ",e
	print ' … restarting player: ',Media[CurrentKC][0]
	logging.info('Error sending controls, restarting player')
	pp.unregister(p.stdout.fileno())
	p.terminate() # p.kill()
	p.wait()
	p = subp.Popen(Media[CurrentKC][-1],
	stdin=subp.PIPE,stdout=subp.PIPE,stderr=subp.STDOUT)
	pp.register(p.stdout.fileno(),select.POLLIN + select.POLLPRI + select.POLLERR)
	print ' … running'
	logging.info(' … running')

	if kc in Media:
	print 'Switching stream to ',Media[kc][0],' -> ',Media[kc][-1][-1]
	logging.info('Switching stream: ' + Media[kc][0] + ' -> ' + Media[kc][-1][-1])
	CurrentKC = kc
	print ' … halting player'
	try:
	p.communicate(input='q')
	except Exception as e:
	print 'Perhaps mplayer died?',e
	print ' … killing it for sure'
	pp.unregister(p.stdout.fileno())
	p.terminate() # p.kill()
	p.wait()
	# print ' … flushing pipes'
	# lw.truncate(0)
	print ' … restarting player: ',Media[CurrentKC][0]
	p = subp.Popen(Media[CurrentKC][-1],
	stdin=subp.PIPE,stdout=subp.PIPE,stderr=subp.STDOUT)
	pp.register(p.stdout.fileno(),select.POLLIN + select.POLLPRI + select.POLLERR)
	print ' … running'
	logging.info(' … running')

	print 'Out of loop!'
	logging.shutdown()

view raw Streamer.py hosted with ❤ by GitHub

2016-10-21

HP 7475A Plotter: Coordinate Pruning

The original SuperFormula equation produces points in polar coordinates, which the Chiplotle library converts to the rectilinear format more useful with Cartesian plotters. I’ve been feeding the equation with 10001 angular values (10 passes around the paper, with 1000 points per pass, plus one more point to close the pattern), which means the angle changes by 3600°/10000 = 0.36° per point. Depending on the formula’s randomly chosen parameters, each successive point can move the plotter pen by almost nothing to several inches.

On the “almost nothing” end of the scale, the plotter slows to a crawl while the serial interface struggles to feed the commands. Given that you can’t see the result, why send the commands?

Computing point-to-point distances goes more easily in rectilinear coordinates, so I un-tweaked my polar-modified superformula function to return the points in rectangular coordinates. I’d originally thought a progressive scaling factor would be interesting, but it never happened.

The coordinate pruning occurs in the supershape function, which now contains a loop to scan through the incoming list of points from the superformula function and add a point to the output path only when it differs by enough from the most recently output point:

    path = []
    path.append(Coordinate(width * points[0][0], height * points[0][1]))
    outi = 0
    xp, yp = points[outi][0], points[outi][1]
    for i in range(len(points))[1:]:
        x,y = width * points[i][0], height * points[i][1]
        dist = sqrt(pow(x - xp,2) + pow(y - yp,2))
        if dist > 60 :
          path.append(Coordinate(x, y))
          outi = i
          xp, yp = x, y

    path.append(Coordinate(width * points[-1][0], height * points[-1][1]))
    print "Pruned",len(points),"to",len(path),"points"

The first and last points always go into the output list; the latter might be duplicated, but that doesn’t matter.

Note that you can’t prune the list by comparing successive points, because then you’d jump directly from the start of a series of small motions to their end. The idea is to step through the small motions in larger units that, with a bit of luck, won’t be too ugly.

The width and height values scale the XY coordinates to fill either A or B paper sheets, with units of “Plotter Units” = 40.2 PU/mm = 1021 PU/inch. You can scale those in various ways to fit various output sizes within the sheets, but I use the defaults that fill the entire sheets with a reasonable margin. As a result, the magic number 60 specifies 60 Plotter Units; obviously, it should have a suitable name.

Pruning to 40 PU = 1.0 mm (clicky for more dots, festooned with over-compressed JPEG artifacts):

Pruning to 60 PU = 1.5 mm:

Pruning to 80 PU = 2.0 mm:

Pruning to 120 PU = 3.0 mm:

All four of those plots have the same pens in the same order, although I refilled a few of them in flight.

By and large, up through 80 PU there’s not much visual difference, although you can definitely see the 3 mm increments at 120 PU. However, the plotting time drops from just under an hour for each un-pruned plot to maybe 15 minutes with 120 PU pruning, with 60 PU producing very good results at half an hour.

Comparing the length of the input point lists to the pruned output path lists, including some pruning values not shown above:

Prune 20
1 - m: 5.3, n1: 0.15, n2=n3: 0.80
Pruned 10001 to 4856 points
2 - m: 5.3, n1: 0.23, n2=n3: 0.75
Pruned 10001 to 5545 points
3 - m: 5.3, n1: 1.15, n2=n3: 0.44
Pruned 10001 to 6218 points
4 - m: 5.3, n1: 0.41, n2=n3: 1.50
Pruned 10001 to 7669 points
5 - m: 5.3, n1: 0.29, n2=n3: 0.95
Pruned 10001 to 6636 points
6 - m: 5.3, n1: 0.95, n2=n3: 0.16
Pruned 10001 to 5076 points

Prune 40
1 - m: 3.1, n1: 0.23, n2=n3: 0.26
Pruned 10001 to 2125 points
2 - m: 3.1, n1: 1.05, n2=n3: 0.44
Pruned 10001 to 5725 points
3 - m: 3.1, n1: 0.25, n2=n3: 0.32
Pruned 10001 to 2678 points
4 - m: 3.1, n1: 0.43, n2=n3: 0.34
Pruned 10001 to 4040 points
5 - m: 3.1, n1: 0.80, n2=n3: 0.40
Pruned 10001 to 5380 points
6 - m: 3.1, n1: 0.55, n2=n3: 0.56
Pruned 10001 to 5424 points

Prune 60
1 - m: 1.1, n1: 0.45, n2=n3: 0.40
Pruned 10001 to 2663 points
2 - m: 1.1, n1: 0.41, n2=n3: 0.14
Pruned 10001 to 1706 points
3 - m: 1.1, n1: 1.20, n2=n3: 0.75
Pruned 10001 to 4446 points
4 - m: 1.1, n1: 0.33, n2=n3: 0.80
Pruned 10001 to 3036 points
5 - m: 1.1, n1: 0.90, n2=n3: 1.40
Pruned 10001 to 4723 points
6 - m: 1.1, n1: 0.61, n2=n3: 0.65
Pruned 10001 to 3601 points

Prune 80
1 - m: 3.7, n1: 0.95, n2=n3: 0.58
Pruned 10001 to 3688 points
2 - m: 3.7, n1: 0.49, n2=n3: 0.22
Pruned 10001 to 2258 points
3 - m: 3.7, n1: 0.57, n2=n3: 0.90
Pruned 10001 to 3823 points
4 - m: 3.7, n1: 0.25, n2=n3: 0.40
Pruned 10001 to 2161 points
5 - m: 3.7, n1: 0.47, n2=n3: 0.30
Pruned 10001 to 2532 points
6 - m: 3.7, n1: 0.45, n2=n3: 0.14
Pruned 10001 to 1782 points

Prune 120
1 - m: 1.9, n1: 0.33, n2=n3: 0.48
Pruned 10001 to 1561 points
2 - m: 1.9, n1: 0.51, n2=n3: 0.18
Pruned 10001 to 1328 points
3 - m: 1.9, n1: 1.80, n2=n3: 0.16
Pruned 10001 to 2328 points
4 - m: 1.9, n1: 0.21, n2=n3: 1.10
Pruned 10001 to 1981 points
5 - m: 1.9, n1: 0.63, n2=n3: 0.24
Pruned 10001 to 1664 points
6 - m: 1.9, n1: 0.45, n2=n3: 0.22
Pruned 10001 to 1290 points

Eyeballometrically, 60 PU pruning halves the number of plotted points, so the average data rate jumps from 9600 b/s to 19.2 kb/s. Zowie!

Most of the pruning occurs near the middle of the patterns, where the pen slows to a crawl. Out near the spiky rim, where the points are few & far between, there’s no pruning at all. Obviously, quantizing a generic plot to 1.5 mm would produce terrible results; in this situation, the SuperFormula produces smooth curves (apart from those spikes) that look just fine.

The Python source code as a GitHub Gist:

	# Adapted from Chiplotle plotter library:
	# http://cmc.music.columbia.edu/chiplotle/

	from chiplotle import *
	from math import *
	from datetime import *
	from time import *
	from types import *
	import random

	def supershape(width, height, m, n1, n2, n3,
	point_count=10 * 1000, percentage=1.0, a=1.0, b=1.0, travel=None):
	'''Supershape, generated using the superformula first proposed
	by Johan Gielis.

	– `points_count` is the total number of points to compute.
	– `travel` is the length of the outline drawn in radians.
	3.1416 * 2 is a complete cycle.
	modified to prune short plotter motions – Ed Nisley KE4ZNU – October 2016
	'''
	travel = travel or (10 * 2 * pi)

	# compute points…
	phis = [i * travel / point_count
	for i in range(1 + int(point_count * percentage))]
	points = [tools.mathtools.superformula(a, b, m, n1, n2, n3, x) for x in phis]

	# scale and prune short motions
	path = []
	path.append(Coordinate(width * points[0][0], height * points[0][1]))
	outi = 0
	xp, yp = points[outi][0], points[outi][1]
	for i in range(len(points))[1:]:
	x,y = width * points[i][0], height * points[i][1]
	dist = sqrt(pow(x – xp,2) + pow(y – yp,2))
	if dist > 60 :
	path.append(Coordinate(x, y))
	outi = i
	xp, yp = x, y

	path.append(Coordinate(width * points[-1][0], height * points[-1][1]))
	print " Pruned",len(points),"to",len(path),"points"

	return Path(path)


	# Run Superformula plots

	if __name__ == '__main__':

	override = False

	plt = instantiate_plotters()[0]
	# plt.write('IN;')

	if plt.margins.soft.width < 11000: # A=10365 B=16640
	maxplotx = (plt.margins.soft.width / 2) – 100
	maxploty = (plt.margins.soft.height / 2) – 150
	legendx = maxplotx – 2900
	legendy = -(maxploty – 750)
	tscale = 0.45
	numpens = 4
	# prime/10 = number of spikes
	m_values = [n / 10.0 for n in [11, 13, 17, 19, 23]]
	# ring-ness 0.1 to 2.0, higher is larger
	n1_values = [
	n / 100.0 for n in range(55, 75, 2) + range(80, 120, 5) + range(120, 200, 10)]
	else:
	maxplotx = plt.margins.soft.width / 2
	maxploty = plt.margins.soft.height / 2
	legendx = maxplotx – 3000
	legendy = -(maxploty – 900)
	tscale = 0.45
	numpens = 6
	m_values = [n / 10.0 for n in [11, 13, 17, 19, 23, 29, 31,
	37, 41, 43, 47, 53, 59]] # prime/10 = number of spikes
	# ring-ness 0.1 to 2.0, higher is larger
	n1_values = [
	n / 100.0 for n in range(15, 75, 2) + range(80, 120, 5) + range(120, 200, 10)]

	print " Max: ({},{})".format(maxplotx, maxploty)

	# spiky-ness 0.1 to 2.0, higher is spiky-er (mostly)
	n2_values = [
	n / 100.0 for n in range(10, 60, 2) + range(65, 100, 5) + range(110, 200, 10)]

	plt.write(chr(27) + '.H200:') # set hardware handshake block size
	plt.set_origin_center()
	# scale based on B size characters
	plt.write(hpgl.SI(tscale * 0.285, tscale * 0.375))
	# slow speed for those abrupt spikes
	plt.write(hpgl.VS(10))

	while True:

	# standard loadout has pen 1 = fine black
	plt.write(hpgl.PA([(legendx, legendy)]))
	pen = 1
	plt.select_pen(pen)
	plt.write(hpgl.PA([(legendx, legendy)]))
	plt.write(hpgl.LB("Started " + str(datetime.today())))

	if override:
	m = 4.1
	n1_list = [1.15, 0.90, 0.25, 0.59, 0.51, 0.23]
	n2_list = [0.70, 0.58, 0.32, 0.28, 0.56, 0.26]
	else:
	m = random.choice(m_values)
	n1_list = random.sample(n1_values, numpens)
	n2_list = random.sample(n2_values, numpens)

	pen = 1
	for n1, n2 in zip(n1_list, n2_list):
	n3 = n2
	print "{0} – m: {1:.1f}, n1: {2:.2f}, n2=n3: {3:.2f}".format(pen, m, n1, n2)
	plt.select_pen(pen)
	plt.write(hpgl.PA([(legendx, legendy – 100 * pen)]))
	plt.write(
	hpgl.LB("Pen {0}: m={1:.1f} n1={2:.2f} n2=n3={3:.2f}".format(pen, m, n1, n2)))
	e = supershape(maxplotx, maxploty, m, n1, n2, n3)
	plt.write(e)
	pen = pen + 1 if (pen % numpens) else 1

	pen = 1
	plt.select_pen(pen)
	plt.write(hpgl.PA([(legendx, legendy – 100 * (numpens + 1))]))
	plt.write(hpgl.LB("Ended " + str(datetime.today())))
	plt.write(hpgl.PA([(legendx, legendy – 100 * (numpens + 2))]))
	plt.write(hpgl.LB("More at https://softsolder.com/?s=7475a"))
	plt.select_pen(0)
	plt.write(hpgl.PA([(-maxplotx,maxploty)]))

	print "Waiting for plotter… ignore timeout errors!"
	sleep(40)
	while NoneType is type(plt.status):
	sleep(5)

	print "Load more paper, then …"
	print " … Press ENTER on the plotter to continue"
	plt.clear_digitizer()
	plt.digitize_point()

	plotstatus = plt.status
	while (NoneType is type(plotstatus)) or (0 == int(plotstatus) & 0x04):
	plotstatus = plt.status

	print "Digitized: " + str(plt.digitized_point)

view raw PlotterShapes.py hosted with ❤ by GitHub

2016-10-20

Tag: Improvements

Storm Door Brace: Now, With Inserts!

Kitchen Spatula Search

Monthly Science: CR2032 Lithium Cell Life

Inside Another 9 V Battery

Hand Sprayer Hose Kink Prevention

Raspberry Pi Streaming Radio Player: Improved Pipe Handling

HP 7475A Plotter: Coordinate Pruning