PostGIS事實(shí)上算是筆者開始寫博客的第一篇內(nèi)容。而事實(shí)上那篇博文的內(nèi)容并不豐富，筆者對(duì)PostGIS的了解仍然不多，然而17年在OSGeo課程學(xué)習(xí)時(shí)對(duì)PostGIS又有了進(jìn)一步了解，并逐步發(fā)現(xiàn)它的強(qiáng)大。剛好最近又遇上一個(gè)問題，也使我萌發(fā)了開坑PostGIS的想法。
1 PostGIS簡(jiǎn)介
PostGIS是對(duì)象關(guān)系型數(shù)據(jù)庫系統(tǒng)PostgreSQL的一個(gè)擴(kuò)展，PostGIS提供如下空間信息服務(wù)功能:空間對(duì)象、空間索引、空間操作函數(shù)和空間操作符。同時(shí)，PostGIS遵循OpenGIS的規(guī)范。PostGIS的版權(quán)被納入到GNU的GPL中，也就是說任何人可以自由得到PostGIS的源碼并對(duì)其做研究和改進(jìn)。正是由于這一點(diǎn)，PostGIS得到了迅速的發(fā)展，越來越多的愛好者和研究機(jī)構(gòu)參與到PostGIS的應(yīng)用開發(fā)和完善當(dāng)中。
以上引自百度百科。
下面是我的第一篇博客。
安裝可以在網(wǎng)上搜索教程，在連接ArcGIS、QGIS拓展內(nèi)容可見上文。
2 關(guān)鍵問題
最近需要的一個(gè)工作是基于土地覆被數(shù)據(jù)和其他數(shù)據(jù)做一個(gè)簡(jiǎn)單的適宜用地提取。其實(shí)總結(jié)起來就是各種基礎(chǔ)的空間疊加分析，但是卻遇上了一個(gè)關(guān)鍵問題。土地覆被數(shù)據(jù)分辨率為300 m，下圖即為土地覆被數(shù)據(jù)（已提取了需要的土地覆被類型）。為了面積準(zhǔn)確性，進(jìn)行柵格轉(zhuǎn)矢量的時(shí)候并沒有選擇簡(jiǎn)化面。結(jié)果導(dǎo)致數(shù)據(jù)量非常巨大。
打開屬性表可以發(fā)現(xiàn)，研究區(qū)一共有3146163條數(shù)據(jù)。
然后需要和另外一個(gè)結(jié)果做空間疊加（Intersection，相交)，得到適宜用地的空間分布。結(jié)果由于數(shù)據(jù)量太大運(yùn)行時(shí)間較長(zhǎng)。在ArcGIS平臺(tái)測(cè)試結(jié)果如圖?；?1分鐘27秒。
筆者同時(shí)嘗試著使用了ArcGIS Pro來進(jìn)行大數(shù)據(jù)量的矢量相交計(jì)算。計(jì)算時(shí)間為13分鐘6秒。
也有人給的建議是選用PostGIS進(jìn)行計(jì)算。于是進(jìn)行了測(cè)試。
需要設(shè)置SRID和geometry。接下來在PostGIS的SQL查詢里運(yùn)行如下的SQL查詢。
SELECT ST_Intersection(suitablewgs84.geometry, ecooutwgs84.geometry)
FROM public.suitablewgs84 INNER JOIN public.ecooutwgs84 on ST_Intersects(suitablewgs84.geometry, ecooutwgs84.geometry)
可能是目前學(xué)習(xí)不夠深入，速度似乎不夠快。而這也是本部分學(xué)習(xí)筆記的開篇。
先放一些PostGIS的資料。
