%ffp

// This template gets you started when you plan to enhance one of 
// your FilterFactory plugins with FilterMeister. For this purpose it is
// important to be put the FF source code into FilterMeister's 
// ForEverTile handler as it gives you a lot of new possibilities.



//------------------------------------------------------------------------------------------
// Filter Infos

Category   :"AstroPlugins"
Title      :"Star Repair"
Copyright  :"@2005"
Author     :"James Ryan"
Organization:"Astroplugins"
URL:"http://www.grekalova.com/plugins/"
Filename   :"Star Repair.8bf"
Description:"Star Repair"
Version    :"v1.01"
About      :"!T !V\n!D\n"
            "!c\n!U"
Dialog: Text = "AstroPlugins Star Repair !V  !c !A"


//Only apply filter to RGB and Grayscale images
SupportedModes: RGBMode




//------------------------------------------------------------------------------------------
// Filter Control Definitions


ctl(0):"Red",val=30
ctl(1):"Green",val=30
ctl(2):"Blue",val=30
ctl(3):"Radius",range = (0,32), val=8
ctl(4):"Bloat",range = (10,100), val=70
ctl(5):"Magnitude", val=100
ctl(6):"Blending",range = (0,100), val=70


//ctl(7):"Counter"
ctl(10):statictext,"Please move the mouse over the preview!",size=(150,*),pos=(300,100)
ctl(9):OWNERDRAW(staticedge),size=(20,20),pos=(300,10),
        range=(0, 0xffffff),
        color=fgColor,val=fgColor,
        action=preview,
        tooltip="Choose foreground color"
ctl(8):checkbox, "Pick Color",pos=(250,15)


//------------------------------------------------------------------------------------------
// Here comes the filter code



OnCtl(n):{ 

int i,j,r;

   setCtlColor(9, RGB(ctl(0),ctl(1),ctl(2)));
   if (n == 8  && e == FME_CLICKED)
     {
        j3 = 0;
        
     
     }
   if (n == CTL_PREVIEW && e == FME_MOUSEMOVE && j3==1 && ctl(8) == 0)
     {
       for (z=0;z<Z;z++)
         for (y=0;y<Y;y++)
           for (x=0;x<X;x++)
             pset(x,y,z, src(x,y,z));



  
       
       j1 = getPreviewCoordX();
       j2 = getPreviewCoordY();
      
       r = ctl(3)/2;
        
                  
               
       for (z=0;z<Z;z++)
         {
           for (i=j1-r;i<j1+r;i++)
             {
               pset(i,j2+r,z,255);
               pset(i,j2-r,z,255);
             } 
         }
           
       for (z=0;z<Z;z++)
         {
           for (j=j2-r;j<j2+r;j++)
             {
               pset(j1+r,j,z,255);
               pset(j1-r,j,z,255);
             } 
         }
       
     }     
        
       if (n==CTL_PREVIEW && e == FME_LEFTCLICKED_UP && ctl(8) == 0)
         { 
            j3=0;
            setCtlTextv (10, "Adjust Sliders");
            j4 = getPreviewCoordX();
            j5 = getPreviewCoordY();

            doAction(CA_PREVIEW);

          }
       if ( n==CTL_PREVIEW && e == FME_RIGHTCLICKED_DOWN && ctl(8) == 0 )
         { 
                             
             j3=1;
             setCtlTextv (10, "Move Cursor over star and leftclick");
          }

       if ( n==CTL_PREVIEW && e == FME_LEFTCLICKED_UP && ctl(8) == 1 )
         { 
             j1 = getPreviewCoordX();
             j2 = getPreviewCoordY();
                
             j3=1;
             setCtlColor(9, RGB(src(j1,j2,0),src(j1,j2,1),src(j1,j2,2)));
             setCtlVal(0, src(j1,j2,0));
             setCtlVal(1, src(j1,j2,1));
             setCtlVal(2, src(j1,j2,2));
             setCtlVal(8, 0);

          }

     updatePreview(n);
   
  
   
   return false;
}

OnFilterStart:
{
    setCtlColor(9, RGB(ctl(0),ctl(1),ctl(2)));
    setClickDrag(2);   // enable left mouse clicks
    setPreviewCursor(32515); // cross
    setCtlTextv (10, "Right Click in the Preview to Start");

    

    return false;
}

OnFilterEnd:
{
    return false;
}

ForEveryTile:{

	int r,g,b,count,counter;
    word cubed,xbox, i,j,k,l,m,n,radius;
    float sumx, sumy, avex, avey, posx, posy,totalr,totalg,totalb,sigmarun,sigmas;
    float mexp, e,r1,g1,b1, twopi,sigmainv,fx;
    float valuex,valuey, power,dx,dy,dxy,bloat,scalem, blend; 
    cubed = ctl(0) + ctl(1) + ctl(2);
    radius = ctl(3);
    counter = ctl(7);
    count = 0;
    bloat = (float)ctl(4)/100;
    scalem = (float)ctl(5)/50;
    blend = (float)ctl(6);

    e = 2.71828;
    twopi = 2.2675;
    

       
            posx = 0;
            posy = 0;
            sumx = 0;
            sumy = 0;
            totalr = 0;
            totalb = 0;
            totalg = 0;
            m = 0;
            n = y;
            m=0;
                

  //if (count < counter){Info("%f %f %f %d",sigmarun,sigmas,posx,xbox);count = count +1;} 
                 
                for (j=y_start; j<y_end;j++)
                  { 
                     for (i=x_start; i<x_end;i++)
                        {
                           if (bCircle(i,j,j4,j5,ctl(3))==1)
                              {

                              r = src(i,j,0);
                              g = src(i,j,1);
                              b = src(i,j,2);
                              sumx = sumx + (r + g + b) * i;
                              sumy = sumy + (r + g + b) * j;
                              totalr = totalr + r;
                              totalg = totalg + g;
                              totalb = totalb + b;
                                m = m + 1;
                           }                      
                        }
                    }
                  avex = totalr/m;
                  posx = sumx/(totalr + totalg + totalb);
                  posy = sumy/(totalr + totalg + totalb);
                  sigmarun = 0;
                    
                 for (j=y_start; j<y_end;j++)
                  {
                     for (i=x_start; i<x_end;i++)
                       {
                         r = src(i,j,0);
                         g = src(i,j,1);
                	     b = src(i,j,2);

                        if (bCircle(i,j,(int)posx,(int)posy,ctl(3))==1)
                           {
                                r = src(i,j,0);
                                sigmarun = (i - posx) * (i - posx) + sigmarun;

                            }
                        }
                    }
                     
                    sigmas = sqrt(sigmarun/m)*bloat;



             for (j=y_start; j<y_end;j++)
                  {
                     for (i=x_start; i<x_end;i++)
                       {
              
                        r = src(i,j,0);
	                    g = src(i,j,1);
                        b = src(i,j,2);

                         if (bCircle(i,j,(int)posx,(int)posy,ctl(3))==1)
                            {
                               pset(i,j,0, ctl(0) ); 
	                           pset(i,j,1, ctl(1)); 
                               pset(i,j,2, ctl(2)); 

                                dx = (i-posx);
                                dy = (j-posy);
                                dxy = sqrt(dx*dx+dy*dy);
                                dx = 1- (dxy*dxy*dxy/((float)radius*(float)radius*(float)radius));

if (count < counter){Info("t %f  %f  %f  %d r %d",blend,dx,dxy,radius,r);count = count +1;}
                                if (dx<0) dx=0;
                                blend = (float)ctl(6) * dx/100;
                               
                                if (dx<0) dx=0;
                                mexp = -dxy*dxy/(2 * sigmas*sigmas);
                                sigmainv = 1/(sigmas * twopi);
                                power = pow(e,mexp);
                                valuex = sigmainv * power;
                                r1 = valuex  * totalr * scalem/sqrt((float)m);
                                g1 = valuex  * totalg * scalem/sqrt((float)m);
                                b1 = valuex  * totalb * scalem/sqrt((float)m);
                                r1 = (blend * r1 + (1-blend) * src(i,j,0));
                                g1 = (blend * g1 + (1-blend) * src(i,j,1));
                                b1 = (blend * b1 + (1-blend) * src(i,j,2));


//if (count < counter){Info("t %f  %f  %f  %d r %d",blend,dx,dxy,radius,r);count = count +1;}
                              if (r1 > ctl(0))                               
                                pset(i,j,0, (int)r1); 
                             if (g1 > ctl(1))
                	            pset(i,j,1, (int)g1 ); 
                              if (b1 > ctl(2))
	                            pset(i,j,2, (int)b1); 
                             }
                            else
                             {
                               pset(i,j,0, r ); 
	                           pset(i,j,1, g ); 
                               pset(i,j,2, b); 
                         
                              }
                       }
                    }
                j1=(int)posx;
                j2=(int)posy;


	
	
	return true;	

}//ForEveryTile