1 V. INSTALASI POSTGRESQL DAN SEKILAS TENTANG POSTGIS Bagian ini akan mengulas secara sederhana tentang POSTGIS (ekstensi spasial untuk database Postg...
V. INSTALASI POSTGRESQL DAN SEKILAS TENTANG POSTGIS
Bagian ini akan mengulas secara sederhana tentang POSTGIS (ekstensi spasial untuk database PostgreSQL). Database yang digunakan : postgresql-8.2.3-1 ; pada platform Windows
INSTALASI POSTGRESQL - POSTGIS Ada berbagai cara untuk membangun database PostgreSQL yang mendukung data spasial (postgis). Akan tetapi cara mudahnya adalah menggunakan postgresql-8.2.msi (installer windows dengan fitur PostgreSQL lengkap, termasuk PostGIS) yang dapat di download di situsnya : http://www.postgresql.org. Kemudian, jangan lupa kunjungi situs PostGIS untuk informasi lebih lengkap, termasuk download dokumentasi/manual posGIS (http://www.postgis.org)
OBJEK SPASIAL PostGIS mendukung semua objek spasial yang di spesifikasikan oleh OpenGIS Consortium (http://www.opengis.org) pada dokumen Simple Features for Specification for SQL (http://www.opengis.org/techno/specs/99-049.pdf). Postgis juga mengembangkan kapabilitasnya dengan 3DZ, 3DM dan 4D koordinat.
WKB and WKT OpenGIS menspesifikasikan ada2 cara standard untuk mengexpresikan objek spasial : format Well-Known Text (WKT) dan Well-Known Binary (WKB).
WebgGIS dengan MapServer – MapLab – Chameleon - PostGIS
Input data pada database biasanya menggunakan format WKT seperti diatas. Akan tetapi data spasial disimpan pada database dengan format heksadesimal. Kita dapat memanipulasi format Input/Output, sebagai berikut:
bytea WKB = asBinary(geometry); text WKT = asText(geometry); geometry = GeomFromWKB(bytea WKB, SRID); geometry = GeometryFromText(text WKT, SRID); NB : SRID = spatial referencing system identifier.
contoh query insert :
INSERT INTO SPATIALTABLE ( THE_GEOM, THE_NAME ) VALUES ( GeomFromText(’POINT(-126.4 45.32)’, 312), ’A Place’ ) ;
MEMBUAT DATABASE SPASIAL Cara paling mudah, adalah menggunakan pgAdmin, yang shortcut nya sudah tersedia pada start menu Windows, pada folder PostgreSQL. Setelah pgAdmin dijalankan, pilih menu Edit Æ New Object Æ New Database. Masukkan pilihan database template : template_postgis
TABEL SPATIAL_REF_SYS DAN TABEL GEOMETRY_COLUMNS Kita akan melihat pada database kita ada beberapa table yang telah terbuat. Salah satunya SPATIAL_REF_SYS dan GEOMETRY_COLUMNS. Yang digunakan untuk menyimpan informasi meta dari objek spasial yang kita buat.
select * from spatial_ref_sys ; select * from geometry_columns ;
MEMBUAT TABEL SPASIAL Langkah pertama adalah membuat database seperti biasa. Misal:
CREATE TABLE ROADS_GEOM ( ID int4, NAME varchar(25) )
www.hatma.info
1
WebgGIS dengan MapServer – MapLab – Chameleon - PostGIS
Langkah kedua adalah menambahkan kolom spasial dengan fungsi AddGeometryColumn. Dengan format : “SELECT AddGeometryColumn(, , <srid>, , )”. Misal :
Perhatikan bahwa OGC geometries tidak memperbolehkan nilai Z atau M. Sehingga bagian dimension cukup kita beri nilai 2.
MENJAGA KONSISTENSI DATA PADA KOLOM SPASIAL Pastikan kolom geometri kita memenuhi spesifikasi standard OpenGIS. Salah satunya dengan fungsi IsValid(). Misal :
Ada cara untuk mempermudah kita menjaga data supaya tetap memenuhi standar, yaitu dengan membuat satu check constraint . Misal : ALTER TABLE mytable ADD CONSTRAINT geometry_valid_check CHECK (isvalid(the_geom));
check constraint juga akan dibuat secara otomatis apabila kita menambahkan kolom geometris dengan perintah AddGeometryColumn().
MEMASUKKAN DATA SPATIAL Ada beberapa metode yang bisa kita gunakan. Cara pertama adalah dengan menggunakan perintah query insert seperti biasa. Misal kita membuat file (roads.sql) yang berisi :
www.hatma.info
2
WebgGIS dengan MapServer – MapLab – Chameleon - PostGIS
BEGIN; INSERT INTO ROADS_GEOM (ID,GEOM,NAME ) VALUES (1,GeomFromText(’LINESTRING(191232 243118,191108 243242)’,-1),’Jeff Rd’); INSERT INTO ROADS_GEOM (ID,GEOM,NAME ) VALUES (2,GeomFromText(’LINESTRING(189141 244158,189265 244817)’,-1),’Geordie Rd’); INSERT INTO ROADS_GEOM (ID,GEOM,NAME ) VALUES (3,GeomFromText(’LINESTRING(192783 228138,192612 229814)’,-1),’Paul St’); INSERT INTO ROADS_GEOM (ID,GEOM,NAME ) VALUES (4,GeomFromText(’LINESTRING(189412 252431,189631 259122)’,-1),’Graeme Ave’); INSERT INTO ROADS_GEOM (ID,GEOM,NAME ) VALUES(5,GeomFromText(’LINESTRING(190131 224148,190871 228134)’,-1),’Phil Tce’); INSERT INTO ROADS_GEOM (ID,GEOM,NAME ) VALUES(6,GeomFromText(’LINESTRING(198231 263418,198213 268322)’,-1),’Dave Cres’); COMMIT;
Kemudian eksekusi dengan perintah pada "psql" SQL terminal monitor:
psql -d [database] -f roads.sql
MENGGUNAKAN SHP2PGSQL Jika kita telah mempunyai file shp yang berisi data-data spasial kita, lebih mudah menggunakan shp2pgsql data loader yang akan mengkonversi file ESRI shp ke dalam SQL. Shp2pgsql dapat kita jalankan pada "psql" SQL terminal monitor. Misal :
Atau : shp2pgsql shaperoads myschema.roadstable | psql -d roadsdb
MENAMPILKAN DATA SPASIAL Untuk menampilkan data spatial, kita gunakan perintah SELECT. Misal :
SELECT id, AsText(geom) AS geom, name FROM ROADS_GEOM;
Kita bisa menggunakan beberapa operator pada query kita :
www.hatma.info
3
WebgGIS dengan MapServer – MapLab – Chameleon - PostGIS
o && Mengatakan apakah batas luar geometry sebuah objek intersects / beririsan dengan batas luar geometry objek lain. o ~= Mengetes apakah 2 buah geometry adalah identik. Misal, apakah ’POLYGON((0 0,1 1,1 0,0 0))’ identik dengan ’POLYGON((0 0,1 1,1 0,0 0))’ Æ(ya). o = Mengetes apakah batas luar 2 buah geometry adalah identik
Contoh : SELECT ID, NAME FROM ROADS_GEOM WHERE GEOM ~= GeomFromText(’LINESTRING(191232 243118,191108 243242)’,-1);
SELECT ID, NAME FROM ROADS_GEOM WHERE GEOM && GeomFromText(’POLYGON((191232 243117,191232 243119,191234 243117,191232 243117))’,-1);
MENGGUNAKAN PGSQL2SHP Kita dapat menggunakan pgsql2shp untuk mengkonversi tabel (atau data spasial tertentu pada database kita) ke dalam format file ESRI .shp
pgsql2shp [] [<schema>.]
pgsql2shp []
MEMBUAT INDEX SPATIAL Index digunakan untuk mempercepat pencarian data. PostGIS Menggunakan GiST (Generalized Search Trees) indexes.
CREATE INDEX [indexname] ON [tablename] USING GIST ( [geometryfield] GIST_GEOMETRY_OPS );
Tabel dengan 1 juta barisdata, pada mesin 300MHz Solaris, membutuhkan waktgu komputasi 1 jam.
www.hatma.info
4
WebgGIS dengan MapServer – MapLab – Chameleon - PostGIS
Kemudian eksekusi perintah untuk membuat / mengupdate geometry statistic yang sangat penting, untuk optimasi eksekusi query oleh database.
VACUUM ANALYZE [table_name] ( [column_name] );
Sebaiknya Vacuum Database dilakukan secara rutin, seperti pada kebanyakan PostgreSQL DBA melakukan VACUUM dengan cron job secara periodic, untuk menjaga agar index dapat berfungsi secara baik dalam mempercepat query.
COMPLEX QUERIES Kita dapat menggunakan fungsi-fungsi geometry dalam perintah query kita. Misalkan fungsi distance() yang digunakan untuk mengukur jarak antar 2 buah objek geometry.
SELECT the_geom FROM geom_table WHERE distance( the_geom, GeomFromText( ’POINT(100000 200000)’, -1 ) ) < 100
Kita dapat menggunakan fungsi expand() untuk menemukan objek geometri dalam radius dari suatu objek:
SELECT * FROM GEOTABLE WHERE GEOCOLUMN && Expand(GeomFromText(’POINT(1000 1000)’,-1),100) AND Distance(GeomFromText(’POINT(1000 1000)’,-1),GEOCOLUMN) < 100;
Fungsi spasial selengkapnya dapat dilihat pada manual/postGIS reference.
CONTOH IMPLEMENTASI SQL SPATIAL
Buat table dan isinya dengan perintah sebagai berikut.
1. Mencari total panjang jalan (dalam kilometer)? SELECT sum(length(the_geom))/1000 AS PANJANG_JALAN FROM bc_roads ;
2. Letak perpotongan semua jalan dengan titik (1 1) SELECT ASTEXT( INTERSECTION(the_geom, GeomFromText('POINT( 1 1 )')) ) FROM bc_roads;
3. Luas Kota bernama ‘AREA PERTAMA’ (dalam hektar) SELECT area(the_geom)/10000 AS LUAS_HEKTAR FROM bc_municipality WHERE name='AREA PERTAMA' ;
4. Apakah kota terbesar ? SELECT name, area(the_geom)/10000 AS LUAS_HEKTAR FROM bc_municipality LIMIT 1 ;
www.hatma.info
6
WebgGIS dengan MapServer – MapLab – Chameleon - PostGIS
Atau SELECT max(area(the_geom)/10000) AS LUAS_HEKTAR FROM bc_municipality ;
5. Berapa panjang jalan pada masing-masing kota? SELECT m.name, SUM(LENGTH(r.the_geom))/1000 AS KM_JALAN FROM bc_roads r, bc_municipality m WHERE CONTAINS(m.the_geom, r.the_geom) AND m.the_geom && r.the_geom GROUP BY m.name ;
6. Buat table baru dengan isi : semua jalan pada kota ’AREA PERTAMA’ dengan tambahan kolom panjang jalan tersebut disemua kota , dan titik perpotongannya dengan kota AREA PERTAMA. CREATE TABLE pg_roads AS SELECT
r.*,
ASTEXT( INTERSECTION(r.the_geom, m.the_geom)) AS "PERPOTONGAN DI :", LENGTH(r.the_geom)/1000 AS "PANJANG ASLI JALAN (KM)" FROM
bc_roads r, bc_municipality m
WHERE
r.the_geom && m.the_geom
AND
CONTAINS(m.the_geom, r.the_geom)
AND m.name = 'AREA PERTAMA' ;
7. Berapa kilometer panjang jalan “JALAN PERTAMA” di kota “AREA PERTAMA” ? SELECT SUM(LENGTH(r.the_geom))/1000 AS "KM PANJANGNYA =" FROM bc_roads r, bc_municipality m WHERE m.name = 'AREA PERTAMA' AND r.name = 'JALAN PERTAMA' AND m.the_geom && r.the_geom
8. Cari kota terbesar yang memiliki lubang SELECT gid, area(the_geom)/1000 AS "LUAS KOTA" FROM bc_municipality WHERE nrings(the_geom) > 1 ORDER BY "LUAS KOTA" DESC LIMIT 1