| 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 long[][] getBasicStatsFromHistogram(boolean[] bands){

	  	if(histogram == null)

	  		return null;

	  	return getBasicStatsFromHistogram(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 long[][] getBasicStatsFromHistogram(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;	

	  }