Saya seorang pemula di postgis dan saya memiliki masalah dalam kinerja permintaan.

Ini pertanyaan saya:

SELECT DISTINCT ON (userid) userid ,ST_AsText(position), timestamp  
FROM table1 
WHERE ST_Intersects ( ST_GeomFromText('a multiypolygon geom goes here',4326),position) 
ORDER BY userid, timestamp desc

dan masalahnya adalah multipoligon saya termasuk poligon SANGAT besar (panjang 600 halaman dalam kata doc!) dan butuh lebih dari 2 jam untuk mengeksekusi!

Apakah ada cara untuk mengoptimalkan permintaan saya atau menggunakan cara lain?

Tolong bantuan Anda sangat dihargai!

Yang harus Anda lakukan adalah meletakkan multipolygon besar Anda dalam sebuah tabel sebagai poligon tunggal (dengan ST_Dump) dan meletakkan indeks di atasnya. Sesuatu seperti:

CREATE TABLE big_polygon as
SELECT (ST_Dump( ST_GeomFromText('a multiypolygon geom goes here',4326))).geom as geom;

-- It is always great to put a primary key on the table
ALTER table big_polygon ADD Column gid serial PRIMARY KEY;

-- Create the index
CREATE INDEX idx_big_polygon_geom
on big_polygon
USING gist(geom);

-- To give the database some information about how the index looks
analyze big_polygon;

-- Then you go:
SELECT DISTINCT ON (userid) userid ,ST_AsText(position), timestamp  
FROM table1, big polygon WHERE ST_Intersects ( big_polygon.geom,position) 
ORDER BY userid, timestamp desc;

Itu harus lebih cepat karena beberapa alasan.

Tergantung jenis kualitas - ketelitian yang Anda butuhkan. Anda jelas dapat menyederhanakan poligon dengan menggunakan: http://postgis.net/docs/ST_Simplify.html

Apa yang sering saya lakukan selama pengembangan aplikasi GIS saya adalah memikirkan cara terbaik untuk meminimalkan data. Misalnya. pilih dulu poligon di dalam kotak batas misalnya. - Bergantung pada tingkat zoom Anda tidak perlu hasil yang sangat presisi (st_simplify) begitu ...

Semoga itu sedikit membantu Anda!

Bergantung pada keahlian postgres dan / atau sql Anda, Anda memiliki beberapa opsi:

1. menganalisis kueri melalui perintah EXPLAIN untuk mencari tahu apakah Anda menekan hambatan tertentu. Peringatan: kadang-kadang output EXPLAIN mungkin sulit dipahami

2. jika Anda berharap bahwa sebagian besar atau sebagian besar geometri dalam tabel1 JANGAN memotong berpoligon, Anda bisa mencoba menerapkan kondisi pendahuluan terhadap poligon yang lebih sederhana (yaitu dengan memecah multiploygon dalam potongan yang lebih kecil) dan kemudian jalankan persimpangan multipoligon yang lebih berat hanya pada hasil tersebut. Lihat di bawah untuk contoh.

3. jika dan hanya jika CPU adalah hambatan (yaitu server macet di persimpangan komputer) Saya sepenuhnya menyarankan Anda mendapatkan CPU yang lebih besar, lebih cepat, lebih kuat atau menyewa Instance CPU Tinggi satu kali dari EC2 Amazon dan menghancurkannya ketika Anda selesai

Contoh permintaan untuk item 2:

SELECT DISTINCT ON (st1.userid) st1.userid ,ST_AsText(st1.position), st1.timestamp  
FROM (
    select userid, position, timestamp from table1 
    WHERE ST_Intersects ( YOUR_MULTIPOL_BOUNDS_HERE,position)
) as st1 
WHERE ST_Intersects ( ST_GeomFromText('a multiypolygon geom goes     here',4326),st1.position) 
ORDER BY st1.userid, st1.timestamp desc

Untuk meningkatkan kinerja, Anda juga dapat secara sementara mematerialisasi subselect st1 sebagai tabel sehingga Anda dapat mengindeksnya.

@Nicklas benar untuk menunjukkan dalam komentar bahwa contoh untuk saran 2 seharusnya tidak membantu. Dia benar, tapi saya pikir saya juga (sebagian) benar.

Bahkan sepertinya pertanyaan yang sangat mirip ditanyakan (dan dijawab) pada November lalu di postgis ML:

http://postgis.refractions.net/pipermail/postgis-users/2011-November/031344.html

dan ternyata sarannya adalah untuk benar-benar memecah poligon sehingga indeks paling efektif dapat memfilter persimpangan palsu yang akan dipicu oleh pemeriksaan batas sederhana.

Menggunakan ST_SubDivide()

Untuk Postgis versi 2.2, Anda dapat menggunakan ST_SubDivide.

ST_Subdivide - Mengembalikan set geometri di mana tidak ada geometri dalam set memiliki lebih dari jumlah simpul yang ditentukan.

setof geometry ST_Subdivide(geometry geom, integer max_vertices=256);

Anda juga bisa

• gunakan tabel temp
• sebuah indeks

Di sini kita gunakan ST_SubDivideuntuk memecah poligon menjadi subpoligon dengan 10 atau lebih sedikit simpul.

CREATE TEMP TABLE divided AS
SELECT ST_SubDivide(bigmultipolygon,10)::geometery AS t(geom);

CREATE INDEX divided_idx ON divided USING gist(geom);

Kemudian

SELECT DISTINCT ON (userid) userid ,ST_AsText(position), timestamp  
FROM table1
JOIN divided AS d
  ON ST_Intersects( d.geom, position )
ORDER BY userid, timestamp desc;

Jangan lakukan hal di atas, itu menyebabkan kesalahan pembulatan

Tuning Umum

Lihat juga bagian yang berjudul Tip Kinerja dalam dokumen. Pastikan Anda disetel dengan tepat. Pertimbangkan meningkatkan max_parallel_workers_per_gatheruntuk memanfaatkan paralelisasi (saat ini default tidak aktif).