【文章內(nèi)容簡介】 MM語句定義了一些插值幾何對象：CIRCULARSTRING(0 0, 1 1, 1 0) ——插值圓弧COMPOUNDCURVE(CIRCULARSTRING(0 0, 1 1, 1 0),(1 0, 0 1)) ——插值復(fù)合曲線CURVEPOLYGON(CIRCULARSTRING(0 0, 4 0, 4 4, 0 4, 0 0),(1 1, 3 3, 3 1, 1 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, 11 , 11 11))) ——多曲面 PostGIS中空間信息處理的實現(xiàn)1. spatial_ref_sys表在基于PostGIS模板創(chuàng)建的數(shù)據(jù)庫的public模式下，有一個spatial_ref_sys表，它存放的是OGC規(guī)范的空間參考。我們?nèi)∥覀冏钍煜さ?326參考看一下：它的srid存放的就是空間參考的WellKnown ID，對這個空間參考的定義主要包括兩個字段，srtext存放的是以字符串描述的空間參考， 投影定義（）。4326空間參考的srtext內(nèi)容：GEOGCS[WGS 84,DATUM[WGS_1984,SPHEROID[WGS 84,6378137,AUTHORITY[EPSG,7030]],TOWGS84[0,0,0,0,0,0,0],AUTHORITY[EPSG,6326]],PRIMEM[Greenwich,0,AUTHORITY[EPSG,8901]],UNIT[degree,AUTHORITY[EPSG,9122]],AUTHORITY[EPSG,4326]]4326空間參考的proj4text內(nèi)容：+proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs2. geometry_columns表geometry_columns表存放了當前數(shù)據(jù)庫中所有幾何字段的信息，比如我當前的庫里面有兩個空間表，在geometry_columns表中就可以找到這兩個空間表中幾何字段的定義：其中f_table_schema字段表示的是空間表所在的模式，f_table_name字段表示的是空間表的表名，f_geometry_column字段表示的是該空間表中幾何字段的名稱，srid字段表示的是該空間表的空間參考。3. 在PostGIS中創(chuàng)建一個空間表在PostGIS中創(chuàng)建一個包含幾何字段的空間表分為2步：第一步創(chuàng)建一個一般表，第二步給這個表添加幾何字段。以下先在test模式下創(chuàng)建一個名為cities的一般表：create table (id int4, name varchar(20))再給cities添加一個名為shape的幾何字段（二維點）：select AddGeometryColumn(39。test39。, 39。cities39。, 39。shape39。, 4326, 39。POINT39。, 2)4. PostGIS對幾何信息的檢查PostGIS可以檢查幾何信息的正確性，這主要是通過IsValid函數(shù)實現(xiàn)的。以下語句分辨檢查了2個幾何對象的正確性，顯然，(0, 0)點和(1,1)點可以構(gòu)成一條線，但是(0, 0)點和(0, 0)點則不能構(gòu)成，這個語句執(zhí)行以后的得出的結(jié)果是TRUE,FALSE。select IsValid(39。LINESTRING(0 0, 1 1)39。), IsValid(39。LINESTRING(0 0,0 0)39。)默認PostGIS并不會使用IsValid函數(shù)檢查用戶插入的新數(shù)據(jù)，因為這會消耗較多的CPU資源（特別是復(fù)雜的幾何對象）。當你需要使用這個功能的時候，你可以使用以下語句為表新建一個約束：ALTER TABLE citiesADD CONSTRAINT geometry_validCHECK (IsValid(shape))這時當我們往這個表試圖插入一個錯誤的空間對象的時候，會得到一個錯誤：INSERT INTO ( shape, name )VALUES ( GeomFromText(39。LINESTRING(0 0,0 0)39。, 4326), 39。北京39。)。ERROR: new row for relation cities violates check constraint geometry_validSQL 狀態(tài): 235145. PostGIS中的空間索引數(shù)據(jù)庫對多維數(shù)據(jù)的存取有兩種索引方案，RTree和GiST（Generalized Search Tree），在PostgreSQL中的GiST比RTree的健壯性更好，因此PostGIS對空間數(shù)據(jù)的索引一般采用GiST實現(xiàn)。以下的語句給sde模式中的cities表添加了一個空間索引shape_index_cities，在pgAdmin中也可以通過圖形界面完成相同的功能。CREATE INDEX shape_index_citiesON USING gist(shape)。另外要注意的是，空間索引只有在進行基于邊界范圍的查詢時才起作用，比如“amp。amp。”操作。 PostGIS中的常用函數(shù)以下內(nèi)容包括比較多的尖括號，發(fā)布到blogger的時候會顯示不正常，內(nèi)容太多我也無暇一個個手動改代碼，因此如有問題就去參考PostGIS官方文檔。首先需要說明一下，這里許多函數(shù)是以ST_[X]yyy形式命名的，事實上很多函數(shù)也可以通過xyyy的形式訪問，在PostGIS的函數(shù)庫中我們可以看到這兩種函數(shù)定義完全一樣。1. OGC標準函數(shù)管理函數(shù)：添加幾何字段 AddGeometryColumn(, , , , , )刪除幾何字段 DropGeometryColumn(, , )檢查數(shù)據(jù)庫幾何字段并在geometry_columns中歸檔 Probe_Geometry_Columns()給幾何對象設(shè)置空間參考（在通過一個范圍做空間查詢時常用） ST_SetSRID(geometry, integer)幾何對象關(guān)系函數(shù)：獲取兩個幾何對象間的距離 ST_Distance(geometry, geometry)如果兩個幾何對象間距離在給定值范圍內(nèi)，則返回TRUE ST_DWithin(geometry, geometry, float)判斷兩個幾何對象是否相等（比如LINESTRING(0 0, 2 2)和LINESTRING(0 0, 1 1, 2 2)是相同的幾何對象） ST_Equals(geometry, geometry)判斷兩個幾何對象是否分離 ST_Disjoint(geometry, geometry)判斷兩個幾何對象是否相交 ST_Intersects(geometry, geometry)判斷兩個幾何對象的邊緣是否接觸 ST_Touches(geometry, geometry)判斷兩個幾何對象是否互相穿過 ST_Crosses(geometry, geometry)判斷A是否被B包含 ST_Within(geometry A, geometry B)判斷兩個幾何對象是否是重疊 ST_Overlaps(geometry, geometry)判斷A是否包含B ST_Contains(geometry A, geometry B)判斷A是否覆蓋 B ST_Covers(geometry A, geometry B)判斷A是否被B所覆蓋 ST_CoveredBy(geometry A, geometry B)通過DE9IM 矩陣判斷兩個幾何對象的關(guān)系是否成立 ST_Relate(geometry, geometry, intersectionPatternMatrix)獲得兩個幾何對象的關(guān)系（DE9IM矩陣） ST_Relate(geometry, geometry)幾何對象處理函數(shù)：獲取幾何對象的中心 ST_Centroid(geometry)面積量測 ST_Area(geometry)長度量測 ST_Length(geometry)返回曲面上的一個點 ST_PointOnSurface(geometry)獲取邊界 ST_Boundary(geometry)獲取緩沖后的幾何對象 ST_Buffer(geometry, double, [integer])獲取多幾何對象的外接對象 ST_ConvexHull(geometry)獲取兩個幾何對象相交的部分 ST_Intersection(geometry, geometry)將經(jīng)度小于0的值加360使所有經(jīng)度值在0360間 ST_Shift_Longitude(geometry)獲取兩個幾何對象不相交的部分（A、B可互換） ST_SymDifference(geometry A, geometry B)從A去除和B相交的部分后返回 ST_Difference(geometry A, geometry B)返回兩個幾何對象的合并結(jié)果 ST_Union(geometry, geometry)返回一系列幾何對象的合并結(jié)果 ST_Union(geometry set)用較少的內(nèi)存和較長的時間完成合并操作，結(jié)果和ST_Union相同 ST_MemUnion(geometry set)幾何對象存取函數(shù)：獲取幾何對象的WKT描述 ST_AsText(geometry)獲取幾何對象的WKB描述 ST_AsBinary(geometry)獲取幾何對象的空間參考ID ST_SRID(geometry)獲取幾何對象的維數(shù) ST_Dimension(geometry)獲取幾何對象的邊界范圍 ST_Envelope(geometry)判斷幾何對象是否為空 ST_IsEmpty(geometry)判斷幾何對象是否不包含特殊點（比如自相交） ST_IsSimple(geometry)判斷幾何對象是否閉合 ST_IsClosed(geometry)判斷曲線是否閉合并且不包含特殊點 ST_IsRing(geometry)獲取多幾何對象中的對象個數(shù) ST_NumGeometries(geometry)獲取多幾何對象中第N個對象 ST_GeometryN(geometry,int)獲取幾何對象中的點個數(shù) ST_NumPoints(geometry)獲取幾何對象的第N個點 ST_PointN(geometry,integer)獲取多邊形的外邊緣 ST_ExteriorRing(geometry)獲取多邊形內(nèi)邊界個數(shù) ST_NumInteriorRings(geometry)同上 ST_NumInteriorRing(geometry)獲取多邊形的第N個內(nèi)邊界 ST_InteriorRingN(geometry,integer)獲取線的終點 ST_EndPoint(geometry)獲取線的起始點 ST_StartPoint(geometry)獲取幾何對象的類型 GeometryType(geometry)類似上，但是不檢查M值，即POINTM對象會被判斷為point ST_GeometryType(geometry)獲取點的X坐標 ST_X(geometry)獲取點的Y坐標 ST_Y(geometry)獲取點的Z坐標 ST_Z(geometry)獲取點的M值 ST_M(geometry)幾何對象構(gòu)造函數(shù)：參考語義：Text：WKTWKB：WKBGeom:GeometryM:MultiBd:BuildAreaColl:Collection ST_GeomFromText(text,[])ST_PointFromText(text,[])ST_LineFromText(text,[])ST_LinestringFromText(text,[])ST_PolyFromText(text,[])ST_PolygonFromText(text,[])ST_MPointFromText(text,[])ST_MLineFromText(text,[])ST_MPolyFromText(text,[])ST_GeomCollFromText(text,[])ST_GeomFromWKB(bytea,[])ST_GeometryFromWKB(bytea,[])ST_PointFromWKB(bytea,[])ST_LineFromWKB(bytea,[])ST_LinestringFromWKB(bytea,[])ST_PolyFromWKB(bytea,[])ST_PolygonFromWKB(bytea,[])ST_MPointFromWKB(bytea,[])ST_MLineFromWKB(bytea,[])ST_MPolyFromWKB(bytea,[])ST_GeomCollFromWKB(bytea,[])ST_BdPolyFromText(text WKT, integer SRID)ST_BdMPolyFromText(text WKT, integer SRID) 2. PostGIS擴展函數(shù)管理函數(shù)：刪除一個空間表（包括geometry_columns中的記錄） DropGeometryTable([], )更新空間表的空間參考 UpdateGeometrySRID([], , , )更新空間表的統(tǒng)計信息 update_geometry_stats([, ])參考語義：Geos：GEOS庫Jts：JTS庫Proj：PROJ4庫 postgis_version()post