{

pixbuf = Gdk::Pixbuf::create(Gdk::COLORSPACE_RGB, true, 8, width, height);

cairo_surface = Cairo::ImageSurface::create(pixbuf->get_pixels(),

Cairo::FORMAT_ARGB32, width, height, pixbuf->get_rowstride());

ivalue_t max_y, string name, string outputdir) :

dpi(dpi), min_x(min_x), max_x(max_x), min_y(min_y), max_y(max_y), zero_x(

-min_x * (ivalue_t) dpi + (ivalue_t) procmargin), zero_y(

-min_y * (ivalue_t) dpi + (ivalue_t) procmargin),

name(name), outputdir(outputdir), fill(false), clr(32)

{

guint8* pixels;

int stride;

make_the_surface((max_x - min_x) * dpi + 2 * procmargin,

(max_y - min_y) * dpi + 2 * procmargin);

usedcolors.push_back(BLACK);

usedcolors.push_back(WHITE);

/* "Note that the buffer is not cleared; you will have to fill it completely yourself." */

//printf("clearing\n");

pixels = cairo_surface->get_data();

stride = cairo_surface->get_stride();

for (int y = 0; y < pixbuf->get_height(); y++)

{

for (int x = 0; x < pixbuf->get_width(); x++)

{

PRC(pixels + x*4 + y*stride) = BLACK;

}

}

throw (import_exception)

{

importer->render(cairo_surface, dpi,

min_x - static_cast<ivalue_t>(procmargin) / dpi,

min_y - static_cast<ivalue_t>(procmargin) / dpi);

if (fill)

fill_outline();

const icoordpair &lb)

{

icoordpair nab; //normal vector to a-b= {-ab_y,ab_x}

nab.first = -(la.second - lb.second);

nab.second = (la.first - lb.first);

double lnab = sqrt(nab.first * nab.first + nab.second * nab.second);

double skalar; //product

skalar = nab.first * (x.first - la.first)

+ nab.second * (x.second - la.second);

return fabs(skalar / lnab);

bool mirrored)

{

Isolator* iso = dynamic_cast<Isolator*>(mill.get());

int extra_passes = iso ? iso->extra_passes : 0;

coords components = fill_all_components();

int added = -1;

int contentions = 0;

int grow = mill->tool_diameter / 2 * dpi;

ivalue_t mirror_axis = mill->mirror_absolute ? min_x : ((min_x + max_x) / 2);

vector<shared_ptr<icoords> > toolpath;

for (int pass = 0; pass <= extra_passes && added != 0; pass++)

{

for (int i = 0; i < grow && added != 0; i++)

{

added = 0;

for ( coordpair c : components )

{

added += grow_a_component(c.first, c.second, contentions);

}

}

coords inside, outside;

for ( coordpair c : components )

{

calculate_outline(c.first, c.second, outside, inside);

inside.clear();

shared_ptr<icoords> outline(new icoords());

shared_ptr<icoords> outline_optimised(new icoords());

// i'm not sure wheter this is the right place to do this...

// that "mirrored" flag probably is a bad idea.

for ( coordpair c : outside )

{

outline->push_back(

icoordpair(

// tricky calculations

mirrored ?

(2 * mirror_axis - xpt2i(c.first)) :

xpt2i(c.first),

min_y + max_y - ypt2i(c.second)));

}

outside.clear();

if (mill->optimise)

{

//Use Boost's Douglas-Peucker simplification algorithm

boost::geometry::simplify( *outline, *outline_optimised, 1.0 / dpi );

toolpath.push_back(outline_optimised);

}

else

toolpath.push_back(outline);

}

}

if (contentions)

{

cerr << "\nWarning: pcb2gcode hasn't been able to fulfill all"

<< " clearance requirements and tried a best effort approach"

<< " instead. You may want to check the g-code output and"

<< " possibly use a smaller milling width.\n";

}

tsp_solver::nearest_neighbour( toolpath, std::make_pair(0, 0), 1.0 / dpi );

save_debug_image("traced_" + name);

return toolpath;

{

bool badcol;

do

{

badcol = false;

clr = rand();

for (unsigned int i = 0; i < usedcolors.size(); i++)

{

if (usedcolors[i] == clr)

{

badcol = true;

break;

}

}

}

while (badcol);

usedcolors.push_back(clr);

return clr;

{

std::vector<pair<int, int> > components;

int max_x = cairo_surface->get_width() - 1;

int max_y = cairo_surface->get_height() - 1;

guint8* pixels = cairo_surface->get_data();

int stride = cairo_surface->get_stride();

for (int y = 0; y <= max_y; y++)

{

for (int x = 0; x <= max_x; x++)

{

if ((PRC(pixels + x*4 + y*stride) | OPAQUE) == WHITE)

{

components.push_back(pair<int, int>(x, y));

fill_a_component(x, y, get_an_unused_color());

}

}

}

return components;

{

guint32 newclr = argb;

guint8* pixels = cairo_surface->get_data();

int stride = cairo_surface->get_stride();

guint8* here = pixels + x * 4 + y * stride;

guint8* maxhere = pixels + (cairo_surface->get_width() - 1) * 4

+ (cairo_surface->get_height() - 1) * stride;

guint32 ownclr = PRC(here);

std::stack<pair<int, int> > queued_pixels;

queued_pixels.push(pair<int, int>(x, y));

while (queued_pixels.size())

{

pair<int, int> current_pixel = queued_pixels.top();

x = current_pixel.first;

y = current_pixel.second;

queued_pixels.pop();

here = pixels + x * 4 + y * stride;

PRC(here) = newclr;

if (here + 4 <= maxhere && PRC(here+4) == ownclr)

queued_pixels.push(pair<int, int>(x + 1, y));

if (pixels <= here - 4 && PRC(here-4) == ownclr)

queued_pixels.push(pair<int, int>(x - 1, y));

if (here + stride <= maxhere && PRC(here+stride) == ownclr)

queued_pixels.push(pair<int, int>(x, y + 1));

if (pixels <= here - stride && PRC(here-stride) == ownclr)

queued_pixels.push(pair<int, int>(x, y - 1));

if (here + 4 + stride <= maxhere && PRC(here+4+stride) == ownclr)

queued_pixels.push(pair<int, int>(x + 1, y + 1));

if (here - 4 + stride <= maxhere && PRC(here-4+stride) == ownclr)

queued_pixels.push(pair<int, int>(x - 1, y + 1));

if (pixels <= here + 4 - stride && PRC(here+4-stride) == ownclr)

queued_pixels.push(pair<int, int>(x + 1, y - 1));

if (pixels <= here - 4 - stride && PRC(here-4-stride) == ownclr)

queued_pixels.push(pair<int, int>(x - 1, y - 1));

}

cairo_surface->mark_dirty();

{

guint8* pixels = cairo_surface->get_data();

int stride = cairo_surface->get_stride();

guint32 start_color = PRC(pixels + x*4 + y * stride);

if (start_color == 0)

{

PRC(pixels + x*4 + y * stride) = RED;

save_debug_image("error_runtoborder");

std::stringstream msg;

msg << "run_to_border: start_color == 0 at (" << x << "," << y << ")\n";

throw std::logic_error(msg.str());

}

while (PRC(pixels + x*4 + y*stride) == start_color)

x++;

{

if (x <= 0 || y <= 0)

return false;

if (x >= cairo_surface->get_width() - 1)

return false;

if (y >= cairo_surface->get_height() - 1)

return false;

guint8* pixels = cairo_surface->get_data();

int stride = cairo_surface->get_stride();

for (int i = 7; i >= 0; i--)

{

int cx = x + offset8[i][0];

int cy = y + offset8[i][1];

guint8* pixel = pixels + cx * 4 + cy * stride;

// surrounding pixel != own color, not black -> other component!

if (PRC(pixel) != ownclr && (PRC(pixel) | OPAQUE) != BLACK)

return false;

}

return true;

vector<pair<int, int> >& outside, vector<pair<int, int> >& inside)

{

guint8* pixels = cairo_surface->get_data();

int stride = cairo_surface->get_stride();

int max_y = cairo_surface->get_height();

guint32 owncolor = PRC(pixels + x*4 + y*stride);

int xstart = x;

int ystart = y;

run_to_border(xstart, ystart); //change xstart so that xstart++ would be outside of the component

int xout = xstart;

int yout = ystart;

int xin = xout - 1;

int yin = yout;

outside.push_back(pair<int, int>(xout, yout));

int blasts = 0;

while (true)

{

int i;

int steps = 0; // number of steps done in 1 iteration of the while loop

// step outside

for (i = 0; i < 8; i++)

{

int xoff = xout - xin + 1;

int yoff = yout - yin + 1;

int xnext = xin + growoff_o[xoff][yoff][0];

int ynext = yin + growoff_o[xoff][yoff][1];

if (xnext == xstart && ynext == ystart)

{

outside.push_back(pair<int, int>(xout, yout));

outside.push_back(pair<int, int>(xstart, ystart));

return;

}

if (xnext < 0 || ynext < 0 || xnext > stride || ynext > max_y)

{

save_debug_image("error_outerpath");

std::stringstream msg;

msg << "Outside path reaches image margins at " << xin << ","

<< yin << ")\n";

throw std::logic_error(msg.str());

}

guint8* next = pixels + xnext * 4 + ynext * stride;

if (PRC(next) != owncolor)

{

outside.push_back(pair<int, int>(xout, yout));

xout = xnext;

yout = ynext;

}

else

break;

}

if (i == 8)

{

save_debug_image("error_outsideoverstepping");

std::stringstream msg;

msg << "Outside over-stepping at in(" << xin << "," << yin << ")\n";

throw std::logic_error(msg.str());

}

steps += i;

// step inside

for (i = 0; i < 8; i++)

{

int xoff = xin - xout + 1;

int yoff = yin - yout + 1;

int xnext = xout + growoff_i[xoff][yoff][0];

int ynext = yout + growoff_i[xoff][yoff][1];

// next pixel is checked clockwise

guint8* next = pixels + xnext * 4 + ynext * stride;

if (xnext < 0 || ynext < 0 || xnext > stride || ynext > max_y)

{

save_debug_image("error_innerpath");

std::stringstream msg;

msg << "Inside path reaches image margins at " << xin << ","

<< yin << ")\n";

throw std::logic_error(msg.str());

}

if (PRC(next) == owncolor)

{

inside.push_back(pair<int, int>(xin, yin));

xin = xnext;

yin = ynext;

}

else

break;

}

if (i == 8)

{

save_debug_image("error_insideoverstepping");

std::stringstream msg;

msg << "Inside over-stepping at out(" << xout << "," << yout

<< ")\n";

throw std::logic_error(msg.str());

}

steps += i;

// check whether we made any progress calculating the trace outline.

// if we haven't, our algorithm is deadlocked by stray pixels

// we try to resolve this by enforcing the algorithm's constraints

// for the components

if (steps == 0)

{

int changes = 0;

// test constraints for surrounding pixels, enforce if necessary

for (i = 0; i < 8; i++)

{

int cx = xin + offset8[i][0];

int cy = yin + offset8[i][1];

guint8* pixel = pixels + cx * 4 + cy * stride;

if (allow_grow(cx, cy, owncolor))

{

PRC(pixel) = owncolor;

changes++;

}

// if a component pixel can't be reached non-diagonally, clear it

// even if it was set just now

int j;

for (j = 0; j < 4; j++)

{

if (PRC( pixels + (cx + offset4[j][0]) * 4

+ (cy + offset4[j][1]) * stride ) != BLACK)

break;

}

if (j == 4)

{

PRC(pixel) = BLACK;

changes++;

}

}

if (allow_grow(xstart, ystart, owncolor))

PRC(pixels+xstart*4+ystart*stride) = owncolor;

if (changes == 0)

{

PRC(pixels + xin*4 + yin*stride) |= RED;

PRC(pixels + xout*4 + yout*stride) |= BLUE;

save_debug_image("failed_repair");

std::stringstream msg;

msg << "Failed repairing @ (" << xin << "," << yin << ")\n";

throw std::logic_error(msg.str());

}

else

blasts++;

// start right at the beginning. still more efficient than keeping

// the history necessary to be able to continue next to the problem.

inside.clear();

outside.clear();

xstart = x;

ystart = y;

run_to_border(xstart, ystart);

xout = xstart;

yout = ystart;

xin = xout - 1;

yin = yout;

outside.push_back(pair<int, int>(xout, yout));

continue;

}

}

fprintf(stderr, "blasts: %d", blasts);

{

if (x < 0 || x >= cairo_surface->get_width() || y < 0

|| y >= cairo_surface->get_height())

{

std::stringstream msg;

msg << "grow_a_component(): invalid starting point: (" << x << "," << y

<< ")";

throw std::logic_error(msg.str());

}

contentions = 0;

vector<pair<int, int> > outside, inside;

calculate_outline(x, y, outside, inside);

guint8* pixels = cairo_surface->get_data();

int stride = cairo_surface->get_stride();

unsigned int pixels_changed = 0;

guint32 ownclr = PRC(pixels + x*4 + y*stride);

for (unsigned int i = 0; i < outside.size(); i++)

{

pair<int, int> coord = outside[i];

if (allow_grow(coord.first, coord.second, ownclr))

{

PRC(pixels + coord.first*4 + coord.second*stride) = ownclr;

pixels_changed++;

}

else

{

contentions++;

}

}

return pixels_changed;

{

shared_ptr<Surface> mask = dynamic_pointer_cast<Surface>(mask_surface);

if (mask)

{

Cairo::RefPtr<Cairo::ImageSurface> mask_cairo_surface = mask->cairo_surface;

int max_x = cairo_surface->get_width();

int max_y = cairo_surface->get_height();

int stride = cairo_surface->get_stride();

if (max_x != mask_cairo_surface->get_width()

|| max_y != mask_cairo_surface->get_height()

|| stride != mask_cairo_surface->get_stride())

{

throw std::logic_error("Surface shapes don't match.");

}

guint8* pixels = cairo_surface->get_data();

guint8* mask_pixels = mask_cairo_surface->get_data();

for (int y = 0; y < max_y; y++)

{

for (int x = 0; x < max_x; x++)

{

PRC(pixels + x*4 + y*stride) &=

PRC(mask_pixels + x*4 + y*stride); /* engrave only on the surface area */

PRC(pixels + x*4 + y*stride) |=

(~PRC(mask_pixels + x*4 + y*stride) & (RED | BLUE)); /* tint the outiside in an own color to block extension */

}

}

}

else

throw std::logic_error("Can't cast Core to Surface");

{

static unsigned int debug_image_index = 0;

opacify(pixbuf);

pixbuf->save( build_filename(outputdir,

(boost::format("outp%1%_%2%.png") % debug_image_index % message).str() ),

"png");

debug_image_index++;

{

int stride = pixbuf->get_rowstride();

guint8* pixels = pixbuf->get_pixels();

for (int y = 0; y < pixbuf->get_height(); y++)

{

for (int x = 0; x < pixbuf->get_width(); x++)

{

PRC(pixels + x*4 + y*stride) |= OPAQUE;

}

}

{

fill = true;

this->linewidth = linewidth;

{

/* paint everything white that can not be reached from outside the image */

int stride = pixbuf->get_rowstride();

guint8* pixels = pixbuf->get_pixels();

/* in order to find out what is "outside", we need to walk "around' the image */

for (int x = 0; x < pixbuf->get_width(); x++)

{

if (PRC(pixels + x*4 + 0*stride) != BLACK)

throw std::logic_error("Non-black pixel at top border");

if (PRC(pixels + x*4 + (pixbuf->get_height() - 1)*stride) != BLACK)

throw std::logic_error("Non-black pixel at bottom border");

}

for (int y = 0; y < pixbuf->get_height(); y++)

{

if (PRC(pixels + 0*4 + y*stride) != BLACK)

throw std::logic_error("Non-black pixel at left border");

if (PRC(pixels + (pixbuf->get_width() - 1)*4 + y*stride) != BLACK)

throw std::logic_error("Non-black pixel at right border");

}

fill_a_component(0, 0, BLUE);

/* everything else (that is, the area of the board) will be black

int first_line_with_black = 0;

for (int y = 0; y < pixbuf->get_height(); y++)

{

for (int x = 0; x < pixbuf->get_width(); x++)

{

if (PRC(pixels + x*4 + y*stride) != BLUE)

{

PRC(pixels + x*4 + y*stride) = BLACK;

if (first_line_with_black == 0)

{

first_line_with_black = y;

}

}

}

}

/* compensate for growth induced by line thicknesses.

int grow = linewidth / 2 * dpi;

int contentions = 0;

int added = 0;

for (int i = 0; i < grow; ++i)

{

// starting at the very left

added = grow_a_component(0, first_line_with_black + grow, contentions);

}

// if you can think of a sane situation in which either of this could

// occur and nevertheless give a meaningful result, change it to a

// warning.

if (!added)

throw std::logic_error(

"Shrinking the outline by half the line width came to a halt.");

if (contentions)

throw std::logic_error(

"Shrinking the outline collided with something while there should not be anything.");

for (int y = 0; y < pixbuf->get_height(); y++)

{

for (int x = 0; x < pixbuf->get_width(); x++)

{

if (PRC(pixels + x*4 + y*stride) == BLUE)

PRC(pixels + x*4 + y*stride) = BLACK;

else

PRC(pixels + x*4 + y*stride) = WHITE;

}

}