[postgis-users] I need to select all featcher touch one polygon

Paul Ramsey pramsey at cleverelephant.ca
Thu Aug 7 12:45:38 PDT 2008


SELECT a.the_geom
FROM "IGFPLAN" a, "IGFPLAN" b
WHERE
  a."Id" = '2' AND
  b."Id" = '2' AND
  ST_Intersects(a.the_geom, b.the_geom)

Leaving your columns in with case preserved really makes writing SQL picky, IMO.

P.

On Thu, Aug 7, 2008 at 10:55 AM, eehab hamzeh <eehab40 at hotmail.com> wrote:
>
> Hello every body,
>
> can any body help me to write sql statment that select all the polygons that
> touch polygon in the same feature class.
> i try to do it in as below and find all polygon that are within 100 m from
> the polygon 2.
>
> SELECT the_geom FROM "public"."IGFPLAN"
> WHERE "Id" = '2' and distance(the_geom, the_geom) < 100
>
>
> Thanks
>
> ihab
>
>
> ________________________________
>> From: postgis-users-request at postgis.refractions.net
>> Subject: postgis-users Digest, Vol 70, Issue 6
>> To: postgis-users at postgis.refractions.net
>> Date: Wed, 6 Aug 2008 12:01:19 -0700
>>
>> Send postgis-users mailing list submissions to
>> postgis-users at postgis.refractions.net
>>
>> To subscribe or unsubscribe via the World Wide Web, visit
>> http://postgis.refractions.net/mailman/listinfo/postgis-users
>> or, via email, send a message with subject or body 'help' to
>> postgis-users-request at postgis.refractions.net
>>
>> You can reach the person managing the list at
>> postgis-users-owner at postgis.refractions.net
>>
>> When replying, please edit your Subject line so it is more specific
>> than "Re: Contents of postgis-users digest..."
>>
>>
>> Today's Topics:
>>
>> 1. Newbie Question (Adam Vande More)
>> 2. Re: ST_Intersects (Paul Ramsey)
>> 3. Re: Newbie Question (Pedro Doria Meunier)
>> 4. Re: Newbie Question (Paul Ramsey)
>> 5. Re: Newbie Question (Pedro Doria Meunier)
>> 6. Update Error (Donald Ijams)
>> 7. RE: Update Error (Obe, Regina)
>> 8. Still cannot display postgis layers using ms4w
>> (kreshna_iceheart at yahoo.com)
>> 9. RE: Still cannot display postgis layers using ms4w
>> (Paragon Corporation)
>> 10. RE: Still cannot display postgis layers using ms4w< br>>
>> (kreshna_iceheart at yahoo.com)
>> 11. Configuration and getting started issues (Khavindra Sivenandan)
>> 12. area calculation problem (danny)
>> 13. filing the holes in the swiss cheese? (Burgholzer,Robert)
>> 14. Re: filing the holes in the swiss cheese? (Paul Ramsey)
>>
>>
>> ----------------------------------------------------------------------
>>
>> Message: 1
>> Date: Tue, 05 Aug 2008 14:24:44 -0500
>> From: Adam Vande More <adam at imedmobility.com>
>> Subject: [postgis-users] Newbie Question
>> To: postgis-users at postgis.refractions.net
>> Message-ID: <4898A8FC.40201 at imedmobility.com>
>> Content-Type: text/plain; charset=ISO-8859-1; format=flowed
>>
>> Hello,
>>
>> Perhaps my questions were addressed earlier or in some documentation I
>> missed. If so, I apologize in advance as I am completely new to GIS and
>> PostGIS. Judging f rom the documentation I have read, my needs are
>> extremely simple vs PostGIS capabilities however I am unsure if the
>> functionality of what I want actually exists.
>>
>> I deal with areas in the US, territories that are regional in nature but
>> not bounded by common things eg state lines, water ways, counties, etc.
>> They are irregular shapes of varying land area. I also have sets of
>> latitude and longitude points. I would like some type of system where I
>> enter in a perimeter of territories, and find out whether a particular
>> latitude and longitude point lies within that perimeter.
>>
>> My understanding of this is that I would need to either manually create
>> a GML type file with the perimeters in it, or use some type of map
>> authoring software to accomplish this goal. Are my understandings
>> correct and if so is there a recommendation for a particular piece of
>> software which can fulfill this?
>>
>> Also, are the libraries/functions available for determining whether a
>> point is in a perimeter?
>>
>> Much appreciated,
>>
>> --
>> Adam Vandemore
>> Systems Administrator
>> IMED Mobility
>> (605) 498-1610
>>
>>
>>
>> ------------------------------
>>
>> Message: 2
>> Date: Tue, 5 Aug 2008 12:39:30 -0700
>> From: "Paul Ramsey" <pramsey at cleverelephant.ca>
>> Subject: Re: [postgis-users] ST_Intersects
>> To: pcreso at pcreso.com, "PostGIS Users Discussion"
>> <postgis-users at postgis.refractions.net>
>> Message-ID:
>> <30fe546d0808051239v54067bcw6335074fe89d2fa6 at mail.gmail.com>
>> Content-Type: text/plain; charset=ISO-8859-1
>>
>> At core, this is not an indexing problem (of the usual sort). When we
>> can fix the memory leaks on the preparedgeometry code, and get it out
>> to users, it will make this case go much much faster, but until then,
>> it's going to be slow.
>>
>> P.
>>
>> On Tue, Aug 5, 2008 at 11:47 AM, <pcreso at pcreso.com> wrote:
>> > Hi Danny,
>> >
>> > As a maritime PostGIS user I have had issues with large & complex
>> > polygons/multipoygons in PostGIS (& GEOS), but not to the extent you seem
>> > to.
>> >
>> > I'm not an expert at this, but here are some things you might try.
>> >
>> > It looks like it is not using the index on sites, perhaps because
>> > postgres is applying a conversion to text for your constant & the attribute.
>> > Can you try modifying your SQL to
>> >
>> > .... where a.sitecode = 'xxxx'::<same spec as the sitecode column>
>> >
>> > (eg: where a.sitecode = 'xxxx'::varchar(10) if sitecode is varchar(10) )
>> >
>> > To see if that enables the index to work better.
>> >
>> > Alternatively, perhaps you could create a new index on a
>> > (sitecode::text)
>> > Though I'm not sure if that's possible....
>> >
>> >
>> > Where I have large complex polygons, I have found tiling them really
>> > enables the spatial index to become useful. I generate a grid in postgis,
>> > build a table of the intersections of the grid & base layer & query against
>> > that.
>> >
>> > Essentially an index helps with long (deep) tables, ie, lots of small
>> > records, much more than fat (broad) tables, ie, few large records. Tiling
>> > geometries converts fat to deep, making things much faster.
>> >
>> >
>> > Cheers,
>> >
>> > Brent Wood
>> >
>> >
>> > --- On Tue, 8/5/08, danny <whatevar89 at gmail.com> wrote:
>> >
>> >> From: danny <whatevar89 at gmail.com>
>> >> Subject: Re: [postgis-users] ST_Intersects
>> >> To: "PostGIS Users Discussion" <postgis-users at postgis.refractions.net>
>> >> Date: Tuesday, August 5, 2008, 11:14 PM
>> >> Guido, it's a bufferzone around the contour of europe.
>> >>
>> >> On Tue, Aug 5, 2008 at 1:02 PM, Guido Lemoine
>> >> <guido.lemoine at jrc.it> wrote:
>> >>
>> >> > Erik,
>> >> >
>> >> > I don't know why european_waters would be a single
>> >> complex polygon?
>> >> > Do you mean it is a MULTIPOLYGON in Postgis? If so, it
>> >> is easy to
>> >> > break it apart in single polygons. If it includes
>> >> rivers, I can't imagine
>> >> > why it
>> >> > is a polygon in the first place. I would expect a
>> >> LINESTRING (or
>> >> > MULTILINESTRING in your case).
>> >> >
>> >> > This explains, of course, w hy performance is not much
>> >> improved. The &&
>> >> > operation does not really affect that.
>> >> >
>> >> > GL
>> >> >
>> >> >
>> >> >
>> >> > danny wrote:
>> >> >
>> >> >> Thanks Guido!
>> >> >>
>> >> >> The sample you gave performed slightly better: 304
>> >> seconds. I'll still
>> >> >> have to wait a few months to let it run on all my
>> >> data though. Yes, sitecode
>> >> >> has a "normal" index, not a spatial one.
>> >> Both the_geom's have spatial ones.
>> >> >> The european_waters polygon is one single complex
>> >> polygon, scale 1/100000.
>> >> >> What's annoying is that when I do an intersect
>> >> in Arc View it's almost
>> >> >> instantaneou s. Of course I don't get exaclty
>> >> what I wan't but with some
>> >> >> tweaking I can maybe get closer. That would be a
>> >> great disappointment
>> >> >> though. I was really looking forward to using
>> >> ST_Within intensively!
>> >> >>
>> >> >> Thanks Jean David, no performance gain though by
>> >> using it in the WHERE
>> >> >> clause, where I intended to use it in the first
>> >> place. I put it in the
>> >> >> Select for testing purposes.
>> >> >>
>> >> >> explain select sitecode from sites a,
>> >> europe_waters b
>> >> >> where st_intersects(a.the_geom,b.the_geom) and
>> >> a.sitecode = 'xxxx';
>> >> >>
>> >> >> "Nested Loop (cost=0.00..461.71 rows=1
>> >> width=10)"
>> >> >&g t; " Join Filter: _st_intersects(a.the_geom,
>> >> b.the_geom)"
>> >> >> " -> Seq Scan on sites a
>> >> (cost=0.00..453.43 rows=1 width=35764)"
>> >> >> " Filter: ((sitecode)::text =
>> >> 'xxxx'::text)"
>> >> >> " -> Index Scan using idx_europe_waters
>> >> on europe_waters b
>> >> >> (cost=0.00..8.27 rows=1 width=32)"
>> >> >> " Index Cond: (a.the_geom &&
>> >> b.the_geom)"
>> >> >> " Filter: (a.the_geom &&
>> >> b.the_geom)"
>> >> >>
>> >> >>
>> >> >>
>> >> >> On Tue, Aug 5, 2008 at 10:53 AM, Guido Lemoine
>> >> <guido.lemoine at jrc.it<mailto:
>> >> >> guido.lemoine at jrc.it>> wrote:
>> >> >>
>> >> >> Erik
>> >&g t; >>
>> >> >> Try this, and see if it is any faster:
>> >> >>
>> >> >> select st_intersects(a.the_geom,b.the_geom)
>> >> from sites a,
>> >> >> europe_waters b where a.the_geom &&
>> >> b.the_geom and
>> >> >> a.sitecode = 'xxxx';
>> >> >>
>> >> >> I guess you mean that sitecode is indexed (not
>> >> spatially indexed),
>> >> >> because it seems to be a varchar.
>> >> >> Both a and b should have a spatial index on
>> >> the_geom, but I reckon
>> >> >> that is already the case. Also,
>> >> >> b (europe_layer) should not be a single
>> >> polygon, but rather a set
>> >> >> of polygons.
>> >> >>
>> >> >> This is lesson 1 in the PostGIS tutor ial, more
>> >> or less. In return
>> >> >> for my 2 minutes effort, you are
>> >> >> obliged to report the new performance report,
>> >> so that future users
>> >> >> will benefit.
>> >> >> After all, I just saved you half a year...
>> >> >>
>> >> >> GL
>> >> >>
>> >> >>
>> >> >>
>> >> >> danny wrote:
>> >> >>
>> >> >> Hello,
>> >> >>
>> >> >> I'm wondering if it's normal to
>> >> have a 420 second response
>> >> >> time for the following query.
>> >> >> A spatial index has been set on the spatial
>> >> field (the_geom)
>> >> >> and other important fields (like sitecode).
>> >> I've g enerously
>> >> >> tweaked the memory options for postgresql.
>> >> >> With such a response time I would have to
>> >> let my query run for
>> >> >> half a year before getting the answer
>> >> I'm interested in! :)
>> >> >>
>> >> >> Anybody know how I can boost up the process
>> >> or is it doomed to
>> >> >> always be so slow? The europe_layer is
>> >> indeed a complex
>> >> >> polygon....
>> >> >>
>> >> >> select
>> >> st_intersects(a.the_geom,b.the_geom) from sites a,
>> >> >> europe_waters b where a.sitecode =
>> >> 'xxxx';
>> >> >>
>> >> >> "Nested Loop (cost=0.00..496.18
>> >> rows=1310 width=35786)"
>> >> >> " -> Seq Scan on sites a
>> >> (cost=0.00..453.43 rows=1
>> >> >> width=35754)"
>> >> >> " Filter: ((sitecode)::text =
>> >> 'xxxx'::text)"
>> >> >> " -> Seq Scan on europe_waters b
>> >> (cost=0.00..23.10
>> >> >> rows=1310 width=32)"
>> >> >>
>> >> >> My true objective would be to join thses
>> >> two tables through an
>> >> >> st_intersects but for the time it is
>> >> unconceivable.
>> >> >> Many thanks,
>> >> >>
>> >> >> Erik
>> >> >>
>> >> >>
>> >>
>> >> ------------------------------------------------------------------------
>> >> >>
>> >> >>
>> >> _________________________________________ ______
>> >> >> postgis-users mailing list
>> >> >> postgis-users at postgis.refractions.net
>> >> >>
>> >> <mailto:postgis-users at postgis.refractions.net>
>> >> >>
>> >> http://postgis.refractions.net/mailman/listinfo/postgis-users
>> >> >>
>> >> >>
>> >> >> --
>> >> ----------------------------------------------------------------
>> >> >> Guido Lemoine
>> >> >> Joint Research Centre, European Commission
>> >> >> Institute for the Protection and Security of
>> >> the Citizen (IPSC)
>> >> >> Support to External Security
>> >> >> Via E. Fermi, 2749 TP 267 Ispra 21027 (VA),
>> >> Italy
>> >> >> Tel. +39 0332 786239 (direct line) Fax. +39
>> >> 0332 785154
>> >> >> WWW: http://ses.jrc.it
>> >> >>
>> >> ----------------------------------------------------------------
>> >> >> Disclaimer:
>> >> >> Views expressed are those of the individual and
>> >> do not represent
>> >> >> the views of the European Commission
>> >> >>
>> >> >>
>> >> >> _______________________________________________
>> >> >> postgis-users mailing list
>> >> >> postgis-users at postgis.refractions.net
>> >> >>
>> >> <mailto:postgis-users at postgis.refractions.net>
>> >> >>
>> >> http://postgis.refractions.net/mailman/listinfo/postgis-users
>> >> >>
>> >> >>
>> >> >>
>> >> ------------------------------------------------------------------
>> >> ------
>> >> >>
>> >> >> _______________________________________________
>> >> >> postgis-users mailing list
>> >> >> postgis-users at postgis.refractions.net
>> >> >>
>> >> http://postgis.refractions.net/mailman/listinfo/postgis-users
>> >> >>
>> >> >>
>> >> >
>> >> > --
>> >> >
>> >> ----------------------------------------------------------------
>> >> > Guido Lemoine
>> >> > Joint Research Centre, European Commission
>> >> > Institute for the Protection and Security of the
>> >> Citizen (IPSC)
>> >> > Support to External Security
>> >> > Via E. Fermi, 2749 TP 267 Ispra 21027 (VA), Italy
>> >> > Tel. +39 0332 786239 (direct line) Fax. +39 0332
>> >> 785154
>> >> > WWW: http://s es.jrc.it
>> >> >
>> >> ----------------------------------------------------------------
>> >> > Disclaimer:
>> >> > Views expressed are those of the individual and do not
>> >> represent the views
>> >> > of the European Commission
>> >> >
>> >> >
>> >> > _______________________________________________
>> >> > postgis-users mailing list
>> >> > postgis-users at postgis.refractions.net
>> >> >
>> >> http://postgis.refractions.net/mailman/listinfo/postgis-users
>> >> >
>> >> _______________________________________________
>> >> postgis-users mailing list
>> >> postgis-users at postgis.refractions.net
>> >> http://postgis.refractions.net/mailman/listinfo/postgis-users
>> > _______________________________________________
>> > postgis-users m ailing list
>> > postgis-users at postgis.refractions.net
>> > http://postgis.refractions.net/mailman/listinfo/postgis-users
>> >
>>
>>
>> ------------------------------
>>
>> Message: 3
>> Date: Tue, 05 Aug 2008 20:45:01 +0100
>> From: Pedro Doria Meunier <pdoria at netmadeira.com>
>> Subject: Re: [postgis-users] Newbie Question
>> To: PostGIS Users Discussion <postgis-users at postgis.refractions.net>
>> Message-ID:
>> <1217965501.4389.8.camel at 213-190-214-90-cmts01.netmadeira.com>
>> Content-Type: text/plain; charset="us-ascii"
>>
>> Hi Adam,
>> This is what I use when I want to see which map to use for a vehicle's
>> current position:
>>
>> SELECT myfield1 FROM maps AS m
>> WHERE m.b_box && geomfromtext('POINT(-16.3 32.7)',4326) LIMIT 1
>>
>> (Note: there might be a newer/better function in PostGIS 1.3 ... )
>>
>> The top-level maps table consists of polygons defining countries
>> boundaries.
>> This query returns the correspondent map for the given location so this
>> should be applicable for your perimeters, which I guess are polygons :)
>>
>> HTH,
>> Pedro Doria Meunier
>>
>>
>> On Tue, 2008-08-05 at 14:24 -0500, Adam Vande More wrote:
>> > Hello,
>> >
>> > Perhaps my questions were addressed earlier or in some documentation I
>> > missed. If so, I apologize in advance as I am completely new to GIS and
>> > PostGIS. Judging from the documentation I have read, my needs are
>> > extremely simple vs PostGIS capabilities however I am unsure if the
>> > functionality of what I want actually exists.
>> >
>> > I deal with areas in the US, territories that are regional in nature but
>> > not bounded by common things eg state lines, water ways, countie s, etc.
>> > They are irregular shapes of varying land area. I also have sets of
>> > latitude and longitude points. I would like some type of system where I
>> > enter in a perimeter of territories, and find out whether a particular
>> > latitude and longitude point lies within that perimeter.
>> >
>> > My understanding of this is that I would need to either manually create
>> > a GML type file with the perimeters in it, or use some type of map
>> > authoring software to accomplish this goal. Are my understandings
>> > correct and if so is there a recommendation for a particular piece of
>> > software which can fulfill this?
>> >
>> > Also, are the libraries/functions available for determining whether a
>> > point is in a perimeter?
>> >
>> > Much appreciated,
>> >
>> -------------- next part --------------
>> A non-tex t attachment was scrubbed...
>> Name: not available
>> Type: application/pgp-signature
>> Size: 189 bytes
>> Desc: This is a digitally signed message part
>> Url :
>> http://lists.refractions.net/pipermail/postgis-users/attachments/20080805/356f2fed/attachment-0001.bin
>>
>> ------------------------------
>>
>> Message: 4
>> Date: Tue, 5 Aug 2008 12:46:14 -0700
>> From: "Paul Ramsey" <pramsey at cleverelephant.ca>
>> Subject: Re: [postgis-users] Newbie Question
>> To: "PostGIS Users Discussion" <postgis-users at postgis.refractions.net>
>> Message-ID:
>> <30fe546d0808051246p6386bee5u2877020e9392a36d at mail.gmail.com>
>> Content-Type: text/plain; charset=ISO-8859-1
>>
>> Well, that's the opposite of the usual point-in-polygon use case
>> (given an existing set of polygons and a point input, tell me what
>> polygons it falls within), but it's equally tractable.
> &gt ;
>> You need an input interface, and that will be the hardest part of your
>> task. You can do something in Google Maps, if you like, then stuff the
>> geometry down into a script (php, asp, whatever) that in turn
>> inscribes it into a SQL statement ala:
>>
>> select points.id, points.name from points where
>> st_contains('POLYGON(....)', points.geometry);
>>
>> Where the 'POLYGON(...)' is constructed by your script from your user
>> interface input.
>>
>> Enjoy,
>>
>> P.
>>
>> On Tue, Aug 5, 2008 at 12:24 PM, Adam Vande More <adam at imedmobility.com>
>> wrote:
>> > Hello,
>> >
>> > Perhaps my questions were addressed earlier or in some documentation I
>> > missed. If so, I apologize in advance as I am completely new to GIS and
>> > PostGIS. Judging from the documentation I have read, my needs are
>> > extremely
>> > simple vs PostGIS capabilit ies however I am unsure if the functionality
>> > of
>> > what I want actually exists.
>> >
>> > I deal with areas in the US, territories that are regional in nature but
>> > not
>> > bounded by common things eg state lines, water ways, counties, etc. They
>> > are irregular shapes of varying land area. I also have sets of latitude
>> > and
>> > longitude points. I would like some type of system where I enter in a
>> > perimeter of territories, and find out whether a particular latitude and
>> > longitude point lies within that perimeter.
>> > My understanding of this is that I would need to either manually create
>> > a
>> > GML type file with the perimeters in it, or use some type of map
>> > authoring
>> > software to accomplish this goal. Are my understandings correct and if
>> > so
>> > is there a recommendation for a particular piece of software which can
>> > fulfill this?
>> >
> & gt; > Also, are the libraries/functions available for determining whether
> a point
>> > is in a perimeter?
>> >
>> > Much appreciated,
>> >
>> > --
>> > Adam Vandemore
>> > Systems Administrator
>> > IMED Mobility
>> > (605) 498-1610
>> >
>> > _______________________________________________
>> > postgis-users mailing list
>> > postgis-users at postgis.refractions.net
>> > http://postgis.refractions.net/mailman/listinfo/postgis-users
>> >
>>
>>
>> ------------------------------
>>
>> Message: 5
>> Date: Tue, 05 Aug 2008 20:57:39 +0100
>> From: Pedro Doria Meunier <pdoria at netmadeira.com>
>> Subject: Re: [postgis-users] Newbie Question
>> To: PostGIS Users Discussion <postgis-users at postgis.refractions.net>
>> Message-ID:
>> <1217966259.4389.12.camel at 213-190-214-90-cmts01.netmadeira.com&gt ;
>> Content-Type: text/plain; charset="us-ascii"
>>
>> Shame on me... I should have read Adam's question more attentively ...
>>
>> Adam:
>> Paul is absolutely right.
>>
>> Regards,
>> Pedro Doria Meunier
>>
>>
>> On Tue, 2008-08-05 at 12:46 -0700, Paul Ramsey wrote:
>> > Well, that's the opposite of the usual point-in-polygon use case
>> > (given an existing set of polygons and a point input, tell me what
>> > polygons it falls within), but it's equally tractable.
>> >
>> > You need an input interface, and that will be the hardest part of your
>> > task. You can do something in Google Maps, if you like, then stuff the
>> > geometry down into a script (php, asp, whatever) that in turn
>> > inscribes it into a SQL statement ala:
>> >
>> > select points.id, points.name from points where
>> > st_contains('POLYGON(....)', poi nts.geometry);
>> >
>> > Where the 'POLYGON(...)' is constructed by your script from your user
>> > interface input.
>> >
>> > Enjoy,
>> >
>> > P.
>> >
>> > On Tue, Aug 5, 2008 at 12:24 PM, Adam Vande More <adam at imedmobility.com>
>> > wrote:
>> > > Hello,
>> > >
>> > > Perhaps my questions were addressed earlier or in some documentation I
>> > > missed. If so, I apologize in advance as I am completely new to GIS
>> > > and
>> > > PostGIS. Judging from the documentation I have read, my needs are
>> > > extremely
>> > > simple vs PostGIS capabilities however I am unsure if the
>> > > functionality of
>> > > what I want actually exists.
>> > >
>> > > I deal with areas in the US, territories that are regional in nature
>> > > but not
>> > > bounded by common things eg state lines, water ways, counties, etc.
>> > > They
>> > > are irregular shapes of varying land area. I also have sets of
>> > > latitude and
>> > > longitude points. I would like some type of system where I enter in a
>> > > perimeter of territories, and find out whether a particular latitude
>> > > and
>> > > longitude point lies within that perimeter.
>> > > My understanding of this is that I would need to either manually
>> > > create a
>> > > GML type file with the perimeters in it, or use some type of map
>> > > authoring
>> > > software to accomplish this goal. Are my understandings correct and if
>> > > so
>> > > is there a recommendation for a particular piece of software which can
>> > > fulfill this?
>> > >
>> > > Also, are the libraries/functions available for determining whether a
>> > > point
>> > > is in a perimeter?
>> > >
>> > > Much appreciated,
>> > >
>> > > --
>> > > Adam Vandemore
>> > > Systems Administrator
>> > > IMED Mobility
>> > > (605) 498-1610
>> > >
>> > > _______________________________________________
>> > > postgis-users mailing list
>> > > postgis-users at postgis.refractions.net
>> > > http://postgis.refractions.net/mailman/listinfo/postgis-users
>> > >
>> > _______________________________________________
>> > postgis-users mailing list
>> > postgis-users at postgis.refractions.net
>> > http://postgis.refractions.net/mailman/listinfo/postgis-users
>> -------------- next part --------------
>> A non-text attachment was scrubbed...
>> Name: not available
>> Type: application/pgp-signature
>> Size: 189 bytes
>> Desc: This is a digitally signed message part
>> Url :
>> http://lists.refractions.net/pipermail/postgis-users/attachments/20080805/98ad00f1/attachment-0001.bin
>>
>> ------------------------------
>>
>> Message: 6
>> Date: Tue, 05 Aug 2008 15:00:35 -0700
>> From: Donald Ijams <dsijams at gmail.com>
>> Subject: [postgis-users] Update Error
>> To: postgis-users at postgis.refractions.net
>> Message-ID: <4898CD83.6070704 at gmail.com>
>> Content-Type: text/plain; charset="iso-8859-1"
>>
>> When I attempt to update a geometry column from a field in another table
>> with:
>>
>> UPDATE recent
>> SET xypoint = ST_GeomFromText((raw.x, raw.y),102249)
>> from raw
>> where recent.casenum=raw.casenum;
>>
>> I get:
>>
>> *SQL error:*
>>
>> ERROR: function st_geomfromtext(record, integer) does not exist
>> LINE 3: xypoint = ST_GeomFromText((raw.x, raw.y),102249)
>> ^
>> HINT: No function matches the given name and argument types. You might
>> need to add explicit type casts.
>>
>>
>> Probably something simple, but I'm not seeing it.
>>
>> Don
>>
>> -------------- next part --------------
>> An HTML attachment was scrubbed...
>> URL:
>> http://lists.refractions.net/pipermail/postgis-users/attachments/20080805/7f4711de/attachment-0001.html
>>
>> ------------------------------
>>
>> Message: 7
>> Date: Tue, 5 Aug 2008 22:59:56 -0400
>> From: "Obe, Regina" <robe.dnd at cityofboston.gov>
>> Subject: RE: [postgis-users] Update Error
>> To: "PostGIS Users Discussion" <postgis-users at postgis.refractions.net>
>> Message-ID:
>> <53F9CF533E1AA14EA1F8C5C08ABC08D20197A16F at ZDND.DND.boston.cob>
>> Content-Type: text/plain; charset="iso-8859-1"
>>
>>
>> ST_GeomFromText takes the form ST_GeomFromText('POINT(' || raw.x || ' ' ||
>> raw.y || ')', 102249)
>>
>> Anyrate for point geometries its more efficient to use ST_MakePoin t so
>> something like
>>
>>
>> UPDATE recent
>> SET xypoint = ST_SetSRID(ST_MakePoint(raw.x, raw.y),102249)
>> from raw
>> where recent.casenum=raw.casenum;
>>
>> Hope that helps,
>> Regina
>>
>>
>> -----Original Message-----
>> From: postgis-users-bounces at postgis.refractions.net on behalf of Donald
>> Ijams
>> Sent: Tue 8/5/2008 6:00 PM
>> To: postgis-users at postgis.refractions.net
>> Subject: [postgis-users] Update Error
>>
>> When I attempt to update a geometry column from a field in another table
>> with:
>>
>> UPDATE recent
>> SET xypoint = ST_GeomFromText((raw.x, raw.y),102249)
>> from raw
>> where recent.casenum=raw.casenum;
>>
>> I get:
>>
>> *SQL error:*
>>
>> ERROR: function st_geomfromtext(record, integer) does not exist
>> LINE 3: xypoint = ST_GeomFromText((raw.x, raw.y),102249)
> &gt ; ^
>> HINT: No function matches the given name and argument types. You might
>> need to add explicit type casts.
>>
>>
>> Probably something simple, but I'm not seeing it.
>>
>> Don
>>
>>
>>
>>
>>
>>
>> -----------------------------------------
>> The substance of this message, including any attachments, may be
>> confidential, legally privileged and/or exempt from disclosure
>> pursuant to Massachusetts law. It is intended
>> solely for the addressee. If you received this in error, please
>> contact the sender and delete the material from any computer.
>> -------------- next part --------------
>> An HTML attachment was scrubbed...
>> URL:
>> http://lists.refractions.net/pipermail/postgis-users/attachments/20080805/61d5ef9c/attachment-0001.html
>>
>> ------------------------------
>>
>> Message: 8
>> Date: Wed, 6 Aug 2008 00:17:05 -0700 (PDT)
>> From: "kreshna_iceheart at yahoo.com" <kreshna_iceheart at yahoo.com>
>> Subject: [postgis-users] Still cannot display postgis layers using
>> ms4w
>> To: postgis-users at postgis.refractions.net
>> Message-ID: <295797.42389.qm at web56101.mail.re3.yahoo.com>
>> Content-Type: text/plain; charset="us-ascii"
>>
>> Help. Somebody please help. Anybody. :-(
>>
>> I am using postgreSQL 8.2.5 for Windows with postGIS. The database was
>> installed using Windows .msi installation package. I am also using ms4w
>> 2.2.7 that comes with MapServer 5.0.2. All installed on Windows XP Service
>> Pack 2.
>>
>> The MapServer does support postGIS. When I executed mapserv -v, I got the
>> following:
>> MapServer version 5.0.2 OUTPUT=GIF OUTPUT=PNG OUTPUT=JPEG OUTPUT=WBMP
>> OUTPUT=PDF OUTPUT=SWF OUTPUT=SVG SUPPORTS=PROJ SUPPORTS=AGG
>> SUPPORTS=FREETYPE SUPPORTS=WMS_SERVER SUPPORTS=WMS_CLIENT
>> SUPPORTS=WFS_SERVER SUPPOR TS=WFS_CLIENT SUPPORTS=WCS_SERVER
>> SUPPORTS=SOS_SERVER SUPPORTS=FASTCGI SUPPORTS=THREADS SUPPORTS=GEOS
>> INPUT=JPEG INPUT=POSTGIS INPUT=OGR INPUT=GDAL INPUT=SHAPEFILE
>>
>> In the postgreSQL database, I have created a database named test_mapserv.
>> I also created a schema named rencana, and a table named kawasanbencana. The
>> name of the geometry column is the_geom.
>> I have created a spatial index (kawasanbencanaspix) on the geometry
>> column. I also use lowercase to write the 'from' statement.
>>
>> The problem is: MapServer always gives the "..Query error. Error executing
>> POSTGIS DECLARE.." error message, no matter what I do. And I have tried
>> everything ad nauseam.
>>
>> I have tried specifying the DATA line without USING UNIQUE ... USING SRID
>> line. I have tried specifying the DATA line with USING UNIQUE ...USING SRID.
>> I have tried using double apostrophe (") and single apostrophe ('). It
>> doesn't matter. Everything I have tried has failed, a nd I still get the
>> error message no matter what I did.
>>
>> My mapfiles are attached in this mail. Each reflect my futile attempt, and
>> all of them generated the error message.
>>
>>
>> My first attempt:
>> =================
>>
>> This is my first attempt. I have both the mapfile and the error message
>> attached on this mail.
>> CONNECTIONTYPE postgis
>> CONNECTION 'host=localhost user=rtrw password=123456 port=5432
>> dbname=test_mapserver'
>> DATA 'the_geom from rencana.kawasanbencana'
>>
>> And I got the following error message:
>> msDrawMap(): Image handling error. Failed to draw layer named
>> 'kawasanbencana'. prepare_database(): Query error. Error executing POSTGIS
>> DECLARE (the actual query) statement: 'DECLARE mycursor BINARY CURSOR FOR
>> SELECT asbinary(force_collection(force_2d(the_geom)),'NDR'),gid::text from
>> rencana.kawasanbencana WHERE () and (the_geom && setSRID( 'BOX3D(528635
>> .525378392 9074209,609092.162121608
>> 9134520)'::BOX3D,find_srid('','rencana.kawasanbencana','the_geom') ))'
>> Postgresql reports the error as 'ERROR: syntax error at or near ")" LINE 1:
>> ...DR'),gid::text from rencana.kawasanbencana WHERE () and (the... ^ ' More
>> Help: Error with POSTGIS data variable. You specified 'check your .map
>> file'. Standard ways of specifiying are : (1) 'geometry_column from
>> geometry_table' (2) 'geometry_column from (sub query) as foo using unique
>> column name using SRID=srid#' Make sure you put in the 'using unique column
>> name' and 'using SRID=#' clauses in. For more help, please see
>> http://postgis.refractions.net/documentation/ Mappostgis.c - version of
>> Jan 23/2004.
>>
>>
>> My second attempt:
>> ==================
>>
>> I tried adding USING UNIQUE to my DATA line. There is something odd with
>> my spatial table though: whenever I view the data contents with pgAdmin, I
>> can always see the oid column (first column to the left, before the PK
>> column). However, I cannot find the oid column when viewing the table's
>> column definition. The question is: does my table have oid column or not?
>> And does it have something to do with the error?
>>
>> Nonethless, here's my lines:
>> CONNECTIONTYPE postgis
>> CONNECTION 'host=localhost user=rtrw password=123456 port=5432
>> dbname=test_mapserver'
>> DATA 'the_geom from rencana.kawasanbencana using unique oid'
>>
>>
>> And I still get the same error message as above. The difference is that
>> the error message now says 'oid' instead of 'gid'. For instance,
>> "(the_geom)),'NDR'),gid" becomes "(the_geom)),'NDR'),oid", but everything
>> else remains the same.
>>
>>
>> My third attempt:
>> =================
>>
>> Like I said in my second attempt, I'm not sure whether my table has oid
>> column or not. Thus, in my third attempt I used the gid column instead,
>> especially since the gid column is my primary key column anyway. This is my
>> lines:
>> CONNECTIONTYPE postgis
>> CONNECTION 'host=localhost user=rtrw password=123456 port=5432
>> dbname=test_mapserver'
>> DATA 'the_geom from rencana.kawasanbencana using unique gid'
>>
>> And I got exactly the same error message as my first attempt.
>>
>>
>> My fourth attempt:
>> ==================
>>
>> In my fourth attempt, I used the gid column for my unique id, and I am
>> using SRID=1. Here is the lines:
>> CONNECTIONTYPE postgis
>> CONNECTION 'host=localhost user=rtrw password=123456 port=5432
>> dbname=test_mapserver'
>> DATA 'the_geom from rencana.kawasanbencana using unique gid using SRID=1'
>>
>> Well I still got the error message. This time, the error message is
>> slightly different. Instead of getting "::BOX3D,find_srid", I got
>> "::BOX3D,1".
>>
>> Here is the complete error message:
>> msDrawMap (): Image handling error. Failed to draw layer named
>> 'kawasanbencana'. prepare_database(): Query error. Error executing POSTGIS
>> DECLARE (the actual query) statement: 'DECLARE mycursor BINARY CURSOR FOR
>> SELECT asbinary(force_collection(force_2d(the_geom)),'NDR'),gid::text from
>> rencana.kawasanbencana WHERE () and (the_geom && setSRID(
>> 'BOX3D(528635.525378392 9074209,609092.162121608 9134520)'::BOX3D,1) )'
>> Postgresql reports the error as 'ERROR: syntax error at or near ")" LINE 1:
>> ...DR'),gid::text from rencana.kawasanbencana WHERE () and (the... ^ ' More
>> Help: Error with POSTGIS data variable. You specified 'check your .map
>> file'. Standard ways of specifiying are : (1) 'geometry_column from
>> geometry_table' (2) 'geometry_column from (sub query) as foo using unique
>> column name using SRID=srid#' Make sure you put in the 'using unique column
>> name' and 'using SRID=#' clauses in. For more help, please see
>> http://postgis.refractions.net/documentation/
>> Mappostgis.c - version of Jan 23/2004.
>>
>>
>> My fifth attempt:
>> =================
>>
>> Now I had been pretty much desperate, I just wrote using SRID=#, closed my
>> eyes, and hoped for the best. Here's the lines:
>> CONNECTIONTYPE postgis
>> CONNECTION 'host=localhost user=rtrw password=123456 port=5432
>> dbname=test_mapserver'
>> DATA 'the_geom from rencana.kawasanbencana using unique gid using SRID=#'
>>
>> I still get an error message. This time, the error message is different
>> altogether. Here's the message:
>> msDrawMap(): Image handling error. Failed to draw layer named
>> 'kawasanbencana'. msPOSTGISLayerParseData(): Query error. Error parsing
>> POSTGIS data variable: You specified 'using SRID=#' but didnt have any
>> numbers! More Help: Error with POSTGIS data variable. You specified
>> 'the_geom from rencana.kawasanbencana using unique gid using SRID=#'.
>> Standard ways of specifiying are : (1) 'geometry_column from ge
>> ometry_table' (2) 'geometry_column from (sub query) as foo using unique
>> column name using SRID=srid#' Make sure you put in the 'using unique column
>> name' and 'using SRID=#' clauses in. For more help, please see
>> http://postgis.refractions.net/documentation/ Mappostgis.c - version of Jan
>> 23/2004.
>>
>>
>> I'm pretty much desperate. Could somebody please tell what's wrong on my
>> mapfile? Or maybe ms4w just cannot display postGIS data? Do I need to
>> install postGIS and MapServer on Linux, in order to get them work together?
>> Please help. :-(
>>
>>
>> Thanks,
>> -Kresh
>>
>>
>>
>>
>>
>>
>>
>>
>>
>>
>>
>> -------------- next part --------------
>> A non-text attachment was scrubbed...
>> Name: test_mapserv_01.map
>> Type: application/octet-stream
>> Size: 1005 bytes
>> Desc: not available
>> Url : http://lists.refractions.net/pipermail/pos
>> tgis-users/attachments/20080806/ff428038/test_mapserv_01.dll
>> -------------- next part --------------
>> msDrawMap(): Image handling error. Failed to draw layer named
>> 'kawasanbencana'. prepare_database(): Query error. Error executing POSTGIS
>> DECLARE (the actual query) statement: 'DECLARE mycursor BINARY CURSOR FOR
>> SELECT asbinary(force_collection(force_2d(the_geom)),'NDR'),gid::text from
>> rencana.kawasanbencana WHERE () and (the_geom && setSRID(
>> 'BOX3D(528635.525378392 9074209,609092.162121608
>> 9134520)'::BOX3D,find_srid('','rencana.kawasanbencana','the_geom') ))'
>> Postgresql reports the error as 'ERROR: syntax error at or near ")" LINE 1:
>> ...DR'),gid::text from rencana.kawasanbencana WHERE () and (the... ^ ' More
>> Help: Error with POSTGIS data variable. You specified 'check your .map
>> file'. Standard ways of specifiying are : (1) 'geometry_column from
>> geometry_table' (2) 'geometry_column from (sub query) as foo using unique
>> column name using SRID=srid#' Make sure you put in the 'using unique column
>> name' and 'using SRID=#' clauses in. For more help, please see
>> http://postgis.refract
>> ions.net/documentation/ Mappostgis.c - version of Jan 23/2004.
>> -------------- next part --------------
>> A non-text attachment was scrubbed...
>> Name: test_mapserv_02.map
>> Type: application/octet-stream
>> Size: 1022 bytes
>> Desc: not available
>> Url :
>> http://lists.refractions.net/pipermail/postgis-users/attachments/20080806/ff428038/test_mapserv_02.bin
>> -------------- next part --------------
>> msDrawMap(): Image handling error. Failed to draw layer named
>> 'kawasanbencana'. prepare_database(): Query error. Error executing POSTGIS
>> DECLARE (the actual query) statement: 'DECLARE mycursor BINARY CURSOR FOR
>> SELECT asbinary(force_collection(force_2d(the_geom)),'NDR'),oid::text from
>> rencana.kawasanbencana WHERE () and (the_geom && setSRID(
>> 'BOX3D(528635.525378392 9074209,609092.16212 1608
>> 9134520)'::BOX3D,find_srid('','rencana.kawasanbencana','the_geom') ))'
>> Postgresql reports the error as 'ERROR: syntax error at or near ")" LINE 1:
>> ...DR'),oid::text from rencana.kawasanbencana WHERE () and (the... ^ ' More
>> Help: Error with POSTGIS data variable. You specified 'check your .map
>> file'. Standard ways of specifiying are : (1) 'geometry_column from
>> geometry_table' (2) 'geometry_column from (sub query) as foo using unique
>> column name using SRID=srid#' Make sure you put in the 'using unique column
>> name' and 'using SRID=#' clauses in. For more help, please see
>> http://postgis.refract
>> ions.net/documentation/ Mappostgis.c - version of Jan 23/2004.
>> -------------- next part --------------
>> A non-text attachment was scrubbed...
>> Name: test_mapserv_03.map
>> Type: application/octet-stream
>> Size: 1022 bytes
>> Desc: not available
>> Url :
>> http://lists.refractions.net/pipermail/postgis-users/attachments/20080806/ff4280
>> 38/test_mapserv_03.dll
>> -------------- next part --------------
>> msDrawMap(): Image handling error. Failed to draw layer named
>> 'kawasanbencana'. prepare_database(): Query error. Error executing POSTGIS
>> DECLARE (the actual query) statement: 'DECLARE mycursor BINARY CURSOR FOR
>> SELECT asbinary(force_collection(force_2d(the_geom)),'NDR'),gid::text from
>> rencana.kawasanbencana WHERE () and (the_geom && setSRID(
>> 'BOX3D(528635.525378392 9074209,609092.162121608
>> 9134520)'::BOX3D,find_srid('','rencana.kawasanbencana','the_geom') ))'
>> Postgresql reports the error as 'ERROR: syntax error at or near ")" LINE 1:
>> ...DR'),gid::text from rencana.kawasanbencana WHERE () and (the... ^ ' More
>> Help: Error with POSTGIS data variable. You specified 'check your .map
>> file'. Standard ways of specifiying are : (1) 'geometry_column from
>> geometry_table' (2) 'geometry_column from (sub query) as foo using unique
>> column name using SRID=srid#' Make sure you put in the 'using unique colu mn
>> name' and 'using SRID=#' clauses in. For more help, please see
>> http://postgis.refract
>> ions.net/documentation/ Mappostgis.c - version of Jan 23/2004.
>> -------------- next part --------------
>> A non-text attachment was scrubbed...
>> Name: test_mapserv_04.map
>> Type: application/octet-stream
>> Size: 1035 bytes
>> Desc: not available
>> Url :
>> http://lists.refractions.net/pipermail/postgis-users/attachments/20080806/ff428038/test_mapserv_04.bin
>> -------------- next part --------------
>> msDrawMap(): Image handling error. Failed to draw layer named
>> 'kawasanbencana'. prepare_database(): Query error. Error executing POSTGIS
>> DECLARE (the actual query) statement: 'DECLARE mycursor BINARY CURSOR FOR
>> SELECT asbinary(force_collection(force_2d(the_geom)),'NDR'),gid::text from
>> rencana.kawasanbencana WHERE () and (the_geom && setSRID(
>> 'BOX3D(528635.525378392 9074209,609092.162121608 9134520)'::BOX3D,1) )'
>> Postgresql reports the error as 'ERROR: syntax error at or near ")" LINE 1:
>> ...DR'),gid::text from rencana.kawasanbencana WHERE () and (the... ^ ' More
>> Help: Error with POSTGIS data variable. You specified 'check your .map
>> file'. Standard ways of specifiying are : (1) 'geometry_column from
>> geometry_table' (2) 'geometry_column from (sub query) as foo using unique
>> column name using SRID=srid#' Make sure you put in the 'using unique column
>> name' and 'using SRID=#' clauses in. For more help, please see
>> http://postgis.refractions.net/documentation/ Mappostgis.c - version o
>> f Jan 23/2004.
>> -------------- next part --------------
>> A non-text attachment was scrubbed...
>> Name: test_mapserv_05.map
>> Type: application/octet-stream
>> Size: 1035 bytes
>> Desc: not available
>> Url :
>> http://lists.refractions.net/pipermail/postgis-users/attachments/20080806/ff428038/test_mapserv_05.dll
>> -------------- next part --------------
>> msDraw Map(): Image handling error. Failed to draw layer named
>> 'kawasanbencana'. msPOSTGISLayerParseData(): Query error. Error parsing
>> POSTGIS data variable: You specified 'using SRID=#' but didnt have any
>> numbers! More Help: Error with POSTGIS data variable. You specified
>> 'the_geom from rencana.kawasanbencana using unique gid using SRID=#'.
>> Standard ways of specifiying are : (1) 'geometry_column from geometry_table'
>> (2) 'geometry_column from (sub query) as foo using unique column name using
>> SRID=srid#' Make sure you put in the 'using unique column name' and 'using
>> SRID=#' clauses in. For more help, please see
>> http://postgis.refractions.net/documentation/ Mappostgis.c - version of Jan
>> 23/2004.
>>
>> ------------------------------
>>
>> Message: 9
>> Date: Wed, 6 Aug 2008 04:10:32 -0400
>> From: "Paragon Corporation" <lr at pcorp.us>
>> Subject: RE: [postgis-users] Still cannot display postgis layers using
>> ms4w
>> To: "'PostGIS Users Discussion'"
>> <postgis-users at postgis.refractions.net>
>> Message-ID: <001501c8f79b$e5fd1860$4d812e40 at l>
>> Content-Type: text/plain; charset="us-ascii"
>>
>> Just a guess. Get rid of your FILTER clause or actually put a where
>> condition in your FILTER clause. I think the fact that your FILTER clause
>> is empty is causing PostgreSQL to generate an invalid query that looks
>> something like this
>>
>> SELECT *
>> FROM sometable
>> WHERE
>>
>> Can't have where without a where condition.
>>
>> Hope that helps,
>> Regina
>>
>>
>>
>> -----Original Message-----
>> From: postgis-users-bounces at postgis.refractions.net
>> [mailto:postgis-users-bounces at postgis.refractions.net] On Behalf Of
>> kreshna_iceheart at yahoo.com
>> Sent: Wednesday, August 06, 2008 3:17 AM
>> To: postgis-users at postgis.refractions.net
>> Subject: [postgis-users] Still cannot display postgis layers using ms4w
>>
>> Help. Somebody please help. Anybody. :-(
>>
>> I am using postgreSQL 8.2.5 for Windows with postGIS. The database was
>> installed using Windows .msi installation package. I am also using ms4w
>> 2.2.7 that comes with MapServer 5.0.2. All installed on Windows XP Service
>> Pack 2.
>>
>> The MapServer does support postGIS. When I executed mapserv -v, I got the
>> following:
>> MapServer version 5.0.2 OUTPUT=GIF OUTPUT=PNG OUTPUT=JPEG OUTPUT=WBMP
>> OUTPUT=PDF OUTPUT=SWF OUTPUT=SVG SUPPORTS=PROJ SUPPORTS=AGG
>> SUPPORTS=FREETYPE SUPPORTS=WMS_SERVER SUPPORTS=WMS_CLIENT
>> SUPPORTS=WFS_SERVER SUPPORTS=WFS_CLIENT SUPPORTS=WCS_SERVER
>> SUPPORTS=SOS_SERVER SUPPORTS=FASTCGI SUPPORTS=THREADS SUPPORTS=GEOS
>> INPUT=JPEG INPUT=POSTGIS INPUT=OGR INPUT=GDAL INPUT=SHAPEFILE
>>
>> In the postgreSQL database, I have created a dat abase named test_mapserv.
>> I
>> also created a schema named rencana, and a table named kawasanbencana. The
>> name of the geometry column is the_geom.
>> I have created a spatial index (kawasanbencanaspix) on the geometry
>> column.
>> I also use lowercase to write the 'from' statement.
>>
>> The problem is: MapServer always gives the "..Query error. Error executing
>> POSTGIS DECLARE.." error message, no matter what I do. And I have tried
>> everything ad nauseam.
>>
>> I have tried specifying the DATA line without USING UNIQUE ... USING SRID
>> line. I have tried specifying the DATA line with USING UNIQUE ...USING
>> SRID.
>> I have tried using double apostrophe (") and single apostrophe ('). It
>> doesn't matter. Everything I have tried has failed, and I still get the
>> error message no matter what I did.
>>
>> My mapfiles are attached in this mail. Each reflect my futile attempt, and
>> a ll of them generated the error message.
>>
>>
>> My first attempt:
>> =================
>>
>> This is my first attempt. I have both the mapfile and the error message
>> attached on this mail.
>> CONNECTIONTYPE postgis
>> CONNECTION 'host=localhost user=rtrw password=123456 port=5432
>> dbname=test_mapserver'
>> DATA 'the_geom from rencana.kawasanbencana'
>>
>> And I got the following error message:
>> msDrawMap(): Image handling error. Failed to draw layer named
>> 'kawasanbencana'. prepare_database(): Query error. Error executing POSTGIS
>> DECLARE (the actual query) statement: 'DECLARE mycursor BINARY CURSOR FOR
>> SELECT asbinary(force_collection(force_2d(the_geom)),'NDR'),gid::text from
>> rencana.kawasanbencana WHERE () and (the_geom && setSRID(
>> 'BOX3D(528635.525378392 9074209,609092.162121608
>> 9134520)'::BOX3D,find_srid('','rencana.kawasanben cana','the_geom') ))'
>> Postgresql reports the error as 'ERROR: syntax error at or near ")" LINE
>> 1:
>> ...DR'),gid::text from rencana.kawasanbencana WHERE () and (the... ^ '
>> More
>> Help: Error with POSTGIS data variable. You specified 'check your .map
>> file'. Standard ways of specifiying are : (1) 'geometry_column from
>> geometry_table' (2) 'geometry_column from (sub query) as foo using unique
>> column name using SRID=srid#' Make sure you put in the 'using unique
>> column
>> name' and 'using SRID=#' clauses in. For more help, please see
>> http://postgis.refractions.net/documentation/ Mappostgis.c - version of
>> Jan
>> 23/2004.
>>
>>
>> My second attempt:
>> ==================
>>
>> I tried adding USING UNIQUE to my DATA line. There is something odd with
>> my
>> spatial table though: whenever I view the data contents with pgAdmin, I
>> can
>> always see the oid column (first column to the l eft, before the PK
>> column).
>> However, I cannot find the oid column when viewing the table's column
>> definition. The question is: does my table have oid column or not? And
>> does
>> it have something to do with the error?
>>
>> Nonethless, here's my lines:
>> CONNECTIONTYPE postgis
>> CONNECTION 'host=localhost user=rtrw password=123456 port=5432
>> dbname=test_mapserver'
>> DATA 'the_geom from rencana.kawasanbencana using unique oid'
>>
>>
>> And I still get the same error message as above. The difference is that
>> the
>> error message now says 'oid' instead of 'gid'. For instance,
>> "(the_geom)),'NDR'),gid" becomes "(the_geom)),'NDR'),oid", but everything
>> else remains the same.
>>
>>
>> My third attempt:
>> =================
>>
>> Like I said in my second attempt, I'm not sure whether my table has oid
>> column or not. Thus, in my third attemp t I used the gid column instead,
>> especially since the gid column is my primary key column anyway. This is
>> my
>> lines:
>> CONNECTIONTYPE postgis
>> CONNECTION 'host=localhost user=rtrw password=123456 port=5432
>> dbname=test_mapserver'
>> DATA 'the_geom from rencana.kawasanbencana using unique gid'
>>
>> And I got exactly the same error message as my first attempt.
>>
>>
>> My fourth attempt:
>> ==================
>>
>> In my fourth attempt, I used the gid column for my unique id, and I am
>> using
>> SRID=1. Here is the lines:
>> CONNECTIONTYPE postgis
>> CONNECTION 'host=localhost user=rtrw password=123456 port=5432
>> dbname=test_mapserver'
>> DATA 'the_geom from rencana.kawasanbencana using unique gid using
>> SRID=1'
>>
>> Well I still got the error message. This time, the error message is
>> slightly
>> different. Instead of gett ing "::BOX3D,find_srid", I got "::BOX3D,1".
>>
>> Here is the complete error message:
>> msDrawMap(): Image handling error. Failed to draw layer named
>> 'kawasanbencana'. prepare_database(): Query error. Error executing POSTGIS
>> DECLARE (the actual query) statement: 'DECLARE mycursor BINARY CURSOR FOR
>> SELECT asbinary(force_collection(force_2d(the_geom)),'NDR'),gid::text from
>> rencana.kawasanbencana WHERE () and (the_geom && setSRID(
>> 'BOX3D(528635.525378392 9074209,609092.162121608 9134520)'::BOX3D,1) )'
>> Postgresql reports the error as 'ERROR: syntax error at or near ")" LINE
>> 1:
>> ...DR'),gid::text from rencana.kawasanbencana WHERE () and (the... ^ '
>> More
>> Help: Error with POSTGIS data variable. You specified 'check your .map
>> file'. Standard ways of specifiying are : (1) 'geometry_column from
>> geometry_table' (2) 'geometry_column from (sub query) as foo using unique
>> column name using SRID=srid#' Make sure you put in the 'using unique
>> column
>> name' and 'using SRID=#' clauses in. For more help, please see
>> http://postgis.refractions.net/documentation/
>> Mappostgis.c - version of Jan 23/2004.
>>
>>
>> My fifth attempt:
>> =================
>>
>> Now I had been pretty much desperate, I just wrote using SRID=#, closed my
>> eyes, and hoped for the best. Here's the lines:
>> CONNECTIONTYPE postgis
>> CONNECTION 'host=localhost user=rtrw password=123456 port=5432
>> dbname=test_mapserver'
>> DATA 'the_geom from rencana.kawasanbencana using unique gid using
>> SRID=#'
>>
>> I still get an error message. This time, the error message is different
>> altogether. Here's the message:
>> msDrawMap(): Image handling error. Failed to draw layer named
>> 'kawasanbencana'. msPOSTGISLayerParseData(): Query error. Error parsing
>> POSTGI S data variable: You specified 'using SRID=#' but didnt have any
>> numbers! More Help: Error with POSTGIS data variable. You specified
>> 'the_geom from rencana.kawasanbencana using unique gid using SRID=#'.
>> Standard ways of specifiying are : (1) 'geometry_column from
>> geometry_table'
>> (2) 'geometry_column from (sub query) as foo using unique column name
>> using
>> SRID=srid#' Make sure you put in the 'using unique column name' and 'using
>> SRID=#' clauses in. For more help, please see
>> http://postgis.refractions.net/documentation/ Mappostgis.c - version of
>> Jan
>> 23/2004.
>>
>>
>> I'm pretty much desperate. Could somebody please tell what's wrong on my
>> mapfile? Or maybe ms4w just cannot display postGIS data? Do I need to
>> install postGIS and MapServer on Linux, in order to get them work
>> together?
>> Please help. :-(
>>
>>
>> Thanks,
>> -Kresh
>>
>>
>> < br>>
>>
>>
>>
>>
>>
>>
>>
>>
>>
>>
>>
>> ------------------------------
>>
>> Message: 10
>> Date: Wed, 6 Aug 2008 01:33:38 -0700 (PDT)
>> From: "kreshna_iceheart at yahoo.com" <kreshna_iceheart at yahoo.com>
>> Subject: RE: [postgis-users] Still cannot display postgis layers using
>> ms4w
>> To: PostGIS Users Discussion <postgis-users at postgis.refractions.net>
>> Message-ID: <801791.19065.qm at web56103.mail.re3.yahoo.com>
>> Content-Type: text/plain; charset=us-ascii
>>
>> Dear Regina....
>>
>> YOU ARE CORRECT!!! I got rid of the empty FILTER clause and now the map is
>> displayed properly! :-)
>>
>> Many many many (times infinity) thanks!
>> -Kresh
>>
>>
>> --- On Wed, 8/6/08, Paragon Corporation <lr at pcorp.us> wrote:
>>
>> > From: Paragon Corporation <lr at pcorp.us&g t;
>> > Subject: RE: [postgis-users] Still cannot display postgis layers using
>> > ms4w
>> > To: "'PostGIS Users Discussion'" <postgis-users at postgis.refractions.net>
>> > Date: Wednesday, August 6, 2008, 4:10 AM
>> > Just a guess. Get rid of your FILTER clause or actually put
>> > a where
>> > condition in your FILTER clause. I think the fact that
>> > your FILTER clause
>> > is empty is causing PostgreSQL to generate an invalid query
>> > that looks
>> > something like this
>> >
>> > SELECT *
>> > FROM sometable
>> > WHERE
>> >
>> > Can't have where without a where condition.
>> >
>> > Hope that helps,
>> > Regina
>> >
>> >
>> >
>> > -----Original Message-----
>> > From: postgis-users-bounces at postgis.refractions.net
>> > [mailto:postgis-users-bounces at postgis.refrac tions.net] On
>> > Behalf Of
>> > kreshna_iceheart at yahoo.com
>> > Sent: Wednesday, August 06, 2008 3:17 AM
>> > To: postgis-users at postgis.refractions.net
>> > Subject: [postgis-users] Still cannot display postgis
>> > layers using ms4w
>> >
>> > Help. Somebody please help. Anybody. :-(
>> >
>> > I am using postgreSQL 8.2.5 for Windows with postGIS. The
>> > database was
>> > installed using Windows .msi installation package. I am
>> > also using ms4w
>> > 2.2.7 that comes with MapServer 5.0.2. All installed on
>> > Windows XP Service
>> > Pack 2.
>> >
>> > The MapServer does support postGIS. When I executed mapserv
>> > -v, I got the
>> > following:
>> > MapServer version 5.0.2 OUTPUT=GIF OUTPUT=PNG OUTPUT=JPEG
>> > OUTPUT=WBMP
>> > OUTPUT=PDF OUTPUT=SWF OUTPUT=SVG SUPPORTS=PROJ SUPPORTS=AGG
>> > SUPPORTS=FREETYPE SUPPORTS=WMS_SERVER SUPPORTS=WMS_CLIENT
>> > SUPPORTS=WFS_SERVER SUPPORTS=WFS_CLIENT SUPPORTS=WCS_SERVER
>> > SUPPORTS=SOS_SERVER SUPPORTS=FASTCGI SUPPORTS=THREADS
>> > SUPPORTS=GEOS
>> > INPUT=JPEG INPUT=POSTGIS INPUT=OGR INPUT=GDAL
>> > INPUT=SHAPEFILE
>> >
>> > In the postgreSQL database, I have created a database named
>> > test_mapserv. I
>> > also created a schema named rencana, and a table named
>> > kawasanbencana. The
>> > name of the geometry column is the_geom.
>> > I have created a spatial index (kawasanbencanaspix) on the
>> > geometry column.
>> > I also use lowercase to write the 'from' statement.
>> >
>> > The problem is: MapServer always gives the "..Query
>> > error. Error executing
>> > POSTGIS DECLARE.." error message, no matter what I do.
>> > And I have tried
>> &g t; everything ad nauseam.
>> >
>> > I have tried specifying the DATA line without USING UNIQUE
>> > ... USING SRID
>> > line. I have tried specifying the DATA line with USING
>> > UNIQUE ...USING SRID.
>> > I have tried using double apostrophe (") and single
>> > apostrophe ('). It
>> > doesn't matter. Everything I have tried has failed, and
>> > I still get the
>> > error message no matter what I did.
>> >
>> > My mapfiles are attached in this mail. Each reflect my
>> > futile attempt, and
>> > all of them generated the error message.
>> >
>> >
>> > My first attempt:
>> > =================
>> >
>> > This is my first attempt. I have both the mapfile and the
>> > error message
>> > attached on this mail.
>> > CONNECTIONTYPE postgis
>> > CONNECTION 'host=localhost user=rtrw< br>> > password=123456 port=5432
>> > dbname=test_mapserver'
>> > DATA 'the_geom from rencana.kawasanbencana'
>> >
>> > And I got the following error message:
>> > msDrawMap(): Image handling error. Failed to draw layer
>> > named
>> > 'kawasanbencana'. prepare_database(): Query error.
>> > Error executing POSTGIS
>> > DECLARE (the actual query) statement: 'DECLARE mycursor
>> > BINARY CURSOR FOR
>> > SELECT
>> > asbinary(force_collection(force_2d(the_geom)),'NDR'),gid::text
>> > from
>> > rencana.kawasanbencana WHERE () and (the_geom &&
>> > setSRID(
>> > 'BOX3D(528635.525378392 9074209,609092.162121608
>> > 9134520)'::BOX3D,find_srid('','rencana.kawasanbencana','the_geom')
>> > ))'
>> > Postgresql reports the error as 'ERROR: syntax error at
>> > or near ")" LINE 1:
>> > ...DR'),gid::text from rencana.kawasanbencana WHERE ()
>> > and (the... ^ ' More
>> > Help: Error with POSTGIS data variable. You specified
>> > 'check your .map
>> > file'. Standard ways of specifiying are : (1)
>> > 'geometry_column from
>> > geometry_table' (2) 'geometry_column from (sub
>> > query) as foo using unique
>> > column name using SRID=srid#' Make sure you put in the
>> > 'using unique column
>> > name' and 'using SRID=#' clauses in. For more
>> > help, please see
>> > http://postgis.refractions.net/documentation/ Mappostgis.c
>> > - version of Jan
>> > 23/2004.
>> >
>> >
>> > My second attempt:
>> > ==================
>> >
>> > I tried adding USING UNIQUE to my DATA line. There is
>> > something odd with my
>> > spatial table though: whenever I view the data contents
>> > with pgAdmin, I can> > always see the oid column (first column to the
>> > left, before
>> > the PK column).
>> > However, I cannot find the oid column when viewing the
>> > table's column
>> > definition. The question is: does my table have oid column
>> > or not? And does
>> > it have something to do with the error?
>> >
>> > Nonethless, here's my lines:
>> > CONNECTIONTYPE postgis
>> > CONNECTION 'host=localhost user=rtrw
>> > password=123456 port=5432
>> > dbname=test_mapserver'
>> > DATA 'the_geom from rencana.kawasanbencana using
>> > unique oid'
>> >
>> >
>> > And I still get the same error message as above. The
>> > difference is that the
>> > error message now says 'oid' instead of
>> > 'gid'. For instance,
>> > "(the_geom)),'NDR'),gid" becomes
>> > "(the_geom)),'NDR'),oid", but everything
>> > else remains the same.
>> >
>> >
>> > My third attempt:
>> > =================
>> >
>> > Like I said in my second attempt, I'm not sure whether
>> > my table has oid
>> > column or not. Thus, in my third attempt I used the gid
>> > column instead,
>> > especially since the gid column is my primary key column
>> > anyway. This is my
>> > lines:
>> > CONNECTIONTYPE postgis
>> > CONNECTION 'host=localhost user=rtrw
>> > password=123456 port=5432
>> > dbname=test_mapserver'
>> > DATA 'the_geom from rencana.kawasanbencana using
>> > unique gid'
>> >
>> > And I got exactly the same error message as my first
>> > attempt.
>> >
>> >
>> > My fourth attempt:
>> > ==================
>> >
>> > In my fourth attempt, I used the gid column for my un ique
>> > id, and I am using
>> > SRID=1. Here is the lines:
>> > CONNECTIONTYPE postgis
>> > CONNECTION 'host=localhost user=rtrw
>> > password=123456 port=5432
>> > dbname=test_mapserver'
>> > DATA 'the_geom from rencana.kawasanbencana using
>> > unique gid using
>> > SRID=1'
>> >
>> > Well I still got the error message. This time, the error
>> > message is slightly
>> > different. Instead of getting
>> > "::BOX3D,find_srid", I got "::BOX3D,1".
>> >
>> > Here is the complete error message:
>> > msDrawMap(): Image handling error. Failed to draw layer
>> > named
>> > 'kawasanbencana'. prepare_database(): Query error.
>> > Error executing POSTGIS
>> > DECLARE (the actual query) statement: 'DECLARE mycursor
>> > BINARY CURSOR FOR
>> > SELECT
>> > asbinary(force_collection(fo rce_2d(the_geom)),'NDR'),gid::text
>> > from
>> > rencana.kawasanbencana WHERE () and (the_geom &&
>> > setSRID(
>> > 'BOX3D(528635.525378392 9074209,609092.162121608
>> > 9134520)'::BOX3D,1) )'
>> > Postgresql reports the error as 'ERROR: syntax error at
>> > or near ")" LINE 1:
>> > ...DR'),gid::text from rencana.kawasanbencana WHERE ()
>> > and (the... ^ ' More
>> > Help: Error with POSTGIS data variable. You specified
>> > 'check your .map
>> > file'. Standard ways of specifiying are : (1)
>> > 'geometry_column from
>> > geometry_table' (2) 'geometry_column from (sub
>> > query) as foo using unique
>> > column name using SRID=srid#' Make sure you put in the
>> > 'using unique column
>> > name' and 'using SRID=#' clauses in. For more
>> > help, please see
>> > http://postgis.refractions.net/documentation/> > Mappostgis.c - version
>> > of Jan 23/2004.
>> >
>> >
>> > My fifth attempt:
>> > =================
>> >
>> > Now I had been pretty much desperate, I just wrote using
>> > SRID=#, closed my
>> > eyes, and hoped for the best. Here's the lines:
>> > CONNECTIONTYPE postgis
>> > CONNECTION 'host=localhost user=rtrw
>> > password=123456 port=5432
>> > dbname=test_mapserver'
>> > DATA 'the_geom from rencana.kawasanbencana using
>> > unique gid using
>> > SRID=#'
>> >
>> > I still get an error message. This time, the error message
>> > is different
>> > altogether. Here's the message:
>> > msDrawMap(): Image handling error. Failed to draw layer
>> > named
>> > 'kawasanbencana'. msPOSTGISLayerParseData(): Query
>> > error. Error parsing
>> > POSTGIS data variable: You s pecified 'using SRID=#'
>> > but didnt have any
>> > numbers! More Help: Error with POSTGIS data variable. You
>> > specified
>> > 'the_geom from rencana.kawasanbencana using unique gid
>> > using SRID=#'.
>> > Standard ways of specifiying are : (1) 'geometry_column
>> > from geometry_table'
>> > (2) 'geometry_column from (sub query) as foo using
>> > unique column name using
>> > SRID=srid#' Make sure you put in the 'using unique
>> > column name' and 'using
>> > SRID=#' clauses in. For more help, please see
>> > http://postgis.refractions.net/documentation/ Mappostgis.c
>> > - version of Jan
>> > 23/2004.
>> >
>> >
>> > I'm pretty much desperate. Could somebody please tell
>> > what's wrong on my
>> > mapfile? Or maybe ms4w just cannot display postGIS data? Do
>> > I need to
>> > install postGIS a nd MapServer on Linux, in order to get
>> > them work together?
>> > Please help. :-(
>> >
>> >
>> > Thanks,
>> > -Kresh
>> >
>> >
>> >
>> >
>> >
>> >
>> >
>> >
>> >
>> >
>> >
>> >
>> >
>> > _______________________________________________
>> > postgis-users mailing list
>> > postgis-users at postgis.refractions.net
>> > http://postgis.refractions.net/mailman/listinfo/postgis-users
>>
>>
>>
>>
>>
>> ------------------------------
>>
>> Message: 11
>> Date: Wed, 6 Aug 2008 16:04:32 +0200
>> From: Khavindra Sivenandan <khavindra at sefsa.co.za>
>> Subject: [postgis-users] Configuration and getting started issues
>> To: "postgis-users at postgis.refractions.net"
>> <postgis-users at postgis.refractions.n et>
>> Message-ID:
>> <1677E63198782240825AEF9E248718DE764FB4613E at sefsaserver.sefsa.local>
>> Content-Type: text/plain; charset="us-ascii"
>>
>> Hi There
>>
>> I am new to PostGIS/PostgreSQL. Can somebody help me with the following
>> support information.
>>
>> I went as far as installing PostgreSQL & PostGIS. I require to configure
>> the Geometry_Columns Table as in the manual. It is somewhat vague as to the
>> inputs into the columns F_Table_Catalogue, F_Table_Schema & F_Table_Name. I
>> have managed to sort out the other columns as indicated in the manual
>> (4.2.2)
>>
>> I have managed to connect to the database from various opensource desktop
>> GIS software, but my database is not configured correctly. Is there any
>> detailed material, other than the online manual I can look up or can
>> somebody please help me to move on further from here. Ultimately, I would
>> like to move data (shapefiles) into the database within some classification
>> system (e.g. by Province, state, etc...) which can then be utilised by
>> desktop and web users.
>>
>>
>> Thanks
>>
>> Khavin Sivenandan
>> Khavindra at sefsa.co.za
>> -------------- next part --------------
>> An HTML attachment was scrubbed...
>> URL:
>> http://lists.refractions.net/pipermail/postgis-users/attachments/20080806/34cc54ab/attachment-0001.html
>>
>> ------------------------------
>>
>> Message: 12
>> Date: Wed, 6 Aug 2008 16:06:14 +0200
>> From: danny <whatevar89 at gmail.com>
>> Subject: [postgis-users] area calculation problem
>> To: "PostGIS Users Discussion" <postgis-users at postgis.refractions.net>
>> Message-ID:
>> <5c218180808060706k3bdd3a0ar5db0e09c6da759cb at mail.gmail.com>
>> Content-Type: text/plain; charset="iso-8859-1"
>>
>> Hi, I have a problem.
>> I'm trying to calculate the area of a great number of si tes on a specific
>> territory.
>> The problem is that these sites can overlap (as they do very frequently)
>> and
>> I do not want to count the same portion of land twice.
>> A st_union would work with this right? Unfortunately my experience with
>> this
>> function is that it is way too slow for the numerous polygons I would have
>> to merge. As a matter fact it would take years.
>> Does anybody know of a hack around this? (In PostGIS or ArcGIS)
>>
>> Many thanks!
>> -------------- next part --------------
>> An HTML attachment was scrubbed...
>> URL:
>> http://lists.refractions.net/pipermail/postgis-users/attachments/20080806/3ce94bd4/attachment-0001.html
>>
>> ------------------------------
>>
>> Message: 13
>> Date: Wed, 6 Aug 2008 11:19:19 -0400
>> From: "Burgholzer,Robert" <rwburgholzer at deq.virginia.gov>
>> Subject: [postgis-users] filing the holes in the swiss cheese?
> &g t; To: "PostGIS Users Discussion" <postgis-users at postgis.refractions.net>
>> Message-ID:
>> <6C097DA58429B743A67070F98BE73A370379B67B at deqex01.deq.local>
>> Content-Type: text/plain; charset="us-ascii"
>>
>> I am trying to do a query of cities that are contained by counties (in
>> Virginia, US), so that I can have a list that cross-references by FIPS
>> these relationships.
>>
>>
>>
>> My query looks like this:
>>
>>
>>
>> select b.poli1 as fips, b.name, a.poli1 as contained_by_fips, a.name
>>
>> from poli_bounds as a , poli_bounds as b
>>
>> where a.poli1 <> b.poli1 and a.the_geom && b.the_geom
>>
>> and contains(convexhull(a.the_geom), b.the_geom)
>>
>> and a.projectid = 1
>>
>> and b.projectid = 1;
>>
>>
>>
>>
>>
>> This works fairly well, but omits a handful of cities that are most
>> definitely within the boundaries of the county.
>>
>>
>>
>> Is "convexhull()" the correct function for this? Are my geometries
>> goofed up?
>>
>>
>>
>> Thanks!
>>
>>
>>
>> Robert W. Burgholzer
>>
>> Surface Water Modeler
>>
>> Office of Water Supply and Planning
>>
>> Virginia Department of Environmental Quality
>>
>> rwburgholzer at deq.virginia.gov
>>
>> 804-698-4405
>>
>> Open Source Modeling Tools:
>>
>> http://sourceforge.net/projects/npsource/
>>
>>
>>
>> -------------- next part --------------
>> An HTML attachment was scrubbed...
>> URL:
>> http://lists.refractions.net/pipermail/postgis-users/attachments/20080806/1d85ac47/attachment-0001.html
>>
>> ------------------------------
>>
>> Message: 14
>> Date: Wed, 6 Aug 2008 11:15:43 -0700> From: "Paul Ramsey"
>> <pramsey at cleverelephant.ca>
>> Subject: Re: [postgis-users] filing the holes in the swiss cheese?
>> To: "PostGIS Users Discussion" <postgis-users at postgis.refractions.net>
>> Message-ID:
>> <30fe546d0808061115y273ba5b1xc178f9e3cd9f9c05 at mail.gmail.com>
>> Content-Type: text/plain; charset=ISO-8859-1
>>
>> try
>>
>> select b.poli1 as fips, b.name, a.poli1 as contained_by_fips, a.name
>> from poli_bounds as a , poli_bounds as b
>> where a.poli1 <> b.poli1
>> and st_contains(a.the_geom, st_pointonsurface(b.the_geom))
>> and a.projectid = 1
>> and b.projectid = 1;
>>
>> that should rid you of the boundary conditions plaguing st_contains.
>>
>> P.
>>
>> ps - note the "modern" st_contains(), with implicit index call.
>>
>> On Wed, Aug 6, 2008 at 8:19 AM, Burgholzer,Robert
>> <rwburgholzer at deq.virginia.gov> wrote:
>> > I am trying to do a query of cities that are contained by counties (in
>> > Virginia, US), so that I can have a list that cross-references by FIPS
>> > these
>> > relationships.
>> >
>> >
>> >
>> > My query looks like this:
>> >
>> >
>> >
>> > select b.poli1 as fips, b.name, a.poli1 as contained_by_fips, a.name
>> >
>> > from poli_bounds as a , poli_bounds as b
>> >
>> > where a.poli1 <> b.poli1 and a.the_geom && b.the_geom
>> >
>> > and contains(convexhull(a.the_geom), b.the_geom)
>> >
>> > and a.projectid = 1
>> >
>> > and b.projectid = 1;
>> >
>> >
>> >
>> >
>> >
>> > This works fairly well, but omits a handful of cities that are most
>> > definitely within the boundaries of the county.
>> >
>> >
>> >
> &g t; > Is "convexhull()" the correct function for this? Are my geometries
> goofed
>> > up?
>> >
>> >
>> >
>> > Thanks!
>> >
>> >
>> >
>> > Robert W. Burgholzer
>> >
>> > Surface Water Modeler
>> >
>> > Office of Water Supply and Planning
>> >
>> > Virginia Department of Environmental Quality
>> >
>> > rwburgholzer at deq.virginia.gov
>> >
>> > 804-698-4405
>> >
>> > Open Source Modeling Tools:
>> >
>> > http://sourceforge.net/projects/npsource/
>> >
>> >
>> >
>> > _______________________________________________
>> > postgis-users mailing list
>> > postgis-users at postgis.refractions.net
>> > http://postgis.refractions.net/mailman/listinfo/postgis-users
>> >
>> >
>>
>>
>> ------------------------------
>> < br>> _______________________________________________
>> postgis-users mailing list
>> postgis-users at postgis.refractions.net
>> http://postgis.refractions.net/mailman/listinfo/postgis-users
>>
>>
>> End of postgis-users Digest, Vol 70, Issue 6
>> ********************************************
>
> ________________________________
> Connect to the next generation of MSN Messenger  Get it now!
> _______________________________________________
> postgis-users mailing list
> postgis-users at postgis.refractions.net
> http://postgis.refractions.net/mailman/listinfo/postgis-users
>
>



More information about the postgis-users mailing list