Funcion Oracle Spatial a Posgis

Buen dia a todos

Alguien me puede ayudar para pasar la funcion de lineas abajo de Oracle
spatial a postgis necesito un resultado similiar en mi BD esto lo vi en
este blog

https://spatialdbadvisor.com/oracle_spatial_tips_tricks/307/computing-cardinal-directions-to-nearby-geometries

1. CREATE OR REPLACE
2. FUNCTION CardinalDirection(p_start_point IN sdo_geometry,
3. p_end_point IN sdo_geometry,
4. p_tolerance IN NUMBER DEFAULT 0.05,
5. p_projected IN INTEGER DEFAULT 1)
6. RETURN varchar2 deterministic
7. AS
8. c_PI CONSTANT NUMBER := acos(-1);
9. v_bearing NUMBER;
10. v_tilt NUMBER;
11. v_vertex1 mdsys.vertex_type;
12. v_vertex2 mdsys.vertex_type;
13. FUNCTION ST_Azimuth(p_dE1 IN NUMBER, p_dN1 IN NUMBER,
14. p_dE2 IN NUMBER, p_dN2 IN NUMBER)
15. RETURN NUMBER
16. IS
17. dBearing NUMBER;
18. dEast NUMBER;
19. dNorth NUMBER;
20. BEGIN
21. IF (p_dE1 IS NULL OR p_dN1 IS NULL OR
22. p_dE2 IS NULL OR p_dE2 IS NULL ) THEN
23. RETURN NULL;
24. END IF;
25. IF ( (p_dE1 = p_dE2) AND (p_dN1 = p_dN2) ) THEN
26. RETURN NULL;
27. END IF;
28. dEast := p_dE2 - p_dE1;
29. dNorth := p_dN2 - p_dN1;
30. IF ( dEast = 0 ) THEN
31. dBearing := CASE WHEN ( dNorth < 0 )
32. THEN c_PI
33. ELSE 0
34. END;
35. ELSE
36. dBearing := -aTan(dNorth / dEast) + c_PI / 2;
37. END IF;
38. RETURN CASE WHEN ( dEast < 0 )
39. THEN dBearing + c_PI
40. ELSE dBearing
41. END;
42. END ST_Azimuth;
43. BEGIN
44. IF (p_start_point IS NULL OR p_end_point IS NULL) THEN
45. RETURN NULL;
46. END IF;
47. -- If Geodetic....
48. IF (NVL(ABS(p_projected),1) = 0 ) THEN
49. SDO_UTIL.BEARING_TILT_FOR_POINTS(
50. p_start_point,
51. p_end_point,
52. p_tolerance,v_bearing,v_tilt);
53. ELSE
54. v_vertex1 := sdo_util.getVertices(p_start_point)(1);
55. v_vertex2 := sdo_util.getVertices(p_end_point)(1);
56. v_bearing := ST_Azimuth(v_vertex1.X, v_vertex1.Y,
57. v_vertex2.X, v_vertex2.Y);
58. END IF;
59. v_bearing := round(v_bearing * (180.0 / acos(-1)),2);
60. RETURN CASE
61. WHEN v_bearing IS NULL THEN 'NULL'
62. WHEN /* Cardinal Point 0 */ v_bearing BETWEEN 348.75 AND
0.0
63. OR v_bearing BETWEEN 0.0 AND
11.25 THEN 'N'
64. WHEN /* Cardinal Point 22.5 */ v_bearing BETWEEN 11.25 AND
33.75 THEN 'NNE'
65. WHEN /* Cardinal Point 45 */ v_bearing BETWEEN 33.75 AND
56.25 THEN 'NE'
66. WHEN /* Cardinal Point 67.5 */ v_bearing BETWEEN 56.25 AND
78.75 THEN 'ENE'
67. WHEN /* Cardinal Point 90 */ v_bearing BETWEEN 78.75 AND
101.25 THEN 'E'
68. WHEN /* Cardinal Point 112.5 */ v_bearing BETWEEN 101.25 AND
123.75 THEN 'ESE'
69. WHEN /* Cardinal Point 135 */ v_bearing BETWEEN 123.75 AND
146.25 THEN 'SE'
70. WHEN /* Cardinal Point 157.5 */ v_bearing BETWEEN 146.25 AND
168.75 THEN 'SSE'
71. WHEN /* Cardinal Point 180 */ v_bearing BETWEEN 168.75 AND
191.25 THEN 'S'
72. WHEN /* Cardinal Point 202.5 */ v_bearing BETWEEN 191.25 AND
213.75 THEN 'SSW'
73. WHEN /* Cardinal Point 225 */ v_bearing BETWEEN 213.75 AND
236.25 THEN 'SW'
74. WHEN /* Cardinal Point 247.5 */ v_bearing BETWEEN 236.25 AND
258.75 THEN 'WSW'
75. WHEN /* Cardinal Point 270 */ v_bearing BETWEEN 258.75 AND
281.25 THEN 'W'
76. WHEN /* Cardinal Point 292.5 */ v_bearing BETWEEN 281.25 AND
303.75 THEN 'WNW'
77. WHEN /* Cardinal Point 315 */ v_bearing BETWEEN 303.75 AND
326.25 THEN 'NW'
78. WHEN /* Cardinal Point 337.5 */ v_bearing BETWEEN 326.25 AND
348.75 THEN 'NNW'
79. ELSE to_char(v_bearing)
80. END;
81. END CardinalDirection;
82. /
83. SHOW errors

--
José Mercedes Venegas Acevedo
cel Mov RPC 964185205