Table of Contents Author Guidelines Submit a Manuscript
The Scientific World Journal
Volume 2014, Article ID 689496, 9 pages
Research Article

Video Texture Synthesis Based on Flow-Like Stylization Painting

Department of Computer Science and Engineering, School of Information Science and Engineering, Yunnan University, Kunming 650091, China

Received 19 March 2014; Accepted 7 June 2014; Published 15 July 2014

Academic Editor: Zhongmei Zhou

Copyright © 2014 Qian Wenhua et al. 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.


The paper presents an NP-video rendering system based on natural phenomena. It provides a simple nonphotorealistic video synthesis system in which user can obtain a flow-like stylization painting and infinite video scene. Firstly, based on anisotropic Kuwahara filtering in conjunction with line integral convolution, the phenomena video scene can be rendered to flow-like stylization painting. Secondly, the methods of frame division, patches synthesis, will be used to synthesize infinite playing video. According to selection examples from different natural video texture, our system can generate stylized of flow-like and infinite video scenes. The visual discontinuities between neighbor frames are decreased, and we also preserve feature and details of frames. This rendering system is easy and simple to implement.

1. Introduction

Most natural phenomena can use image or video to demonstrate directly. However, both natural and artificial phenomena cannot adequately be captured by a single static photo or image. On the other hand, though video is considered as the best media to show natural scene, it also has some shortcomings, such as if video scene should be stored on a computer or other storage devices, we always apply finite duration video clip. Therefore, video needs a beginning, a processing, and the ending frame. Further, though video captures the time-varying behavior of the phenomenon, it lacks the “timeless” quality of the photograph or image. In addition, it needs too much storage space.

In order to solve the above problems, Schödl and Essa proposed a new type of medium called “video texture” [1]. Video texture has qualities between photograph and video. This medium provides a continuous, infinitely varying stream of video images, and it can show an arbitrary length of video. The natural short video can be captured to synthesize infinite scene based on the technique of video texture synthesis. Then video texture works well for motions not only repetitive but also quasi-repetitive, such as swaying trees, waterfalls, and flickering flames. So, video texture is the best media to represent natural phenomena of our nature. Because of its wide application in the entertainment and industry and so forth, computer-based phenomena video texture becomes hot topic.

Generally, video texture may be regarded as expansion of image not only in the time domain but also in natural phenomena simulation. Schödl and Essa also employed video texture technique to animation. They used L2 distance to measure similarity between neighbor frames and maintained playing continuity among video frames. The proposed method has realized highly third dimension scenery animation successfully. Based on image cutting and synthesizing technique, Kuwahara improved video texture synthesis method. To ensure continuity, the input frames were divided into some parts, and then they would be synthesized to video texture [2]. Bhat et al. put forward an editing nature phenomena method to generate the final video texture. Their method first analyzed the motion of texture particles and segmented frames into patches, then they reorganized the patches to reduce error dither [3]. Their method can merge an arbitrary length video based on natural scene. Agarwala et al. described an automatic method for taking the output of a single panning video camera and creating a panoramic video texture [4]. Their method applied a dynamic programming step through hierarchical min-cut optimization process to improve synthesis speed.

Though there are so many video synthesis methods to synthesize natural video phenomena, the synthesis results are traditionally striven photorealism scenes. That is to say photorealism effects can not produce artistic sensation. Sometimes, people want to enjoy artistic video with cartoon, oil, and water painting effects. In addition, photorealistic images often include more details than necessary to communicate intended information. It is known that nonphotorealistic rendering (NPR) is the subject of intense debate for representative methods of artistic style creation. Artistic expression can often convey a specific mood which is difficult to imbue in a photorealistic scene. NPR also can focus viewer’s attention on the important information while downplaying extraneous or unimportant features. Additionally, a NPR look is often more engaging than the traditional, photorealism computer graphic rendering [5]. So, if video texture can be rendered to some nonphotorealistic effects, we can obtain artistic video with infinite playing scene.

Some NPR methods address the computer generated abstract stylization artistic effects, and this process fit perceptual task performance and evaluation. Decaudin applied modeling tools to generate abstract effect from natural scene [6], and his method kept details very well and shadow stayed coherent. Based on mean-shift, DeCarlo and Santella extended three-dimensional video volumes to automatically generate abstract artistic [7]. Particularly in the presence of occlusions and camera movement, contour tracking required substantial user correction of the segmentation results.

This paper presents an improved method to synthesize phenomena video with flow-like abstract painting. The new NP-video model can convert natural phenomena to a flow-like stylization video fast, and we also can achieve the infinite video texture easily. The amalgamations of the NPR and video texture synthesis can make the user enjoy artistic phenomena work, looking on the original phenomena in different artistic style, immersing vision and spirit into the artistic kingdom, and avoid some critical problems in photorealism systems.

2. NP-Video Phenomena Model System

We will introduce the basic architecture of our NP-video phenomena model in this section. As shown in Figure 1, our system proceeds in two steps. (1) Flow-like effects rendering: the natural phenomena video should be rendered to NPR artistic scene. (2) Video texture synthesis: the phenomena with flow-like effect will be synthesized to infinite playing video. We define the final play sequence of video texture to eliminate visual discontinuities. In addition, frames dividing and recombination method is used to arrive coherent playing frames. Our architecture supports a common NPR technique merging to video texture, and the final NP-video can give user flexibility and enjoyment.

Figure 1: NP-video phenomena of our system.

3. Flow-Like Artistic Style Simulating

Many arts create flow-like artistic work, and people always enjoy this work. Van Gogh and Munch emphasized flow-like structures and directional features in their paintings. The paintings are so famous because they are harmonic, interesting, and pleasant. Directional coherence and flow structure in this artistic work can help salience region features and boundaries. It also helps to evoke viewer’s imagination. So, before synthesizing video texture, we transfer the input phenomena to flow-like art first, and our method provides a versatile NPR flow-like rendering interface.

As same as the abstract effects, flow-like stylization artistic work can generate using edge-preserving filter method. Bilateral filter [8] and mean-shift method [9] are very famous examples of edge-preserving filter. Kang et al. improved the filter shapes, which is based on the vector field derived from the salient image features [10]. However, because edge-preserving filter cannot obtain good result with weak vector, this method may lose some details such as weak edge and structure. Kuwahara filter provides an overall stereoscopic painting visual and maintains a uniform level of abstraction across the image [11]. Unfortunately, Kuwahara filter can result in clustering artifacts during filter process. Papari et al. introduced a new weighting filter window, and the shape of the filter window can be changed based on local feature directions [12]. Kyprianidis also generated an abstract effect along the local feature directions [13], but the stylized surface of their methods cannot preserve both spatial and temporal coherence. Like bilateral filter and Kuwahara filter, Kang and Lee used Shock filter and mean curvature flow to generate abstraction effect, but the boundaries are not salient and simplified [14]. Recently, Kim et al. put forward bristle maps to generate aggregation, abstraction, and stylization of spatiotemporal data [15]. Isenberg explored visual map technique for generating of stylized renderings of 2D map data [16]. His method can render different stylization artistic work.

Though there are many techniques that can be used to obtain flow-like artistic work, these methods should be improved to adopt phenomena video. After capturing phenomena video, we will resolve this video and generate flow-like stylization painting of every frame. Our method begins with the anisotropic Kuwahara filtering. Firstly, eigenvalues and eigenvectors can be used to calculate structure tensor. Secondly, local orientation and anisotropy structure are used to guide Kuwahara filter process. Then the flow-like effects can be generated using improved line integral convolution.

3.1. Local Orientation Calculate

Brox et al. obtained anisotropy through calculating eigenvalues and eigenvectors of structure tensor [17]. We use Brox’s method to calculate structure tensor first. Let and be the spatial derivatives in - and -direction and let be the input frame with standard deviation ; structure tensor can be calculated as

The pixel’s gradient vector is . Let denote convolution operation; the partial derivatives of can be calculated:

The structure tensor of each video frame can be defined as

Let and be the minimum and the maximum eigenvalue of , which represents gray scale variation degree of vector and where parameter denotes filter kernel of Gauss function variance, denotes filter radius which is calculated using . If and are close, this filter field tends to be circular. The parameter represents edge direction of pixel , which reflects curve direction of eigenvalue. Then anisotropic structure tensor can be expressed as :

Because the local vector field maybe discontinuous, we use Gaussian filter technique to smooth this vector field, and eigenvalues of structure tensor are given:

Then, the eigenvector in local curve direction is calculated as

The local orientation is , and we calculate anisotropy based on Yang’s method [18]: where parameter ranges between 0 and 1, which denotes isotropic and anisotropic regions, respectively.

3.2. Anisotropic Filter

Because pixel of filter result is determined by the spatial distance from the center pixel , as well as its relative intensity difference, the filtering output at pixel is calculated: where parameter denotes filter spatial space. is the local intensity or color similarity with weighted , and denotes the squared standard deviations, and can be defined:

Let be the local orientation, be anisotropy, and be matrix defining a rotation, the filter spatial space is where parameter . Parameter is coefficient to adjust anisotropy . So, is a linear coordinate transform which maps to a disc with radius . This anisotropic filter ensures that if local region has lower standard deviation, the more weight value is given during filter process. Figure 2 shows our filter results. Figure 2(a) is input image, and Figures 2(b) and 2(d) are our results with different . If is smaller, more details will be retained. Figures 2(e) and 2(h) are filter results with other methods. Figure 3 also shows another result with blowing water. Figures 3(a) and 3(d) are different frame of natural texture, and Figures 3(e) and 3(h) are filter results. From these results, we can find that our method can maintain more details and features, such as edge, color, and local structure. Our method also eliminates artifacts of shape-simplifying and Shock filter method.

Figure 2: Comparison of different filter.
Figure 3: Video texture filter of blowing water.
3.3. Integral Convolution

We benefit from our efficient stylization painting not only in anisotropic filtering, but also with flow-like effect. Papari and Petkov applied integral convolution technique to produce nonphotorealistic rendering glass patterns texture effect [19]. We refer their integral convolution method to generate final artistic work. During convolution process, our method calculates a local stream line which moves out from the positive and negative directions. Let be flow-like output image, and let be a stream line with length of . Parameter denotes each pixel in . Let be current pixel; we can obtain . Parameter denotes convolution kernel, and we utilize Hanning window function in this paper. denotes all the pixels in this line, the convolution implement is defined as follows [20]: where parameter denotes convolution length. If is small, an insufficient amount of filtering flow occurs. However, too large will result in close convolution value and regular output. Because variance can well reflect local statistic information, it can be used to calculate local convolution length automatically: where and are the largest and the smallest integral length which can be assigned. Parameter is variance value of local area. Based on anisotropic filtering and integral convolution, input video can be rendered to flow-like effects. Next we will synthesize nonphotorealistic video texture to infinite video scene.

4. Phenomena Video Texture Synthesis

After phenomena video is rendered to flow-like effects, we will synthesize video texture clip to infinite video scene through video texture synthesis method. There are two significant challenges for video texture synthesis. (1) Because infinite video scene is generated through playing period video sequence repeat, the end frame in one sequence has large difference with the beginning frame in the next loop sequence. So, there are visual breaks between different loop arrays. (2) During video texture playing, visual discontinuities exist in different frames in one sequence. When frames dispose to nonphotorealistic artistic work, the discontinuities will be more prominent. So, we will determinate loop sequence of input video frames first. Then, methods of frame division and recombination, patches synthesis will be applied to eliminate visual discontinuities.

4.1. Loop Sequence Determinate

Many natural phenomena present playing cycle in video texture. If we play video cycle of natural phenomena repeat, people can enjoy infinite scene, and visual discontinuities will be reduced. So, we should find the loop cycle of a certain video texture. We use Schödl’s method of L2 distance to measure the similarity among different video frames [1].

Because human’s vision is sensitive to the change of luminance in frames, during the process of L2 distance calculation, we transfer the color space of frame from RGB to YCbCr. So L2 distance can be calculated only in luminance channel Y. If , denote two different frames in one video texture, similarity can be acquired as where , are luminance information of and . Parameters , represent pixel positions of and . So, if is smaller than the threshold , frames and can be regarded as the same loop sequence. That is to say, and are neighbor frames. If denotes frame number of input natural video texture and denotes each frame in ; threshold can be calculated:

Based on similarity, video texture can be divided into different sequences. Because L2 distance of neighbor frames in one sequence is small, we consider using different sequences to combine final infinite video scene. If L2 distance satisfies formula (16), then there only exists one sequence, and this sequence of input video texture has frames:

Figure 4 shows one natural video texture of stream water. There are fifteen frames to in this video, and we can obtain similarity in Table 1 using formula (14).

Table 1: L2 distance between the first frame and others.
Figure 4: Natural input video texture.

From Table 1, L2 distance satisfied . Then frames to belong to the one sequence, and fourteen is the cycle length of this video. If we want to obtain infinite video texture, frame should be played when frame is ended, and these fourteen frames should play repeatedly.

4.2. Frame Division and Recombination

Though L2 distance is smaller than threshold in one sequence, when playing video sequence repeatedly, there also exists visual discontinuity between neighbor frames, such as ending and beginning frames. Therefore, we utilize Bhat’s method to eliminate these discontinuities. Most natural phenomena video texture is similar particle moving, such as waterfall, fountain, flame, and stream water. Figure 5(a) shows a particle moving procedure along a flow line, and we can find that a patch texture of one frame in Figure 5(b) will shift to neighbor frame in next playing moment.

Figure 5: Natural phenomena moving like particle.

One frame in video texture can be divided into some texture patches, such as Figure 5(c). If the sequence of input video is four, then each frame should be divided into four texture parts, and each patch will shift along flow line like Figure 5(a). The patch will move to different position in neighbor frames. For example, when , the patch moves to the third part in the fourth frame, whose beginning is the first part in the second frame. If frames are divided into three parts, Figure 5(d) shows the result after frames recombined together, and the same color stands for the same texture patch. We can find that the same patch is in different frames when this video scene is playing.

Table 2 shows original playing sequence without frame dividing and recombination. The different column expresses different playing video frames. When this video texture is playing repeatedly, there is a visual discontinuity between the end frame of front cycle and the beginning frame of next cycle. If sequence length of this input video texture is three, then we divide each frame into three texture parts. After using frames recombination method, we obtain new group frames like Table 3. From Table 3, we can find that the discontinuity is scattered to different frames and forms a ladder-shaped discontinuity. Because this technique disperses visual dither to the whole playing sequence, visual break between different loop arrays will be reduced.

Table 2: Original frame loop sequence.
Table 3: Frame division and recombination.

The method of frames division and recombination is suitable to natural phenomena moving downwards, such as rivers, waterfalls, and fallen leaves. If input video texture moves upwards, such as fountains or flames, the frames should be divided into patches and recombination retroflex like Table 3. However, if a natural phenomenon is complicated and its sequence is more than ten, it is impossible to divide every frame into ten patches. Then, if video texture’s sequence is more than appointed threshold, each frame of this video can be divided using this threshold.

4.3. Patches Synthesis Method

Though visual discontinuity can be removed applying frame division and recombination, some distinct boundaries will be produced between different patches in one frame. Texture synthesis method can be used to eliminate these boundaries.

Texture synthesis method can be used to synthesize texture sample to an unlimited size image. Efros and Freeman proposed a texture quilting method to synthesize surprisingly good results for a wide range of textures [21]. Efros and Freeman chose some blocks from the sample texture to fill in the final texture image. The blocks satisfy the overlap constraints within some error tolerance. As shown in Figure 6, there have been some overlap regions between neighbor blocks.

Figure 6: Overlap regions between patches.

To eliminate discontinuity between overlap regions, based on the error cost path in overlap region, Efros found the best seam line to synthesize neighbor blocks. Hence, when frames of input video texture are divided into patches, we apply Efros’ method to recombine these patches. Let be error cost function and width of overlap region. denotes distance error of color, denotes structure error, and the minimum error boundary can be calculated as

Figure 7 shows the slide window when we calculate . Let be error cost of previous slider window and error cost of region A; is error cost of region B. Equation (18) is used to calculate the error cost between New and Old windows, and this method can optimize synthesis process between different patches:

Figure 7: The slide of matching window.

5. Experiment Results

To verify the feasibility and effectiveness of our proposed methods in this paper, we implemented the relevant algorithms and carried out our corresponding experiments. Figures 7 and 8 show the experimental results of video texture about blowing water and fountain scene. Our system works on Microsoft Visual C and the OpenGL texture mapping routines.

Figure 8: Video texture of blowing water.

Figure 8 shows blowing water video scene, and this input video has seventeen frames. Figure 9 shows fountain video scene, and this input video has fifty-eight frames. Figures 8(a) and 8(d) and Figures 9(a) and 9(d) are input frames from natural texture samples. The size of input frames is . Based on L2 distance calculation, we choose eight frames as loop sequence from blowing water and choose eighteen frames as loop sequence from fountain video scene. Figures 8(e) and 8(h) and Figures 9(e) and 9(h) show final frames with flow-like stylization effects based on our nonphotorealistic rendering and video texture synthesis methods. We can obtain the infinite video scenes playing loop sequence repeat, and this infinite video scenes has flow-like artistic effect. We also can conclude that the experimental results are realistic and they all have some esthetics.

Figure 9: Video texture of fountain.

6. Conclusions

In this paper, we have presented an NP-video texture synthesis system based on natural phenomena. Based on a simple and effective anisotropic filtering and integral convolution, flow-like stylization video effects can be generated. In addition, through loop sequence determination, frame division, and recombination, our algorithm removes visual discontinuities between different frames during video texture synthesis process. The experiment results indicate that our algorithm is easy to synthesize infinite video scene with nonphotorealistic artistic work.

An obvious limitation of our method is that the synthesis speed should be improved. Some strategies should be used to protect highly important area in the frames. We may further extend the GPU implementation of our algorithm to process video in real time.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.


This research was funded by the Grants (nos. 60663010, 61163019, and 61271361) from the Research Natural Science Foundation of China, (no. 2010CD024) from the Research Foundation of Yunnan Province, (no. 2012Z047) from the Research Foundation of the Educational Department of Yunnan Province, and (no. 20125301120008) from the Research Foundation of New Teacher Fund for Doctor Station, the Ministry of Education.


  1. A. Schödl and I. Essa, “Controlled animation of videosprites,” in Proceedings of ACM SIGGRAPH, pp. 121–127, San Antonio, Tex, USA, 2002.
  2. V. Kwatra, A. Schödl, I. Essa, G. Turk, and A. Bobick, “Graphcut textures: image and video synthesis using graph cuts,” ACM Transactions on Graphics (TOG)—Proceedings of ACM SIGGRAPH 2003, vol. 22, no. 3, pp. 277–286. View at Publisher · View at Google Scholar
  3. D. H. Bhat, S. M. Seitz, K. J. Hodgins et al., “Flow based video synthesis and editing,” Proceedings of ACM SIGGRAPH, vol. 23, no. 3, pp. 360–363, 2004. View at Google Scholar
  4. A. Agarwala, K. C. Zheng, C. Pal et al., “Panoramic video textures,” in Proceedings of the ACM SIGGRAPH (SIGGRAPH '05), pp. 821–827, August 2005. View at Publisher · View at Google Scholar · View at Scopus
  5. A. W. Klein, W. Li, M. M. Kazhdan, W. T. Corrêa, A. Finkelstein, and T. A. Funkhouser, “Non-photorealistic virtual environments,” in Proceedings of the 27th International Conference on Computer Graphics and Interactive Techniques (SIGGRAPH '00), pp. 527–534, July 2000. View at Scopus
  6. D. Philippe, Cartoon-looking rendering of 3D-scenes, [Ph.D. thesis], Universite de Techonologie de Compeigne, Inria, France, 1996.
  7. D. DeCarlo and A. Santella, “Stylization and abstraction of photographs,” in Proceedings of the ACM Transactions on Graphics (ACM SIGGRAPH '02), pp. 769–776, July 2002. View at Scopus
  8. C. Tomasi and R. Manduchi, “Bilateral filtering for gray and color images,” in Proceedings of the 6th International Conference on Computer Vision (ICCV '98), pp. 839–846, Bombay, India, January 1998. View at Scopus
  9. D. Comaniciu and P. Meer, “Mean shift: a robust approach toward feature space analysis,” IEEE Transactions on Pattern Analysis and Machine Intelligence, vol. 24, no. 5, pp. 603–619, 2002. View at Publisher · View at Google Scholar · View at Scopus
  10. H. Kang, S. Lee, and C. K. Chui, “Flow-based image abstraction,” IEEE Transactions on Visualization and Computer Graphics, vol. 15, no. 1, pp. 62–76, 2009. View at Publisher · View at Google Scholar · View at Scopus
  11. M. Kuwahara, K. Hachimura, S. Eiho, and M. Kinoshita, Digital Processing of Biomedical Images, Plenum Press, 1976.
  12. G. Papari, N. Petkov, and P. Campisi, “Artistic edge and corner enhancing smoothing,” IEEE Transactions on Image Processing, vol. 16, no. 10, pp. 2449–2462, 2007. View at Publisher · View at Google Scholar · View at MathSciNet · View at Scopus
  13. J. E. Kyprianidis, H. Kang, and J. Döllner, “Image and video abstraction by anisotropic Kuwahara filtering,” Computer Graphics Forum, vol. 28, no. 7, pp. 1955–1963, 2009. View at Publisher · View at Google Scholar · View at Scopus
  14. H. Kang and S. Lee, “Shape-simplifying image abstraction,” Computer Graphics Forum, vol. 27, no. 7, pp. 1773–1780, 2008. View at Publisher · View at Google Scholar · View at Scopus
  15. S. Kim, R. Maciejewski, A. Malik, Y. Jang, D. S. Ebert, and T. Isenberg, “Bristle maps: a multivariate abstraction technique for geovisualization,” IEEE Transactions on Visualization and Computer Graphics, vol. 19, no. 9, pp. 1438–1454, 2013. View at Publisher · View at Google Scholar · View at Scopus
  16. T. Isenberg, “Visual abstraction and stylisation of maps,” The Cartographic Journal, vol. 50, no. 1, pp. 8–18, 2013. View at Publisher · View at Google Scholar · View at Scopus
  17. T. Brox, J. Weickert, B. Burgeth, and P. Mrázek, “Nonlinear structure tensors,” Image and Vision Computing, vol. 24, no. 1, pp. 41–55, 2006. View at Publisher · View at Google Scholar · View at Scopus
  18. G. Z. Yang, P. Burger, D. N. Firmin, and S. R. Underwood, “Structure adaptive anisotropic image filtering,” Image and Vision Computing, vol. 14, no. 2, pp. 135–145, 1996. View at Publisher · View at Google Scholar · View at Scopus
  19. G. Papari and N. Petkov, “Continuous glass patterns for painterly rendering,” IEEE Transactions on Image Processing, vol. 18, no. 3, pp. 652–664, 2009. View at Publisher · View at Google Scholar · View at MathSciNet · View at Scopus
  20. B. Cabral and L. C. Leedom, “Imaging vector fields using line integral convolution,” in Proceedings of the 20th Annual Conference on Computer Graphics and Interactive Techniques (SIGGRAPH '93), pp. 263–270, August 1993. View at Scopus
  21. A. A. Efros and W. T. Freeman, “Image quilting for texture synthesis and transfer,” in Proceedings of the Computer Graphics Annual Conference (SIGGRAPH '01), pp. 341–346, Los Angeles, Calif, USA, August 2001. View at Scopus