/trunk/extensions/extRemoteSensing/src/org/gvsig/remotesensing/georeferencing/geotransform/GeoTransformProcess.java - Diff - | gvSIG desktop 1 - gvSIG

Revision 18389 trunk/extensions/extRemoteSensing/src/org/gvsig/remotesensing/georeferencing/geotransform/GeoTransformProcess.java

     	// Lista con los valores de las potencias de x e y
     	private int exp[][] = null;
     	// error de la funcion en x
     	// vector con error de la funcion en x para cada punto
     	private double[] xError= null;
     	// error de la funcion para las y
     	//vector con el error de la funcion en y para cada punto
     	private double [] yError=null;
     	// rms
     	// rms vector con los rms para cada punto
     	private double rms[]=null;
     	// rms total en las x
     	private double rmsXTotal=0;
     	// rms total en las y
     	private double rmsYTotal=0;
     	// rms total
     	private double rmsTotal=0;
     	/** Metodo que recoge los parametros del proceso de clasificacion de
     	* m?xima probabilidad
     	* <LI>geo_points: lista de puntos con coordenadas reales</LI>
-...
     		getPolinomialCoefY();
     		// calculo de los resultados
     		getRMSerror();
+    	}
     	/**
     	 * @return coeficientes para el polinomio de transformaci?n en x.
     	 **/
     	public double[] getPolinomyalCoefX(){
     		if(coefX==null)
     			setDxGeneral();
     		return coefX;
+    	}
     		System.out.print(evaluate(1382,2893.5));
     	/**
     	 * @return coeficientes para el polinomio de transformaci?n en y.
     	 * */
     	public double[] getPolinomialCoefY(){
     		if(coefY==null)
     			setDyGeneral();
     		return coefY;
+    	}
     /////////!!!!!!! Obtiene matriz general a partir de los puntos de control y el orden del polinomio aproximador
     	public void setDxGeneral(){
         /**
          *  Calculo de los coeficientes del polinimio aproximador.
          *  @return array con los coeficientes para las x.
+         *
          * */
     	public double [] getPolinomyalCoefX(){
     		double matrixDx[][]= new double [minGPC][minGPC];
     		double result[]= new double[minGPC];
     		int exp[][]=new int[minGPC][2];
     		int k=-1;
     		// Obtencion de la primera fila de la matriz
     		int k=-1;
     		for (int filas=0; filas<minGPC;filas++)
     		{	k=-1;
     		for (int i=0; i<=orden; i++)
-...
     						result[k]+=(Math.pow(geo_points[v][0],i-j)* Math.pow(geo_points[v][0],exp[filas][0]))
     						* (Math.pow(geo_points[v][1],j)* Math.pow(geo_points[v][1],exp[filas][1]))*image_points[v][0];
+    				}
+    			}
+    		}
     		Matrix matrixResult= new Matrix(matrixDx);
     		coefX=solveSistem(matrixResult,result);
     		return coefX;
+    	}
     	// TO DO: Ver la manera de unificar con setDxGeneral(Parametrizar un metodo general)..... Estudiar optimizaciones!!!
     	// Cambios necesarios para el caolculo de coeficientes para coordenadas y
     	public void setDyGeneral(){
     	 /**
          *  Calculo de los coeficientes del polinimio aproximador.
          *  @return array con los coeficientes para las x.
+         *
          * */
     	public double[] getPolinomialCoefY(){
     		double matrixDy[][]= new double [minGPC][minGPC];
     		double result[]= new double[minGPC];
     		int k=-1;
     		// Obtencion de la primera fila de la matriz
     		for (int filas=0; filas<minGPC;filas++)
     		{	k=-1;
     		for (int i=0; i<=orden; i++)
-...
     						result[k]+=(Math.pow(geo_points[v][0],i-j)* Math.pow(geo_points[v][0],exp[filas][0]))
     						* (Math.pow(geo_points[v][1],j)* Math.pow(geo_points[v][1],exp[filas][1]))*image_points[v][1];
+    				}
+    			}
+    		}
     		Matrix matrixResult= new Matrix(matrixDy);
     		coefY=solveSistem(matrixResult,result);
     		return coefY;
+    	}
     	/**
-...
     		rms = new double [image_points.length];
     		xError= new double [image_points.length];
     		yError= new double[image_points.length];
     		double rmsxTotal=0;
     		double rmsyTotal=0;
     		double rmsTotal=0;
     		for(int im_point=0; im_point<image_points.length;im_point++){
     			for(int i=0; i<minGPC;i++)
-...
+    			}
     			xError[im_point]= Math.pow(xEvaluate[im_point] -image_points[im_point][0], 2);
     			rmsxTotal+= xError[im_point];
     			yError[im_point]= Math.pow(yEvaluate[im_point] -image_points[im_point][1], 2);
     			rmsyTotal+= yError[im_point];
     			rms[im_point]=Math.sqrt
+    			(
     					xError[im_point]+ yError[im_point]
     			);
     			rmsTotal+=rms[im_point];
+    		}
     		rmsTotal= rmsTotal/image_points.length;
     		rmsxTotal= rmsTotal/image_points.length;
     		rmsyTotal= rmsyTotal/image_points.length;
     		System.out.print("Base X\t\t");
     		System.out.print("Base Y\t\t");
     		System.out.print("WarpX\t\t");
-...
+    	}
     	/**
     	 *  @return error total para la coordenada X
     	 * */
     	public double getRMSx(){
     		return rmsXTotal;
+    	}
     	/**
     	 *  @return error total para la coordenada y
     	 * */
     	public double getRMSy(){
     		return rmsYTotal;
+    	}
     	/**
     	 *  @return error total para la coordenada y
     	 * */
     	public double getRMSTotal(){
     		return rmsTotal;
+    	}
     	/**
     	 *  Funci?n que evalua el polinomio de transformaci?n para un punto especifico
     	 * @param x coordenada x del punto
     	 * @param y coordenada y del punto

Also available in: Unified diff