PostGIS: Tetapkan ID titik di lapisan A ke titik terdekat di Lapisan B

Ini seharusnya merupakan prekursor yang jelas (yang tidak saya tanyakan) untuk pertanyaan saya yang lain: Bagaimana cara membuat diagram laba-laba (garis hub) di PostGIS?

Jika saya tidak tahu hubungan antara titik di lapisan A (toko) dan titik di lapisan B (pelanggan), saya biasanya ingin mengatakan "Pelanggan 1 dilayani oleh toko terdekat." Sementara saya menyadari fakta ini mungkin tidak benar, itu bisa menjadi pengganti yang layak.

Menggunakan PostGIS, apa cara paling efisien untuk menetapkan ID titik terdekat di lapisan A (toko) ke setiap titik di lapisan B (pelanggan). Output yang saya cari adalah sesuatu seperti di bawah ini.

Customer | Store
    1    |   A
    2    |   A
    3    |   B
    4    |   C

juga:

pilih A.ID sebagai CUST_ID, (pilih B.ID dari B order dengan st_distance (A.geom, B.geom) batas 1) sebagai STORE_ID dari A

Sepertinya jika Anda memiliki lebih banyak pelanggan daripada toko, maka mungkin akan lebih efisien untuk membuat lapisan poligon voronoi untuk toko, lalu melakukan penggabungan spasial pelanggan terhadap poligon toko.

Diskusi tentang beberapa solusi tetangga terdekat dasar dapat ditemukan di sini: http://www.bostongis.com/?content_name=postgis_nearest_neighbor#120

/ Nicklas

Dari http://www.bostongis.com/?content_name=postgis_nearest_neighbor :

Jika Anda perlu mendapatkan tetangga terdekat untuk semua catatan dalam sebuah tabel, tetapi Anda hanya perlu tetangga terdekat pertama untuk masing-masing, maka Anda dapat menggunakan sintaks berbeda DISTINCT ON PostgreSQL. Yang akan terlihat seperti ini:

SELECT DISTINCT ON(g1.gid)  g1.gid As gref_gid, 
       g1.description As gref_description, 
       g2.gid As gnn_gid, 
       g2.description As gnn_description  
FROM sometable As g1, sometable As g2   
WHERE g1.gid <> g2.gid 
      AND ST_DWithin(g1.the_geom, g2.the_geom, 300)   
ORDER BY g1.gid, ST_Distance(g1.the_geom,g2.the_geom) 

Ini akan menemukan jarak minimum hingga 300 unit. Jadi, Anda harus terlebih dahulu memeriksa data Anda dan mencari tahu seberapa besar jarak minimum Anda akan dapatkan.

Anda harus menemukan pasangan terpendek di kotak pencarian, dan jika kotak itu tidak ada di dalamnya, perluas. Itu tidak cantik tapi berhasil. Ada contoh kode PL / PgSQL di sini http://trac.osgeo.org/postgis/wiki/UsersWikiNearest

Terima kasih atas masukan semua orang. Saya akhirnya pergi dengan kombinasi saran eprand dan underdark. Kode terakhir yang saya gunakan adalah:

CREATE TABLE closest_point as
SELECT DISTINCT ON (A.GID) A.GID AS CUST_ID, 
      (SELECT B.GID FROM "STORES" as B 
       ORDER BY ST_Distance(A.the_geom, B.the_geom) limit 1) as STORE_ID, 
       A.the_geom 
FROM "CUSTOMERS" as A, "STORES" as B;

Saya kemudian membuat diagram voronoi pada lapisan toko untuk mengkonfirmasi hasil yang bekerja dengan benar, yang tentu saja mereka lakukan. Terima kasih atas kerja kerasnya semua!