%ffp
//-------------------------------------
// FM_TEMPLATE_1
//-------------------------------------
Category: "AstroPlugins"
Title: "Low End Smoother"
Copyright: "Copyright ©2005
Author: "James Ryan"
Organization: "Astroplugins"
Filename: "Low End Smoother"
Description: ""
Version: "v1.0"
URL: "www.grekalova.com"

About: "AstroPlugins !V\n" "!c\n!A\n!U"

Dialog: Text = "Astroplugins Artifical Flat !V  !c !A"
//-------------------------------------
	
ctl(0):"Red Cutoff",val=50
ctl(1):"Green Cutoff",val=50
ctl(2):"Blue Cutoff", val=50
ctl(4):"Amount",range = (1,30),val=10
ctl(5):"Segments",range = (4,35), val=20
ctl(6):"Red Range",range = (1,40), val=10
ctl(7):"Green Range",range = (1,40), val=10
ctl(8):"Blue Range",range = (1,40), val=10

//    ctl[6]: "count", val = 0, range=(0,255)

	
OnFilterStart:
{

//    Info("%d",imageMode);
 
	if( imageMode != 3) {
		ErrorOk("This filter only works on RGB 8bit images.");
		doAction(CA_CANCEL);
	}
	return false;
}//OnFilterStart
//-------------------------------------
// ForEveryTile
//-------------------------------------
ForEveryTile:
{
    int count, counter, reject, ave,ave1,sigmareject,sigmareject1;
    float sigma, sigmarun,fave1;
    int minval, lastmin, interval;
    int xval,yval, ix, iy, ixval,iyval, yseg,xseg;
    int x,y,r,g,b,r1,g1,b1,rctl,gctl,bctl,amount;
    int intensity, xbox, i,j,k,l,m;
    dword cubed;
    int rave, gave, bave;
    int rbar, gbar, bbar;



    rctl = ctl(6);
    gctl = ctl(7);
    bctl = ctl(8); 
    xbox = ctl(5);
    amount = ctl(4);

    sigmareject = ctl(5);
    sigmareject1 = (float) sigmareject/10;

    counter = 0;
    
    xseg = X/ctl(5);
    yseg = Y/ctl(5);
    xval = X/xseg;
    yval = Y/yseg;
    
    iy = 0;
    ix = 0; 


for (z=0; z<Z; Z++)
{
   if ( updateProgress(z,Z) ) break;
   for (y=y_start; y<y_end; y++)
     { 
       for (x=x_start; x<x_end; x++)
          { 		
            pset(x,y,0, 0 );         
	    }
	 }
}

	
for (z = 0; z<Z; z++)
  {
    for (iyval = 0; iyval < yseg + 1; iyval++)
      {
	  for (ixval = 0; ixval< xseg + 1; ixval++)
          {
               j = 0;
               ave1 = 0; 
               for (y= iyval * yval; y < iyval * yval + yval; y++)
                {
                  for  (x = ixval * xval; x < ixval * xval + xval; x++)
                    {
                       k = src(x,y,z);

                      if (k < ctl(z))
                        {
                          ave1 = ave1 + k;
                          j = j+1;
                        }
                    }
                }
      
                ave1  = ave1/j;
                j = 0;  
                sigmarun = 0;

                for (y= iyval * yval; y < iyval * yval + yval; y++)
              	{
                    for  (x = ixval * xval; x < ixval * xval + xval; x++)
                      {
                         k = src(x,y,z);
                        if (k < ctl(z))
                          {
                             sigmarun =(k-ave1) * (k-ave1) + sigmarun;
                             j = j+1;
                          }
                      }
                  }
                 if (j == 0) j = 1; 
                 sigma = (sqrt(sigmarun /j));

//if (counter < count){ Info("sigma  %f %f %d %d",sigma,sigmarun,j,ave1);counter = counter +1;}
          
                 for (y= iyval * yval; y < iyval * yval + yval; y++)
                	 {
                     for  (x = ixval * xval; x < ixval * xval + xval; x++)
                        {
                           minval = src(x,y,z);	

// if (y = Y && x = X){ Info("ave1  %d %d %d %d",x,y,minval,ave1);counter = counter +1;} 

                           if (minval  < ctl(z))
                             { 
                                amount = abs (minval - ave1);

            		        if (minval < ave1)
                                   {
                                      pset (x, y, z, 233);
                                   }else
                                   {
                                      pset (x, y, z, 166);
                                   }
                             }
                               
                        }
                   }
           }
      if (updateProgress(y,Y)) abort();
      } //end iyval
  }//end z
 	return true;
}//ForEveryTile