Compute partial polygon area

Stephen Woodbridge woodbri at SWOODBRIDGE.COM
Sat Dec 31 15:53:14 PST 2005


Krung Saengpole wrote:
> Hi all,
>  
> I'd like to compute partial area of polygon that intersects with another 
> polygon. Is there a way to do with PHPMapscript? For example, I want to 
> know how many square kilometers of forests that reside in the area of 
> 50km-radius circled area. Or where can I get to the example?

If you use PostGIS this is a straight forward SQL query like:

select polyid, area2d(intersection(the_geom, circleAsPolygon(point(x, 
y), radius, npoints)) as area from mytable where the_geom && 
expand(point(x, y), radius);

where the_geom is the geometry column for the polygons
and x,y is the center of your circle radius = 50000
and npoints is the number of points to generate for the polygonal circle.

you would need to write a function circleAsPolygon() in pgsql, but that 
is pretty trivial and most of us that would with PostGIS could post that 
if you needed us to.

-Steve



More information about the MapServer-users mailing list