http://blog.gmane.org/gmane.comp.gis.postgis hourly 1 1901-01-01T00:00+00:00 http://gmane.org/img/gmane-25t.png http://gmane.org http://comments.gmane.org/gmane.comp.gis.postgis/19775 I'm trying to fix some polygons-with-holes imported from shapefiles where the main polygon is not the first in the list, one of the hole polygons is, making them invalid. These are 3D multipolygons - imported from the SRTM SWBD shapefiles (with GDAL). (file w162n55n if anyone wants to try) The cleanGeometry funtion at: http://postgis.refractions.net/support/wiki/index.php?CleanPolygons isn't working, I get an error: =>select astext(cleangeometry(wkb_geometry)) from swbd where cell = 145018 and wb = 353; NOTICE: Hole lies outside shell at or near point -161.201 55.9992 0 CONTEXT: PL/pgSQL function "cleangeometry" line 15 at IF ERROR: Dimensions mismatch in lwcurve CONTEXT: PL/pgSQL function "cleangeometry" line 19 at assignment I tried forcing the geometry to 2D (st_force_2d()) and I get a different error: NOTICE: Hole lies outside shell at or near point -161.201 55.9992 CONTEXT: PL/pgSQL function "cleangeometry" line 15 at IF ERROR: Exception in LWGEOM2GEOS: curved geometry not supported. CONTEXT: PL/pgSQL function "cleangeometry" line 19 at assignment These shapefiles don't have curved geometry (I'm not sure it's supported in shapefiles). I wonder if there is a simpler operation that could reorder sub- polygons to put the one that contains all the others first? And would it be faster or slower? ----- William Kyngesburye <kyngchaos*at*kyngchaos*dot*com> http://www.kyngchaos.com/ "History is an illusion caused by the passage of time, and time is an illusion caused by the passage of history." - Hitchhiker's Guide to the Galaxy William Kyngesburye 2008-12-04T00:08:45 http://comments.gmane.org/gmane.comp.gis.postgis/19770 This problem occurs on versions... POSTGIS="1.3.1" GEOS="3.0.0rc4-CAPI-1.3.3" PROJ="Rel. 4.5.0, 22 Oct 2006" POSTGIS="1.3.2" GEOS="3.0.0-CAPI-1.4.1" PROJ="Rel. 4.6.0, 21 Dec 2007" The first query gives... select st_isvalid( st_GeomFromText('MULTILINESTRING((491181.9816 532686.6037,491209.6059 532664.6521),(491215.9577 532660.8051,491215.9577 532660.8051))')); NOTICE: IllegalArgumentException: Number of points must be 0 or >3 st_isvalid ------------ f ... while the second gives... select st_isvalid( st_GeomFromText('MULTILINESTRING((491181.9816 532686.6037,491209.6059 532664.6521),(491215.957 532660.8051,491215.9577 532660.8051))')); st_isvalid ------------ t ...and the only difference between the two is that the first point of the second line has an x value of 491215.9577 in the first query and 491215.957 in the second. In other words, the last seven is just stripped and it works. What gives here?! Any ideas? Ken Southerland Ken Southerland 2008-12-03T20:22:17 http://comments.gmane.org/gmane.comp.gis.postgis/19762 Hi all, I'm trying to figure out how to explode a set of polygons to display them with space between. In my case they are municipal boundary polys but could be anything. So far the only approach I can think of is using translate such that: select st_translate(st_translate(poly, -x(st_centroid(poly)), -y(st_centroid(poly))),x(p2) + (x(p2)- x(p1)), y(p2) + (y(p2) - y(p1))) from (select st_centroid(st_collect(geom))as p1 from mypolytable) as s1, (select st_centroid(geom) as p2, geom as poly from mypolytable) as s2 The problem with this is that small polys (Boroughs, etc) relative larger ones often don't get their own space until a large enough multiplier is applied, i.e. x(p2) + (x(p2)- x(p1))*5, y(p2) + (y(p2) - y(p1))*5 , but then there is too much space around the larger polys. Anyone have a suggestion or done this another way? Thanks. think i'm gonna splode..., Eric Randall, Eric 2008-12-02T22:27:31 http://comments.gmane.org/gmane.comp.gis.postgis/19759 Hi list! In my file .map I wrote something like select the_geom .......... (end_date - now() >= 10) as green and (end_date - now() < 0) as red .... from .... The idea is using different coloured markers for expired points or valid points After upgrading Postgres/Gis with last available versions, the query above doesn't work anymore because I cannot compare date types with integers ....... Any idea? Thanks in advance Vinc Vincenzo Patruno 2008-12-02T11:04:16 http://comments.gmane.org/gmane.comp.gis.postgis/19753 Many thanks Regina for your help, I wondered why I was using OGR - I guess I got confused as I was trying this earlier with joining a shapefile dbf to postgresql, which is where my problem arose! I will figure this out over the next few days. Thanks again Mark Mark 2008-11-28T11:36:09 http://comments.gmane.org/gmane.comp.gis.postgis/19751 Hi all, I'v got a validity issue with one polygon. As you can see in the attached picture, i've got a polygon with a linear ring that touch itself at one point : <gml:Polygon srsName="urn:ogc:def:derivedCRSType:OGC:1.0:image"> <gml:exterior> <gml:LinearRing> <gml:posList>33.000000 889.000000 51.000000 894.000000 64.000000 912.000000 80.000000 924.000000 119.000000 923.000000 125.000000 929.000000 157.000000 928.000000 160.000000 915.000000 179.000000 914.000000 186.000000 907.000000 186.000000 901.000000 189.000000 901.000000 190.000000 913.000000 211.000000 918.000000 211.000000 936.000000 255.500000 938.500000 300.000000 941.000000 295.000000 945.000000 249.000000 947.000000 223.000000 973.000000 210.000000 977.000000 211.000000 936.000000 184.000000 934.500000 157.000000 933.000000 156.000000 929.000000 125.250000 930.500000 94.500000 932.000000 63.750000 933.500000 33.000000 935.000000 33.000000 889.000000</gml:posList> </gml:LinearRing> </gml:exterior> </gml:Polygon> Or in WKT MULTIPOLYGON(((889.000000 33.000000,894.000000 51.000000,912.000000 64.000000,924.000000 80.000000,923.000000 119.000000,929.000000 125.000000,928.000000 157.000000,915.000000 160.000000,914.000000 179.000000,907.000000 186.000000,901.000000 186.000000,901.000000 189.000000,913.000000 190.000000,918.000000 211.000000,936.000000 211.000000,938.500000 255.500000,941.000000 300.000000,945.000000 295.000000,947.000000 249.000000,973.000000 223.000000,977.000000 210.000000,936.000000 211.000000,934.500000 184.000000,933.000000 157.000000,929.000000 156.000000,930.500000 125.250000,932.000000 94.500000,933.500000 63.750000,935.000000 33.000000,889.000000 33.000000))) This polygon seems to be not valid for postgis since it detects a Self-Intersection select polygon_id from polygon where isvalid(the_geom) = false; INFO: Ring Self-Intersection The OGC Simple Feature for SQL says : "A Polygon has no rings that cross. The rings in the boundary of a Polygon may intersect at a Point, but only as a tangent" If i well understand this definition, the polygon seems to be correct and should be valid for postgis... Am i wrong ? Regards Jerome Gasperi Nicolas Ribot 2008-11-28T10:45:03 http://comments.gmane.org/gmane.comp.gis.postgis/19747 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 function then exploding the multi features to single part features CREATE TABLE omnr_glcbp_geomunion AS SELECT upgis_cascadeunion(the_geom) as the_geom from omnr_glcbp_watersheds; CREATE TABLE omnr_glcbp_boundary AS SELECT ST_GeometryN(the_geom, generate_series(1, ST_NumGeometries(the_geom))) AS the_geom FROM omnr_glcbp_geomunion ; I then checked that that the geometry was valid and a multipolygon and everything looks fine. When the data is dumped to a shp file I get a error when opening the file with arcgis "Number of shapes does not match number of table records". Also when I check the shapefile with ogrinfo it does not report the geometry (returns Unknown). I thought that upgis_cascadeunion might be causing the problem however I had the same results when using the ST_Buffer(Collect(the_geom), 0) to dissolve the polygons. Travis Kirstine 2008-11-27T19:33:10 http://comments.gmane.org/gmane.comp.gis.postgis/19746 Hi All I am currently trying to create a query to pass through to mapserver but am having a problem. Here is my scenario: Two tables in postgresql A quarter degree grid with 273600 records (spatial geometry) with GiST index on the geometry and a btree index on FID A species table with 10 million records (btree index on FID1) contains species ids with related grid cell ids (essentially an intersect) In mapserver I have defined the layer thus LAYER CONNECTIONTYPE postgis NAME "SpcDensity" CONNECTION "user="" password="" dbname="" host=localhost" DATA "the_geom FROM (select spcdensity.fid_1, spcdensity.the_geom AS the_geom, gid AS gid, spcgrid.spcid from spcdensity, spcgrid where spcdensity.fid_1 = spcgrid.fid) as myquery using SRID=4326 using unique gid" STATUS OFF TYPE POLYGON TEMPLATE "templates/spc_header.html" CLASS NAME "density" OUTLINECOLOR 128 128 128 END JOIN NAME "species" CONNECTION "host=localhost port=5432 user="" password="" dbname=""" CONNECTIONTYPE ogr TABLE "spcgrid" FROM "fid_1" TO "fid" TYPE ONE-TO-ONE END END When I run this it works fine, click in a quarter degree grid cell and it returns a species ID using the one-to-one join but I need to get it to return the many species ids. Changing the join type to ONE-TO-MANY results in the error: loadJoin(): Premature End-of-File. One-to-many joins must define template and name properties If anyone has any suggestions as to how I can get this to work I would be most grateful. Thanks in advance Mark Mark 2008-11-27T19:01:22 http://comments.gmane.org/gmane.comp.gis.postgis/19740 _______________________________________________ postgis-users mailing list postgis-users< at >postgis.refractions.net http://postgis.refractions.net/mailman/listinfo/postgis-users lisek lichu 2008-11-27T00:23:16 http://comments.gmane.org/gmane.comp.gis.postgis/19737 Hi List, i have troubles transforming from WGS84 to Googles Mercator and back. On the commandline everthing is fine. WGS84 to google: cs2cs +proj=latlong +datum=WGS84 +to +proj=merc +a=6378137 +b=6378137 +lat_ts=0.0 +lon_0=0.0 +x_0=0.0 +y_0=0 +k=1.0 +units=m +nadgrids=< at >null +no_defs 56 46 6233891.485780349.22 0.00 Reverse cs2cs -I +proj=latlong +datum=WGS84 +to +proj=merc +a=6378137 +b=6378137 +lat_ts=0.0 +lon_0=0.0 +x_0=0.0 +y_0=0 +k=1.0 +units=m +nadgrids=< at >null +no_defs 6233891.485780349.22 56dE46dN 0.000 same in postgres forward select st_astext(st_transform(geomfromtext('POINT(56 46)',4326),900913)); st_astext ------------------------------------------ POINT(6233891.48442332 5780349.22025635) reverse select st_astext(st_transform(geomfromtext('POINT(6233891.48442332 5780349.22025635)',900913),4326)); WARNING: transform: -38 (failed to load NAD27-83 correction file) st_astext -------------------------------------------------- POINT(8.77999328016943e-06 7.21213733809987e-06) (1 row) There is a posting in google (german): http://groups.google.com/group/postgis-users/browse_thread/thread/158c4ea06d710010?hl=de&ie=UTF-8&q=ST_Transform+failed+to+load+NAD27-83+correction+file#e7adf1f4291550cd Explaining that some files are missing on ubuntu. But as transformation on the commandline are working that shouldn't be the case. all my packages are stock ubuntu 8.04. proj: 4.6.0 postgresql: 8.3 postgis: 1.3.3 Did i miss something? Any hints? Greetings Christoph Christoph Handel 2008-11-26T17:38:58 http://comments.gmane.org/gmane.comp.gis.postgis/19723 _______________________________________________ postgis-users mailing list postgis-users< at >postgis.refractions.net http://postgis.refractions.net/mailman/listinfo/postgis-users naina susen 2008-11-26T09:01:44 http://comments.gmane.org/gmane.comp.gis.postgis/19722 Hello, with the new release postgis will support Compound Curves in a Curve polygon? I realy need this feature. Thanks mario< at >geosar.ch 2008-11-26T07:43:41 http://comments.gmane.org/gmane.comp.gis.postgis/19708 Hi all, I am developing an application that uses google maps as a front-end and postgis in the backend. Currently, we are using google mercator projection (SRID = 4326) for all geometries in the database. But, since we perform lots of geometric operations with these objects within the database, I am concerned about performance issues. Wouldn't it be better if instead of a mercator projection, we stored all information in metric-based values? If so, is there any metric-based projection that works in any part of the world? Thanks for any help. Davi Pires Davi Pires 2008-11-25T12:58:15 http://comments.gmane.org/gmane.comp.gis.postgis/19702 After a lengthy development and testing process, it is with some trepidation that we bring PostGIS 1.3.4 out into a harsh and unforgiving world. http://postgis.refractions.net/download/postgis-1.3.4.tar.gz - New Features - Add new ST_AsGeoJSON() function to allow conversion of geometries 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 comments to be installed for PostGIS functions. This attaches a description and associated parameter information to each function so that this information is available directly within psql/pgAdmin (Regina Obe) - Enhancements and Bug Fixes - General documentation improvements (Regina Obe, Kevin Neufield) - Improve PiP code by removing substantial memory leaks and adding multi-polygon support (Paul Ramsey) - Improve GiST consistent performance (Paul Ramsey) - GBT#20: Improved algorithm for ST_DWithin (Paul Ramsey) - GBT#21: locate_along_measure: wrong values, invalid data (Mark Cave-Ayland) - GBT#23: wrong proj4text in spatial_ref_sys for SRID 31300 and 31370 (Paul Ramsey) - GBT#43: Heap over-read in compute_geometry_stats() (Landon Fuller) - GBT#44: ST_Relate(geom,geom,pattern) is case sensitive (Paul Ramsey) - GBT#45: bug in postgis_restore.pl - not passing options to createdb (Paul Ramsey) - GBT#58: bounding box of circular string is incorrect (Mark Leslie) - GBT#65: ST_AsGML kills the backend when fed a CIRCULAR STRING (Mark Cave-Ayland) - GBT#66: ST_Dump kills backend when fed CIRCULAR STRING (Mark Cave-Ayland) - GBT#68: ST_Shift_Longitude doesn't work with MULTIPOINT (Mark Cave-Ayland) - GBT#69: ST_Translate crashes when fed circular string (Mark Cave-Ayland) - GBT#70: ST_NPoints Crashes with Curves (Mark Cave-Ayland) Paul Ramsey 2008-11-25T06:30:36 http://comments.gmane.org/gmane.comp.gis.postgis/19700 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) 3.3.1 (Mandrake Linux 9.2 3.3.1-2mdk) Postgis 1.3.4rc3 - 2008-11-25 03:24:13 GEOS: 3.0.0rc4-CAPI-1.3.3 PROJ: Rel. 4.5.0, 22 Oct 2006 Running tests . . regress_lrs. failed (diff expected obtained: /var/lib/pgsql/tmp/test_17_diff) Diff file has: 17c17 < LINEZM_6|POINT(9.5 0.5 0.5 2) --- 60c60 < LINEZM_6|POINT(9.5 0.5 0.5 2) --- All else seems AOK. Cheers, Stephen Davies Stephen Davies 2008-11-25T05:09:33 http://comments.gmane.org/gmane.comp.gis.postgis/19697 Hi guys :) -- a *special* wink to David Techer ;-) I'm having a bit of trouble figuring out what can I do with the wrapper functions of pgRouting... (blush) The thing is: I have everything perfectly in place, pgRouting working as it's supposed to, but! I have a dataset that describes the road type. I want to filter the results provided by pgRouting in a way that it *doesn't* return return pedestrian segments (btw: road_type=22). So, as to the query itself (basic) 1st try: SELECT gid, AsBinary(the_geom) AS the_geom FROM dijkstra_sp_directed('pt_madeira_roads', 11780, 424, true, true); Works like a charm except that it returns pedestrian segments (of course). 2nd try: SELECT rt.gid, AsText(rt.the_geom) AS wkt, length(rt.the_geom) AS length, roads.gid FROM my_roads as roads, (SELECT gid, the_geom FROM dijkstra_sp_directed( 'my_roads', 11780, 424, true,true) ) as rt WHERE roads.road_type<>22 and roads.gid=rt.gid; Leaves me with *holes* in the route. All of this is, of course, a consequence of my own stupidity and late hours but I'd surely appreciate if some brave soul could look into it. ;-) Already thankful in advance, Pedro Doria Meunier. Pedro Doria Meunier 2008-11-25T02:16:48 http://comments.gmane.org/gmane.comp.gis.postgis/19691 _______________________________________________ postgis-users mailing list postgis-users< at >postgis.refractions.net http://postgis.refractions.net/mailman/listinfo/postgis-users John Zhang 2008-11-25T00:14:25 http://comments.gmane.org/gmane.comp.gis.postgis/19674 _______________________________________________ postgis-users mailing list postgis-users< at >postgis.refractions.net http://postgis.refractions.net/mailman/listinfo/postgis-users John Zhang 2008-11-22T20:42:35 http://comments.gmane.org/gmane.comp.gis.postgis/19667 _______________________________________________ postgis-users mailing list postgis-users< at >postgis.refractions.net http://postgis.refractions.net/mailman/listinfo/postgis-users John Zhang 2008-11-22T02:22:31 http://comments.gmane.org/gmane.comp.gis.postgis/19666 hi list! does postgis provide a method to reduce/simplify a linestring to a specific number of points? i'm aware of st_simplify() but i have to set an lower limit. thanks, //richard richard -rw- weinberger 2008-11-22T01:34:52 http://comments.gmane.org/gmane.comp.gis.postgis/19662 I'm having some trouble loading PostGIS data into MapWindow GIS http://www.mapwindow.org/phorum/read.php?4,9561,10771#msg-10771 Wondering if anyone here may be able to lend a hand to me and others on this thread. I would normally use QGIS / Grass here but making edits to the attribute table is a bit smoother in MapWindow. Thank you! ~ Joe Joe Larson 2008-11-21T22:06:51 Search the mailing list at Gmane query http://search.gmane.org/?group=$group=gmane.comp.gis.postgis