PostgreSQL+PostGIS 的使用

from PostgreSQL+PostGIS 的使用

一、PostGIS中的幾何類型

PostGIS支持所有OGC規(guī)范的“Simple Features”類型昂拂,同時在此基礎(chǔ)上擴(kuò)展了對3DZ、3DM抛猖、4D坐標(biāo)的支持格侯。

1. OGC的WKB和WKT格式

OGC定義了兩種描述幾何對象的格式鼻听,分別是WKB(Well-Known Binary)和WKT(Well-Known Text)。

在SQL語句中联四,用以下的方式可以使用WKT格式定義幾何對象:

  • POINT(0 0) ——點(diǎn)
  • 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) ——多點(diǎn)
  • 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))) ——幾何集合

以下語句可以使用WKT格式插入一個點(diǎn)要素到一個表中撑碴,其中用到的GeomFromText等函數(shù)在后面會有詳細(xì)介紹:
INSERT INTO table ( SHAPE, NAME )
VALUES ( GeomFromText('POINT(116.39 39.9)', 4326), '北京');

2. EWKT、EWKB和Canonical格式

EWKT和EWKB相比OGC WKT和WKB格式主要的擴(kuò)展有3DZ朝墩、3DM醉拓、4D坐標(biāo)和內(nèi)嵌空間參考支持。

以下以EWKT語句定義了一些幾何對象:
POINT(0 0 0) ——3D點(diǎn)
SRID=32632;POINT(0 0) ——內(nèi)嵌空間參考的點(diǎn)
POINTM(0 0 0) ——帶M值的點(diǎn)
POINT(0 0 0 0) ——帶M值的3D點(diǎn)
SRID=4326;MULTIPOINTM(0 0 0,1 2 1) ——內(nèi)嵌空間參考的帶M值的多點(diǎn)

以下語句可以使用EWKT格式插入一個點(diǎn)要素到一個表中:
INSERT INTO table ( SHAPE, NAME )
VALUES ( GeomFromEWKT('SRID=4326;POINTM(116.39 39.9 10)'), '北京' )

Canonical格式是16進(jìn)制編碼的幾何對象收苏,直接用SQL語句查詢出來的就是這種格式廉嚼。

3. SQL-MM格式

SQL-MM格式定義了一些插值曲線,這些插值曲線和EWKT有點(diǎn)類似倒戏,也支持3DZ怠噪、3DM、4D坐標(biāo)杜跷,但是不支持嵌入空間參考傍念。

以下以SQL-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.5 11, 11 11.5, 11 11))) ——多曲面

二、 PostGIS中空間信息處理的實現(xiàn)

1. spatial_ref_sys表

在基于PostGIS模板創(chuàng)建的數(shù)據(jù)庫的public模式下葛闷,有一個spatial_ref_sys表憋槐,它存放的是OGC規(guī)范的空間參考。我們?nèi)∥覀冏钍煜さ?326參考看一下:

它的srid存放的就是空間參考的Well-Known ID淑趾,對這個空間參考的定義主要包括兩個字段阳仔,srtext存放的是以字符串描述的空間參考,proj4text存放的則是以字符串描述的PROJ.4 投影定義(PostGIS使用PROJ.4實現(xiàn)投影)扣泊。

4326空間參考的srtext內(nèi)容:
GEOGCS["WGS 84",DATUM["WGS_1984",SPHEROID["WGS 84",6378137,298.257223563,AUTHORITY["EPSG","7030"]],TOWGS84[0,0,0,0,0,0,0],AUTHORITY["EPSG","6326"]],PRIMEM["Greenwich",0,AUTHORITY["EPSG","8901"]],UNIT["degree",0.01745329251994328,AUTHORITY["EPSG","9122"]],AUTHORITY["EPSG","4326"]]

4326空間參考的proj4text內(nèi)容:
+proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs

2. geometry_columns表

geometry_columns表存放了當(dāng)前數(shù)據(jù)庫中所有幾何字段的信息近范,比如我當(dāng)前的庫里面有兩個空間表,在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 test.cities (id int4, name varchar(20))

再給cities添加一個名為shape的幾何字段(二維點(diǎn)):
select AddGeometryColumn('test', 'cities', 'shape', 4326, 'POINT', 2)

4. PostGIS對幾何信息的檢查

PostGIS可以檢查幾何信息的正確性沥匈,這主要是通過IsValid函數(shù)實現(xiàn)的蔗喂。
以下語句分辨檢查了2個幾何對象的正確性,顯然高帖,(0, 0)點(diǎn)和(1,1)點(diǎn)可以構(gòu)成一條線缰儿,但是(0, 0)點(diǎn)和(0, 0)點(diǎn)則不能構(gòu)成,這個語句執(zhí)行以后的得出的結(jié)果是TRUE,FALSE棋恼。

select IsValid('LINESTRING(0 0, 1 1)'), IsValid('LINESTRING(0 0,0 0)')
默認(rèn)PostGIS并不會使用IsValid函數(shù)檢查用戶插入的新數(shù)據(jù)返弹,因為這會消耗較多的CPU資源(特別是復(fù)雜的幾何對象)。當(dāng)你需要使用這個功能的時候爪飘,你可以使用以下語句為表新建一個約束:
ALTER TABLE cities
ADD CONSTRAINT geometry_valid
CHECK (IsValid(shape))

這時當(dāng)我們往這個表試圖插入一個錯誤的空間對象的時候义起,會得到一個錯誤:
INSERT INTO test.cities ( shape, name )
VALUES ( GeomFromText('LINESTRING(0 0,0 0)', 4326), '北京');

ERROR: new row for relation "cities" violates check constraint "geometry_valid"
SQL 狀態(tài): 23514

5. PostGIS中的空間索引

數(shù)據(jù)庫對多維數(shù)據(jù)的存取有兩種索引方案,R-Tree和GiST(Generalized Search Tree)师崎,在PostgreSQL中的GiST比R-Tree的健壯性更好默终,因此PostGIS對空間數(shù)據(jù)的索引一般采用GiST實現(xiàn)。

以下的語句給sde模式中的cities表添加了一個空間索引shape_index_cities犁罩,在pgAdmin中也可以通過圖形界面完成相同的功能齐蔽。
CREATE INDEX shape_index_cities
ON sde.cities
USING gist
(shape);

另外要注意的是,空間索引只有在進(jìn)行基于邊界范圍的查詢時才起作用床估,比如“&&”操作含滴。

三、 PostGIS中的常用函數(shù)

以下內(nèi)容包括比較多的尖括號丐巫,發(fā)布到blogger的時候會顯示不正常谈况,內(nèi)容太多我也無暇一個個手動改代碼,因此如有問題就去參考PostGIS官方文檔递胧。

首先需要說明一下碑韵,這里許多函數(shù)是以ST_[X]yyy形式命名的,事實上很多函數(shù)也可以通過xyyy的形式訪問缎脾,在PostGIS的函數(shù)庫中我們可以看到這兩種函數(shù)定義完全一樣祝闻。

1. OGC標(biāo)準(zhǔn)函數(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)
通過DE-9IM 矩陣判斷兩個幾何對象的關(guān)系是否成立 ST_Relate(geometry, geometry, intersectionPatternMatrix)
獲得兩個幾何對象的關(guān)系(DE-9IM矩陣) ST_Relate(geometry, geometry)

幾何對象處理函數(shù):
獲取幾何對象的中心 ST_Centroid(geometry)
面積量測 ST_Area(geometry)
長度量測 ST_Length(geometry)
返回曲面上的一個點(diǎn) ST_PointOnSurface(geometry)
獲取邊界 ST_Boundary(geometry)
獲取緩沖后的幾何對象 ST_Buffer(geometry, double, [integer])
獲取多幾何對象的外接對象 ST_ConvexHull(geometry)
獲取兩個幾何對象相交的部分 ST_Intersection(geometry, geometry)
將經(jīng)度小于0的值加360使所有經(jīng)度值在0-360間 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)
判斷幾何對象是否不包含特殊點(diǎn)(比如自相交) ST_IsSimple(geometry)
判斷幾何對象是否閉合 ST_IsClosed(geometry)
判斷曲線是否閉合并且不包含特殊點(diǎn) ST_IsRing(geometry)
獲取多幾何對象中的對象個數(shù) ST_NumGeometries(geometry)
獲取多幾何對象中第N個對象 ST_GeometryN(geometry,int)
獲取幾何對象中的點(diǎn)個數(shù) ST_NumPoints(geometry)
獲取幾何對象的第N個點(diǎn) ST_PointN(geometry,integer)
獲取多邊形的外邊緣 ST_ExteriorRing(geometry)
獲取多邊形內(nèi)邊界個數(shù) ST_NumInteriorRings(geometry)
同上 ST_NumInteriorRing(geometry)
獲取多邊形的第N個內(nèi)邊界 ST_InteriorRingN(geometry,integer)
獲取線的終點(diǎn) ST_EndPoint(geometry)
獲取線的起始點(diǎn) ST_StartPoint(geometry)
獲取幾何對象的類型 GeometryType(geometry)
類似上,但是不檢查M值辙纬,即POINTM對象會被判斷為point ST_GeometryType(geometry)
獲取點(diǎn)的X坐標(biāo) ST_X(geometry)
獲取點(diǎn)的Y坐標(biāo) ST_Y(geometry)
獲取點(diǎn)的Z坐標(biāo) ST_Z(geometry)
獲取點(diǎn)的M值 ST_M(geometry)

幾何對象構(gòu)造函數(shù) :
參考語義:
Text:WKT
WKB:WKB
Geom:Geometry
M:Multi
Bd:BuildArea
Coll: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)

?著作權(quán)歸作者所有,轉(zhuǎn)載或內(nèi)容合作請聯(lián)系作者
  • 序言:七十年代末耸袜,一起剝皮案震驚了整個濱河市,隨后出現(xiàn)的幾起案子牲平,更是在濱河造成了極大的恐慌堤框,老刑警劉巖,帶你破解...
    沈念sama閱讀 216,651評論 6 501
  • 序言:濱河連續(xù)發(fā)生了三起死亡事件纵柿,死亡現(xiàn)場離奇詭異蜈抓,居然都是意外死亡,警方通過查閱死者的電腦和手機(jī)昂儒,發(fā)現(xiàn)死者居然都...
    沈念sama閱讀 92,468評論 3 392
  • 文/潘曉璐 我一進(jìn)店門沟使,熙熙樓的掌柜王于貴愁眉苦臉地迎上來,“玉大人渊跋,你說我怎么就攤上這事腊嗡∽徘悖” “怎么了?”我有些...
    開封第一講書人閱讀 162,931評論 0 353
  • 文/不壞的土叔 我叫張陵燕少,是天一觀的道長卡者。 經(jīng)常有香客問我,道長客们,這世上最難降的妖魔是什么崇决? 我笑而不...
    開封第一講書人閱讀 58,218評論 1 292
  • 正文 為了忘掉前任,我火速辦了婚禮底挫,結(jié)果婚禮上恒傻,老公的妹妹穿的比我還像新娘。我一直安慰自己建邓,他們只是感情好盈厘,可當(dāng)我...
    茶點(diǎn)故事閱讀 67,234評論 6 388
  • 文/花漫 我一把揭開白布。 她就那樣靜靜地躺著官边,像睡著了一般扑庞。 火紅的嫁衣襯著肌膚如雪。 梳的紋絲不亂的頭發(fā)上拒逮,一...
    開封第一講書人閱讀 51,198評論 1 299
  • 那天罐氨,我揣著相機(jī)與錄音,去河邊找鬼滩援。 笑死栅隐,一個胖子當(dāng)著我的面吹牛,可吹牛的內(nèi)容都是我干的玩徊。 我是一名探鬼主播租悄,決...
    沈念sama閱讀 40,084評論 3 418
  • 文/蒼蘭香墨 我猛地睜開眼,長吁一口氣:“原來是場噩夢啊……” “哼恩袱!你這毒婦竟也來了泣棋?” 一聲冷哼從身側(cè)響起,我...
    開封第一講書人閱讀 38,926評論 0 274
  • 序言:老撾萬榮一對情侶失蹤畔塔,失蹤者是張志新(化名)和其女友劉穎潭辈,沒想到半個月后,有當(dāng)?shù)厝嗽跇淞掷锇l(fā)現(xiàn)了一具尸體澈吨,經(jīng)...
    沈念sama閱讀 45,341評論 1 311
  • 正文 獨(dú)居荒郊野嶺守林人離奇死亡把敢,尸身上長有42處帶血的膿包…… 初始之章·張勛 以下內(nèi)容為張勛視角 年9月15日...
    茶點(diǎn)故事閱讀 37,563評論 2 333
  • 正文 我和宋清朗相戀三年,在試婚紗的時候發(fā)現(xiàn)自己被綠了谅辣。 大學(xué)時的朋友給我發(fā)了我未婚夫和他白月光在一起吃飯的照片修赞。...
    茶點(diǎn)故事閱讀 39,731評論 1 348
  • 序言:一個原本活蹦亂跳的男人離奇死亡,死狀恐怖桑阶,靈堂內(nèi)的尸體忽然破棺而出柏副,到底是詐尸還是另有隱情勾邦,我是刑警寧澤,帶...
    沈念sama閱讀 35,430評論 5 343
  • 正文 年R本政府宣布割择,位于F島的核電站眷篇,受9級特大地震影響,放射性物質(zhì)發(fā)生泄漏锨推。R本人自食惡果不足惜铅歼,卻給世界環(huán)境...
    茶點(diǎn)故事閱讀 41,036評論 3 326
  • 文/蒙蒙 一公壤、第九天 我趴在偏房一處隱蔽的房頂上張望换可。 院中可真熱鬧,春花似錦厦幅、人聲如沸沾鳄。這莊子的主人今日做“春日...
    開封第一講書人閱讀 31,676評論 0 22
  • 文/蒼蘭香墨 我抬頭看了看天上的太陽译荞。三九已至,卻和暖如春休弃,著一層夾襖步出監(jiān)牢的瞬間吞歼,已是汗流浹背。 一陣腳步聲響...
    開封第一講書人閱讀 32,829評論 1 269
  • 我被黑心中介騙來泰國打工塔猾, 沒想到剛下飛機(jī)就差點(diǎn)兒被人妖公主榨干…… 1. 我叫王不留篙骡,地道東北人。 一個月前我還...
    沈念sama閱讀 47,743評論 2 368
  • 正文 我出身青樓丈甸,卻偏偏與公主長得像糯俗,于是被迫代替她去往敵國和親。 傳聞我的和親對象是個殘疾皇子睦擂,可洞房花燭夜當(dāng)晚...
    茶點(diǎn)故事閱讀 44,629評論 2 354

推薦閱讀更多精彩內(nèi)容