Elsevier

Graphical Models

Volume 76, Issue 5, September 2014, Pages 484-495
Graphical Models

A graph-based optimization algorithm for fragmented image reassembly

https://doi.org/10.1016/j.gmod.2014.03.001Get rights and content

Abstract

We propose a graph-based optimization framework for automatic 2D image fragment reassembly. First, we compute the potential matching between each pair of the image fragments based on their geometry and color. After that, a novel multi-piece matching algorithm is proposed to reassemble the overall image fragments. Finally, the reassembly result is refined by applying the graph optimization algorithm. We perform experiments to evaluate our algorithm on multiple torn real-world images, and demonstrate the robustness of this new assembly framework outperforms the existing algorithms in both reassembly accuracy (in handling accumulated pairwise matching error) and robustness (in handling small image fragments).

Introduction

Fragmented image reassembly recomposes a group of picture fragments into the original complete image. It is a geometric processing problem that has important applications in many fields such as archaeology and forensics. For example, archaeologists need to spend a lot of efforts to recompose fractured ancient paintings to restore their original appearances; forensic specialists manually compose damaged documents or pictures to recover the original evidence that are broken by humans. Developing robust geometric algorithms for automatic image reassembly can reduce the expensive human labors and improve the restoration efficiency, and hence is highly desirable.

Existing automatic reassembly algorithms can be generally divided into two categories: color-based approaches and geometry-based approaches. Color-based methods mainly use the color information to predict the adjacency relationship of the fragments and guide the matching [1], [2], [3]. These algorithms are usually efficient, but they sometimes suffer from composition accuracy and may fail when the textures of different fragments are similar. Geometry-based methods reassemble the fragments by matching their boundary curves [4], [5], [6], [7]. They can more accurately align adjacent fragments along their breaking regions, but they are sometimes slow and can easily get trapped in local optima then fail to reach the correct result. Therefore, incorporating both the geometry and color information in the reassembly computation could make process more efficient, accurate, and reliable. This paper presents a novel 3-step composition algorithm, as illustrated in Fig. 1.

The first step is called the pairwise matching, which aims to identify adjacent pairs and compute their initial alignments. Geometry-based pairwise matching methods rely on analyzing the shape of the boundary curve contours; color-based pairwise matching methods match fragments using their color information. Our integrated algorithm (1) extracts the border of each fragment and represents it as a curve contour, (2) clusters such a curve contour into multiple segments based on both color and geometry, and (3) matches contour segments to suggest the potential alignment between adjacent fragments. This pairwise matching algorithm will result in a set of possible matching between identified pairs of fragments, some of which are correct while some are not.

The pairwise matching produces redundant matches, which we intentionally generate in order to tolerate erroneous adjacency identification due to noise and local minima. A global groupwise matching can better filter out these false matches. Extensive examination of this is a combinatorial optimization problem and is NP-hard. Many previous works use a best-first search strategy (i.e., to explore a graph by expanding the most promising node) with backtracking to solve this problem. In our work, the global composition is formulated as a novel graph-based searching problem, solving which will give us a more reliable global reassembly of multiple fragments.

Groupwise matching can compose the fragments globally. However, the result is dependent on the composition sequence, and the alignment errors will be accumulated which may even cause a failure in the reassembly (i.e., fragment intersection, see Section 6 for details). We derive a graph optimization algorithm to reduce the global matching errors between adjacent fragments.

The main contributions of this work are as follows:

  • 1.

    We integrate both geometry and color information in pairwise matching computation to obtain more reliable alignment between adjacent fragments.

  • 2.

    We propose a graph-based algorithm that performs groupwise matching to better handle the errors resulted from pairwise alignments and obtain correct adjacency information of all the fragments.

  • 3.

    We develop a variational graph optimization in the end to reduce the accumulated errors to refine the reassembly and achieve a global optimal result.

Section snippets

Local pairwise matching

The essence of fragment reassembly is to find the relative transformations between adjacent fragments. In the procedure of 2D images composition, fragments can be modeled and aligned using their boundaries, which are 2D curve contours. Therefore, pairwise matching often reduces to a partial curve matching problem, which is mainly solved in existing literatures via either geometry-based or color-based approaches.

Wolfson [8] approximates the curves as the polygonal shapes and solves their

Problem formulation and experimental settings

Given a set of fragments F={Fi},i=1,,n,FiR2, which are from a torn image I, we want to solve the transformation for each fragment Fi, i.e., a 3×3 matrix Ti, so that these transformations compose all the fragments back into the original image I=Ti(Fi).

In our experiments, each image is torn into multiple fragmented pieces. We scan all the fragments. For each fragment Fi we obtain a digital image scan Ii which has a colorful region appeared on a white background. We can segment such a region Fi

Pairwise matching between fragments

This section elaborates our algorithm of computing the pairwise matching between two image fragments. First, the boundary of each fragment Fi is extracted from the scan image Ii; the boundary is a 2D curve loop. Then, the pairwise matching between two image fragments can reduce to a partial curve matching problem. Our pairwise matching algorithm utilizes both the color and geometry information. Unlike existing algorithms that directly use dynamic programming [5] or geometric hashing [8] for

Global image reassembly

The pairwise matching step provides us a set of possible matchings between each pair of the image fragments; each matching is associated with a matching score. This score is usually a good indicator to whether this alignment is correct. A straightforward strategy is to iteratively pick the matching following the scores from high to low. However, this may not always work as the highest score is not a guarantee to the correct match, especially when there are many small fragments: small fragment’s

Refinement of overall reassembly

In the previous section, we obtain a global reassembly of all the fragmented images. But small error may be inevitable in the transformations computed in the previous two steps, caused during e.g. polygon approximation, the inconsistency of the pixelization in image scan. In Step 2, we start with one fragment and iteratively glue other fragments. Transformation errors will be accumulated and sometimes this could even cause a failure in the reassembly (Fig. 7(a)). The error becomes especially

Experimental results

We perform multiple experiments to evaluate our reassembly algorithm. In our experiments, we print out randomly selected digital images on papers. Most images, after printed, are with the size of about 20 cm × 25 cm. Then we randomly tear an image into multiple pieces. The size of the fragments is ranging from 4 cm × 4 cm–7 cm × 7 cm. Then each image fragment is scanned (in our experiments, with a scanner of 150-dpi resolution). Our reassembly algorithm is then used to recompose the scanned digital image

Conclusion and future work

We present a novel computational pipeline for the automatic reassembly of fragmented images. It consists of three main steps: pairwise matching between two image fragments, graph-based global fragment reassembly, and refinement of the reassembly via graph optimization. In Step 1, we present a better curve-matching algorithm: each pair of image fragments are aligned via a pairwise matching integrating both geometry and color; In step 2, we present a reliable graph-based global search algorithm,

Acknowledgments

This research is partially supported by National Science Foundation – United States (IIS-1320959), IBM Faculty Award, and National Natural Science Foundation of China (No. 61170323).

References (25)

  • L. Zhu et al.

    Globally consistent reconstruction of ripped-up documents

    IEEE Trans. Pattern Anal. Mach. Intell.

    (2008)
  • H. Wolfson

    On curve matching

    IEEE Trans. Pattern Anal. Mach. Intell.

    (1990)
  • Cited by (53)

    • Solving archaeological puzzles

      2021, Pattern Recognition
      Citation Excerpt :

      Handling 2D general-shaped pieces and gaps. Though not in archaeology, the most closely related works to ours are [31,49], which assume general-shaped pieces of natural images. However, in [31] there are no gaps between the pieces.

    • DAFNE: A dataset of fresco fragments for digital anastlylosis

      2020, Pattern Recognition Letters
      Citation Excerpt :

      More precisely, matching curves and contour-shape based methods were used to recompose 2D and 3D puzzles [6], pottery [7], wall paintings [8] and archaeological fragments [9]. Considering fragmented image reconstruction, the available solutions include pure color-based approaches [10], the combination of color and geometric features [11,12] and, more recently, also deep learning [13]. All these methods, however, deal with a situation similar to a puzzle, where all (or quite all) the pieces are available, and there is an almost perfect match between fragment edges.

    View all citing articles on Scopus
    View full text