tag:old.nabble.com,2006:forum-1223Nabble - PostGIS - User2009-11-23T22:44:06Ztag:old.nabble.com,2006:post-26491153Re: identifying blank geometry fields2009-11-23T22:44:06Z2009-11-23T22:44:06ZNicklas Avén<html>
<head>
<title></title>
<meta name="GENERATOR" content="MSHTML 8.00.6001.18852"></meta>
<meta name="SKYPE_FRAMEID" content="HBROAMAKDT"></meta>
<meta id="skype_v3_tb_marker_id" name="SKYPE_PARSING_HAS_FINISHED" content="metacontent"></meta>
</head>
<body>
<div align="left">Hallo </div>
<div align="left"> </div>
<div align="left">At least in pgAdmin the cells is showed as empty cells if the geometries is above some size. But they are there.</div>
<div align="left">An easy way to confirm they are there is to use som function on them. If it is polygons for instance you can run st_area or st_length for lines. </div>
<div align="left"> </div>
<div align="left">Hope that helps</div>
<div align="left">/Nicklas<br />
<br />
2009-11-24 Bob Hume wrote:<br />
<br />
></div>
<div>select count(*) from dum where the_geom is null gives a result of 0 <br />
><br />
> select count(the_geom) from dum gives a result is 1676, which is the total number of rows in the table <br />
><br />
> This is consistent with my previous description. The apparently blank the_geom fields apparently have some data, but the data doesn't display. This seems odd, and makes it difficult (impossible for me, so far) to identify the blank fields. <br />
><br />
><br />
> (Sorry this seems to be a duplicate. I couldn't figure out how to reply to a message in a thread and assumed the mail program would figure it out. It didn't. Just for the heck of it I clicked on the name of the person I was replying to with the idea that I'd send a message directly to him, and a mail window opened to send a message to the list. Who knew? I searched, and couldn't find any information explaining this.)<br />
></div>
</body>
</html>
<br />_______________________________________________
<br>postgis-users mailing list
<br><a href="http://old.nabble.com/user/SendEmail.jtp?type=post&post=26491153&i=0" target="_top" rel="nofollow">postgis-users@...</a>
<br><a href="http://postgis.refractions.net/mailman/listinfo/postgis-users" target="_top" rel="nofollow">http://postgis.refractions.net/mailman/listinfo/postgis-users</a><br>tag:old.nabble.com,2006:post-26490812Re: identifying blank geometry fields2009-11-23T21:49:56Z2009-11-23T21:49:56ZBob Hume<!DOCTYPE html PUBLIC "-//W3C//DTD HTML 4.01 Transitional//EN">
<html>
<head>
</head>
<body bgcolor="#ffffff" text="#000000">
<div class="moz-text-flowed" style="font-family: -moz-fixed; font-size: 13px;" lang="x-western">select
count(*) from dum where the_geom is null gives a result of 0
<br>
<br>
select count(the_geom) from dum gives a result is 1676, which is the
total number of rows in the table
<br>
<br>
This is consistent with my previous description. The apparently blank
the_geom fields apparently have some data, but the data doesn't
display. This seems odd, and makes it difficult (impossible for me, so
far) to identify the blank fields.
<br>
<br>
<br>
(Sorry this seems to be a duplicate. I couldn't figure out how to
reply to a message in a thread and assumed the mail program would
figure it out. It didn't. Just for the heck of it I clicked on the
name of the person I was replying to with the idea that I'd send a
message directly to him, and a mail window opened to send a message to
the list. Who knew? I searched, and couldn't find any information
explaining this.)<br>
</div>
</body>
</html>
<br />_______________________________________________
<br>postgis-users mailing list
<br><a href="http://old.nabble.com/user/SendEmail.jtp?type=post&post=26490812&i=0" target="_top" rel="nofollow">postgis-users@...</a>
<br><a href="http://postgis.refractions.net/mailman/listinfo/postgis-users" target="_top" rel="nofollow">http://postgis.refractions.net/mailman/listinfo/postgis-users</a><br>tag:old.nabble.com,2006:post-26490796Re: identifying blank geometry fields2009-11-23T21:46:51Z2009-11-23T21:46:51ZStephen WoodbridgeBob Hume wrote:
<br>> select count(*) from dum where the_geom is null gives a result of 0
<br>>
<br>> select count(the_geom) from dum gives a result is 1676, which is the
<br>> total number of rows in the table
<br>>
<br>> This is consistent with my previous description. The apparently blank
<br>> the_geom fields apparently have some data, but the data doesn't
<br>> display. This seems odd, and makes it difficult (impossible for me, so
<br>> far) to identify the blank fields.
<br><br>Have you run:
<br><br>select count(*) from dum where not isvalid(the_geom);
<br><br>-Steve
<br>_______________________________________________
<br>postgis-users mailing list
<br><a href="http://old.nabble.com/user/SendEmail.jtp?type=post&post=26490796&i=0" target="_top" rel="nofollow">postgis-users@...</a>
<br><a href="http://postgis.refractions.net/mailman/listinfo/postgis-users" target="_top" rel="nofollow">http://postgis.refractions.net/mailman/listinfo/postgis-users</a><br>tag:old.nabble.com,2006:post-26490704identifying blank geometry fields2009-11-23T21:34:53Z2009-11-23T21:34:53ZBob Humeselect count(*) from dum where the_geom is null gives a result of 0
<br><br>select count(the_geom) from dum gives a result is 1676, which is the
<br>total number of rows in the table
<br><br>This is consistent with my previous description. The apparently blank
<br>the_geom fields apparently have some data, but the data doesn't
<br>display. This seems odd, and makes it difficult (impossible for me, so
<br>far) to identify the blank fields.
<br><br>_______________________________________________
<br>postgis-users mailing list
<br><a href="http://old.nabble.com/user/SendEmail.jtp?type=post&post=26490704&i=0" target="_top" rel="nofollow">postgis-users@...</a>
<br><a href="http://postgis.refractions.net/mailman/listinfo/postgis-users" target="_top" rel="nofollow">http://postgis.refractions.net/mailman/listinfo/postgis-users</a><br>tag:old.nabble.com,2006:post-26490593Re: loading many similar shapefiles2009-11-23T21:14:02Z2009-11-23T21:14:02ZPaul Ramsey-3Just use a temp table and interleave the load steps with psql -c
<br>"insert into realtable select a,b,c from tmptable"
<br><br>P
<br><br>On Mon, Nov 23, 2009 at 8:27 PM, P Kishor <<a href="http://old.nabble.com/user/SendEmail.jtp?type=post&post=26490593&i=0" target="_top" rel="nofollow">punk.kish@...</a>> wrote:
<div class='shrinkable-quote'><br>> I have many shapefiles that are conceptually identical in structure,
<br>> except a couple of columns in each file are named differently. These
<br>> columns are unimportant.
<br>>
<br>> How can I load all these files in the same table without the offending columns?
<br>>
<br>> # generate the create table and insert commands
<br>> shp2pgsql -s 4269 -c file1 table > load.sql
<br>>
<br>> # generate the insert commands for subsequent files
<br>> shp2pgsql -s 4269 -a file2 table >> load.sql
<br>> shp2pgsql -s 4269 -a file3 table >> load.sql
<br>> ..
<br>>
<br>> The above commands generate INSERT statements like so --
<br>>
<br>> INSERT INTO "table" (
<br>> "area",
<br>> "perimeter",
<br>> "file1_",
<br>> "file1_i",
<br>> the_geom
<br>> ) VALUES (
<br>> '2.29080994500086E-04',
<br>> '6.34477469815224E-02',
<br>> '0106000020AD100000010000000103000020AD100000010000000D000000A4552DE928EA55C0C6F99B5088443E404582A966D6E955C050C8CEDBD8443E40187D0569C6E955C02B508BC1C3443E403DEFC682C2E955C0BB4560AC6F443E40EB8D5A61FAE955C0A6D24F38BB413E402907B30930EA55C0503750E09D403E40C28A53AD85EA55C0A8716F7EC3403E403D997FF44DEB55C058FFE7305F423E402C4833164DEB55C00F4240BE84423E40DAAD65321CEB55C0AFCF9CF529433E40D48041D2A7EA55C046274BADF7433E40FB58C16F43EA55C045F295404A443E40A4552DE928EA55C0C6F99B5088443E40'
<br>> );
<br>>
<br>>
<br>> I don't want file1_ and file1_i and their corresponding values. Is
<br>> there a way to specify which columns to load?
<br>>
<br>>
<br>>
<br>> --
<br>> Puneet Kishor
<br>> _______________________________________________
<br>> postgis-users mailing list
<br>> <a href="http://old.nabble.com/user/SendEmail.jtp?type=post&post=26490593&i=1" target="_top" rel="nofollow">postgis-users@...</a>
<br>> <a href="http://postgis.refractions.net/mailman/listinfo/postgis-users" target="_top" rel="nofollow">http://postgis.refractions.net/mailman/listinfo/postgis-users</a><br>>
</div>_______________________________________________
<br>postgis-users mailing list
<br><a href="http://old.nabble.com/user/SendEmail.jtp?type=post&post=26490593&i=2" target="_top" rel="nofollow">postgis-users@...</a>
<br><a href="http://postgis.refractions.net/mailman/listinfo/postgis-users" target="_top" rel="nofollow">http://postgis.refractions.net/mailman/listinfo/postgis-users</a><br>tag:old.nabble.com,2006:post-26490333loading many similar shapefiles2009-11-23T20:27:31Z2009-11-23T20:27:31ZP Kishor-3I have many shapefiles that are conceptually identical in structure,
<br>except a couple of columns in each file are named differently. These
<br>columns are unimportant.
<br><br>How can I load all these files in the same table without the offending columns?
<br><br># generate the create table and insert commands
<br>shp2pgsql -s 4269 -c file1 table > load.sql
<br><br># generate the insert commands for subsequent files
<br>shp2pgsql -s 4269 -a file2 table >> load.sql
<br>shp2pgsql -s 4269 -a file3 table >> load.sql
<br>..
<br><br>The above commands generate INSERT statements like so --
<br><br>INSERT INTO "table" (
<br> "area",
<br> "perimeter",
<br> "file1_",
<br> "file1_i",
<br> the_geom
<br>) VALUES (
<br> '2.29080994500086E-04',
<br> '6.34477469815224E-02',
<br> '0106000020AD100000010000000103000020AD100000010000000D000000A4552DE928EA55C0C6F99B5088443E404582A966D6E955C050C8CEDBD8443E40187D0569C6E955C02B508BC1C3443E403DEFC682C2E955C0BB4560AC6F443E40EB8D5A61FAE955C0A6D24F38BB413E402907B30930EA55C0503750E09D403E40C28A53AD85EA55C0A8716F7EC3403E403D997FF44DEB55C058FFE7305F423E402C4833164DEB55C00F4240BE84423E40DAAD65321CEB55C0AFCF9CF529433E40D48041D2A7EA55C046274BADF7433E40FB58C16F43EA55C045F295404A443E40A4552DE928EA55C0C6F99B5088443E40'
<br>);
<br><br><br>I don't want file1_ and file1_i and their corresponding values. Is
<br>there a way to specify which columns to load?
<br><br><br><br>--
<br>Puneet Kishor
<br>_______________________________________________
<br>postgis-users mailing list
<br><a href="http://old.nabble.com/user/SendEmail.jtp?type=post&post=26490333&i=0" target="_top" rel="nofollow">postgis-users@...</a>
<br><a href="http://postgis.refractions.net/mailman/listinfo/postgis-users" target="_top" rel="nofollow">http://postgis.refractions.net/mailman/listinfo/postgis-users</a><br>tag:old.nabble.com,2006:post-26490170Re: Using a rectangular selection2009-11-23T20:00:41Z2009-11-23T20:00:41ZBrian Modra2009/11/23 Henri De Feraudy <<a href="http://old.nabble.com/user/SendEmail.jtp?type=post&post=26490170&i=0" target="_top" rel="nofollow">feraudyh@...</a>>
<br>>
<br>> Hi,
<br>> I'm working in planar coordinates. I have a Postgis database and
<br>> for starters I want to limit my attention to all geometries that intersect a given rectangle (the "selection").
<br>> I can imagine that one could enumerate all geometries and ask in each case if they
<br>> intersect the rectangle, but that doesnt seem very efficient, so what do you suggest?
<br><br>Use ST_Intersects
<br><div class='shrinkable-quote'><br>> I only need general guidance so I can look at the various parts of the doc.
<br>> Later I'll be wanting to clip the geometies with this selection.
<br>> Still later I'll want to consider polygonal selections.
<br>> By the way I'm accessing the PostGis database via LibPQ and have written my application
<br>> in Qt.
<br>>
<br>> Thanking you in advance
<br>> Henri
<br>>
<br>> _______________________________________________
<br>> postgis-users mailing list
<br>> <a href="http://old.nabble.com/user/SendEmail.jtp?type=post&post=26490170&i=1" target="_top" rel="nofollow">postgis-users@...</a>
<br>> <a href="http://postgis.refractions.net/mailman/listinfo/postgis-users" target="_top" rel="nofollow">http://postgis.refractions.net/mailman/listinfo/postgis-users</a><br>>
</div><br><br><br>--
<br>Brian Modra Land line: +27 23 5411 462
<br>Mobile: +27 79 69 77 082
<br>5 Jan Louw Str, Prince Albert, 6930
<br>Postal: P.O. Box 2, Prince Albert 6930
<br>South Africa
<br><a href="http://www.zwartberg.com/" target="_top" rel="nofollow">http://www.zwartberg.com/</a><br>_______________________________________________
<br>postgis-users mailing list
<br><a href="http://old.nabble.com/user/SendEmail.jtp?type=post&post=26490170&i=2" target="_top" rel="nofollow">postgis-users@...</a>
<br><a href="http://postgis.refractions.net/mailman/listinfo/postgis-users" target="_top" rel="nofollow">http://postgis.refractions.net/mailman/listinfo/postgis-users</a><br>tag:old.nabble.com,2006:post-26487571Re: C++ and Postgis2009-11-23T14:53:41Z2009-11-23T14:53:41ZPedro Doria MeunierPaul,
<br><br>I could just *kiss* you mate! :D
<br>(Right now I'm feeling pretty ashamed for not thinking of it right away
<br>... :] guess it's been too many hours coding... )
<br><br>All the best,
<br>Pedro.
<br><br><br>On 11/23/2009 10:45 PM, Paul Ramsey wrote:
<div class='shrinkable-quote'><br>> If it's just points, do it all on the server side, skip encoding formats
<br>>
<br>> select st_x(pt) as x, st_y(pt) as y from mytable where myclause();
<br>>
<br>> P.
<br>>
<br>> On Mon, Nov 23, 2009 at 2:43 PM, Pedro Doria Meunier
<br>> <<a href="http://old.nabble.com/user/SendEmail.jtp?type=post&post=26487571&i=0" target="_top" rel="nofollow">pdoria@...</a>> wrote:
<br>>
<br>>> Thank you Paul and John for your ideas. :)
<br>>>
<br>>> Actually what I'm trying to do is to decode simple POINT geoms...
<br>>>
<br>>> This is what I've done in PHP (using the astext() approach)
<br>>>
<br>>> // extract the geometry
<br>>> $geometry = str_replace("POINT(","", $loc);
<br>>> $geometry = str_replace(")","",$geometry);
<br>>> $coordinate = explode(" ",$geometry);
<br>>> $lon=$coordinate[0];
<br>>> $lat=$coordinate[1];
<br>>>
<br>>> Now this was for the web interface... ( the user can wait another 100ms or
<br>>> so ... :] )
<br>>> As far as a server solution is concerned this *is* a time-consuming
<br>>> procedure... especially so when you've got *thousands* of connections per
<br>>> minute...
<br>>>
<br>>> Furthermore, I'm developing the "beast" in Qt4 and - *ahem* - had never
<br>>> needed to use libpq so I'm a bit lost here... :)
<br>>>
<br>>> Could someone be so kind as to point to a code example (outside using
<br>>> astext() method) ? :]
<br>>>
<br>>> P.S. @John - found this to be interesting as I could store it in QByteArray
<br>>> and go from there... ;) -- the real question here is: how to interpret the
<br>>> byte array?
<br>>>
<br>>> BR,
<br>>> Pedro.
<br>>>
<br>>>
<br>>> On 11/23/2009 09:50 PM, Paul Ramsey wrote:
<br>>>
<br>>>> That depends entirely on what you're planning to do with them once you
<br>>>> get them. You'll probably want to serialize them into a known form so
<br>>>> you can use them on the client side. ST_AsText() is one, very heavy,
<br>>>> option. Mapserver, to take an example, uses
<br>>>> encode(ST_AsBinary(ST_Force2D(geom)), 'hex'), so what goes over the
<br>>>> wire is hex strings, but those are easy to turn into binary arrays, in
<br>>>> WKB format, suitable for reading by a number of C++ libraries (like
<br>>>> OGR or GEOG or etc).
<br>>>>
<br>>>> Best,
<br>>>>
<br>>>> P.
<br>>>>
<br>>>> On Mon, Nov 23, 2009 at 1:19 PM, Pedro Doria Meunier
<br>>>> <<a href="http://old.nabble.com/user/SendEmail.jtp?type=post&post=26487571&i=1" target="_top" rel="nofollow">pdoria@...</a>> wrote:
<br>>>>
<br>>>>
<br>>>>> Hi All,
<br>>>>>
<br>>>>> First off sorry for cross-posting. I'm hoping for a bigger audience... :]
<br>>>>>
<br>>>>> As for my question:
<br>>>>>
<br>>>>> How does one deal with Postgis geom fields in C++ ?
<br>>>>> Google hasn't offered much help so far.
<br>>>>>
<br>>>>> I know I can make a "SELECT astext(my_geom) FROM my_table" and after
<br>>>>> parsing
<br>>>>> the resulting text.
<br>>>>> This seems a bit awkward and time-consuming...
<br>>>>>
<br>>>>> Is there a better way to deal with it?
<br>>>>>
<br>>>>> Already thankful for any ideas,
<br>>>>> BR,
<br>>>>> Pedro.
<br>>>>> _______________________________________________
<br>>>>> postgis-users mailing list
<br>>>>> <a href="http://old.nabble.com/user/SendEmail.jtp?type=post&post=26487571&i=2" target="_top" rel="nofollow">postgis-users@...</a>
<br>>>>> <a href="http://postgis.refractions.net/mailman/listinfo/postgis-users" target="_top" rel="nofollow">http://postgis.refractions.net/mailman/listinfo/postgis-users</a><br>>>>>
<br>>>>>
<br>>>>>
<br>>>> _______________________________________________
<br>>>> postgis-users mailing list
<br>>>> <a href="http://old.nabble.com/user/SendEmail.jtp?type=post&post=26487571&i=3" target="_top" rel="nofollow">postgis-users@...</a>
<br>>>> <a href="http://postgis.refractions.net/mailman/listinfo/postgis-users" target="_top" rel="nofollow">http://postgis.refractions.net/mailman/listinfo/postgis-users</a><br>>>>
<br>>>>
<br>>>>
<br>>> _______________________________________________
<br>>> postgis-users mailing list
<br>>> <a href="http://old.nabble.com/user/SendEmail.jtp?type=post&post=26487571&i=4" target="_top" rel="nofollow">postgis-users@...</a>
<br>>> <a href="http://postgis.refractions.net/mailman/listinfo/postgis-users" target="_top" rel="nofollow">http://postgis.refractions.net/mailman/listinfo/postgis-users</a><br>>>
<br>>>
<br>> _______________________________________________
<br>> postgis-users mailing list
<br>> <a href="http://old.nabble.com/user/SendEmail.jtp?type=post&post=26487571&i=5" target="_top" rel="nofollow">postgis-users@...</a>
<br>> <a href="http://postgis.refractions.net/mailman/listinfo/postgis-users" target="_top" rel="nofollow">http://postgis.refractions.net/mailman/listinfo/postgis-users</a><br>>
<br>>
</div>_______________________________________________
<br>postgis-users mailing list
<br><a href="http://old.nabble.com/user/SendEmail.jtp?type=post&post=26487571&i=6" target="_top" rel="nofollow">postgis-users@...</a>
<br><a href="http://postgis.refractions.net/mailman/listinfo/postgis-users" target="_top" rel="nofollow">http://postgis.refractions.net/mailman/listinfo/postgis-users</a><br>tag:old.nabble.com,2006:post-26487478Re: C++ and Postgis2009-11-23T14:45:01Z2009-11-23T14:45:01ZPaul Ramsey-4If it's just points, do it all on the server side, skip encoding formats
<br><br>select st_x(pt) as x, st_y(pt) as y from mytable where myclause();
<br><br>P.
<br><br>On Mon, Nov 23, 2009 at 2:43 PM, Pedro Doria Meunier
<br><<a href="http://old.nabble.com/user/SendEmail.jtp?type=post&post=26487478&i=0" target="_top" rel="nofollow">pdoria@...</a>> wrote:
<div class='shrinkable-quote'><br>> Thank you Paul and John for your ideas. :)
<br>>
<br>> Actually what I'm trying to do is to decode simple POINT geoms...
<br>>
<br>> This is what I've done in PHP (using the astext() approach)
<br>>
<br>> // extract the geometry
<br>> $geometry = str_replace("POINT(","", $loc);
<br>> $geometry = str_replace(")","",$geometry);
<br>> $coordinate = explode(" ",$geometry);
<br>> $lon=$coordinate[0];
<br>> $lat=$coordinate[1];
<br>>
<br>> Now this was for the web interface... ( the user can wait another 100ms or
<br>> so ... :] )
<br>> As far as a server solution is concerned this *is* a time-consuming
<br>> procedure... especially so when you've got *thousands* of connections per
<br>> minute...
<br>>
<br>> Furthermore, I'm developing the "beast" in Qt4 and - *ahem* - had never
<br>> needed to use libpq so I'm a bit lost here... :)
<br>>
<br>> Could someone be so kind as to point to a code example (outside using
<br>> astext() method) ? :]
<br>>
<br>> P.S. @John - found this to be interesting as I could store it in QByteArray
<br>> and go from there... ;) -- the real question here is: how to interpret the
<br>> byte array?
<br>>
<br>> BR,
<br>> Pedro.
<br>>
<br>>
<br>> On 11/23/2009 09:50 PM, Paul Ramsey wrote:
<br>>>
<br>>> That depends entirely on what you're planning to do with them once you
<br>>> get them. You'll probably want to serialize them into a known form so
<br>>> you can use them on the client side. ST_AsText() is one, very heavy,
<br>>> option. Mapserver, to take an example, uses
<br>>> encode(ST_AsBinary(ST_Force2D(geom)), 'hex'), so what goes over the
<br>>> wire is hex strings, but those are easy to turn into binary arrays, in
<br>>> WKB format, suitable for reading by a number of C++ libraries (like
<br>>> OGR or GEOG or etc).
<br>>>
<br>>> Best,
<br>>>
<br>>> P.
<br>>>
<br>>> On Mon, Nov 23, 2009 at 1:19 PM, Pedro Doria Meunier
<br>>> <<a href="http://old.nabble.com/user/SendEmail.jtp?type=post&post=26487478&i=1" target="_top" rel="nofollow">pdoria@...</a>> wrote:
<br>>>
<br>>>>
<br>>>> Hi All,
<br>>>>
<br>>>> First off sorry for cross-posting. I'm hoping for a bigger audience... :]
<br>>>>
<br>>>> As for my question:
<br>>>>
<br>>>> How does one deal with Postgis geom fields in C++ ?
<br>>>> Google hasn't offered much help so far.
<br>>>>
<br>>>> I know I can make a "SELECT astext(my_geom) FROM my_table" and after
<br>>>> parsing
<br>>>> the resulting text.
<br>>>> This seems a bit awkward and time-consuming...
<br>>>>
<br>>>> Is there a better way to deal with it?
<br>>>>
<br>>>> Already thankful for any ideas,
<br>>>> BR,
<br>>>> Pedro.
<br>>>> _______________________________________________
<br>>>> postgis-users mailing list
<br>>>> <a href="http://old.nabble.com/user/SendEmail.jtp?type=post&post=26487478&i=2" target="_top" rel="nofollow">postgis-users@...</a>
<br>>>> <a href="http://postgis.refractions.net/mailman/listinfo/postgis-users" target="_top" rel="nofollow">http://postgis.refractions.net/mailman/listinfo/postgis-users</a><br>>>>
<br>>>>
<br>>>
<br>>> _______________________________________________
<br>>> postgis-users mailing list
<br>>> <a href="http://old.nabble.com/user/SendEmail.jtp?type=post&post=26487478&i=3" target="_top" rel="nofollow">postgis-users@...</a>
<br>>> <a href="http://postgis.refractions.net/mailman/listinfo/postgis-users" target="_top" rel="nofollow">http://postgis.refractions.net/mailman/listinfo/postgis-users</a><br>>>
<br>>>
<br>>
<br>> _______________________________________________
<br>> postgis-users mailing list
<br>> <a href="http://old.nabble.com/user/SendEmail.jtp?type=post&post=26487478&i=4" target="_top" rel="nofollow">postgis-users@...</a>
<br>> <a href="http://postgis.refractions.net/mailman/listinfo/postgis-users" target="_top" rel="nofollow">http://postgis.refractions.net/mailman/listinfo/postgis-users</a><br>>
</div>_______________________________________________
<br>postgis-users mailing list
<br><a href="http://old.nabble.com/user/SendEmail.jtp?type=post&post=26487478&i=5" target="_top" rel="nofollow">postgis-users@...</a>
<br><a href="http://postgis.refractions.net/mailman/listinfo/postgis-users" target="_top" rel="nofollow">http://postgis.refractions.net/mailman/listinfo/postgis-users</a><br>tag:old.nabble.com,2006:post-26487456Re: C++ and Postgis2009-11-23T14:43:07Z2009-11-23T14:43:07ZPedro Doria MeunierThank you Paul and John for your ideas. :)
<br><br>Actually what I'm trying to do is to decode simple POINT geoms...
<br><br>This is what I've done in PHP (using the astext() approach)
<br><br> // extract the geometry
<br> $geometry = str_replace("POINT(","", $loc);
<br> $geometry = str_replace(")","",$geometry);
<br> $coordinate = explode(" ",$geometry);
<br> $lon=$coordinate[0];
<br> $lat=$coordinate[1];
<br><br>Now this was for the web interface... ( the user can wait another 100ms
<br>or so ... :] )
<br>As far as a server solution is concerned this *is* a time-consuming
<br>procedure... especially so when you've got *thousands* of connections
<br>per minute...
<br><br>Furthermore, I'm developing the "beast" in Qt4 and - *ahem* - had never
<br>needed to use libpq so I'm a bit lost here... :)
<br><br>Could someone be so kind as to point to a code example (outside using
<br>astext() method) ? :]
<br><br>P.S. @John - found this to be interesting as I could store it in
<br>QByteArray and go from there... ;) -- the real question here is: how to
<br>interpret the byte array?
<br><br>BR,
<br>Pedro.
<br><br><br>On 11/23/2009 09:50 PM, Paul Ramsey wrote:
<div class='shrinkable-quote'><br>> That depends entirely on what you're planning to do with them once you
<br>> get them. You'll probably want to serialize them into a known form so
<br>> you can use them on the client side. ST_AsText() is one, very heavy,
<br>> option. Mapserver, to take an example, uses
<br>> encode(ST_AsBinary(ST_Force2D(geom)), 'hex'), so what goes over the
<br>> wire is hex strings, but those are easy to turn into binary arrays, in
<br>> WKB format, suitable for reading by a number of C++ libraries (like
<br>> OGR or GEOG or etc).
<br>>
<br>> Best,
<br>>
<br>> P.
<br>>
<br>> On Mon, Nov 23, 2009 at 1:19 PM, Pedro Doria Meunier
<br>> <<a href="http://old.nabble.com/user/SendEmail.jtp?type=post&post=26487456&i=0" target="_top" rel="nofollow">pdoria@...</a>> wrote:
<br>>
<br>>> Hi All,
<br>>>
<br>>> First off sorry for cross-posting. I'm hoping for a bigger audience... :]
<br>>>
<br>>> As for my question:
<br>>>
<br>>> How does one deal with Postgis geom fields in C++ ?
<br>>> Google hasn't offered much help so far.
<br>>>
<br>>> I know I can make a "SELECT astext(my_geom) FROM my_table" and after parsing
<br>>> the resulting text.
<br>>> This seems a bit awkward and time-consuming...
<br>>>
<br>>> Is there a better way to deal with it?
<br>>>
<br>>> Already thankful for any ideas,
<br>>> BR,
<br>>> Pedro.
<br>>> _______________________________________________
<br>>> postgis-users mailing list
<br>>> <a href="http://old.nabble.com/user/SendEmail.jtp?type=post&post=26487456&i=1" target="_top" rel="nofollow">postgis-users@...</a>
<br>>> <a href="http://postgis.refractions.net/mailman/listinfo/postgis-users" target="_top" rel="nofollow">http://postgis.refractions.net/mailman/listinfo/postgis-users</a><br>>>
<br>>>
<br>> _______________________________________________
<br>> postgis-users mailing list
<br>> <a href="http://old.nabble.com/user/SendEmail.jtp?type=post&post=26487456&i=2" target="_top" rel="nofollow">postgis-users@...</a>
<br>> <a href="http://postgis.refractions.net/mailman/listinfo/postgis-users" target="_top" rel="nofollow">http://postgis.refractions.net/mailman/listinfo/postgis-users</a><br>>
<br>>
</div>_______________________________________________
<br>postgis-users mailing list
<br><a href="http://old.nabble.com/user/SendEmail.jtp?type=post&post=26487456&i=3" target="_top" rel="nofollow">postgis-users@...</a>
<br><a href="http://postgis.refractions.net/mailman/listinfo/postgis-users" target="_top" rel="nofollow">http://postgis.refractions.net/mailman/listinfo/postgis-users</a><br>tag:old.nabble.com,2006:post-26487416Re: Corner coordinates from Geometry2009-11-23T14:39:50Z2009-11-23T14:39:50ZSimon Greener-3Robert,
<br><br>What do you mean by "corner coordinates"?
<br><br>Do you mean the 4 corners of the MBR of each MULTIPOLYGON or do you mean that each MULTIPOLYGON is some sort of a collection of rectangular polygons and you want the corners of each of the rectangles?
<br><br>regards
<br>SImon
<br>On Tue, 24 Nov 2009 05:53:33 +1100, Robert Hicks <<a href="http://old.nabble.com/user/SendEmail.jtp?type=post&post=26487416&i=0" target="_top" rel="nofollow">robhyx@...</a>> wrote:
<br><br>> Hi all, I have a table full of MULTIPOLYGON features and I am looking
<br>> for a way to parse the geometry from a given row and get back the
<br>> corner coordinates in decimal degrees. Is there a stored procedure for
<br>> this?
<br>>
<br>> Thanks!
<br>>
<br><br><br>--
<br>SpatialDB Advice and Design, Solutions Architecture and Programming,
<br>Oracle Database 10g Administrator Certified Associate; Oracle Database 10g SQL Certified Professional
<br>Oracle Spatial, SQL Server, PostGIS, MySQL, ArcSDE, Manifold GIS, FME, Radius Topology and Studio Specialist.
<br>39 Cliff View Drive, Allens Rivulet, 7150, Tasmania, Australia.
<br>Website: www.spatialdbadvisor.com
<br> Email: <a href="http://old.nabble.com/user/SendEmail.jtp?type=post&post=26487416&i=1" target="_top" rel="nofollow">simon@...</a>
<br> Voice: +61 362 396397
<br>Mobile: +61 418 396391
<br>Skype: sggreener
<br>Longitude: 147.20515 (147° 12' 18" E)
<br>Latitude: -43.01530 (43° 00' 55" S)
<br>GeoHash: r22em9r98wg
<br>NAC:W80CK 7SWP3
<br>_______________________________________________
<br>postgis-users mailing list
<br><a href="http://old.nabble.com/user/SendEmail.jtp?type=post&post=26487416&i=2" target="_top" rel="nofollow">postgis-users@...</a>
<br><a href="http://postgis.refractions.net/mailman/listinfo/postgis-users" target="_top" rel="nofollow">http://postgis.refractions.net/mailman/listinfo/postgis-users</a><br>tag:old.nabble.com,2006:post-26486679Re: ST_* versus no ST_*2009-11-23T13:55:13Z2009-11-23T13:55:13ZPaul Ramsey-4Use the ST_* variants, in case we ever decide to delete the non-ST_ versions.
<br><br>P
<br><br>On Mon, Nov 23, 2009 at 1:27 PM, Brian Peschel <<a href="http://old.nabble.com/user/SendEmail.jtp?type=post&post=26486679&i=0" target="_top" rel="nofollow">brianp@...</a>> wrote:
<div class='shrinkable-quote'><br>> Probably an obvious thing, but.....
<br>>
<br>> Why would I use an ST_ fucntion versus a non ST_ function? For example (in
<br>> Postgis 1.3.3, I haven't had a chance to upgrade yet) there is a AsGeom()
<br>> and a ST_AsGeom() function. Both seem to return the same results. Which
<br>> one should I use?
<br>>
<br>> Thanks!
<br>> - Brian
<br>> _______________________________________________
<br>> postgis-users mailing list
<br>> <a href="http://old.nabble.com/user/SendEmail.jtp?type=post&post=26486679&i=1" target="_top" rel="nofollow">postgis-users@...</a>
<br>> <a href="http://postgis.refractions.net/mailman/listinfo/postgis-users" target="_top" rel="nofollow">http://postgis.refractions.net/mailman/listinfo/postgis-users</a><br>>
</div>_______________________________________________
<br>postgis-users mailing list
<br><a href="http://old.nabble.com/user/SendEmail.jtp?type=post&post=26486679&i=2" target="_top" rel="nofollow">postgis-users@...</a>
<br><a href="http://postgis.refractions.net/mailman/listinfo/postgis-users" target="_top" rel="nofollow">http://postgis.refractions.net/mailman/listinfo/postgis-users</a><br>tag:old.nabble.com,2006:post-26486607Re: C++ and Postgis2009-11-23T13:50:39Z2009-11-23T13:50:39ZPaul Ramsey-4That depends entirely on what you're planning to do with them once you
<br>get them. You'll probably want to serialize them into a known form so
<br>you can use them on the client side. ST_AsText() is one, very heavy,
<br>option. Mapserver, to take an example, uses
<br>encode(ST_AsBinary(ST_Force2D(geom)), 'hex'), so what goes over the
<br>wire is hex strings, but those are easy to turn into binary arrays, in
<br>WKB format, suitable for reading by a number of C++ libraries (like
<br>OGR or GEOG or etc).
<br><br>Best,
<br><br>P.
<br><br>On Mon, Nov 23, 2009 at 1:19 PM, Pedro Doria Meunier
<br><<a href="http://old.nabble.com/user/SendEmail.jtp?type=post&post=26486607&i=0" target="_top" rel="nofollow">pdoria@...</a>> wrote:
<div class='shrinkable-quote'><br>> Hi All,
<br>>
<br>> First off sorry for cross-posting. I'm hoping for a bigger audience... :]
<br>>
<br>> As for my question:
<br>>
<br>> How does one deal with Postgis geom fields in C++ ?
<br>> Google hasn't offered much help so far.
<br>>
<br>> I know I can make a "SELECT astext(my_geom) FROM my_table" and after parsing
<br>> the resulting text.
<br>> This seems a bit awkward and time-consuming...
<br>>
<br>> Is there a better way to deal with it?
<br>>
<br>> Already thankful for any ideas,
<br>> BR,
<br>> Pedro.
<br>> _______________________________________________
<br>> postgis-users mailing list
<br>> <a href="http://old.nabble.com/user/SendEmail.jtp?type=post&post=26486607&i=1" target="_top" rel="nofollow">postgis-users@...</a>
<br>> <a href="http://postgis.refractions.net/mailman/listinfo/postgis-users" target="_top" rel="nofollow">http://postgis.refractions.net/mailman/listinfo/postgis-users</a><br>>
</div>_______________________________________________
<br>postgis-users mailing list
<br><a href="http://old.nabble.com/user/SendEmail.jtp?type=post&post=26486607&i=2" target="_top" rel="nofollow">postgis-users@...</a>
<br><a href="http://postgis.refractions.net/mailman/listinfo/postgis-users" target="_top" rel="nofollow">http://postgis.refractions.net/mailman/listinfo/postgis-users</a><br>tag:old.nabble.com,2006:post-26486190ST_* versus no ST_*2009-11-23T13:27:15Z2009-11-23T13:27:15ZBrian PeschelProbably an obvious thing, but.....
<br><br>Why would I use an ST_ fucntion versus a non ST_ function? For example
<br>(in Postgis 1.3.3, I haven't had a chance to upgrade yet) there is a
<br>AsGeom() and a ST_AsGeom() function. Both seem to return the same
<br>results. Which one should I use?
<br><br>Thanks!
<br>- Brian
<br>_______________________________________________
<br>postgis-users mailing list
<br><a href="http://old.nabble.com/user/SendEmail.jtp?type=post&post=26486190&i=0" target="_top" rel="nofollow">postgis-users@...</a>
<br><a href="http://postgis.refractions.net/mailman/listinfo/postgis-users" target="_top" rel="nofollow">http://postgis.refractions.net/mailman/listinfo/postgis-users</a><br>tag:old.nabble.com,2006:post-26486182Re: C++ and Postgis2009-11-23T13:26:54Z2009-11-23T13:26:54ZHenri De Feraudy<table cellspacing="0" cellpadding="0" border="0"><tr><td valign="top" style="font: inherit;"><DIV>Well Pedro, your suggestion is what I'm doing.</DIV>
<DIV>I'm using LibPQ.</DIV>
<DIV>So I'm interested too.</DIV>
<DIV>Henri</DIV>
<DIV><BR><BR>--- On <B>Mon, 11/23/09, Pedro Doria Meunier <I><<a href="http://old.nabble.com/user/SendEmail.jtp?type=post&post=26486182&i=0" target="_top" rel="nofollow">pdoria@...</a>></I></B> wrote:<BR></DIV>
<BLOCKQUOTE style="BORDER-LEFT: rgb(16,16,255) 2px solid; PADDING-LEFT: 5px; MARGIN-LEFT: 5px"><BR>From: Pedro Doria Meunier <<a href="http://old.nabble.com/user/SendEmail.jtp?type=post&post=26486182&i=1" target="_top" rel="nofollow">pdoria@...</a>><BR>Subject: [postgis-users] C++ and Postgis<BR>To: "PostGIS Users Discussion" <<a href="http://old.nabble.com/user/SendEmail.jtp?type=post&post=26486182&i=2" target="_top" rel="nofollow">postgis-users@...</a>><BR>Cc: "Postgresql Mailing List" <<a href="http://old.nabble.com/user/SendEmail.jtp?type=post&post=26486182&i=3" target="_top" rel="nofollow">pgsql-general@...</a>><BR>Date: Monday, November 23, 2009, 10:19 PM<BR><BR>
<DIV class=plainMail>Hi All,<BR><BR>First off sorry for cross-posting. I'm hoping for a bigger audience... :]<BR><BR>As for my question:<BR><BR>How does one deal with Postgis geom fields in C++ ?<BR>Google hasn't offered much help so far.<BR><BR>I know I can make a "SELECT astext(my_geom) FROM my_table" and after parsing the resulting text.<BR>This seems a bit awkward and time-consuming...<BR><BR>Is there a better way to deal with it?<BR><BR>Already thankful for any ideas,<BR>BR,<BR>Pedro.<BR>_______________________________________________<BR>postgis-users mailing list<BR><A href="http://us.mc1108.mail.yahoo.com/mc/compose?to=postgis-users@postgis.refractions.net" ymailto="mailto:postgis-users@postgis.refractions.net" target="_top" rel="nofollow">postgis-users@...</A><BR><A href="http://postgis.refractions.net/mailman/listinfo/postgis-users" target=_blank rel="nofollow">http://postgis.refractions.net/mailman/listinfo/postgis-users</A><BR></DIV></BLOCKQUOTE></td></tr></table><br />_______________________________________________
<br>postgis-users mailing list
<br><a href="http://old.nabble.com/user/SendEmail.jtp?type=post&post=26486182&i=4" target="_top" rel="nofollow">postgis-users@...</a>
<br><a href="http://postgis.refractions.net/mailman/listinfo/postgis-users" target="_top" rel="nofollow">http://postgis.refractions.net/mailman/listinfo/postgis-users</a><br>tag:old.nabble.com,2006:post-26486109Using a rectangular selection2009-11-23T13:23:02Z2009-11-23T13:23:02ZHenri De Feraudy<table cellspacing="0" cellpadding="0" border="0"><tr><td valign="top" style="font: inherit;"><DIV>Hi,</DIV>
<DIV> I'm working in planar coordinates. I have a Postgis database and </DIV>
<DIV>for starters I want to limit my attention to all geometries that intersect a given rectangle (the "selection").</DIV>
<DIV>I can imagine that one could enumerate all geometries and ask in each case if they </DIV>
<DIV>intersect the rectangle, but that doesnt seem very efficient, so what do you suggest?</DIV>
<DIV>I only need general guidance so I can look at the various parts of the doc.</DIV>
<DIV>Later I'll be wanting to clip the geometies with this selection.</DIV>
<DIV>Still later I'll want to consider polygonal selections.</DIV>
<DIV>By the way I'm accessing the PostGis database via LibPQ and have written my application </DIV>
<DIV>in Qt.</DIV>
<DIV> </DIV>
<DIV>Thanking you in advance</DIV>
<DIV>Henri</DIV>
<DIV> </DIV></td></tr></table><br />_______________________________________________
<br>postgis-users mailing list
<br><a href="http://old.nabble.com/user/SendEmail.jtp?type=post&post=26486109&i=0" target="_top" rel="nofollow">postgis-users@...</a>
<br><a href="http://postgis.refractions.net/mailman/listinfo/postgis-users" target="_top" rel="nofollow">http://postgis.refractions.net/mailman/listinfo/postgis-users</a><br>tag:old.nabble.com,2006:post-26486046C++ and Postgis2009-11-23T13:19:44Z2009-11-23T13:19:44ZPedro Doria MeunierHi All,
<br><br>First off sorry for cross-posting. I'm hoping for a bigger audience... :]
<br><br>As for my question:
<br><br>How does one deal with Postgis geom fields in C++ ?
<br>Google hasn't offered much help so far.
<br><br>I know I can make a "SELECT astext(my_geom) FROM my_table" and after
<br>parsing the resulting text.
<br>This seems a bit awkward and time-consuming...
<br><br>Is there a better way to deal with it?
<br><br>Already thankful for any ideas,
<br>BR,
<br>Pedro.
<br>_______________________________________________
<br>postgis-users mailing list
<br><a href="http://old.nabble.com/user/SendEmail.jtp?type=post&post=26486046&i=0" target="_top" rel="nofollow">postgis-users@...</a>
<br><a href="http://postgis.refractions.net/mailman/listinfo/postgis-users" target="_top" rel="nofollow">http://postgis.refractions.net/mailman/listinfo/postgis-users</a><br>tag:old.nabble.com,2006:post-26484382Re: identifying blank geometry fields2009-11-23T11:30:13Z2009-11-23T11:30:13ZStephen WoodbridgeDid you try:
<br><br>select count(*) from yourtable where the_geom is null;
<br><br>delete from yourtable where the_geom is null:
<br><br><br>-Steve
<br><br>Bob Hume wrote:
<div class='shrinkable-quote'><br>>
<br>> I loaded a shapefile into my Postgresql/Postgis table. I'm
<br>> trying to put data in some records that appear to have no data and
<br>> eventually delete records with no geometry data. One problem is that I
<br>> can't seem to identify records with no geometry data. If I review the
<br>> the_geom field for all records, there are several records that seem to
<br>> be blank. Nothing is shown in the field when I do a listing of all
<br>> fields and all records. However, using ST_Npoints shows, for example,
<br>> 438 points for the “blank” record and St_IsValidReason says that it is a
<br>> Valid Geometry. St_GeometryType says it is ST_MultiString. However,
<br>> St_AsEWKT for the record is blank. It doesn't show any coordinates like
<br>> the other records that do have data in the the_geom field. St_IsEmpty
<br>> shows "f" for all records, including the ones that are apparently empty.
<br>>
<br>> So, it seems that the geometry data is missing, but it isn't.
<br>>
<br>> How can identify the records with apparently missing geometry
<br>> data? Thanks for the advice.
<br>>
<br>>
<br>>
<br>>
<br>> ------------------------------------------------------------------------
<br>>
<br>> _______________________________________________
<br>> postgis-users mailing list
<br>> <a href="http://old.nabble.com/user/SendEmail.jtp?type=post&post=26484382&i=0" target="_top" rel="nofollow">postgis-users@...</a>
<br>> <a href="http://postgis.refractions.net/mailman/listinfo/postgis-users" target="_top" rel="nofollow">http://postgis.refractions.net/mailman/listinfo/postgis-users</a></div><br>_______________________________________________
<br>postgis-users mailing list
<br><a href="http://old.nabble.com/user/SendEmail.jtp?type=post&post=26484382&i=1" target="_top" rel="nofollow">postgis-users@...</a>
<br><a href="http://postgis.refractions.net/mailman/listinfo/postgis-users" target="_top" rel="nofollow">http://postgis.refractions.net/mailman/listinfo/postgis-users</a><br>tag:old.nabble.com,2006:post-26483855Corner coordinates from Geometry2009-11-23T10:53:33Z2009-11-23T10:53:33ZdnbphysicistHi all, I have a table full of MULTIPOLYGON features and I am looking
<br>for a way to parse the geometry from a given row and get back the
<br>corner coordinates in decimal degrees. Is there a stored procedure for
<br>this?
<br><br>Thanks!
<br><br>--
<br>web <a href="http://www.hyxspace.com" target="_top" rel="nofollow">http://www.hyxspace.com</a><br>aim hyx1138
<br>_______________________________________________
<br>postgis-users mailing list
<br><a href="http://old.nabble.com/user/SendEmail.jtp?type=post&post=26483855&i=0" target="_top" rel="nofollow">postgis-users@...</a>
<br><a href="http://postgis.refractions.net/mailman/listinfo/postgis-users" target="_top" rel="nofollow">http://postgis.refractions.net/mailman/listinfo/postgis-users</a><br>tag:old.nabble.com,2006:post-26484271identifying blank geometry fields2009-11-23T10:21:31Z2009-11-23T10:21:31ZBob Hume<!DOCTYPE html PUBLIC "-//W3C//DTD XHTML 1.0 Transitional//EN" "http://www.w3.org/TR/xhtml1/DTD/xhtml1-transitional.dtd">
<html xmlns="http://www.w3.org/1999/xhtml" xml:lang="en" lang="en"><head>
<title></title>
<meta http-equiv="content-type" content="text/html;charset=utf-8" />
<meta http-equiv="Content-Style-Type" content="text/css" />
</head>
<body>
<div align="left"><font face="Arial" size="2"><span style=" font-size:10pt">           I loaded a shapefile into my Postgresql/Postgis table. I'm trying to put data in some
records that appear to have no data and eventually delete records with no geometry data.
One problem is that I can't seem to identify records with no geometry data.  If I review the
the_geom field for all records, there are several records that seem to be blank.  Nothing is
shown in the field when I do a listing of all fields and all records. However, using ST_Npoints
shows, for example, 438 points for the “blank” record and St_IsValidReason says that it is a
Valid Geometry.  St_GeometryType says it is ST_MultiString.  However, St_AsEWKT for the
record is blank. It doesn't show  any coordinates like the other records that do have data in
the the_geom field.  St_IsEmpty shows "f" for all records, including the ones that are
apparently empty.</span></font></div>
<div align="left"><font face="Arial" size="2"><span style=" font-size:10pt"><br />
</span></font></div>
<div align="left"><font face="Arial" size="2"><span style=" font-size:10pt">           So, it seems that the geometry data is  missing, but it isn't.</span></font></div>
<div align="left"><font face="Arial" size="2"><span style=" font-size:10pt"><br />
</span></font></div>
<div align="left"><font face="Arial" size="2"><span style=" font-size:10pt">           How can identify the records with apparently missing geometry data?  Thanks for the
advice.</span></font></div>
<div align="left"><font face="Arial" size="2"><span style=" font-size:10pt"><br />
</span></font></div>
<div align="left"> </div>
</body>
</html>
<br />_______________________________________________
<br>postgis-users mailing list
<br><a href="http://old.nabble.com/user/SendEmail.jtp?type=post&post=26484271&i=0" target="_top" rel="nofollow">postgis-users@...</a>
<br><a href="http://postgis.refractions.net/mailman/listinfo/postgis-users" target="_top" rel="nofollow">http://postgis.refractions.net/mailman/listinfo/postgis-users</a><br>tag:old.nabble.com,2006:post-26482916Re: Small benchmark with PostGIS and MS SQL-ServerSpatial: Adjustment2009-11-23T09:52:33Z2009-11-23T09:52:33ZSufficool, Stanley-2#2 will always be faster on SQL Server because it is not *entirely* ACID
<br>compliant. It has a cache rather than a full table scan for row counts
<br>IIRC.
<br><div class='shrinkable-quote'><br>> -----Original Message-----
<br>> From: <a href="http://old.nabble.com/user/SendEmail.jtp?type=post&post=26482916&i=0" target="_top" rel="nofollow">postgis-users-bounces@...</a>
<br>> [mailto:<a href="http://old.nabble.com/user/SendEmail.jtp?type=post&post=26482916&i=1" target="_top" rel="nofollow">postgis-users-bounces@...</a>] On
<br>> Behalf Of <a href="http://old.nabble.com/user/SendEmail.jtp?type=post&post=26482916&i=2" target="_top" rel="nofollow">rkrummen@...</a>
<br>> Sent: Monday, November 23, 2009 9:46 AM
<br>> To: <a href="http://old.nabble.com/user/SendEmail.jtp?type=post&post=26482916&i=3" target="_top" rel="nofollow">postgis-users@...</a>
<br>> Subject: [postgis-users] Small benchmark with PostGIS and MS
<br>> SQL-ServerSpatial: Adjustment
<br>>
<br>>
<br>> Two weeks ago, my advisor Stefan Keller posted some queries
<br>> about database benchmarking. We would like to thank you for
<br>> the useful suggestions. Wherever possible we adapted our
<br>> queries. Due to restrictions in time, data and the limited
<br>> spatial functionality of MS SQL Server, we could not
<br>> implement any proposed refinement, but we like to show you
<br>> what we got now. Feel free to comment.
<br>>
<br>> - Roland
<br>>
<br>> 1. Loading data and generating indices
<br>>
<br>> SELECT * INTO {dataset} FROM {original dataset}
<br>> WHERE ST_Intersects(@geom, the_geom) ;
<br>>
<br>> ALTER TABLE {dataset} ADD PRIMARY KEY (gid);
<br>>
<br>> CREATE INDEX gidx ON {dataset} USING gist (the_geom);
<br>>
<br>> 2. Non-spatial selection query
<br>>
<br>> SELECT Count(*) FROM {dataset lines} l
<br>> WHERE l.roadflg='Y';
<br>>
<br>> 3. Spatial query I:
<br>>
<br>> a)
<br>> SELECT Count(*) FROM {dataset points} p
<br>> WHERE ST_Intersects(@poly, p.the_geom) ;
<br>>
<br>> b)
<br>> SELECT Count(*) FROM {dataset polygons} pg
<br>> WHERE ST_Intersects(@geom, pg.the_geom) ;
<br>>
<br>> c)
<br>> SELECT Count(*) FROM {dataset lines} l
<br>> WHERE ST_Intersects(@geom, l.the_geom) ;
<br>>
<br>> 4. Spatial query II:
<br>>
<br>> a)
<br>> SELECT Count(*) FROM {dataset points} p
<br>> WHERE ST_DWITHIN(@point, p.the_geom, x)
<br>>
<br>> b)
<br>> SELECT Count(*) FROM {dataset polygons} pg
<br>> WHERE ST_DWITHIN(@point, pg.the_geom, x)
<br>>
<br>> c)
<br>> SELECT Count(*) FROM {dataset lines} l
<br>> WHERE ST_DWITHIN(@point, l.the_geom, x)
<br>>
<br>> 5. Join of a linestring and a polygon table
<br>>
<br>> SELECT SUM(ST_Length(pg.the_geom)) FROM {dataset lines} l
<br>> JOIN {dataset polygons} pg ON ST_DWITHIN(pg.the_geom,
<br>> l.the_geom, 10) WHERE l.railflg = 'Y'
<br>> _______________________________________________
<br>> postgis-users mailing list <a href="http://old.nabble.com/user/SendEmail.jtp?type=post&post=26482916&i=4" target="_top" rel="nofollow">postgis-users@...</a>
<br>> <a href="http://postgis.refractions.net/mailman/listinfo/postgis-users" target="_top" rel="nofollow">http://postgis.refractions.net/mailman/listinfo/postgis-users</a><br>>
</div>_______________________________________________
<br>postgis-users mailing list
<br><a href="http://old.nabble.com/user/SendEmail.jtp?type=post&post=26482916&i=5" target="_top" rel="nofollow">postgis-users@...</a>
<br><a href="http://postgis.refractions.net/mailman/listinfo/postgis-users" target="_top" rel="nofollow">http://postgis.refractions.net/mailman/listinfo/postgis-users</a><br>tag:old.nabble.com,2006:post-26482805Small benchmark with PostGIS and MS SQL-Server Spatial: Adjustment2009-11-23T09:45:36Z2009-11-23T09:45:36ZrkrummenTwo weeks ago, my advisor Stefan Keller posted some queries about database benchmarking. We would like to thank you for the useful suggestions. Wherever possible we adapted our queries. Due to restrictions in time, data and the limited spatial functionality of MS SQL Server, we could not implement any proposed refinement, but we like to show you what we got now. Feel free to comment.
<br><br>- Roland
<br><br>1. Loading data and generating indices
<br><br>SELECT * INTO {dataset} FROM {original dataset}
<br>WHERE ST_Intersects(@geom, the_geom) ;
<br><br>ALTER TABLE {dataset} ADD PRIMARY KEY (gid);
<br><br>CREATE INDEX gidx ON {dataset} USING gist (the_geom);
<br><br>2. Non-spatial selection query
<br><br>SELECT Count(*) FROM {dataset lines} l
<br>WHERE l.roadflg='Y';
<br><br>3. Spatial query I:
<br><br>a)
<br>SELECT Count(*) FROM {dataset points} p
<br>WHERE ST_Intersects(@poly, p.the_geom) ;
<br><br>b)
<br>SELECT Count(*) FROM {dataset polygons} pg
<br>WHERE ST_Intersects(@geom, pg.the_geom) ;
<br><br>c)
<br>SELECT Count(*) FROM {dataset lines} l
<br>WHERE ST_Intersects(@geom, l.the_geom) ;
<br><br>4. Spatial query II:
<br><br>a)
<br>SELECT Count(*) FROM {dataset points} p
<br>WHERE ST_DWITHIN(@point, p.the_geom, x)
<br><br>b)
<br>SELECT Count(*) FROM {dataset polygons} pg
<br>WHERE ST_DWITHIN(@point, pg.the_geom, x)
<br><br>c)
<br>SELECT Count(*) FROM {dataset lines} l
<br>WHERE ST_DWITHIN(@point, l.the_geom, x)
<br><br>5. Join of a linestring and a polygon table
<br><br>SELECT SUM(ST_Length(pg.the_geom)) FROM {dataset lines} l
<br>JOIN {dataset polygons} pg ON ST_DWITHIN(pg.the_geom, l.the_geom, 10)
<br>WHERE l.railflg = 'Y'
<br>_______________________________________________
<br>postgis-users mailing list
<br><a href="http://old.nabble.com/user/SendEmail.jtp?type=post&post=26482805&i=0" target="_top" rel="nofollow">postgis-users@...</a>
<br><a href="http://postgis.refractions.net/mailman/listinfo/postgis-users" target="_top" rel="nofollow">http://postgis.refractions.net/mailman/listinfo/postgis-users</a><br>tag:old.nabble.com,2006:post-26481356Re: [Fossgis-talk] problems using epsg:37852009-11-23T08:27:50Z2009-11-23T08:27:50ZAhmed DassoukiHopefully my bad day is over now that i spilled yogurt on my lap lool. i
<br>laughed. So how's the pomegranate ? how's your day going so far?
<br><br>Jan Hartmann wrote:
<div class='shrinkable-quote'><br>> Hi Peter,
<br>>
<br>> The two projection definitions you give for Google Mercator differ
<br>> only in the datum shift:
<br>>
<br>> +nadgrids=@null (epsg and PostGIS)
<br>> +towgs84=0,0,0,0,0,0,0 (spatialreference.org)
<br>>
<br>> Number one works OK for me, number two is wrong. You can find out
<br>> about nadgrids at:
<br>>
<br>> <a href="http://trac.osgeo.org/proj/wiki/FAQ" target="_top" rel="nofollow">http://trac.osgeo.org/proj/wiki/FAQ</a>, section "Changing Ellipsoid / Why
<br>> can't I convert from WGS84 to Virtual Globe Mercator?"
<br>>
<br>> As to the two gdalwarp's: *both* projection definitions need a
<br>> +datum, +towgs84 or a +nadgrids parameter, else the
<br>> WGS84-datum-transformation won't work. In your case, the first warp
<br>> looks OK to me (both projections are OK), the second one gets its
<br>> source projection from the source file. I would guess that something
<br>> is wrong with the datum information there, and that the first warp is
<br>> the correct one. Don't you have any ground information?
<br>>
<br>> Jan
<br>>
<br>> Peter Freimuth wrote:
<br>>> Hi all,
<br>>>
<br>>> i am struggling around with a projection problem. I want to harmonize
<br>>> several hundreds of satellite images from different UTM Zones to the
<br>>> Google Projection Spherical Mercator. I am using a python script which
<br>>> handles the job for me using "gdalwarp" from the 1.6.1 and 1.7.0dev
<br>>> release (same result).
<br>>>
<br>>> Find attached the gdalinfo output of the original file
<br>>> (gdalinfo_2009-03-02T140133_RE3_3A-NAC_642351_34785_enhancedRGB.txt):
<br>>>
<br>>> Upper Left ( 691500.000,-2471500.000) ( 49d 8'25.94"W, 22d20'19.37"S)
<br>>> Lower Left ( 691500.000,-2496500.000) ( 49d 8'15.09"W, 22d33'52.01"S)
<br>>> Upper Right ( 716500.000,-2471500.000) ( 48d53'52.47"W, 22d20'8.69"S)
<br>>> Lower Right ( 716500.000,-2496500.000) ( 48d53'40.21"W, 22d33'41.21"S)
<br>>>
<br>>> echo "warp giving a source epsg code"
<br>>> gdalwarp -s_srs "EPSG:32622" -t_srs "EPSG:3785"
<br>>> 2009-03-02T140133_RE3_3A-NAC_642351_34785_enhancedRGB.pix
<br>>> c:\tmp\2009-03-02T140133_RE3_3A-NAC_642351_34785_enhancedRGB_32622-pix-3785.tif
<br>>> this results in
<br>>> (gdalinfo_2009-03-02T140133_RE3_3A-NAC_642351_34785_32622-pix-3785.txt)
<br>>>
<br>>> Upper Left (-5470299.856,-2535649.342)
<br>>> Lower Left (-5470299.856,-2563039.918)
<br>>> Upper Right (-5442909.280,-2535649.342)
<br>>> Lower Right (-5442909.280,-2563039.918)
<br>>>
<br>>> echo "warp without giving a source epsg code"
<br>>> gdalwarp -t_srs "EPSG:3785"
<br>>> 2009-03-02T140133_RE3_3A-NAC_642351_34785_enhancedRGB.pix
<br>>> c:\tmp\2009-03-02T140133_RE3_3A-NAC_642351_34785_enhancedRGB_pix-3785.tif
<br>>>
<br>>> This results in
<br>>> (gdalinfo_2009-03-02T140133_RE3_3A-NAC_642351_34785_pix-3785.txt)
<br>>>
<br>>> Upper Left (-5470299.856,-2551883.715)
<br>>> Lower Left (-5470299.856,-2579432.525)
<br>>> Upper Right (-5442913.703,-2551883.715)
<br>>> Lower Right (-5442913.703,-2579432.525)
<br>>>
<br>>> Somehow i get big shift in the Y values of the images.
<br>>>
<br>>> As a kind of reference i have a polygon layer in epsg:4326 stored in a
<br>>> postgis 1.4.0 database table which contains the exact definitions of the
<br>>> tiling grid which is used to create the satallite images.
<br>>>
<br>>> The same shift can be reproduce in postgis when using a somehow
<br>>> different epsg:3785 definition:
<br>>> As you can see, there is a significant difference in the Y value.
<br>>>
<br>>> select
<br>>> astext(envelope(transform(GeomFromEWKT('SRID=4326;POLYGON((-49.1405410766602
<br>>> -22.5644474029541,-49.1405410766602 -22.3357467651367,-48.8945007324219
<br>>> -22.3357467651367,-48.8945007324219 -22.5644474029541,-49.1405410766602
<br>>> -22.5644474029541))'::TEXT),3785)))
<br>>> =>>"POLYGON((-5470300 -2579430.25,-5470300 -2551883.5,-5442911
<br>>> -2551883.5,-5442911 -2579430.25,-5470300 -2579430.25))"
<br>>>
<br>>> select
<br>>> astext(envelope(transform(GeomFromEWKT('SRID=4326;POLYGON((-49.1405410766602
<br>>> -22.5644474029541,-49.1405410766602 -22.3357467651367,-48.8945007324219
<br>>> -22.3357467651367,-48.8945007324219 -22.5644474029541,-49.1405410766602
<br>>> -22.5644474029541))'::TEXT),903785)))
<br>>> =>>"POLYGON((-5470300 -2563038,-5470300 -2535649.25,-5442911
<br>>> -2535649.25,-5442911 -2563038,-5470300 -2563038))"
<br>>>
<br>>> Here are the two definitons as they are in my PostGIS spatial_ref_sys table.
<br>>>
<br>>> "3785" as defined in postgis 1.4.0 and proj-4.7.0 epsg file
<br>>> +proj=merc +a=6378137 +b=6378137 +lat_ts=0.0 +lon_0=0.0 +x_0=0.0 +y_0=0
<br>>> +units=m +k=1.0 +nadgrids=@null +no_defs
<br>>> PROJCS["Popular Visualisation CRS / Mercator (deprecated)",
<br>>> GEOGCS["Popular Visualisation CRS",
<br>>> DATUM["Popular_Visualisation_Datum",
<br>>> SPHEROID["Popular Visualisation Sphere",6378137,0,
<br>>> AUTHORITY["EPSG","7059"]],
<br>>> TOWGS84[0,0,0,0,0,0,0],
<br>>> AUTHORITY["EPSG","6055"]],
<br>>> PRIMEM["Greenwich",0,
<br>>> AUTHORITY["EPSG","8901"]],
<br>>> UNIT["degree",0.01745329251994328,
<br>>> AUTHORITY["EPSG","9122"]],
<br>>> AUTHORITY["EPSG","4055"]],
<br>>> UNIT["metre",1,
<br>>> AUTHORITY["EPSG","9001"]],
<br>>> PROJECTION["Mercator_1SP"],
<br>>> PARAMETER["central_meridian",0],
<br>>> PARAMETER["scale_factor",1],
<br>>> PARAMETER["false_easting",0],
<br>>> PARAMETER["false_northing",0],
<br>>> AUTHORITY["EPSG","3785"],
<br>>> AXIS["X",EAST],
<br>>> AXIS["Y",NORTH]]"
<br>>>
<br>>> I added a second definition which i took from spatialreference.org
<br>>> (<a href="http://spatialreference.org/ref/epsg/3785" target="_top" rel="nofollow">http://spatialreference.org/ref/epsg/3785</a>) with the srid 903785
<br>>> "903785"
<br>>> +proj=merc +lon_0=0 +k=1 +x_0=0 +y_0=0 +a=6378137 +b=6378137
<br>>> +towgs84=0,0,0,0,0,0,0 +units=m +no_defs
<br>>> PROJCS["Popular Visualisation CRS / Mercator",
<br>>> GEOGCS["Popular Visualisation CRS",
<br>>> DATUM["Popular_Visualisation_Datum",
<br>>> SPHEROID["Popular Visualisation Sphere",6378137,0,
<br>>> AUTHORITY["EPSG","7059"]],
<br>>> TOWGS84[0,0,0,0,0,0,0],
<br>>> AUTHORITY["EPSG","6055"]],
<br>>> PRIMEM["Greenwich",0,
<br>>> AUTHORITY["EPSG","8901"]],
<br>>> UNIT["degree",0.01745329251994328,
<br>>> AUTHORITY["EPSG","9122"]],
<br>>> AUTHORITY["EPSG","4055"]],
<br>>> UNIT["metre",1,
<br>>> AUTHORITY["EPSG","9001"]],
<br>>> PROJECTION["Mercator_1SP"],
<br>>> PARAMETER["central_meridian",0],
<br>>> PARAMETER["scale_factor",1],
<br>>> PARAMETER["false_easting",0],
<br>>> PARAMETER["false_northing",0],
<br>>> AUTHORITY["EPSG","3785"],
<br>>> AXIS["X",EAST],
<br>>> AXIS["Y",NORTH]]
<br>>>
<br>>>
<br>>> Now the big questions:
<br>>> Which of the above given definitions is the correct one?? So which of
<br>>> the polygons is the correct result when i want epsg:3785?
<br>>>
<br>>> How can it be explained that gdalwarp seems to use two different
<br>>> projections when warping the same file to the same output projection?
<br>>>
<br>>> Any ideas or comments??
<br>>>
<br>>> Thanks for any kind of help,
<br>>> Peter
<br>>>
<br>>>
<br>>> ------------------------------------------------------------------------
<br>>>
<br>>> _______________________________________________
<br>>> Fossgis-talk-liste mailing list
<br>>> <a href="http://old.nabble.com/user/SendEmail.jtp?type=post&post=26481356&i=0" target="_top" rel="nofollow">Fossgis-talk-liste@...</a>
<br>>> <a href="https://lists.fossgis.de/mailman/listinfo/fossgis-talk-liste" target="_top" rel="nofollow">https://lists.fossgis.de/mailman/listinfo/fossgis-talk-liste</a><br>> ------------------------------------------------------------------------
<br>>
<br>> _______________________________________________
<br>> postgis-users mailing list
<br>> <a href="http://old.nabble.com/user/SendEmail.jtp?type=post&post=26481356&i=1" target="_top" rel="nofollow">postgis-users@...</a>
<br>> <a href="http://postgis.refractions.net/mailman/listinfo/postgis-users" target="_top" rel="nofollow">http://postgis.refractions.net/mailman/listinfo/postgis-users</a><br>>
</div>_______________________________________________
<br>postgis-users mailing list
<br><a href="http://old.nabble.com/user/SendEmail.jtp?type=post&post=26481356&i=2" target="_top" rel="nofollow">postgis-users@...</a>
<br><a href="http://postgis.refractions.net/mailman/listinfo/postgis-users" target="_top" rel="nofollow">http://postgis.refractions.net/mailman/listinfo/postgis-users</a><br>tag:old.nabble.com,2006:post-26481105Re: [Fossgis-talk] problems using epsg:37852009-11-23T08:16:52Z2009-11-23T08:16:52ZJan Hartmann<!DOCTYPE html PUBLIC "-//W3C//DTD HTML 4.01 Transitional//EN">
<html>
<head>
<meta content="text/html;charset=ISO-8859-15" http-equiv="Content-Type">
</head>
<body bgcolor="#ffffff" text="#000000">
Hi Peter,<br>
<br>
The two projection definitions you give for Google Mercator differ only
in the datum shift:<br>
<br>
+nadgrids=@null (epsg and PostGIS)<br>
+towgs84=0,0,0,0,0,0,0 (spatialreference.org)<br>
<br>
Number one works OK for me, number two is wrong. You can find out
about nadgrids at:<br>
<br>
<a class="moz-txt-link-freetext" href="http://trac.osgeo.org/proj/wiki/FAQ" target="_top" rel="nofollow">http://trac.osgeo.org/proj/wiki/FAQ</a>, section "Changing Ellipsoid / Why
can't I convert from WGS84 to Virtual Globe Mercator?"<br>
<br>
As to the two gdalwarp's: *both* projection definitions need a
+datum, +towgs84 or a +nadgrids parameter, else the
WGS84-datum-transformation won't work. In your case, the first warp
looks OK to me (both projections are OK), the second one gets its
source projection from the source file. I would guess that something is
wrong with the datum information there, and that the first warp is the
correct one. Don't you have any ground information?<br>
<br>
Jan<br>
<br>
Peter Freimuth wrote:
<blockquote cite="mid:4B0AA744.1070202@arcor.de" type="cite">
<pre wrap="">Hi all,
i am struggling around with a projection problem. I want to harmonize
several hundreds of satellite images from different UTM Zones to the
Google Projection Spherical Mercator. I am using a python script which
handles the job for me using "gdalwarp" from the 1.6.1 and 1.7.0dev
release (same result).
Find attached the gdalinfo output of the original file
(gdalinfo_2009-03-02T140133_RE3_3A-NAC_642351_34785_enhancedRGB.txt):
Upper Left ( 691500.000,-2471500.000) ( 49d 8'25.94"W, 22d20'19.37"S)
Lower Left ( 691500.000,-2496500.000) ( 49d 8'15.09"W, 22d33'52.01"S)
Upper Right ( 716500.000,-2471500.000) ( 48d53'52.47"W, 22d20'8.69"S)
Lower Right ( 716500.000,-2496500.000) ( 48d53'40.21"W, 22d33'41.21"S)
echo "warp giving a source epsg code"
gdalwarp -s_srs "EPSG:32622" -t_srs "EPSG:3785"
2009-03-02T140133_RE3_3A-NAC_642351_34785_enhancedRGB.pix
c:\tmp\2009-03-02T140133_RE3_3A-NAC_642351_34785_enhancedRGB_32622-pix-3785.tif
this results in
(gdalinfo_2009-03-02T140133_RE3_3A-NAC_642351_34785_32622-pix-3785.txt)
Upper Left (-5470299.856,-2535649.342)
Lower Left (-5470299.856,-2563039.918)
Upper Right (-5442909.280,-2535649.342)
Lower Right (-5442909.280,-2563039.918)
echo "warp without giving a source epsg code"
gdalwarp -t_srs "EPSG:3785"
2009-03-02T140133_RE3_3A-NAC_642351_34785_enhancedRGB.pix
c:\tmp\2009-03-02T140133_RE3_3A-NAC_642351_34785_enhancedRGB_pix-3785.tif
This results in
(gdalinfo_2009-03-02T140133_RE3_3A-NAC_642351_34785_pix-3785.txt)
Upper Left (-5470299.856,-2551883.715)
Lower Left (-5470299.856,-2579432.525)
Upper Right (-5442913.703,-2551883.715)
Lower Right (-5442913.703,-2579432.525)
Somehow i get big shift in the Y values of the images.
As a kind of reference i have a polygon layer in epsg:4326 stored in a
postgis 1.4.0 database table which contains the exact definitions of the
tiling grid which is used to create the satallite images.
The same shift can be reproduce in postgis when using a somehow
different epsg:3785 definition:
As you can see, there is a significant difference in the Y value.
select
astext(envelope(transform(GeomFromEWKT('SRID=4326;POLYGON((-49.1405410766602
-22.5644474029541,-49.1405410766602 -22.3357467651367,-48.8945007324219
-22.3357467651367,-48.8945007324219 -22.5644474029541,-49.1405410766602
-22.5644474029541))'::TEXT),3785)))
=>>"POLYGON((-5470300 -2579430.25,-5470300 -2551883.5,-5442911
-2551883.5,-5442911 -2579430.25,-5470300 -2579430.25))"
select
astext(envelope(transform(GeomFromEWKT('SRID=4326;POLYGON((-49.1405410766602
-22.5644474029541,-49.1405410766602 -22.3357467651367,-48.8945007324219
-22.3357467651367,-48.8945007324219 -22.5644474029541,-49.1405410766602
-22.5644474029541))'::TEXT),903785)))
=>>"POLYGON((-5470300 -2563038,-5470300 -2535649.25,-5442911
-2535649.25,-5442911 -2563038,-5470300 -2563038))"
Here are the two definitons as they are in my PostGIS spatial_ref_sys table.
"3785" as defined in postgis 1.4.0 and proj-4.7.0 epsg file
+proj=merc +a=6378137 +b=6378137 +lat_ts=0.0 +lon_0=0.0 +x_0=0.0 +y_0=0
+units=m +k=1.0 +nadgrids=@null +no_defs
PROJCS["Popular Visualisation CRS / Mercator (deprecated)",
GEOGCS["Popular Visualisation CRS",
DATUM["Popular_Visualisation_Datum",
SPHEROID["Popular Visualisation Sphere",6378137,0,
AUTHORITY["EPSG","7059"]],
TOWGS84[0,0,0,0,0,0,0],
AUTHORITY["EPSG","6055"]],
PRIMEM["Greenwich",0,
AUTHORITY["EPSG","8901"]],
UNIT["degree",0.01745329251994328,
AUTHORITY["EPSG","9122"]],
AUTHORITY["EPSG","4055"]],
UNIT["metre",1,
AUTHORITY["EPSG","9001"]],
PROJECTION["Mercator_1SP"],
PARAMETER["central_meridian",0],
PARAMETER["scale_factor",1],
PARAMETER["false_easting",0],
PARAMETER["false_northing",0],
AUTHORITY["EPSG","3785"],
AXIS["X",EAST],
AXIS["Y",NORTH]]"
I added a second definition which i took from spatialreference.org
(<a class="moz-txt-link-freetext" href="http://spatialreference.org/ref/epsg/3785" target="_top" rel="nofollow">http://spatialreference.org/ref/epsg/3785</a>) with the srid 903785
"903785"
+proj=merc +lon_0=0 +k=1 +x_0=0 +y_0=0 +a=6378137 +b=6378137
+towgs84=0,0,0,0,0,0,0 +units=m +no_defs
PROJCS["Popular Visualisation CRS / Mercator",
GEOGCS["Popular Visualisation CRS",
DATUM["Popular_Visualisation_Datum",
SPHEROID["Popular Visualisation Sphere",6378137,0,
AUTHORITY["EPSG","7059"]],
TOWGS84[0,0,0,0,0,0,0],
AUTHORITY["EPSG","6055"]],
PRIMEM["Greenwich",0,
AUTHORITY["EPSG","8901"]],
UNIT["degree",0.01745329251994328,
AUTHORITY["EPSG","9122"]],
AUTHORITY["EPSG","4055"]],
UNIT["metre",1,
AUTHORITY["EPSG","9001"]],
PROJECTION["Mercator_1SP"],
PARAMETER["central_meridian",0],
PARAMETER["scale_factor",1],
PARAMETER["false_easting",0],
PARAMETER["false_northing",0],
AUTHORITY["EPSG","3785"],
AXIS["X",EAST],
AXIS["Y",NORTH]]
Now the big questions:
Which of the above given definitions is the correct one?? So which of
the polygons is the correct result when i want epsg:3785?
How can it be explained that gdalwarp seems to use two different
projections when warping the same file to the same output projection?
Any ideas or comments??
Thanks for any kind of help,
Peter
</pre>
<pre wrap=""><pre wrap=""><pre wrap=""><pre wrap="">
<hr size="4" width="90%">
_______________________________________________
Fossgis-talk-liste mailing list
<a href="http://old.nabble.com/user/SendEmail.jtp?type=post&post=26481105&i=0" target="_top" rel="nofollow">Fossgis-talk-liste@...</a>
<a class="moz-txt-link-freetext" href="https://lists.fossgis.de/mailman/listinfo/fossgis-talk-liste" target="_top" rel="nofollow">https://lists.fossgis.de/mailman/listinfo/fossgis-talk-liste</a>
</pre></pre></pre></pre>
</blockquote>
</body>
</html>
<br />_______________________________________________
<br>postgis-users mailing list
<br><a href="http://old.nabble.com/user/SendEmail.jtp?type=post&post=26481105&i=1" target="_top" rel="nofollow">postgis-users@...</a>
<br><a href="http://postgis.refractions.net/mailman/listinfo/postgis-users" target="_top" rel="nofollow">http://postgis.refractions.net/mailman/listinfo/postgis-users</a><br>tag:old.nabble.com,2006:post-26479900problems using epsg:37852009-11-23T07:16:20Z2009-11-23T07:16:20ZPeter FreimuthHi all,
<br><br>i am struggling around with a projection problem. I want to harmonize
<br>several hundreds of satellite images from different UTM Zones to the
<br>Google Projection Spherical Mercator. I am using a python script which
<br>handles the job for me using "gdalwarp" from the 1.6.1 and 1.7.0dev
<br>release (same result).
<br>
<br>Find attached the gdalinfo output of the original file
<br>(gdalinfo_2009-03-02T140133_RE3_3A-NAC_642351_34785_enhancedRGB.txt):
<br><br>Upper Left ( 691500.000,-2471500.000) ( 49d 8'25.94"W, 22d20'19.37"S)
<br>Lower Left ( 691500.000,-2496500.000) ( 49d 8'15.09"W, 22d33'52.01"S)
<br>Upper Right ( 716500.000,-2471500.000) ( 48d53'52.47"W, 22d20'8.69"S)
<br>Lower Right ( 716500.000,-2496500.000) ( 48d53'40.21"W, 22d33'41.21"S)
<br><br>echo "warp giving a source epsg code"
<br>gdalwarp -s_srs "EPSG:32622" -t_srs "EPSG:3785"
<br>2009-03-02T140133_RE3_3A-NAC_642351_34785_enhancedRGB.pix
<br>c:\tmp\2009-03-02T140133_RE3_3A-NAC_642351_34785_enhancedRGB_32622-pix-3785.tif
<br>this results in
<br>(gdalinfo_2009-03-02T140133_RE3_3A-NAC_642351_34785_32622-pix-3785.txt)
<br><br>Upper Left (-5470299.856,-2535649.342)
<br>Lower Left (-5470299.856,-2563039.918)
<br>Upper Right (-5442909.280,-2535649.342)
<br>Lower Right (-5442909.280,-2563039.918)
<br><br>echo "warp without giving a source epsg code"
<br>gdalwarp -t_srs "EPSG:3785"
<br>2009-03-02T140133_RE3_3A-NAC_642351_34785_enhancedRGB.pix
<br>c:\tmp\2009-03-02T140133_RE3_3A-NAC_642351_34785_enhancedRGB_pix-3785.tif
<br><br>This results in
<br>(gdalinfo_2009-03-02T140133_RE3_3A-NAC_642351_34785_pix-3785.txt)
<br><br>Upper Left (-5470299.856,-2551883.715)
<br>Lower Left (-5470299.856,-2579432.525)
<br>Upper Right (-5442913.703,-2551883.715)
<br>Lower Right (-5442913.703,-2579432.525)
<br><br>Somehow i get big shift in the Y values of the images.
<br><br>As a kind of reference i have a polygon layer in epsg:4326 stored in a
<br>postgis 1.4.0 database table which contains the exact definitions of the
<br>tiling grid which is used to create the satallite images.
<br><br>The same shift can be reproduce in postgis when using a somehow
<br>different epsg:3785 definition:
<br>As you can see, there is a significant difference in the Y value.
<br><br>select
<br>astext(envelope(transform(GeomFromEWKT('SRID=4326;POLYGON((-49.1405410766602
<br>-22.5644474029541,-49.1405410766602 -22.3357467651367,-48.8945007324219
<br>-22.3357467651367,-48.8945007324219 -22.5644474029541,-49.1405410766602
<br>-22.5644474029541))'::TEXT),3785)))
<br>=>>"POLYGON((-5470300 -2579430.25,-5470300 -2551883.5,-5442911
<br>-2551883.5,-5442911 -2579430.25,-5470300 -2579430.25))"
<br><br>select
<br>astext(envelope(transform(GeomFromEWKT('SRID=4326;POLYGON((-49.1405410766602
<br>-22.5644474029541,-49.1405410766602 -22.3357467651367,-48.8945007324219
<br>-22.3357467651367,-48.8945007324219 -22.5644474029541,-49.1405410766602
<br>-22.5644474029541))'::TEXT),903785)))
<br>=>>"POLYGON((-5470300 -2563038,-5470300 -2535649.25,-5442911
<br>-2535649.25,-5442911 -2563038,-5470300 -2563038))"
<br><br>Here are the two definitons as they are in my PostGIS spatial_ref_sys table.
<br><br>"3785" as defined in postgis 1.4.0 and proj-4.7.0 epsg file
<br>+proj=merc +a=6378137 +b=6378137 +lat_ts=0.0 +lon_0=0.0 +x_0=0.0 +y_0=0
<br>+units=m +k=1.0 +nadgrids=@null +no_defs
<br>PROJCS["Popular Visualisation CRS / Mercator (deprecated)",
<br> GEOGCS["Popular Visualisation CRS",
<br> DATUM["Popular_Visualisation_Datum",
<br> SPHEROID["Popular Visualisation Sphere",6378137,0,
<br> AUTHORITY["EPSG","7059"]],
<br> TOWGS84[0,0,0,0,0,0,0],
<br> AUTHORITY["EPSG","6055"]],
<br> PRIMEM["Greenwich",0,
<br> AUTHORITY["EPSG","8901"]],
<br> UNIT["degree",0.01745329251994328,
<br> AUTHORITY["EPSG","9122"]],
<br> AUTHORITY["EPSG","4055"]],
<br> UNIT["metre",1,
<br> AUTHORITY["EPSG","9001"]],
<br> PROJECTION["Mercator_1SP"],
<br> PARAMETER["central_meridian",0],
<br> PARAMETER["scale_factor",1],
<br> PARAMETER["false_easting",0],
<br> PARAMETER["false_northing",0],
<br> AUTHORITY["EPSG","3785"],
<br> AXIS["X",EAST],
<br> AXIS["Y",NORTH]]"
<br>
<br>I added a second definition which i took from spatialreference.org
<br>(<a href="http://spatialreference.org/ref/epsg/3785" target="_top" rel="nofollow">http://spatialreference.org/ref/epsg/3785</a>) with the srid 903785
<br>"903785"
<br>+proj=merc +lon_0=0 +k=1 +x_0=0 +y_0=0 +a=6378137 +b=6378137
<br>+towgs84=0,0,0,0,0,0,0 +units=m +no_defs
<br>PROJCS["Popular Visualisation CRS / Mercator",
<br> GEOGCS["Popular Visualisation CRS",
<br> DATUM["Popular_Visualisation_Datum",
<br> SPHEROID["Popular Visualisation Sphere",6378137,0,
<br> AUTHORITY["EPSG","7059"]],
<br> TOWGS84[0,0,0,0,0,0,0],
<br> AUTHORITY["EPSG","6055"]],
<br> PRIMEM["Greenwich",0,
<br> AUTHORITY["EPSG","8901"]],
<br> UNIT["degree",0.01745329251994328,
<br> AUTHORITY["EPSG","9122"]],
<br> AUTHORITY["EPSG","4055"]],
<br> UNIT["metre",1,
<br> AUTHORITY["EPSG","9001"]],
<br> PROJECTION["Mercator_1SP"],
<br> PARAMETER["central_meridian",0],
<br> PARAMETER["scale_factor",1],
<br> PARAMETER["false_easting",0],
<br> PARAMETER["false_northing",0],
<br> AUTHORITY["EPSG","3785"],
<br> AXIS["X",EAST],
<br> AXIS["Y",NORTH]]
<br><br><br>Now the big questions:
<br>Which of the above given definitions is the correct one?? So which of
<br>the polygons is the correct result when i want epsg:3785?
<br><br>How can it be explained that gdalwarp seems to use two different
<br>projections when warping the same file to the same output projection?
<br><br>Any ideas or comments??
<br><br>Thanks for any kind of help,
<br>Peter
<br><br>--
<br>Peter Freimuth
<br><br><br />Driver: GTiff/GeoTIFF
<br>Files: 2009-03-02T140133_RE3_3A-NAC_642351_34785_enhancedRGB_pix-3785.tif
<br>Size is 5051, 5081
<br>Coordinate System is:
<br>PROJCS["unnamed",
<br> GEOGCS[,
<br> DATUM["unknown",
<br> SPHEROID["unretrievable - using WGS84",6378137,298.257223563]],
<br> PRIMEM["Greenwich",0],
<br> UNIT["degree",0.0174532925199433]],
<br> UNIT["metre",1,
<br> AUTHORITY["EPSG","9001"]],
<br> AUTHORITY["EPSG","3785"]]
<br>Origin = (-5470299.855895809800000,-2551883.714705614400000)
<br>Pixel Size = (5.421926883630149,-5.421926883630149)
<br>Metadata:
<br> AREA_OR_POINT=Area
<br>Image Structure Metadata:
<br> INTERLEAVE=PIXEL
<br>Corner Coordinates:
<br>Upper Left (-5470299.856,-2551883.715)
<br>Lower Left (-5470299.856,-2579432.525)
<br>Upper Right (-5442913.703,-2551883.715)
<br>Lower Right (-5442913.703,-2579432.525)
<br>Center (-5456606.780,-2565658.120)
<br>Band 1 Block=5051x1 Type=Byte, ColorInterp=Red
<br>Band 2 Block=5051x1 Type=Byte, ColorInterp=Green
<br>Band 3 Block=5051x1 Type=Byte, ColorInterp=Blue
<br><br />Driver: GTiff/GeoTIFF
<br>Files: 2009-03-02T140133_RE3_3A-NAC_642351_34785_enhancedRGB_32622-pix-3785.tif
<br>Size is 5066, 5066
<br>Coordinate System is:
<br>PROJCS["unnamed",
<br> GEOGCS[,
<br> DATUM["unknown",
<br> SPHEROID["unretrievable - using WGS84",6378137,298.257223563]],
<br> PRIMEM["Greenwich",0],
<br> UNIT["degree",0.0174532925199433]],
<br> UNIT["metre",1,
<br> AUTHORITY["EPSG","9001"]],
<br> AUTHORITY["EPSG","3785"]]
<br>Origin = (-5470299.855895809800000,-2535649.341730865200000)
<br>Pixel Size = (5.406746173490638,-5.406746173490638)
<br>Metadata:
<br> AREA_OR_POINT=Area
<br>Image Structure Metadata:
<br> INTERLEAVE=PIXEL
<br>Corner Coordinates:
<br>Upper Left (-5470299.856,-2535649.342)
<br>Lower Left (-5470299.856,-2563039.918)
<br>Upper Right (-5442909.280,-2535649.342)
<br>Lower Right (-5442909.280,-2563039.918)
<br>Center (-5456604.568,-2549344.630)
<br>Band 1 Block=5066x1 Type=Byte, ColorInterp=Red
<br>Band 2 Block=5066x1 Type=Byte, ColorInterp=Green
<br>Band 3 Block=5066x1 Type=Byte, ColorInterp=Blue
<br><br />Driver: PCIDSK/PCIDSK Database File
<br>Files: 2009-03-02T140133_RE3_3A-NAC_642351_34785_enhancedRGB.pix
<br>Size is 5000, 5000
<br>Coordinate System is:
<br>PROJCS["UTM Zone 22, Northern Hemisphere",
<br> GEOGCS["Unknown datum based upon the WGS 84 ellipsoid",
<br> DATUM["Not specified (based on WGS 84 spheroid)",
<br> SPHEROID["WGS 84",6378137,298.257223563,
<br> AUTHORITY["EPSG","7030"]]],
<br> PRIMEM["Greenwich",0],
<br> UNIT["degree",0.0174532925199433]],
<br> PROJECTION["Transverse_Mercator"],
<br> PARAMETER["latitude_of_origin",0],
<br> PARAMETER["central_meridian",-51],
<br> PARAMETER["scale_factor",0.9996],
<br> PARAMETER["false_easting",500000],
<br> PARAMETER["false_northing",0],
<br> UNIT["Meter",1]]
<br>Origin = (691500.000000000000000,-2471500.000000000000000)
<br>Pixel Size = (5.000000000000000,-5.000000000000000)
<br>Corner Coordinates:
<br>Upper Left ( 691500.000,-2471500.000) ( 49d 8'25.94"W, 22d20'19.37"S)
<br>Lower Left ( 691500.000,-2496500.000) ( 49d 8'15.09"W, 22d33'52.01"S)
<br>Upper Right ( 716500.000,-2471500.000) ( 48d53'52.47"W, 22d20'8.69"S)
<br>Lower Right ( 716500.000,-2496500.000) ( 48d53'40.21"W, 22d33'41.21"S)
<br>Center ( 704000.000,-2484000.000) ( 49d 1'3.43"W, 22d27'0.49"S)
<br>Band 1 Block=5000x1 Type=Byte, ColorInterp=Undefined
<br> Description = Contents Not Specified
<br>...
<br>Band 2 Block=5000x1 Type=Byte, ColorInterp=Undefined
<br> Description = Contents Not Specified
<br>...
<br>Band 3 Block=5000x1 Type=Byte, ColorInterp=Undefined
<br> Description = Contents Not Specified
<br>...<br />_______________________________________________
<br>postgis-users mailing list
<br><a href="http://old.nabble.com/user/SendEmail.jtp?type=post&post=26479900&i=0" target="_top" rel="nofollow">postgis-users@...</a>
<br><a href="http://postgis.refractions.net/mailman/listinfo/postgis-users" target="_top" rel="nofollow">http://postgis.refractions.net/mailman/listinfo/postgis-users</a><br>tag:old.nabble.com,2006:post-26449125Re: Server Crash on Invalid 3d Polygon2009-11-20T11:36:35Z2009-11-20T11:36:35ZPaul Ramsey-4Love to see a bug with such a succinct test case. It crashes for me
<br>too, I've ticked it and will be looking now.
<br><br><a href="http://trac.osgeo.org/postgis/ticket/316" target="_top" rel="nofollow">http://trac.osgeo.org/postgis/ticket/316</a><br><br>P
<br><br>On Fri, Nov 20, 2009 at 11:31 AM, Brad Hecox <<a href="http://old.nabble.com/user/SendEmail.jtp?type=post&post=26449125&i=0" target="_top" rel="nofollow">brad@...</a>> wrote:
<div class='shrinkable-quote'><br>> I am having a problem were I get a server crash any time I try to validate or generate a geometry on an invalid 3D polygon. I have a set of polygons in EWKT format that I need to validate and convert to geometries if valid. In most cases the st_isvalid function returns a true or false, however of the polygon is 3D and is invalid rather than returning a false, the system crashes. I have been able to reproduce this on Solaris x86, Solaris SPARC and Max OSX. We are running PostgreSQL 8.4.1 with PostGIS 1.4.0. The following is the pstack from the core dump:
<br>>
<br>> ore 'core' of 4774: /opt/postgres/pgserver/8.4.1/bin/postgres
<br>> fef68610 strncat (4d8e38, 0, 40923bff, 28, 0, 0) + 20
<br>> fedd3360 pg_parser_errhint (ffbfd9e8, 4d8e38, 7, fee02a08, fee026d0, fee34224) + 24
<br>> fedd9c74 LWGEOM_in (ffbfda70, 5167bc, d3280f32, 406000, 0, 3) + 70
<br>> 002a0ef4 InputFunctionCall (ffbfdce8, 4d8e38, 4010, ffffffff, 516688, 4010) + 58
<br>> 002a1120 OidInputFunctionCall (4012, 4d8e38, 4010, ffffffff, 0, 445000) + 20
<br>> 000e43e8 coerce_type (4d9150, 4d90b0, 2c1, 4010, ffffffff, 0) + 3d4
<br>> 000da760 make_fn_arguments (4d8d60, 4d90d8, ffbfdef0, 4d9118, 0, 0) + 5c
<br>> 000db534 ParseFuncOrColumn (4d8d60, 4d8e20, 4d90e8, 2, 0, 0) + 3d8
<br>> 000d7d00 transformExpr (0, 4d8f00, 4d9020, 31, 1518, 4da538) + 150
<br>> 000e6cb8 transformTargetEntry (4d8d60, 4d8f00, 0, 0, 0, 2b65c4) + 80
<br>> 000e6d0c transformTargetList (0, 4d8f50, 398ea0, fee62a00, 1e70, 0) + 44
<br>> 000bc6a4 transformStmt (4d9028, 4d8f78, 0, 11ab4, 529f2c, 529f2c) + e0
<br>> 000bed78 parse_analyze (4d8f78, 4d8550, 0, 0, 41e400, 529f00) + 20
<br>> 001e8664 pg_analyze_and_rewrite (4d8f78, 4d8550, 0, 0, 4627d0, 462748) + 24
<br>> 001ea244 PostgresMain (2, 4627d0, 4d9000, 0, 0, 1000000) + 14bc
<br>> 001bc2f0 ServerLoop (47d0c8, ffbff278, 38e800, 4b06afa6, 41d400, 3ccc00) + 9e8
<br>> 001be9ec PostmasterMain (0, 45c180, 3cd800, 461698, 45a400, 3ccc00) + 11c0
<br>> 00164550 main (1, 45c180, ffbff8fc, 459bac, fee60100, fee60140) + 39c
<br>> 0004ace4 _start (0, 0, 0, 0, 0, 0) + 5c
<br>>
<br>> I am able to reproduce the crash consistently with a query such as the following:
<br>>
<br>> select st_isvalid('POLYGON ((-97.78351581 36.37358796 1167.172,-97.78303556 36.37319923 1167.172,-98.88351581 38.37358796 1167.172))');
<br>>
<br>> Any help would be greatly appreciated!!
<br>>
<br>> Brad
<br>>
<br>> _______________________________________________
<br>> postgis-users mailing list
<br>> <a href="http://old.nabble.com/user/SendEmail.jtp?type=post&post=26449125&i=1" target="_top" rel="nofollow">postgis-users@...</a>
<br>> <a href="http://postgis.refractions.net/mailman/listinfo/postgis-users" target="_top" rel="nofollow">http://postgis.refractions.net/mailman/listinfo/postgis-users</a><br>>
</div>_______________________________________________
<br>postgis-users mailing list
<br><a href="http://old.nabble.com/user/SendEmail.jtp?type=post&post=26449125&i=2" target="_top" rel="nofollow">postgis-users@...</a>
<br><a href="http://postgis.refractions.net/mailman/listinfo/postgis-users" target="_top" rel="nofollow">http://postgis.refractions.net/mailman/listinfo/postgis-users</a><br>tag:old.nabble.com,2006:post-26449052Server Crash on Invalid 3d Polygon2009-11-20T11:31:20Z2009-11-20T11:31:20ZBrad HecoxI am having a problem were I get a server crash any time I try to validate or generate a geometry on an invalid 3D polygon. I have a set of polygons in EWKT format that I need to validate and convert to geometries if valid. In most cases the st_isvalid function returns a true or false, however of the polygon is 3D and is invalid rather than returning a false, the system crashes. I have been able to reproduce this on Solaris x86, Solaris SPARC and Max OSX. We are running PostgreSQL 8.4.1 with PostGIS 1.4.0. The following is the pstack from the core dump:
<br><br>ore 'core' of 4774: /opt/postgres/pgserver/8.4.1/bin/postgres
<br> fef68610 strncat (4d8e38, 0, 40923bff, 28, 0, 0) + 20
<br> fedd3360 pg_parser_errhint (ffbfd9e8, 4d8e38, 7, fee02a08, fee026d0, fee34224) + 24
<br> fedd9c74 LWGEOM_in (ffbfda70, 5167bc, d3280f32, 406000, 0, 3) + 70
<br> 002a0ef4 InputFunctionCall (ffbfdce8, 4d8e38, 4010, ffffffff, 516688, 4010) + 58
<br> 002a1120 OidInputFunctionCall (4012, 4d8e38, 4010, ffffffff, 0, 445000) + 20
<br> 000e43e8 coerce_type (4d9150, 4d90b0, 2c1, 4010, ffffffff, 0) + 3d4
<br> 000da760 make_fn_arguments (4d8d60, 4d90d8, ffbfdef0, 4d9118, 0, 0) + 5c
<br> 000db534 ParseFuncOrColumn (4d8d60, 4d8e20, 4d90e8, 2, 0, 0) + 3d8
<br> 000d7d00 transformExpr (0, 4d8f00, 4d9020, 31, 1518, 4da538) + 150
<br> 000e6cb8 transformTargetEntry (4d8d60, 4d8f00, 0, 0, 0, 2b65c4) + 80
<br> 000e6d0c transformTargetList (0, 4d8f50, 398ea0, fee62a00, 1e70, 0) + 44
<br> 000bc6a4 transformStmt (4d9028, 4d8f78, 0, 11ab4, 529f2c, 529f2c) + e0
<br> 000bed78 parse_analyze (4d8f78, 4d8550, 0, 0, 41e400, 529f00) + 20
<br> 001e8664 pg_analyze_and_rewrite (4d8f78, 4d8550, 0, 0, 4627d0, 462748) + 24
<br> 001ea244 PostgresMain (2, 4627d0, 4d9000, 0, 0, 1000000) + 14bc
<br> 001bc2f0 ServerLoop (47d0c8, ffbff278, 38e800, 4b06afa6, 41d400, 3ccc00) + 9e8
<br> 001be9ec PostmasterMain (0, 45c180, 3cd800, 461698, 45a400, 3ccc00) + 11c0
<br> 00164550 main (1, 45c180, ffbff8fc, 459bac, fee60100, fee60140) + 39c
<br> 0004ace4 _start (0, 0, 0, 0, 0, 0) + 5c
<br><br>I am able to reproduce the crash consistently with a query such as the following:
<br><br>select st_isvalid('POLYGON ((-97.78351581 36.37358796 1167.172,-97.78303556 36.37319923 1167.172,-98.88351581 38.37358796 1167.172))');
<br><br>Any help would be greatly appreciated!!
<br><br>Brad
<br><br>_______________________________________________
<br>postgis-users mailing list
<br><a href="http://old.nabble.com/user/SendEmail.jtp?type=post&post=26449052&i=0" target="_top" rel="nofollow">postgis-users@...</a>
<br><a href="http://postgis.refractions.net/mailman/listinfo/postgis-users" target="_top" rel="nofollow">http://postgis.refractions.net/mailman/listinfo/postgis-users</a><br>tag:old.nabble.com,2006:post-26445935Re: how many vertices in a linestring2009-11-20T08:07:18Z2009-11-20T08:07:18ZMaxime van NoppenMarco Lechner - FOSSGIS e.V. wrote:
<br>> exactly. thnx. by the way: what's the difference between st_numpoints()
<br>> and st_npoints()? The results in my case are the same
<br><br>From the doc of ST_NumPoints :
<br><br>« From 1.4 forward this is an alias for ST_NPoints which returns number
<br>of vertexes for not just line strings. Consider using ST_NPoints instead
<br>which is multi-purpose and works with many geometry types. »
<br><br>--
<br>Maxime
<br>_______________________________________________
<br>postgis-users mailing list
<br><a href="http://old.nabble.com/user/SendEmail.jtp?type=post&post=26445935&i=0" target="_top" rel="nofollow">postgis-users@...</a>
<br><a href="http://postgis.refractions.net/mailman/listinfo/postgis-users" target="_top" rel="nofollow">http://postgis.refractions.net/mailman/listinfo/postgis-users</a><br>tag:old.nabble.com,2006:post-26445872Re: how many vertices in a linestring2009-11-20T08:07:16Z2009-11-20T08:07:16ZMarco Lechner - FOSSGIS e.V.Maxime van Noppen schrieb:
<div class='shrinkable-quote'><br>> Marco Lechner - FOSSGIS e.V. wrote:
<br>>
<br>>> is there a function that returns the number of points of a geometry?
<br>>>
<br>>> LINESTRING(9.6567389 47.8442843,9.6572511 47.8443736,9.6570058
<br>>> 47.8449885,9.6565007 47.8449207,9.6567378 47.8442853) => 5
<br>>> LINESTRING(8.9289313 48.6910326,8.9289313 48.6910326) => 2
<br>>> ...
<br>>>
<br>>
<br>> Hi,
<br>>
<br>> You might want to look at ST_NPoints and ST_Numpoints for linestrings :
<br>>
<br>> <a href="http://postgis.refractions.net/documentation/manual-1.4/ST_NumPoints.html" target="_top" rel="nofollow">http://postgis.refractions.net/documentation/manual-1.4/ST_NumPoints.html</a><br>> <a href="http://postgis.refractions.net/documentation/manual-1.4/ST_NPoints.html" target="_top" rel="nofollow">http://postgis.refractions.net/documentation/manual-1.4/ST_NPoints.html</a><br>>
<br>> However for geometries in general it's more complicated
</div>exactly. thnx. by the way: what's the difference between st_numpoints()
<br>and st_npoints()? The results in my case are the same
<br>_______________________________________________
<br>postgis-users mailing list
<br><a href="http://old.nabble.com/user/SendEmail.jtp?type=post&post=26445872&i=0" target="_top" rel="nofollow">postgis-users@...</a>
<br><a href="http://postgis.refractions.net/mailman/listinfo/postgis-users" target="_top" rel="nofollow">http://postgis.refractions.net/mailman/listinfo/postgis-users</a><br>tag:old.nabble.com,2006:post-26445600Re: how many vertices in a linestring2009-11-20T07:48:43Z2009-11-20T07:48:43ZMaxime van NoppenMarco Lechner - FOSSGIS e.V. wrote:
<br>> is there a function that returns the number of points of a geometry?
<br>>
<br>> LINESTRING(9.6567389 47.8442843,9.6572511 47.8443736,9.6570058
<br>> 47.8449885,9.6565007 47.8449207,9.6567378 47.8442853) => 5
<br>> LINESTRING(8.9289313 48.6910326,8.9289313 48.6910326) => 2
<br>> ...
<br><br>Hi,
<br><br>You might want to look at ST_NPoints and ST_Numpoints for linestrings :
<br><br><a href="http://postgis.refractions.net/documentation/manual-1.4/ST_NumPoints.html" target="_top" rel="nofollow">http://postgis.refractions.net/documentation/manual-1.4/ST_NumPoints.html</a><br><a href="http://postgis.refractions.net/documentation/manual-1.4/ST_NPoints.html" target="_top" rel="nofollow">http://postgis.refractions.net/documentation/manual-1.4/ST_NPoints.html</a><br><br>However for geometries in general it's more complicated.
<br><br>--
<br>Maxime
<br>_______________________________________________
<br>postgis-users mailing list
<br><a href="http://old.nabble.com/user/SendEmail.jtp?type=post&post=26445600&i=0" target="_top" rel="nofollow">postgis-users@...</a>
<br><a href="http://postgis.refractions.net/mailman/listinfo/postgis-users" target="_top" rel="nofollow">http://postgis.refractions.net/mailman/listinfo/postgis-users</a><br>tag:old.nabble.com,2006:post-26445406how many vertices in a linestring2009-11-20T07:41:29Z2009-11-20T07:41:29ZMarco Lechner - FOSSGIS e.V.is there a function that returns the number of points of a geometry?
<br><br>LINESTRING(9.6567389 47.8442843,9.6572511 47.8443736,9.6570058
<br>47.8449885,9.6565007 47.8449207,9.6567378 47.8442853) => 5
<br>LINESTRING(8.9289313 48.6910326,8.9289313 48.6910326) => 2
<br>...
<br><br>Marco
<br>_______________________________________________
<br>postgis-users mailing list
<br><a href="http://old.nabble.com/user/SendEmail.jtp?type=post&post=26445406&i=0" target="_top" rel="nofollow">postgis-users@...</a>
<br><a href="http://postgis.refractions.net/mailman/listinfo/postgis-users" target="_top" rel="nofollow">http://postgis.refractions.net/mailman/listinfo/postgis-users</a><br>tag:old.nabble.com,2006:post-26440921Re: TopologyException: no outgoing dirEdge found2009-11-20T04:18:50Z2009-11-20T04:18:50ZJohn PowellI was having this error and in my case it was to do with consecutive identical polygons.
<br><br>I ran
<br><br>delete from geom_table using (select t1.id as delid from geom_table t1, geom_table t2 where t1.id<>t2.id and st_equals(t1.g,t2.g)='t' group by t1.id) as dupes where geom_table.id=dupes.delid;
<br><br>to remove all the duplicates. After that ST_Multi(ST_Union(geom)) worked as expected.
<br><br><br><blockquote class="quote light-black dark-border-color"><div class="quote light-border-color">
<div class="quote-author" style="font-weight: bold;">Andreas Laggner wrote:</div>
<div class="quote-message shrinkable-quote">Hi,
<br><br>does someone know what´s behind this Error?
<br>NOTICE: TopologyException: no outgoing dirEdge found
<br>(3.5888e+006,5.82215e+006)
<br>ERROR: GEOS union() threw an error!
<br><br>I tried the functions Intersection() and Union() with different datasets
<br>from differnet regions (65 to 6000 polygons) in many combinations and i
<br>always get this error (even if the data isvalid and issimple). With
<br>google search i found some hints: Maybe to high precision is the reason
<br>for this???
<br>I want to intersect whole germany with postgis and i did not realize a
<br>small municipality so far! I tried WinXP with PostgreSQL8.1 and Linux
<br>with PostgreSQL8.2 (new installation with highest version of postgis
<br>available).
<br><br>My SQL statements:
<br>create table intersect_gf2(gid integer) with oids;
<br>select
<br>addgeometrycolumn('','intersect_gf2','the_geom','31467','MULTIPOLYGON',2);
<br>insert into intersect_gf2 select t1.gid,geomunion(t1.the_geom,
<br>t2.the_geom) from gf_veg1 t1, gf_corine t2;
<br>insert into intersect_gf2 select t1.gid,intersection(t1.the_geom,
<br>t2.the_geom) from gf_buk t1, gf_veg2 t2;
<br><br>hasta luego Andreas
<br>_______________________________________________
<br>postgis-users mailing list
<br>postgis-users@postgis.refractions.net
<br><a href="http://postgis.refractions.net/mailman/listinfo/postgis-users" target="_top" rel="nofollow">http://postgis.refractions.net/mailman/listinfo/postgis-users</a></div>
</div></blockquote>
tag:old.nabble.com,2006:post-26437075Re: st_union exception2009-11-20T02:30:12Z2009-11-20T02:30:12ZJohn PowellI have also seen this error and I believe it is to do with identical geometries being unioned. I cleaned my dataset using someting like
<br><br>select t1.id from points t1, points t2 where t1.id<>t2.id and st_equals(t1.geom,t2.geom)='t' group by t1.id
<br><br>to find all the duplicates which i then deleted. The st_union(geom) function then worked as expected.
<br><br><br><br><blockquote class="quote light-black dark-border-color"><div class="quote light-border-color">
<div class="quote-author" style="font-weight: bold;">Mark Mitchell-5 wrote:</div>
<div class="quote-message shrinkable-quote">Anyone get an error similar to this when unioning a bunch of polygons
<br>with st_union?
<br><br>NOTICE: TopologyException: no outgoing dirEdge found (-149.073,61.6101,61.6101)
<br>ERROR: GEOS union() threw an error!
<br><br>I must say before I get into the details only a few of you will read
<br>that portion that I've used postgis for a while now (since 1.1.5) and
<br>it is *PHENOMENAL*. I don't want my only posts to this list in that
<br>time frame to be negative -- can't tell you how much
<br>faster/better/intuitive I've found postgis over very expensive
<br>commercial alternatives :-).
<br><br><br>My query is similar to the following:
<br><br>SELECT st_multi(st_union(poly)) from (select ... (a bunch of polys) ... ) a;
<br><br>The subselect is operating on a fairly dynamic table, and I've gone
<br>months without any problems (with many different sets of polygons too)
<br><br>I tried st_isvalid() on all polygons returned from the subselect, and
<br>all are valid, and are all non-null...
<br>If I have duplicate (or close to duplicate) polygons returned, is that
<br>going to potentially screw up st_union()?
<br><br><br>I'm running Ubuntu server (x64) 8.04 LTS, and am using the packages
<br>for GEOS and postgis from their universe repo:
<br><br>postgis 1.3.3
<br>libgeos 2.2.3
<br><br>(postgresql 8.3)
<br><br><br>Thanks for the help!
<br><br>-Mark
<br>_______________________________________________
<br>postgis-users mailing list
<br>postgis-users@postgis.refractions.net
<br><a href="http://postgis.refractions.net/mailman/listinfo/postgis-users" target="_top" rel="nofollow">http://postgis.refractions.net/mailman/listinfo/postgis-users</a></div>
</div></blockquote>
tag:old.nabble.com,2006:post-26430745Re: Postgis and arcgis with data interoperability extension2009-11-19T09:41:37Z2009-11-19T09:41:37ZParagon Corporation-2Thomas,
<br>Try changing the client encoding of your db to latin1. That's what I do.
<br><br>ALTER DATABASE mydb SET client_encoding=latin1;
<br><br><br>Leo
<br><br>-----Original Message-----
<br>From: <a href="http://old.nabble.com/user/SendEmail.jtp?type=post&post=26430745&i=0" target="_top" rel="nofollow">postgis-users-bounces@...</a>
<br>[mailto:<a href="http://old.nabble.com/user/SendEmail.jtp?type=post&post=26430745&i=1" target="_top" rel="nofollow">postgis-users-bounces@...</a>] On Behalf Of Thomas
<br>Sent: Thursday, November 19, 2009 10:47 AM
<br>To: PostGIS Users Discussion
<br>Subject: Re: [postgis-users] Postgis and arcgis with data interoperability
<br>extension
<br><br><br><br>Maybe I have an idea how. to solve the Problem: The error message that i get
<br>from the database is that the query is using the wrong encoding: IBM/latin
<br>but in the database we are using UTF8 Does anybody know how to change that
<br><br><br>Regards
<br>Thomas
<br><br><br>On 19.11.09 10:50, "Andreas Neumann" <<a href="http://old.nabble.com/user/SendEmail.jtp?type=post&post=26430745&i=2" target="_top" rel="nofollow">a.neumann@...</a>> wrote:
<br><div class='shrinkable-quote'><br>> Hi Thomas,
<br>>
<br>> did you have a look at ziggis? It enables data connection and editing of
<br>> native Postgis geometries without the need for ArcSDE.
<br>>
<br>> It is commercial, but much cheaper than SDE: <a href="http://pub.obtusesoft.com/" target="_top" rel="nofollow">http://pub.obtusesoft.com/</a><br>>
<br>> I don't have ESRI myself, so I don't have experience with it, but it is
<br>> definitely worth a try.
<br>>
<br>> Andreas
<br>>
<br>>
<br>> On Thu, November 19, 2009 10:19 am, Thomas wrote:
<br>>>
<br>>> Ok thanks for the OLE DB info.
<br>>> I am still new to working with databases so maybe I am wrong, but ESRI
<br>>> told
<br>>> me that I can use ArcSDE only with the Server Version that I until now
</div>don
<div class='shrinkable-quote'><br>>> t
<br>>> have. But still it should be possible to connect through the Data
<br>>> Interoperability Extension.
<br>>> The strange thing is that it works with polygons but not with point
<br>>> geometry
<br>>>
<br>>> Thanks for your ideas
<br>>>
<br>>> thomas
<br>>>
<br>>>
<br>>> On 18.11.09 22:51, "George Silva" <<a href="http://old.nabble.com/user/SendEmail.jtp?type=post&post=26430745&i=3" target="_top" rel="nofollow">georger.silva@...</a>> wrote:
<br>>>
<br>>>> This is very weird. Very weird. I have experience with Esri SDE and if
<br>>>> you are
<br>>>> working with OGC standarts everything should work fine.
<br>>>>
<br>>>> Did you registered all your tables trough ArcSDE? If not, you need to do
<br>>>> it.
<br>>>>
<br>>>> OLE Connections will only let you handle tabular information, so let
<br>>>> that go.
<br>>>> It's not the way.
<br>>>>
<br>>>> AFAIK ESRI implements OGC standarts. You should NOT have problems with
<br>>>> point
<br>>>> geometry.
<br>>>>
<br>>>> Att
<br>>>>
<br>>>> George
<br>>>>
<br>>>> On Wed, Nov 18, 2009 at 6:39 PM, Ragi Y. Burhum <<a href="http://old.nabble.com/user/SendEmail.jtp?type=post&post=26430745&i=4" target="_top" rel="nofollow">rburhum@...</a>>
<br>>>> wrote:
<br>>>>>
<br>>>>>> From: Thomas Jakubicka <<a href="http://old.nabble.com/user/SendEmail.jtp?type=post&post=26430745&i=5" target="_top" rel="nofollow">crazyh@...</a>>
<br>>>>>> Subject: [postgis-users] Postgis and arcgis with data interoperability
<br>>>>>> extension
<br>>>>>> To: <<a href="http://old.nabble.com/user/SendEmail.jtp?type=post&post=26430745&i=6" target="_top" rel="nofollow">postgis-users@...</a>>
<br>>>>>> Message-ID: <<a href="http://old.nabble.com/user/SendEmail.jtp?type=post&post=26430745&i=7" target="_top" rel="nofollow">C729D2DF.2936%crazyh@...</a>
<br>>>>>> <mailto:<a href="http://old.nabble.com/user/SendEmail.jtp?type=post&post=26430745&i=8" target="_top" rel="nofollow">C729D2DF.2936%25crazyh@...</a>> >
<br>>>>>> Content-Type: text/plain; charset="US-ASCII"
<br>>>>>>
<br>>>>>>
<br>>>>>> Hello,
<br>>>>>>
<br>>>>>> I am trying to access my postgis/postgreSQL database through the
<br>>>>>> datainteroperability extension and have problems with all my point
<br>>>>>> data. I
<br>>>>>> can load and export polygons but no points. The connection to the
<br>>>>>> database
<br>>>>>> works fine but when i try to load point data the attribute table is
<br>>>>>> empty
<br>>>>>> and only the header row is shown
<br>>>>>
<br>>>>> That seems bizarre. FYI, ESRI provides their own implementation of a
<br>>>>> Geometry
<br>>>>> type on PostgreSQL. You can also use PostGIS. My "psychic debugging"
<br>>>>> skills
<br>>>>> tell me there is a chance you may be mixing the two?
<br>>>>>
<br>>>>> Also, did you load the PostGIS functions as described in the
<br>>>>> documentation
<br>>>>> for the database you are working with *before* loading the data?
<br>>>>>
<br>>>>>>
<br>>>>>> I also tried to access the database through an OLE DB connection witch
<br>>>>>> is
<br>>>>>> working so far that I get access to all tables stored in the database
<br>>>>>> but
<br>>>>>> still can not show the spatial information. With this connection all
<br>>>>>> datasets are shown with the table-symbol and the geometry field stays
<br>>>>>> empty.
<br>>>>>
<br>>>>> This is how the ESRI OLE DB provider works - it doesn't load and
<br>>>>> spatial
<br>>>>> information. Anything you load from there will be "just a table"
<br>>>>>
<br>>>>>> We are using the Postgis/PostgreSQL version 8.3 and ArcEditor 9.3.1
<br>>>>>>
<br>>>>>> Thanks for your help
<br>>>>>>
<br>>>>>
<br>>>>> The best thing I could recommend to you now is to do some logging from
<br>>>>> the
<br>>>>> PostgreSQL side of things. See what SQL is being executed during the
<br>>>>> load and
<br>>>>> it should show you your problem right away.
<br>>>>>
<br>>>>> Hope this helps,
<br>>>>>
<br>>>>> - Ragi
<br>>>>>
<br>>>>> _______________________________________________
<br>>>>> postgis-users mailing list
<br>>>>> <a href="http://old.nabble.com/user/SendEmail.jtp?type=post&post=26430745&i=9" target="_top" rel="nofollow">postgis-users@...</a>
<br>>>>> <a href="http://postgis.refractions.net/mailman/listinfo/postgis-users" target="_top" rel="nofollow">http://postgis.refractions.net/mailman/listinfo/postgis-users</a><br>>>>>
<br>>>>
<br>>>>
<br>>>
<br>>> _______________________________________________
<br>>> postgis-users mailing list
<br>>> <a href="http://old.nabble.com/user/SendEmail.jtp?type=post&post=26430745&i=10" target="_top" rel="nofollow">postgis-users@...</a>
<br>>> <a href="http://postgis.refractions.net/mailman/listinfo/postgis-users" target="_top" rel="nofollow">http://postgis.refractions.net/mailman/listinfo/postgis-users</a><br>>>
<br>>
</div><br><br>_______________________________________________
<br>postgis-users mailing list
<br><a href="http://old.nabble.com/user/SendEmail.jtp?type=post&post=26430745&i=11" target="_top" rel="nofollow">postgis-users@...</a>
<br><a href="http://postgis.refractions.net/mailman/listinfo/postgis-users" target="_top" rel="nofollow">http://postgis.refractions.net/mailman/listinfo/postgis-users</a><br><br><br>_______________________________________________
<br>postgis-users mailing list
<br><a href="http://old.nabble.com/user/SendEmail.jtp?type=post&post=26430745&i=12" target="_top" rel="nofollow">postgis-users@...</a>
<br><a href="http://postgis.refractions.net/mailman/listinfo/postgis-users" target="_top" rel="nofollow">http://postgis.refractions.net/mailman/listinfo/postgis-users</a><br>tag:old.nabble.com,2006:post-26428576Re: Postgis and arcgis with data interoperability extension2009-11-19T07:47:09Z2009-11-19T07:47:09ZThomas Jakubicka<br><br>Maybe I have an idea how. to solve the Problem: The error message that i get
<br>from the database is that the query is using the wrong encoding: IBM/latin
<br>but in the database we are using UTF8
<br>Does anybody know how to change that
<br><br><br>Regards
<br>Thomas
<br><br><br>On 19.11.09 10:50, "Andreas Neumann" <<a href="http://old.nabble.com/user/SendEmail.jtp?type=post&post=26428576&i=0" target="_top" rel="nofollow">a.neumann@...</a>> wrote:
<br><div class='shrinkable-quote'><br>> Hi Thomas,
<br>>
<br>> did you have a look at ziggis? It enables data connection and editing of
<br>> native Postgis geometries without the need for ArcSDE.
<br>>
<br>> It is commercial, but much cheaper than SDE: <a href="http://pub.obtusesoft.com/" target="_top" rel="nofollow">http://pub.obtusesoft.com/</a><br>>
<br>> I don't have ESRI myself, so I don't have experience with it, but it is
<br>> definitely worth a try.
<br>>
<br>> Andreas
<br>>
<br>>
<br>> On Thu, November 19, 2009 10:19 am, Thomas wrote:
<br>>>
<br>>> Ok thanks for the OLE DB info.
<br>>> I am still new to working with databases so maybe I am wrong, but ESRI
<br>>> told
<br>>> me that I can use ArcSDE only with the Server Version that I until now don
<br>>> t
<br>>> have. But still it should be possible to connect through the Data
<br>>> Interoperability Extension.
<br>>> The strange thing is that it works with polygons but not with point
<br>>> geometry
<br>>>
<br>>> Thanks for your ideas
<br>>>
<br>>> thomas
<br>>>
<br>>>
<br>>> On 18.11.09 22:51, "George Silva" <<a href="http://old.nabble.com/user/SendEmail.jtp?type=post&post=26428576&i=1" target="_top" rel="nofollow">georger.silva@...</a>> wrote:
<br>>>
<br>>>> This is very weird. Very weird. I have experience with Esri SDE and if
<br>>>> you are
<br>>>> working with OGC standarts everything should work fine.
<br>>>>
<br>>>> Did you registered all your tables trough ArcSDE? If not, you need to do
<br>>>> it.
<br>>>>
<br>>>> OLE Connections will only let you handle tabular information, so let
<br>>>> that go.
<br>>>> It's not the way.
<br>>>>
<br>>>> AFAIK ESRI implements OGC standarts. You should NOT have problems with
<br>>>> point
<br>>>> geometry.
<br>>>>
<br>>>> Att
<br>>>>
<br>>>> George
<br>>>>
<br>>>> On Wed, Nov 18, 2009 at 6:39 PM, Ragi Y. Burhum <<a href="http://old.nabble.com/user/SendEmail.jtp?type=post&post=26428576&i=2" target="_top" rel="nofollow">rburhum@...</a>>
<br>>>> wrote:
<br>>>>>
<br>>>>>> From: Thomas Jakubicka <<a href="http://old.nabble.com/user/SendEmail.jtp?type=post&post=26428576&i=3" target="_top" rel="nofollow">crazyh@...</a>>
<br>>>>>> Subject: [postgis-users] Postgis and arcgis with data interoperability
<br>>>>>> extension
<br>>>>>> To: <<a href="http://old.nabble.com/user/SendEmail.jtp?type=post&post=26428576&i=4" target="_top" rel="nofollow">postgis-users@...</a>>
<br>>>>>> Message-ID: <<a href="http://old.nabble.com/user/SendEmail.jtp?type=post&post=26428576&i=5" target="_top" rel="nofollow">C729D2DF.2936%crazyh@...</a>
<br>>>>>> <mailto:<a href="http://old.nabble.com/user/SendEmail.jtp?type=post&post=26428576&i=6" target="_top" rel="nofollow">C729D2DF.2936%25crazyh@...</a>> >
<br>>>>>> Content-Type: text/plain; charset="US-ASCII"
<br>>>>>>
<br>>>>>>
<br>>>>>> Hello,
<br>>>>>>
<br>>>>>> I am trying to access my postgis/postgreSQL database through the
<br>>>>>> datainteroperability extension and have problems with all my point
<br>>>>>> data. I
<br>>>>>> can load and export polygons but no points. The connection to the
<br>>>>>> database
<br>>>>>> works fine but when i try to load point data the attribute table is
<br>>>>>> empty
<br>>>>>> and only the header row is shown
<br>>>>>
<br>>>>> That seems bizarre. FYI, ESRI provides their own implementation of a
<br>>>>> Geometry
<br>>>>> type on PostgreSQL. You can also use PostGIS. My "psychic debugging"
<br>>>>> skills
<br>>>>> tell me there is a chance you may be mixing the two?
<br>>>>>
<br>>>>> Also, did you load the PostGIS functions as described in the
<br>>>>> documentation
<br>>>>> for the database you are working with *before* loading the data?
<br>>>>>
<br>>>>>>
<br>>>>>> I also tried to access the database through an OLE DB connection witch
<br>>>>>> is
<br>>>>>> working so far that I get access to all tables stored in the database
<br>>>>>> but
<br>>>>>> still can not show the spatial information. With this connection all
<br>>>>>> datasets are shown with the table-symbol and the geometry field stays
<br>>>>>> empty.
<br>>>>>
<br>>>>> This is how the ESRI OLE DB provider works - it doesn't load and
<br>>>>> spatial
<br>>>>> information. Anything you load from there will be "just a table"
<br>>>>>
<br>>>>>> We are using the Postgis/PostgreSQL version 8.3 and ArcEditor 9.3.1
<br>>>>>>
<br>>>>>> Thanks for your help
<br>>>>>>
<br>>>>>
<br>>>>> The best thing I could recommend to you now is to do some logging from
<br>>>>> the
<br>>>>> PostgreSQL side of things. See what SQL is being executed during the
<br>>>>> load and
<br>>>>> it should show you your problem right away.
<br>>>>>
<br>>>>> Hope this helps,
<br>>>>>
<br>>>>> - Ragi
<br>>>>>
<br>>>>> _______________________________________________
<br>>>>> postgis-users mailing list
<br>>>>> <a href="http://old.nabble.com/user/SendEmail.jtp?type=post&post=26428576&i=7" target="_top" rel="nofollow">postgis-users@...</a>
<br>>>>> <a href="http://postgis.refractions.net/mailman/listinfo/postgis-users" target="_top" rel="nofollow">http://postgis.refractions.net/mailman/listinfo/postgis-users</a><br>>>>>
<br>>>>
<br>>>>
<br>>>
<br>>> _______________________________________________
<br>>> postgis-users mailing list
<br>>> <a href="http://old.nabble.com/user/SendEmail.jtp?type=post&post=26428576&i=8" target="_top" rel="nofollow">postgis-users@...</a>
<br>>> <a href="http://postgis.refractions.net/mailman/listinfo/postgis-users" target="_top" rel="nofollow">http://postgis.refractions.net/mailman/listinfo/postgis-users</a><br>>>
<br>>
</div><br><br>_______________________________________________
<br>postgis-users mailing list
<br><a href="http://old.nabble.com/user/SendEmail.jtp?type=post&post=26428576&i=9" target="_top" rel="nofollow">postgis-users@...</a>
<br><a href="http://postgis.refractions.net/mailman/listinfo/postgis-users" target="_top" rel="nofollow">http://postgis.refractions.net/mailman/listinfo/postgis-users</a><br>