File:Algol-type variable binary star animation 4.gif

magick convert  -delay 10 -loop 0 *.png algol_type_binary_animation_1.gif

1. 3. generate light curve+img animation images from double star images
2. 1. python3 script
   3. 28.9.2023 0000.0000

import math import numpy as np import matplotlib.pyplot as plt import matplotlib.image as mpimg

1. import PIL

inputbase1=".\\pngs\\" outputbase1=".\\pngs2\\"

imagebase1="doublestar1"

animsize=72

animcoeff1=72/360

1. xx=np.linspace(0,359,360)
2. yy=np.sin(xx*math.pi/180.)*200+300

x0=2*np.array([00,30,60,95,125,160,180,220,250,280,300,330,360]) y0=300+100*(np.array([4,4,4,4,6,4,4,4,5,4,4,4,4])-4) x2 = np.linspace(0,360,360)*2 y2=np.interp(x2, x0, y0)

1. print(y2)

for n in range(1, 360,5): angel1=n nummero1=int(animcoeff1*angel1)+1 in_imagename1=inputbase1+imagebase1+str(nummero1).zfill(2)+".png" out_imagename1=outputbase1+imagebase1+str(nummero1).zfill(2)+".png" x3=angel1*2 y3=y2[int(angel1)] print(n, angel1,y3) img = mpimg.imread(in_imagename1) imgplot = plt.imshow(img) plt.scatter(x3,y3,color="red",s=400) plt.yticks(ticks=np.array([0,200,400,600]),labels=["0","3","6","9"]) plt.xticks(ticks=np.array([0,90,180,270,360])*2,labels=["0","90","180","270","360"]) plt.plot(x2,y2, color="yellow", lw=5) plt.savefig(out_imagename1) plt.show() #plt.scatter(x3,y3,color="black",s=400)

POV-Ray animation ini file

Antialias=Off Antialias_Threshold=0.1 Antialias_Depth=2

Input_File_Name="doublestar1.pov"

Initial_Frame=1 Final_Frame=72 Initial_Clock=0 Final_Clock=360

1. include "functions.inc"

1. declare taunolla=clock+90;

1. declare tau1=radians(taunolla);
2. declare tau2=radians(taunolla+180.0);

1. declare ra=1.2;
2. declare rb=1.2;
3. declare ax = sin(tau1)*ra+0.5;
4. declare ay = cos(tau1)*ra;

1. declare bx = sin(tau2)*rb-0.5;
2. declare by = cos(tau2)*rb;

1. declare Orange=<1,0.5,0>;
2. declare Yellow=<1,1,0>;

1. declare White=<1,1,1>;
2. declare Blue=<0.8,0.8,1>;
3. declare Black=<0.0,0.0,0.0>;
4. declare Grey=<0.5,0.5,0.5>;

1. declare atm_thickness1=0.1 ;

camera { location <0,10,-80> look_at 0 angle 4.5 }

1. macro convection_cells(c)

crackle pigment_map{ [0.0 rgb c*0.2] [0.2 rgb c] [0.9 rgb c] [1.0 rgb 1.0]

} warp{turbulence 0.3}

1. end

// yellowish star

#declare stara1= union {

  sphere
  { 
  0, 1

// translate <ax,0,ay>

  texture{
        pigment{
         // ripples
           granite scale 1/10
          // wrinkles
            scale 1 turbulence 0.3
          color_map {
                 [0.0 color <1,0.5,0> ]
               //  [0.1 color Grey   ]
               //  [0.2 color Yellow    ]
                 [1.0 color <1,1,0>   ]
                } // cmap
   
 
      } // pig
 
     finish {
     ambient 1
     diffuse 0.01
     reflection 0.01
      } // fin
  
  } // tex
    
    
    
  } // sphere

// halo 1

sphere

  { 
  
  0, 1 hollow no_shadow
  

    pigment {

    rgbt 1 }
    interior
    { media
      { emission 1/1
        density
        { 
            function {
          //    1*exp(-1*(sqrt(x*x+(y)*(y)+z*z)- 1 - 0.00001)/atm_thickness1)    
          // exp(-sqrt((x*x+y*y+z*z)-0) )*0.3
            pow( 2*sqrt(x*x+y*y+z*z),-6 )
           }
        
       density_map  
          
          { 
               [0 rgb 0.0001*<1,1,0>] 
               [0.01 rgb 0.003*<1,1,0>] 
               [.5 rgb <1,1>*1/2]
               [1 rgb <1,1,0>*1.5]
          }
       

      }
        samples 30,30 intervals 1 confidence .1 // quite fast settings
        method 3
      }
      
      
    }

       scale 1.75

  }  

}

// bluish star

1. declare stara2=union {

  sphere
  { 
  0, 1

 texture {
 pigment {
	crackle

pigment_map{ [0.0 rgb <1/100,1/100,1>*2]

[1.0 rgb <8/10,8/10,1>*2]

} warp{turbulence 0.75}

 }
 scale 0.1
 }
 
     finish {
     ambient 10
     diffuse 0.01
     reflection 0.01
      } // fin

     scale 1/2

  } // sphere

   sphere
   { 
  
  0, 1 hollow no_shadow

    pigment {

    rgbt 1 }
    interior
    { media
      { emission 1.5
        density
        {

              function {
            

// exp(-1.5*(sqrt(x*x+y*y+z*z)-0) )

  pow( 2*sqrt(x*x+y*y+z*z),-4 )

   }
   density_map
          { 
           [0 rgb 0.01*<0.5,0.5,1>] 
           [0.125 rgb 0.1*<0.5,0.5,1>] 
           [.5 rgb <0.5,0.5,1>*1/2]
           [1 rgb <0.5,0.5,1>*1.5]
          }

        }
        samples 30,30 intervals 2 confidence .001 // quite fast settings
        method 3
      }
    }
    
       scale 1/2
       scale 2

  }

}

// stars

object {stara1

       translate <ax,0,ay>

}

object {stara2
       translate <bx,0,by>

}