################################################################################
### Creating networks from imported OSM data (8 / 11)
################################################################################
### Description:
### - This script creates a relation containing junctions between two sections
### belonging to two Streets.
### (the number 1 in the relation name stands for incoming, the number 2 means
### outgoing)
### Please, see tu_shp_import.sh for further details. 
###
### Preconditions:
### - activated necessary algebras
### - existing open database with successfully imported osm data
### - FullSectionsView-relation
###
### Postconditions:
### - twoSecJuncBtwTwoSts-relation
### twoSecJuncBtwTwoSts: rel{Road1: int, Pos1: double, Road2: int,
###                          Pos2: double, JunctionType: int}
###
### Author:
### - Thomas Uchdorf, thomas.uchdorf(at)fernuni-hagen.de
################################################################################

#      |
#      |
#      |
#      v
#----->xxxxxx
#      x
#      x
#      x
#      x
let TwoSecInZeroSecOut11JuncTmp =
   ((NodesRel feed
       FullSectionsView feed
       sortby [SectionId]
       filter [(.StreetEndNodeId = .EndNodeId) and
          not(.StreetStartNodeId = .StreetEndNodeId)]
       project [StreetId,SectionId,EndNodeId,StreetData,SectionData,
          OneWay] {a}
       hashjoin [NodeId, endNodeId_a,99997])
      FullSectionsView feed
      sortby [SectionId]
      filter [(.StreetEndNodeId = .EndNodeId) and
         not(.StreetStartNodeId = .StreetEndNodeId)]
      project [StreetId,SectionId,EndNodeId,StreetData,SectionData,
         OneWay] {b}
      hashjoin [NodeId,EdNodeId_b,99997]
      filter [.SectionId_a < .SectionId_b])
   projectextend [; 
      Road1: .StreetId_a,
      Pos1: atpoint(.StreetData_a,.Node,TRUE),
      Road2: .StreetId_b,
      Pos2: atpoint(.StreetData_b,.Node,TRUE),
      JunctionType:
         getconnectivitycode(1,2,0,0,.OneWay_a = 1,.OneWay_b = 1,FALSE,FALSE),
      Sec1: .SectionId_a,
      Sec2: .SectionId_b,
      Node: .Node,
      SecData1: .SectionData_a,
      SecData2: .SectionData_b]
   filter [not(.Road1 = .Road2)]
   consume;

#      ^
#      |
#      |
#      |
#----->xxxxxx
#      x
#      x
#      x
#      x
let OneSecInOneSecOut12JuncTmp =
   ((NodesRel feed
       FullSectionsView feed
       sortby [SectionId]
       filter [(.StreetEndNodeId = .EndNodeId) and
          not(.StreetStartNodeId = .StreetEndNodeId)]
       project [StreetId,SectionId,EndNodeId,StreetData,SectionData,
          OneWay] {a}
       hashjoin [NodeId, EndNodeId_a,99997])
      FullSectionsView feed
      sortby [SectionId]
      filter [(.StreetStartNodeId = .StartNodeId) and
         not(.StreetStartNodeId = .StreetEndNodeId)]
      project [StreetId,SectionId,StartNodeId,StreetData,SectionData,
         OneWay] {b}
      hashjoin [NodeId,StartNodeId_b,99997]
      filter [.SectionId_a < .SectionId_b])
   projectextend [; 
      Road1: .StreetId_a,
      Pos1: atpoint(.StreetData_a,.Node,TRUE),
      Road2: .StreetId_b,
      Pos2: atpoint(.StreetData_b,.Node,TRUE),
      JunctionType:
         getconnectivitycode(1,1,0,0,.OneWay_a = 1,.OneWay_b = 1,FALSE,FALSE),
      Sec1: .SectionId_a,
      Sec2: .SectionId_b,
      Node: .Node,
      SecData1: .SectionData_a,
      SecData2: .SectionData_b]
   filter [not(.Road1 = .Road2)]
   consume;

#      |
#      |
#      |
#      v
#<-----xxxxxx
#      x
#      x
#      x
#      x
let OneSecInOneSecOut21JuncTmp =
   ((NodesRel feed
       FullSectionsView feed
       sortby [SectionId]
       filter [(.StreetStartNodeId = .StartNodeId) and
          not(.StreetStartNodeId = .StreetEndNodeId)]
       project [StreetId,SectionId,StartNodeId,StreetData,SectionData,
          OneWay] {a}
       hashjoin [NodeId, StartNodeId_a,99997])
      FullSectionsView feed
      sortby [SectionId]
      filter [(.StreetEndNodeId = .EndNodeId) and
         not(.StreetStartNodeId = .StreetEndNodeId)]
      project [StreetId,SectionId,EndNodeId,StreetData,SectionData,
         OneWay] {b}
      hashjoin [NodeId,EndNodeId_b,99997]
      filter [.SectionId_a < .SectionId_b])
   projectextend [; 
      Road1: .StreetId_a,
      Pos1: atpoint(.StreetData_a,.Node,TRUE),
      Road2: .StreetId_b,
      Pos2: atpoint(.StreetData_b,.Node,TRUE),
      JunctionType:
         getconnectivitycode(2,2,0,0,.OneWay_a = 1,.OneWay_b = 1,FALSE,FALSE),
      Sec1: .SectionId_a,
      Sec2: .SectionId_b,
      Node: .Node,
      SecData1: .SectionData_a,
      SecData2: .SectionData_b]
   filter [not(.Road1 = .Road2)]
   consume;

#      ^
#      |
#      |
#      |
#<-----xxxxxx
#      x
#      x
#      x
#      x
let ZeroSecInTwoSecOut22JuncTmp =
   ((NodesRel feed
       FullSectionsView feed
       sortby [SectionId]
       filter [(.StreetStartNodeId = .StartNodeId) and
          not(.StreetStartNodeId = .StreetEndNodeId)]
       project [StreetId,SectionId,StartNodeId,StreetData,SectionData,
          OneWay] {a}
       hashjoin [NodeId, StartNodeId_a,99997])
      FullSectionsView feed
      sortby [SectionId]
      filter [(.StreetStartNodeId = .StartNodeId) and
         not(.StreetStartNodeId = .StreetEndNodeId)]
      project [StreetId,SectionId,StartNodeId,StreetData,SectionData,
         OneWay] {b}
      hashjoin [NodeId,startNodeId_b,99997]
      filter [.SectionId_a < .SectionId_b])
   projectextend [; 
      Road1: .StreetId_a,
      Pos1: atpoint(.StreetData_a,.Node,TRUE),
      Road2: .StreetId_b,
      Pos2: atpoint(.StreetData_b,.Node,TRUE),
      JunctionType:
         getconnectivitycode(2,1,0,0,.OneWay_a = 1,.OneWay_b = 1,FALSE,FALSE),
      Sec1: .SectionId_a,
      Sec2: .SectionId_b,
      Node: .Node,
      SecData1: .SectionData_a,
      SecData2: .SectionData_b]
   filter [not(.Road1 = .Road2)]
   consume;

# Creating the union
let TwoSecJuncBtwTwoSts =
   ((TwoSecInZeroSecOut11JuncTmp feed
     OneSecInOneSecOut12JuncTmp feed
     concat)
    OneSecInOneSecOut21JuncTmp feed
    concat)
   ZeroSecInTwoSecOut22JuncTmp feed
   concat
   consume;