In my first post I asked where am i? The gps in my phone said I was standing at 43.73066°N, 79.26482°W. In real life, I was standing at the junction of Kenmark Blvd and Chevron Cres.
With the Open Toronto centreline data, I can check the location of road intersections:
select distinct ( astext ( transform ( intersection ( r1.geometry, r2.geometry ), 4326 ) ) ) from centreline as r1, centreline as r2 where r1.lf_name = 'KENMARK BLVD' and r2.lf_name = 'CHEVRON CRES'
NULL POINT(-79.262423 43.730019) POINT(-79.264815 43.730591)
Hmm; three answers. NULL I can’t answer; I’ll put it down as an exhortation to Be Here Now. The second is given a clue by one of the street names: Chevron Crescent – first thing I learned when I had my newspaper round is that a crescent’s going to end you up back on the same road you started from. The last one, though, agrees to 5(ish) decimal places – well within the accuracy of my simple GPS.
Update: This query gets really, really slow on long streets. Both Chevron and Kenmark only have two road segments. Kennedy and Steeles East each have over 90.