import logging import matplotlib.pyplot as plt import time from threading import Thread import numpy from PIL import Image class Tracer(object): def __init__(self, fname): super(Tracer, self).__init__() self.threshold = 64 self.orig = Image.open(fname) self.conv1 = self.remove_transparency(self.orig) self.conv1 = self.conv1.convert("L") self.src = numpy.asarray(self.conv1) self.dst = numpy.copy(self.src) self.dst.fill(False) self.mark = numpy.copy(self.dst) # start in center self.height, self.width = self.dst.shape[0:2] self.posy = self.height / 2 self.posx = self.width / 2 def remove_transparency(self, im, bg_colour=(255, 255, 255)): # from https://stackoverflow.com/questions/35859140/remove-transparency-alpha-from-any-image-using-pil # Only process if image has transparency (http://stackoverflow.com/a/1963146) if im.mode in ('RGBA', 'LA') or (im.mode == 'P' and 'transparency' in im.info): # Need to convert to RGBA if LA format due to a bug in PIL (http://stackoverflow.com/a/1963146) alpha = im.convert('RGBA').split()[-1] # Create a new background image of our matt color. # Must be RGBA because paste requires both images have the same format # (http://stackoverflow.com/a/8720632 and http://stackoverflow.com/a/9459208) bg = Image.new("RGBA", im.size, bg_colour + (255,)) bg.paste(im, mask=alpha) return bg else: return im def trace(self): while self._has_unchecked_pixels(): # go circles to find a pixel in src if not self._spiral_till_pixel(): break # move until we find new pixels self._move_while_you_can() def _has_unchecked_pixels(self): ix, iy = numpy.where(self.mark == False) return len(ix) or len(iy) def _spiral_till_pixel(self): radius = 1 direction = 0 offset = 0 while self._has_unchecked_pixels(): in_lower = self.posy < self.height and self.posx < self.width in_upper = self.posy >= 0 and self.posx >= 0 if in_lower and in_upper and not self.mark[self.posy][self.posx]: if self.src[self.posy][self.posx] < self.threshold: return True self.mark[self.posy][self.posx] = True if direction == 0: self.posx += 1 self.posy += 0 elif direction == 1: self.posx += 0 self.posy += 1 elif direction == 2: self.posx += -1 self.posy += 0 elif direction == 3: self.posx += 0 self.posy += -1 else: pass offset += 1 if offset >= radius: # time.sleep(0.01) offset = 0 direction += 1 if direction > 3: direction = 0 if direction in (0, 2): radius += 1 logging.debug("End of image") return False def _move_while_you_can(self): # time.sleep(0.1) logging.debug("%s:%s=%s", self.posy, self.posx, self.src[self.posy][self.posx]) self.mark[self.posy][self.posx] = True self.dst[self.posy][self.posx] = True # self.posx += 1 # self.posy += 1 class TracerVisualizer(object): def __init__(self, tracer): """ :type tracer: Tracer """ self.tracer = tracer def run(self): tracer = self.tracer fig, ((ax1, ax2), (ax3, ax4)) = plt.subplots(nrows=2, ncols=2) ax1.imshow(tracer.orig) ax2.imshow(tracer.src) plt.show(block=False) thr = Thread(target=tracer.trace) thr.setDaemon(True) thr.start() while plt.get_fignums(): # weird trick to react on close ax3.set_title("%s:%s" % (tracer.posx, tracer.posy)) ax3.imshow(tracer.mark, cmap='gray') ax4.imshow(tracer.dst, cmap='gray') plt.pause(0.5) if __name__ == "__main__": logging.basicConfig(level=logging.DEBUG) trc = Tracer("test1.png") TracerVisualizer(trc).run() time.sleep(5)