| gvSIG desktop 1

  /**

   * Calculo de estad?sticas a partir de un histograma. El resultado de la funci?n es un array

   * bidimensional donde el primer ?ndice inndica la estadistica y el segundo el n?mero de banda.

   * 

   * <UL>

   * <LI>m?nimo</LI>

   * <LI>m?ximo</LI>

   * <LI>media</LI>

   * <LI>mediana</LI>

   * <LI>N?mero de pixels</LI>

   * </UL>

   * @param histogram

   * @param bandas solicitadas. Cada elemento del vector representa una banda. Si est? a true se calcula la 

   * estadistica para esa banda y si est? a false no se calcular?.

   * @return

   */

  public static long[][] getBasicStatsFromHistogram(int[][] histogram, boolean[] bands){

  	if(histogram == null)

  		return null;

  	return getBasicStatsFromHistogram(histogram, 0, histogram[0].length - 1, bands);

  }

  /**

   * Calculo de estad?sticas a partir de un histograma. El resultado de la funci?n es un array

   * bidimensional donde el primer ?ndice inndica la estadistica y el segundo el n?mero de banda.

   * 

   * <UL>

   * <LI>m?nimo</LI>

   * <LI>m?ximo</LI>

   * <LI>media</LI>

   * <LI>mediana</LI>

   * <LI>N?mero de pixels</LI>

   * </UL>

   * @param histogram

   * @param beginPos Posici?n de inicio del histograma para contabilizar estadisticas

   * @param endPos Posici?n de fin del histograma para contabilizar estadisticas

   * @param bandas solicitadas. Cada elemento del vector representa una banda. Si est? a true se calcula la 

   * estadistica para esa banda y si est? a false no se calcular?.

   * @return

   */

  public static long[][] getBasicStatsFromHistogram(int[][] histogram, int beginPos, int endPos, boolean[] bands){

  	if(histogram == null)

  		return null;

  	//Contamos el n?mero de bandas para las cuales se calcula la estad?stica

  	int bandCount = 0;

  	for(int iBand = 0; iBand < bands.length; iBand ++)

  		if(bands[iBand])

  			bandCount ++;

  	int values = 5;

  	long[][] res = new long[values][];

  	long[] min = new long[bandCount];//M?nimo

  	long[] max = new long[bandCount];//M?ximo

  	for(int iBand = 0; iBand < bandCount; iBand ++){

  		max[iBand] = beginPos;

  		min[iBand] = endPos;

  	}

  	long[] average = new long[bandCount]; //Valor de pixel medio (Media)

  	long[] middle = new long[bandCount]; //Mediana

  	long[] nPixelsBand = new long[bandCount];//N?mero de pixels por banda

  	int showBandCounter = 0;  //Contador de bandas de las que hay calcular la estadistica

  	for(int iBand = 0; iBand < histogram.length; iBand ++){

  		if(bands[iBand]){

    		int pixels = 0; //N?mero de valores por banda (entre 0 y 255)

			for(int i = beginPos; i <= endPos; i ++){

				//Calculo del m?nimo

				if(histogram[iBand][i] != 0 && i < min[showBandCounter])

					min[showBandCounter] = i;

				//Calculo del m?ximo										

				if(histogram[iBand][i] != 0 && i > max[showBandCounter])

					max[showBandCounter] = i;

				//Calculo del n?mero de pixeles

				nPixelsBand[showBandCounter] += (long)histogram[iBand][i];

				if(histogram[iBand][i] != 0)

					pixels ++;

				average[showBandCounter] += histogram[iBand][i] * i;

			}

			//Calculo de la media

			try{

				average[showBandCounter] /= nPixelsBand[showBandCounter];

			}catch(ArithmeticException exc){

				average[showBandCounter] = 0;

			}

			//Calculo de mediana

			long middlePos = nPixelsBand[showBandCounter] >> 1;

			int aux = 0;

			int i = beginPos;

			for(i = beginPos; aux < middlePos; i++)

				aux += histogram[iBand][i];

			middle[showBandCounter] = i - 1;

			showBandCounter ++;

  		}

  	}

  	res[0] = min;

  	res[1] = max;

  	res[2] = average;

  	res[3] = middle;

  	res[4] = nPixelsBand;

  	return res;	

  }