gmane.comp.gis.postgis

Re: geom null

Chris Hermansen — 2008-11-28T17:18:46

Then if you can't add to a null geometry, you will probably be forced to digitize your places into another table and join the two after the data entry phase. Which means you will need identifiers on the geometry sub table that can be linked 1-1 with identifiers on the attribute sub table. You can use a correlated subquery to find rows in the geometry sub table that do not connect with the attribute sub table. ju wrote:

Re: geom null

Paul Ramsey — 2008-11-28T17:05:50

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 espaces-naturels.fr> wrote:

Re: geom null

ju — 2008-11-28T16:12:04

thanks but my table contains data for 1500 geographical place we are going to digitize some will be point, other line and other polygon so i can not make a false object and then modify it. Chris Hermansen a écrit :

Re: geom null

Chris Hermansen — 2008-11-28T15:32:43

Interesting question. What if, instead of null, you put in a simple piece of geometry that you can modify? (a point or a line or a 3-point linestring or a rectangle)? ju wrote:

Help with query

Mark — 2008-11-28T11:36:09

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

Re: Polygon validity with one common point

Nicolas Ribot — 2008-11-28T11:20:01

Sorry about this post, reading an old spec leads to this misinterpretation of the polygon validity. I attach anyway the picture describing the polygon. The common point in the given polygon is ("936.000000 211.000000") 2008/11/28 Nicolas Ribot gmail.com>: _______________________________________________ postgis-users mailing list postgis-users< at >postgis.refractions.net http://postgis.refractions.net/mailman/listinfo/postgis-users

Polygon validity with one common point

Nicolas Ribot — 2008-11-28T10:45:03

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 : 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 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

geom null

ju — 2008-11-28T08:33:11

Hi list, i have imported a table in postgres with only data fields. i add a geom field and i would like to digize object, for example with Qgis or GV sig. my table record is visible in qgis or gv sig but impossible to digitize... i think the trouble is not qgis or gvsig but a trouble of concept : is it possible to digitize a geom in postgis table if the record ( with geom null) allready exist ? thanks for helping,

RE: Help with query

Paragon Corporation — 2008-11-27T20:42:07

Mark, I guess my skills are a bit rusty in Mapserver. My first question is if they are both in PostgreSQL, why are you using OGR JOIN for this? Would seem better and more efficient to just do a JOIN in SQL which it seems you are doing already. You may actually want to separate into 2 layers one that has your grid and one that returns your species. Then what you have in TEMPLATE you put in HEADER instead And then have a TEMPLATE file with just the attributes as you want them laid out. Nowadays I just do Mapserver WFS/WMS calls for this kind of thing and parse the WFS XML in OpenLayers and overlay the WMS images and skip the whole template exercise altogether. Anyrate might be worthwhile to ask this question in Mapserver mailing list. Hope that helps, Regina -----Original Message----- From: postgis-users-bounces< at >postgis.refractions.net [mailto:postgis-users-bounces< at >postgis.refractions.net] On Behalf Of Mark Sent: Thursday, November 27, 2008 2:01 PM To: postgis-users< at >postgis.refractions.net Subject: [postgis-users] Help with query 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 _______________________________________________ postgis-users mailing list postgis-users< at >postgis.refractions.net http://postgis.refractions.net/mailman/listinfo/postgis-users

Re: Invalid geometry

Paul Ramsey — 2008-11-27T19:40:06

Try including one non-spatial field in the table you dump. P On Thu, Nov 27, 2008 at 11:33 AM, Travis Kirstine gmail.com> wrote:

Invalid geometry

Travis Kirstine — 2008-11-27T19:33:10

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.

Help with query

Mark — 2008-11-27T19:01:22

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

Re: maps google problem with SRID

lisek lichu — 2008-11-27T14:53:45

_______________________________________________ postgis-users mailing list postgis-users< at >postgis.refractions.net http://postgis.refractions.net/mailman/listinfo/postgis-users

Re: Integer (PostGIS) fields become real in QGIS!

Vito Meuli — 2008-11-27T08:59:45

Thank you Paul. They say it's an (old) OGR problem. I'll try and ask them a solution. Regards, Vito Alle 17:47, mercoledì 26 novembre 2008, Paul Ramsey ha scritto:

Re: maps google problem with SRID

Dirk Jesko — 2008-11-27T08:28:34

_______________________________________________ postgis-users mailing list postgis-users< at >postgis.refractions.net http://postgis.refractions.net/mailman/listinfo/postgis-users

RE: Arc catalog (Arc GIS) import

Paragon Corporation — 2008-11-27T02:32:09

see an ESRI frontend, they tend to look at you like you have three heads. Having 3 heads I imagine is more useful than having one assuming they are not all fighting with each other :) -----Original Message----- From: postgis-users-bounces< at >postgis.refractions.net [mailto:postgis-users-bounces< at >postgis.refractions.net]On Behalf Of John Abraham Sent: Tuesday, November 25, 2008 7:00 PM To: postgis-users< at >postgis.refractions.net Subject: [postgis-users] Arc catalog (Arc GIS) import A colleague told me that ArcGIS could import postgis layers directly, without using Ziggis. I've been googling around and playing with ArcMAP and ArcCatalog trying to figure it out, but no luck so far. I know I could use pgsql2shp. And I know I could use Ziggis commercial (2.x) or Ziggis open source (1.x). But for our needs a direct import menu item within ArcMAP or ArcCatalog would be best, as they would like to frequently view snapshots of the data. Any other suggestions also appreciated, thanks. -- John Abraham jabraham< at >ucalgary.ca _______________________________________________ postgis-users mailing list postgis-users< at >postgis.refractions.net http://postgis.refractions.net/mailman/listinfo/postgis-users _______________________________________________ postgis-users mailing list postgis-users< at >postgis.refractions.net http://postgis.refractions.net/mailman/listinfo/postgis-users

Re: maps google problem with SRID

lisek lichu — 2008-11-27T00:36:34

_______________________________________________ postgis-users mailing list postgis-users< at >postgis.refractions.net http://postgis.refractions.net/mailman/listinfo/postgis-users

maps google problem with SRID

lisek lichu — 2008-11-27T00:23:16

_______________________________________________ postgis-users mailing list postgis-users< at >postgis.refractions.net http://postgis.refractions.net/mailman/listinfo/postgis-users

Re: Union of MultiPolygons

Paul Ramsey — 2008-11-26T19:10:40

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. CREATE TABLE states_from_couties AS select st_multi(st_union(the_geom)) from county group by state; P. On Wed, Nov 26, 2008 at 11:02 AM, C. Y. W. dls.net> wrote:

Union of MultiPolygons

C. Y. W. — 2008-11-26T19:02:12

Hi All, I have a table with counties as Multipolygons. I need to define a new multipolygon with some of these counties. I tried to run this : INSERT INTO county (polygon_nm, the_geom) VALUES ('X', st_union((select the_geom from county where polygon_nm='LAKE'), (select the_geom from county where polygon_nm='COOK'))) And I got constraint "enforce_geotype_the_geom" violations. I guess st_union returns something that's not a Multipolygon. What should I do? Anyway to merge more than 2 Multipolygons at a time? Thanks, cyw

Transforming WG84 to google Mercator and back

Christoph Handel — 2008-11-26T17:38:58

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