<!DOCTYPE html PUBLIC "-//W3C//DTD HTML 4.01 Transitional//EN">
<html>
<head>
<meta content="text/html;charset=ISO-8859-1" http-equiv="Content-Type">
</head>
<body bgcolor="#ffffff" text="#000000">
Hello!<br>
Kevin, you are growing to be my personal hero. I am getting a good
understanding here. I used to work with Spatialware from MapInfo a lot,
so I know the matter well. I had figured out that the geometry_columns
table would be something like mapinfo's mapcatalog.<br>
<br>
I would love to write an article on the proces I now go to for the
official user documentation once I feel I controle the proces. The
article would be called something like:<br>
"How to - copy spatial tables with options in postGIS"<br>
<br>
Kevin Neufeld schreef:
<blockquote cite="mid:46A101BA.4020301@refractions.net" type="cite">
<meta content="text/html;charset=ISO-8859-1" http-equiv="Content-Type">
Ah, ok. Getting a clearer picture.<br>
<br>
It looks like your production table is currently empty which is why it
returned nothing. <br>
</blockquote>
That is right!<br>
<blockquote cite="mid:46A101BA.4020301@refractions.net" type="cite"><br>
In response to your comment regarding the geometry_columns table...<br>
<br>
PostGIS is an add-on to PostgreSQL that simply adds a few new datatypes
(ie. geometry, box3d), a few functions (ie. intersects(...),
distance(...)) and a couple of tables (geometry_columns and
spatial_ref_sys) to a database instance. <br>
<br>
In postgres, it's entirely possible to simply create a table with a
geometry column by simply typing:<tt><br>
CREATE TABLE public.tmp (id integer, the_geom geometry);</tt><br>
<br>
However, creating a table this way would permit someone to insert
geometries of varying types (ie. POINT, LINESTRING,...), of varying
dimensions (2, 3 or 4D), and of varying projections (different SRIDs).
What would be ideal is to add constraints to your table that would
ensure that geometries are of the same type.<br>
<br>
One of the many functions added by PostGIS is addgeometrycolumn(...).
Use this function whenever you want to add a geometry column to a
table. So our previous example would be:<br>
<tt>CREATE TABLE public.tmp (id integer);<br>
SELECT addgeometrycolumn('public','tmp','the_geom',4326,'LINESTRING',2);</tt><br>
</blockquote>
That is clear! Thank you, I will drop my hand made tables and recreate
them in this manner. It seems that OGR does it correct during import,
so the staging table can remain as is.<br>
<blockquote cite="mid:46A101BA.4020301@refractions.net" type="cite"><br>
This function simply adds a column called '<tt>the_geom</tt>' of type '<tt>geometry</tt>'
to the <tt>public.tmp</tt> table, adds three constraints to the table (<i>enforce_dims_the_geom,
enforce_geotype_the_geom, </i>and <i>enforce_srid_the_geom</i>), and
inserts a row into the geometry_colums table.<br>
This ensures that the public.tmp table will only have geometries of
type 'LINESTRING', are of 2 dimensions, and have an SRID of 4326.<br>
</blockquote>
Perfect!<br>
<blockquote cite="mid:46A101BA.4020301@refractions.net" type="cite"><br>
See (<a moz-do-not-send="true" class="moz-txt-link-freetext"
href="http://postgis.refractions.net/docs/ch06.html#id2997535">http://postgis.refractions.net/docs/ch06.html#id2997535</a>)<br>
<br>
Whenever you drop a table, you must first use the
dropgeometrycolumn(...) function that will automatically remove the
necessary row from the geometry_columns table and drop the geometry
column. Then you can drop the table.<br>
</blockquote>
That's clear and that was a function I was looking for. I will check if
it has been discussed to create a trigger for this, because I think
that would make life a bit easier. <br>
<blockquote cite="mid:46A101BA.4020301@refractions.net" type="cite"><br>
If you drop a table without calling dropgeometrycolumn first, you will
have to manually delete the corresponding row in the geometry_columns
table. Not a big deal, but a pain.<br>
</blockquote>
clear<br>
<blockquote cite="mid:46A101BA.4020301@refractions.net" type="cite"><br>
There is nothing special at all about the geometry_columns table. It's
just a simple table that stores metadata about the tables with geometry
columns in the database, and as you've already discovered, it's easy to
get out of date. It's really only there to be OGC compliant. Some
applications, like mapserver, use this table to quickly find all the
registered geometry columns in a database ... although a complete list
can be obtained by looking at the system tables.<br>
</blockquote>
I understand. On the things that follow below, I have to get back on
that next week. It seems the person who contracted me closes it's VPN
connection in the weekend, so I cannot access the machine I was working
on! Hahaha, that is a way to force people to take weekends off.<br>
<br>
<br>
<blockquote cite="mid:46A101BA.4020301@refractions.net" type="cite"><br>
Hope that makes things clearer.<br>
<br>
<br>
Anywho, back to your problem, it sounds like there is an SRID
constraint on your production table that only allows geometries with a
certain SRID to be inserted into the table.<br>
<br>
You mentioned that you were using pgAdmin, right? When you select your
table in pgAdmin, what is the exact table definition? <br>
If you use postgres's command-line terminal, psql, you could simply
issue the following command: "<tt>\d infra_line2</tt>".<br>
<br>
This will show you the table definition and all the constraints on the
table. <br>
<br>
I suppose you could also view the table's constraints which this query:
(I think i got this right)<br>
<tt>SELECT c.relname as table_name, conname as constraint_name,
pg_catalog.pg_get_constraintdef(r.oid, true) as constraint_def<br>
FROM pg_catalog.pg_constraint r, pg_class c<br>
WHERE r.conrelid = c.oid<br>
AND r.contype = 'c'<br>
AND c.relname = 'infra_line2</tt><tt>';</tt><br>
<br>
What you need to do is wrap your geometries with a
setsrid(wkb_geometry, XXXX) when you insert them into the infra_line2
table, where XXXX is the <i>enforce_srid_the_geom</i> constraint you
have on that table.<br>
<br>
Clear as mud?<br>
<br>
Have a great weekend!<br>
-- Kevin<br>
<br>
<br>
Milo van der Linden wrote:
<blockquote cite="mid:46A0EF63.5060600@zeelandnet.nl" type="cite">
<meta content="text/html;charset=ISO-8859-1"
http-equiv="Content-Type">
<title></title>
Thank you Kevin, One step closer to goal!<br>
<br>
I ran the 2 select statements you mentioned below and found out the
following;<br>
- For infra_line2, it returns nothing!<br>
- bm_staging returns the right value: 32767. <br>
<br>
I then tried the following; <br>
UPDATE infra_line2 SET the_geom = SetSRID(the_geom, (SELECT srid FROM
geometry_columns WHERE f_table_name = 'infra_line2'));<br>
SELECT distinct( getsrid( the_geom ) ) FROM infra_line2;<br>
SELECT distinct( getsrid( wkb_geometry ) ) FROM bm_staging;<br>
<br>
But still infra_line2 returns nothing.<br>
<br>
I noticed something else; when I drop a table, the entry in the
geometry_columns tables remains. And opposite, when I create a table,
no entry is set in the geometry_columns.<br>
<br>
I am doing a project initiated by someone else, could it be that I have
to reinstall? Is there a way to check what version of postGIS is the
installed one?<br>
<br>
I also read on the OGR forum that when there is just a tiny mismatch
between the way EPSG:4326 is described in the conversion script and the
way it is implemented in postGIS, OGR creates a new SRID during import,
that's why my coordsys is in 32767. For this I am going to use
SetSRID(the_geom, 4326)<br>
<br>
<br>
Help apreciated! Special thanks to you, Kevin.<br>
<br>
<br>
Kevin Neufeld schreef:
<blockquote cite="mid:46A0E620.6090607@refractions.net" type="cite">
<meta content="text/html;charset=ISO-8859-1"
http-equiv="Content-Type">
<title></title>
The error is being generated because the geometries you are inserting
into your production table have a different SRID than the constraint
you have on your table. <br>
<br>
You can either transform the geometries to the SRID of your table as
you insert them, or you can simply set the SRID of your inserted
geometries to be the same. <br>
<br>
For your insert query to work, these numbers should both be the same:<br>
<tt>SELECT distinct( getsrid( the_geom ) ) FROM
public.infra_line2;<font color="#006600"><br>
</font>SELECT distinct( getsrid( wkb_geometry ) ) FROM bm_staging;</tt><br>
<tt><font color="#006600"><br>
</font></tt>Cheers,<br>
-- Kevin<br>
<br>
Milo van der Linden wrote:
<blockquote cite="mid:46A0E0C1.40107@zeelandnet.nl" type="cite">
<meta content="text/html;charset=ISO-8859-1"
http-equiv="Content-Type">
<title></title>
<meta content="text/html;charset=ISO-8859-1"
http-equiv="Content-Type">
<title></title>
Hello Kevin, thank you for the response!<br>
<br>
Kevin Neufeld schreef:
<blockquote cite="mid:46A0D632.9040408@refractions.net"
type="cite">
<meta content="text/html;charset=ISO-8859-1"
http-equiv="Content-Type">
<br>
Milo van der Linden wrote:
<blockquote cite="mid:46A0B40C.6000704@zeelandnet.nl"
type="cite"><br>
The original data is in shape files. With ogr2ogr I do a conversion
from the Dutch coordinate system to WGS84 and import it in a temp_table
in PostGIS.<br>
</blockquote>
Did you know that PostGIS has a shape loader? You may find loading
shape files easier with the shp2pgsql utility. Also, PostGIS can
perform transformations on geometries ... you may not even need ogr.
Just a thought.<br>
</blockquote>
I knew that, but I am very familiar with ogr and not with shp2pgsql.
Besides that, it is fired from a php script that is currently well
under control. Also, shp2pgsql is from the same development team as
ogr, they probably use the same code base.<br>
<blockquote cite="mid:46A0D632.9040408@refractions.net"
type="cite">
<blockquote cite="mid:46A0B40C.6000704@zeelandnet.nl"
type="cite"><br>
When I execute the statement in pgAdmin, I get an error message from
postGIS:<br>
<br>
<font color="#ff0000">ERROR: new row for relation
"infra_line2"
violates check constraint "enforce_srid_the_geom"</font><br>
<br>
I tried to add USING SRID=32767 to the insert statement, but it fails
with another error message<br>
</blockquote>
It sounds like the geometries in your temp table either doesn't have,
or has a different SRID than your production table. Is the data in your
staging table supposed to be the same projection as your production
table? (ogr possibly didn't add an SRID to your imported geometries and
they may be defaulting to -1). If this is the case, then you simply
need to wrap your geometries with a setsrid function call, ie
"SetSRID(the_geom, 32767)", assuming of course that 32767 is the SRID
of your production table.<br>
</blockquote>
When I check my geometry meta table, the SRID for the staging table is
correct, even though it should be 4326 (wgs84), during the import ogr
decided to create a new srid that is a perfect copy of the standard
wgs84 SRID.<br>
<blockquote cite="mid:46A0D632.9040408@refractions.net"
type="cite">
<blockquote><tt><font color="#006600">INSERT
INTO public.infra_line2 </font></tt><br>
<tt><font color="#006600">
(infraline_type, </font></tt><br>
<tt><font color="#006600">
infraline_polution, </font></tt><br>
<tt><font color="#006600">
infraline_comment, </font></tt><br>
<tt><font color="#006600">
infraline_startdate, </font></tt><br>
<tt><font color="#006600">
infraline_enddate, </font></tt><br>
<tt><font color="#006600">
infraline_top, </font></tt><br>
<tt><font color="#006600">
infraline_base, </font></tt><br>
<tt><font color="#006600">
the_geom, </font></tt><br>
<tt><font color="#006600">
idproject, </font></tt><br>
<tt><font color="#006600">
idsh0, </font></tt><br>
<tt><font color="#006600">
idshp) </font></tt><br>
<tt><font color="#006600">SELECT </font></tt><br>
<tt><font color="#006600">'type' as infraline_type,</font></tt><br>
<tt><font color="#006600">CASE WHEN ver_inr='0' THEN FALSE
WHEN
ver_inr='1' THEN TRUE ELSE FALSE
END as infraline_polution,</font></tt><br>
<tt><font color="#006600">opmerking as infraline_comment,</font></tt><br>
<tt><font color="#006600">d_start as infraline_startdate, </font></tt><br>
<tt><font color="#006600">d_stop as infraline_stopdate, </font></tt><br>
<tt><font color="#006600">hoogte_top as infraline_top, </font></tt><br>
<tt><font color="#006600">hoogte_bas as infraline_base, </font></tt><br>
<tt><font color="#006600">setsrid(wkb_geometry, 32767) as
the_geom, </font></tt><br>
<tt><font color="#006600">1 as idproject,</font></tt><br>
<tt><font color="#006600">6 as indsh0,</font></tt><br>
<tt><font color="#006600">id as idshp</font></tt><br>
<tt><font color="#006600">FROM bm_staging;</font></tt><br>
</blockquote>
Hope this helps,<br>
</blockquote>
Thank you for giving it a try, but no, it doesn't. I need to know why
the error is generated. Perhaps some one from postGIS development can
tell me?<br>
<blockquote cite="mid:46A0D632.9040408@refractions.net"
type="cite">--
Kevin<br>
<br>
<br>
--------<br>
Kevin Neufeld<br>
Refractions Research Inc.<br>
<a moz-do-not-send="true" class="moz-txt-link-abbreviated"
href="mailto:kneufeld@refractions.net">kneufeld@refractions.net</a><br>
(250) 383-3022<br>
<br>
</blockquote>
<br>
<br>
<div class="moz-signature">-- <br>
<title>3DSite</title>
<meta http-equiv="Content-Type" content="text/html; ">
<style type="text/css">
<!--
body {
background-color: #FFFFFF;
margin-left: 5px;
margin-top: 5px;
margin-right: 5px;
margin-bottom: 5px;
}
body,td,th {
font-family: Arial, Helvetica, sans-serif;
font-size: 12px;
color: #333333;
}
.style4 {font-size: 9px; }
.style5 {font-size: 9px; color: #CCCCCC; }
-->
</style>
<meta content="MSHTML 6.00.2900.2912" name="GENERATOR">
<br>
<table border="0" cellpadding="0" cellspacing="0" width="400">
<tbody>
<tr>
<td rowspan="3" align="left" height="0" valign="bottom"
width="15"><br>
</td>
<td colspan="2" align="left" height="78" valign="top"
width="0">
<p class="style4">Milo van der Linden <br>
<a moz-do-not-send="true"
href="mailto:mlinden@zeelandnet.nl">mlinden@zeelandnet.nl</a><br>
<a moz-do-not-send="true"
href="mailto:milovanderlinden@gmail.com">milovanderlinden@gmail.com</a><br>
<a moz-do-not-send="true" href="mailto:milo@3dsite.nl">milo@3dsite.nl</a><br>
<a moz-do-not-send="true" href="http://www.3dsite.nl">http://www.3dsite.nl</a><br>
</p>
</td>
<td rowspan="3" align="left" height="0" valign="top"
width="15"> </td>
<td valign="top" width="300">
<p class="style5"><span lang="NL">De informatie in dit
bericht
reflecteerd mijn persoonlijke mening en niet die van een bedrijf of
instantie. Aan de informatie kunnen geen rechten worden ontleend.
Indien dit bericht onderdeel is van een forum, mailing-list of
community dan gelden automatisch de bijbehorende voorwaarden.</span></p>
</td>
</tr>
</tbody>
</table>
</div>
<br>
<div class="moz-signature">-- <br>
<title>3DSite</title>
<meta http-equiv="Content-Type" content="text/html; ">
<style type="text/css">
<!--
body {
background-color: #FFFFFF;
margin-left: 5px;
margin-top: 5px;
margin-right: 5px;
margin-bottom: 5px;
}
body,td,th {
font-family: Arial, Helvetica, sans-serif;
font-size: 12px;
color: #333333;
}
.style4 {font-size: 9px; }
.style5 {font-size: 9px; color: #CCCCCC; }
-->
</style>
<meta content="MSHTML 6.00.2900.2912" name="GENERATOR">
<br>
<table border="0" cellpadding="0" cellspacing="0" width="400">
<tbody>
<tr>
<td rowspan="3" align="left" height="0" valign="bottom"
width="15"><br>
</td>
<td colspan="2" align="left" height="78" valign="top"
width="0">
<p class="style4">Milo van der Linden <br>
<a moz-do-not-send="true"
href="mailto:mlinden@zeelandnet.nl">mlinden@zeelandnet.nl</a><br>
<a moz-do-not-send="true"
href="mailto:milovanderlinden@gmail.com">milovanderlinden@gmail.com</a><br>
<a moz-do-not-send="true" href="mailto:milo@3dsite.nl">milo@3dsite.nl</a><br>
<a moz-do-not-send="true" href="http://www.3dsite.nl">http://www.3dsite.nl</a><br>
</p>
</td>
<td rowspan="3" align="left" height="0" valign="top"
width="15"> </td>
<td valign="top" width="300">
<p class="style5"><span lang="NL">De informatie in dit
bericht
reflecteerd mijn persoonlijke mening en niet die van een bedrijf of
instantie. Aan de informatie kunnen geen rechten worden ontleend.
Indien dit bericht onderdeel is van een forum, mailing-list of
community dan gelden automatisch de bijbehorende voorwaarden.</span></p>
</td>
</tr>
</tbody>
</table>
</div>
<pre wrap=""><hr size="4" width="90%">
_______________________________________________
postgis-users mailing list
<a moz-do-not-send="true" class="moz-txt-link-abbreviated"
href="mailto:postgis-users@postgis.refractions.net">postgis-users@postgis.refractions.net</a>
<a moz-do-not-send="true" class="moz-txt-link-freetext"
href="http://postgis.refractions.net/mailman/listinfo/postgis-users">http://postgis.refractions.net/mailman/listinfo/postgis-users</a>
</pre>
</blockquote>
</blockquote>
<br>
<br>
<div class="moz-signature">-- <br>
<title>3DSite</title>
<meta http-equiv="Content-Type" content="text/html; ">
<style type="text/css">
<!--
body {
background-color: #FFFFFF;
margin-left: 5px;
margin-top: 5px;
margin-right: 5px;
margin-bottom: 5px;
}
body,td,th {
font-family: Arial, Helvetica, sans-serif;
font-size: 12px;
color: #333333;
}
.style4 {font-size: 9px; }
.style5 {font-size: 9px; color: #CCCCCC; }
-->
</style>
<meta content="MSHTML 6.00.2900.2912" name="GENERATOR">
<br>
<table border="0" cellpadding="0" cellspacing="0" width="400">
<tbody>
<tr>
<td rowspan="3" align="left" height="0" valign="bottom"
width="15"><br>
</td>
<td colspan="2" align="left" height="78" valign="top"
width="0">
<p class="style4">Milo van der Linden <br>
<a moz-do-not-send="true" href="mailto:mlinden@zeelandnet.nl">mlinden@zeelandnet.nl</a><br>
<a moz-do-not-send="true"
href="mailto:milovanderlinden@gmail.com">milovanderlinden@gmail.com</a><br>
<a moz-do-not-send="true" href="mailto:milo@3dsite.nl">milo@3dsite.nl</a><br>
<a moz-do-not-send="true" href="http://www.3dsite.nl">http://www.3dsite.nl</a><br>
</p>
</td>
<td rowspan="3" align="left" height="0" valign="top"
width="15"> </td>
<td valign="top" width="300">
<p class="style5"><span lang="NL">De informatie in dit
bericht
reflecteerd mijn persoonlijke mening en niet die van een bedrijf of
instantie. Aan de informatie kunnen geen rechten worden ontleend.
Indien dit bericht onderdeel is van een forum, mailing-list of
community dan gelden automatisch de bijbehorende voorwaarden.</span></p>
</td>
</tr>
</tbody>
</table>
</div>
</blockquote>
</blockquote>
<br>
<br>
<div class="moz-signature">-- <br>
<title>3DSite</title>
<meta http-equiv="Content-Type" content="text/html; ">
<style type="text/css">
<!--
body {
background-color: #FFFFFF;
margin-left: 5px;
margin-top: 5px;
margin-right: 5px;
margin-bottom: 5px;
}
body,td,th {
font-family: Arial, Helvetica, sans-serif;
font-size: 12px;
color: #333333;
}
.style4 {font-size: 9px; }
.style5 {font-size: 9px; color: #CCCCCC; }
-->
</style>
<meta content="MSHTML 6.00.2900.2912" name="GENERATOR">
<br>
<table border="0" cellpadding="0" cellspacing="0" width="400">
<tbody>
<tr>
<td rowspan="3" align="left" height="0" valign="bottom" width="15"><br>
</td>
<td colspan="2" align="left" height="78" valign="top" width="0">
<p class="style4">Milo van der Linden
<br>
<a href="mailto:mlinden@zeelandnet.nl">mlinden@zeelandnet.nl</a><br>
<a href="mailto:milovanderlinden@gmail.com">milovanderlinden@gmail.com</a><br>
<a href="mailto:milo@3dsite.nl">milo@3dsite.nl</a><br>
<a href="http://www.3dsite.nl">http://www.3dsite.nl</a><br>
</p>
</td>
<td rowspan="3" align="left" height="0" valign="top" width="15"> </td>
<td valign="top" width="300">
<p class="style5"><span lang="NL">De informatie in dit bericht
reflecteerd mijn persoonlijke mening en niet die van een bedrijf of
instantie. Aan de informatie kunnen geen rechten worden ontleend.
Indien dit bericht onderdeel is van een forum, mailing-list of
community dan gelden automatisch de bijbehorende voorwaarden.</span></p>
</td>
</tr>
</tbody>
</table>
</div>
</body>
</html>