batch/includes/partenaires/classMMap.php

<?
if ( !function_exists('htmlspecialchars_decode') )
{
    function htmlspecialchars_decode($text)
    {
        return strtr($text, array_flip(get_html_translation_table(HTML_SPECIALCHARS)));
    }
}

include_once(FWK_PATH.'common/curl.php');

function supprDecimales($dec) {
	if ($dec>0 )
		return floor($dec);
	else 
		return ceil($dec);
}

function dec2dms($dec) {
	$d = supprDecimales($dec);
	$m = supprDecimales(($dec - $d) * 60);
	$s = abs(round(((($dec - $d) * 60) - $m) * 60));		
	$m = abs($m);
	return $d.'<27>'.$m."'".$s.'"';
}

/**/
function ALG0001($phi,$e) {
	$temp = ( 1 - ( $e * sin( $phi ) ) ) / ( 1 + ( $e * sin( $phi ) ) );
	$L = log ( tan ( (pi()/4) + ($phi/2) ) * pow ($temp, ($e/2) ));
	return $L;
}

/** Calcul de la latitude <EFBFBD> partir de la latitude isom<EFBFBD>trique
 **/
function ALG0002($L,$e,$epsilon) {
	$phi[0] = 2 * atan(exp($L)) - (pi()/2);

	$i=0;
	do {
		$i++;
		$temp = ( 1 + ( $e * sin( $phi[$i-1] ) ) ) / ( 1 - ( $e * sin( $phi[$i-1] ) ) );
		$phi[$i] = 2 * atan ( pow ($temp, ($e/2)) * exp ($L) ) - pi()/2;
	}
	while (abs($phi[$i] - $phi[$i - 1]) >= $epsilon);

	return $phi[$i];
}

/** Transformation de coordonn<EFBFBD>es en projection conique conforme de Lambert, en coordonn<EFBFBD>es g<EFBFBD>ographiques
 ** @param double $X		Coordonn<EFBFBD>e X en projection conique conforme de Lambert du point
 ** @param double $Y		Coordonn<EFBFBD>e Y en projection conique conforme de Lambert du point
 ** @param double $n		Exposant de la projection
 ** @param double $c		Constante de la projection
 ** @param double $Xs		Coordonn<EFBFBD>e Xs en projection du p<EFBFBD>le
 ** @param double $Ys		Coordonn<EFBFBD>e Ys en projection du p<EFBFBD>le
 ** @param double $lambdac	Longitude de l'origine par rapport au m<EFBFBD>ridien origine	
 ** @param double $e		Premi<EFBFBD>re excentricit<EFBFBD> de l'ellipso<EFBFBD>de
 ** @param double $epsilon	Tol<EFBFBD>rance de convergence
 ** @return array	lambda	Longitude par rapport au m<EFBFBD>ridien origine
 **					phi		Latitude
 **/
function ALG0004($X,$Y,$n,$c,$Xs,$Ys,$lambdac,$e,$epsilon) {
        $R = sqrt( pow(($X - $Xs),2) + pow(($Y - $Ys),2) );
        $gamma = atan(($X - $Xs)/($Ys - $Y));

        $lambda = $lambdac + ($gamma / $n);

        $L = (-1 / $n) * log(abs($R/$c));

        $phi = ALG0002($L,$e,$epsilon);

        $coords['lambda']=$lambda;
        $coords['phi']=$phi;

        return $coords;
}

/** Transformation des coordonn<EFBFBD>es g<EFBFBD>ographiques ellipso<EFBFBD>dales en coordonn<EFBFBD>es cart<EFBFBD>siennes
 **  
 **/
function ALG0009($lambda,$phi,$he,$a,$e) {
        $N = ALG0021($phi,$a,$e);

        $X = ($N + $he) * cos($phi) * cos($lambda);

        $Y = ($N + $he) * cos($phi) * sin($lambda);

        $Z = ($N * (1 - $e*$e)  + $he) * sin ($phi);

        $coords['X']=$X;
        $coords['Y']=$Y;
        $coords['Z']=$Z;

        return $coords;
}

/**/
function ALG0012($X,$Y,$Z,$a,$e,$epsilon) {
        $lambda = atan ($Y/$X);

        $P = sqrt($X*$X + $Y*$Y);
        $phi[0] = atan ($Z/ ($P * (1 - ( ($a*$e*$e)/sqrt($X*$X + $Y*$Y + $Z*$Z) ) ) ) );

        $i = 0;
        do
        {
                $i++;
                $temp =  pow((1 - ( $a * $e*$e * cos($phi[$i - 1] )/( $P * sqrt(1 - $e*$e * sin($phi[$i - 1])*sin($phi[$i - 1]))))),-1);
                $phi[$i] = atan( $temp * $Z / $P );
        }
        while (abs($phi[$i] - $phi[$i - 1]) >= $epsilon);

        $phix = $phi[$i];

        $he = ($P/cos($phix)) - ($a/sqrt(1 - $e*$e * sin($phix)*sin($phix)));

        $coords['lambda']=$lambda;
        $coords['phi']=$phix;
        $coords['he']=$he;

        return $coords;
}

/** Transformation d'un jeu de 7 param<EFBFBD>tres entre 2 syst<EFBFBD>mes g<EFBFBD>od<EFBFBD>siques
 ** @param	double $Tx	Translation suivant l'axe des X (de 1 vers 2)
 ** @param	double $Ty	Translation suivant l'axe des Y (de 1 vers 2)
 ** @param	double $Tz	Translation suivant l'axe des Z (de 1 vers 2)
 ** @param	double $D	Facteur d'<EFBFBD>chelle de 1 vers 2
 ** @param	double $Rx	Angle de rotation autour de l'axe des X en radian (de 1 vers 2)
 ** @param	double $Ry	Angle de rotation autour de l'axe des Y en radian (de 1 vers 2)
 ** @param	double $Rz	Angle de rotation autour de l'axe des Z en radian (de 1 vers 2)
 ** @param	double $U	Vecteur de coordonn<EFBFBD>es cart<EFBFBD>siennes tridimensionnelles dans le syst<EFBFBD>me 1 : U=(Ux,Uy,Uz)
 ** @return array Vecteur de coordonn<EFBFBD>es cart<EFBFBD>siennes tridimensionnelles dans le syst<EFBFBD>me 2, V=(Vx,Vy,Vz)
 */
function ALG0013($Tx,$Ty,$Tz,$D,$Rx,$Ry,$Rz,$U) {
	$V=array();
	$V['X'] = $Tx + $U['X'] * (1 + $D) + $U['Z'] * $Ry - $U['Y'] * $Rz;
	$V['Y'] = $Ty + $U['Y'] * (1 + $D) + $U['X'] * $Rz - $U['Z'] * $Rx;
	$V['Z'] = $Tz + $U['Z'] * (1 + $D) + $U['Y'] * $Rx - $U['X'] * $Ry;
	return $V;
}

/** D<EFBFBD>termination des param<EFBFBD>tres de calcul d'une projection Lambert conique
 ** dans le cas tangent, avec ou sans facteur d'<EFBFBD>chelle en fonction des param<EFBFBD>tres de d<EFBFBD>finition usuels
 ** @return array n, C, lambdac, Xs, Ys
 **/
function ALG0019($lambda0,$phi0,$k0,$X0,$Y0,$a,$e) {
        $lambdac = $lambda0;
        $n = sin($phi0);
        $C = $k0 * ALG0021($phi0,$a,$e) * tan (pi()/2 - $phi0) * exp ( $n * ALG0001($phi0,$e) );
        $Xs = $X0;
        $Ys = $Y0 + $k0 * ALG0021($phi0,$a,$e) * tan (pi()/2 - $phi0) ;

        $tab ['e'] = $e;
        $tab ['n'] = $n;
        $tab ['C'] = $C;
        $tab ['lambdac'] = $lambdac;
        $tab ['Xs'] = $Xs;
        $tab ['Ys'] = $Ys;

        return $tab;

}

/** Calcul de la grande normale de l'ellipso<EFBFBD>de
 ** @param  double $phi		Latitude
 ** @param  double $a		Demi-grand axe de l'ellipso<EFBFBD>de 
 ** @param  double $e		Premi<EFBFBD>re excentricit<EFBFBD> de l'ellipso<EFBFBD>de
 ** @return double
 **/
function ALG0021($phi,$a,$e) {
	$N = $a/sqrt( 1 - $e * $e * sin($phi) * sin($phi) );
	return $N;
}

/** Calcul des constantes d'une projection Lambert conique conforme dans le cas s<EFBFBD>cant
 ** @param	double $lambda0	Longitude origine en radian par rapport au m<EFBFBD>ridien origine	
 ** @param	double $phi0	Latitude origine	
 ** @param	double $X0		Coordonn<EFBFBD>e X en projection du point origine
 ** @param	double $Y0		Coordonn<EFBFBD>e Y en projection du point origine
 ** @param	double $phi1	Latitude en radian du 1er parall<EFBFBD>le autom<EFBFBD>co<EFBFBD>que	
 ** @param	double $phi2	Latitude en radian du 2<EFBFBD>me parall<EFBFBD>le autom<EFBFBD>co<EFBFBD>que
 ** @param  double $a		Demi-grand axe de l'ellipso<EFBFBD>de 
 ** @param	double $e		Premi<EFBFBD>re excentricit<EFBFBD> de l'ellipso<EFBFBD>de
 ** @return array n, C, lambdac, Xs, Ys
 **/
function ALG0054($lambda0,$phi0,$X0,$Y0,$phi1,$phi2,$a,$e) {
        $lambdac = $lambda0;
        $n = ( (log( (ALG0021($phi2,$a,$e)*cos($phi2))/(ALG0021($phi1,$a,$e)*cos($phi1)) ))/(ALG0001($phi1,$e) - ALG0001($phi2,$e) ));
        $C = ((ALG0021($phi1,$a,$e)* cos($phi1))/$n) * exp($n * ALG0001($phi1,$e));

        if ($phi0 == (pi()/2)) {
                $Xs = $X0;
                $Ys = $Y0;
        } else {
			$Xs = $X0;
			$Ys = $Y0 + $C * exp(-1 * $n * ALG0001($phi0,$e));
        }
		$tab=array(	'e' => $e,
        			'n' => $n,
        			'C' => $C,
        			'lambdac' => $lambdac,
        			'Xs'=> $Xs,
        			'Ys'=> $Ys);
        return $tab;
}

function Lambert2WGS84($X,$Y,$orig='L93') {
        $epsilon = 0.00000000001;

        $lambdac = 0.04079234433;	// pour greenwich
        $e = 0.08248325676; 		// premi<6D>re excentricit<69> de l ellipso<73>de Clarke 1880 fran<61>ais
        $he = 100;
        $a = 6378249.2;				// demi-grand axe de l'ellipsoide

        $Tx = -168;
        $Ty = -60;
        $Tz = +320;
        $D  = 0;
        $Rx = $Ry = $Rz = 0;

        $orig=strtoupper($orig);
        switch ($orig) {
			case 'LI':
			case 'L1':		$n = 0.7604059656;
							$c = 11603796.98;
							$Xs = 600000;
							$Ys = 5657616.674;
				break;
        	case 'LII':
        	case 'LIIE':
        	case 'L2E':
        	case 'L2':		$n = 0.7289686274;
							$c = 11745793.39;
							$Xs = 600000;
							if ($orig=='L2E' || $orig=='LIIE')
									$Ys = 8199695.768;
							else	$Ys = 6199695.768;
				break;
        	case 'LIII':
			case 'L3':		$n = 0.6959127966;
							$c = 11947992.52;
							$Xs = 600000;
							$Ys = 6791905.085;
				break;
			case 'LIV':
			case 'L4':		$n = 0.6712679322;
							$c = 12136281.99;
							$Xs = 234.358;
							$Ys = 7239161.542;
				break;
			case 'L93':
			default:		$n = 0.7256077650;
							$c = 11745255.426;
							$Xs = 700000;
							$Ys = 12655612.050;
				break;
        }

        $coords = ALG0004($X,$Y,$n,$c,$Xs,$Ys,$lambdac,$e,$epsilon);

        $coords = ALG0009($coords['lambda'],$coords['phi'],$he,$a,$e);

        $coords = ALG0013($Tx,$Ty,$Tz,$D,$Rx,$Ry,$Rz,$coords);

        $a = 6378137.0; // ellipso<73>des WGS84
        $f = 1/298.257223563;
        $b = $a*(1-$f);
        $e = sqrt(($a*$a - $b*$b)/($a*$a));

        $X = $coords['X'];
        $Y = $coords['Y'];
        $Z = $coords['Z'];
        $coords = ALG0012($X,$Y,$Z,$a,$e,$epsilon);

        $xy['long']  = rad2deg($coords['lambda']);
        $xy['lat']   = rad2deg($coords['phi']);
        return $xy;
}

  /**
  * Function to convert geographic coordinates to "lambert 2 etendue" coordinates
  * from: http://www.forumsig.org/showthread.php?p=64050#post64050
  **/
function geos2lambert($latitude,$longitude) {
    //0)degres-minutes-secondes + orientation (d,m,s,o) en radian
    $lambda_w = $longitude * pi()/180 ;
    $phi_w =  $latitude * pi()/180 ;

    //1) coordonn<6E>es g<>ographiques WGS84 (phi_w,lambda_w) en coordonn<6E>es cart<72>siennes WGS84 (X_w,Y_w,Z_w)
    $a_w = 6378137.0;
    $b_w = 6356752.314;

    $e2_w = ($a_w*$a_w-$b_w*$b_w)/($a_w*$a_w);

    //et on a la grande normale de l'ellipsoide WGS84
    $N = $a_w/sqrt(1-$e2_w*pow(sin($phi_w),2));

    //ainsi on obtient
    $X_w = $N * cos($phi_w) * cos($lambda_w);
    $Y_w = $N * cos($phi_w) * sin($lambda_w);
    $Z_w = $N * (1-$e2_w) * sin($phi_w);
    //Ref.: http://www.ign.fr/telechargement/MPr...RCE/NTG_80.pdf and http://de.wikipedia.org/wiki/WGS84

    //2) coordonn<6E>es cart<72>siennes WGS84 (X_w,Y_w,Z_w) en coordonn<6E>es cart<72>siennes NTF (X_n,Y_n,Z_n)

    $dX = 168.0;
    $dY = 60.0;
    $dZ = -320.0;

    $X_n = $X_w + $dX;
    $Y_n = $Y_w + $dY;
    $Z_n = $Z_w + $dZ;
    //ref.: http://support.esrifrance.fr/Documents/Generalites/Projections/Generalites/Generalites.htm#2

    //3) coordonn<6E>es cart<72>siennes NTF (X_n,Y_n,Z_n) en coordonn<6E>es g<>ographiques NTF (phi_n,lambda_n)
    $a_n = 6378249.2;
    $b_n = 6356515.0;

    $e2_n = ($a_n*$a_n-$b_n*$b_n)/($a_n*$a_n);
    //on d<>finit une tol<6F>rance de convergence
    $epsilon = pow(10,-10);

    //puis on amorce une boucle de calcul
    $p0=atan($Z_n/sqrt($X_n*$X_n+$Y_n*$Y_n)*(1-($a_n*$e2_n)/(sqrt($X_n*$X_n+$Y_n*$Y_n+$Z_n*$Z_n))));
    $p1=atan(($Z_n/sqrt($X_n*$X_n+$Y_n*$Y_n))/(1-($a_n*$e2_n*cos($p0))/(sqrt(($X_n*$X_n+$Y_n*$Y_n)*(1-$e2_n*pow(sin($p0),2))))));
    while(!(abs($p1-$p0)<$epsilon)){
        $p0 = $p1;
        $p1 = atan(($Z_n/sqrt($X_n*$X_n+$Y_n*$Y_n))/(1-($a_n*$e2_n*cos($p0))/(sqrt(($X_n*$X_n+$Y_n*$Y_n)*(1-$e2_n*pow(sin($p0),2))))));
    }
    $phi_n = $p1;
    $lambda_n = atan($Y_n/$X_n);

    //4) coordonn<6E>es g<>ographiques NTF (phi_n,lambda_n) en coordonn<6E>es projet<65>es en Lambert II <20>tendu (X_l2e, Y_l2e)
    $n = 0.7289686274;
    $c = 11745793.39;
    $Xs = 600000.0;
    $Ys = 8199695.768;

    $e_n = sqrt($e2_n);
    $lambda0 = 0.04079234433198;   //correspond <20> la longitude en radian de Paris (2<>20'14.025" E) par rapport <20> Greenwich
    //puis on calcule la latitude isom<6F>trique
    $L = log(tan(pi()/4 + $phi_n/2) * pow(((1-$e_n*sin($phi_n))/(1+$e_n*sin($phi_n))),($e_n/2)));

    //enfin on projette
    $X_l2e = $Xs + $c*exp((-$n*$L))*sin($n*($lambda_n-$lambda0));
    $Y_l2e = $Ys - $c*exp((-$n*$L))*cos($n*($lambda_n-$lambda0));

    return array("x_l2e"=>$X_l2e,"y_l2e"=>$Y_l2e);
}


/** Conversion WGS84 en Lambert93
 ** NTG_71.pdf
 ** http://www.ign.fr/affiche_rubrique.asp?rbr_id=1700&lng_id=FR 
 **/
function geos2lambert93($latitude,$longitude) {

	// Syst<73>me WGS84
	$a=6378137; 		// demi grand axe de l'ellipsoide (m)
	$e=0.08181919106;	// premi<6D>re excentricit<69> de l'ellipsoide

	// Param<61>tres de projections
    $l0=$lc=deg2rad(3);		// longitude de r<>f<EFBFBD>rence
    $phi0=deg2rad(46.5);	// latitude d'origine en radian
    $phi1=deg2rad(44);		// 1er parallele automécoïque
    $phi2=deg2rad(49);		// 2eme parallele automécoïque
    
    $x0=700000;		//coordonn<6E>es <20> l'origine
    $y0=6600000;	//coordonn<6E>es <20> l'origine
    
	// coordonn<6E>es du point <20> traduire
    $phi=deg2rad($latitude);
    $l=deg2rad($longitude);
    
 	// calcul des grandes normales
    $gN1=$a/sqrt(1-$e*$e*sin($phi1)*sin($phi1));
    $gN2=$a/sqrt(1-$e*$e*sin($phi2)*sin($phi2));
    
 	//calculs de slatitudes isom<6F>triques
    $gl1=log(tan(pi()/4+$phi1/2)*pow((1-$e*sin($phi1))/(1+$e*sin($phi1)),$e/2));
    $gl2=log(tan(pi()/4+$phi2/2)*pow((1-$e*sin($phi2))/(1+$e*sin($phi2)),$e/2));
    $gl0=log(tan(pi()/4+$phi0/2)*pow((1-$e*sin($phi0))/(1+$e*sin($phi0)),$e/2));
    $gl=log(tan(pi()/4+$phi/2)*pow((1-$e*sin($phi))/(1+$e*sin($phi)),$e/2));
    
	//calcul de l'exposant de la projection
    $n=(log(($gN2*cos($phi2))/($gN1*cos($phi1))))/($gl1-$gl2);//ok
    
	//calcul de la constante de projection
    $c=(($gN1*cos($phi1))/$n)*exp($n*$gl1);//ok
    
	//calcul des coordonn<6E>es
    $ys=$y0+$c*exp(-1*$n*$gl0);
    
	//calcul des coordonn<6E>es lambert
    $x93=$x0+$c*exp(-1*$n*$gl)*sin($n*($l-$lc));
    $y93=$ys-$c*exp(-1*$n*$gl)*cos($n*($l-$lc));
    
    return array("x_93"=>$x93,"y_93"=>$y93);   
}


//http://www.ign.fr/telechargement/MPro/geodesie/CIRCE/transfo.pdf
//http://www.ign.fr/telechargement/MPro/geodesie/CIRCE/NTG_80.pdf
//http://www.ign.fr/telechargement/MPro/geodesie/CIRCE/NTG_71.pdf
/**
print_r(Lambert2WGS84(591647.56,2426659.65));
print_r(geos2lambert(48.838245134184 ,2.2227849815878));
**/

/** The point-in-polygon algorythm allows you to programmatically check if a particular point is inside a polygon or outside of it. A common way to tackle the problem is to count how many times a line drawn from the point (in any direction) intersects with the polygon boundary. If they intersect an even number of times, then the point is outside.
	I used that approach in this PHP code, which doesn't contain detailed comments yet. A few people asked if I could post it anyway, so there it is... I'll try to add some more comment as soon as I have some spare time.

    The returned values are:
    "vertex" if the point sits exactly on a vertex AND you left true as the value for $pointOnVertex.
    "boundary" if the point sits on the boundary. If $pointOnVertex is false, then "boundary" is also returned if the point is on a vertex.
    "inside" if the point is inside the polygon.
    "outside" if, you guessed it, the point is outside of the polygon.
*/
class pointLocation {
	
	var $pointOnVertex = true; // Check if the point sits exactly on one of the vertices

	function pointLocation() {
	}

	function pointInPolygon($point, $polygon, $pointOnVertex = true) {
		$this->pointOnVertex = $pointOnVertex;

		// Transform string coordinates into arrays with x and y values
		$point = $this->pointStringToCoordinates($point);
		$vertices = array();
		foreach ($polygon as $vertex) {
			$vertices[] = $this->pointStringToCoordinates($vertex);
		}
//print_r($vertices);
//die();
		// Check if the point sits exactly on a vertex
		if ($this->pointOnVertex == true && $this->pointOnVertex($point, $vertices) == true) {
			return "vertex";
		}

		// Check if the point is inside the polygon or on the boundary
		$intersections = 0;
		$vertices_count = count($vertices);
		echo "vertices_count=$vertices_count".EOL;
		
		for ($i=1; $i < $vertices_count; $i++) {
			$vertex1 = $vertices[$i-1];
			$vertex2 = $vertices[$i];
			if ($vertex1['y'] == $vertex2['y'] && $vertex1['y'] == $point['y'] && $point['x'] > min($vertex1['x'], $vertex2['x']) && $point['x'] < max($vertex1['x'], $vertex2['x'])) {
				// Check if point is on an horizontal polygon boundary
				return "boundary";
			}
			if ($point['y'] > min($vertex1['y'], $vertex2['y']) && $point['y'] <= max($vertex1['y'], $vertex2['y']) && $point['x'] <= max($vertex1['x'], $vertex2['x']) && $vertex1['y'] != $vertex2['y']) {
				$xinters = ($point['y'] - $vertex1['y']) * ($vertex2['x'] - $vertex1['x']) / ($vertex2['y'] - $vertex1['y']) + $vertex1['x'];
				if ($xinters == $point['x']) {
					// Check if point is on the polygon boundary (other than horizontal)
					return "boundary";
				}
				if ($vertex1['x'] == $vertex2['x'] || $point['x'] <= $xinters) {
					$intersections++;
				}
			}
		}
		// If the number of edges we passed through is even, then it's in the polygon.
		echo "intersections=$intersections".EOL;
		if ($intersections % 2 != 0) {
			return "inside";
		} else {
			return "outside";
		}
	}

	function pointOnVertex($point, $vertices) {
		foreach($vertices as $vertex) {
			if ($point == $vertex) {
				return true;
			}
		}
	}
	
	function pointStringToCoordinates($pointString) {
		$coordinates = explode(' ', $pointString);
		return array('x'=>$coordinates[0], 'y'=>$coordinates[1]);
	}

}

class MMap {

	private $accesDist = true;
	private $referer ='';
	private $body = '';
	private $header = '';
	private $cookie = '';
	private $codeRetour = 0;
	public $codeRetourG = 0;	// Code retour de l'api Google
	public $codeRetourY = 0;	// Code retour de l'api Yahoo
	public $codeRetourM = 0;	// Code retour de l'api Mappy
	public $codeRetourC = 0;	// Code retour de l'api Cloudmade
	
	public $geocodeur='';
	
	public $latitudeDec= 0;		// Latitude en D<>cimal
	public $latitudeDeg= 0;		// Latitude en D<>gr<67>s
	public $longitudeDec= 0;	// Longitude en D<>cimal
	public $longitudeDeg= 0;	// Longitude en D<>gr<67>s
	public $altitude = NULL;	// Altitude en m<>tres
	public $l93x = NULL;		// Lambert 93, coordonn<6E>e X
	public $l93y = NULL;		// Lambert 93, coordonn<6E>e Y

	public $precision = 0;
	public $adresseValidee='';
	public $enCache=false;

//	private $apiKeyG='ABQIAAAAuKBtUyFonYJBl1fqfc78tRQvADPcxwXf3Q2QIE-M32vuSkrxiBRLUHDB_YSLeTscTDeWRKM_wILaaw';
	private $apiKeyG='AIzaSyBfyIMaVwm1nn8AupBf-asC816R71LZl8g';
//	private $apiKey='ABQIAAAAuKBtUyFonYJBl1fqfc78tRT4Qe_2x2fax4Crd4sBzQI-tDkl3BRm2l2feTV5Lnx0Ah-aTIlwHtglpA';
	private $apiKeyY='klII0gzV34Ho2TCsN3XiXeh4No033AifxMeDXfFXj8TDCXO3WAtOLj9c74VsV45HcII-';
	private $apiKeyC='ae5216123e414fcd84b8ca9a0f43a68c';

	function __construct($accesDist=true) {// $adresse, $cp, $ville, $pays='France') {
		$this->accesDist=$accesDist;
	}

/*
http://platform.beta.mapquest.com/geocoding/v1/address?key=YOUR_KEY_HERE&location=Lancaster,PA&callback=renderGeocode
http://platform.beta.mapquest.com/geocoding/v1/address?key=YOUR_KEY_HERE&callback=renderOptions&inFormat=kvp&outFormat=json&location=77 rue lost al lann, 29000 quimper,france

Client ID:		83955 
Email Address: 		buzuk77@aim.com
Registry Password:	tU2cA4qV 
Forums User Name: 
JS/AS3 Key:		Fmjtd%7Clu61nuuzlu%2Cr0%3Do5-5a1sd 
*/

/*
GeoCodingHandler.ashx?addr=
<wsdl:import namespace="http://dev.virtualearth.net/webservices/v1/imagery/contracts" location="http://dev.virtualearth.net/webservices/v1/metadata/imageryservice/imageryservice1.wsdl" />
API keybing = 56D6CBA671C986D3EA11B1B48F97507BE9B46999
*/	
	public function geoCodeAdrCpVilleM($adresse,$cp,$ville,$pays='') {

		$this->geocodeur='Mappy';
		$this->precision=0;
		
		$fileName=REP_TEMP.'geoCodeur'.$this->geocodeur;
		$tmp=@explode('|',file_get_contents($fileName));
		$nbQuery=@$tmp[1]*1;
		
		if ($pays=='')	$pays='France';
		if (preg_match('/^([0-9]{1,4} )(.*)$/', $adresse, $matches)) {
			$adresseVoie=trim($matches[2]);
			$adresseNum=trim($matches[1]);
		} else {
			$adresseVoie=trim($adresse);
			$adresseNum='';
		}
		$url='http://axe.mappy.com/services/loc/get.aspx?countryName='.urlencode($pays).'&townName='.urlencode($ville).'&wayName='.urlencode($adresseVoie).'&number='.urlencode($adresseNum).'&namedPlaceSearch=0&interactive=1&auth=maiaTest&format=json';
		
		$page=getUrl($url, '', '', '', false, '', '', 3);
		
		$res=json_decode($page['body']);
		$result=$res->geocode->answer;
		$nbQuery++;
		file_put_contents($fileName, $this->geocodeur.'|'.$nbQuery.'|'.$this->codeRetour."|$url|".date('Y-m-d H:i:s').'|'.$this->body.EOL);
		
		
		print_r($page);
		print_r($result);
		die($url);
		
		switch ($result->geocode_level->code) {
			case 0:		$this->precision=0;	break;	// No geocoding
			case 1:		$this->precision=1;	break;	//  Country level geocoding
			case 2:		$this->precision=4;	break;	//  Town level geocoding
			case 3:		$this->precision=6;	break;	//  Way level geocoding
			case 4:		$this->precision=8;	break;	//  Road element level geocoding
			default:	die('Precisions Inconnue : "'.$result->geocode_level->code.'"='.$result->geocode_level->label.EOL.$page['body']);	break;
		}
		
		$this->adresseValidee=@preg_replace('/ +/', ' ', strtoupper(str_replace('&apos;', "'", trim(htmlspecialchars_decode($result->road_element->number.' '.$result->way->name,ENT_QUOTES)))));
		
		$this->codeRetourM=@$result->global_score;
/*
        $result->postal_code->postal_code
        $result->subcountry->iso_code	// 11
        $result->subcountry->language
		$result->subcountry->Ile-de-France

		$result->town->language] => FRE
		$result->town->name] => Nanterre
        $result->town->official_code] => 92050
        $result->town->postal_codes->postal_code] => 92000           

0 	 Unknown location. (Since 2.59)
1 	Country level accuracy. (Since 2.59)
2 	Region (state, province, prefecture, etc.) level accuracy. (Since 2.59)
3 	Sub-region (county, municipality, etc.) level accuracy. (Since 2.59)
4 	Town (city, village) level accuracy. (Since 2.59)
5 	Post code (zip code) level accuracy. (Since 2.59)
6 	Street level accuracy. (Since 2.59)
7 	Intersection level accuracy. (Since 2.59)
8 	Address level accuracy. (Since 2.59)
*/
		$this->latitudeDec=@$result->location->coordinates->y;
		$this->longitudeDec=@$result->location->coordinates->x;
		$this->latitudeDeg=dec2dms($this->latitudeDec);
		$this->longitudeDeg=dec2dms($this->longitudeDec);

		//if ($this->codeRetourM>0)
			return true;

		return false;
	}
	
	public function geoCodeAdrCpVilleC($adresse,$cp,$ville,$pays='') {
	
		$this->geocodeur='Cloudmade';
		$this->precision=0;
		
		$fileName=REP_TEMP.'geoCodeur'.$this->geocodeur;
		$tmp=@explode('|',file_get_contents($fileName));
		$nbQuery=@$tmp[1]*1;
		
		if ($pays=='')	$pays='France';
		/*if (preg_match('/^([0-9]{1,4} )(.*)$/', $adresse, $matches)) {
			$adresseVoie=trim($matches[2]);
			$adresseNum=trim($matches[1]);
		} else {
			$adresseVoie=trim($adresse);
			$adresseNum='';
		}*/
		$url='http://beta.geocoding.cloudmade.com/v3/'.$this->apiKeyC.'/api/geo.location.search.2?format=json&source=OSM&enc=UTF-8&limit=10&locale=fr&q='.urlencode("$adresse $cp $ville $pays");
		$page=getUrl($url, '', '', '', false, '', '', 15);
		$this->body=$page['body'];
		$this->codeRetour=$page['code'];
		$this->header=$page['header'];
		
		$nbQuery++;
		file_put_contents($fileName, $this->geocodeur.'|'.$nbQuery.'|'.$this->codeRetour."|$url|".date('Y-m-d H:i:s').'|'.$this->body.EOL);
				
		$res=@json_decode(utf8_decode($page['body']));
		if ($this->codeRetour==200 && $res->status->success==1) {
			$geoCodage=$res->places[0];
			$this->codeRetourC=1;
			switch($geoCodage->addressType) {
				case 'housenumber':		$this->precision=8;				break;
				case 'street':			$this->precision=6;				break;
				case 'zip':				$this->precision=5;				break;
				case 'city':			$this->precision=4;				break;
				case 'district':		$this->precision=3;				break;
				case 'sight':			$this->precision=0;				break;
				default:
					print_r($res);
					die('addressType inconnue : '.$geoCodage->addressType.EOL);	
				break;
			}
		
			$this->adresseValidee=@trim(strtoupper(preg_replace('/\s+/',' ',$geoCodage->houseNumber.' '.$geoCodage->street.', '.
								$geoCodage->zip.' '.$geoCodage->city.', '.str_replace(',',' ',$geoCodage->country))));
				
			$this->latitudeDec=$geoCodage->position->lat;
			$this->longitudeDec=$geoCodage->position->lon;
			$this->latitudeDeg=dec2dms($this->latitudeDec);
			$this->longitudeDeg=dec2dms($this->longitudeDec);
		
			/*print_r($res);
			echo $this->adresseValidee.EOL;
			echo 'lat='.$this->latitudeDec.EOL;
			echo 'lon='.$this->longitudeDec.EOL;
			echo 'precis='.$this->precision.EOL;
			die($url);*/
			if ($this->precision>0)
				return true;
		}
		/*print_r($page);
		print_r($res);
		die($url);*/
		return false;
	}
	
	private function geoCodeAdrCpVilleY($address) {
		$this->geocodeur='Yahoo';
		$this->precision=0;
		
		$fileName=REP_TEMP.'geoCodeur'.$this->geocodeur;
		$tmp=@explode('|',file_get_contents($fileName));
		$nbQuery=@$tmp[1]*1;
		
		$url='http://query.yahooapis.com/v1/public/yql?q=select+*+from+geo.placefinder+where+text%3D%22'.rawurlencode($address).'%22&format=xml';
		$_result = false;
		$page=getUrl($url, '', '', '', false, '', '', 3);
		$this->body=$page['body'];
		$this->codeRetour=$page['code'];
		$this->header=$page['header'];
		
		$nbQuery++;
		file_put_contents($fileName, $this->geocodeur.'|'.$nbQuery.'|'.$this->codeRetour."|$url|".date('Y-m-d H:i:s').'|'.$this->body.EOL);
				
		if ($this->codeRetour==200) {
			
			/*R<EFBFBD>ponse Yahoo = <?xml version="1.0" encoding="UTF-8"?>
			<query xmlns:yahoo="http://www.yahooapis.com/v1/base.rng" yahoo:count="1" yahoo:created="2013-09-25T13:52:43Z" yahoo:lang="en-US"><results><Result><quality>71</quality>
			<latitude>45.782864</latitude><longitude>4.828002</longitude><offsetlat>45.782864</offsetlat><offsetlon>4.828002</offsetlon><radius>400</radius>
			<name/><line1>Rue Jacques Louis-Henon</line1><line2>69004 Lyon</line2><line3/><line4>France</line4><house/><street>Rue Jacques Louis-Henon</street><xstreet/>
			<unittype/><unit/><postal>69004</postal><neighborhood>4E Arrondissement</neighborhood><city>Lyon</city><county>Rhône</county><state>Rhone-Alpes</state>
			<country>France</country><countrycode>FR</countrycode><statecode/><countycode/><uzip>69004</uzip><hash/><woeid>20068131</woeid><woetype>11</woetype></Result></results>
			</query><!-- total: 90 -->
			<!-- engine8.yql.bf1.yahoo.com -->
			
			<query xmlns:yahoo="http://www.yahooapis.com/v1/base.rng" yahoo:count="1" yahoo:created="2014-05-20T18:13:30Z" yahoo:lang="en-US"><results><Result><quality>87</quality>
			<addressMatchType>POINT_ADDRESS</addressMatchType>
			<latitude>48.281769</latitude><longitude>4.09041</longitude><offsetlat>48.281731</offsetlat><offsetlon>4.09047</offsetlon><radius>400</radius>
			<name/><line1>38B, Rue Maurice Romagon</line1><line2>10800 Saint-Julien-les-Villas</line2><line3/><line4>France</line4><house>38B</house><street>Rue Maurice Romagon</street><xstreet/>
			<unittype/><unit/><postal>10800</postal><neighborhood/><city>Saint-Julien-les-Villas</city><county>Aube</county><state>Champagne-Ardenne</state><country>France</country><countrycode>FR</countrycode><statecode/><countycode>10</countycode><uzip>10800</uzip><hash>D49A3CF6AD6DEB1C</hash><woeid>12722975</woeid><woetype>11</woetype></Result></results></query><!-- total: 49 -->
<!-- engine4.yql.bf1.yahoo.com -->

			*/
			
			//echo "Yahoo :'$address'".EOL.$this->body.EOL;
			if (preg_match('! yahoo:count="0" !U', $this->body, $_match)) {
				$this->codeRetourY='ZERO_RESULTS';
				//echo $this->codeRetourY.EOL;
				return false;
			} elseif (preg_match('!<quality>(.*)</quality>(?:.*)<latitude>(.*)</latitude><longitude>(.*)</longitude>!U', $this->body, $_match)) {
				if (preg_match('!<quality>(.*)</quality><addressMatchType>(.*)</addressMatchType><latitude>(.*)</latitude><longitude>(.*)</longitude>!U', $this->body, $_match2)) {
					$quality=$_match2[1];
					$qualityTxt=$_match2[2];
					$this->codeRetourY='OK';
					if     ($quality>87 || $qualityTxt=='POINT_ADDRESS')	$this->precision=9;	// POI or Coordinate
					elseif ($quality>82 || $qualityTxt=='INTERPOLATED')		$this->precision=8; // Address (N<> compris dans une Tranche d'adresse)
					elseif ($quality>75)	$this->precision=-7;	// Intersection
					elseif ($quality>70)	$this->precision=-6;	// Street match
					elseif ($quality>40)	$this->precision=-5; // Zip Code
					elseif ($quality>30)	$this->precision=-4; // Town
					elseif ($quality>20)	$this->precision=-3; // Sous/R<>gion
					elseif ($quality>10)	$this->precision=-2; // Region/state
					elseif ($quality> 8)	$this->precision=-1;	// Country
					else	{				$this->precision=0;	// Not an address
											$this->codeRetourY='ZERO_RESULTS';
					}
					if ($this->precision<=0)
						die('Code Retour Yahoo 200 mais addressMatchType inconnu :'.$this->body.EOL);
				} else {
					$quality=$_match[1];
					$this->codeRetourY='OK';
					if     ($quality>87)	$this->precision=9;	// POI or Coordinate
					elseif ($quality>82)	$this->precision=8; // Address
					elseif ($quality>75)	$this->precision=7;	// Intersection
					elseif ($quality>70)	$this->precision=6;	// Street match
					elseif ($quality>40)	$this->precision=5; // Zip Code
					elseif ($quality>30)	$this->precision=4; // Town
					elseif ($quality>20)	$this->precision=3; // Sous/R<>gion
					elseif ($quality>10)	$this->precision=2; // Region/state
					elseif ($quality> 8)	$this->precision=1;	// Country
					else	{				$this->precision=0;	// Not an address
											$this->codeRetourY='ZERO_RESULTS';
					}
				}
				$this->latitudeDec=$_match[2];
				$this->longitudeDec=$_match[3];
				$this->latitudeDeg=dec2dms($this->latitudeDec);
				$this->longitudeDeg=dec2dms($this->longitudeDec);
				if (preg_match('!<line1>(.*)</line1>(?:|.*)<line2>(.*)</line2>(?:|.*)<line4>(.*)</line4>!Ui', $this->body, $_match)) {
					$adrValide=utf8_decode($_match[1].', '.$_match[2].', '.$_match[3]);
					$this->adresseValidee=strtoupper(str_replace('&apos;', "'", htmlspecialchars_decode($adrValide,ENT_QUOTES)));
					$this->codeRetourY='OK';
				} elseif (preg_match('!<line2>(.*)</line2>(?:|.*)<line4>(.*)</line4>!Ui', $this->body, $_match)) {
					$adrValide=utf8_decode($_match[2].', '.$_match[3]);
					$this->adresseValidee=strtoupper(str_replace('&apos;', "'", htmlspecialchars_decode($adrValide,ENT_QUOTES)));
					$this->codeRetourY='OK';
				} else 
					$this->codeRetourY='INVALID_REQUEST';
			} else
				die('Code Retour Yahoo 200 mais :'.$this->body.EOL);
			
			if ($this->latitudeDec==0 && ($this->longitudeDec==0 || $this->precision==0))	die('R<>ponse Yahoo = '.$this->body.EOL);
			return true;//array("lat"=>$lat,"lng"=>$lng,"address"=>$address);
		} else
			$this->codeRetourY='INVALID_REQUEST';
		
		return false;
	}
	
	private function geoCodeAdrCpVilleG($adresse,$cp,$ville,$pays, $proxy='') {
		$this->geocodeur='Google';
		$this->precision=0;
		
		$fileName=REP_TEMP.'geoCodeur'.$this->geocodeur;
		$tmp=@explode('|',file_get_contents($fileName));
		$nbQuery=@$tmp[1]*1;

		//$url='http://maps.googleapis.com/maps/api/geocode/xml?address='.urlencode($adresse.', '.$cp.' '. $ville.', '.$pays).'&sensor=false';
		$url='https://maps.googleapis.com/maps/api/geocode/xml?address='.urlencode($adresse.', '.$cp.' '. $ville.', '.$pays).'&sensor=false&key='.$this->apiKeyG;
		
		$page=getUrl($url, '', '', $this->referer, false, 'maps.googleapis.com', $proxy,3);
		//print_r($page);
		//die($url);
		//getUrl($url, $strCookies='', $postData='', $referer='', $debug=false, $host='', $proxy='', $timeout=0) {
		$this->body=$page['body'];
		$this->codeRetour=$page['code'];
		$this->header=$page['header'];

		$nbQuery++;
		file_put_contents($fileName, $this->geocodeur.'|'.$nbQuery.'|'.$this->codeRetour."|$url|".date('Y-m-d H:i:s').'|'.$this->body.EOL);		
		
		if ($this->codeRetour==200) {
		// REVERSE GEOCODING
		// http://maps.google.com/maps/geo?oe=utf-8&ll=48.808955,2.34227&output=xml&key=ABQIAAAAuKBtUyFonYJBl1fqfc78tRQvADPcxwXf3Q2QIE-M32vuSkrxiBRLUHDB_YSLeTscTDeWRKM_wILaaw
		
/*		
0 	 Unknown location. (Since 2.59)
1 	Country level accuracy. (Since 2.59)
2 	Region (state, province, prefecture, etc.) level accuracy. (Since 2.59)
3 	Sub-region (county, municipality, etc.) level accuracy. (Since 2.59)
4 	Town (city, village) level accuracy. (Since 2.59)
5 	Post code (zip code) level accuracy. (Since 2.59)
6 	Street level accuracy. (Since 2.59)
7 	Intersection level accuracy. (Since 2.59)
8 	Address level accuracy. (Since 2.59)
*/
/*			<?xml version="1.0" encoding="UTF-8"?>
			<GeocodeResponse>
			 <status>OK</status>
			 <result>
			  <type>street_address</type>
			  <formatted_address>0 Impasse Saint-Germier, 31270 Frouzins, France</formatted_address>
			  <address_component>
			   <long_name>0</long_name>
			   <short_name>0</short_name>
			   <type>street_number</type>
			  </address_component>
			  <address_component>
			   <long_name>Impasse Saint-Germier</long_name>
			   <short_name>Imp. Saint-Germier</short_name>
			   <type>route</type>
			  <geometry>
			   <location>
			    <lat>43.5160071</lat>
			    <lng>1.3252185</lng>
			   </location>
			   <location_type>RANGE_INTERPOLATED</location_type>
			  </geometry>
			  <partial_match>true</partial_match>
			 </result>
			</GeocodeResponse>
 */					
		$this->codeRetourG=@getTextInHtml($this->body, '<GeocodeResponse>', '<status>', '</status>');
		$location_type=@getTextInHtml($this->body, '<geometry>', '<location_type>', '</location_type>');
		switch($location_type) {
			case 'ROOFTOP':				// Precise geocode for which we have location information accurate down to street address precision.
				$this->precision=8;
				break;
			case 'RANGE_INTERPOLATED':	// Approximation (usually on a road) interpolated between two precise points (such as intersections). Interpolated results are generally returned when rooftop geocodes are unavailable for a street address.
				$this->precision=7;
				break;
			case 'GEOMETRIC_CENTER':	// Geometric center of a result such as a polyline (for example, a street) or polygon (region).
				$this->precision=6;
				break;
    		case 'APPROXIMATE':			// indicates that the returned result is approximate.
    			$this->precision=4;
    			break;
    		default:
    			if ($this->codeRetourG=='ZERO_RESULTS')	// Indicates that the geocode was successful but returned no results.
    				$this->precision=0;					// This may occur if the geocode was passed a non-existent address or a latlng in a remote location.
    			break;
		}
	
		$this->adresseValidee=strtoupper(str_replace('&apos;', "'", htmlspecialchars_decode(utf8_decode(@getTextInHtml($this->body, '<result', '<formatted_address>', '</formatted_address>')),ENT_QUOTES)));
		if (preg_match('!<geometry>(?:.*)<location>(?:.*)<lat>(.*)</lat>(?:.*)<lng>(.*)</lng>(?:.*)</location>!Uis', $this->body, $matches2)) {
			$this->latitudeDec=@$matches2[1];
			$this->longitudeDec=@$matches2[2];
			$this->latitudeDeg=dec2dms($this->latitudeDec);
			$this->longitudeDeg=dec2dms($this->longitudeDec);
		}
/*
200 G_GEO_SUCCESS				No errors occurred; the address was successfully parsed and its geocode has been returneds->geocodeur='Google';

400 G_GEO_BAD_REQUEST			A directions request could not be successfully parsed.
500 G_GEO_SERVER_ERROR			A geocoding or directions request could not be successfully processed, yet the exact reason for the failure is not known.
601 G_GEO_MISSING_QUERY			The HTTP q parameter was either missing or had no value. For geocoding requests, this means that an empty address was
								specified as input. For directions requests, this means that no query was specified in the input.
602 G_GEO_UNKNOWN_ADDRESS		No corresponding geographic location could be found for the specified address. This may be due to the fact that the address
								is relatively new, or it may be incorrect.
603 G_GEO_UNAVAILABLE_ADDRESS	The geocode for the given address or the route for the given directions query cannot be returned due to legal or contractual
								reasons.
604 G_GEO_UNKNOWN_DIRECTIONS	The GDirections object could not compute directions between the points mentioned in the query. This is usually because
								there is no route available between the two points, or because we do not have data for routing in that region.
610 G_GEO_BAD_KEY				The given key is either invalid or does not match the domain for which it was given.
620 G_GEO_TOO_MANY_QUERIES		The given key has gone over the requests limit in the 24 hour period.

<GeocodeResponse><status>OK</status><result><type>street_address</type><formatted_address>2 Route du Cure, 92410 Ville-d'Avray, France</formatted_address><address_component><long_name>2</long_name><short_name>2</short_name><type>street_number</type></address_component><address_component><long_name>Route du Cure</long_name><short_name>Route du Cure</short_name><type>route</type></address_component><address_component><long_name>Ville-d'Avray</long_name><short_name>Ville-d'Avray</short_name><type>locality</type><type>political</type></address_component><address_component><long_name>Hauts-de-Seine</long_name><short_name>92</short_name><type>administrative_area_level_2</type><type>political</type></address_component><address_component><long_name><EFBFBD>le-de-France</long_name><short_name>IDF</short_name><type>administrative_area_level_1</type><type>political</type></address_component><address_component><long_name>France</long_name><short_name>FR</short_name><type>country</type><type>political</type></address_component><address_component><long_name>92410</long_name><short_name>92410</short_name><type>postal_code</type></address_component><geometry><location><lat>48.8252802</lat><lng>2.1821693</lng></location>

		/*
"OK" indicates that no errors occurred; the address was successfully parsed and at least one geocode was returned.
    
    "OVER_QUERY_LIMIT" indicates that you are over your quota.
    "REQUEST_DENIED" indicates that your request was denied, generally because of lack of a sensor parameter.
    "" generally indicates that the query (address or latlng) is missing.
    UNKNOWN_ERROR indicates that the request could not be processed due to a server error. The request may succeed if you try again. 
		 */
		
		/*
		<?xml version="1.0" encoding="UTF-8"?>
		<GeocodeResponse>
		 <status>OVER_QUERY_LIMIT</status>
		 <error_message>You have exceeded your daily request quota for this API.</error_message>
		</GeocodeResponse>
				*/
				
		if ($this->codeRetourG=='OK') {
			//echo  $this->body.EOL;
			//echo 'Adr='.$this->adresseValidee.EOL;
			//echo 'Lat='.$this->latitudeDec.EOL;
			//echo 'Lon='.$this->longitudeDec.EOL;
			//echo 'Pre='.$this->precision.EOL;
			//die();				
			return true;
		}
		elseif ($this->codeRetourG=='ZERO_RESULTS' ||				// http://maps.googleapis.com/maps/api/geocode/xml?address=0044+RUE+JACQUES+-+LOUIS+HENON%2C+69004+LYON+4EME%2C+France&sensor=false
				$this->codeRetourG=='INVALID_REQUEST' ||		// http://maps.googleapis.com/maps/api/geocode/xml?address=0028+B+RUE+MGR+DE+BEAUMONT%2C+97400+SAINT+DENIS%2C+La+R%E9union&sensor=false
				$this->codeRetourG=='OVER_QUERY_LIMIT') {
			return false;
		} else  {
			echo $url.' ('.$this->codeRetour.')'.EOL;
			die('Code Retour Google ='.$this->codeRetourG.EOL);
		}
	}
	
	return false;
}
	
/** G<EFBFBD>ocodage d'une adresse		
0	Unknown location. (Since 2.59)
1 	Country level accuracy. (Since 2.59)
2 	Region (state, province, prefecture, etc.) level accuracy. (Since 2.59)
3 	Sub-region (county, municipality, etc.) level accuracy. (Since 2.59)
4 	Town (city, village) level accuracy. (Since 2.59)
5 	Post code (zip code) level accuracy. (Since 2.59)
6 	Street level accuracy. (Since 2.59)
7 	Intersection level accuracy. (Since 2.59)
8 	Address level accuracy. (Since 2.59)
9   Address level accuracy +++
 **/	
	function geoCodeAdresse($adrNum, $adrIndRep, $adrTypeVoieCourt,  $adrTypeVoieLong, $adrLibVoie, $cp, $ville='', $pays='France', $codeRivoli='') {
		$iDb=new WDB();
		$tabRep=array();
		$adresse=addslashes(trim(preg_replace('/ +/',' ', "$adrNum $adrIndRep $adrTypeVoieLong $adrLibVoie")));
		$ville=addslashes($ville);
		$ret=$iDb->select(	'zonageXY', 'lat, lon, 	l93_x, l93_y, alt, precis, adresseValidee, dateInsert', "address='$adresse' AND adr_cp='$cp' AND adr_ville='$ville'",false, MYSQL_ASSOC);
		if (count($ret)>0) {
			$zonage=$ret[0];
			$this->precision=$zonage['precis']*1;
			$this->adresseValidee=strtoupper($zonage['adresseValidee']);
			$this->latitudeDec=$zonage['lat']*1;
			$this->longitudeDec=$zonage['lon']*1;
			$this->l93x=$zonage['l93_x'];
			$this->l93y=$zonage['l93_y'];
			$this->altitude=$zonage['alt'];
			if ($this->altitude==NULL && $this->precision>5) {
				if ($this->accesDist) {
					$alt=$this->getAltitude($this->latitudeDec,$this->longitudeDec);
					$this->altitude=$alt['alt'];
				}
			}
			if ($this->l93x==NULL && $this->precision>5) {
				$tmp=geos2lambert93($this->latitudeDec,$this->longitudeDec);
				$this->l93x=$tmp['x_93'];
				$this->l93y=$tmp['y_93'];
			}
			if ($this->altitude<>NULL || $this->l93x<>NULL) {
					$tabUpdate=array(	'alt'=>$this->altitude,
										'l93_x'=>$this->l93x,
										'l93_y'=>$this->l93y);
										//'dateUpdate'=>DATETIME,
					$iDb->update('zonageXY', $tabUpdate, "address='$adresse' AND adr_cp='$cp' AND adr_ville='$ville'");					
			}
				
			$this->latitudeDeg=dec2dms($this->latitudeDec);
			$this->longitudeDeg=dec2dms($this->longitudeDec);
			$this->enCache=true;
			$ligne='En base : '.print_r($zonage, true);	
		} elseif ($this->accesDist) {
			$this->enCache=false;
			$retG=$retY=$retC=$retM=false;
			if (substr($codeRivoli,0,3)*1==971)		$pays='Guadeloupe';
			elseif (substr($codeRivoli,0,3)*1==972)	$pays='Martinique';
			elseif (substr($codeRivoli,0,3)*1==973)	$pays='Guyane';
			elseif (substr($codeRivoli,0,3)*1==974)	$pays='La R<>union';
			elseif (substr($codeRivoli,0,3)*1==975)	$pays='Saint-Pierre-et-Miquelon';
			elseif (substr($codeRivoli,0,3)*1==976)	$pays='Mayotte';
			elseif (substr($codeRivoli,0,3)*1==977)	$pays='Saint-Martin';
			elseif (substr($codeRivoli,0,3)*1==978)	$pays='Saint-Barth<74>l<EFBFBD>my';
			elseif (substr($codeRivoli,0,3)*1==985)	$pays='Saint-Pierre-et-Miquelon';//98501 <20> 98502		St Pierre et Miquelon (anciens codes non compatibles 97501 <20> 97502)
			elseif (substr($codeRivoli,0,3)*1==986)	$pays='Wallis-et-Futuna';		// 98601 <20> 98617 		Mayotte (anciens codes non compatibles : 98501 <20> 98517)
			elseif (substr($codeRivoli,0,3)*1==987)	$pays='Polyn<79>sie fran<61>aise';	// 98711 <20> 98758,98799 	Polyn<79>sie fran<61>aise (sans changement)
			elseif (substr($codeRivoli,0,3)*1==988)	$pays='Nouvelle Cal<61>donie';		// 98801 <20> 98832 		Nouvelle Cal<61>donie (sans changement)
			elseif (substr($codeRivoli,0,3)*1==989)	$pays='Wallis-et-Futuna';		// 98911 <20> 98913 		Wallis et Futuna (anciens codes non compatibles : 98611 <20> 98613)
			elseif (substr($codeRivoli,0,2)*1>96 && substr($codeRivoli,0,2)*1<99)
				die("Code Rivoli non g<>r<EFBFBD> : '$codeRivoli'".EOL);
				
			$retG=$this->geoCodeAdrCpVilleG($adresse, $cp, $ville, $pays);
			if ($this->codeRetourG=='OVER_QUERY_LIMIT' || $this->precision==0) {
				$retY=$this->geoCodeAdrCpVilleY("$adresse, $cp $ville, $pays");
				if ($this->codeRetourY=='OVER_QUERY_LIMIT' || $this->codeRetourY<>'OK' || $this->precision==0)
					$retC=$this->geoCodeAdrCpVilleC($adresse, $cp, $ville, $pays);
				//	$retM=$this->geoCodeAdrCpVilleM($adresse, $cp, $ville, $pays);
			}
			
			// On r<>cup<75>re l'altitude
			//echo ('Pr<50>cision : '. $this->precision.EOL);
			if ($this->precision>5 && $this->latitudeDec<>0 && $this->longitudeDec<>0) {	
				$alt=$this->getAltitude($this->latitudeDec,$this->longitudeDec);
				$this->altitude=$alt['alt'];
			}
			$tmp=geos2lambert93($this->latitudeDec,$this->longitudeDec);
			$this->l93x=$tmp['x_93'];
			$this->l93y=$tmp['y_93'];
				
			//echo ('$retG='.$retG.EOL.'$retY='.$retY.EOL.'$retC='.$retC.EOL.'$retM='.$retM.EOL);
			
			if ($retG || $retY || $retM || $retC) {
				$adresse=stripslashes(strtoupper(trim(preg_replace('/ +/',' ', "$adrNum $adrIndRep $adrTypeVoieLong $adrLibVoie"))));
				$ville=stripslashes($ville);
				$tabInsert=array(	'address'=>$adresse,
									'adr_cp'=>$cp,
									'adr_ville'=>$ville,
									'adrNum'=>$adrNum,
									'adrIndRep'=>$adrIndRep,
									'adrTypeVoie'=>$adrTypeVoieCourt,
									'adrLibVoie'=>$adrLibVoie,
									'rivoli'=>$codeRivoli,
									'adresseValidee'=>$this->adresseValidee,
									'lat'=>$this->latitudeDec,
									'lon'=>$this->longitudeDec,
									'l93_x'=>$this->l93x,
									'l93_y'=>$this->l93y,
									'alt'=>$this->altitude,
									'precis'=>$this->precision,
									'source'=>$this->geocodeur,
									'dateInsert'=>DATETIME,
								);
				/** Insertion de l'adresse que si elle est possible
				 **/
				//echo ('Avant Insert '.$this->latitudeDec.EOL.$this->longitudeDec.EOL);
				if ($this->latitudeDec<>0 && $this->longitudeDec<>0 && $this->precision>0) {
				//	print_r($tabInsert);
					//echo ($this->latitudeDec.EOL.$this->longitudeDec.EOL);
					$iDb->insert('zonageXY', $tabInsert);
					if (mysql_errno()>0)	echo mysql_error().die(EOL);
				} else {
					$this->precision=$this->latitudeDec=$this->longitudeDec=0;
					$this->adresseValidee='';
					$this->latitudeDeg=$this->longitudeDeg='';
					$ligne="Erreur Geocodage : $adrNum, $adrIndRep, $adrTypeVoieCourt, $adrTypeVoieLong, $adrLibVoie, $cp, $ville, $pays, $codeRivoli, ".$this->body;
					//die('Latitude et longitude <20> zero :'.$this->latitudeDec.EOL.$this->longitudeDec.EOL);
				}
				//echo ('Apr<70>s Insert : '.$this->latitudeDec.EOL.$this->longitudeDec.EOL);
				$ligne='Acc<63>s Google : '.print_r($tabInsert, true);	
			}
			else {
				$this->precision=$this->latitudeDec=$this->longitudeDec=0;
				$this->adresseValidee='';
				$this->latitudeDeg=$this->longitudeDeg='';
				$ligne="Erreur : $adrNum, $adrIndRep, $adrTypeVoieCourt, $adrTypeVoieLong, $adrLibVoie, $cp, $ville, $pays, $codeRivoli, ".$this->body;
			}
		}
	}
/*      if (latnorth == 'S') latdir='-'
      if (lonwest == 'W') longdir='-'
      latdec=Math.round(lat1)+lat2/60
      londec=Math.round(long1)+long2/60
      gmdatalat=latdir+latdec.toFixed(6)
      gmdatalon=longdir+londec.toFixed(6)
      //GM_log( 'Latdec:'+gmdatalat+' LongDec:'+gmdatalon)
      gmx=gmdatalat*Math.pow(10,6)
      gmy=gmdatalon*Math.pow(10,6)
      if (gmx < 0) gmx=gmx+Math.pow(2,32)*/


/** Retourne la distance en kilom<EFBFBD>tres entre 2 coordonn<EFBFBD>es GPS <EFBFBD> la surface de la terre
 **
 ** @param double $latA		Latitude  du point A en d<EFBFBD>cimal
 ** @param double $lonA		Longitude du point A en d<EFBFBD>cimal
 ** @param double $latB		Latitude  du point B en d<EFBFBD>cimal
 ** @param double $lonB		Longitude du point B en d<EFBFBD>cimal
 ** @param string $sphere	Type de sph<EFBFBD>re
 ** @return unknown
 **/
function distance($latA=0, $lonA=0, $latB=0, $lonB=0, $sphere='GRS80') {
	switch ($sphere) {
		case 'HAYFORD':	$R=6378.388;	// Demi grand axe ou Rayon de la sph<70>re International Hayford 1909 en Kms
						$f=1/297;		// Aplatissement
						break;
		case 'PICARD':	$R=6371.598;	// Demi grand axe ou Rayon de la sph<70>re Picard en Kms
						break;
		case 'GRS80':
		default:		$R=6378.137;	// Demi grand axe ou Rayon de la sph<70>re GRS80 en Kms
						$f=1/298.257222101;
						break;	
	}
	$a=pi()/180;
	$e=$latA*$a;
	$f=$lonA*$a; 
	$g=$latB*$a;
	$h=$lonB*$a;

	$j=acos(cos($e)*cos($g)*cos($f)*cos($h) + cos($e)*sin($f)*cos($g)*sin($h) + sin($e)*sin($g));

	return round($R*$j,3); // div par 1.852 ==> r<>sultat en miles nautiques
}

/** Retourne la distance en m<EFBFBD>tres entre 2 points sur une surface plane
 ** @param	double $x1	Coordonn<EFBFBD>e X du point 1 en d<EFBFBD>cimal
 ** @param	double $y1	Coordonn<EFBFBD>e Y du point 1 en d<EFBFBD>cimal
 ** @param	double $x2	Coordonn<EFBFBD>e X du point 2 en d<EFBFBD>cimal
 ** @param	double $y2	Coordonn<EFBFBD>e Y du point 2 en d<EFBFBD>cimal
 ** @return double
 **/
function distanceLambert($x1,$y1,$x2,$y2) {
	return sqrt ( pow($x1-$x2,2) + (pow($y1-$y2,2) ) );
}

/** Retourne l'altitude en m<>tres d'un point GPS. Shuttle Radar Topography Mission (SRTM) elevation data (resolution 90mx90m)
 ** @url	http://ws.geonames.org/export/web-services.html 
 ** @param	double $tabLatLon	Latitude du point ou Tableau de Points avec latitudes longitudes array('lat'=>x,'lon'=>y)
 ** @param	double $lon			Longitude ou false si tableau
 ** @return double
 **/
function getAltitude($tabLatLon,$lon=false) {
	if ($lon)	{
		/*Elevation - SRTM3
Shuttle Radar Topography Mission (SRTM) elevation data. SRTM consisted of a specially modified radar system that flew onboard the Space Shuttle Endeavour during an 11-day mission in February of 2000. The dataset covers land areas between 60 degrees north and 56 degrees south.
This web service is using SRTM3 data with data points located every 3-arc-second (approximately 90 meters) on a latitude/longitude grid. Documentation : Nasa

Webservice Type : REST
Url : api.geonames.org/srtm3?
Parameters : lat,lng;
sample area: ca 90m x 90m Result : a single number giving the elevation in meters according to srtm3, ocean areas have been masked as "no data" and have been assigned a value of -32768
Example http://api.geonames.org/srtm3?lat=50.01&lng=10.2&username=demo

This service is also available in XML and JSON format :api.geonames.org/srtm3XML?lat=50.01&lng=10.2&username=demo api.geonames.org/srtm3JSON?lat=50.01&lng=10.2&username=demo

The text version of the service also accepts a list of lat/lng for the parameters 'lats' and 'lngs'. On the free server the number of points per call is limited to 20, for the premium service the limit is 2000:

    [0] => the daily limit of 30000 credits 78.31.45.206 has been exceeded. Please throttle your requests or use the commercial service
    [1] => daily
    [2] => 30000
    [3] =>  78.31.45.206

*/
		$url="http://ws.geonames.org/srtm3?lat=$tabLatLon&lng=$lon&username=scores_et_decisions&password=bj10sx";
		$page=getUrl($url, '', '', '', false, '', '', 3);
		$alt=trim($page['body']);
		if (preg_match('/the (.*) limit of (.*) credits(.*) has been exceeded\. Please throttle your requests or use the commercial service/Uis', $alt, $matches)) {
			print_r($matches);
			$alt=NULL;
		} elseif ($page['code']<>200 || $alt==-32768)	{
			if ($page['code']<>200) echo '==========================>'. $page['code']. '<=================================='.EOL;
			$alt=NULL;
		}
		return array('lat'=>$tabLatLon,'lon'=>$lon,'alt'=>$alt);
	} else {
		$LatLon='';
		$Lats='';
		$Lons='';
		foreach ($tabLatLon as $i=>$latlon) {
			$lat=$latlon['lat'];
			$lon2=$latlon['lon'];
			if ($i>0) {
				$LatLon.="%0D%0A";
				$Lats.=',';
				$Lons.=',';
			}
			$LatLon.="$lat%2C$lon2";
			$Lats.=$lat;
			$Lons.=$lon;	
		}
		//	http://api.geonames.org/srtm3?lats=50.01,51.01&lngs=10.2,11.2&username=demo
	/*	$url="http://api.geonames.org/srtm3?lats=$Lats&lng=$Lons&username=demo";
		$page=getUrl($url, '', '', '', false, '', '', 3);
		$alt=trim($page['body']);
		die($alt);
		if ($page['code']<>200 || $alt==-32768)	{
			if ($page['code']<>200) echo '==========================>'. $page['code']. '<=================================='.EOL;
			$alt=NULL;
		}
		return array('lat'=>$tabLatLon,'lon'=>$lon,'alt'=>$alt);
		*/
	}
	
	$url='http://www.fredorando.fr/Pages/altitude.php';
	$post=array('LatLon'=>$LatLon,'qui'=>'');
	$page=getUrl($url, '', $post, '', false, '', '', 3);
	$body=$page['body'];
	$tabRet=false;
	if (preg_match('/<BODY>\s+IP\:(?:.*)<br\/>\s+url\:(?:.*)<br\/>(?:.*)<br\/>(?:.*)<br\/>(.*)<br\/>(?:.*)<a href="http/Uis', $body, $matches)) {
		$tabRet=array();
		$tabTmp=explode('<br />', $matches[1]);
		foreach($tabTmp as $i=>$ret) {
			$ret2=explode(',',trim($ret));
			$alt=$ret2[2];
			if ($alt==-32768)	$alt=NULL;
			$tabRet[$i]=@array('lat'=>$ret2[0],'lon'=>$ret2[1],'alt'=>$alt);
		}
		unset($tabRet[$i]);
	}
	if ($lon)	return $tabRet[0];
	return $tabRet;
}

}
if ( !function_exists('json_decode') ){
    function json_decode($content, $assoc=false){
                require_once 'Services/JSON.php';
                if ( $assoc ){
                    $json = new Services_JSON(SERVICES_JSON_LOOSE_TYPE);
        } else {
                    $json = new Services_JSON;
                }
        return $json->decode($content);
    }
}

if ( !function_exists('json_encode') ){
    function json_encode($content){
                require_once 'Services/JSON.php';
                $json = new Services_JSON;
               
        return $json->encode($content);
    }
}
?>