Wie finde ich den nächsten Punkt mithilfe der PostGIS-Funktion?

8

Tabelle A:

lat | long | the_geom | code | sign 
13.8433095 | 100.6360357 | 0101000020E61.... | ABC | start_point
13.7544738 | 100.5459646 | 0101000020E6..... | ABC | end_point
13.4124215 | 100.6232332 | 0101000020E61.... | DEF | start_point
13.2423438 | 100.2324426 | 0101000020E6..... | DEF | end_point

Tabelle B:

lat | long | the_geom | code    
13.7546285 | 100.5458729 | 0101000020E.... | ABC    
13.7546698 | 100.5458513 | 0101000020E.... | ABC
13.7547107 | 100.5458233 | 0101000020E.... | DEF    
...

Ich möchte die nächsten Punkte von jedem Punkt (Start- und Endpunkt von Tabelle A) finden, die mit jedem Punkt mit demselben 'Code' in Tabelle B verglichen werden.

Was ist die beste PostGIS-Funktion / PostgreSQL-Abfrage, um dies zu lösen? Was ich versucht habe sind ..

ST_Distance_Spheroid und ST_DWithin und ST_Distance

Aber ich habe keine Ahnung, wie ich sie in einer einzigen Abfrage kombinieren soll.

user9404
quelle

Antworten:

9

Vielleicht so etwas (ich gehe davon aus, dass Sie in jeder Tabelle eine Primärschlüsselspalte "id" haben):

SELECT A.id, A.code AS Code, A.sign AS Sign, B.id, 
   ST_Distance(A.geom, B.geom) AS Distance 
FROM Table_A AS A, Table_B AS B
WHERE A.id IN ( 
   SELECT X.id
   FROM TableA as X, TableB as Y
-- Here's the important part: refer to the A table **outside** of the subquery
   WHERE X.code=Y.code AND X.id=A.id
-- and use ORDER BY ... LIMIT 1 to get the closest point
   ORDER BY ST_Distance(X.geom, Y.geom) ASC LIMIT 1
);
Micha
quelle