Advances in Fuzzy Systems

Volume 2013 (2013), Article ID 593694, 6 pages

http://dx.doi.org/10.1155/2013/593694

Research Article

Generating Suitable Basic Functions Used in Image Reconstruction by F-Transform

Pavel Vlašánek

Department of Informatics and Computers, University of Ostrava, 30. Dubna 22, 701 03 Ostrava, Czech Republic

Received 5 July 2013; Accepted 12 July 2013

Academic Editor: Salvatore Sessa

Copyright © 2013 Pavel Vlašánek. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

Image reconstruction technique based on F-transform uses clearly defined basic functions. These functions have strong impact on the quality of reconstruction. We can use some predefined shape and radius, but also we can create a new one from the scratch. The aim of this paper is to analyze the creating process and based on that find best basic function for input set of damaged testing images.

1. Introduction

Image reconstruction aims at recovering damaged parts. In real situation this damage can be caused by various sources on various surfaces. Input information has to be digitalized as a first step and after that divided on damaged and nondamaged parts. Typical situation lays on photography. We can distinguish among many types of damage like stains, scratches, text, or noise. We can say that every kind of damage covers typical way of damaging process. Unwanted time or date stamp is covered by text. Scratches, cracks, and folds are covered by inpaint. These are shown in Figure 1 which includes also noise.

(a)

(b)

(c)

Figure 1: (a) Inpaint damage; (b) text damage; (c) noise damage.

Target of reconstruction is removing damaged parts from input image and replacing them by the parts with recomputed values. These values are computed from the neighborhood ones. The technique mentioned in this paper is based on F-transform which brings particular way of valuation neighborhood pixels and their usage in the computation [1, 2]. Other methods [3, 4] differ by the choice at the technique.

In this paper, we will focus on the valuation part of the computation. Quality of reconstruction will be measured by RMSE value (RMSE stands for the root mean square error). In Figure 2 you can see damaged input images and reconstructed ones with usage of the ideal basic function described later. We will analyse influence of changing parameter on target quality.

(a)

(b)

(c)

Figure 2: (a) Original image; (b) damaged image; (c) reconstructed image.

2. F-Transform

In image reconstruction, a discrete version of the F-transform is used. Details can be seen in [5, 6]. In this section two-dimensional (2D) variant and also conditions for proper functioning will be briefly introduced.

2.1. Fuzzy Partition with Ruspini Condition

Let be fixed nodes within such that , and . We say that the fuzzy sets , identified with their membership functions defined on , establish a fuzzy partition with Ruspini condition of if they fulfill the following conditions for :

(1),  ;(2) if , where for uniformity of notation, we set and ;(3) is continuous;(4), for , increases on , and , for , strictly decreases on ;(5)for all ,

The condition (1) is known as the Ruspini condition. The membership functions are called basic functions. A point is covered by basic function if .

2.2. Shape of the Basic Function

The shape of the basic functions is not predetermined, and therefore, it can be chosen according to additional requirements (e.g., smoothness). Let us give examples of various fuzzy partitions with the Ruspini condition. In Figure 3, two such partitions with triangular and cosine basic functions are shown. The formulas given below represent generic fuzzy partitions with the Ruspini condition and triangular functions:

where and .

(a)

(b)

Figure 3: (a) Ruspini partitions with triangular basic function; (b) Ruspini partitions with cosine basic function.

We say that a Ruspini partition of is-uniform if its nodes , where , are -equidistant; that is, , for , where , and two additional properties are met:

(6), for all , ;(7), for all and , and , for all and .

An -uniform fuzzy partition of can be determined by the so called generating function , which is assumed to be even (The function is even if for all ,  ), continuous, have a bell shape, and fulfill . Basic functions of an -uniform fuzzy partition with generating function are shifted copies of in the sense that

and for ,

As an example, we notice that the function is a generating function for any -uniform triangular partition. In the sequel, we will be using -uniform fuzzy partitions only and refer to as to a radius of partition [7].

2.3. Discrete 2D F-Transform

We say that the -matrix of real numbers is called the (discrete) F-transform ofwith respect to and if for all ,

The elements are called components of the F-transform.

The inverse F-transform is defined as follows:

We use discrete basic functions and . Our results show that sufficient image reconstruction is enough to use maximal radius . Discrete basic function with respect to Ruspini condition is building on top of template in Figure 4.

Figure 4: Template for basic function definition.

You can see white squares in Figure 4. These squares can be marked as parts of the basic function. One marked per column. After that is automatically computed mirrored part of the basic function as can be seen in Figure 5. We choose for the column . Value in column is the same, and because of Ruspini condition, values in columns 5 and 7 are easily computable as

(a)

(b)

Figure 5: (a) Result after selection of fourth pixel in column; (b) result after selection of fifth pixel in column.

Because of using Ruspini condition, we can use copies of the basic function for covering the whole range. Height of the template determines values on the axis. For better usability we choose value 0 for first row, 1 for last one, and twice 0.5 on the middle. Two basic functions with marked values can be seen in Figure 6.

Figure 6: Template for basic function definition.

The overlapping part of basic functions have functional values equal 1 in each column. For better visualization, black vertical lines are also plotted.

3. Influence of Basic Functions

Shapes of the basic functions can be nondecreasing, oscillating, or nonincreasing. Differences and influences on computation are described in [8]. We choose basic functions with radius 2 in the first step of reconstruction and radius 4 in the second step because of fully sufficient usage on the input set of testing images. In this paper, we focus on building shape of basic function step by step based on results provided by RMSE. We identified two ways of this process which will be demonstrated on template in Figure 4 where .

Column by Column. Find the best value for current column and continue with next column:

(1)choose radius , (2)choose active column as (3)choose current row as , (4)based on those values create a basic function by mirroring and with respect to Ruspini condition formula (7), (5)use the basic function for image reconstruction, (6)compare reconstructed image with original undamaged one by RMSE, (7)if , then , and continue by step 3, (8)change current row as , (9)continue by step 4.

Radius by Radius. Find the best value for current radius and continue with next radius:

(1)choose current radius as , (2)choose current column as (3)choose current row as , (4)based on those values create a basic function, (5)use the basic function for image reconstruction, (6)compare reconstructed image with the original undamaged by RMSE, (7)if and , then and continue by step 3, (8)if , then and continue by step 3, (9)choose current row as , (10)continue by step 4.

First way column by column is computed in Table 1. You can see that there is for column best RMSE value for . It means that we choose sixth row for column . Next column is the best RMSE value in row for the Lena and Nature images and for the Picture image. We choose basic function with values for column and for column . This function is in Figure 8.

Table 1: RMSE values of the column by column basic function creation.

Lena
Picture
Nature

1
62.99
5.16
46.72
4.01
71.25
7.55

2
49.99
5.04
36.63
3.89
56.39
7.39

3
49.98
4.97
36.62
3.82
56.37
7.30

4
49.97
4.93
36.62
3.79
56.36
7.27

5
49.97
4.92
36.61
3.80
56.35
7.27

6
49.96
4.95
36.61
3.84
56.35
7.32

Second way of basic function creation radius by radius is computed in Table 2. From the table comes that the best value is for and for column . This basic function is shown in Figure 7.

Table 2: RMSE values of the radius by radius basic function creation.

1
2
3
4
5
6

Lena

1
6.12
5.71
5.58
5.51
5.50
5.52

2
5.16
4.97
4.85
4.79
4.78
4.81

3
5.12
4.96
4.86
4.81
4.80
4.83

4
5.11
4.97
4.88
4.83
4.83
4.86

5
5.13
5.00
4.92
4.87
4.87
4.90

6
5.16
5.04
5.97
4.93
4.92
4.95

Picture

1
4.93
4.32
4.21
4.16
4.15
4.19

2
3.93
3.76
3.65
3.61
3.61
3.65

3
3.91
3.76
3.67
3.63
3.64
3.68

4
3.92
3.79
3.71
3.67
3.68
3.72

5
3.96
3.83
3.76
3.73
3.73
3.77

6
4.01
3.89
3.82
3.80
3.80
3.84

Nature

1
8.50
7.68
7.45
7.33
7.30
7.33

2
7.63
7.33
7.16
7.07
7.06
7.10

3
7.53
7.30
7.15
7.09
7.09
7.13

4
7.50
7.30
7.18
7.13
7.14
7.19

5
7.51
7.34
7.23
7.19
7.20
7.25

6
7.55
7.39
7.30
7.27
7.27
7.32

(a)

(b)

Figure 7: (a) Basic function for radius ; (b) basic function for radius .

(a)

(b)

Figure 8: (a) Basic function for radius; (b) basic function for radius .

4. Conclusion

As a result we can say that step by step process converges to oscillating or linear shape. Column by column way is the best solution oscillating basic function. Value for column and for column provides RMSE values for Lena, for Picture, and for Nature. Because of very close results for Picture image between and , we choose value for all of them. Radius by radius is the best linear function with for column and for column . For Lena RMSE is equal to , for Picture is equal to , and for Nature is equal to .

Acknowledgment

This work was supported by SGS14/PRF/2013 (advanced techniques of applications of soft computing methods in image processing).

References

P. Vlasanek, I. Perfilieva, and M. Wrublova, “Fuzzy transform for image reconstruction, Uncertainty modeling in knowledge engineering and decision making,” in Proceedings of the 10th International FLINS Conference, pp. 615–620, World Scientific, Istanbul, Turkey, August 2012.

P. Vlasanek and I. Perfilieva, “Image reconstruction with usage of the F-Transform,” in Proceedings of the Advances in Intelligent Systems and Computing, Interna- tional Joint Conference (CISIS'12-ICEUTE'12-SOCO '12), A. Herrero, V. Snášel, A. Abraham, et al., Eds., vol. 189, pp. 507–514, Springer, 2012.

J. A. Parker, R. V. Kenyon, and D. E. Troxel, “Comparison of interpolating methods for image resampling,” IEEE Transactions on Medical Imaging, vol. MI-2, no. 1, pp. 31–39, 1983.

K. Uhlir and V. Skala, “Radial basis function use for the restoration of damaged images,” in Computer Vision and Graphics, K. Wojciechowski, B. Smolka, H. Palus, et al., Eds., vol. 32, pp. 839–844, Kluwer Academic, 2006.

F. Di Martino, V. Loia, I. Perfilieva, and S. Sessa, “An image coding/decoding method based on direct and inverse fuzzy transforms,” International Journal of Approximate Reasoning, vol. 48, no. 1, pp. 110–131, 2008.

I. Perfilieva, “Fuzzy transforms: theory and applications,” Fuzzy Sets and Systems, vol. 157, no. 8, pp. 993–1023, 2006.

I. Perfilieva and P. Vlasanek, Image Reconstruction by means of F-transform, Elsevier Editorial System for Knowledge-Based Systems, 2013

P. Vlasanek and I. Perfilieva, “Influence of various types of basic functions on image reconstruction using F-transform,” European Society for Fuzzy Logic and Technology. In press