Revision 15999 trunk/extensions/extRemoteSensing/src/org/gvsig/remotesensing/classification/ClassificationMaximumLikelihoodProcess.java

View differences:

ClassificationMaximumLikelihoodProcess.java
257 257
			Y = new Matrix(y);
258 258
			result= (Y.transpose().times(inverseVarCovMatrix[clase])).times(Y);
259 259
			// Obtencion probabilidad de pertenencia del pixel a la clase clase
260
//			probability= Math.log(detS[clase])+ result.get(0, 0);
261
			/*
262
			 * Me ahorro el c?lculo de log (no lo necesito para comparar). ????????????????????
263
			 */
264
			probability = detS[clase] + result.get(0, 0);
260
			probability= Math.log(detS[clase])+ result.get(0, 0);
265 261
			if(clase==0)
266 262
				finalProbability=probability;
267 263
			else if(probability<finalProbability){
......
323 319
			Y = new Matrix(y);
324 320
			result= (Y.transpose().times(inverseVarCovMatrix[clase])).times(Y);
325 321
			// Obtencion probabilidad de pertenencia del pixel a la clase clase
326
//			probability= Math.log(detS[clase])+ result.get(0, 0);
327
			/*
328
			 * Me ahorro el c?lculo de log (no lo necesito para comparar). ????????????????????
329
			 */
330
			probability = detS[clase] + result.get(0, 0);
322
			probability= Math.log(detS[clase])+ result.get(0, 0);
331 323
			if(clase==0)
332 324
				finalProbability=probability;
333 325
			else if(probability<finalProbability){
......
358 350
			Y = new Matrix(y);
359 351
			result= (Y.transpose().times(inverseVarCovMatrix[clase])).times(Y);
360 352
			// Obtencion probabilidad de pertenencia del pixel a la clase clase
361
//			probability= Math.log(detS[clase])+ result.get(0, 0);
362
			/*
363
			 * Me ahorro el c?lculo de log (no lo necesito para comparar). ????????????????????
364
			 */
365
			probability = detS[clase] + result.get(0, 0);
353
			probability= Math.log(detS[clase])+ result.get(0, 0);
366 354
			if(clase==0)
367 355
				finalProbability=probability;
368 356
			else if(probability<finalProbability){
......
393 381
			Y = new Matrix(y);
394 382
			result= (Y.transpose().times(inverseVarCovMatrix[clase])).times(Y);
395 383
			// Obtencion probabilidad de pertenencia del pixel a la clase clase
396
			//probability= Math.log(detS[clase])+ result.get(0, 0);
397
			/*
398
			 * Me ahorro el c?lculo de log (no lo necesito para comparar). ????????????????????
399
			 */
400
			probability = detS[clase] + result.get(0, 0);
384
			probability= Math.log(detS[clase])+ result.get(0, 0);
401 385
			if(clase==0)
402 386
				finalProbability=probability;
403 387
			else if(probability<finalProbability){

Also available in: Unified diff