PostGIS - National Education Center for GIS: Open Source GIS

2012년 8월 20일 ~ 21일
국토해양부 공간정보 거점대학
오픈 소스 GIS 과정

OSGeo 한국어 지부
망고시스템㈜
이민파 (mapplus@gmail.com)

목차

1. PostGIS 개요
2. PostGIS 설치 및 환경 설정
3. 공간 데이터베이스의 생성
4. GIS 자료의 Import
5. 좌표계 설정 및 변경
6. PostGIS Objects & Function
7. PostGIS 함수 활용
8. Backup & Restore
9. pgRouting 소개
2

 1.1 What is a Spatial Database?
 1.2 PostGIS
 1.3 PostGIS를 지원하는 GIS 프로그램들

3

 Spatial databases store/manipulate
spatial objects
 – data types, indexes, and functions
 Spatial data types
 shapes - point, line, polygon
 Spatial indexing
 efficient processing of spatial operations
 Spatial functions,
 querying of spatial properties and relationships.

4

http://www.postgresql.org/
http://www.postgis.org/
5

Open/Free Closed/Proprietary
•Loading/Extracting •Loading/Extracting
• Shp2Pgsql • Safe FME Desktop Translator/Converter
• ogr2ogr •Web-Based
• Dxf2PostGIS • Ionic Red Spider (now ERDAS)
•Web-Based • Cadcorp GeognoSIS
• Mapserver • Iwan Mapserver
• GeoServer (Java-based WFS / WMS -se • MapDotNet Server
rver ) • MapGuide Enterprise (using FDO)
• SharpMap SDK - for ASP.NET 2.0 • ESRI ArcGIS Server 9.3+
• MapGuide Open Source (using FDO) •Desktop
•Desktop • Cadcorp SIS
• uDig • Microimages TNTmips GIS
• QGIS • ESRI ArcGIS 9.3+
• mezoGIS • Manifold
• OpenJUMP • GeoConcept
• OpenEV • MapInfo (v10)
• SharpMap SDK for Microsoft.NET 2.0 • AutoCAD Map 3D (using FDO)
• ZigGIS for ArcGIS/ArcObjects.NET
• GvSIG
• GRASS

6

HTTP TMS, WMTS WMS, WFS, WFS-T

JDBC

JDBC
ODBC

7

 2.1 PostGIS 설치
 2.2 Dashboard
 2.3 Starting and Stopping
 2.4 환경 설정
 2.5 PostGIS Web Administration

8

 1. PostgreSQL 단독 설치
 Application Stack Builder
 Binaries 수동 설치
 http://postgis.refractions.net/download
 pgRouting을 사용하는 경우 8.4.x 버전 설치할 것
 http://www.pgrouting.org/download.html

 2. OpenGeoSuite 통합 설치
 PostGIS와 함께 설치
 OpenGeoSuite 2.5 = PostgreSQL 8.4.9 + PostGIS
1.5.4 + GeoServer + GeoWebCache + GeoExplorer +
Client SDK
 http://opengeo.org/technology/suite/download/

9

PostgreSQL + Application Stack Builder

PostgreSQL: http://www.postgresql.org/download/
10

PostgreSQL + PostGIS Binaries

 아래 예는 postgreSQL 9.0 버전에 PostGIS 1.5.4.2 버전을
설치하는 과정
 http://postgis.refractions.net/download/windows/ 이동
 설치된 PostgreSQL 버전과 호환하는 PostGIS 1.5.x 바이
너리 버전 다운로드, ex) postgis-pg90-binaries-1.5.4-
2.zip
 압축 해제 후 postgis-pg90-binaries-1.5.4-2 폴더로 이동
 makepostgisdb.bat 파일의 연결정보 및 경로를 수정
 makepostgisdb.bat 파일 실행하여 PostGIS 설치

11

makepostgisdb.bat 파일
 set PGPORT=5432
 set PGHOST=localhost
 set PGUSER=postgres
 set PGPASSWORD=postgis
 set THEDB=template_postgis
 set PGBIN=C:Program FilesPostgreSQL9.0bin
 set PGLIB=C:Program FilesPostgreSQL9.0lib
 set POSTGISVER=1.5
 xcopy bin*.* "%PGBIN%"
 xcopy /I /S binpostgisgui* "%PGBIN%postgisgui"
 xcopy lib*.* "%PGLIB%"
 "%PGBIN%psql" -c "CREATE DATABASE %THEDB%"
 "%PGBIN%psql" -d "%THEDB%" -c "CREATE LANGUAGE plpgsql"
 "%PGBIN%psql" -d "%THEDB%" -f "sharecontribpostgis-%POSTGISVER%postgis.sql"
 "%PGBIN%psql" -d "%THEDB%" -f "sharecontribpostgis-%POSTGISVER%spatial_ref_sys.sql"
 "%PGBIN%psql" -d "%THEDB%" -f "sharecontribpostgis-%POSTGISVER%postgis_comments.sql"

 REM Uncomment the below line if this is a template database
 REM "%PGBIN%psql" -d "%THEDB%" -c "UPDATE pg_database SET datistemplate = true WHERE datname =
'%THEDB%';GRANT ALL ON geometry_columns TO PUBLIC; GRANT ALL ON spatial_ref_sys TO PUBLIC―
 pause

12

OpenGeoSuite 2.5 = PostgreSQL 8.4.9 + PostGIS 1.5.4 + GeoServer 2.1.x

다운로드: http://opengeo.org/technology/suite/download/
13

컴퓨터의 등록정보를 확인하여 사용자 이름이 한글일
경우 오류가 발생할 수 있으므로 변경 후 설치

14

도움말을 적극 활용하자
16

GeoServer의 Username: admin, Password:
geoserver 를 반드시 기억하십시오.

17

설치 디렉터리(C:Program FilesOpenGeoOpenGeo Suite)로 이동 후 다음의 명령을 실행.
시작하기
> opengeo-suite start
종료하기
> opengeo-suite stop
20

메뉴
툴바

객체 브라우저 View

25

도구 설명
pgAdmin 도구에 새로운 서버를 추가합니다.
객체의 정보를 생성, 수정, 갱신, 삭제 후 새로고침 합니다.
테이블 등의 객체의 등록정보를 확인하거나 수정합니다.
선택된 객체와 같은 형태의 새로운 객체를 생성합니다.
선택된 객체를 삭제합니다.
SQL을 실행할 수 있는 Query 도구를 불러옵니다.
선택된 테이블을 편집할 수 있는 테이블 편집기를 불러옵니다.
필터를 적용해서 테이블을 편집할 수 있는 테이블 편집기를 불러옵니다.

Vacuum, Analyze 등 테이블 유지에 필요한 도구를 불러옵니다.
26

 3.1 공간 데이터베이스의 생성
 3.2 공간 데이터베이스의 탐색

27

이름: seoul
오너: postgres
인코딩: UTF8
Template: template_postgis

CREATE DATABASE seoul
WITH ENCODING='UTF8'
TEMPLATE=template_postgis
CONNECTION LIMIT=-1;
28

geometry_columns 테이블은 PostGIS 내에 등록된
공간데이터의 메타 정보를 포함

29

SELECT postgis_full_version();

SELECT version();

SELECT 필드 FROM 테이블 WHERE 조건 30

 4.1 샘플 데이터셋의 준비
 4.2 GIS 자료의 Import
 4.3 레이어 확인

34

Shapefile 레이어명 유형 SRID 비고
admin_emd 읍면동경계 MultiPolygon 2097
admin_sgg 시군구경계 MultiPolygon 2097
admin_sid 시도경계 MultiPolygon 2097
river 실폭하천 MultiPolygon 2097
road_link2 도로 MultiLinestring 2097
firestation 소방서 Point 2097
healthcenter 보건소 Point 2097
policestation 경찰서 Point 2097
stores 대형마트 Point 2097
wardoffice 시군구청 Point 2097
subway 지하철노선 MultiLinestring 2097
Subway_station 지하철역 Point 2097
35

 PostGIS
 Using the Loader
 psql -d [database] -f roads.sql
 Using the Loader
 shp2pgsql -c -D -s 4269 -i -I shaperoads.shp public.roadstable >
roads.sql # psql -d roadsdb -f roads.sql
 pgsql2shp
 PostGIS Shapefile and DBF loader
 GUI – 실습에 사용

 Third Party
 Quantum GIS, uDIG…
 GDAL/OGR
 ogr2ogr -f "ESRI Shapefile" admin.shp PG:"host=localhost
user=postgres dbname=seoul password=postgis" ―admin" 39

OpenGeo Suite의 PostGIS를 사용한다면 54321, 독립형으로
PostgreSQL을 설치했을 경우 5432가 기본 포트입니다.
40

레이어에 대한 이름, geometry 컬럼명,
dimension, srid, geometry 유형 등의 정보
를 확인

43

 5.1 Unknown 좌표계 변환
 5.2 좌표계 설정
 5.3 좌표계 확인

44

PROJCS["Korean 1985 / Korea Central Belt",GEOGCS["Korean
1985",DATUM["Korean_Datum_1985",SPHEROID["Bessel
1841",6377397.155,299.1528128,AUTHORITY["EPSG","7004"]],AUTHORITY["EPSG","6162"]],PRIMEM["Greenwich",
0,AUTHORITY["EPSG","8901"]],UNIT["degree",0.01745329251994328,AUTHORITY["EPSG","9122"]],AUTHORITY["E
PSG","4162"]],UNIT["metre",1,AUTHORITY["EPSG","9001"]],PROJECTION["Transverse_Mercator"],PARAMETER["lati
tude_of_origin",38],PARAMETER["central_meridian",127],PARAMETER["scale_factor",1],PARAMETER["false_easting"
,200000],PARAMETER["false_northing",500000],AUTHORITY["EPSG","2097"],AXIS["Y",EAST],AXIS["X",NORTH]]

46

•첫째, srid 제약(Constraint)을 삭제.
•둘째, geometry_columns 테이블의 정의하고자 하는 레이어의
srid를 업데이트.
•셋째, 정의하고자 하는 레이어의 geometry column(일반적으로
the_geom)을 업데이트.
•넷째, srid 제약을 새로운 srid로 생성.

47

ALTER TABLE wardoffice DROP CONSTRAINT enforce_srid_the_geom;
UPDATE geometry_columns SET srid = 2097 WHERE f_table_name = 'wardoffice';
UPDATE wardoffice SET the_geom = ST_SetSRID(the_geom, 2097);
ALTER TABLE wardoffice ADD CONSTRAINT enforce_srid_the_geom CHECK (ST_Srid(the_geom) =
(2097));

48

 http://postgis.org/docs/UpdateGeometrySRID.html

 SELECT
UpdateGeometrySRID('public', 'admin_sgg', 'the_geo
m', 2097);

 Name
 UpdateGeometrySRID — Updates the SRID of all features in a geometry
column, geometry_columns metadata and srid table constraint
 Synopsis
 text UpdateGeometrySRID(varchar table_name, varchar column_name, integer srid);
 text UpdateGeometrySRID(varchar schema_name, varchar table_name, varchar column_
name, integer srid);
 text UpdateGeometrySRID(varchar catalog_name, varchar schema_name, varchar table_
name, varchar column_name, integer srid);
50

한국 주요 좌표계 EPSG코드 및 proj4 인자 정리
현재 PostGIS는 2096, 2097, 2098을 제외하고는 미등록상태
<4927> GRS80 Korea 2000 / Geocentric datum of Korea +proj=longlat +ellps=GRS80 +towgs84=0.0,0.0,0.0,0.0,0.0,0.0,0.0 +no_defs

<2096> Bessel Korean 1985 / East Belt +proj=tmerc +lat_0=38 +lon_0=129 +k=1 +x_0=200000 +y_0=500000 +ellps=bessel +units=m +no_defs
<2097> Bessel Korean 1985 / Central Belt +proj=tmerc +lat_0=38 +lon_0=127 +k=1 +x_0=200000 +y_0=500000 +ellps=bessel +units=m +no_defs
<2098> Bessel Korean 1985 / West Belt +proj=tmerc +lat_0=38 +lon_0=125 +k=1 +x_0=200000 +y_0=500000 +ellps=bessel +units=m +no_defs

<5167> Bessel Korean 1985 / East Sea Belt +proj=tmerc +lat_0=38 +lon_0=131 +k=1 +x_0=200000 +y_0=500000 +ellps=bessel +units=m +no_defs
<5168> Bessel Korean 1985 / Central Belt Jeju +proj=tmerc +lat_0=38 +lon_0=127 +k=1 +x_0=200000 +y_0=550000 +ellps=bessel +units=m +no_defs

<5173> Bessel Korean 1985 / Modified West Belt +proj=tmerc +lat_0=38 +lon_0=125.00289 +k=1 +x_0=200000 +y_0=500000 +ellps=bessel +units=m +no_defs
<5174> Bessel Korean 1985 / Modified Central Belt +proj=tmerc +lat_0=38 +lon_0=127.00289 +k=1 +x_0=200000 +y_0=500000 +ellps=bessel +units=m +no_defs
<5175> Bessel Korean 1985 / Modified Central Belt Jeju +proj=tmerc +lat_0=38 +lon_0=127.00289 +k=1 +x_0=200000 +y_0=550000 +ellps=bessel +units=m +no_defs
<5176> Bessel Korean 1985 / Modified East Belt +proj=tmerc +lat_0=38 +lon_0=129.00289 +k=1 +x_0=200000 +y_0=500000 +ellps=bessel +units=m +no_defs
<5177> Bessel Korean 1985 / Modified East Sea Belt +proj=tmerc +lat_0=38 +lon_0=131.00289 +k=1 +x_0=200000 +y_0=500000 +ellps=bessel +units=m +no_defs
<5178> Bessel Korean 1985 / Unified CS +proj=tmerc +lat_0=38 +lon_0=127.5 +k=0.9996 +x_0=1000000 +y_0=2000000 +ellps=bessel +units=m +no_defs

<5179> GRS80 Korea 2000 / Unified CS +proj=tmerc +lat_0=38 +lon_0=127.5 +k=0.9996 +x_0=1000000 +y_0=2000000 +ellps=GRS80 +units=m +no_defs
<5180> GRS80 Korea 2000 / West Belt +proj=tmerc +lat_0=38 +lon_0=125 +k=1 +x_0=200000 +y_0=500000 +ellps=GRS80 +units=m +no_defs
<5181> GRS80 Korea 2000 / Central Belt +proj=tmerc +lat_0=38 +lon_0=127 +k=1 +x_0=200000 +y_0=500000 +ellps=GRS80 +units=m +no_defs
<5182> GRS80 Korea 2000 / Central Belt Jeju +proj=tmerc +lat_0=38 +lon_0=127 +k=1 +x_0=200000 +y_0=550000 +ellps=GRS80 +units=m +no_defs
<5183> GRS80 Korea 2000 / East Belt +proj=tmerc +lat_0=38 +lon_0=129 +k=1 +x_0=200000 +y_0=500000 +ellps=GRS80 +units=m +no_defs
<5184> GRS80 Korea 2000 / East Sea Belt +proj=tmerc +lat_0=38 +lon_0=131 +k=1 +x_0=200000 +y_0=500000 +ellps=GRS80 +units=m +no_defs

<5185> GRS80 Korea 2000 / West Belt 2010 +proj=tmerc +lat_0=38 +lon_0=125 +k=1 +x_0=200000 +y_0=600000 +ellps=GRS80 +units=m +no_defs
<5186> GRS80 Korea 2000 / Central Belt 2010 +proj=tmerc +lat_0=38 +lon_0=127 +k=1 +x_0=200000 +y_0=600000 +ellps=GRS80 +units=m +no_defs
<5187> GRS80 Korea 2000 / East Belt 2010 +proj=tmerc +lat_0=38 +lon_0=129 +k=1 +x_0=200000 +y_0=600000 +ellps=GRS80 +units=m +no_defs
<5188> GRS80 Korea 2000 / East Sea Belt 2010 +proj=tmerc +lat_0=38 +lon_0=131 +k=1 +x_0=200000 +y_0=600000 +ellps=GRS80 +units=m +no_defs

51

한국 주요 좌표계 등록 전 활용
2097=PROJCS["Korean 1985 / Korea Central Belt",
GEOGCS["Korean 1985",
DATUM["Korean Datum 1985",
SPHEROID["Bessel 1841", 6377397.155, 299.1528128, AUTHORITY["EPSG","7004"]],
TOWGS84[-146.43, 507.89, 681.46, 0.0, 0.0, 0.0, 0.0], AUTHORITY["EPSG","6162"]],
PRIMEM["Greenwich", 0.0, AUTHORITY["EPSG","8901"]],
UNIT["degree", 0.017453292519943295],
AXIS["Geodetic longitude", EAST],
AXIS["Geodetic latitude", NORTH], AUTHORITY["EPSG","4162"]],
PROJECTION["Transverse_Mercator"],
PARAMETER["central_meridian", 127.00289027777001],
PARAMETER["latitude_of_origin", 38.0],
PARAMETER["scale_factor", 1.0],
PARAMETER["false_easting", 200000.0],
PARAMETER["false_northing", 500000.0],
UNIT["m", 1.0],
AXIS["Easting", EAST],
AXIS["Northing", NORTH],
AUTHORITY["EPSG","2097"]]
52

한국 주요 좌표계 등록 전 활용
DELETE FROM spatial_ref_sys WHERE srid =2097;

INSERT INTO spatial_ref_sys (srid, auth_name, auth_srid, srtext, proj4text)
VALUES (2097, 'EPSG', 2097,
'PROJCS["Korean 1985 / Central Belt", GEOGCS["Korean 1985", DATUM["Korean Datum 1985",
SPHEROID["Bessel 1841", 6377397.155, 299.1528128, AUTHORITY["EPSG","7004"]], TOWGS84[-
146.43, 507.89, 681.46, 0.0, 0.0, 0.0, 0.0], AUTHORITY["EPSG","6162"]],
PRIMEM["Greenwich", 0.0, AUTHORITY["EPSG","8901"]], UNIT["degree", 0.017453292519943295],
AXIS["Geodetic longitude", EAST], AXIS["Geodetic latitude", NORTH], AUTHORITY["EPSG","4162"]],
PROJECTION["Transverse_Mercator", AUTHORITY["EPSG","9807"]], PARAMETER["central_meridian",
127.00289027777001], PARAMETER["latitude_of_origin", 38.0], PARAMETER["scale_factor", 1.0],
PARAMETER["false_easting", 200000.0], PARAMETER["false_northing", 500000.0], UNIT["m", 1.0],
AXIS["Easting", EAST], AXIS["Northing", NORTH], AUTHORITY["EPSG","2097"]]',
'+proj=tmerc +lat_0=38 +lon_0=127 +k=1 +x_0=200000 +y_0=500000 +towgs84=-
146.43,507.89,681.46 +ellps=bessel +units=m +no_defs');

53

앞으로는 이렇게…
샘플데이터셋 폴더의 postgis_korea_epsg_towgs84.sql에 정리되어 있음

Ex)EPSG:5174

INSERT INTO spatial_ref_sys (srid, auth_name, auth_srid, srtext, proj4text)
VALUES (5174, 'EPSG', 5174,
'PROJCS["Korean 1985 / Modified Central Belt", GEOGCS["Korean 1985", DATUM["Korean Datum
1985", SPHEROID["Bessel
1841", 6377397.155, 299.1528128, AUTHORITY["EPSG","7004"]], TOWGS84[-
146.43, 507.89, 681.46, 0.0, 0.0, 0.0, 0.0], AUTHORITY["EPSG","6162"]], PRIMEM["Greenwic
h", 0.0, AUTHORITY["EPSG","8901"]], UNIT["degree", 0.017453292519943295], AXIS["Geodetic
longitude", EAST], AXIS["Geodetic
latitude", NORTH], AUTHORITY["EPSG","4162"]], PROJECTION["Transverse_Mercator", AUTHORITY["E
PSG","9807"]], PARAMETER["central_meridian", 127.00289027777775], PARAMETER["latitude_of_o
rigin", 38.0], PARAMETER["scale_factor", 1.0], PARAMETER["false_easting", 200000.0], PARAMETER["f
alse_northing", 500000.0], UNIT["m", 1.0], AXIS["Easting", EAST], AXIS["Northing", NORTH], AUTHOR
ITY["EPSG","5174"]]',
'+proj=tmerc +lat_0=38 +lon_0=127.00289 +k=1 +x_0=200000 +y_0=500000 +towgs84=-
146.43,507.89,681.46 +ellps=bessel +units=m +no_defs');

54

 6.1 GIS Objects
 6.2 PostGIS Functions

55

 4.1.1. OpenGIS WKB and WKT
 4.1.2. PostGIS EWKB, EWKT and Canonical Forms
 4.1.3. SQL-MM Part 3

 The GIS objects supported by PostGIS are a superset of the "Simple
Features" defined by the OpenGIS Consortium (OGC). As of version 0.9,
PostGIS supports all the objects and functions specified in the OGC
"Simple Features for SQL" specification.
 PostGIS extends the standard with support for 3DZ,3DM and 4D
coordinates.

56

 POINT(0 0)
 LINESTRING(0 0,1 1,1 2)
 POLYGON((0 0,4 0,4 4,0 4,0 0),(1 1, 2 1, 2 2, 1 2,1 1))
 MULTIPOINT(0 0,1 2)
 MULTILINESTRING((0 0,1 1,1 2),(2 3,3 2,5 4))
 MULTIPOLYGON(((0 0,4 0,4 4,0 4,0 0),(1 1,2 1,2 2,1 2,1
1)), ((-1 -1,-1 -2,-2 -2,-2 -1,-1 -1)))
 GEOMETRYCOLLECTION(POINT(2 3),LINESTRING(2 3,3 4))

58

 POINT(0 0 0) -- XYZ
 SRID=32632;POINT(0 0) -- XY with SRID
 POINTM(0 0 0) -- XYM
 POINT(0 0 0 0) -- XYZM
 SRID=4326;MULTIPOINTM(0 0 0,1 2 1) -- XYM with SRID
 MULTILINESTRING((0 0 0,1 1 0,1 2 1),(2 3 1,3 2 1,5 4 1))
 POLYGON((0 0 0,4 0 0,4 4 0,0 4 0,0 0 0),(1 1 0,2 1 0,2 2 0,1
2 0,1 1 0))
 MULTIPOLYGON(((0 0 0,4 0 0,4 4 0,0 4 0,0 0 0),(1 1 0,2 1
0,2 2 0,1 2 0,1 1 0)),((-1 -1 0,-1 -2 0,-2 -2 0,-2 -1 0,-1 -1 0)))
 GEOMETRYCOLLECTIONM(POINTM(2 3 9), LINESTRINGM(2 3
4, 3 4 5))
59

 CIRCULARSTRING(0 0, 1 1, 1 0)
 CIRCULARSTRING(0 0, 4 0, 4 4, 0 4, 0 0)
 COMPOUNDCURVE(CIRCULARSTRING(0 0, 1 1, 1 0),(1 0, 0 1))
 CURVEPOLYGON(CIRCULARSTRING(0 0, 4 0, 4 4, 0 4, 0 0),(1 1, 3 3, 3 1,
1 1))
 CURVEPOLYGON(COMPOUNDCURVE(CIRCULARSTRING(0 0,2 0, 2 1, 2 3,
4 3),(4 3, 4 5, 1 4, 0 0)), CIRCULARSTRING(1.7 1, 1.4 0.4, 1.6 0.4, 1.6
0.5, 1.7 1) )
 MULTICURVE((0 0, 5 5),CIRCULARSTRING(4 0, 4 4, 8 4))
 MULTISURFACE(CURVEPOLYGON(CIRCULARSTRING(0 0, 4 0, 4 4, 0 4, 0
0),(1 1, 3 3, 3 1, 1 1)),((10 10, 14 12, 11 10, 10 10),(11 11, 11.5 11, 11
11.5, 11 11)))

60

윈도우  프로그램 OpenGeo Suite 2.5 
Documentation  PostGIS Documentation 61

 AddGeometryColumn
 테이블을 레이어로
 DropGeometryColumn
 레이어를 테이블로
 DropGeometryTable
 레이어 삭제
 PostGIS_Full_Version
 라이브러리 버전정보 확인
 UpdateGeometrySRID
 좌표체계 정의
 Populate_Geometry_Columns
 레이어 제약조건 확인 및 업데이트
 Probe_Geometry_Columns
 레이어 제약조건 확인 및 업데이트(View는 제외)

 http://postgis.org/docs/reference.html#Management_Functions
62

 ST_GeometryFromText ST_MakeBox2D
 ST_GeomFromGML ST_MakeLine
 ST_GeomFromKML
ST_MakeEnvelope
ST_MakePolygon
 ST_GeomFromText
ST_MakePoint
 ST_GeomFromWKB ST_Point = ST_MakePoint.

 http://postgis.org/docs/reference.html#Geometry_Constructors

63

 GeometryType  ST_IsSimple
 ST_Boundary  ST_IsValid
 ST_Dimension  ST_NPoints
 ST_EndPoint  ST_NRings
 ST_Envelope  ST_NumGeometries
 ST_ExteriorRing  ST_NumInteriorRings
 ST_GeometryN  ST_NumPoints
 ST_GeometryType  ST_PointN
 ST_InteriorRingN  ST_SRID
 ST_IsClosed  ST_StartPoint
 ST_IsEmpty  ST_X
 ST_IsRing  ST_Y

http://postgis.org/docs/reference.html#Geometry_Accessors
64

 ST_AddPoint ST_RotateY
 ST_Affine ST_RotateZ
 ST_LineMerge
ST_Scale
ST_SetPoint
 ST_RemovePoint
ST_SetSRID
 ST_Reverse ST_Transform
 ST_Rotate
 ST_RotateX

 http://postgis.org/docs/reference.html#Geometry_Editors

65

 ST_AsBinary ST_AsGeoJSON
 ST_AsEWKB ST_AsGML
ST_AsHEXEWKB
 ST_AsText
ST_AsKML
 ST_AsEWKT ST_AsSVG
ST_GeoHash

 http://postgis.org/docs/reference.html#Geometry_Outputs

66

 ST_Contains
ST_Equals
 ST_Covers ST_Intersects
 ST_CoveredBy ST_Overlaps
 ST_Crosses ST_Relate
 ST_Disjoint ST_Touches

 http://postgis.org/docs/reference.html#Spatial_Relationships_Measure
ments

67

 ST_Area ST_MaxDistance
 ST_Centroid ST_Dwithin
 ST_ClosestPoint
ST_Length
ST_LongestLine
 ST_Distance
ST_ShortestLine
 ST_HausdorffDistance

 http://postgis.org/docs/reference.html#Spatial_Relationships_Measure
ments

68

ST_Buffer
 ST_Difference ST_ConvexHull
 ST_Intersection ST_MinimumBoundingCircle
 ST_SymDifference
ST_Polygonize
 ST_Union
ST_Simplify
 ST_MemUnion ST_SimplifyPreserveTopology
ST_Collect

 http://postgis.org/docs/reference.html#Geometry_Processing

69

 ST_Accum ST_XMax
 ST_Estimated_Extent ST_XMin
ST_YMax
 ST_Expand
ST_YMin
 ST_Extent ST_ZMax
 ST_Extent3D ST_ZMin

 http://postgis.org/docs/reference.html#Miscellaneous_Functions

70

 7.1 테이블 관리
 7.2 공간함수 활용 실습

71

 공간테이블 생성하기
 공간인덱싱 생성하기
 공간테이블 삭제하기
 공간 뷰(View) 생성 후 등록하기
 SQL을 이용하여 공간테이블 생성 후 등록하기

72

 cityhall 이라는 공간 테이블을 생성하고 다음의 정보를 입력
 X=197977.3885 , Y=451601.1444, city_name =서울

 Function
 AddGeometryColumn( <schema_name>, <table_name>, <column_name>, <srid>, <type>,
<dimension> )
 AddGeometryColumn( <table_name>, <column_name>, <srid>, <type>, <dimension> )

CREATE TABLE cityhall (city_name VARCHAR);
SELECT AddGeometryColumn('cityhall', 'the_geom', 2097, 'POINT', 2 );

INSERT INTO cityhall (the_geom, city_name) VALUES
(ST_GeomFromText('POINT(197977.3885 451601.1444)', 2097), '서
울');

 http://postgis.org/docs/ch04.html#OpenGISWKBWKT

73

 위에서 생성한 cityhall 테이블에 공간인덱싱을 생성

 SQL
 CREATE INDEX [indexname] ON [tablename] USING GIST ( [geometryfield] );

CREATE INDEX spatial_cityhall_the_geom ON public.cityhall USING
GIST(the_geom);

 http://postgis.org/docs/ch04.html#id2628148
74

 위에서 생성한 cityhall 테이블을 삭제

 일반적인 과정
 Remove refs from geometry_columns table
 Remove table

 Function
 DropGeometryTable(varchar table_name);
 DropGeometryTable(varchar schema_name, varchar table_name);
 DropGeometryTable(varchar catalog_name, varchar schema_name,
varchar table_name);

SELECT DropGeometryTable('public', 'cityhall');

75

 읍면동 행정경계(admin_emd)를 이용하여 종로구에 해당하는 읍면동의 중
심점을 admin_emd_point View로 등록
 Function
 CREATE OR REPLACE VIEW, DROP VIEW
 ST_Centroid(geometry g1)
 Populate_Geometry_Columns()

CREATE OR REPLACE VIEW admin_emd_point AS
SELECT
gid, ST_Centroid(the_geom) as the_geom, sgg_nm,
emd_cd, emd_nm, pop2008, pop_den
FROM admin_emd
WHERE sgg_nm = '종로구';

SELECT Populate_Geometry_Columns();
76

 시군구 행정경계(admin_sgg)를 이용하여 종로구에 해당하는 시군구 폴리
곤과 Intersects되는 도로(road_link2)를 Clip하여 새로운 테이블로 생성

CREATE TABLE clip_road_link AS
SELECT
ST_Intersection(r.the_geom, a.the_geom) as the_geom,
ST_Length(r.the_geom) as new_length,
a.sgg_cd as sgg_cd, a.sgg_nm as sgg_nm, r.road_name as
road_name, r.link_id as link_id
FROM
road_link2 as r,
admin_sgg as a
WHERE a.sgg_nm = '종로구' AND ST_Intersects(r.the_geom,
a.the_geom);

SELECT Populate_Geometry_Columns();; 77

 Measurement
 Proximity & Distance
 Spatial Join
 Overlay
 Projection
 Advanced Exercises

79

 1. 서울특별시 중구의 행정경계(admin_sgg) 면
적은 얼마인가? 단위는 ㎢로 표현하시오.

SELECT ST_Area(the_geom) / 1000000 AS area
FROM admin_sgg
WHERE sgg_nm = '중구'

80

 2. 도로 레이어(road_link2) 중 6차선(lanes) 이
상의 도로의 길이는 얼마인가? 단. ㎞로 표현하
시오.

SELECT SUM(ST_Length(the_geom)) / 1000 AS
km_roads
FROM road_link2
WHERE lanes >= 6
81

 3. 서울시의 읍면동(admin_emd) 중 면적이 가
장 큰 상위 10개를 선택하여 읍면동이름 및 면
적값을 조회하시오.

SELECT emd_nm, ST_Area(the_geom) as area
FROM admin_emd
ORDER BY area DESC
LIMIT 10
82

 4. 강(river) 레이어 중 hole을 가진 강은?

SELECT *
FROM river
WHERE ST_NumInteriorRings(the_geom) > 0

83

 1. 좌표 X = 197215 Y = 447711 지점에서 1KM
반경 내에 있는 대형매장(stores)은 무엇인가?
 ST_Distance, ST_Dwithin

SELECT *
FROM stores
WHERE
ST_Distance(the_geom, ST_GeomFromTex
t('POINT(197215 447711)', 2097)) < 1000
84

 2. 좌표 X = 197215 Y = 447711 지점에서 1KM
반경 내에 있는 대형매장(stores)은 무엇인가?

SELECT *
FROM stores
WHERE ST_DWithin(the_geom,
ST_GeomFromText('POINT(197215 447711)',
2097), 1000)
85

 2. 좌표 X = 197215 Y = 447711 지점에서 1KM 반
경 내에 있는 대형매장(stores)은 무엇인가?
 단. ST_Buffer, ST_Intersects 함수를 사용해서 구하
시오.

SELECT *
FROM stores
WHERE ST_Intersects(the_geom,
ST_Buffer(ST_GeomFromText('POINT(197215
447711)', 2097), 1000))

86

 1. 한강(river, river_cd = '1')과 인접한 서울시의
구(admin_sgg)는?

SELECT a.sgg_nm
FROM admin_sgg As a, river r
WHERE r.river_cd = '1' AND
ST_Intersects(a.the_geom, r.the_geom)
GROUP BY a.sgg_nm
87

 2. 영등포구(admin_sgg)에 속한 대형매장
(stores)은 무엇인가?

SELECT *
FROM stores s, admin_sgg a
WHERE a.sgg_nm = '영등포구' AND
ST_Within(s.the_geom, a.the_geom)

88

 3. 소방서(firestation)에서 500미터 반경 내의 도로
중 가장 가까운 도로(road_link2)와 거리를 계산하
시오.

SELECT DISTINCT ON(f.nam) f.nam, r.roadname_a,
r.lanes, ST_Distance(r.the_geom, f.the_geom)
As dist
FROM firestation AS f LEFT JOIN road_link2 As r ON
ST_DWithin(r.the_geom, f.the_geom, 500)
ORDER BY f.nam, ST_Distance(r.the_geom,
f.the_geom)
89

 1. 서울시 구별(admin_sgg) 6차선 이상 도로(road_link2, lanes)의 길
이를 구하시오.

SELECT
a.sgg_nm as sgg_nm,
SUM(ST_Length(ST_Intersection(r.the_geom,
a.the_geom))) / 1000 as road_length
FROM
road_link2 as r,
admin_sgg as a
WHERE r.lanes >= 6 AND ST_Intersects(r.the_geom,
a.the_geom)
GROUP BY sgg_nm
ORDER BY road_length DESC
90

 1. 서울특별시 중구의 중심점은 어디인가?
단, 경위도로 표시하라

SELECT
ST_AsText(ST_Transform(ST_Centroid(the_
geom), 4326))
FROM admin_sgg
WHERE sgg_nm = '중구'
91

 8.1 Backup
 8.2 Restore
 8.3 주의사항

92

•PostgreSQL 9.0 버전 이후는 백업 옵션이 강화됨

93

 절차
 템플릿 데이터베이스 생성
 Restore(되돌리기)

94

 상위버전과 하위버전 호환안됨
 하위버전 백업  상위버전 Restore 가능
 상위버전 백업  하위버전 Restore 오류

  버전충돌시 테이블/레이어 분리백업/복
원

95

 9.1 pgRuting

96

 pgRouting 은 PostGIS/PostgreSQL에 라우팅 기능을 추한 PostGIS의
Extension
 Camptocamp SA 에 의해 시작된 pgDijkstra의 확장판이며 Orkney 에
의해 확장
 현재 Georepublic 에 의해 개발 및 유지 관리되고 있음
 pgRouting는 GPLv2 라이선스 하에 사용 가능하며 개인, 기업 및 조직
의 커뮤니티에 의해 지원

 데이터베이스 기반 라우팅 접근 방법의 장점
 데이터와 속성은 Quantum GIS, uDig 그리고 JDBC, ODBC 또는
Pl/pgSQL 등 다양한 클라이언트에 의해 직접 갱신 가능하며, 클라이언트
는 PC나 모바일 기기가 될 수 있음
 데이터 변경 내용은 라우팅 엔진을 통해 즉시 반영되며, 재계산이 필요
없음
 ―cost‖ 파라미터는 SQL을 통해 동적으로 계산되고 그 값은 다중 필드 및
테이블의 다양한 값을 사용할 수 있음

97

 Shortest Path Dijkstra
 휴리스틱 기법(heuristic method)을 사용하지 않는 최단거리 알고리즘, 알고리즘을 개
발한 Dr. Edsger Wybe Dijkstra 교수의 이름을 사용.

 Shortest Path A-Star(A*)
 휴리스틱 기법을 사용하며 대용량 데이터셋에 적합한 최단거리 알고리즘

 Shortest Path Shooting-Star(Shooting*)
 휴리스틱 기법을 사용하며 turn restrictions(U-turn, P-turn, left-turn 등), 신호등, 편도
차선 등의 시제 도로 네트워크를 위한 최단거리 알고리즘

 Traveling Salesperson Problem (TSP)
 최대 40개의 포인트를 지원하는 Traveling Salesman Problem(TSP, 외판원 문제) 알고리
즘

 Driving Distance calculation (Isolines)
 도로네트워크에 기반을 둔 특정지점에서 특정 시간내에 도달 할 수 있는 영역. Service
area

98

 도로 네트워크 데이터셋은 최소한 다음의 정보를
포함

 Road link ID (gid)
 Road class (class_id): primary roads, secondary roads,
and local roads 등 Hierarchies
 Road link length (length): Cost - Length, Travel Time
등
 Road name (name)
 Road geometry (the_geom)
 Road restrictions & rule (optional): Lanes, Speed Limit,
traffic light, one-way streets 등

99

 공공기관에서 생산하는 도로네트워크로 활용 가능한
데이터는 다음의 4가지가 가능
 1. ITS 전국표준노드링크: http://nodelink.its.go.kr/
 2. UTIS 통합노드링크: http://www.utis.go.kr/
 3. 국가교통DB센터 노드링크: http://www.ktdb.go.kr/
 4. 새주소 도로구간
 1. 2. 3. 데이터는 1:5000 축척을 기반으로 제작되었으
며 주요 도로(고속국도/일반국도/지방도 수준) 중심으
로 상세한 지역에서의 활용은 어려움
 4. 새주소 도로구간은 도로중심선으로 구축되어 있으
나 네트워크 분석자료로 활용하기 위해서는 위상구조
편집이 필요하며 도로에 대한 상세정보(Restrictions &
rules)는 없음

100

사례지역: 제주도, 배경은 주거지(건물) 밀집도

ITS 전국표준 노드링크
새주소 도로중심선

101

 PostgreSQL 8.4 버전을 대상으로 하며 2012년 8월 현재 Windows에서
PostgreSQL 9.0 이상은 지원하지 않음
 설치에 대한 더 자세한 내용은 다음을 참고
 PgRouting 1.02 on Win32
 Installing PostgreSQL 8.4, PostGIS 1.4.1, and pgRouting 1.0.3 on Ubuntu 9.10

 다운로드
 http://www.pgrouting.org/download.html 이동
 Source Packages에서 pgrouting-1.05.tar.gz 다운로드
 Windows Binaries에서 pgRouting-1.03_pg-8.4.2.zip 다운로드
 특정폴더(C:OpenGeoSuitepgrouting)에 pgRouting-1.03_pg-8.4.2.zip 파일 압축 해제
 doc : 설치 도움말 포함
 lib : pgRouting 라이브러리
 share : 설치 sql
 share 폴더에는 Shortest Path(Dijkstra, A-Star, Shooting-Star) 알고리즘만 포함
 위 3 폴더를 C:Program FilesOpenGeoOpenGeo Suitepgsql8.4 폴더에 복사

 pgrouting-1.05.tar.gz 파일 열어 core/sql/matching.sql 압축풀기
 matching.sql은 Utility 함수를 포함하고 있음

102

 template_routing 이라는 라우팅 템플릿 데이터베이스를
생성
 Query 윈도우에서 다음을 실행

CREATE DATABASE template_routing
TEMPLATE=template_postgis

103

 PgAdmin III를 사용하는 경우 각 SQL Editor에서 각각의 sql 파일들을 불러와서 실행

 Core function
 psql -U postgres -d template_routing -f routing_core.sql
 psql -U postgres -d template_routing -f routing_core_wrappers.sql
 psql -U postgres -d template_routing -f routing_topology.sql
 psql -U postgres -d template_routing -f matching.sql
 TSP
 psql -U postgres -d template_routing -f routing_tsp.sql
 psql -U postgres -d template_routing -f routing_tsp_wrappers.sql
 Driving Distance
 psql -U postgres -d template_routing -f routing_dd.sql
 psql -U postgres -d template_routing -f routing_dd_wrappers.sql

 Optimize template database
 VACUUM FULL;
 VACUUM FREEZE;
 UPDATE pg_database SET datistemplate='true' WHERE datname='template_routing';

 UPDATE pg_database SET datallowconn='false' WHERE datname='template_routing';

104

 다음은 라우팅 데이터베이스 생성 후 도로 shapefile을 이용해 네트워크
데이터를 생성하는 과정
 Routing 데이터베이스 생성
CREATE DATABASE routing
TEMPLATE=template_routing

 샘플 데이터 변환
 샘플 데이터셋의 road_link2.shp를 road_link로 변환

 토폴로지 생성
 ALTER TABLE road_link ADD COLUMN source integer;
 ALTER TABLE road_link ADD COLUMN target integer;
 SELECT assign_vertex_id('road_link', 0.1, 'the_geom', 'gid');

 CREATE INDEX road_link_source_idx ON road_link("source");
 CREATE INDEX road_link_target_idx ON road_link("target");

105

 Dijkstra 알고리즘은 pgRouting에 구현된 최초의 알
고리즘
 source 와 target ID, id 속성과 cost 이외의 속성을
필요로 하지 않음
 directed와 undirected 그래프를 구분할 수 있으며,
네트워크의 역방향 cost(reverse cost) 유무를 설정
할 수 있음
 Cost를 도로의 길이로 사용하는 경우
 ALTER TABLE road_link ADD COLUMN length double precision;
 UPDATE road_link SET length = ST_Length(the_geom);

 ALTER TABLE road_link ADD COLUMN reverse_cost double precision;
 UPDATE road_link SET reverse_cost = length;

106

Dijkstra shortest_path 함수 선언
CREATE OR REPLACE FUNCTION shortest_path(
sql text,
source_id integer,
target_id integer,
directed boolean,
has_reverse_cost boolean)
RETURNS SETOF path_result

 sql: 다음의 컬럼을 포함한 결과를 생성하는 SQL query
 SELECT id, source, target, cost FROM edge_table
 * id: integer, edge 테이블의 identifier, 일반적으로 gid
 * source: integer, source vertex의 identifier
 * target: integer, target vertex의 identifier
 * cost: float8 value, 경로 탐색 비용(cost)에 사용되는 edge의 속성(-값의 cost를 가지는 edge는 그
래프 탐색에서 제외).
 * reverse_cost (옵션): 역방향 탐색(reverse traversal)에 사용되는 edge의 속성, 방향성이 true이고
has_reverse_cost가 true일 경우에만 사용
 source_id: integer 경로탐색 시점의 노드 id
 target_id: integer 경로탐색 종점의 노드 id
 directed: 그래프의 방향성이 정해진 경우 true
 has_reverse_cost: true이면 역방향 탐색에 SQL의 reverse_cost 컬럼을 사용

107

 다음은 네트워크의 길이(length)를 cost로 사용하며, 방향성을 무시하여 노드 id
3에서 7까지의 경로를 탐색하는 예
SELECT *
FROM shortest_path('
SELECT
gid AS id,
source::integer,
target::integer,
length::double precision AS cost
FROM road_link',
3,
7,
false,
false);

 Wrapper WITHOUT bounding box
SELECT gid AS edge_id, ST_Length(the_geom) AS cost FROM dijkstra_sp('road_link', 3, 7);

 Wrapper WITH bounding box
SELECT gid AS edge_id, ST_Length(the_geom) AS cost FROM dijkstra_sp_delta('road_link', 3, 7, 500);

108

 QIS에서 확인하기
 QGIS를 실행합니다.
 Add PostGIS Layer -> PosotGIS 연결정보 생성 후 road_link 및
vertices_tmp 레이어를 배경 레이어로 추가합니다.
 Add PostGIS Layer -> road_link 레이어를 선택한 후 Build query
버튼을 눌러 다음을 SQL where clause에 추가합니다.
 ―gid‖ IN (SELECT gid FROM dijkstra_sp('road_link', 3, 7))

109

 유용한 함수 – matching.sql

 POINT(211894.692 447070.0273)에서 50미터 반경 내 가장 가까운 링
크의 soruce, target 노드 정보를 검색
 SELECT gid, source, target
 FROM road_link
 WHERE gid IN (SELECT
find_nearest_link_within_distance(ST_AsText(ST_GeomFromText('P
OINT(211894.692 447070.0273)', 2097)), 300, 'road_link'));

 POINT(211894.692 447070.0273)에서 50미터 반경 내 가장 가까운 노
드의 ID와 좌표를 검색
 SELECT id, ST_AsText(the_geom) AS wkt
 FROM vertices_tmp
 WHERE id IN (SELECT
find_nearest_node_within_distance(ST_AsText(ST_GeomFromText('
POINT(211894.692 447070.0273)', 2097)), 50, 'road_link'))

110

 탐색 경로의 도로 정보 확인

 SELECT edge_id, vertex_id, cost, roadname_a, ST_AsText(the_geom) as
wkt
 FROM road_link
 JOIN
 (SELECT vertex_id, edge_id, cost FROM shortest_path('
 SELECT gid AS id,
 source::integer AS source,
 target::integer AS target,
 length::double precision AS cost
 FROM road_link',
 1385,
 2059,
 false,
 false)) AS route
 ON
 road_link.gid = route.edge_id;

111

 탐색 경로 주변 100미터 이내의 stores 찾기

 SELECT *
 FROM stores
 WHERE ST_Intersects(stores.the_geom,
 (SELECT ST_Buffer(ST_Union(road_link.the_geom), 2000)
 FROM shortest_path('SELECT gid as id,
 source::integer,
 target::integer,
 length::double precision as cost
 FROM road_link', 3, 7, false, false) AS sp
 INNER JOIN road_link ON road_link.gid = sp.edge_id));

112

 pgRouting Demo
 http://websi.openvrp.com
 http://map.veloland.ch
 http://www.ridethecity.com
 http://www.pgrouting.org/gallery.html

 Reference
 http://www.pgrouting.org/
 http://workshop.pgrouting.org/
 http://underdark.wordpress.com/2011/02/07/a-beginners-guide-to-
pgrouting/
 https://github.com/pgRouting/pgrouting/
 http://2011.foss4g.org/sites/2011.foss4g.org/files/slides/FOSS4G2011_
pgrouting2.pdf
 http://2010.foss4g.org/presentations/3284.pdf

113

 PostGIS - http://postgis.org/
 Docs - http://postgis.org/docs/
 PostgreSQL - http://www.postgresl.org/
 Downloads - http://www.postgresql.org/download/
 Docs - http://www.postgresql.org/docs/
 JDBC Driver - http://jdbc.postgresql.org/
 .Net Driver - http://npgsql.projects.postgresql.org/
 Python Driver - http://www.pygresql.org/
 C/C++ Driver - http://www.postgresql.org/docs/8.4/static/libpq.html
 PgAdmin III - http://www.pgadmin.org/
 Open Source Desktop Clients
 uDig - http://udig.refractions.net/
 QGIS - http://qgis.org/
 OpenJUMP - http://openjump.org/

114

PostGIS - National Education Center for GIS: Open Source GIS

Recommended

More Related Content

What's hot (20)

Viewers also liked (13)

Similar to PostGIS - National Education Center for GIS: Open Source GIS (20)

More from MinPa Lee (7)

PostGIS - National Education Center for GIS: Open Source GIS