I will be needing the database to integrate data from the other tables. I want to create a field containing all the bus lines that serve this stop, but I was able to do that before performing joins. That's why I was trying to do all of it in PostGIS. I have not been very succesful with GDAL in Python, when I was trying to do this before. <br>
I started with fields named x and y, but since the target is to have lon and lat in there, I renamed them. COPY won't allow me to import a csv file if the number of columns in the DB doesn't correspond to the number of columns in the csv file.<br>
<br>Here is a sample line from the csv file:<br><br>5349;308999;Dorp;Dorp;Moerzeke;Hamme;134858;194785;false;true<br><br>134858;194785 is what I also find in my lon and lat columns in the DB after COPY ran successfully. I know the trigger was invoked, as I had to debug it before it would work...<br>
<br>I had expected something like 2.34 and 50.1 in those columns. (those are not the real corresponding coordinates).<br><br>I found an example doing the reverse, but they had to use string concatenation to  achieve what they wanted to do. Should I also let it create a string and then parse that? That seems kind of backwards to me...<br>
<br>Jo<br><br><br>
<div class="gmail_extra"><br><br><div class="gmail_quote">2012/11/14 Nicolas Ribot <span dir="ltr"><<a href="mailto:nicolas.ribot@gmail.com" target="_blank">nicolas.ribot@gmail.com</a>></span><br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Hi,<div><br></div><div>If the lat and lon fields of the stops table are already in WGS84, as their names suggest (latitude and longitude), then you don't have to transform them.</div><div><br></div><div>Could you give us a sample of the table, to see the coordinates in lat and lon fields ?</div>

<div>In a first step, I would add 2 new columns to the stops table storing transformed coordinates (WGS84) and keep columns storing Lambert72 values, to check.</div><div><br></div><div>Otherwise, if the sole purpose is to produce an XML file, Paolo's approach would be the most straightforward: use ogr to both reproject the data and change their format.</div>
<span class="HOEnZb"><font color="#888888">
<div><br></div><div>Nicolas</div></font></span><div class="HOEnZb"><div class="h5"><div><br></div><div><br></div><div><br></div><div class="gmail_extra"><br><br><div class="gmail_quote">On 14 November 2012 01:36, Jo <span dir="ltr"><<a href="mailto:winfixit@gmail.com" target="_blank">winfixit@gmail.com</a>></span> wrote:<br>

<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Hi Nicolas and Paolo,<br><br>This is what I have been able to come up with:<br><br><blockquote><span></span>DROP EXTENSION postgis;<br>

CREATE EXTENSION postgis;<br><br>DROP TABLE IF EXISTS stops;<br>
<br>CREATE TABLE stops<br>(<br>  stopid bigint NOT NULL,<br>  stopidentifier bigint,<br>  description text,<br>  street text,<br>  municipality text,<br>  parentmunicipality text,<br>  lat double precision,<br>  lon double precision,<br>


  stopisaccessible boolean,<br>  stopispublic boolean,<br>  PRIMARY KEY (stopid)<br>)<br>WITH (<br>  OIDS=FALSE<br>);<br>ALTER TABLE stops<br>  OWNER TO postgres;<br><br>CREATE OR REPLACE FUNCTION addwgs84() RETURNS trigger AS<br>


  $body$<br>  DECLARE<br>    node geometry;<br>  BEGIN<br>    node := st_transform(st_setSRID(st_makePoint(NEW.lat, NEW.lon), 31370),4326);<br>    NEW.lat := st_x(node);<br>    NEW.lon := st_y(node);   <br>    RETURN NEW;<br>


  END;<br>  $body$<br>  LANGUAGE 'plpgsql';<br><br>CREATE TRIGGER wgs84<br>  AFTER INSERT OR UPDATE ON stops<br>  FOR EACH ROW EXECUTE PROCEDURE addwgs84();<br><br>COPY stops FROM 'C:/stops.csv' DELIMITERS ';' CSV HEADER;<span></span><br>


</blockquote><br><span></span><span>Unfortunately, this doesn't result in WGS84 lat and lon, but gives me lambert72 once again. I must be missing something crucial somewhere. I don't need the geometry in the table, so I took it out again. I'm not going to use it to calculate distances. This is only an intermediate step in converting all of it to an XML-format suitable for Openstreetmap.org.<br>


<br>Thanks for your help. Hopefully I'm near to the goal I wanted to reach, now and you can nudge me a bit further in the right direction.<span><font color="#888888"><br><br>Jo<br></font></span></span><div>
<div><div class="gmail_extra"><br><br><div class="gmail_quote">
2012/11/13 Nicolas Ribot <span dir="ltr"><<a href="mailto:nicolas.ribot@gmail.com" target="_blank">nicolas.ribot@gmail.com</a>></span><br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">


Hi,<div><br></div><div>Use the  COPY command to load CSV to a table.<br><br></div><div>The rest depends if you want to display/use the geometry, or you're only interested in the coordinates:</div><div>• Add a geometry column</div>



<div>• st_makePoint to build a point from coordinates</div><div>• A trigger to create the WGS84 point from coordinates with st_transform: st_transform(st_setSRID(st_makePoint(<span style="font-family:arial,sans-serif;font-size:13px">coordinateX, </span><span style="font-family:arial,sans-serif;font-size:13px">coordinateY), 31300), 4326)</span></div>



<div><span style="font-family:arial,sans-serif;font-size:13px"><br></span></div><div><span style="font-family:arial,sans-serif;font-size:13px">Nicolas</span></div><div class="gmail_extra"><br><br><div class="gmail_quote">


<div><div>
On 13 November 2012 13:54, Jo <span dir="ltr"><<a href="mailto:winfixit@gmail.com" target="_blank">winfixit@gmail.com</a>></span> wrote:<br></div></div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">


<div><div>
Hi,<br>
<br>
I've been reading docs, but I'm afraid I can't figure it out on my own.<br>
<br>
I have a csv-file:<br>
POI_ID;POI_name;coordinateX;coordinateY<br>
<br>
<br>
Where X and Y are in Lambert72.<br>
<br>
<br>
How can I get these in a table that includes 2 columns with those same<br>
coordinates in WGS84? Should I have a column with a geometry type and<br>
a point in it? Can I create a trigger function which updates the WGS84<br>
as the Lambert72 are inserted?<br>
</div></div><span><font color="#888888"><br>
Jo<br>
_______________________________________________<br>
postgis-users mailing list<br>
<a href="mailto:postgis-users@lists.osgeo.org" target="_blank">postgis-users@lists.osgeo.org</a><br>
<a href="http://lists.osgeo.org/cgi-bin/mailman/listinfo/postgis-users" target="_blank">http://lists.osgeo.org/cgi-bin/mailman/listinfo/postgis-users</a><br>
</font></span></blockquote></div><br></div>
</blockquote></div><br></div>
</div></div></blockquote></div><br></div>
</div></div></blockquote></div><br></div>