No subject
Tue Oct 7 10:08:30 PDT 2008
columns
as per the argument group_=3D'GONDWANA'
But the resultant table coal is not having a 'primary key'.
table geology has gid as primary key.
Pl give me a way how I can get a primary key on table coal.
Table geology has (copy pasted from pgadmin)
CREATE TABLE public.geology
(
gid integer NOT NULL DEFAULT nextval('geology_gid_seq'::regclass),
objectid bigint,
group_ character varying(35),
geometry geometry,
CONSTRAINT geology_pkey PRIMARY KEY (gid),
CONSTRAINT enforce_dims_geometry CHECK (ndims(geometry) =3D 2),
CONSTRAINT enforce_geotype_geometry CHECK (geometrytype(geometry) =3D =
'MULTIPOLYGON'::text OR geometry IS NULL),
CONSTRAINT enforce_srid_geometry CHECK (srid(geometry) =3D (-1))
)
WITH (OIDS=3DFALSE);
ALTER TABLE public.geology OWNER TO postgres;
The query results table coal without primary key (copy pasted from =
pgadmin)
CREATE TABLE public.coal
(
gid integer,
objectid bigint,
group_ character varying(35),
geometry geometry
)
WITH (OIDS=3DFALSE);
ALTER TABLE public.coal OWNER TO postgres;
Looking for local information? Find it on Yahoo! Local =
http://in.local.yahoo.com/
------------------------------
Message: 9
Date: Mon, 14 Sep 2009 18:12:38 -0500
From: P Kishor <punk.kish at gmail.com>
Subject: Re: [postgis-users] No Primary Key
To: PostGIS Users Discussion <postgis-users at postgis.refractions.net>
Message-ID:
<cdf6db500909141612p5962b097nff813ff0955d5bb4 at mail.gmail.com>
Content-Type: text/plain; charset=3DISO-8859-1
On Mon, Sep 14, 2009 at 5:34 PM, Ravi <ravivundavalli at yahoo.com> wrote:
> Hi,
> have a problem with a query resulting in a table without a primary =
key.
> Pl suggest a proper SQL
> Thanks in anticipation
> Ravi
>
> create table coal AS SELECT * from geology where group_=3D'GONDWANA';
>
> From this SQL a table coal is created from table geology with all its =
columns
> as per the argument group_=3D'GONDWANA'
>
> But the resultant table coal is not having a 'primary key'.
> table geology has gid as primary key.
>
> Pl give me a way how I can get a primary key on table coal.
>
>
> Table geology has (copy pasted from pgadmin)
> CREATE TABLE public.geology
> (
> ?gid integer NOT NULL DEFAULT nextval('geology_gid_seq'::regclass),
> ?objectid bigint,
> ?group_ character varying(35),
> ?geometry geometry,
> ?CONSTRAINT geology_pkey PRIMARY KEY (gid),
> ?CONSTRAINT enforce_dims_geometry CHECK (ndims(geometry) =3D 2),
> ?CONSTRAINT enforce_geotype_geometry CHECK (geometrytype(geometry) =3D =
'MULTIPOLYGON'::text OR geometry IS NULL),
> ?CONSTRAINT enforce_srid_geometry CHECK (srid(geometry) =3D (-1))
> )
> WITH (OIDS=3DFALSE);
> ALTER TABLE public.geology OWNER TO postgres;
>
> The query results table coal without primary key (copy pasted from =
pgadmin)
>
> CREATE TABLE public.coal
> (
> ?gid integer,
> ?objectid bigint,
> ?group_ character varying(35),
> ?geometry geometry
> )
> WITH (OIDS=3DFALSE);
> ALTER TABLE public.coal OWNER TO postgres;
>
>
>
Create table "coal" with the correct primary key definition (make the
table defintion the same as the "geology" table) and then
INSERT INTO coal SELECT FROM geology WHERE group_=3D'GONDWANA';
Step 1:
> ? ? ?Looking for local information? Find it on Yahoo! Local =
http://in.local.yahoo.com/
> _______________________________________________
> postgis-users mailing list
> postgis-users at postgis.refractions.net
> http://postgis.refractions.net/mailman/listinfo/postgis-users
>
--=20
-----------------------------------------------------------------------
Assertions are politics; backing up assertions with evidence is science
=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=
=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=
=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D=3D
Sent from Madison, WI, United States
------------------------------
Message: 10
Date: Mon, 14 Sep 2009 18:25:49 -0700 (PDT)
From: pcreso at pcreso.com
Subject: Re: [postgis-users] Center of Points Collection
To: PostGIS Users Discussion <postgis-users at postgis.refractions.net>
Message-ID: <536453.95341.qm at web33208.mail.mud.yahoo.com>
Content-Type: text/plain; charset=3Dus-ascii
Hi guys,
A bit more difficult, & way out in left field, but if you use PL/R to =
create a median function for Postgres, you could build your point from =
the median(X) & (Y) values instead of the average.
Where this would actually lie obviously depends on the distribution of =
the points. The centroid is most affected by (actually defined by) =
outlying points, the avg somewhat less & the median less still.
Of course once you have PL/R to play with, you have much more =
flexibility to look at returning statistics from datasets than just the =
median.
Cheers,
Brent Wood
--- On Tue, 9/15/09, Kevin Neufeld <kneufeld at refractions.net> wrote:
> From: Kevin Neufeld <kneufeld at refractions.net>
> Subject: Re: [postgis-users] Center of Points Collection
> To: "PostGIS Users Discussion" <postgis-users at postgis.refractions.net>
> Date: Tuesday, September 15, 2009, 8:51 AM
> Paul Ramsey wrote:
> > Faster than creating a multipoint is to recognize that
> ST_Centroid()
> > is just going to return the center of the bbox of the
> collection
> > anyways...
>=20
> Unfortunately, Paul, ST_Centroid returns the center of
> mass, not the center of the bbox.
>=20
> SELECT astext(st_centroid(st_collect(column1))),
> FROM (values ('POINT(0 0)'),
> =20
> ('POINT(0 1)'),
> =20
> ('POINT(0 2)'),
> =20
> ('POINT(1 0)')) as foo;
> astext
> ------------------
> POINT(0.25 0.75)
> (1 row)
>=20
> Your second post, taking the avg of the x,y does seem to be
> the nice approach, and produces the same results as
> ST_Centroid - the center of mass.
>=20
> SELECT astext(st_makepoint(avg(st_x(column1)),
> avg(st_y(column1))))
> FROM (values ('POINT(0 0)'),
> =20
> ('POINT(0 1)'),
> =20
> ('POINT(0 2)'),
> =20
> ('POINT(1 0)')) as foo;
> astext
> ------------------
> POINT(0.25 0.75)
> (1 row)
>=20
> If Dustin is after the center of the collection, then
> something along your first suggestion might be more
> appropriate.
> (taking the center of the extents)
>=20
> Cheers,
> Kevin
> _______________________________________________
> postgis-users mailing list
> postgis-users at postgis.refractions.net
> http://postgis.refractions.net/mailman/listinfo/postgis-users
>=20
------------------------------
Message: 11
Date: Mon, 14 Sep 2009 21:25:13 -0600
From: Richard Greenwood <richard.greenwood at gmail.com>
Subject: Re: [postgis-users] multipoint2point
To: PostGIS Users Discussion <postgis-users at postgis.refractions.net>
Message-ID:
<ae9185aa0909142025l179d5822t888f5fd38bea4f79 at mail.gmail.com>
Content-Type: text/plain; charset=3DISO-8859-1
You could use st_multi() to cast them al to multi point:
ST_Multi(geometry)
Returns the geometry as a MULTI* geometry. If the geometry is
already a MULTI*, it is returned unchanged.
SELECT ST_MULTI(wkb_geometry) FROM table;
On Mon, Sep 14, 2009 at 6:22 AM, <jj.wag at gmx.de> wrote:
> Hi,
> I have a postgistable with different geometrytyps in one table: points =
and
> multipoints (get this information via ST_GeometryType).
> UMN Mapserver have no problems to show the geometrys, but if I try to =
export
> this table via ogr2ogr there is an error:
> Attempt to write non-poin(MULTIPOINT) to point shapefile.
> How can I cast all my geometrys to point?
> Jo
> _______________________________________________
> postgis-users mailing list
> postgis-users at postgis.refractions.net
> http://postgis.refractions.net/mailman/listinfo/postgis-users
>
--=20
Richard Greenwood
richard.greenwood at gmail.com
www.greenwoodmap.com
------------------------------
Message: 12
Date: Tue, 15 Sep 2009 11:52:56 +0530
From: poonam jalwan <poonam.jalwan at orkash.com>
Subject: [postgis-users] Pgrouting assign_vertex_id function finding
lat lon values of vertexes problem
To: postgis-users at postgis.refractions.net
Message-ID: <4AAF32C0.6060809 at orkash.com>
Content-Type: text/plain; charset=3D"iso-8859-1"; Format=3D"flowed"
Hi All,
I have roads data i run the following function for creating source and=20
target and cost values.
SELECT assign_vertex_id(table_name, snapping_range,=20
geometry_column_name, edge_id_column_name)
now i want lat lon values of each and every start and end(vrtexes) with=20
respect to roads data.Can anyone help me to resolve this problem.
Thanks,
Poonam
-------------- next part --------------
An HTML attachment was scrubbed...
URL: =
<http://postgis.refractions.net/pipermail/postgis-users/attachments/20090=
915/d647efea/attachment-0001.html>
------------------------------
Message: 13
Date: Tue, 15 Sep 2009 08:43:10 +0200
From: D?ster Horst <Horst.Duester at bd.so.ch>
Subject: [postgis-users] st_intersection error
To: postgis-users <postgis-users at postgis.refractions.net>
Message-ID: <H00002eb058418d3.1252996988.srsofaioi6145.ktso.ch at MHS>
Content-Type: text/plain; charset=3D"iso-8859-1"
I do have 2 simple linestrings which partly have the same geometry (see
attached dump and image). My aim is to detect the common part of these
two linestrings with the following query:
select st_intersection(a.the_geom, b.the_geom), a.myid
from aline1, aline2
But as the result I get a geometrycollection with one POINT and one
LINESTING. I expected one linestring. Is it a bug??
I'll appreciate any hints.
My system:
Postgis 1.4.0
Geos 3.1.1
Postgres 8.3.6
Best regards
=20
Dr. Horst D?ster
Stv. Amtschef / GIS-Koordinator=20
Kanton Solothurn
Bau- und Justizdepartement
Amt f?r Geoinformation
SO!GIS Koordination
R?tistrasse 4
CH-4501 Solothurn
Telefon ++41(0)32 627 25 32
Telefax ++41(0)32 627 22 14
mailto:horst.duester at bd.so.ch
http://www.agi.so.ch
-------------- next part --------------
An HTML attachment was scrubbed...
URL: =
<http://postgis.refractions.net/pipermail/postgis-users/attachments/20090=
915/d3b3b03b/attachment-0001.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: aline1.sql
Type: application/octet-stream
Size: 1199 bytes
Desc: not available
URL: =
<http://postgis.refractions.net/pipermail/postgis-users/attachments/20090=
915/d3b3b03b/attachment-0003.obj>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: aline2.sql
Type: application/octet-stream
Size: 1268 bytes
Desc: not available
URL: =
<http://postgis.refractions.net/pipermail/postgis-users/attachments/20090=
915/d3b3b03b/attachment-0004.obj>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: intersecterror.png
Type: application/octet-stream
Size: 13093 bytes
Desc: not available
URL: =
<http://postgis.refractions.net/pipermail/postgis-users/attachments/20090=
915/d3b3b03b/attachment-0005.obj>
------------------------------
Message: 14
Date: Tue, 15 Sep 2009 10:03:58 +0200
From: <jj.wag at gmx.de>
Subject: Re: [postgis-users] multipoint2point
To: "PostGIS Users Discussion" <postgis-users at postgis.refractions.net>
Message-ID: <89B0DFA22A7B49EEA95F836FAF78EAB4 at SVEN>
Content-Type: text/plain; format=3Dflowed; charset=3D"iso-8859-1";
reply-type=3Doriginal
thank you very much ST_MULTI works fine and solves my problem, is there =
a=20
function to cast a MULTIPOINT to POINT?
----- Original Message -----=20
From: "Richard Greenwood" <richard.greenwood at gmail.com>
To: "PostGIS Users Discussion" <postgis-users at postgis.refractions.net>
Sent: Tuesday, September 15, 2009 5:25 AM
Subject: Re: [postgis-users] multipoint2point
> You could use st_multi() to cast them al to multi point:
> ST_Multi(geometry)
>
> Returns the geometry as a MULTI* geometry. If the geometry is
> already a MULTI*, it is returned unchanged.
>
> SELECT ST_MULTI(wkb_geometry) FROM table;
>
>
> On Mon, Sep 14, 2009 at 6:22 AM, <jj.wag at gmx.de> wrote:
>> Hi,
>> I have a postgistable with different geometrytyps in one table: =
points=20
>> and
>> multipoints (get this information via ST_GeometryType).
>> UMN Mapserver have no problems to show the geometrys, but if I try to =
>> export
>> this table via ogr2ogr there is an error:
>> Attempt to write non-poin(MULTIPOINT) to point shapefile.
>> How can I cast all my geometrys to point?
>> Jo
>> _______________________________________________
>> postgis-users mailing list
>> postgis-users at postgis.refractions.net
>> http://postgis.refractions.net/mailman/listinfo/postgis-users
>>
>
>
>
> --=20
> Richard Greenwood
> richard.greenwood at gmail.com
> www.greenwoodmap.com
> _______________________________________________
> postgis-users mailing list
> postgis-users at postgis.refractions.net
> http://postgis.refractions.net/mailman/listinfo/postgis-users
>=20
------------------------------
Message: 15
Date: Tue, 15 Sep 2009 10:07:27 +0200
From: nicklas.aven at jordogskog.no
Subject: Re: [postgis-users] multipoint2point
To: PostGIS Users Discussion <postgis-users at postgis.refractions.net>
Message-ID: <200909150807.n8F87Ri6011714 at mail-core.space2u.com>
Content-Type: text/plain; charset=3D"iso-8859-1"
You can use st_dump for that pupose. Then you will get one row per point =
in your multipoint. select (st_dump(the_multipoint)).geom from thetable =
/Nicklas
2009-09-15 wrote:
thank you very much ST_MULTI works fine and solves my problem, is there =
a=20
>function to cast a MULTIPOINT to POINT?
>
>
>----- Original Message -----=20
>From: "Richard Greenwood"
>To: "PostGIS Users Discussion"
>Sent: Tuesday, September 15, 2009 5:25 AM
>Subject: Re: [postgis-users] multipoint2point
>
>
>> You could use st_multi() to cast them al to multi point:
>> ST_Multi(geometry)
>>
>> Returns the geometry as a MULTI* geometry. If the geometry is
>> already a MULTI*, it is returned unchanged.
>>
>> SELECT ST_MULTI(wkb_geometry) FROM table;
>>
>>
>> On Mon, Sep 14, 2009 at 6:22 AM,wrote:
>>> Hi,
>>> I have a postgistable with different geometrytyps in one table: =
points=20
>>> and
>>> multipoints (get this information via ST_GeometryType).
>>> UMN Mapserver have no problems to show the geometrys, but if I try =
to=20
>>> export
>>> this table via ogr2ogr there is an error:
>>> Attempt to write non-poin(MULTIPOINT) to point shapefile.
>>> How can I cast all my geometrys to point?
>>> Jo
>>> _______________________________________________
>>> postgis-users mailing list
>>> postgis-users at postgis.refractions.net
>>> http://postgis.refractions.net/mailman/listinfo/postgis-users
>>>
>>
>>
>>
>> --=20
>> Richard Greenwood
>> richard.greenwood at gmail.com
>> www.greenwoodmap.com
>> _______________________________________________
>> postgis-users mailing list
>> postgis-users at postgis.refractions.net
>> http://postgis.refractions.net/mailman/listinfo/postgis-users
>>=20
>
>_______________________________________________
>postgis-users mailing list
>postgis-users at postgis.refractions.net
>http://postgis.refractions.net/mailman/listinfo/postgis-users
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: =
<http://postgis.refractions.net/pipermail/postgis-users/attachments/20090=
915/097326a9/attachment-0001.html>
------------------------------
Message: 16
Date: Tue, 15 Sep 2009 10:27:01 +0200
From: strk <strk at keybit.net>
Subject: Re: [postgis-users] st_intersection error
To: PostGIS Users Discussion <postgis-users at postgis.refractions.net>
Message-ID: <20090915082701.GB30456 at keybit.net>
Content-Type: text/plain; charset=3Diso-8859-1
On Tue, Sep 15, 2009 at 08:43:10AM +0200, D??ster Horst wrote:
> I do have 2 simple linestrings which partly have the same geometry =
(see
> attached dump and image). My aim is to detect the common part of these
> two linestrings with the following query:
>=20
> select st_intersection(a.the_geom, b.the_geom), a.myid
> from aline1, aline2
>=20
> But as the result I get a geometrycollection with one POINT and one
> LINESTING. I expected one linestring. Is it a bug??
> I'll appreciate any hints.
Could be a bug, or a dimensional collapse.
Try posting the geometries, Martin Davis might have a try
with JTS.
Try also upgrading GEOS.
--strk;
>=20
> My system:
> Postgis 1.4.0
> Geos 3.1.1
> Postgres 8.3.6
>=20
> Best regards
> =20
> Dr. Horst D?ster
> Stv. Amtschef / GIS-Koordinator=20
>=20
> Kanton Solothurn
> Bau- und Justizdepartement
> Amt f?r Geoinformation
> SO!GIS Koordination
> R?tistrasse 4
> CH-4501 Solothurn
>=20
> Telefon ++41(0)32 627 25 32
> Telefax ++41(0)32 627 22 14
>=20
> mailto:horst.duester at bd.so.ch
> http://www.agi.so.ch
>=20
>=20
> _______________________________________________
> postgis-users mailing list
> postgis-users at postgis.refractions.net
> http://postgis.refractions.net/mailman/listinfo/postgis-users
--=20
Free GIS & Flash consultant/developer () ASCII Ribbon Campaign
http://foo.keybit.net/~strk/services.html /\ Keep it simple!=20
------------------------------
Message: 17
Date: Tue, 15 Sep 2009 11:32:48 +0200
From: D?ster Horst <Horst.Duester at bd.so.ch>
Subject: Re: [postgis-users] st_intersection error
To: postgis-users <postgis-users at postgis.refractions.net>, strk
<strk at keybit.net>
Message-ID: <H00002eb058418dd.1253007168.srsofaioi6145.ktso.ch at MHS>
Content-Type: text/plain; charset=3D"iso-8859-1"
strk
I attached the geometries in my initial mail.
Take a look at aline1.sql and aline2.sql
regards
Horst
------------------------------------------------
Dr. Horst D?ster
Stv. Amtschef / GIS-Koordinator=20
Kanton Solothurn
Bau- und Justizdepartement
Amt f?r Geoinformation
SO!GIS Koordination
R?tistrasse 4
CH-4501 Solothurn
Telefon ++41(0)32 627 25 32
Telefax ++41(0)32 627 22 14
mailto:horst.duester at bd.so.ch
http://www.agi.so.ch
-----Urspr?ngliche Nachricht-----
Von: strk [mailto:strk at keybit.net]
Gesendet am: Dienstag, 15. September 2009 10:27
An: PostGIS Users Discussion
Betreff: Re: [postgis-users] st_intersection error
On Tue, Sep 15, 2009 at 08:43:10AM +0200, D??ster Horst wrote:
> I do have 2 simple linestrings which partly have the same geometry
(see
> attached dump and image). My aim is to detect the common part of these
> two linestrings with the following query:
>=20
> select st_intersection(a.the_geom, b.the_geom), a.myid
> from aline1, aline2
>=20
> But as the result I get a geometrycollection with one POINT and one
> LINESTING. I expected one linestring. Is it a bug??
> I'll appreciate any hints.
Could be a bug, or a dimensional collapse.
Try posting the geometries, Martin Davis might have a try
with JTS.
Try also upgrading GEOS.
--strk;
>=20
> My system:
> Postgis 1.4.0
> Geos 3.1.1
> Postgres 8.3.6
>=20
> Best regards
> =20
> Dr. Horst D?ster
> Stv. Amtschef / GIS-Koordinator=20
>=20
> Kanton Solothurn
> Bau- und Justizdepartement
> Amt f?r Geoinformation
> SO!GIS Koordination
> R?tistrasse 4
> CH-4501 Solothurn
>=20
> Telefon ++41(0)32 627 25 32
> Telefax ++41(0)32 627 22 14
>=20
> mailto:horst.duester at bd.so.ch
> http://www.agi.so.ch
>=20
>=20
> _______________________________________________
> postgis-users mailing list
> postgis-users at postgis.refractions.net
> http://postgis.refractions.net/mailman/listinfo/postgis-users
--=20
Free GIS & Flash consultant/developer () ASCII Ribbon Campaign
http://foo.keybit.net/~strk/services.html /\ Keep it simple!=20
_______________________________________________
postgis-users mailing list
postgis-users at postgis.refractions.net
http://postgis.refractions.net/mailman/listinfo/postgis-users
-------------- next part --------------
An HTML attachment was scrubbed...
URL: =
<http://postgis.refractions.net/pipermail/postgis-users/attachments/20090=
915/202c8a98/attachment-0001.html>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: aline1.sql
Type: application/octet-stream
Size: 1199 bytes
Desc: not available
URL: =
<http://postgis.refractions.net/pipermail/postgis-users/attachments/20090=
915/202c8a98/attachment-0002.obj>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: aline2.sql
Type: application/octet-stream
Size: 1268 bytes
Desc: not available
URL: =
<http://postgis.refractions.net/pipermail/postgis-users/attachments/20090=
915/202c8a98/attachment-0003.obj>
------------------------------
Message: 18
Date: Tue, 15 Sep 2009 11:53:54 +0200
From: strk <strk at keybit.net>
Subject: Re: [postgis-users] st_intersection error
To: PostGIS Users Discussion <postgis-users at postgis.refractions.net>
Message-ID: <20090915095354.GC30456 at keybit.net>
Content-Type: text/plain; charset=3Diso-8859-1
On Tue, Sep 15, 2009 at 11:32:48AM +0200, D??ster Horst wrote:
> strk
>=20
> I attached the geometries in my initial mail.
>=20
> Take a look at aline1.sql and aline2.sql
I confirm a POINT and a LINESTRING in a COLLECTION
is returned by ST_Intersection() with latest version:
POSTGIS=3D"1.5.0SVN" GEOS=3D"3.2.0-CAPI-1.6.0" PROJ=3D"Rel. 4.6.0, 21 =
Dec 2007" USE_STATS
Martin, how about JTS ?
--strk;=20
Free GIS & Flash consultant/developer () ASCII Ribbon Campaign
http://foo.keybit.net/~strk/services.html /\ Keep it simple!=20
------------------------------
Message: 19
Date: Tue, 15 Sep 2009 12:50:09 +0200
From: nicklas.aven at jordogskog.no
Subject: Re: [postgis-users] st_intersection error
To: PostGIS Users Discussion <postgis-users at postgis.refractions.net>
Message-ID: <200909151050.n8FAo9DN018385 at mail-core.space2u.com>
Content-Type: text/plain; charset=3D"iso-8859-1"
The problem is the vertex represented only in line 2 between your =
resulting line and point 'POINT(618971.85092967 251568.179431664) . =
Because of precision aspects it will never be exactly on the line =
causing st_intersection find only the point and line to be exactly the =
same place. /Nicklas
2009-09-15 D??ster Horst wrote:
>
I do have 2 simple linestrings which partly have the same geometry (see =
attached dump and image). My aim is to detect the common part of these =
two linestrings with the following query:
>>
select st_intersection(a.the_geom, b.the_geom), a.myid
>from aline1, aline2
>>
But as the result I get a geometrycollection with one POINT and one =
LINESTING. I expected one linestring. Is it a bug??
>I'll appreciate any hints.
>>
My system:
>Postgis 1.4.0
>Geos 3.1.1
>Postgres 8.3.6
>>
Best regards
>>
Dr. Horst D?ster
>Stv. Amtschef / GIS-Koordinator=20
>>
Kanton Solothurn
>Bau- und Justizdepartement
>Amt f?r Geoinformation
>SO!GIS Koordination
>R?tistrasse 4
>CH-4501 Solothurn
>>
Telefon ++41(0)32 627 25 32
>Telefax ++41(0)32 627 22 14
>>
mailto:horst.duester at bd.so.ch
>www.agi.so.ch
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: =
<http://postgis.refractions.net/pipermail/postgis-users/attachments/20090=
915/ee7cc43b/attachment-0001.html>
------------------------------
Message: 20
Date: Tue, 15 Sep 2009 08:23:52 -0700
From: Kevin Neufeld <kneufeld at refractions.net>
Subject: Re: [postgis-users] No Primary Key
To: punkish at eidesis.org, PostGIS Users Discussion
<postgis-users at postgis.refractions.net>
Message-ID: <4AAFB188.6080303 at refractions.net>
Content-Type: text/plain; charset=3DISO-8859-1; format=3Dflowed
P Kishor wrote:
>> create table coal AS SELECT * from geology where group_=3D'GONDWANA';
>>
>> But the resultant table coal is not having a 'primary key'.
>> table geology has gid as primary key.
>>
>> Pl give me a way how I can get a primary key on table coal.
>>
ALTER TABLE coal ADD PRIMARY KEY (gid);
ANALYZE coal;
------------------------------
Message: 21
Date: Tue, 15 Sep 2009 08:54:03 -0700
From: Martin Davis <mbdavis at refractions.net>
Subject: Re: [postgis-users] st_intersection error
To: PostGIS Users Discussion <postgis-users at postgis.refractions.net>,
Martin Davis <mbdavis at refractions.net>
Message-ID: <4AAFB89B.9080709 at refractions.net>
Content-Type: text/plain; charset=3DISO-8859-1; format=3Dflowed
Same result in JTS.
As Nicklas points out, this is a problem of precision. The software is=20
working as designed, but it's not designed to use it to perform=20
intersections with a tolerance.
Some more sophisticated approach is required, such as vertex snapping=20
the lines to a given tolerance, or perhaps using a narrow buffer and=20
some sort of segment-based containment tests.
strk wrote:
> On Tue, Sep 15, 2009 at 11:32:48AM +0200, D??ster Horst wrote:
> =20
>> strk
>>
>> I attached the geometries in my initial mail.
>>
>> Take a look at aline1.sql and aline2.sql
>> =20
>
> I confirm a POINT and a LINESTRING in a COLLECTION
> is returned by ST_Intersection() with latest version:
>
> POSTGIS=3D"1.5.0SVN" GEOS=3D"3.2.0-CAPI-1.6.0" PROJ=3D"Rel. 4.6.0, 21 =
Dec 2007" USE_STATS
>
> Martin, how about JTS ?
>
> --strk;=20
>
> Free GIS & Flash consultant/developer () ASCII Ribbon Campaign
> http://foo.keybit.net/~strk/services.html /\ Keep it simple!=20
>
> =20
--=20
Martin Davis
Senior Technical Architect
Refractions Research, Inc.
(250) 383-3022
------------------------------
Message: 22
Date: Tue, 15 Sep 2009 19:20:38 +0200
From: Stefan Keller <sfkeller at gmail.com>
Subject: Re: [postgis-users] PostGIS to DXF output
To: PostGIS Users Discussion <postgis-users at postgis.refractions.net>
Message-ID:
<25bc040b0909151020r2199b87o79d85310c645f60 at mail.gmail.com>
Content-Type: text/plain; charset=3DISO-8859-1
I think that Frank Wamerdam is working on a DXF in/out driver for OGR
without dependencies.
Not sure if it's out before end of year.
If you anyone has any DXF test samples that would be helpful.
Yours, S.
2009/9/9 Paul Ramsey <pramsey at cleverelephant.ca>:
> Not that I know of opensource-wise. FMEServer from Safe Software
> (www.safe.com) will give you every format under the sun and can read
> from PostGIS.
>
> P.
>
> On Tue, Sep 8, 2009 at 5:50 PM, GIS<GIS at asiaq.gl> wrote:
>>
>> Hi, I have a question concerning DXF/DWG output.
>> On our mapping website http://en.nunagis.gl users can download our =
postgis
>> data as TAB, DGN, SHP, GML and other formats supported by ogr. =
However, we
>> would like them to be able to download dxf files, too. As far as I =
know,
>> this is not possible through ogr.
>> Any other possibilities at hand?
>>
>> Karl
>>
>>
>> _______________________________________________
>> 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 mailing list
postgis-users at postgis.refractions.net
http://postgis.refractions.net/mailman/listinfo/postgis-users
End of postgis-users Digest, Vol 85, Issue 15
*********************************************