Re: [postgis-users] Determining file's bit encoding

2008-09-26 Thread Paul Ramsey
You don't need iconv, you can use the -W flag in shp2pgsql: shp2pgsql -W LATIN1 -s 4366 -D foo.shp foo If you're in a western country, your encoding is likely either LATIN1 or WIN1252 (they are almost identical, the WIN one has some Windows special characters in it). Here's the list of possibili

Re: [postgis-users] ST_SPHEROID_LENGTH SLOW

2008-09-27 Thread Paul Ramsey
That's deeply odd, and worth investigating. Can you try 8.3+1.3.3, so we can eliminate the PostGIS version as a variable first? Then I'll see if I can duplicate the result on my system. P. On Sat, Sep 27, 2008 at 12:59 PM, David Vaz <[EMAIL PROTECTED]> wrote: > Hi, > > I am a debian user, several

Re: [postgis-users] ST_SPHEROID_LENGTH SLOW

2008-09-27 Thread Paul Ramsey
so I stayed > with 8.2+1.3.2, this week I tried again, given that debian testing > already had postgis packages, but the speed was the same. > > Paul Ramsey wrote: >> That's deeply odd, and worth investigating. Can you try 8.3+1.3.3, so >> we can eliminate the PostGIS v

Re: [postgis-users] Re: Comparing Geometries with Different SRIDs

2008-09-28 Thread Paul Ramsey
There's sort of a grey area between RDBMS and spatial semantics if A = B does st_equals(A,B) ? same thing with null and GEOMETRY EMPTY, there is similarity but the ideas are not identical I would expect A = B to return false if SRIDs are not the same, but ST_Equals(A,B) to error out if SRIDs

Re: [postgis-users] GeometryCollection(Empty) is SRID dependent (orsupporting nulls in st_union)

2008-09-29 Thread Paul Ramsey
I think we all buy that. File it. P On Mon, Sep 29, 2008 at 4:48 PM, Charlie Savage <[EMAIL PROTECTED]> wrote: >> So - do we have agreement that NULL's should be ignored in functions that > act on aggregates? > > To recap: > > * NULL values cause functions such as ST_UNION to return NULL > * This

Re: [postgis-users] Bug Converting from WGS 84 to UTM ?

2008-09-29 Thread Paul Ramsey
Use 26910. The reference you chose is for the southern hemisphere (hence the "S"). P On Mon, Sep 29, 2008 at 7:55 PM, Oliver Monson <[EMAIL PROTECTED]> wrote: > Hi All, > > I'm not sure if I'm seeing some random bug, or just not doing/understanding > something correctly? > > I have a list of poi

Re: [postgis-users] lwpostgis.sql and liblwgeom.so

2008-10-02 Thread Paul Ramsey
You've got a mismatch between your GEOS version and your postgis version. Try upgrading your geos version. P On Thu, Oct 2, 2008 at 7:18 AM, ROBERT HOLLINGSWORTH <[EMAIL PROTECTED]> wrote: > I'm new to building postgresql and postgis, and have run into a problem > that, from googles on the subjec

[postgis-users] Re: [postgis-devel] Dropped DM (Time) dimension with intersections

2008-10-07 Thread Paul Ramsey
3-D geometries without > spending a fortune on proprietary stuff which is not a PostGIS issue, but > has to gain in momentum to make PostGIS 3D more powerful (e.g. uDig for 3D > or OpenJump for 3D or OpenLayers for 3D?) > > Thanks, > Regina > >

Re: [postgis-users] Conterra and ESRI 9.3

2008-10-08 Thread Paul Ramsey
Are you suggesting doing something like PostGIS => WFS Server => Conterra WFS Editor => ArcMap? Wouldn't using the new SDE-on-PostGIS functionality from ESRI be a little more direct? Paul On Wed, Oct 8, 2008 at 12:59 PM, <[EMAIL PROTECTED]> wrote: > > Is anyone using the Conterra.de wfsEditor

Re: [postgis-users] How to find out SRID

2008-10-09 Thread Paul Ramsey
There isn't one, you'll have to add a new entry. The PROJ4TEXT will be +proj=tmerc +lat_0=36.66 +lon_0=-88.33 +k=0.75 +x_0=30 +y_0=0 +ellps=GRS80 +datum=NAD83 +to_meter=0.3048006096012192 +no_defs On Thu, Oct 9, 2008 at 9:47 AM, <[EMAIL PROTECTED]> wrote: > Hi A

[postgis-users] GEOS Project Sponsors Wanted

2008-10-09 Thread Paul Ramsey
Do you want to help the PostGIS project get faster and remove nagging bugs? Then become a project sponsor for GEOS, the library that PostGIS uses for topological operations! GEOS project sponsorship will be used for things like: bringing the latest performance improvements in from JTS (new Buffer i

Re: [postgis-users] TRANSFORM failures with some multiploygons - data sample attached

2008-10-09 Thread Paul Ramsey
I'm not seeing this problem on your polygon. What does select postgis_full_version() say? On Thu, Oct 9, 2008 at 7:56 PM, <[EMAIL PROTECTED]> wrote: > Hi All, > > In debugging the problem below, I found that TRANSFORM would fail with some > of my city multipolygons. > > I would think that the T

Re: [postgis-users] TRANSFORM failures with some multiploygons - data sample attached

2008-10-09 Thread Paul Ramsey
>> It ran without any error. But some of the resulted values in the_geom >> column became null or empty. How did you check this? select gid from thetable where the_geom is null? select gid from thetabel where npoints(the_geom) = P ___ postgis-users mai

Re: [postgis-users] TRANSFORM failures with some multiploygons - datasample attached

2008-10-09 Thread Paul Ramsey
t; CYW > > - Original Message - From: "Paul Ramsey" <[EMAIL PROTECTED]> > To: "PostGIS Users Discussion" > Sent: Thursday, October 09, 2008 10:12 PM > Subject: Re: [postgis-users] TRANSFORM failures with some multiploygons - > datasample attached &g

Re: [postgis-users] PostGIS Transform and Quantum GIS

2008-10-10 Thread Paul Ramsey
Try removing the table from qgis and re-adding it. And don't forget to change the SRID value in GEOMETRY_COLUMNS before you do. I think you're altering the coordinates underneath QGIS's understanding of that they are (it things it has a 4326 table, and boom, you turn it into something else). P. O

Re: [postgis-users] Problem with LWGEOM_collect_garray - never returns MultiX from MultiX?

2008-10-10 Thread Paul Ramsey
OK, I agree, the comment header and the behavior differ... I'm a bit afraid to change this behavior though... it's a bit long-standing now. I wonder if doing so would cause and serious problems for people. P. On Wed, Feb 6, 2008 at 5:31 AM, Marcus Jenkins <[EMAIL PROTECTED]> wrote: > Dare I say i

[postgis-users] 1.3.4SVN Testers Wanted

2008-10-13 Thread Paul Ramsey
1.3.4 is very close to release, but we can release it faster, with less intestinal distress, if we get some solid pre-release testing from folks. You can download the current SVN version directly from the downloads page, no SVN client required: http://postgis.refractions.net/download/postgis-1.3.4

Re: [postgis-users] 1.3.4SVN Testers Wanted

2008-10-14 Thread Paul Ramsey
Got it, sorry, fell into gaps. Applied now. (BTW, unified diff (-u) is preferred, if possible in future.) On Tue, Oct 14, 2008 at 1:45 AM, Mark Cave-Ayland <[EMAIL PROTECTED]> wrote: > Paul Ramsey wrote: >> >> 1.3.4 is very close to release, but we can release it faster, wi

Re: [postgis-users] 1.3.4SVN Testers Wanted

2008-10-14 Thread Paul Ramsey
ory `/home/clh/src/postgis-1.3.svn/utils' > chmod +x postgis_restore.pl create_undef.pl postgis_proc_upgrade.pl > profile_intersects.pl test_estimation.pl test_joinestimation.pl > make[1]: Leaving directory `/home/clh/src/postgis-1.3.svn/utils' > [EMAIL PROTECTED]:~/src/postgi

Re: [postgis-users] Indices for ST_distance_sphere

2008-10-14 Thread Paul Ramsey
What about 3112 (# GDA94 / Geoscience Australia Lambert) or 3577 (Australian Albers), those should cover things fine, no? P. On Tue, Oct 14, 2008 at 3:46 PM, Stephen Baillie <[EMAIL PROTECTED]> wrote: > G'day Regina, > >> It really depends your area of interest, how far out you span out when >> l

Re: [postgis-users] Indices for ST_distance_sphere

2008-10-14 Thread Paul Ramsey
If you're in a planar projection, then select * from footable where st_dwithin(the_geom, 'POINT(foo bar)', DISTANCE) will return based on an indexed search for arbitrary foo and bar. You need to be in a planar projection so the spatial index works in the same cartesian domain as the distance calc

Re: [postgis-users] Indices for ST_distance_sphere

2008-10-14 Thread Paul Ramsey
larger variance than close ones, but whether you care depends, again, on your actual use cases. P On Tue, Oct 14, 2008 at 5:40 PM, Paul Ramsey <[EMAIL PROTECTED]> wrote: > If you're in a planar projection, then > > select * from footable where st_dwithin(the_geom, '

Re: [postgis-users] Indices for ST_distance_sphere

2008-10-14 Thread Paul Ramsey
d behavior: CREATE INDEX mygeoidx ON mytable USING GIST (mygeocolumn); Paul On Tue, Oct 14, 2008 at 6:43 PM, Stephen Baillie <[EMAIL PROTECTED]> wrote: > Paul Ramsey wrote: > >> If you're in a planar projection, then >> >> select * from footable where st_dw

Re: [postgis-users] Contains failing

2008-10-15 Thread Paul Ramsey
There's a contains bug in GEOS 3.0.0 We need to cut 3.0.1 P On Wed, Oct 15, 2008 at 9:23 AM, Fred Lehodey <[EMAIL PROTECTED]> wrote: > Hi, > I have the same results as David: (false) > postgis_full_version(): > "POSTGIS="1.3.3" GEOS="3.0.0-CAPI-1.4.1" PROJ="Rel. 4.6.0, 21 Dec 2007" > USE_STATS"

[postgis-users] GEOS 3.0.1

2008-10-15 Thread Paul Ramsey
I'm pleased to report, GEOS 3.0.1 is available: http://download.osgeo.org/geos/geos-3.0.1.tar.bz2 This removes a contains bug that has been reported here a couple time recently. P. ___ postgis-users mailing list postgis-users@postgis.refractions.net ht

Re: [postgis-users] GEOS 3.0.1

2008-10-16 Thread Paul Ramsey
Does adding #include help? What is the value of HAVE_MEMCPY in your source/headers/config..h file? P. On Thu, Oct 16, 2008 at 10:52 AM, Dylan Beaudette <[EMAIL PROTECTED]> wrote: > On Wednesday 15 October 2008, Paul Ramsey wrote: >> I'm pleased to report, GEOS 3.0.1 is a

[postgis-users] 1.3.4rc1

2008-10-19 Thread Paul Ramsey
The PostGIS team is pleased to provide release version 1.3.4rc1 for your testing delectation. You may download it from: http://postgis.refractions.net/download For extra fun, you can also download the SVN version of GEOS, which includes support for "prepared geometry" acceleration of ST_Intersect

Re: [postgis-users] Out of Memory problem for large table by ST_Contains(..)

2008-10-21 Thread Paul Ramsey
This should be fixed when 1.3.4 comes out. An rc2 will be out tomorrow you can try. P. On Tue, Oct 21, 2008 at 3:28 PM, John Zhang <[EMAIL PROTECTED]> wrote: > Hello list, > > I am writing to seek your input on how to handle such an issue: > > I have a large table containing over 3 million polygo

[postgis-users] 1.3.4rc2 Released

2008-10-24 Thread Paul Ramsey
We continue to inch closer to a 1.3.4 release. Please compile and test 1.3.4rc2, available from http://postgis.refractions.net/download/postgis-1.3.4rc2.tar.gz Thanks! Paul ___ postgis-users mailing list postgis-users@postgis.refractions.net http://pos

Re: [postgis-users] Roeck test

2008-10-24 Thread Paul Ramsey
Wouldn't the ratio of area to perimeter be a similar metric? Also, wouldn't it be easier to just stop letting politicians draw their own boundaries? :) I wonder what the algorithm for minimum enclosing circle is... (Short answer, no, because we don't have a min-circle routine.) P. On Fri, Oct 2

Re: [postgis-users] Roeck test

2008-10-24 Thread Paul Ramsey
It's an inscribed circle? I thought it was a containing one. P. On Fri, Oct 24, 2008 at 12:50 PM, Lee Meilleur <[EMAIL PROTECTED]> wrote: > > > The area of the convex hull is actually larger than the area of the polygon, > so the ratio ends up being greater than 1. With the Roeck test, the rati

Re: [postgis-users] Using ST_ (doesn't work)

2008-10-31 Thread Paul Ramsey
Did you upgrade an existing database? Create a fresh one and load postgis and then see if it works. The SQL signatures in old databases don't get updated when you update the dll/so underneath. P On Fri, Oct 31, 2008 at 9:54 AM, lisek lichu <[EMAIL PROTECTED]> wrote: > Hello, > > For last week i

Re: [postgis-users] Live PostGIS Edit tool?

2008-10-31 Thread Paul Ramsey
uDig doesn't checkin/checkout, it works directly on the records one at time. P. On Fri, Oct 31, 2008 at 9:58 AM, Sufficool, Stanley <[EMAIL PROTECTED]> wrote: > So far all the tools I have worked with (OpenJUMP, uDIG) use a check out, > check in process on PostGIS to update geometries. This is im

Re: [postgis-users] Using ST_ (doesn't work)

2008-10-31 Thread Paul Ramsey
Either you are working on an old database with the new library underneath, or you accidentally loaded the old lwpostgis.sql file into a new database. Go back and review your install procedure, you'll find the problem there. P. On Fri, Oct 31, 2008 at 12:08 PM, lisek lichu <[EMAIL PROTECTED]> wrot

Re: [postgis-users] Union, but need to exclude any void

2008-11-02 Thread Paul Ramsey
A time-honored approach to tricking flood-fill routines into honoring holes is the "zero width corridor" between the outer and inner rings, using clockwise orientation on the outer ring and counter-clockwise on the inner. So, instead of this: POLYGON((0 0, 0 10, 10 10, 10 0, 0 0), (2 2, 2 3, 3 3,

Re: [postgis-users] WKT polygon and interior ring orientation

2008-11-03 Thread Paul Ramsey
Since PostGIS explicitly keeps track of whether a ring is exterior or interior, the orientation is neither here nor there. Internally I think it matters not at all. For ST_As*() functions, I can see enforcing orientation as a useful thing, since that would remove the need for clients to do so post-

[postgis-users] PostgreSQL 8.3 GiST Bug

2008-11-12 Thread Paul Ramsey
Important information for PostGIS users! If you are on the PostgreSQL 8.3 series, upgrade to 8.3.5 as soon as possible, as earlier releases contain a GiST index bug which may cause data to occasionally not show up in spatially indexed searches. This bug will only affect tables that have had record

Re: [postgis-users] PostGIS 1.3.4 RC3 released

2008-11-13 Thread Paul Ramsey
If you re-run ./autogen yourself, do things get better? P On Thu, Nov 13, 2008 at 9:33 AM, Kevin Neufeld <[EMAIL PROTECTED]> wrote: > I'm trying to compile it with PostgreSQL 8.2.4 with no luck. > > I did this (configure and make had no errors): > ./configure --with-pgsql=/opt/pgsql-test/bin/pg_c

Re: [postgis-users] PostGIS 1.3.4 RC3 released

2008-11-13 Thread Paul Ramsey
Werks fer me on rackmount. SUMMARY --- HOST_OS: linux-gnu PGSQL: /opt/pgsql/8.2/bin/pg_config GEOS: /usr/local/bin/geos-config (with C-API) (ldflags: -L/usr/local/lib) PROJ: prefix=/usr/local libdir=/usr/local/lib ICO

Re: [postgis-users] PostGIS 1.3.4 RC3 released

2008-11-21 Thread Paul Ramsey
Tim, what did the diff file look like? tmp/pgis_reg_15519/test_38_diff P On Fri, Nov 21, 2008 at 5:50 AM, Tim Bowden <[EMAIL PROTECTED]> wrote: > > On Fri, 2008-11-21 at 12:25 +, Mark Cave-Ayland wrote: >> Tim Bowden wrote: >> >> > On Fri, 2008-11-21 at 00:42 +0900, Tim Bowden wrote: >> >> D

Re: [postgis-users] PostGIS 1.3.4 RC3 released

2008-11-21 Thread Paul Ramsey
Thanks. I guess it's possible the syntax in my regression doesn't agree with pgsql 8.1... something to check. P On Fri, Nov 21, 2008 at 4:36 PM, Tim Bowden <[EMAIL PROTECTED]> wrote: > On Fri, 2008-11-21 at 08:42 -0800, Paul Ramsey wrote: >> Tim, what did the diff fil

Re: [postgis-users] PostGIS 1.3.4 RC3 released

2008-11-24 Thread Paul Ramsey
I have fixed the regression failure on pre-8.2 pgsql, which was unrelated to PostGIS, just the use of SELECT FROM VALUES that I used in the regress_ogc_prep.sql file, and was not supported in earlier versions. P. On Mon, Nov 24, 2008 at 3:23 AM, Mark Cave-Ayland <[EMAIL PROTECTED]> wrote: > Peter

Re: [postgis-users] 1.3.4rc3 regression failure

2008-11-24 Thread Paul Ramsey
Excellent timing, I just pushed to release. An issue for 1.3.5 :) P On Mon, Nov 24, 2008 at 9:09 PM, Stephen Davies <[EMAIL PROTECTED]> wrote: > I have just built 1.3.4rc3 and run the regression tests. > > Test 17 failed with: > > PostgreSQL 8.3.5 on i686-pc-linux-gnu, compiled by GCC gcc (GCC)

[postgis-users] PostGIS 1.3.4 Released

2008-11-24 Thread Paul Ramsey
to GeoJSON format within the database. (Olivier Courtin) - Add forthcoming PostgreSQL 8.4 support (Paul Ramsey, Mark Cave-Ayland) - Improved CSS for PostGIS documentation (Dane Springmeyer) - Inclusion of new "postgis_comments.sql" file to enable detailed function comm

Re: [postgis-users] Integer (PostGIS) fields become real in QGIS!

2008-11-26 Thread Paul Ramsey
Tell the QGIS developers to expand their notion of an int to 11 digits. The whole problem is the mismatch between formal typing in C/C++ (integer, float, double) and the "typing" used in DBF headers (number(n,m)). QGIS interprets anything over 10 digits as float so that it doesn't accidentally bust

Re: [postgis-users] Union of MultiPolygons

2008-11-26 Thread Paul Ramsey
Probably your multipolygons are merging into a proper polygon. INSERT INTO county (polygon_nm, the_geom) VALUES ('X', st_multi(st_union((select the_geom from county where polygon_nm='LAKE'), (select the_geom from county where polygon_nm='COOK' Yes, st_union can operate as an aggregate. CREAT

Re: [postgis-users] Invalid geometry

2008-11-27 Thread Paul Ramsey
Try including one non-spatial field in the table you dump. P On Thu, Nov 27, 2008 at 11:33 AM, Travis Kirstine <[EMAIL PROTECTED]> wrote: > I'm having some problems with shapefiles dumped out of postgis using > the pgsql2shp utility. I am creating a geometry union using the > upgis_cascadeunion

Re: [postgis-users] geom null

2008-11-28 Thread Paul Ramsey
It's very much a client side issue, not a PostGIS issue. Ask the QGIS/gvSIG/uDig folks how they support your use case. P. On Fri, Nov 28, 2008 at 8:12 AM, ju <[EMAIL PROTECTED]> wrote: > thanks but my table contains data for 1500 geographical place we are going > to digitize > some will be point,

Re: [postgis-users] relationship functions not working well

2008-12-05 Thread Paul Ramsey
add wkb_geometry && SetSRID('BOX3D(-162 55,-161 56)'::box3d,4326) to your where clause to engage the index. P On Fri, Dec 5, 2008 at 10:26 AM, William Kyngesburye <[EMAIL PROTECTED]> wrote: > On Dec 5, 2008, at 11:22 AM, Martin Davis wrote: > >> Good news! I thought about this some more overni

Re: [postgis-users] variant documentation paths, 1.3.4 works on NetBSD-current

2008-12-10 Thread Paul Ramsey
Sorry Greg, it's manual, so variation is seeping in. Another reason to bring the scripted packaging system back up to snuff. P On Wed, Dec 10, 2008 at 3:34 PM, Greg Troxel <[EMAIL PROTECTED]> wrote: > I am updating postgis to 1.3.4 (from 1.3.3) in pkgsrc. It seems the > names of installed docs h

Re: [postgis-users] Changed bahaviour of PostGIS confuses MapServer

2008-12-13 Thread Paul Ramsey
Yes, Peter, I think it should call for a quick 1.3.5. Mapserver users are a huge contingent of PostGIS users, and anyone with a LINE layer now has a crasher. 1.3.5 should come off the line ASAP. There's always more version numbers, IMO. Can we have a quick vote on 1.3.5 release-from-branch, no-wai

Re: [postgis-users] Crossing roads

2008-12-13 Thread Paul Ramsey
Stop thinking lines and start thinking endpoints. If you convert your lines into a set of start/end points then group on x/y you'll end up an aggregation which defines your sets of road intersections. If you want until Monday, I bet Kevin or Regina will drop full SQL into your lap. I'm too lazy (

Re: [postgis-users] Problem loading shape data

2008-12-13 Thread Paul Ramsey
You have to scroll up in the window until you can see what the *actual* error is... increase your command terminal scroll buffers to 5000 then re-run and scroll up to the first error message. On Fri, Dec 12, 2008 at 3:17 PM, John Zhang wrote: > Hi there, > > I used shp2pgsql command to convert s

Re: [postgis-users] Re: Problem loading shape data (Paul Ramsey)

2008-12-15 Thread Paul Ramsey
HaidaMapping.local -d hm -U postgres -f > D:\tmp\AbaloneClosures2.sql > > Based on the error message, the table "hm_vect"."abaloneclosures2" seemed > created. However, it is not. > > Your input on the problem would be much appreciated! > > Thanks >

Re: [postgis-users] Re: Problem loading shape data (Paul Ramsey)

2008-12-15 Thread Paul Ramsey
2" seemed >> created. However, it is not. >> >> Your input on the problem would be much appreciated! >> >> Thanks >> John >> >>> Date: Sat, 13 Dec 2008 16:36:00 -0800 >>> From: "Paul Ramsey" >>> Subject: Re: [postgis-us

[postgis-users] PostGIS 1.3.5

2008-12-15 Thread Paul Ramsey
The 1.3.5 release of PostGIS is available. I am cross-posting to Mapserver because of the particular interest for that community in this release. This release fixes a minor change in 1.3.4 that caused and Mapserver LINE layer to start crashing PostGIS. If you are using Mapserver and PostGIS, do no

Re: [postgis-users] Confusing Distance

2008-12-17 Thread Paul Ramsey
Yeah, north polar stereographic, not an optimal choice for Brazil :) P On Wed, Dec 17, 2008 at 12:53 PM, David William Bitner wrote: > How did you chose to use 32661? You should use an SRID that uses a > projection that works well for the area that you are working in. > > On Wed, Dec 17, 2008 a

Re: [postgis-users] Does Slony-I work with PostGIS data?

2008-12-22 Thread Paul Ramsey
I don't have facts on the ground, just surmise: large objects likely refers to BLOBs, which PostGIS objects are *not* (though PostGIS objects are large, they are large in a confusing, different way than BLOBs). Slony should work for PostGIS. P. On Mon, Dec 22, 2008 at 7:14 AM, Peter N. Schweitzer

Re: [postgis-users] caching distance calculations

2008-12-22 Thread Paul Ramsey
Totally up to you... depends on how it fits your application. It's not a bad design pattern, by any means, and it does provide guaranteed performance, assuming your data doesn't change too much. What are you doing this distance calculation set in aid of? P On Mon, Dec 22, 2008 at 8:52 AM, Matthe

Re: [postgis-users] Strange intersection behavior

2008-12-29 Thread Paul Ramsey
Indeed. If I test the distance of your derived point from the original line, it is: 5.82594422389951e-15 We only have 64 bits of precision to work with (less, really, since the exponent and sign take up space) in a double, so the derived point will not be "on" the original line, in a mathematica

Re: [postgis-users] The Old Who is using PostGIS survey again?

2009-01-07 Thread Paul Ramsey
Fastest way to migrate: drop a couple thousand bucks on an FME license and run it once. :) If you understand Oracle Spatial, you understand PostGIS, conceptually. Just a few function names different. The drudgery of a migration though, that's something FME can cure. P. On Wed, Jan 7, 2009 at 12:

Re: [postgis-users] The Old Who is using PostGIS survey again?

2009-01-08 Thread Paul Ramsey
nect using old SQL*Net clients. P. On Jan 8, 2009, at 10:30 AM, Stephen Carville wrote: On Wednesday 07 January 2009 15:34, Paul Ramsey wrote: Fastest way to migrate: drop a couple thousand bucks on an FME license and run it once. :) If you understand Oracle Spatial, you understa

Re: [postgis-users] Help with Bad Query Plan

2009-01-08 Thread Paul Ramsey
The information contained in the EMail and any attachments is confidential and intended solely and for the attention and use of the named addressee(s). It may not be disclosed to any other person without the express authority of the HPA, or the inten

Re: [postgis-users] The Old Who is using PostGIS survey again?

2009-01-08 Thread Paul Ramsey
Another option is ogr2ogr <http://www.gdal.org/ogr/> which should allow you to do flat table transfers. P. On Jan 8, 2009, at 10:30 AM, Stephen Carville wrote: On Wednesday 07 January 2009 15:34, Paul Ramsey wrote: Fastest way to migrate: drop a couple thousand bucks on an FME l

Re: [postgis-users] Help with Bad Query Plan

2009-01-09 Thread Paul Ramsey
I just tried to load up the geometry portion of the problem only, and I am *not* seeing the same effect. Does the geometry-only part of the query also exhibit the long query time for you? select count(*) from gdors_geography g1, gdors_geography g2 where st_intersects(g1.the_geom, g2.the_geom) and

Re: [postgis-users] Encoding issue in PostGIS/PostgreSQL

2009-01-10 Thread Paul Ramsey
anks in advance. -- Thanh ___ postgis-users mailing list postgis-users@postgis.refractions.net http://postgis.refractions.net/mailman/listinfo/postgis-users -- Paul Ramsey pram...@cleverelephant.ca +1 250 885 0632 __

[postgis-users] Toronto Code Sprint: March 7-10

2009-01-11 Thread Paul Ramsey
Community members, Reminder, there is a Code Sprint event occurring this spring that members of the "C tribe" of open source GIS projects might be interested in attending. http://wiki.osgeo.org/wiki/Toronto_Code_Sprint_2009 We have space for only about five more attendees, and a couple more sp

Re: [postgis-users] How to disable password in postgreSQL (pardon the stupid question)

2009-01-12 Thread Paul Ramsey
abases. HTH, Mark. -- Mark Cave-Ayland Sirius Corporation - The Open Source Experts http://www.siriusit.co.uk T: +44 870 608 0063 ___ postgis-users mailing list postgis-users@postgis.refractions.net http://postgis.refractions.net/mailman/listinfo/postgis-us

Re: [postgis-users] Query Help

2009-01-13 Thread Paul Ramsey
SELECT o.name, o.id, s.name FROM offenders o, schools s WHERE ST_DWithing(o.the_geom, s.the_geom, 500) AND o.status = 'TRUE'; On Tue, Jan 13, 2009 at 7:56 AM, chris brisendine wrote: > I am trying to do a query where I have two point tables one is > offenders and the other is schools I want to

Re: [postgis-users] postgre arcsde9.3

2009-01-14 Thread Paul Ramsey
st_is_registered_spatial_column() is an ESRI function, best ask at an ESRI forum. P On Wed, Jan 14, 2009 at 1:03 AM, P Adji wrote: > dear all, > > ive been messing around with Arcsde + postgre 8.3.0 and postgis 1.3.2. After > i ran a command "select st_is_registered_spatial_column(bla bla)' > i

Re: [postgis-users] Query uses too much memory

2009-01-21 Thread Paul Ramsey
Upgrade to 1.3.5, it looks like you are doing point-in-polygon intersects, and there were some massive leaks there that I fixed in later revisions. P. On Wed, Jan 21, 2009 at 3:15 PM, Reid Priedhorsky wrote: > Dear all, > > I have a fairly simple (or so I thought) PostGIS query which eats my > m

Re: [postgis-users] pgsql2shp does not have an encoding option!

2009-01-25 Thread Paul Ramsey
Setting an environment variable should do the trick, but a commandline option would be an improvement Sent from my iPod On 25-Jan-09, at 10:33 AM, Stephen Woodbridge wrote: Hi all, I tend to keep all my databases in UTF8, and I just had an occasion where I needed to dump a postGIS tab

Re: [postgis-users] Shapes to PG

2009-01-26 Thread Paul Ramsey
select st_area(the_geom) from onb; You *cannot* use the visibility/non-visibility of fields in PgAdmin as any indicator of the presence/absence of geometry data in the record, because PgAdmin tosses things that are too long to display and *many* geometries fall into that category. On Mon, Jan 26,

Re: [postgis-users] Shapes to PG

2009-01-26 Thread Paul Ramsey
On Mon, Jan 26, 2009 at 10:45 AM, Neil.Young wrote: >> select count(*) from my_table where the_geom is null; > The first returns 0 of 5204 records. So you have no null records. >> select distinct st_isvalid(the_geom) from my_table; -- for OGC validity >> etc. > The output of the latter I don't u

Re: [postgis-users] Some questions concerning postgis

2009-01-26 Thread Paul Ramsey
On Mon, Jan 26, 2009 at 3:54 PM, Paragon Corporation wrote: > Though if you want to maintain borders I guess you would probably want to > 1) convert your polygons to linestrings using ST_Boundary, > 2) dump out into edges and grouping by the gid > 3) Simplify and then repolygonize. It is hard. K

Re: [postgis-users] help optimize distance calculation

2009-01-27 Thread Paul Ramsey
Is there a maximum expected mindistance? select min(distance(f.the_geom,g.the_geom)) from features f, grid g where st_dwithin(f,g,[maxmindistance]) group by ao.gid P. On Tue, Jan 27, 2009 at 3:56 PM, Sky wrote: > Hello all, I'm new to postgis and I'm hoping someone might help advise me! > > I h

Re: [postgis-users] help optimize distance calculation

2009-01-28 Thread Paul Ramsey
Judging by your index condition in your query, you are restricting your test to just one grid cell (gid = 1?) and PostgreSQL doesn't use indexes willynilly, it uses them when it thinks they will help. You only have 2000 features in the other candidate table, PostgreSQL is guessing that a seqscan wi

Re: [postgis-users] help optimize distance calculation

2009-01-28 Thread Paul Ramsey
On Wed, Jan 28, 2009 at 10:41 AM, Sky wrote: > Again, the query doesn't seem to use the ao.the_geom index... correct? > Here is the query I used: > > explain analyse select ao.gid, min(distance(f.the_geom,ao.the_geom)) > from features f, ao ao group by ao.gid limit 40 > > This query should u

[postgis-users] Call for Semantics (ST_CollectionTypes())

2009-01-28 Thread Paul Ramsey
I would like to add a function that takes an anonymous GEOMETRYCOLLECTION and returns something more structured: - a GEOMETRYCOLLECTION where the first element is a MULTLIPOLYGON, the second a MULTILINESTRING and the third a MULTIPOINT? - a tuple, with elements 'polygons', 'linestrings', 'points',

Re: [postgis-users] Call for Semantics (ST_CollectionTypes())

2009-01-28 Thread Paul Ramsey
casionally happens. P. On Wed, Jan 28, 2009 at 4:40 PM, Mark Cave-Ayland wrote: > Paul Ramsey wrote: > >> I would like to add a function that takes an anonymous >> GEOMETRYCOLLECTION and returns something more structured: >> - a GEOMETRYCOLLECTION where the first element i

Re: [postgis-users] PostGIS analogues to ArcInfo geoprocessing functions

2009-01-29 Thread Paul Ramsey
On Thu, Jan 29, 2009 at 8:42 AM, Martin Davis wrote: > The basic approach to computing polygon overlays has been well-understood > for a long time (although this does not imply well-documented!). The > implementation however is quite tricky, especially if performance and > robustness is required.

[postgis-users] More Semantics: SRID Matching

2009-01-29 Thread Paul Ramsey
I just did this: aggtest=# select name from tm_world_2 where st_dwithin(the_geom, 'POINT(32.4122 -21.2178)', 0.0001); ERROR: Operation on two geometries with different SRIDs And you know what, that seems a bit harsh to me. Given an operation where one SRID is known (st_srid(the_geom) == 4326), a

Re: [postgis-users] ST_Intersects of a diagonal buffer on a polygon grid

2009-01-29 Thread Paul Ramsey
Bruteforce-wise, upgrading to 1.4 and GEOS 3.1 will bring you prepared geometries which will make things faster. However, given your basic apparent question "what cells are within N distance of this line?" a much faster approach, which doesn't require upgrading, is as follows: First, don't use ju

Re: [postgis-users] PostGIS analogues to ArcInfo geoprocessing functions

2009-01-29 Thread Paul Ramsey
>> I can't even tell if it is memory-only or whether it uses external >> files to hold intermediate results! >> >> I look forward to some cleverer programmer than me figuring out how to >> integrate the Grass code with PostGIS... >> >> Chris Herman

Re: [postgis-users] portable spatial db

2009-01-30 Thread Paul Ramsey
If the versions are the same, you should be able to move the instance from place to place (note, I say "instance" not database, instances are the collections of databases that reside in one data area) initdb on your SATA start your pgsql, with -D to your SATA load your data stop your pgsql remove

Re: [postgis-users] portable spatial db

2009-01-30 Thread Paul Ramsey
---Original Message- > From: postgis-users-boun...@postgis.refractions.net > [mailto:postgis-users-boun...@postgis.refractions.net] On Behalf Of Paul > Ramsey > Sent: Friday, January 30, 2009 11:57 AM > To: PostGIS Users Discussion > Subject: Re: [postgis-users] portable spatial db > > If the

Re: [postgis-users] portable spatial db

2009-01-30 Thread Paul Ramsey
Same endian / hardware architecture required. Things are only transformed to architecture independent forms on the way out of database. On Fri, Jan 30, 2009 at 9:23 AM, Greg Troxel wrote: > Paul Ramsey writes: > >> If the versions are the same, you should be able to move the inst

Re: [postgis-users] portable spatial db

2009-01-30 Thread Paul Ramsey
Bruce > > > -Original Message- > From: postgis-users-boun...@postgis.refractions.net > [mailto:postgis-users-boun...@postgis.refractions.net] On Behalf Of Greg > Troxel > Sent: Friday, January 30, 2009 12:24 PM > To: PostGIS Users Discussion > Subject: Re: [postgis-users] portable

Re: [postgis-users] why isn't st_translate working for points and polys in same dataset?

2009-01-30 Thread Paul Ramsey
Your constraint is an SRID constraint, and you aren't setting your SRID when you construct your new points. Try this: truncate stands; insert into stands select tmp_stands.gid as id, tmp_stands.strata as strata, tmp_stands.initage as initage, st_translate( tmp_stands.boundary, 5

[postgis-users] Test the Future (GEOS 3.1RC2)

2009-02-02 Thread Paul Ramsey
PostGIS users, I've been blogging over the past several months about some big performance improvements in PostGIS, prepared geometries and cascaded union

[postgis-users] Re: Test the Future (GEOS 3.1RC2)

2009-02-03 Thread Paul Ramsey
P. On Mon, Feb 2, 2009 at 10:40 AM, Paul Ramsey wrote: > PostGIS users, > > I've been blogging over the past several months about some big > performance improvements in PostGIS, prepared geometries > <http://blog.cleverelephant.ca/2008/10/postgis-performance-prepared-geome

Re: [postgis-users] Re: Test the Future (GEOS 3.1RC2)

2009-02-03 Thread Paul Ramsey
No, that's as far as the thread got. For the short term, edit your makefile and remove the shp2pgsql-cli target and any other target that errors out on you. The old shp2pgsql is still there and will work fine. P. On Tue, Feb 3, 2009 at 10:58 AM, Bruce Rindahl wrote: > Sorry - I just found the th

Re: [postgis-users] Using a buffer with a LineString covering the whole world

2009-02-11 Thread Paul Ramsey
I find running screaming into the hills helps. Sorry, that's the use case (really large/long things) where there just is no proper solution except to start writing code that works on a sphere. Most things just aren't that long/large. But a few are. I would love it if someone/manyones would expres

Re: [postgis-users] Using a buffer with a LineString covering the whole world

2009-02-11 Thread Paul Ramsey
that's really not a good option. mercator does not preserve length/distance, so the buffer would end up larger at one end than the other, and furthermore would not be the expected width, since it would be built far from the latitude of true scale. P. On Wed, Feb 11, 2009 at 8:52 AM, Pavel Iacovle

Re: [postgis-users] still problems with Mapserver and MULTIPOLYGON layer

2009-02-11 Thread Paul Ramsey
I've brought it over to Mapserver trac. No luck reproducing it though. On Wed, Feb 11, 2009 at 9:48 AM, Mark Cave-Ayland wrote: > Ronald Woita wrote: > >> hi all, >> >> after installing PostGIS-1.3.5 (GEOS-3.0.0) in a PostgreSQL-8.3.6 >> server instance I tried to access on a MULTIPOLYGON Dataset

Re: [postgis-users] Using a buffer with a LineString covering the whole world

2009-02-11 Thread Paul Ramsey
might be > pretty close. If it was say Seattle Paris, that would be ok. > > Paul Ramsey wrote: >> that's really not a good option. mercator does not preserve >> length/distance, so the buffer would end up larger at one end than the >> other, and furthermore would not be t

Re: [postgis-users] design comments

2009-02-11 Thread Paul Ramsey
Pivoting a raster coverage into a table of points is an unbelievably bad idea. Don't do it. To the extent that the wktraster project starts to output some raster *processing* tools that lay on top of the storage abstraction, it will become useful for you. In the meanwhile, there is no reason you

Re: [postgis-users] Mapserver layer selects

2009-02-11 Thread Paul Ramsey
The first query refers to a table 'swip2' and the succeeding ones refer to 'swip3'. Also, the point returned in the third query is *different* that the one returned in the first. P On Wed, Feb 11, 2009 at 8:27 PM, Stephen Davies wrote: > Could somebody please help me explain the results below. >

Re: [postgis-users] points from st_distance

2009-02-12 Thread Paul Ramsey
Anything that you *can* do native in PostGIS, you *should* do. The benefit in speed and memory handling is high. If it was easy to implement all of GEOS functionality in pure C, we'd do that... but it's not :) It's easier to leverage the shared library. On the topic itself, this is doable, but not

Re: [postgis-users] SRID 27700 problem even with datum shift correction

2009-02-17 Thread Paul Ramsey
By saying @nadgrids=null you are implicitly saying "please don't try to do spheroid conversions for me". In cases where you are already in WGS84, this is good. In other cases, it isn't. The @nadgrids hack was invented by folks in North America, where their data is usually already in NAD83, which i

  1   2   3   4   5   6   7   8   9   10   >