Menggambar garis antara titik-titik pada jarak tertentu di PostGIS?

Saya memiliki data titik di sepanjang jalan, saya ingin mengubah titik-titik itu menjadi garis berwarna sederhana. Adakah petunjuk apa masalah ini dapat disebut atau algoritma yang dapat membantu saya menyelesaikan ini?

Saya berharap menggunakan PostGISfungsi untuk melakukan ini tetapi saya terbuka untuk saran, ini adalah data dari .shpfile.

Sunting1: Diperbarui gambar untuk menunjukkan solusi ideal dari masalah ini.

Menggambar garis akan murni berdasarkan jarak antara titik-titik itu, tidak ada lagi yang bisa saya gunakan untuk mengelompokkannya. Idealnya ini adalah titik pada jarak maksimum yang ditentukan di sepanjang garis yang diproyeksikan? Dan dengan memproyeksikan garis yang saya maksud menemukan titik 1 kemudian berikutnya yang terdekat dengan itu memproyeksikan garis dan memeriksa apakah ada titik pada garis ini pada jarak maksimal ke salah satu yang sudah ada di garis.

​

Anda dapat menggunakan kueri rekursif untuk menjelajahi tetangga terdekat dari setiap titik mulai dari setiap ujung garis yang terdeteksi yang ingin Anda buat.

Prasyarat : siapkan layer postgis dengan poin Anda dan poin lain dengan objek multi-linestring tunggal yang berisi jalan Anda. Dua lapisan harus berada pada CRS yang sama. Ini adalah kode untuk set data uji yang saya buat, harap modifikasi sesuai kebutuhan. (Diuji pada postgres 9.2 dan postgis 2.1)

WITH RECURSIVE
points as (SELECT id, st_transform((st_dump(wkb_geometry)).geom,2154) as geom, my_comment as com FROM mypoints),
roads as (SELECT st_transform(ST_union(wkb_geometry),2154) as geom from highway),

Berikut langkah-langkahnya :

1. Buat untuk setiap titik daftar setiap tetangga dan jarak mereka yang memenuhi tiga kriteria ini.

   • Jarak tidak boleh melebihi ambang batas yang ditentukan pengguna (ini akan menghindari menghubungkan ke titik yang terisolasi)

     graph_full as (
     SELECT a.id, b.id as link_id, a.com, st_makeline(a.geom,b.geom) as geom, st_distance(a.geom,b.geom) as distance
     FROM points a
     LEFT JOIN points b ON a.id<>b.id
     WHERE st_distance(a.geom,b.geom) <= 15
     ),

   • Jalur langsung tidak boleh melewati jalan

     graph as (
     SELECt graph_full.*
     FROM graph_full RIGHT JOIN
     roads ON st_intersects(graph_full.geom,roads.geom) = false
     ),

   • Jarak tidak boleh melebihi rasio jarak yang ditentukan pengguna dari tetangga terdekat (ini harus mengakomodasi lebih baik digitalisasi tidak teratur daripada jarak tetap) Bagian ini sebenarnya terlalu sulit untuk diterapkan, melekat pada radius pencarian tetap

   Sebut saja tabel ini "grafik"

2. Pilih ujung titik garis dengan bergabung ke grafik dan pertahankan hanya titik yang memiliki satu entri persis dalam grafik.

   eol as (
   SELECT points.* FROM
   points  JOIN
   (SELECT id, count(*) FROM graph 
   GROUP BY id
   HAVING count(*)= 1) sel
   ON points.id = sel.id),

   Sebut saja tabel ini "eol" (end of line) dengan
   mudah? bahwa hadiah untuk melakukan grafik yang bagus tetapi hal-hal yang bertahan akan menjadi gila pada langkah berikutnya

3. Menyiapkan kueri rekursif yang akan berputar dari tetangga ke tetangga mulai dari setiap hari

   • Inisialisasi kueri rekursif menggunakan tabel eol dan menambahkan penghitung untuk kedalaman, agregator untuk path, dan konstruktor geometri untuk membangun garis
   • Pindah ke iterasi berikutnya dengan beralih ke tetangga terdekat menggunakan grafik dan periksa bahwa Anda tidak pernah mundur menggunakan jalur
   • Setelah iterasi selesai, simpan hanya jalur terpanjang untuk setiap titik awal (jika dataset Anda menyertakan persimpangan potensial antara garis ekspektasi, bagian itu akan membutuhkan lebih banyak kondisi)

   recurse_eol (id, link_id, depth, path, start_id, geom) AS (--initialisation
   SELECT id, link_id, depth, path, start_id, geom FROM (
       SELECT eol.id, graph.link_id,1 as depth,
       ARRAY[eol.id, graph.link_id] as path,
       eol.id as start_id,
       graph.geom as geom,
       (row_number() OVER (PARTITION BY eol.id ORDER BY distance asc))=1 as test
       FROM eol JOIn graph ON eol.id = graph.id 
       ) foo
   WHERE test = true
   
   UNION ALL ---here start the recursive part
   
   SELECT id, link_id, depth, path, start_id, geom  FROM (
       SELECT graph.id, graph.link_id, r.depth+1 as depth,
       path || graph.link_id as path,
       r.start_id,
       ST_union(r.geom,graph.geom) as geom,
       (row_number() OVER (PARTITION BY r.id ORDER BY distance asc))=1 as test
       FROM recurse_eol r JOIN graph ON r.link_id = graph.id AND NOT graph.link_id = ANY(path)) foo
   WHERE test = true AND depth < 1000), --this last line is a safe guard to stop recurring after 1000 run adapt it as needed

   Sebut saja tabel ini "recurse_eol"

4. Pertahankan hanya garis terpanjang untuk setiap titik awal dan hapus setiap jalur duplikat yang tepat Contoh: jalur 1,2,3,5 DAN 5,3,2,1 adalah garis yang sama yang ditemukan oleh dua perbedaan "ujung jalur"

   result as (SELECT start_id, path, depth, geom FROM
   (SELECT *,
   row_number() OVER (PARTITION BY array(SELECT * FROM unnest(path) ORDER BY 1))=1 as test_duplicate,
   (max(depth) OVER (PARTITION BY start_id))=depth as test_depth
   FROM recurse_eol) foo
   WHERE  test_depth = true AND test_duplicate = true)
   
   SELECT * FROM result

5. Secara manual memeriksa kesalahan yang tersisa (titik terisolasi, garis yang tumpang tindih, jalan berbentuk aneh)

Diperbarui seperti yang dijanjikan, saya masih tidak tahu mengapa kadang-kadang permintaan rekursif tidak memberikan hasil yang sama persis ketika mulai dari lawan dari baris yang sama sehingga beberapa duplikat mungkin tetap ada di lapisan hasil seperti sekarang.

Jangan ragu untuk bertanya, saya benar-benar mengerti bahwa kode ini membutuhkan lebih banyak komentar. Ini pertanyaan lengkapnya:

WITH RECURSIVE
points as (SELECT id, st_transform((st_dump(wkb_geometry)).geom,2154) as geom, my_comment as com FROM mypoints),
roads as (SELECT st_transform(ST_union(wkb_geometry),2154) as geom from highway),

graph_full as (
    SELECT a.id, b.id as link_id, a.com, st_makeline(a.geom,b.geom) as geom, st_distance(a.geom,b.geom) as distance
    FROM points a
    LEFT JOIN points b ON a.id<>b.id
    WHERE st_distance(a.geom,b.geom) <= 15
    ),

graph as (
    SELECt graph_full.*
    FROM graph_full RIGHT JOIN
    roads ON st_intersects(graph_full.geom,roads.geom) = false
    ),

eol as (
    SELECT points.* FROM
    points  JOIN
        (SELECT id, count(*) FROM graph 
        GROUP BY id
        HAVING count(*)= 1) sel
    ON points.id = sel.id),


recurse_eol (id, link_id, depth, path, start_id, geom) AS (
    SELECT id, link_id, depth, path, start_id, geom FROM (
        SELECT eol.id, graph.link_id,1 as depth,
        ARRAY[eol.id, graph.link_id] as path,
        eol.id as start_id,
        graph.geom as geom,
        (row_number() OVER (PARTITION BY eol.id ORDER BY distance asc))=1 as test
        FROM eol JOIn graph ON eol.id = graph.id 
        ) foo
    WHERE test = true

UNION ALL
    SELECT id, link_id, depth, path, start_id, geom  FROM (
        SELECT graph.id, graph.link_id, r.depth+1 as depth,
        path || graph.link_id as path,
        r.start_id,
        ST_union(r.geom,graph.geom) as geom,
        (row_number() OVER (PARTITION BY r.id ORDER BY distance asc))=1 as test
        FROM recurse_eol r JOIN graph ON r.link_id = graph.id AND NOT graph.link_id = ANY(path)) foo
    WHERE test = true AND depth < 1000),

result as (SELECT start_id, path, depth, geom FROM
    (SELECT *,
    row_number() OVER (PARTITION BY array(SELECT * FROM unnest(path) ORDER BY 1))=1 as test_duplicate,
    (max(depth) OVER (PARTITION BY start_id))=depth as test_depth
    FROM recurse_eol) foo
WHERE  test_depth = true AND test_duplicate = true)

SELECT * FROM result

Seperti yang ditunjukkan oleh @FelixIP, langkah pertama adalah menemukan poin yang akan membentuk setiap baris. Anda dapat melakukan ini dengan menghubungi ST_ClusterWithin dengan jarak pemisahan maksimum Anda:

SELECT
  row_number() OVER () AS cid, 
  (ST_Dump(geom)).geom 
FROM (
  SELECT unnest(st_clusterwithin(geom, 0.05)) AS geom 
  FROM inputs) sq

Kemudian, Anda harus menggunakan beberapa heuristik untuk membangun garis melalui semua titik di setiap cluster. Misalnya, jika Anda dapat menganggap baris yang diinginkan sebagai Y-monoton, Anda dapat mengurutkan titik di setiap kluster dan mengumpankannya ke dalam ST_MakeLine . Menggabungkan bahwa semuanya akan terlihat seperti ini:

SELECT 
  ST_MakeLine(geom ORDER BY ST_Y(geom)) AS geom
FROM (
  SELECT row_number() OVER () AS cid, 
  (ST_Dump(geom)).geom FROM (
    SELECT unnest(st_clusterwithin(geom, 0.05)) AS geom 
    FROM inputs) sq) ssq 
GROUP BY cid