Oh, yes, sorry... I forgot to do the intersection on the geom field! below
works perfectly, thanks!
select ST_Intersection(TABLE1.geom,
ST_GeomFromText('POLYGON((31.7 -13.1, 31.6
-13.04,31.96 -13.04,31.9 -13.1,31.7
-13.1))',4326) ) as geom,
a_code as code,ST_Area(ST_Transform(ST_Intersection(TABLE1.geom,
ST_GeomFromText('POLYGON((31.7 -13.1, 31.6
-13.04,31.96 -13.04,31.9 -13.1,31.7
-13.1))',4326) ),32736))*0.0001::double precision as area from TABLE1
where
ST_Intersection(TABLE1.geom,
ST_GeomFromText('POLYGON((31.7 -13.1, 31.6
-13.04,31.96 -13.04,31.9 -13.1,31.7
-13.1))',4326) )
--
Sent from: http://osgeo-org.1560.x6.nabble.com/GDAL-Dev-f3742093.html
_______________________________________________
gdal-dev mailing list
gdal-dev@lists.osgeo.org
https://lists.osgeo.org/mailman/listinfo/gdal-dev
