OneBitImageView* to_logical(const T &src)
{
OneBitImageData* data = new OneBitImageData(src.size(), src.origin());
OneBitImageView* view = new OneBitImageView(*data);
transform(src.vec_begin(),
src.vec_end(),
view->vec_begin(),
to_logical_core<typename T::value_type>());
return view;
}
OneBitImageView* border_removal(const T &src,
int win_dil, int win_avg, int win_med,
double threshold1_scale, double threshold1_gradient,
double threshold2_scale, double threshold2_gradient,
double scale_length,
int terminate_time1, int terminate_time2, int terminate_time3,
unsigned int interval2, unsigned int interval3)
{
// image resize
double scalar=sqrt(double(AREA_STANDARD)/(src.nrows()*src.ncols()));
GreyScaleImageView* src_scale=static_cast<GreyScaleImageView*>(scale(src, scalar, 1));
// paper estimation
GreyScaleImageView* blur1=paper_estimation(*src_scale, SMOOTH, win_dil, win_avg, win_med);
GreyScaleImageView* blur2=paper_estimation(*src_scale, DETAIL, win_dil, win_avg, win_med);
// edge detection
OneBitImageView* boundary=edge_detection(*blur1, *blur2,
threshold1_scale, threshold1_gradient,
threshold2_scale, threshold2_gradient,
scale_length);
// boundary reconstruct
OneBitImageView* mask_scale=boundary_reconstruct(*boundary,
terminate_time1, terminate_time2, terminate_time3,
interval2, interval3);
// image resize back
OneBitImageView* mask_temp=static_cast<OneBitImageView*>(scale(*mask_scale, 1.0/scalar, 1));
OneBitImageData* data = new OneBitImageData(src.size(), src.origin());
OneBitImageView* mask = new OneBitImageView(*data);
for (unsigned int m=0; m<mask->nrows(); m++) {
for (unsigned int n=0; n<mask->ncols(); n++) {
mask->set(Point(n, m), mask_temp->get(Point(n, m)));
}
}
delete src_scale->data();
delete src_scale;
delete blur1->data();
delete blur1;
delete blur2->data();
delete blur2;
delete boundary->data();
delete boundary;
delete mask_scale->data();
delete mask_scale;
delete mask_temp->data();
delete mask_temp;
return mask;
}
OneBitImageView* add_edge(const T &src, unsigned int interval)
{
OneBitImageData* data = new OneBitImageData(src.size(), src.origin());
OneBitImageView* view = new OneBitImageView(*data);
copy(src.vec_begin(),
src.vec_end(),
view->vec_begin());
for (unsigned int y=interval; y<view->nrows()-interval; y++) {
view->set(Point(0, y), 1);
view->set(Point(view->ncols()-1, y), 1);
}
for (unsigned int x=interval; x<view->ncols()-interval; x++) {
view->set(Point(x, 0), 1);
view->set(Point(0, view->nrows()-1), 1);
}
return view;
}
OneBitImageView* flatten(const T &src )
{
typedef typename T::value_type value_type;
int x,y;
OneBitImageData* dest_data = new OneBitImageData(src.size(), src.origin());
OneBitImageView* dest = new OneBitImageView(*dest_data);
OneBitPixel blackval = black(*dest);
OneBitPixel whiteval = white(*dest);
for (y = 0; y < (int)src.nrows(); y++)
for (x = 0; x < (int)src.ncols(); x++) {
if (src.get(Point(x,y))!=0)
dest->set(Point(x,y),blackval);
else
dest->set(Point(x,y),whiteval);
}
return dest;
}
Image* labeled_region_edges(const T& src, bool mark_both=false) {
OneBitImageData* edges_data = new OneBitImageData(src.size(), src.origin());
OneBitImageView* edges = new OneBitImageView(*edges_data);
size_t x,y,max_x,max_y;
max_x = src.ncols()-1;
max_y = src.nrows()-1;
// the following mask is sufficient: xx
// xx
// because we assume that no pixel is unlabeled
// check bulk of image
for (y=0; y<max_y; ++y) {
for (x=0; x<max_x; ++x) {
if (src.get(Point(x,y)) != src.get(Point(x+1,y))) {
edges->set(Point(x,y),1);
if (mark_both)
edges->set(Point(x+1,y),1);
}
if (src.get(Point(x,y)) != src.get(Point(x,y+1))) {
edges->set(Point(x,y),1);
if (mark_both)
edges->set(Point(x,y+1),1);
}
if (src.get(Point(x,y)) != src.get(Point(x+1,y+1))) {
edges->set(Point(x,y),1);
if (mark_both)
edges->set(Point(x+1,y+1),1);
}
}
}
// check last row
for (x=0; x<max_x; ++x) {
if (src.get(Point(x,max_y)) != src.get(Point(x+1,max_y))) {
edges->set(Point(x,max_y),1);
if (mark_both)
edges->set(Point(x+1,max_y),1);
}
}
// check last column
for (y=0; y<max_y; ++y) {
if (src.get(Point(max_x,y)) != src.get(Point(max_x,y+1))) {
edges->set(Point(max_x,y),1);
if (mark_both)
edges->set(Point(max_x,y+1),1);
}
}
return edges;
}