分形之美

分形之美

看到有的网站搞分形艺术,才知道这个Z=Z*Z+C 如此的奇妙。
于是干脆在opencv下写出来,看看。效果还不错。呵呵。
这个小例子只能修改C的值,有兴趣的朋友可以自己改改代码,看能不能做出些美妙的图形。

代码提示:要获得高质量的图像,可以增加image_w,不过性能很低。

 

图片
 


参考链接
http://blog.csdn.net/spanzhang/
http://www.fractal.net.cn/

代码

 

//  test.cpp : 定义控制台应用程序的入口点。
//

#include 
" stdafx.h "
#include 
" cv.h "
#include 
" highgui.h "

const   int  image_w  =   100 ;
const   double  a  = 0.0 ;
const   double  b  = 0.0 ;
const   long  Nmax  =   255 ;   // Z值衰减(趋向于零) 的临界
const   double  Rmax  =   20 //  Z值没有界限地增加(趋向无穷)的临界

IplImage
*  image  = 0 ;
int  c_r_value  =   60 ;
int  c_i_value  =   40 ;

// 返回迭代的次数
int  getnum( int  i, int  j)
{
    
double  C_RealPart  =  ( double )(c_r_value - 50 ) / 100 * 2 * 2 ;   // -2 到 +2   // Julia 集合 常数C
     double  C_ImagPart  =  ( double )(c_i_value - 50 ) / 100 * 2 * 2 ;   //  -2 到 +2  // Julia 集合 常数C
     double  RealPart  =  ( double )i / image_w + a; 
    
double  ImagPart  =  ( double )j / image_w + b;
    
int  n = 0 ;
    
while (n < Nmax)
    {
        
// Z=Z*Z+C;
         double  ZM  =   RealPart * RealPart  +  ImagPart * ImagPart;
        
if (ZM > Rmax)  break ;
        
double  temp_RealPart  =  RealPart * RealPart  -  ImagPart * ImagPart + C_RealPart;
        
double  temp_ImagPart  =   2 * ImagPart * RealPart + C_ImagPart;
        RealPart
=  temp_RealPart;
        ImagPart
=  temp_ImagPart;  
        n
++ ;
    }
    
return  n;
}
void  onChangeimg( int  pos)
{
    cvZero(image);
    
for ( int  i  = 0 ;i < image_w;i ++ )
        
for ( int  j  = 0 ;j < image_w;j ++ )
        {
            
int  n1  =  getnum(i,j);
            
int  n2  =  getnum(j,i);
            
int  n3  =  getnum(j,i / 2 );
            
// printf("%d,",n);

            CvScalar sc1;
            
if (n1 > Nmax - 2 ) n1 = 0 ;
            
if (n2 > Nmax - 2 ) n2 = 0 ;
            
if (n3 > Nmax - 2 ) n3 = 0 ;

            sc1.val[
0 =  ( int )abs( 255 * ( float )n1 / Nmax);
            sc1.val[
1 =  ( int )abs( 255 * ( float )n2 / Nmax);
            sc1.val[
2 =  ( int )abs( 255 * ( float )n3 / Nmax);
            cvSet2D(image,i,j,sc1);

        }

        
double  nvalue  =  cvNorm(image, 0 ,CV_C);

        
for ( int  i  = 0 ;i < image_w;i ++ )
        
for ( int  j  = 0 ;j < image_w;j ++ )
        {
            CvScalar sc1;
            CvScalar sc2 
= cvGet2D(image,i,j);
            sc1.val[
0 =  ( int )abs( 255 * sc2.val[ 0 ] / nvalue);
            sc1.val[
1 =  ( int )abs( 255 * sc2.val[ 1 ] / nvalue);;
            sc1.val[
2 =  ( int )abs( 255 * sc2.val[ 2 ] / nvalue);;
            cvSet2D(image,i,j,sc1);
        }

}
int  main(  int  argc,  char **  argv )
{
   

    cvNamedWindow( 
" WIN " 0  );
    cvCreateTrackbar( 
" C real " " WIN " & c_r_value,  100 , onChangeimg );
    cvCreateTrackbar( 
" C imag " " WIN " & c_i_value,  100 , onChangeimg );
    
    image 
=  cvCreateImage( cvSize(image_w,image_w),  8 3  );
    onChangeimg(
0 );
    

        
for (;;)
        {
            
            cvShowImage(
" WIN " ,image);
            
if (cvWaitKey( 40 ) == 27 break ;
        }
        
return   0 ;

}

 

下载

opencv_share@163.com download
 

已标记关键词 清除标记
相关推荐
©️2020 CSDN 皮肤主题: 大白 设计师:CSDN官方博客 返回首页