В PostGIS версии 1.5.3 ST_Azimuth () в качестве входных данных требовалось две точки - насколько я могу судить, для всех версий - которые я просто извлек из строки, используя ST_Line_Interpolate_Point ( ) , вот так:
ST_Azimuth(
ST_Line_Interpolate_Point(wkb_geometry, 0.2), /* point at 20% length */
ST_Line_Interpolate_Point(wkb_geometry, 0.8) /* point at 80% length */
)/(2*pi())*360 as angle /* radians to degrees */