I want to calculate the angle of two lines of a routing network exactly at the "touching point" (the topological vertex).
My previous attempt was using the
ST_ANGLE() function of PostGIS, which gave me a "wrong" result of 76° for the point 4060 between line [16543,4060] and line [4060,3382], since
ST_ANGLE() takes the beginning and end node of the input lines to measure the angle betwen those two.
ST_LineInterpolatePoint to generate auxiliary points close to the node isn't that easy also due to the different orientations of the edges.