// samuel bravo silva, santiago, chile. 2008. http://terreno.wordpress.com 
 
PImage ii; 
PFont font; 
 
float px,pr,pl; 
int i=1;//intervalo de análisis 
int pp=1; 
/*intervalo de puntos por rayo. si es igual al anterior solo veremos líneas, 
 pero eso ayuda a tener un análisis homogeneo..  
 este es el espesor mínimo de elementos detectables*/ 
int xb=0; 
int yb=0; 
int yy=0; 
int p=0; 
 
int ma, mi; //máximo y mínimo 
boolean dos; 
boolean ff; // saveframe 
boolean mm=true;//switch para calcular máximo y mínimo SÓLO UNA VEZ 
 
float  y;//coordenada del haz de rayos 
float rr=PI/20;//*2=cantidad de rayos 
 
void setup(){ 
  ii = loadImage("compostela.gif"); 
  size(ii.width,ii.height); 
  image(ii, 0, 0); 
  frameRate(300); 
  font = loadFont("DIN-Regular-14.vlw");  
} 
 
int trvx [][]= new int [1200][1200]; 
 
void draw(){ 
  smooth(); 
  if(key==' ')ff=true; 
  if(key=='p')ff=false; 
  if(ff){ 
    saveFrame("o4_####.png"); 
    ff=false;   
  } 
 
  if(yy<width){ 
    image(ii, 0, 0); 
    yy++; 
 
    for(int xx=0;xx<width;xx+=i){ 
      pr=brightness(get(xx,yy));  
      if(pr!=0){ 
        for (float r=-PI/4;r<PI/4;r+=rr){ 
          for (float x=pp;x<width/4;x+=pp){ 
            y=x*tan(r); 
            pl=brightness(get(int(x+xx),int(y+yy))); 
            if(pl==0)break; 
            //  point(x+xx,y+yy); 
            trvx[int(x+xx)][int(y+yy)]+=1; 
          } 
        } 
        for (float r=PI/4;r>-PI/4;r-=rr){ 
          for (float x=pp;x<width/4;x+=pp){ 
            y=x*tan(r); 
            pl=brightness(get(int(x+xx),int(y+yy))); 
            if(pl==0)break; 
            //  point(x+xx,y+yy); 
            trvx[int(x+xx)][int(y+yy)]+=1; 
          } 
        } 
        for (float r=-PI/4;r<PI/4;r+=rr){ 
          for (float x=pp;x<width/4;x+=pp){ 
            y=x*tan(r); 
            pl=brightness(get(int(-x+xx),int(-y+yy))); 
            if(pl==0)break; 
            //  point(-x+xx,-y+yy); 
            trvx[int(-x+xx)][int(-y+yy)]+=1; 
          } 
        } 
        for (float r=PI/4;r>-PI/4;r-=rr){ 
          for (float x=pp;x<width/4;x+=pp){ 
            y=x*tan(r); 
            pl=brightness(get(int(-x+xx),int(-y+yy))); 
            if(pl==0)break; 
            // point(-x+xx,-y+yy); 
            trvx[int(-x+xx)][int(-y+yy)]+=1; 
          } 
        } 
      }   
    }        
    progreso(yy); 
  }   
 
  if(yy>=width){  
    //situar rango 
    if(mm){  //hace que esto se lea solo una vez 
      for(int xx=0;xx<width;xx+=1){ 
        for(int yy=0;yy<height;yy+=1){ 
          if(ma<trvx[xx][yy])ma=trvx[xx][yy]; 
          if(mi>trvx[xx][yy]&&0<trvx[xx][yy])mi=trvx[xx][yy]; 
          mm=false; 
        } 
      } 
    } 
 
    //recolectar valores en cuadrados 
    if(yb<height){ 
      println(mi); 
      println(ma); 
      yb++;   
      loadPixels(); 
      for(xb=0;xb<width;xb++){ 
        int cc=0; 
        cc=250*(trvx[xb][yb]-mi)/ma; 
        color aa=color(0,0,0); 
        aa=color (cc+40,2*cc-200,130-cc); 
        //aa=color (cc,cc,cc); 
        if(cc==0) aa=color (0); 
        pixels[p]=aa; 
        p++; 
      } 
      updatePixels(); 
    } 
    if(yb==height) { 
      yb++;  
      escala(ma); 
    }   
  } 
} 
 
void escala(int ma){ 
  int i=ma/200; 
  println(i); 
  fill(180,120); 
  textFont(font);  
  text("0%", width-65, 36);  
  text("100%", width-82, 20+ma/i);  
  for (int a=0;a<ma;a+=i){ 
    float dd=a*200/ma; 
    stroke(dd+40,2*dd-200,130-dd,230); 
    line (width-40,20+a/i,width-20,20+a/i); 
  } 
  noFill(); 
  stroke(180,120); 
  quad(width-40,20,width-20,20,width-20,20+ma/i,width-40,20+ma/i);   
} 
 
void progreso(int cc){ 
  noStroke(); 
  fill(180,80); 
  quad(20,20,20+width/4,20,20+width/4,40,20,40); 
  noFill(); 
  stroke(180,120); 
  quad(20,20,20+cc/4,20,20+cc/4,40,20,40); 
  textFont(font);  
  fill(180,120); 
  text("sight cast", 20, 16);  
}