PostGIS: Porovnání verzí

Z GeoWikiCZ
m (Obsah stránky nahrazen textem „{{freegiswiki|PostGIS}}“)
 
(Není zobrazeno 47 mezilehlých verzí od stejného uživatele.)
Řádek 1: Řádek 1:
[[Image:postgis.png|150px|right]]
{{freegiswiki|PostGIS}}
[http://www.postgis.org PostGIS] je rozšíření objektově-relačního databázového systému [[PostgreSQL]] pro podporu geografických objektů. PostGIS implementuje specifikaci [http://www.opengeospatial.org/standards/sfa Simple Features] konsorcia [http://www.opengeospatial.org Open Geospatial Consortium].
 
Pro experimenty používejte '''[[Cvičná databáze PostGIS|cvičnou databázi]]''' na serveru [http://josef.fsv.cvut.cz josef].
 
== Příklady ==
 
=== [[153GIS1]] ===
 
* [[153GIS1 - 2. cvičení - PostGIS|2. cvičení]]
* [[153GIS1 - 3. cvičení - PostGIS|3. cvičení]]
* [[153GIS1 - 4. cvičení - PostGIS|4. cvičení]]
 
=== [[153UZPD]] ===
 
* http://geo.fsv.cvut.cz/~gin/uzpd/examples/postgis/
 
== Vytvoření databáze ==
 
createdb <databáze>
createlang plpgsql <databáze>
psql -d <databáze> -f cesta/k/postgis.sql
psql -d <databáze> -f cesta/k/spatial_ref_sys.sql
 
Příklad založení PostGIS databáze:
 
<source lang=bash>
export DB=pgisďb
export PGPATH=/usr/share/postgresql/8.4/contrib/postgis-1.5/
 
createdb $DB
createlang plpgsql $DB
psql -d $DB -f $PGPATH/postgis.sql
psql -d $DB -f $PGPATH/spatial_ref_sys.sql
</source>
 
 
'''Poznámka:''' Pro nahrání dávky 'postgis.sql' jsou vyžadována práva "superuživatele" (superuser) databázového systému [[PostgreSQL#Poznámky|PostgreSQL]].
 
== Import dat ==
 
V případě, že jsou atributová v jiném kódování, je nutné nastavit před importem proměnnou prostředí <code>PGCLIENTENCODING</code>, např.
 
<source lang="bash">
export PGCLIENTENCODING=WIN1250
</source>
 
=== SQL ===
 
<source lang="sql">
INSERT INTO zb (cat, zb_geom) VALUES (1, ST_GeomFromText('POINT(13.30150173 49.79076912)', 4326));
</source>
 
nebo
 
<source lang="sql">
INSERT INTO zb (cat, zb_geom) VALUES (1, ST_SetSRID(ST_MakePoint(13.30150173,49.79076912), 4326));
</source>
 
Viz [http://www.postgis.org/documentation/manual-svn/ch07.html#Geometry_Constructors geometrické konstruktory].
 
=== ESRI Shapefile (shp2pgsql, ogr2ogr) ===
 
shp2pgsql -s 4326 -D -I -W latin2 cr.shp cr | psql pgis_yfsg
 
;OGR včetně transformace S-JTSK &rarr; WGS-84
 
ogr2ogr -f PostgreSQL -s_srs '+proj=krovak \
+a=6377397.155 +rf=299.1528128 +no_defs \
+towgs84=570.8,85.7,462.8,4.998,1.587,5.261,3.56 +to_meter=1.0' \
-t_srs EPSG:4326 \
pg:dbname=pgis_yfsg cr.shp
 
=== [[GRASS GIS]] (v.out.ogr) ===
 
{{GrassPrikaz|v.out.ogr}} input=cr type=area format=PostgreSQL dsn="pg:dbname=pgis_yfsg"
 
Příklad skriptu pro import vektorových vrstev z databáze 'pgis_student' schéma 'gis1'.
 
<source lang='python'>
#!/usr/bin/python                                                                                                                                                         
 
import sys
import psycopg2
import grass.script as grass
 
dbname = 'pgis_student'
schema = 'gis1'
 
def import_maps(maps):
    global dbname
    global schema
    for name in maps:
        grass.message("Importuji vektorovou mapu <%s>..." % name)
        grass.run_command('v.in.ogr',
                          quiet = True,
                          overwrite = True,
                          dsn = 'PG:dbname=%s' % dbname,
                          layer = schema + '.' + name,
                          output = name)
        grass.run_command('v.info',
                          flags = 't',
                          map = name)
def main():
    global dbname
    global schema
    conn = psycopg2.connect("dbname=%s" % dbname)
    cur = conn.cursor()
 
    cur.execute("SELECT tablename FROM pg_tables WHERE schemaname = '%s';" % schema)
    maps = list()
    for row in cur.fetchall():
        maps.append(row[0])
 
    cur.close()
    conn.close()
 
    import_maps(maps)
 
    return 0
 
if __name__ == "__main__":
    sys.exit(main())
</source>
 
=== Další formáty podporované knihovnou [[GDAL/OGR|OGR]] (ogr2ogr) ===
 
;GPX
 
ogr2ogr -f PostgreSQL -t_srs EPSG:4326 pg:dbname=pgis_yfsg yfsg-w.gpx waypoints
 
== Poznámky ==
 
=== Atribut geometrie ===
 
Pro manipulaci s atributem geometrie slouží specializované funkce, které aktualizují tabulku <tt>geometry_columns</tt>, vytvoří nad atributem prostorový index.
 
*;Vytvoření prostorového atributu
 
SELECT AddGeometryColumn(
schéma,                        # volitelný atribut
název tabulky,
sloupec,
srid,
typ geometrie,
dimenze);
 
Např.
<source lang="sql">
SELECT AddGeometryColumn('zb', 'zb_geom', 4326, 'POINT', 2);
</source>
 
*;Odstranění atributu geometrie
 
DropGeometryColumn(
schéma,                        # volitelný atribut
název tabulky,
sloupec)
 
*;Odstranění tabulky
 
DropGeometryTable(
schéma,                        # volitelný atribut
název tabulky)
 
Např.
 
<source lang="sql">
SELECT DropGeometryTable('osm', 'czech_roads');
</source>
 
Při manuálním vytvoření prostorového atributu (např. <code>CREATE TABLE tabulka AS ...</code>) je potřeba aktualizovat tabulku <tt>geometry_columns</tt> (INSERT nebo pomocí funkce <tt>Populate_Geometry_Columns</tt> (od verze 1.4.0)). Např.:
 
<source lang="sql">
SELECT Populate_Geometry_Columns('public.zb'::regclass);
</source>
 
=== Definice S-JTSK (ESRI:102067) ===
 
<source lang="sql">
INSERT into spatial_ref_sys (srid, auth_name, auth_srid, proj4text, srtext) values
( 102067, 'esri', 102067, '+proj=krovak +lat_0=49.5 +lon_0=24.83333333333333
+alpha=30.28813975277778 +k=0.9999 +x_0=0 +y_0=0 +ellps=bessel +units=m +no_defs ',
'PROJCS["S-JTSK_Krovak_East_North",GEOGCS["GCS_S_JTSK",
DATUM["Jednotne_Trigonometricke_Site_Katastralni",SPHEROID["Bessel_1841",6377397.155,299.1528128]],
PRIMEM["Greenwich",0],UNIT["Degree",0.017453292519943295]],PROJECTION["Krovak"],
PARAMETER["False_Easting",0],PARAMETER["False_Northing",0],PARAMETER["Pseudo_Standard_Parallel_1",78.5],
PARAMETER["Scale_Factor",0.9999],PARAMETER["Azimuth",30.28813975277778],
PARAMETER["Longitude_Of_Center",24.83333333333333],PARAMETER["Latitude_Of_Center",49.5],
PARAMETER["X_Scale",-1],PARAMETER["Y_Scale",1],PARAMETER["XY_Plane_Rotation",90],UNIT["Meter",1],
AUTHORITY["EPSG","102067"]]');
</source>
 
Viz [http://spatialreference.org/ref/esri/102067/ SpatialReference.org]
 
=== Rozsah dat ===
 
<source lang="sql">
SELECT ST_Extent(the_geom) from points;
-- EPSG:4326 --
SELECT ST_Extent(ST_Transform(the_geom, 4326)) from points;
</source>
 
=== Self-intersection ===
 
Příklad detekce průniku polygonů lesních porostů (data [[OpenStreetMap]]).
 
<source lang="sql">
CREATE VIEW lesy AS SELECT osm_id,way FROM czech_polygon WHERE landuse = 'forest';
 
CREATE AGGREGATE array_accum (
    sfunc = array_append,
    basetype = anyelement,
    stype = anyarray,
    initcond = '{}'
);
 
CREATE TABLE ilesy AS SELECT COUNT(DISTINCT g1.osm_id) AS count1, COUNT(DISTINCT g2.osm_id) AS count2,
        array_accum(DISTINCT g1.osm_id) AS array_accum1, array_accum(DISTINCT g2.osm_id) AS array_accum2,
        st_collect(DISTINCT g1.way) AS st_collect1,
        st_collect(DISTINCT g2.way) AS st_collect2,
        st_intersection(g1.way, g2.way), st_area(st_intersection(g1.way, g2.way))
FROM lesy g1, lesy g2
WHERE g1.osm_id < g2.osm_id AND st_intersects(g1.way, g2.way) AND
        st_isvalid(g1.way) AND st_isvalid(g2.way)
GROUP BY st_intersection(g1.way, g2.way)
ORDER BY COUNT(DISTINCT g1.osm_id), st_area(st_intersection(g1.way, g2.way));
</source>
 
== Procedury ==
 
Příklady uživatelských procedur pro PostGIS najdete např. [http://trac.osgeo.org/postgis/wiki/UsersWikiplpgsqlfunctions zde].
 
Příklad funkce pro výpočet výměry v km<sup>2</sup>.
 
<source lang="sql">
CREATE OR REPLACE FUNCTION area_km(geometry)
RETURNS integer AS
$$
SELECT ROUND(ST_Area($1)/1e6)::integer
$$
LANGUAGE 'sql';
</source>
 
Aplikace funkce <tt>area_km</tt>.
 
<source lang="sql">
SELECT nazev, area_km(the_geom) FROM gis1.obce limit 3;
</source>
 
<pre>
  nazev  | area_km
----------+---------
Abertamy |      9
Adamov  |      1
Adamov  |      3
(3 rows)
</pre>
 
Další příklad je převzat z [[153GIS1 - 3._cvičení - PostGIS#3.|3. cvičení GIS1]].
 
<source lang="sql">
CREATE OR REPLACE FUNCTION boundary_split(geometry, integer)
RETURNS SETOF geometry AS
$$
SELECT ST_Line_Substring(the_geom, $2 * n / length,
  CASE
WHEN $2 * (n + 1) < length THEN $2 * (n + 1) / length
ELSE 1
  END) AS the_geom
FROM
(SELECT ST_LineMerge($1) AS the_geom,
  ST_Length($1) AS length) AS t
CROSS JOIN generate_series(0, $2) AS n
WHERE n * $2 / length < 1;
$$
LANGUAGE 'sql';
</source>
 
Počet liniových segmentů pro délku linie 2km.
 
<source lang="sql">
CREATE TABLE cr AS SELECT ST_Boundary(ST_Union(the_geom)) AS the_geom FROM gis1.obce GROUP BY pomoc;
 
SELECT COUNT(*) FROM (SELECT boundary_split(the_geom, 2000) FROM cr) AS cr_2km;
</source>
 
<pre>
1059
</pre>
 
Počet liniových segmentů pro délku linie 5km.
 
<source lang="sql">
SELECT COUNT(*) FROM (SELECT boundary_split(the_geom, 5000) FROM cr) AS cr_5km;
</source>
 
<pre>
424
</pre>
 
Podívejte se také na:
 
* [http://www.postgres.cz/index.php/PL/pgSQL PL/pgSQL]
* [http://www.root.cz/clanky/jemny-uvod-do-plpgsql/ Jemný úvod do PL/pgSQL]
* [http://www.postgresql.org/docs/8.3/static/plpgsql.html PL/pgSQL - SQL Procedural Language]
 
== Související články ==
 
* [[Instalace PostGIS]]
* [[Cvičná databáze PostGIS]]
* [[Psql a Emacs]]
* [[WKTRaster]]
* [[SpatiaLite]]
* [[OpenStreetMap]]
* [[pgRouting]]
 
== Výuka ==
 
* [[153UZPD|Úvod do zpracování prostorových dat]]
* [[153YFSG|Free software GIS]]
 
== Externí odkazy ==
 
* [http://www.postgis.org/documentation/ PostGIS documentation]
* [http://www.foss4g2007.org/workshops/W-04/ Introduction to PostGIS] by Paul Ramsey
* [http://revenant.ca/www/postgis/workshop/index.html Introduction to PostGIS (OSGeo)] by Mark Leslie
* [http://www.mapbender.org/presentations/Spatial_Data_Management_Arnulf_Christl/Spatial_Data_Management_Arnulf_Christl.pdf Spatial Data Management]
* [http://edndoc.esri.com/arcsde/9.0/general_topics/understand_spatial_relations.htm Understanding spatial relations] by ESRI
* [http://www.bostongis.com/postgis_quickguide.bqg PostGIS Quick Guide - Cheatsheet]
* [http://grass.osgeo.org/wiki/Spatial_SQL Spatial SQL on GRASSWiki]
* [http://www.manning.com/obe/ PostGIS in action]
* [http://pgrouting.postlbs.org pgRouting]
* [http://spatialreference.org Spatial Reference]
* [http://trac.osgeo.org/postgis/wiki/UsersWikiPostgisTopology PostGIS Topology]
* [http://s3.opengeo.org/postgis-power.pdf Tips for the PostGIS Power User] (FOSS4G 2010 - Paul Ramsey)
* [http://workshops.opengeo.org/postgis-spatialdbtips/ Spatial database tips and tricks]
* [http://www.kappasys.ch/cms/index.php?id=23&L=5 PostGIS Versioning - pgVersion]
 
<!-- * [http://gis.zcu.cz/studium/ugi/referaty/05/PostGIS/index.html Představení projektu PostGIS] -->
* [http://geoinformatics.fsv.cvut.cz/gwiki/PostGIS_pro_v%C3%BDvoj%C3%A1%C5%99e PostGIS pro vývojáře]
* [http://www.openweekend.cz/slides/ow_2005/orlik_postgis.pdf Správa časoprostorových dat v prostředí PostgreSQL/PostGIS]
 
{{Databáze}}
{{GIS}}
{{GFOSS}}

Aktuální verze z 8. 1. 2013, 14:32

Stránky přesunuty na Free GIS Portál: http://geo.fsv.cvut.cz/freegis/PostGIS