#include<stdio.h>
#include<stdlib.h>
#include<math.h>
#include<time.h>


FILE *pf1 ,*pf2 ;
#define D 1447680 
#define  n  1040
#define  m  1392
#define star 3000

int i , j , k , d ,c ,b  ,num_star , num_starmaxi ;

double x[star],y[star];
float somme ,somme2 ,somme3 , moy ,carre_moy , ecart_type , seuil ,seuil2 ,max ;

float V[n*m];

int main(){
    
    num_starmaxi=2 ;
    somme3=7671;
    somme=0;
    somme2=0;
    seuil2=7670;
   pf1=fopen("allsky_20150716_232659.txt","r");
  // pf2=fopen("image_ciel.dat","w");
    
   for (i=0;i<n;i++){
  for(j=0;j<m;j++){
 
 
  fscanf(pf1,"%f \n",&V[i*m+j]);
  //fprintf(pf2,"%f \n",V[i*m+j]);
  
   // somme=somme+V[i*m+j];
   //moy=somme/D ;
   
  // printf("%f \n" , moy);
   
  
  
 
  //printf("%f \n",V[i*m+j]);
  
 
  }
   }//
  fclose(pf1);
 // fclose(pf2);
 
//pf1=fopen("image_ciel.dat","r");
//pf2=fopen("seuil_ciel.rtf","w");
    
//somme=0;
    
for(i=0;i<D ;i++){
  //fscanf(pf1,"%f \n ",&V[i]);

   somme=somme+V[i];
   somme2=somme2+V[i]*V[i];


 //printf("%f \n" , V[i]);
}
  moy= somme/D ;
  carre_moy= somme2/D;
  ecart_type= sqrt(carre_moy-(moy*moy));
  seuil= moy + 4*ecart_type ;
    
  printf("%f %f %f %f \n" , moy , carre_moy , ecart_type , seuil);
    printf("%f\n" , somme);
          for(i=0;i<D;i++){
        if(V[i]< seuil){
            V[i]=0.0 ;
//fprintf(pf2,"%f ∖n ",V[i]);
//printf("%f",V[i] );
  }
  }
//fclose(pf1);
//fclose(pf2);
    num_star=0;
   while(  num_star<num_starmaxi){
       
       
    for(i=0;i<n;i++){
        for(j=0;j<m;j++){
            k=i*m+j;

            somme3=somme3+V[k];
           
            
            if(max<V[k]){
                max=V[k];
                x[num_star]=j;
                y[num_star]=i;
            
            num_star=num_star+1;
            
              
                for(b=i-2;i<=i+2;i++){
                    for(c=j-2;c<=j+2;j++){
                        d=b*m+c;
                       V[d]=0;
                     }}
          
                }
            
                       
           
            
        }
           
        

               
            
        }
        printf("%f \n ",somme3);
    }
    
    
   
    for(i=0;i<num_starmaxi;i++){
    
        printf("%lf %lf  \n",x[i],y[i]);
    }
  return 0;
  }