webservice/library/Metier/partenaires/classMMap.php

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

require_once 'framework/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.'°'.$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;
}


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];
}

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;
}

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;
}

function ALG0013($Tx,$Ty,$Tz,$D,$Rx,$Ry,$Rz,$U) {
        $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;
}

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;

}

function ALG0021($phi,$a,$e) {
        $N = $a/sqrt( 1 - $e * $e * sin($phi) * sin($phi) );
        return $N;

}

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
        {
                        echo ('coucou');
                $Xs = $X0;
                $Ys = $Y0 + $C * exp(-1 * $n * ALG0001($phi0,$e));
        }

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

        return $tab;

}

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

        switch ($orig)
        {
                case 'LII' :
                        $n = 0.7289686274;
                        $c = 11745793.39;
                        $Xs = 600000;
                        $Ys = 6199695.768;
                        $lambdac = 0.04079234433; // pour greenwich

                        $e = 0.08248325676; //(première excentricité de l ellipsoïde Clarke 1880 français)

                        $he = 100;
                        $a = 6378249.2;

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

                case 'L93' :
                        $n = 0.7256077650;
                        $c = 11745255.426;
                        $Xs = 700000;
                        $Ys = 12655612.050;
                        $lambdac = 0.04079234433; // pour greenwich

                        $e = 0.08248325676; //(première excentricité de l ellipsoïde Clarke 1880 français)

                        $he = 100;
                        $a = 6378249.2;

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

                case 'LIIe' :
                default:
                        $n = 0.7289686274;
                        $c = 11745793.39;
                        $Xs = 600000;
                        $Ys = 8199695.768;
                        $lambdac = 0.04079234433; // pour greenwich

                        $e = 0.08248325676; //(première excentricité de l ellipsoïde Clarke 1880 français)

                        $he = 100;
                        $a = 6378249.2;

                        $Tx = -168;
                        $Ty = -60;
                        $Tz = +320;
                        $D  = 0;
                        $Rx = $Ry = $Rz = 0;
                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ï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ées géographiques WGS84 (phi_w,lambda_w) en coordonnées carté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ées cartésiennes WGS84 (X_w,Y_w,Z_w) en coordonnées carté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ées cartésiennes NTF (X_n,Y_n,Z_n) en coordonné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é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ées géographiques NTF (phi_n,lambda_n) en coordonnées projetées en Lambert II é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 à la longitude en radian de Paris (2°20'14.025" E) par rapport à Greenwich
    //puis on calcule la latitude isomé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ème WGS84
	$a=6378137; 		// demi grand axe de l'ellipsoide (m)
	$e=0.08181919106;	// première excentricitè de l'ellipsoide

	// Paramètres de projections
    $l0=$lc=deg2rad(3);		// longitude de réfé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ées à l'origine
    $y0=6600000;	//coordonnées à l'origine

	// coordonnées du point à 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é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ées
    $ys=$y0+$c*exp(-1*$n*$gl0);

	//calcul des coordonné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));
**/


class MMap {

	private $accesDist = true;
	private $referer ='';
	private $body = '';
	private $header = '';
	private $cookie = '';
	private $codeRetour = 0;
	private $codeRetourG = 0;	// Code retour de l'api Google
	private $codeRetourY = 0;	// Code retour de l'api Yahoo
	private $codeRetourM = 0;	// Code retour de l'api Mappy

	public $geocodeur='';

	public $latitudeDec= 0;		// Latitude en Décimal
	public $latitudeDeg= 0;		// Latitude en Dégrés
	public $longitudeDec= 0;	// Longitude en Décimal
	public $longitudeDeg= 0;	// Longitude en Dégrés

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

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

	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='') {

		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;
		//print_r($result);die();
		//$this->codeRetourM=

		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);

		$this->geocodeur='Mappy';

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

		return false;
	}


	private function geoCodeAdrCpVilleY($address) {
		$_url = 'http://api.local.yahoo.com/MapsService/V1/geocode';
		$_url .= sprintf('?appid=%s&location=%s',$this->apiKeyY,rawurlencode($address));
		$_result = false;
		$page=getUrl($_url, '', '', '', false, '', '', 3);
		$_result=$page['body'];
//		die("URL='$_url'".EOL.print_r($page));
		$try=1;
		while ($page['code']!=200 && $page['code']!=403) {
			$page=getUrl($_url, '', '', '', false, '', '', 3);
			$_result=$page['body'];
			echo 'Yahoo '.$page['code']." : $address... Attente...".EOL;
			$try++;
			if ($try>=4)	break;
			sleep(1);
		}
		if ($page['code']==400 && preg_match('/<Message>unable to parse location<\/Message>/', $_result)) {
			$this->latitudeDec=$this->longitudeDec=0;
            $this->latitudeDeg=$this->longitudeDeg=0;
			$this->geocodeur='Yahoo';
			$this->precision=0;
			return false;
		} elseif ($page['code']==403 && preg_match('/<Message>limit exceeded<\/Message>/', $_result)) {
			$this->latitudeDec=$this->longitudeDec=0;
            $this->latitudeDeg=$this->longitudeDeg=0;
			$this->geocodeur='Yahoo';
			$this->precision=0;
			$this->codeRetourY=620;	// On adopte les codes Google
		} elseif ($page['code']==408) { // Connexion impossible au site du partenaire
			$this->latitudeDec=$this->longitudeDec=0;
            $this->latitudeDeg=$this->longitudeDeg=0;
			$this->geocodeur='Yahoo';
			$this->precision=0;
			$this->codeRetourY=500;	// On adopte les codes Google
	    } elseif ($page['code']==200) {
			preg_match('!<Latitude>(.*)</Latitude><Longitude>(.*)</Longitude>!U', $_result, $_match);
			$lng = $_match[2];
			$lat = $_match[1];
			preg_match('/warning="(.*)"/Ui', $_result, $_match);
/*
1       Country level accuracy. (Since 2.59)
2       Region (state, province, prefecture, etc.) level accuracy. (Since 2.59)
7       Intersection level accuracy. (Since 2.59)
*/
			$this->precision=0;	// Unknown location
			switch(@$_match[1]) {
				case "The exact location could not be found. Here is a nearby county.":				$this->precision=3;     break;  // Sub-region
				case "The city could not be found. Here is the state/province.":					$this->precision=3;     break;  // Sub-region
				case "The exact location could not be found. Here is a nearby neighorhood.":		$this->precision=4;     break;  // Town level
				case "The exact location could not be found. Here is the center of the ZIP code.":
				case "The street could not be found. Here is the center of the city.":			$this->precision=5;     break;  // Post code (zip code)
				case "The street number could not be found. Here is a nearby location.":
				case "The Street name might have been changed":						$this->precision=6;	break;	// Street level accuracy
				case '':										$this->precision=8;     break;  // Address level
				default:
					if (preg_match('/The exact location could not be found, here is the closest match/', @$_match[1]))
						$this->precision=6;
					else
						die('Precisions Inconnue : "'.$_match[1].'"'.EOL.$_result);
					break;
			}

			preg_match('!<Address>(.*)</Address>(?:|.*)<City>(.*)</City>(?:|.*)<State>(.*)</State>!Ui', $_result, $_match);
                        $adrValide=utf8_decode($_match[1].', '.$_match[2].', '.$_match[3]);
			$this->adresseValidee=strtoupper(str_replace('&apos;', "'", htmlspecialchars_decode($adrValide,ENT_QUOTES)));
                	$this->latitudeDec=$lat;
                	$this->longitudeDec=$lng;
               	 	$this->latitudeDeg=dec2dms($this->latitudeDec);
                	$this->longitudeDeg=dec2dms($this->longitudeDec);
			$this->geocodeur='Yahoo';
			if ($lat==0 && $lng==0 || $this->precision==0)	die('Réponse Yahoo = '.$_result.EOL);
			return true;//array("lat"=>$lat,"lng"=>$lng,"address"=>$address);
		}
		else {die("URL='$_url'".EOL.print_r($page));
			return false;
		}
	}

	private function geoCodeAdrCpVilleG($adresse,$cp,$ville,$pays, $proxy='') {

		$url='http://maps.google.com/maps/geo?q='.urlencode($adresse.', '.$cp.' '. $ville.', '.$pays).'&output=xml&key='.$this->apiKeyG;

		$page=getUrl($url, '', '', $this->referer, false, 'maps.google.com', $proxy,3);
		//getUrl($url, $strCookies='', $postData='', $referer='', $debug=false, $host='', $proxy='', $timeout=0) {
		$this->body=$page['body'];
		$this->codeRetour=$page['code'];
		$this->header=$page['header'];

		// REVERSE GEOCODING
		// http://maps.google.com/maps/geo?oe=utf-8&ll=48.808955,2.34227&output=xml&key=ABQIAAAAuKBtUyFonYJBl1fqfc78tRQvADPcxwXf3Q2QIE-M32vuSkrxiBRLUHDB_YSLeTscTDeWRKM_wILaaw

		//die($this->body=$page['body']);
		//{"name":"3 rue viète, 75017 paris, france","Status":{"code":200,"request":"geocode"},"Placemark":[{"id":"p1","address":"3, Rue Viète, 75017 17ème Arrondissement, Paris, France","AddressDetails":{"Country":{"CountryNameCode":"FR","AdministrativeArea":{"AdministrativeAreaName":"Ile-de-France","SubAdministrativeArea":{"SubAdministrativeAreaName":"Paris","Locality":{"LocalityName":"Paris","DependentLocality":{"DependentLocalityName":"17ème Arrondissement","Thoroughfare":{"ThoroughfareName":"3, Rue Viète"},"PostalCode":{"PostalCodeNumber":"75017"}}}}}},"Accuracy": 8},"Point":{"coordinates":[2.306174,48.883705,0]}}]}
		/*<?xml version="1.0" encoding="UTF-8"?><kml xmlns="http://earth.google.com/kml/2.0"><Response><name>3 rue viète, 75017 paris, france</name><Status><code>200</code><request>geocode</request></Status><Placemark id="p1"><address>3, Rue Viète, 75017 17ème Arrondissement, Paris, France</address><AddressDetails Accuracy="8" xmlns="urn:oasis:names:tc:ciq:xsdschema:xAL:2.0"><Country><CountryNameCode>FR</CountryNameCode><AdministrativeArea><AdministrativeAreaName>Ile-de-France</AdministrativeAreaName><SubAdministrativeArea><SubAdministrativeAreaName>Paris</SubAdministrativeAreaName><Locality><LocalityName>Paris</LocalityName><DependentLocality><DependentLocalityName>17ème Arrondissement</DependentLocalityName><Thoroughfare><ThoroughfareName>3, Rue Viète</ThoroughfareName></Thoroughfare><PostalCode><PostalCodeNumber>75017</PostalCodeNumber></PostalCode></DependentLocality></Locality></SubAdministrativeArea></AdministrativeArea></Country></AddressDetails><Point><coordinates>2.306174,48.883705,0</coordinates></Point></Placemark></Response></kml>*/

/*
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->codeRetourG=@getTextInHtml($this->body, '<Status>', '<code>', '/code>');
		$this->precision=@getTextInHtml($this->body, '<AddressDetails Accuracy', '="', '" ');

		$this->adresseValidee=strtoupper(str_replace('&apos;', "'", htmlspecialchars_decode(@getTextInHtml($this->body, '<Placemark', '<address>', '</address>'),ENT_QUOTES)));
		$strTmp=@getTextInHtml($this->body, '<Point><coordinates>', '<coordinates>', '</coordinates>');
		$tabTmp=explode(',', $strTmp);
		$this->latitudeDec=@$tabTmp[1];
		$this->longitudeDec=$tabTmp[0];
		$this->latitudeDeg=dec2dms($this->latitudeDec);
		$this->longitudeDeg=dec2dms($this->longitudeDec);

		$this->geocodeur='Google';

/*
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.
*/
		if ($this->codeRetourG==200)
			return true;

		return false;
}


	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, 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->latitudeDeg=dec2dms($this->latitudeDec);
			$this->longitudeDeg=dec2dms($this->longitudeDec);
			$this->enCache=true;
			$ligne='En base : '.print_r($zonage, true);
		} else {
			$this->enCache=false;
			$retM=$retY=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élémy';
			elseif (substr($codeRivoli,0,3)*1==986)	$pays='Wallis-et-Futuna';
			elseif (substr($codeRivoli,0,3)*1==987)	$pays='Polynésie française';
			elseif (substr($codeRivoli,0,3)*1==988)	$pays='Nouvelle Calédonie';
			elseif (substr($codeRivoli,0,2)*1>96 && substr($codeRivoli,0,2)*1<99)
				die("Code Rivoli non géré : '$codeRivoli'".EOL);

			$retG=$this->geoCodeAdrCpVilleG($adresse, $cp, $ville, $pays);
			if ($this->codeRetourG==620 || $this->precision==0)
				if ($this->codeRetourG==620)
					$retY=$this->geoCodeAdrCpVilleY("$adresse, $cp $ville");

				if ($this->codeRetourY==620 || $this->precision==0)
					$retM=$this->geoCodeAdrCpVilleM($adresse, $cp, $ville, $pays);

			if ($retG || $retY || $retM) {
				$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,
									'precis'=>$this->precision,
									'source'=>$this->geocodeur,
									'dateInsert'=>DATETIME,
								);
				/** Insertion de l'adresse que si elle est possible
				 **/
				if ($this->latitudeDec<>0 && $this->longitudeDec<>0)
					$iDb->insert('zonageXY', $tabInsert);

				$ligne='Accè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;
			}
		}
		/*$fp=fopen('/var/www/log/accesGMap.log', 'a');
		fwrite($fp,$ligne.EOL);
		fclose($fp);*/
	}
/*      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ètres entre 2 points à la surface de la terre
 ** Calcul effectué avec la sphère « GRS80 » de rayon R = 6378,187 km
 ** Autre sphère possible : 	  « Picard » de rayon R = 6371,598 km
 **
 ** @param double $latA Latitude  du point A en décimal
 ** @param double $lonA Longitude du point A en décimal
 ** @param double $latB Latitude  du point B en décimal
 ** @param double $lonB Longitude du point B en décimal
 ** @return unknown
 **/
function distance($latA=0, $lonA=0, $latB=0, $lonB=0) {
//s(AB) = arc cos (sinjA sinjB + cos jA cosjB cosdl)
//	avec dl = lB - lA
	$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(6378.187*$j,3); // div par 1.852 ==> résultat en miles nautiques
}
/*function distance(lat_a, lon_a, lat_b, lon_b)
{
  a = Math.PI / 180;
  lat1 = lat_a * a;
  lat2 = lat_b * a;
  lon1 = lon_a * a;
  lon2 = lon_b * a;

  t1 = Math.sin(lat1) * Math.sin(lat2);
  t2 = Math.cos(lat1) * Math.cos(lat2);
  t3 = Math.cos(lon1 - lon2);
  t4 = t2 * t3;
  t5 = t1 + t4;
  rad_dist = Math.atan(-t5/Math.sqrt(-t5 * t5 +1)) + 2 * Math.atan(1);

  return (rad_dist * 3437.74677 * 1.1508) * 1.6093470878864446;
}
*/

}
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);
    }
}
?>