Esempio n. 1
0
    def test_00_generate_layers(self):
        """Generate some layers and validate them"""
        print "\nGenerating a set of simple slices on my head."
        stl = helpers.get_head()
        surf = cairo.ImageSurface(cairo.FORMAT_ARGB32, 1024, 1024)
        ctx = cairo.Context(surf)
        ctx.scale(1024, 1024)
        ctx.set_source_rgb(0.8, 0.2, 0.5)
        for facet in stl._facets:
            x, y, z = facet['p'][0]
            ctx.move_to(0.5 + x / 30, 0.5 + y / 30)

            x, y, z = facet['p'][1]
            ctx.line_to(0.5 + x / 30, 0.5 + y / 30)
            x, y, z = facet['p'][2]
            ctx.line_to(0.5 + x / 30, 0.5 + y / 30)
        ctx.set_line_width(0.001)
        ctx.stroke()

        #outpath = os.path.join(
        #os.path.dirname(__file__),
        #'out', 'stl-topview.png')
        #surf.write_to_png(outpath)

        model = Model.from_stl(stl)
        #layers = model.generate_planar_intersections(-0, 1, 10)
        now = time.time()
        layers = model.generate_planar_intersections(0, 0.1, 30, processes=1)
        print "Simple planar intersections: That took %5.5f seconds" % (
            time.time() - now)
        p = Pool(cpu_count() * 4)
    def test_00_generate_layers(self):
        """Generate some layers and validate them"""
        print "\nGenerating a set of simple slices on my head."
        stl = helpers.get_head()
        surf = cairo.ImageSurface(cairo.FORMAT_ARGB32, 1024, 1024)
        ctx = cairo.Context (surf)
        ctx.scale (1024, 1024)
        ctx.set_source_rgb (0.8, 0.2, 0.5)
        for facet in stl._facets:
            x,y,z = facet['p'][0]
            ctx.move_to(0.5+x/30, 0.5+y/30)

            x,y,z = facet['p'][1]
            ctx.line_to(0.5+x/30,0.5+y/30)
            x,y,z = facet['p'][2]
            ctx.line_to(0.5+x/30,0.5+y/30)
        ctx.set_line_width (0.001)
        ctx.stroke()

        #outpath = os.path.join(
            #os.path.dirname(__file__),
            #'out', 'stl-topview.png')
        #surf.write_to_png(outpath)

        model = Model.from_stl(stl)
        #layers = model.generate_planar_intersections(-0, 1, 10)
        now=time.time()
        layers = model.generate_planar_intersections(0, 0.1, 30, processes=1)
        print "Simple planar intersections: That took %5.5f seconds" % (time.time()-now)
        p = Pool(cpu_count()*4)
    def test_02_complex_insets(self):
        """Test a complicated set of insets, on the head model."""
        print "\nPerforming a series of profiling and insetting, " \
                "on a model of rantenki's head."
        stl = helpers.get_head()
        model = Model.from_stl(stl)
        now=time.time()
        layers = [
                (plane[0], Layer.from_CGAL_intersections(plane[1])) for
                    plane in
                    model.generate_planar_intersections(0, 0.5, 30)
                    if len(plane[1])]
        print "Generating slices for %d layers took %3.3f seconds." % (len(layers), time.time()-now)



        now = time.time()
        for layer in layers:
            polys = list(layer[1].polys)

            for erosion in [0.15, 0.3, .45]:
                eroded = [poly.buffer(-erosion) for poly in polys]
        print "Insetting perimeters for 3 loops for %d layers took %3.3f seconds." % (len(layers), time.time()-now)


        return #We don't need to dump files now.
        now = time.time()
        for layer in layers:
            outpath = os.path.join(
                os.path.dirname(__file__),
                'out',
                "layer-%3.3fz.png" % layer[0])
            surf = get_cairo_surf(outpath)
            polys = list(layer[1].polys)
            draw_cairo_polys(polys, surf, (0,0,0.5,1))

            for erosion in [0.15, 0.3, .45]:
                eroded = [poly.buffer(-erosion) for poly in polys]
                draw_cairo_polys(eroded, surf, (0.2, 2*erosion, 0.2, 1))

                for internal in poly.interiors:
                    draw_cairo_polys(internal, surf, (0.2, 0.2, 2*erosion, 1))
            surf.write_to_png(outpath)
        print "Generating hires pngs of the sliced/inset layers took %3.3f seconds." % (time.time()-now)
Esempio n. 4
0
    def test_02_complex_insets(self):
        """Test a complicated set of insets, on the head model."""
        print "\nPerforming a series of profiling and insetting, " \
                "on a model of rantenki's head."
        stl = helpers.get_head()
        model = Model.from_stl(stl)
        now = time.time()
        layers = [(plane[0], Layer.from_CGAL_intersections(plane[1]))
                  for plane in model.generate_planar_intersections(0, 0.5, 30)
                  if len(plane[1])]
        print "Generating slices for %d layers took %3.3f seconds." % (
            len(layers), time.time() - now)

        now = time.time()
        for layer in layers:
            polys = list(layer[1].polys)

            for erosion in [0.15, 0.3, .45]:
                eroded = [poly.buffer(-erosion) for poly in polys]
        print "Insetting perimeters for 3 loops for %d layers took %3.3f seconds." % (
            len(layers), time.time() - now)

        return  #We don't need to dump files now.
        now = time.time()
        for layer in layers:
            outpath = os.path.join(os.path.dirname(__file__), 'out',
                                   "layer-%3.3fz.png" % layer[0])
            surf = get_cairo_surf(outpath)
            polys = list(layer[1].polys)
            draw_cairo_polys(polys, surf, (0, 0, 0.5, 1))

            for erosion in [0.15, 0.3, .45]:
                eroded = [poly.buffer(-erosion) for poly in polys]
                draw_cairo_polys(eroded, surf, (0.2, 2 * erosion, 0.2, 1))

                for internal in poly.interiors:
                    draw_cairo_polys(internal, surf,
                                     (0.2, 0.2, 2 * erosion, 1))
            surf.write_to_png(outpath)
        print "Generating hires pngs of the sliced/inset layers took %3.3f seconds." % (
            time.time() - now)