<?php

$db=pg_pconnect('host=localhost port=5432 dbname=osm user=osm password=osm') or die("Keine DB, sonst alles in Ordnung");

// Erst alle Wehre holen
//

$weirquery="select osm_id,way from osm_line where waterway='weir' and way && st_transform(ST_MakeEnvelope(9,46,15,49,4326),900913);";
$resultweir=pg_query($db,$weirquery);
while ($w=pg_fetch_array($resultweir)){
 $osmidw=$w['osm_id'];
 $wayw=$w['way'];

// Dann den ersten kreuzenden Fluss suchen
//

 $streamquery="select osm_id,waterway,way from osm_line where waterway in ('stream','river','canal','ditch','drain','brook') and st_intersects(way,'$wayw') limit 1;";
 $resultstream=pg_query($db,$streamquery);
 if($resultstream){
  $s=pg_fetch_array($resultstream);
  if($s){
   $osmids=$s['osm_id'];
   $ways=$s['way'];
   $wways=$s['waterway'];

// Zu diesem Fluss die Kreuzungsrichtung ermitteln und nachsehen, wo dieser Fluss vom Wehr geteilt wird
//

   $crossquery="select ST_LineCrossingDirection('$wayw','$ways') as lc, ST_Line_Locate_Point('$ways',(select ST_ClosestPoint('$ways', '$wayw'))) as llp,st_length('$ways') as l;";
   $resultcrossing=pg_query($db,$crossquery);
   $c=pg_fetch_array($resultcrossing);
   $cross=$c['lc'];
   $cpoint=$c['llp'];
   if($cpoint<0.5) {$len=$cpoint*$c['l'];}
   if($cpoint>=0.5){$len=(1-$cpoint)*$c['l'];}

// Kreuzung: -1 heisst, das Wehr kommt von links, beim Wehr gehts LinksRunter
//           +1 Wehr kommt von rechts, Rechts gehts Runter
//           0: keine Kreuzung, +-2, +-3 heisst Mehrfachkreuzung
//
   $dir="??";
   if($cross=="-1"){$dir="LR";}
   if($cross=="1"){$dir="RR";}
   $ae="MID";
   if(($cpoint<0.5)&&($len<0.05)){$ae="START";}
   if(($cpoint>=0.5)&&($len<0.05)){$ae="END";}
   echo "weir $osmidw $wways $osmids $cross $dir $ae\n";
  }else{
   echo "weir $osmidw NC\n";
 }}
}


?>