<?php 


function teilweg($db,$p1,$p2,$po,$stueckel){
//
// Bestimmt den naechstgelegenen Weg zu P1/P2
// Teilt diesen Weg an p1,p2 auf, sofern das nicht sowieso schon der Anfangs- oder Endpunkt ist
// gibt den Punkt zurueck, an dem die Teilung stattfand, also entwerder den schon
// vorhandenen source/target oder den neu entstandenen.
// Die Nummer des evtl neu entstehenden Knotens wird in po uebergeben. In der Datenbank stehen
// naemlich noch nicht alle Wege, manche Knoten werden erst spaeter erstellt. Deshalb muss
// irgendwo uebergeordnet ueber die po's buchgefuehrt werden.
//
// Ist kein Weg in 10cm Umkreis, wird nichts gemacht und der Rueckgabewert is der uebergebene po
// Wurde nicht geteilt, weil da eh schon ein Start/Endpunkt war, wird dieser zurueckgegeben
// Wurde nicht geteilt, weils nicht noetig ist und stueckel das nicht verlangt hat, wird auch -1 
// zurueckgegeben
//
 if($stueckel){$ret=$po;}else{$ret=-1;}
 $query="select gid,source,target,st_distance(the_geom,st_transform(st_geometryfromtext('POINT($p1 $p2)',900913),4326)) as dist,st_astext(the_geom) as weg from routing_ways where the_geom && st_transform(expand(st_geometryfromtext('POINT($p1 $p2)',900913),5),4326) order by dist limit 1;";
// echo "Teile \n";
 $result=pg_query($db,$query);
 $r=pg_fetch_array($result);
 $dist=100;
 if($r){
  $gid=$r['gid'];
  $dist=$r['dist'];
 }
// echo "Teile liefert $dist \n";
 if($dist<0.000001){
  $query="select ST_Line_Locate_Point((select the_geom from routing_ways where gid=$gid),st_transform(st_geometryfromtext('POINT($p1 $p2)',900913),4326)) as ll;";
  $result=pg_query($db,$query);
  $e=pg_fetch_array($result);
  $ll=$startll=$e["ll"];
//  echo "Teile weg $gid an $p1 $p2 bei $ll (stck=$stueckel) und vergebe $po \n";
  if($ll<=0.005)$ret=$r['source'];
  if($ll>=0.995)$ret=$r['target'];
  if(($ll>0.005)&&($ll<0.995)&&($stueckel==0)){$ret=-1;}
  if(($ll>0.005)&&($ll<0.995)&&($stueckel==1)){
//   echo "-- Teile Weg $gid bei $ll am Punkt $p1 $p2\n";
   pg_query($db, "BEGIN WORK");  
   $query="select max(gid) as m from routing_ways;";
   $result=pg_query($db,$query);
   $e=pg_fetch_array($result); 
   $gi=$e['m']+1;
   $query="select gid,class_id,osm_id,source,target,waysource,the_geom,length,costfactor,hdiff from routing_ways where gid=$gid;";
   $result=pg_query($db,$query);
   $e=pg_fetch_array($result);
   $gid=$e['gid'];
   $q="insert into routing_ways (gid,class_id,osm_id,source,target,waysource,the_geom,length,hdiff,costfactor)  values ($gi,".$e['class_id'].",".$e['osm_id'].",".$e['source'].",".$e['target'].",'".$e['waysource']."',(select the_geom from routing_ways where gid=$gid),".$e['length'].",".$e['hdiff'].",".$e['costfactor'].");";
   $r=pg_query($db,$q);
   $q="update routing_ways set target=$po where gid=$gid;";
   $r=pg_query($db,$q);   
   $q="update routing_ways set source=$po where gid=$gi;";
   $r=pg_query($db,$q);   
   $q="update routing_ways set hdiff=hdiff*$ll  where gid=$gid;";
   $r=pg_query($db,$q);   
   $q="update routing_ways set hdiff=hdiff*(1-$ll)  where gid=$gi;"; 
   $r=pg_query($db,$q);   
   $q="update routing_ways set the_geom=st_line_substring(the_geom,$ll,1) where gid=$gi;";
   $r=pg_query($db,$q);   
   $q="update routing_ways set the_geom=st_line_substring(the_geom,0,$ll) where gid=$gid;";
   $r=pg_query($db,$q);
   $q="update routing_ways set length=st_length(st_transform(the_geom,900913))*0.66  where gid=$gid;";
   $r=pg_query($db,$q);
   $q="update routing_ways set length=st_length(st_transform(the_geom,900913))*0.66  where gid=$gi;";
   $r=pg_query($db,$q);
   $q="update routing_ways set x1=st_x(st_startpoint(the_geom))  where gid=$gi;";
   $r=pg_query($db,$q);   
   $q="update routing_ways set y1=st_y(st_startpoint(the_geom))  where gid=$gi;";       
   $r=pg_query($db,$q);   
   $q="update routing_ways set x2=st_x(st_endpoint(the_geom))  where gid=$gi;";       
   $r=pg_query($db,$q);   
   $q="update routing_ways set y2=st_y(st_endpoint(the_geom))  where gid=$gi;";         
   $r=pg_query($db,$q);   
   $q="update routing_ways set x1=st_x(st_startpoint(the_geom))  where gid=$gid;";
   $r=pg_query($db,$q);   
   $q="update routing_ways set y1=st_y(st_startpoint(the_geom))  where gid=$gid;";       
   $r=pg_query($db,$q);   
   $q="update routing_ways set x2=st_x(st_endpoint(the_geom))  where gid=$gid;";       
   $r=pg_query($db,$q);   
   $q="update routing_ways set y2=st_y(st_endpoint(the_geom))  where gid=$gid;";         
   $r=pg_query($db,$q);   
   pg_query($db, "COMMIT"); 
   $ret=$po;
  }
 }
// echo "Teile liefert $ret\n";
 return $ret;
}


function verbinde($db,$osm_id,$source,$target,$x1,$y1,$x2,$y2,$way){
//
// Fuegt den den Weges von source/x1/x2 zu target/x2/y2 way mit der osm_id in den Graphen ein
//
  $query="select max(gid) as m from routing_ways;";
  $result=pg_query($db,$query);
  $e=pg_fetch_array($result);
  $gid=$e['m']+1;
  $q="insert into routing_ways (gid,class_id,osm_id,source,target,waysource,costfactor,hdiff) values ($gid,0,$osm_id,$source,$target,'virtual_area_way',1.0,0);";
  $r=pg_query($db,$q);
  $q="update routing_ways set the_geom=st_transform(ST_GeomFromText('$way',900913),4326) where gid=$gid;";
  $r=pg_query($db,$q);   
  $q="update routing_ways set length=st_length(st_transform(the_geom,900913))*0.66  where gid=$gid;";
  $r=pg_query($db,$q);
  $q="update routing_ways set x1=st_x(st_startpoint(the_geom))  where gid=$gid;";
  $r=pg_query($db,$q);   
  $q="update routing_ways set y1=st_y(st_startpoint(the_geom))  where gid=$gid;";       
  $r=pg_query($db,$q);   
  $q="update routing_ways set x2=st_x(st_endpoint(the_geom))  where gid=$gid;";       
  $r=pg_query($db,$q);   
  $q="update routing_ways set y2=st_y(st_endpoint(the_geom))  where gid=$gid;";         
  $r=pg_query($db,$q);   
}


function winkel($x1,$y1,$x2,$y2,$x3,$y3){
//
// Gibt den Winkel an xy2 zurueck
//
 $w1=rad2deg(atan2(($y2-$y1),($x2-$x1)));if($w1<0){$w1+=360;}
 $w2=rad2deg(atan2(($y3-$y2),($x3-$x2)));if($w2<0){$w2+=360;}
 $w=$w2-$w1;
 if($w>180 ){$w=$w-360;}
 if($w<-180){$w=$w+360;}
// echo "Winkel: $w1 -> $w2 = $w \n";
 return $w;
}


// ------------------------- hier gehts los

 $db=pg_pconnect('host=localhost port=5432 dbname=osm user=osm password=osm') or die("Keine DB, sonst alles in Ordnung");
//
// Letzten Source/targe abholen, brauchen wir dann später
//
 $query="select max(source) as m from routing_ways;";
 $result=pg_query($db,$query);
 $e=pg_fetch_array($result);
 $po=$e['m'];
 $query="select max(target) as m from routing_ways;";
 $result=pg_query($db,$query);
 $e=pg_fetch_array($result);
 if($e['m']>$po){$po=$e['m'];}
 $po++;

//
// Alle Polygone holen 
//
 $placecount=$ringcount=0;
// $query1="select osm_id,st_astext(way) as way from osm_pedestrian where osm_id=-1951581;";
// $query1="select osm_id,st_astext(way) as way from osm_pedestrian where st_contains((select way from osm_polygon where osm_id=-62428 order by way_area desc limit 1),way);";
// $query1="select osm_id,st_astext(way) as way from osm_pedestrian where osm_id='-1934891'";
// $query1="select osm_id,st_astext(way) as way from osm_pedestrian where osm_id='45697424'";
 $query1="select osm_id,st_astext(way) as way from osm_pedestrian;";
// $query1="select osm_id,st_astext(way) as way from osm_pedestrian where (osm_id='-152653');";
//  $query1="select osm_id,st_astext(way) as way from osm_pedestrian where (osm_id='33364851 ');";
 
 $result1=pg_query($db,$query1);
 while($r1=pg_fetch_array($result1)) {
  $placecount++;
  $ringcount=0;
  $osm_id=$r1['osm_id'];
  $flaeche=$r1['way'];  
  $query5="select ipoint from osm_pedestrian where osm_id='$osm_id' limit 1;";
  $result5=pg_query($db,$query5);   
  $r5=pg_fetch_array($result5);     
  $ipo=split(',',$r5['ipoint']);    
  if((count($ipo)>0)&&(count($ipo)<80)){         
//
// von jedem Polygon alle Ringe holen
//
  $query2="select st_astext(st_transform((o.w).geom::geometry,900913)) as polygone,(o.w).path as path from (select ST_DumpRings((g.dump).geom) as w from (select st_dump(way) as dump from osm_pedestrian where osm_id='".$osm_id."') as g) as o;";
  $result2=pg_query($db,$query2);
  $n=0;$i=0;
  while($r2=pg_fetch_array($result2)) {
//
// Aus jedem Ring alle Punkte holen
//
   $polygone=$r2['polygone'];
   if($r2['path']=='{0}'){$exterior=1;}else{$exterior=0;}
   $ringcount++;
   echo "$placecount/$ringcount  (e:$exterior von  $osm_id ".$r2['path'].")\n";
   $p2=split(',',str_ireplace(array('POLYGON','(',')'),'',$polygone));
   for($j=0;$j<count($p2);$j++){
    $p3=split(' ',$p2[$j]);
    $x[$j]=$p3[0];
    $y[$j]=$p3[1];
   }
//
// An jedem Punkt den vorhandenen Graphen auftrennen, oder sich einen neuen
// Knoten geben lassen
//
   for($n=0;$n<count($p2)-1;$n++){
    $stueckel=0;$wi=0;
    if(($n==0)||($n==count($p2)-1)){
     $wi=winkel($x[count($p2)-2],$y[count($p2)-2],$x[0],$y[0],$x[1],$y[1]);
     $rewi=rad2deg(atan2($y[0]-$y[count($p2)-2],$x[0]-$x[count($p2)-2]));
     $liwi=rad2deg(atan2($y[0]-$y[1],$x[0]-$x[1]));              
    }else{
     $wi=winkel($x[$n-1],$y[$n-1],$x[$n],$y[$n],$x[$n+1],$y[$n+1]);
     $rewi=rad2deg(atan2($y[$n]-$y[$n-1],$x[$n]-$x[$n-1]));
     $liwi=rad2deg(atan2($y[$n]-$y[$n+1],$x[$n]-$x[$n+1]));
    }
    if(($exterior==1) && ($wi>0)){$stueckel=1;}
    if(($exterior==0) && ($wi>0)){$stueckel=1;}
//    echo "Versuche Teilung mit $po (n=$n stc=$stueckel ext=$exterior wi=$wi)\n";
    $w[$n]=teilweg($db,$x[$n],$y[$n],$po,$stueckel);
    if($w[$n]>0){
     if($w[$n]==$po){$po++;}
     $ex[$i]=$x[$n];$ey[$i]=$y[$n];$ew[$i]=$w[$n];
     $rewinkel[$i]=$rewi;
     $liwinkel[$i]=$liwi;
//     echo "Nehme ".$w[$n]." in Graphen auf ReWi=$rewi LiWi=$liwi\n";
     $i++;
    }
   }
  }
 }
//
// Notbremse: nicht mehr als 500 neue Knoten und 100 Hindernisse
//
  
  
  if(($i<600)&&($ringcount<120)){
  $query5="select ipoint from osm_pedestrian where osm_id='$osm_id' limit 1;";
  $result5=pg_query($db,$query5);   
  $r5=pg_fetch_array($result5);
  $ipo=split(',',$r5['ipoint']);
//
// von jedem zu jedem angelegten Punkt eine Linie ziehen, falls die innerhalb des Polygons liegt
// 
    for($k=0;$k<$i;$k++){
     $kisipo=0;
     for($ii=0;$ii<count($ipo);$ii++){if($ipo[$ii]==$ew[$k]){$kisipo=1;}}
     for($j=0;$j<$k;$j++){
      $jisipo=0;
      for($ii=0;$ii<count($ipo);$ii++){if($ipo[$ii]==$ew[$j]){$jisipo=1;}}
      $testweg="LINESTRING(".$ex[$j]." ".$ey[$j].",".$ex[$k]." ".$ey[$k].")";
      $winkel=rad2deg(atan2($ey[$k]-$ey[$j],$ex[$k]-$ex[$j]));
      $winkel2=rad2deg(atan2($ey[$j]-$ey[$k],$ex[$j]-$ex[$k]));
      $query3="select st_astext(st_difference(st_geometryfromtext('$testweg',900913),st_geometryfromtext('$flaeche',900913))) as t;";
      $result3=pg_query($db,$query3);
      $e3=pg_fetch_array($result3);
      if($e3['t']=='GEOMETRYCOLLECTION EMPTY'){
//       echo "Verbinde ".$ew[$j]." mit ".$ew[$k]." Winkel: $winkel / $winkel2  (rw=".$rewinkel[$j]." lw=".$liwinkel[$j]." zu rw=".$rewinkel[$k]." lw=".$liwinkel[$k].")\n";
       $qliegtrechts=$qliegtlinks=0;
       $zliegtrechts=$zliegtlinks=0;
       if(($winkel>=0)&&($rewinkel[$j]>=0)&&($winkel< $rewinkel[$j]    )){$qliegtrechts=1;}
       if(($winkel>=0)&&($rewinkel[$j] <0)&&($winkel> $rewinkel[$j]+180)){$qliegtrechts=1;}
       if(($winkel <0)&&($rewinkel[$j]>=0)&&($winkel> $rewinkel[$j]-180)){$qliegtrechts=1;}
       if(($winkel <0)&&($rewinkel[$j] <0)&&($winkel< $rewinkel[$j]    )){$qliegtrechts=1;}
       
       if(($winkel>=0)&&($liwinkel[$j]>=0)&&($winkel> $liwinkel[$j]    )){$qliegtlinks=1;}
       if(($winkel>=0)&&($liwinkel[$j] <0)&&($winkel< $liwinkel[$j]+180)){$qliegtlinks=1;}
       if(($winkel <0)&&($liwinkel[$j]>=0)&&($winkel< $liwinkel[$j]-180)){$qliegtlinks=1;}
       if(($winkel <0)&&($liwinkel[$j] <0)&&($winkel> $liwinkel[$j]))    {$qliegtlinks=1;}
       
       if(($winkel2>=0)&&($rewinkel[$k]>=0)&&($winkel2< $rewinkel[$k]    )){$zliegtrechts=1;}
       if(($winkel2>=0)&&($rewinkel[$k] <0)&&($winkel2> $rewinkel[$k]+180)){$zliegtrechts=1;}
       if(($winkel2 <0)&&($rewinkel[$k]>=0)&&($winkel2> $rewinkel[$k]-180)){$zliegtrechts=1;}
       if(($winkel2 <0)&&($rewinkel[$k] <0)&&($winkel2< $rewinkel[$k]    )){$zliegtrechts=1;}
       
       if(($winkel2>=0)&&($liwinkel[$k]>=0)&&($winkel2> $liwinkel[$k]    )){$zliegtlinks=1;}
       if(($winkel2>=0)&&($liwinkel[$k] <0)&&($winkel2< $liwinkel[$k]+180)){$zliegtlinks=1;}
       if(($winkel2 <0)&&($liwinkel[$k]>=0)&&($winkel2< $liwinkel[$k]-180)){$zliegtlinks=1;}
       if(($winkel2 <0)&&($liwinkel[$k] <0)&&($winkel2> $liwinkel[$k]))    {$zliegtlinks=1;}
       
       if((($qliegtlinks==0)||($qliegtrechts==0)||($jisipo==1)) && (($zliegtlinks==0)||($zliegtrechts==0)||($kisipo==1))){
//        echo "-------- Verbinde  lili=$qliegtlinks lire=$qliegtrechts   zu  lili=$zliegtlinks lire=$zliegtrechts\n";
        verbinde($db,$osm_id,$ew[$j],$ew[$k],$ex[$j],$ey[$j],$ex[$k],$ey[$k],$testweg);
       }
      }
     }
    }
//
// Danach die "weiteren wichtigen Punkte" holen, aber erst schaun,
// ob die schon abgehandelt wurden
//
   $query5="select ipoint from osm_pedestrian where osm_id='$osm_id' limit 1;";
   $result5=pg_query($db,$query5);   
   $r5=pg_fetch_array($result5);
   $ipo=split(',',$r5['ipoint']);
   for($k=0;$k<count($ipo);$k++){
    $tuwas=1;
    for($j=0;$j<$i;$j++){if($ipo[$k]==$ew[$j]){$tuwas=0;}}
    if($tuwas){
     $query6="select st_x(st_transform(ST_SetSRID(st_makepoint(x,y),4326),900913)),st_y(st_transform(ST_SetSRID(st_makepoint(x,y),4326),900913)) from (select x1 as x,y1 as y,source as p from routing_ways where source=".$ipo[$k]." union  select x2 as x,y2 as y,target as p from routing_ways where target=".$ipo[$k]." ) as l limit 1;";
     $result6=pg_query($db,$query6);
     $r6=pg_fetch_array($result6);
     $ex[$i]=$r6['st_x'];
     $ey[$i]=$r6['st_y'];
     $ew[$i]=$ipo[$k];
     for($j=0;$j<$i;$j++){
      $jisipo=0;
      for($ii=0;$ii<count($ipo);$ii++){if($ipo[$ii]==$ew[$j]){$jisipo=1;}}
      $testweg="LINESTRING(".$ex[$j]." ".$ey[$j].",".$ex[$i]." ".$ey[$i].")";
      $query3="select st_astext(st_difference(st_geometryfromtext('$testweg',900913),st_geometryfromtext('$flaeche',900913))) as t;";
      $result3=pg_query($db,$query3);
      $e3=pg_fetch_array($result3);
      if($e3['t']=='GEOMETRYCOLLECTION EMPTY'){
       verbinde($db,$osm_id,$ew[$j],$ew[$i],$ex[$j],$ey[$j],$ex[$i],$ey[$i],$testweg);
      }
     }
     $i++;
    }
   }
   $query5="select st_x(way) as x,st_y(way) as y  from osm_point as P where (layer is null or layer='0') and  (tags->'ford'='yes' or tags->'entrance'!='' or railway like '%entrance' or amenity in ('watering_place','hunting_stand','shelter','drinking_water') or tourism in ('viewpoint','attraction') or man_made='water_well' or nature in ('spring','cave_entrance','col','saddle','peak','mountain_pass') ) and st_contains((select way from osm_pedestrian where osm_id='".$osm_id."' limit 1),P.way) and not exists (select  osm_id from osm_line as L where highway is not null and P.way && L.way and st_distance(L.way,P.way)<2);";
   $result5=pg_query($db,$query5);
   while($r5=pg_fetch_array($result5)){
    $po++;
    $ex[$i]=$r5['x'];
    $ey[$i]=$r5['y'];
    $ew[$i]=$po;
    $query5="update osm_pedestrian set ipoint=ipoint||',$po' where osm_id='".$osm_id."';";
    pg_query($db,$query5);
    for($j=0;$j<$i;$j++){
     $testweg="LINESTRING(".$ex[$j]." ".$ey[$j].",".$ex[$i]." ".$ey[$i].")";
     $query3="select st_astext(st_difference(st_geometryfromtext('$testweg',900913),st_geometryfromtext('$flaeche',900913))) as t;";
     $result3=pg_query($db,$query3);
     $e3=pg_fetch_array($result3);
     if($e3['t']=='GEOMETRYCOLLECTION EMPTY'){
      verbinde($db,$osm_id,$ew[$j],$ew[$i],$ex[$j],$ey[$j],$ex[$i],$ey[$i],$testweg);
     }   
    }
    $i++;
   }
  }else{
   echo " -- Notbremse wegen $i Punkten und $ringcount Ringen\n";
  }
 }
 exit;  
?>