## Services on Demand

## Journal

## Article

## Indicators

- Cited by SciELO
- Access statistics

## Related links

- Cited by Google
- Similars in SciELO
- Similars in Google

## Share

## Materials Research

##
*Print version* ISSN 1516-1439*On-line version* ISSN 1980-5373

### Mat. Res. vol.11 no.4 São Carlos Oct./Dec. 2008

#### https://doi.org/10.1590/S1516-14392008000400006

**REGULAR ARTICLES**

**Hydrophobicity classification of polymeric materials based on fractal dimension**

**Daniel Thomazini ^{I, *}; Maria Virginia Gelfuso^{I}; Ruy Alberto Corrêa Altafim^{II}**

^{I}CCT - Engenharia de Controle e Automação, Universidade de Fortaleza - UNIFOR, Av. Washington Soares, 1321, 60811-905 Fortaleza - CE, Brazil

^{II}Departamento de Engenharia Elétrica, Universidade de São Paulo - USP, Av. Trabalhador Sancarlense, 400, 13566-590 São Carlos - SP, Brazil

**ABSTRACT**

This study proposes a new method to obtain hydrophobicity classification (HC) in high voltage polymer insulators. In the method mentioned, the HC was analyzed by fractal dimension (fd) and its processing time was evaluated having as a goal the application in mobile devices. Texture images were created from spraying solutions produced of mixtures of isopropyl alcohol and distilled water in proportions, which ranged from 0 to 100% volume of alcohol (%AIA). Based on these solutions, the contact angles of the drops were measured and the textures were used as patterns for fractal dimension calculations.

**Keywords:** digital image, hydrophobicity, texture, fractal dimension

**1. Introduction**

In the last decades, new polymeric materials such as silicone rubber, polyurethane and epoxy^{1-3} have been used for high voltage outdoor insulation. The main characteristics of these materials are the low cost, easiness of molding, chemical stability and high resistance to natural conditions.

Among other characteristics, they also present low surface energy which can be explained by the presence of methyl groups on their surface. In fact, these materials are water-repellent, also called hydrophobic^{4}.

The hydrophobicity classification (HC) is defined by the STRI guide^{5}, by measuring the contact angle or the amount of the wetted insulator surface. This method, generally used in the field has an intrinsic inaccuracy due to its subjective image analyses.

About 30 years ago, mathematicians as Von Kock, Peano, Hausdorff, Besicovitch and others ones created pictures which were called "monsters" or "pathologic" by some others mathematicians^{6} because these pictures were created by fractal geometry calculations. Recently, the fractal geometry has been applied in digital image processing methods to image-file compression^{7}, segmentation techniques^{8}, surface roughness characterizations in fractures^{9}, wear and erosion surfaces^{10}.

Some researchers have proposed digital image analyses to classify the hydrophobicity of the electrical insulators. Circular factor^{11,12}, goniometric measurement using Hough transformation^{13}, segmentation, scaled entropy and histogram analyses^{14} and surface energy^{15} were some of the proposed methods, but neither of them were able to obtain a mathematical relationship which may define HC. This paper presents a new approach where this HC identification is possible using the fractal dimension as a parameter to quantify the texture of the images and classify the HC.

**2. Experimental Procedure**

The material employed in the experiment procedure was silicon rubber outdoor insulators. Ten specimen were cut out from each sheet, each one with (10 x 10) cm^{2} and 5 mm thickness. The solutions of water and isopropyl alcohol (%AIA solution) were sprayed on samples, which ranged from 0 to 100% volume in alcohol. In this way, discrete texture images were obtained (%AIA image), providing a possibility of comparison and classification according to the STRI Guide.

At first, ten grayscale VGA images from each solution of the drops and the surfaces (OR) were taken and the contact angles were measured. In order to avoid differences or gradients in the illumination, caused by natural conditions, the original images were modified by histogram equalization (EQU) and/or white top-hat filter (WTH). The EQU promotes the balance in the grayscale distribution; in other words, all range of the grayscale is presented in the image. Another way to suppress the conditions of illumination is the edge detection tool, which produces images with the drop patterns, with the Sobel operator (SOB)^{16} being used in this case. The contact angle was measured from the drops with the aid of a CAD software.

The WTH works as a high-pass filter, reducing the gradient illumination of an image normally obtained in natural conditions. In order to apply the WTH filter, the opening image should be initially defined from Equation 1:

where δ* _{B}* is the dilatation and ε

*is the erosion of an f image.*

_{B}Equation 2 presents the mathematical operator of the WTH. In order to define the best opening size for WTH images, the opening mask size (*B*) has been changed from 3 x 3 to 121 x 121 pixels, and convergence was analyzed. This procedure was taken to reduce the computer processing time.

*2.1. Fractal dimension*

The fractal dimension calculations have been performed by the box-count method in the grayscale images^{17}, with the box-size ranging from 3 to 11. The fractal dimension based on the box-count method can be defined by:

where the linear fitting of D_{n} for many *n* values presents the fractal dimension (df) of an image divided *n* times and N_{n} is the number of boxes crossing the image.

The N_{n} calculation can be obtained based on grayscale values from pixels grid (i,j) as it follows:

The n_{n} values are obtained from the equation 5, where:

where Gray_{max}(i,j) and Gray_{min}(i,j) are the maximum and minimum grayscale image values in the grid (i,j) and S*'* is the size of the box.

In this way, N_{n} is taken for different values of n, that is, from different sizes of grids. This way of counting leads to a better approximation of the boxes that intercept the surface of the grayscale levels in the image. The df is obtained from the average value of D_{n}, or the angular coefficient of the linear fit from n vs. D_{n} dilog plot.

**3. Results and Discussion**

Figure 1 shows pictures of surfaces, in which each one of them was sprayed with %AIA solution. It was possible to observe a discrete evolution from hydrophobic to hydrophilic behavior, where the shape of the drops changes from circular to a water film, passing through irregular shapes, and the increase in the amount of the wetted surface.

In the pictures presented in the Figure 1, it is possible to observe the changing in the texture of the images. As the contact angle increases, the wetted area also increases. Another observed behavior was the decreasing in the circular factor of the drops with the addition of alcohol; this situation can be attributed to the low homogeneity of the surface^{18}.

Differences in the illumination could also be observed in the Figure 1. Based on this, EQU and/or WTH filter should be applied to reduce these illumination gradients. In order to compare these images with the STRI classification, the contact angles were measured as a function of the %AIA solution, thus, the HC was related to each one (see Figure 2), until the HC03. According to the STRI Guide, for an HC level above 04, the value must be defined using not only the contact angle, but also the amount of the wetted surface, which can be observed in the Figure 2. The mathematical fit was performed based on a third degree polynomial function, only for future comparisons between the contact angle and the image texture.

In order to define the opening image to be used as the WTH filter, the behavior of the volume of these images, i.e. the sum of all pixel values in the image, as a function of *B* (see Figure 3) was performed. The derivate curve begins to converge in *B* = 41 x 41 indicating that *B* values above 41 x 41 will not produce higher changes in the image. Thus, this value was adopted to reduce the processing time for each image.

As it follows, it can be observed in the Figure 4a, the opening image process obtained from the original Figure 1a with *B* = 41 x 41. As the image has one area clearer than the other, this part is more visible in this image. If the original image had the same illumination as the previous one, the opening image should have been completely dark. Figure 4b shows the final result obtained after the WTH operation, where it can be noticed that the variation in the illumination was significantly reduced, showing the efficiency of this method then.

Figure 5 presents the a) OR, b) EQU, c) WTH, d) EQU+WTH, e) EQU+WTH+EQU and f) SOB images, showing artificial illumination adjustments produced in each sample of image.

The fractal dimension of these images has been calculated based on Equation 3 and the results are plotted in the Figure 6. It can be observed that dispersion in the results, in general, are around 3% and all the curves tend to converge. As the main goal of this study is to predict the behavior of the hydrophobicity in polymeric surfaces, an injective function was sought in order to provide a reliable mathematical identity.

As a result, this injective function between %AIA image and *df* could only be obtained from the WTH images, with its equation being as it follows:

The Equation 6 relates the image texture from the patterns based on %AIA images and its fractal dimension value. In other words, this method provides a direct quantification of the hydrophobicity due to its texture characteristics. The main advantage of it is the objective classification promoted by this method, avoiding the subjective analysis of the operator.

**4. Conclusions**

A new method to relate the hydrophobicity classification according to STRI is proposed based on the digital image processing. Several samples of the image textures have been produced by spraying water-isopropyl alcohol solutions on polymer surfaces. In order to avoid the gradient of illumination of the pictures taken under natural conditions, different procedures as histogram equalization, white top-hat filter and edge detection have been tested. By using the fractal dimension to analyze the textures of the polymeric wetted surface images, it was possible to determine a mathematical function that can classify the hydrophobicity of the polymeric surfaces by digital image analyses.

**References**

1. Hackam R. Outdoor high voltage composite polymeric insulators. *Dielectrics and Electrical Insulation, IEEE Transactions* 1999; 6(5):557 585. [ Links ]

2. Kikuchi T, Nishimura S, Nagao M, Izumi K, Kubota Y, Sakata M. Survey on the Use on Non-Ceramic Composite Insulator. *Dielectrics and Electrical Insulation, IEEE Transactions* 1999; 6(5):548-556. [ Links ]

3. Sherif EM. *Performance and aging of HVAC and HVDC overhead line insulators.* [Ph. D Thesis]. Sweden: Chalmers Institute of Technology; 1987. [ Links ]

4. Altafim RAC, Silva JFR, Gonzaga DP, Ribeiro C, Godoy J, Basso HC, Bueno B, Calil-Jr C, Sartori JC, Altafim RAP, Silveira A. Wood cross-arms coated with polyurethane resin - tests and numerical simulations. *Materials Research*. 2006; 9(1):77-81. [ Links ]

*5. Hydrophobicity Classification Guide*. Guide 1, 92/1. Ludvika: STRI; 1992. [ Links ]

6. Stevens RT. *Fractal programming in C*. San Mateo: M&T Books; 1989. [ Links ]

7. Fisher Y, Jacobs, EW, Boss, RD. Fractal Image Compression Using Iterated Transforms. In: James A. *Strorer, Image and Text Compression.* Boston: Kluwer Academic Publishers; 1992. p. 35-61. [ Links ]

8. Peleg S, Naor J, Hartley R, Avnir D. Multiple Resolution Texture Analysis and Classification. *IEEE Transactions on Pattern Analysis and Machine Intelligence* 1984; PAMI(6):518-523. [ Links ]

9. Fahmy Y, Russ, JC, Koch CC. Application of Fractal Geometry Measurements to the Evaluation of Fracture Toughness of Brittle Intermetallics. *Journal of Materials Research* 1991; 6(9):1856-1861. [ Links ]

10. Srinivasan S, Russ, JC, Scattergood RO. Fractal Analysis of Erosion Surfaces. *Journal of Materials Research* 1990; 5(11):2616-2619. [ Links ]

11. Tokoro T, Nagao M, Kosaki M. Image Analysis of Hydrophobicity of silicone rubber insulator. In: *Annual Report Conference on Electrical Insulations and Dielectric Phenomena* (CEIDP); 1999; Piscataway, NJ, USA: IEEE, v.2, p. 763-766. [ Links ]

12. Tokoro T, Yanagihara S, Kosakii M. Diagnosis of Hydrophobic Condition of Silicone Rubber using Dielectric Measurement and Image Analysis. In: *Annual Report Conference on Electrical Insulation and Dielectric Phenomena*; 2005; USA, p. 281- 284. [ Links ]

13. Altafim RAC, Santana AMI, Murakami CRI, Basso HC, Chierice GO, Claro Neto S. Hydrophobicity of Polyurethane Resins. In: *International Conference on Solid Dielectrics;* 2004; Toulouse, France, p. 452-455. [ Links ]

14. Berg M, Thottappillil R, Scuka V. Hydrophobicity estimation of HV polymeric Insulating Materials - Development of a digital image processing method. *IEEE Transactions on Dielectric and Electrical Insulation*. 2001; 8(6):1098-2007. [ Links ]

15. Li Z, Liang X, Zhou YG, Tang J, Cui J, Lid Y. Influence of temperature on the hydrophobicity of silicone rubber. In: *Annual Report Conference on Electrical Insulation and Dielectric Phenomena;* 2004; Colorado, USA, p. 679-682. [ Links ]

16. Gonzalez RC, Woods RE. *Digital Image Processing*. 2 ed. New Jersey: Prentice-Hall; 2002. [ Links ]

17. Conci A, Proença CB. A Box Counting Approach for Color Segmentation. In: *International Conference on Image Processes - ICIP'97*; 1997 Oct 26-29; Santa Barbara, USA, p. 228-230. [ Links ]

18. Wu S. *Polymer Interface and Adhesion*. New York: Marcel Dekker, Inc.; 1982. [ Links ]

Received: February 8, 2008

Revised: 3/11/08

* e-mail: thomazini@unifor.br