Changes between Version 5 and Version 6 of UsersWikiplpgsqlfunctions

Timestamp:

04/19/09 00:39:50 (16 years ago)

Author:

robe

Comment:

--

UsersWikiplpgsqlfunctions

@@ -7 +7 @@
  * '''Find UTM (WGS84)  SRID for a point (in any SRID)'''
 
- {{{
+{{{
+#!sql
  -- Function: utmzone(geometry)
  -- DROP FUNCTION utmzone(geometry);
@@ -20 +21 @@
 
  BEGIN
-     geomgeog:=transform($1,4326);
+     geomgeog:= ST_Transform($1,4326);
 
-     IF (y(geomgeog))>0 THEN
+     IF (ST_Y(geomgeog))>0 THEN
         pref:=32600;
      ELSE
@@ -28 +29 @@
      END IF;
 
-     zone:=floor((x(geomgeog)+180)/6)+1;
+     zone:=floor((ST_X(geomgeog)+180)/6)+1;
 
      RETURN zone+pref;
@@ -38 +39 @@
  * '''Find Length of Day for a given location and date'''
 {{{
+#!sql
  CREATE OR REPLACE FUNCTION daylength(geometry, date)
- RETURNS double precision AS 'select 24*(ACOS(1-(1-TAN(radians(Y(transform($1,4326))))*TAN(.409088 * COS(.0172024 * extract(doy from $2)))))/PI());'
+ RETURNS double precision AS 
+$$ 
+SELECT 24*(ACOS(1-(1-TAN(radians(Y(transform($1,4326))))*TAN(.409088 * COS(.0172024 * extract(doy from $2)))))/PI());
+$$
  LANGUAGE 'sql';
 }}}
@@ -45 +50 @@
  * '''Rotate Geometry around a point(geom,x,y)'''
 {{{
+#!sql
 CREATE OR REPLACE FUNCTION RotateAtPoint(geometry, double precision, double precision, double precision)
- RETURNS geometry AS 'SELECT translate( rotate( translate( $1, -1*$2, -1*$3), $4), $2, $3)'
+ RETURNS geometry AS 
+$$
+SELECT ST_Translate( ST_Rotate( ST_Translate( $1, -1*$2, -1*$3), $4), $2, $3)
+$$
    LANGUAGE 'sql';
 }}}
@@ -52 +61 @@
  * '''Create Ellipse (x,y,rx,ry,rotation,#of segments in 1/4 of ellipse)'''
 {{{
+#!sql
 CREATE OR REPLACE FUNCTION Ellipse(double precision, double precision, double precision, double precision, double precision, integer)
  RETURNS geometry AS
- 'SELECT translate( rotate( scale( buffer(makepoint(0,0), 0.5, $6), $3, $4), $5), $1, $2)'
+ $$
+   SELECT ST_Translate( ST_Rotate( ST_Scale( ST_Buffer(ST_Point(0,0), 0.5, $6), $3, $4), $5), $1, $2)
+$$
  LANGUAGE 'sql';
 }}}