1. 從一個(gè)被問爛了的問題說起緩沖區(qū)到底在算什么做GIS這行十來年被新手問得最多的空間分析操作里緩沖區(qū)Buffer絕對排前三。很多人第一次聽到這個(gè)詞腦子里浮現(xiàn)的是化學(xué)實(shí)驗(yàn)室里的緩沖溶液或者視頻播放時(shí)的緩沖進(jìn)度條。放到地理信息里它的意思其實(shí)很樸素給定一個(gè)點(diǎn)、一條線或者一個(gè)面向外擴(kuò)展出指定距離的范圍形成一個(gè)新多邊形。這個(gè)新多邊形就是緩沖區(qū)英文叫buffer zone或者buffer polygon。聽起來簡單到不值一提但真正落到項(xiàng)目里坑多得能讓人懷疑人生。我見過太多人拿著ArcGIS或者QGIS點(diǎn)一下Buffer工具輸入個(gè)500米回車然后拿著結(jié)果去做后續(xù)分析最后發(fā)現(xiàn)結(jié)果偏了幾十米甚至幾百米或者本該合并的區(qū)域碎成了一地雞毛。問題出在哪出在坐標(biāo)系、投影、距離單位、融合參數(shù)、端點(diǎn)樣式這些細(xì)節(jié)上。緩沖區(qū)分析是空間分析里最基礎(chǔ)的操作之一但基礎(chǔ)不等于簡單它背后牽扯的是地理坐標(biāo)系與投影坐標(biāo)系的轉(zhuǎn)換邏輯、歐氏距離與測地線距離的差異、矢量幾何的拓?fù)涮幚磉@一整套東西。這篇文章想做的事情很明確把Buffer這個(gè)操作從“點(diǎn)一下按鈕”的層面拆解到“為什么這么點(diǎn)、參數(shù)怎么選、結(jié)果怎么驗(yàn)”的層面。不管你是剛接觸GIS的學(xué)生還是做了幾年但一直沒深究過緩沖區(qū)的從業(yè)者我都希望你看完之后能對這個(gè)小工具有一個(gè)全新的認(rèn)識(shí)。我會(huì)用QGIS和ArcGIS Pro兩個(gè)平臺(tái)的操作做對照穿插Python的geopandas和shapely代碼實(shí)現(xiàn)把參數(shù)選擇的邏輯、常見錯(cuò)誤的排查、性能優(yōu)化的技巧都講透。全文基于我這些年做城市規(guī)劃、選址分析、環(huán)境評估類項(xiàng)目時(shí)積累的實(shí)際經(jīng)驗(yàn)不是照搬幫助文檔。2. 緩沖區(qū)的核心邏輯距離這件事沒那么簡單2.1 緩沖區(qū)到底在幾何上做了什么先從幾何本質(zhì)說起。對一個(gè)點(diǎn)做緩沖區(qū)就是以該點(diǎn)為圓心、以指定距離為半徑畫一個(gè)圓用這個(gè)圓去近似表示“距離該點(diǎn)不超過指定距離的所有位置”。對一條線做緩沖區(qū)是以線上每一點(diǎn)為圓心畫圓所有這些圓的并集形成一個(gè)“香腸形”的多邊形。對一個(gè)面做緩沖區(qū)是面本身加上面邊界向外擴(kuò)展的部分。這里有個(gè)關(guān)鍵點(diǎn)很多人忽略緩沖區(qū)操作本質(zhì)上是在做“距離場”的等值線提取。你指定的緩沖距離就是這個(gè)距離場的閾值。所有到原始要素的距離小于等于這個(gè)閾值的位置構(gòu)成緩沖區(qū)內(nèi)部大于閾值的在外部。理解這一點(diǎn)后面很多參數(shù)選擇就順理成章了。在矢量數(shù)據(jù)模型下緩沖區(qū)是通過幾何運(yùn)算生成的。點(diǎn)緩沖生成的是正多邊形用有限個(gè)頂點(diǎn)去逼近圓線緩沖生成的是帶圓角或方角的帶狀多邊形面緩沖生成的是原面加上外擴(kuò)環(huán)。在柵格數(shù)據(jù)模型下緩沖區(qū)可以通過歐氏距離工具生成距離柵格再做閾值分割得到。兩種路線各有適用場景后面會(huì)展開講。2.2 投影坐標(biāo)系緩沖區(qū)分析的第一道生死線我敢說新手做緩沖區(qū)出的問題八成以上跟坐標(biāo)系有關(guān)。核心邏輯是這樣的緩沖區(qū)距離是一個(gè)長度量而長度量只有在投影坐標(biāo)系下才有穩(wěn)定的物理意義。地理坐標(biāo)系比如WGS84EPSG:4326的單位是度不是米。你在這種坐標(biāo)系下輸入“500”軟件會(huì)理解為500度而不是500米。500度是什么概念地球一圈才360度這個(gè)緩沖區(qū)能把整個(gè)地球包進(jìn)去還多。所以做緩沖區(qū)之前第一件事永遠(yuǎn)是確認(rèn)數(shù)據(jù)在什么坐標(biāo)系下。如果是地理坐標(biāo)系必須先投影到合適的投影坐標(biāo)系。投影坐標(biāo)系的選擇要看研究區(qū)域的位置和范圍。國內(nèi)常用的有UTM各帶比如EPSG:32650對應(yīng)UTM 50N、高斯-克呂格各帶比如EPSG:4547對應(yīng)CGCS2000 3度帶中央經(jīng)線114度或者Web MercatorEPSG:3857但這個(gè)在高緯度地區(qū)面積變形很大做緩沖區(qū)要謹(jǐn)慎。注意Web Mercator在赤道附近1米就是1米但到了北緯40度左右1米的實(shí)際地面距離在投影上大約是1.3米。如果你在北京做緩沖區(qū)用EPSG:3857輸入500米實(shí)際地面距離只有約380米。這個(gè)誤差在選址分析里是致命的。選投影坐標(biāo)系的原則是研究區(qū)域盡量靠近投影帶的中央經(jīng)線范圍不要跨帶。如果研究區(qū)域跨了兩個(gè)UTM帶要么分帶處理要么選一個(gè)覆蓋范圍更大的投影比如Albers等面積投影但要接受一定的距離變形。我個(gè)人的經(jīng)驗(yàn)是城市級別的分析用當(dāng)?shù)氐母咚?克呂格3度帶投影最穩(wěn)省級以上用Albers跨國或全球分析才考慮用測地線緩沖。2.3 歐氏距離、測地線距離與曼哈頓距離緩沖區(qū)默認(rèn)用的是歐氏距離也就是直線距離。但在實(shí)際地理空間中兩點(diǎn)之間的“距離”可能有多種定義。測地線距離是球面上兩點(diǎn)之間的最短弧長適合大范圍分析。曼哈頓距離是沿網(wǎng)格線走的距離適合城市路網(wǎng)分析。緩沖區(qū)工具通常只支持歐氏距離如果你需要其他距離度量得用別的方法。QGIS的Buffer工具在投影坐標(biāo)系下用歐氏距離在地理坐標(biāo)系下會(huì)自動(dòng)切換到測地線計(jì)算較新版本支持。ArcGIS Pro的Buffer工具有一個(gè)“Method”參數(shù)可以選Planar平面或Geodesic測地線。選Geodesic時(shí)即使數(shù)據(jù)是地理坐標(biāo)系也能算出正確的米制距離但計(jì)算速度會(huì)慢一些而且生成的幾何在投影后可能會(huì)有輕微變形。我的建議是能投影就投影投影后用Planar方法又快又準(zhǔn)。只有在數(shù)據(jù)范圍極大比如跨洲或者實(shí)在不方便投影時(shí)才用Geodesic。測地線緩沖在極地附近和跨越日期變更線時(shí)尤其要注意幾何可能出問題。3. 參數(shù)詳解每個(gè)選項(xiàng)背后都是一道選擇題3.1 緩沖距離單值、多值與字段驅(qū)動(dòng)最簡單的緩沖區(qū)是固定距離比如所有要素都緩沖500米。但實(shí)際項(xiàng)目里不同要素的緩沖距離往往不同。比如做學(xué)校選址小學(xué)的服務(wù)半徑可能是800米中學(xué)是1500米這時(shí)候就需要用字段驅(qū)動(dòng)緩沖。QGIS的Buffer工具里Distance參數(shù)可以輸入固定值也可以點(diǎn)右邊的按鈕選擇字段。ArcGIS Pro的Buffer工具里Distance參數(shù)支持輸入字段名也支持輸入表達(dá)式。Python的geopandas里buffer()方法接受標(biāo)量或數(shù)組數(shù)組長度必須和幾何數(shù)量一致。多環(huán)緩沖區(qū)是另一個(gè)常見需求。比如做商業(yè)選址想同時(shí)看300米、500米、1000米三個(gè)范圍。QGIS里可以用“Multi-ring buffer”工具ArcGIS Pro里用“Multiple Ring Buffer”工具。底層邏輯是多次調(diào)用單環(huán)緩沖然后把結(jié)果合并到一個(gè)圖層里用distance字段區(qū)分。注意多環(huán)緩沖默認(rèn)是每個(gè)環(huán)獨(dú)立生成環(huán)與環(huán)之間會(huì)有重疊如果你想要的是“甜甜圈”形狀外環(huán)減去內(nèi)環(huán)需要額外做擦除操作。3.2 融合類型Dissolve的三種策略這是緩沖區(qū)操作里第二容易出錯(cuò)的地方。融合Dissolve決定了當(dāng)多個(gè)要素的緩沖區(qū)重疊時(shí)結(jié)果怎么處理。通常有三種選項(xiàng)None不融合每個(gè)要素生成獨(dú)立的緩沖區(qū)多邊形重疊部分保留。結(jié)果圖層里要素?cái)?shù)量等于輸入要素?cái)?shù)量。All全部融合所有緩沖區(qū)合并成一個(gè)大多邊形重疊部分被合并。結(jié)果圖層里只有一個(gè)要素。By field按字段融合按指定字段分組融合同一組內(nèi)的緩沖區(qū)合并不同組之間保持獨(dú)立。選哪種取決于你的分析目的。如果你要統(tǒng)計(jì)每個(gè)設(shè)施的服務(wù)范圍用None因?yàn)槟阈枰烂總€(gè)設(shè)施各自覆蓋了哪里。如果你要統(tǒng)計(jì)總覆蓋面積用All避免重疊區(qū)域被重復(fù)計(jì)算。如果你要按行政區(qū)統(tǒng)計(jì)覆蓋情況用By field按行政區(qū)代碼融合。實(shí)操心得用None生成緩沖區(qū)后如果直接做面積統(tǒng)計(jì)重疊區(qū)域的面積會(huì)被重復(fù)計(jì)算。我見過一個(gè)項(xiàng)目用None生成的緩沖區(qū)算總覆蓋面積結(jié)果比實(shí)際大了將近40%因?yàn)樵O(shè)施密集區(qū)域重疊嚴(yán)重。后來改成All重新算數(shù)據(jù)才對上。3.3 端點(diǎn)樣式與連接樣式線的緩沖區(qū)專屬參數(shù)對線要素做緩沖區(qū)時(shí)有兩個(gè)參數(shù)專門控制線的端點(diǎn)和拐角怎么處理。端點(diǎn)樣式End cap style有三種Round圓頭線的兩端是半圓形這是默認(rèn)值也是最符合“距離不超過X”定義的樣式。Flat平頭線的兩端是平的緩沖區(qū)在端點(diǎn)處直接截?cái)嗖幌蛲鈹U(kuò)展。Square方頭線的兩端是方形相當(dāng)于在端點(diǎn)處向外擴(kuò)展了一個(gè)距離等于緩沖距離的矩形。連接樣式Join style控制線拐角處的處理方式Round圓角拐角處是圓弧默認(rèn)值。Mitre尖角拐角處是尖角可以設(shè)置尖角限制Mitre limit超過限制時(shí)自動(dòng)切成平角。Bevel斜角拐角處是斜切角。這些參數(shù)在道路分析、管線分析里很重要。比如做道路噪聲影響范圍用Round端點(diǎn)更符合聲波傳播的實(shí)際情況。做建筑退線分析用Flat端點(diǎn)可能更合適因?yàn)橥司€通常從建筑邊界算起不向外延伸。3.4 分段數(shù)圓滑度與性能的權(quán)衡點(diǎn)緩沖和圓角連接都是用正多邊形逼近圓分段數(shù)Segments決定用多少條邊來逼近。QGIS默認(rèn)是5段每四分之一圓5段整個(gè)圓20段ArcGIS Pro默認(rèn)是100段整個(gè)圓100段。分段數(shù)越高緩沖區(qū)越圓滑但頂點(diǎn)數(shù)越多文件越大后續(xù)分析越慢。怎么選分段數(shù)我的經(jīng)驗(yàn)是看你的分析精度要求和后續(xù)用途。如果只是做可視化展示20到30段足夠了肉眼幾乎看不出棱角。如果要做面積精確計(jì)算分段數(shù)至少要50以上否則多邊形面積會(huì)比真實(shí)圓面積小。如果要做疊加分析分段數(shù)太高會(huì)導(dǎo)致頂點(diǎn)爆炸疊加運(yùn)算慢得讓人想砸鍵盤。有個(gè)計(jì)算公式可以參考正N邊形逼近圓的面積誤差約為(π/N) * 100%。N20時(shí)誤差約15.7%N50時(shí)約6.3%N100時(shí)約3.1%N200時(shí)約1.6%。如果你能接受1%的面積誤差N取300以上。但實(shí)際項(xiàng)目中我通常用50到100段在精度和性能之間取平衡。4. 實(shí)操全流程從數(shù)據(jù)準(zhǔn)備到結(jié)果驗(yàn)證4.1 數(shù)據(jù)準(zhǔn)備與坐標(biāo)系檢查拿到數(shù)據(jù)后第一步永遠(yuǎn)是檢查坐標(biāo)系。在QGIS里右鍵圖層選Properties看Information里的CRS。在ArcGIS Pro里看圖層屬性里的Spatial Reference。如果顯示的是GCS開頭的比如GCS_WGS_1984說明是地理坐標(biāo)系需要投影。投影操作在QGIS里用“Reproject layer”工具在ArcGIS Pro里用“Project”工具。選目標(biāo)坐標(biāo)系時(shí)如果研究區(qū)域在國內(nèi)我通常用CGCS2000的高斯-克呂格投影。具體選哪個(gè)帶看研究區(qū)域的經(jīng)度范圍。3度帶的話中央經(jīng)線從75度開始每3度一帶帶號從25開始。比如北京在東經(jīng)116度左右中央經(jīng)線114度是第38帶EPSG代碼是4547CGCS2000 3度帶中央經(jīng)線114度。投影完之后用“Add geometry attributes”工具算一下要素的坐標(biāo)范圍確認(rèn)投影后的坐標(biāo)值在合理范圍內(nèi)。比如北京的高斯投影坐標(biāo)X應(yīng)該在500000左右?guī)?nèi)偏移Y應(yīng)該在4400000左右。如果X是116、Y是40說明投影沒成功還是經(jīng)緯度。4.2 緩沖區(qū)生成的完整操作以QGIS為例完整流程如下打開Processing Toolbox搜索“Buffer”。選擇輸入圖層。Distance輸入緩沖距離。如果是字段驅(qū)動(dòng)點(diǎn)右邊的按鈕選字段。Segments輸入分段數(shù)我通常填50。End cap style選Round。Join style選Round。Mitre limit保持默認(rèn)2。Dissolve選擇融合策略。勾選“Keep disjoint features separate”如果你想讓不相連的要素保持獨(dú)立。點(diǎn)Run。ArcGIS Pro的流程類似在Geoprocessing里搜“Buffer”參數(shù)名稱略有不同但邏輯一致。注意ArcGIS Pro的Buffer工具默認(rèn)會(huì)做融合如果你不想融合要把Dissolve Type設(shè)為None。Python實(shí)現(xiàn)的話用geopandas大概是這樣import geopandas as gpd gdf gpd.read_file(input.shp) gdf_proj gdf.to_crs(EPSG:4547) gdf_proj[geometry] gdf_proj.geometry.buffer(500, resolution50) gdf_proj.to_file(buffer_output.shp)buffer()方法的resolution參數(shù)對應(yīng)分段數(shù)默認(rèn)是16。注意geopandas的buffer是在要素級別操作的不會(huì)自動(dòng)融合。要融合的話用unary_union或者dissolve。4.3 結(jié)果驗(yàn)證三個(gè)必做的檢查緩沖區(qū)生成后別急著往下做分析先做三個(gè)檢查第一目視檢查。把緩沖區(qū)疊加到原始要素上看看形狀對不對有沒有明顯的變形或缺失。特別注意線要素的端點(diǎn)和拐角以及面要素的內(nèi)外環(huán)。第二距離抽查。隨機(jī)選幾個(gè)點(diǎn)量一下從原始要素邊界到緩沖區(qū)邊界的距離看是否等于設(shè)定的緩沖距離。在QGIS里用Measure工具在ArcGIS Pro里用Measure工具。如果偏差超過1%說明坐標(biāo)系或分段數(shù)有問題。第三面積核對。對點(diǎn)要素緩沖區(qū)面積應(yīng)該接近πr2。比如500米緩沖理論面積是785398平方米。如果算出來是700000說明分段數(shù)太低。對線要素緩沖區(qū)面積約等于線長度乘以2r加上πr2兩端半圓。對面要素緩沖區(qū)面積等于原面積加上周長乘以r加上πr2外角。這些公式可以用來快速估算判斷結(jié)果是否合理。踩過的坑有一次做管線緩沖區(qū)用了一個(gè)自定義的投影坐標(biāo)系結(jié)果緩沖區(qū)生成后距離抽查發(fā)現(xiàn)偏了15%。查了半天才發(fā)現(xiàn)那個(gè)投影坐標(biāo)系的單位是英尺不是米。輸入500以為是500米實(shí)際是500英尺約152米。從那以后我養(yǎng)成了習(xí)慣投影后一定用測量工具量一下已知距離確認(rèn)單位。5. 進(jìn)階話題緩沖區(qū)之外的距離分析5.1 柵格歐氏距離另一種思路矢量緩沖區(qū)生成的是硬邊界多邊形柵格歐氏距離生成的是連續(xù)的距離表面。在ArcGIS Pro里用“Euclidean Distance”工具在QGIS里用“Proximity”下的“Raster distance”工具。輸入是柵格化的要素輸出是每個(gè)像元到最近要素的距離值。柵格距離的優(yōu)勢在于可以做連續(xù)的距離分級比如把距離分成0-100、100-200、200-500、500-1000幾個(gè)帶每個(gè)帶用不同顏色顯示。也可以直接做閾值提取得到和矢量緩沖區(qū)類似的結(jié)果。劣勢在于精度受像元大小限制而且柵格數(shù)據(jù)量大。選擇矢量還是柵格看你的分析需求。如果后續(xù)要做疊加分析、網(wǎng)絡(luò)分析矢量緩沖區(qū)更方便。如果要做表面分析、成本距離分析柵格距離更合適。5.2 網(wǎng)絡(luò)緩沖區(qū)沿路網(wǎng)的距離歐氏緩沖區(qū)假設(shè)你可以從任何方向直線到達(dá)目標(biāo)但現(xiàn)實(shí)中你只能沿著路網(wǎng)走。網(wǎng)絡(luò)緩沖區(qū)Network Buffer就是解決這個(gè)問題的。它基于路網(wǎng)數(shù)據(jù)計(jì)算從設(shè)施點(diǎn)出發(fā)沿路網(wǎng)行進(jìn)指定距離能覆蓋的范圍。ArcGIS Pro里有“Network Analyst”擴(kuò)展模塊可以做這個(gè)QGIS里可以用“QNEAT3”插件。網(wǎng)絡(luò)緩沖區(qū)的結(jié)果不是平滑的多邊形而是沿路網(wǎng)分布的帶狀區(qū)域。做應(yīng)急服務(wù)覆蓋分析、公交站點(diǎn)服務(wù)范圍分析時(shí)網(wǎng)絡(luò)緩沖區(qū)比歐氏緩沖區(qū)準(zhǔn)確得多。我做過一個(gè)急救站覆蓋分析的項(xiàng)目用歐氏緩沖區(qū)算出來覆蓋率95%用網(wǎng)絡(luò)緩沖區(qū)算出來只有78%。差距在哪歐氏緩沖區(qū)把河流、山體、封閉小區(qū)都算進(jìn)去了但實(shí)際上救護(hù)車得繞路走。這個(gè)案例讓我深刻體會(huì)到緩沖區(qū)分析的結(jié)果是否可信取決于你對現(xiàn)實(shí)約束的建模是否到位。5.3 性能優(yōu)化大數(shù)據(jù)量下的緩沖區(qū)處理當(dāng)要素?cái)?shù)量上萬、緩沖距離又大時(shí)緩沖區(qū)操作可能變得很慢。幾個(gè)優(yōu)化思路第一先做空間索引。QGIS和ArcGIS Pro都會(huì)自動(dòng)建空間索引但如果你用Python腳本處理記得先建索引。geopandas的sindex屬性可以查看索引。第二分塊處理。把大范圍數(shù)據(jù)切成小塊逐塊緩沖最后合并。這樣可以避免內(nèi)存溢出也能利用多核并行。第三降低分段數(shù)。如果只是中間結(jié)果分段數(shù)用20就夠了最后出圖時(shí)再重新生成高精度的。第四用PostGIS。如果數(shù)據(jù)在PostgreSQL里用PostGIS的ST_Buffer函數(shù)性能比桌面軟件好很多而且支持并行查詢。語法是ST_Buffer(geometry, distance, quad_segs50)。第五考慮用GeoPandas的buffer配合Dask。對于超大數(shù)據(jù)集Dask-GeoPandas可以把數(shù)據(jù)分片并行處理速度提升明顯。6. 常見問題與排查速查表6.1 緩沖區(qū)結(jié)果偏移或變形這是最常見的問題原因通常有三個(gè)坐標(biāo)系不對、投影帶選錯(cuò)、數(shù)據(jù)本身有拓?fù)溴e(cuò)誤。排查步驟先確認(rèn)坐標(biāo)系再檢查投影參數(shù)最后用拓?fù)錂z查工具查一下原始數(shù)據(jù)有沒有自相交、懸掛點(diǎn)之類的問題。6.2 緩沖區(qū)面積與理論值差距大分段數(shù)太低是主因。把Segments從默認(rèn)值提高到50以上重新生成。如果還不對檢查是不是用了地理坐標(biāo)系或者投影變形太大。6.3 融合后要素?cái)?shù)量不對Dissolve參數(shù)選錯(cuò)了。想要每個(gè)要素獨(dú)立選None想要全部合并選All想要按字段分組選By field并指定字段。注意有些軟件默認(rèn)是All有些默認(rèn)是None操作前看清楚。6.4 線緩沖區(qū)端點(diǎn)形狀不對End cap style選錯(cuò)了。Round是圓頭Flat是平頭Square是方頭。根據(jù)分析目的選擇不確定就用Round。6.5 緩沖區(qū)在極地或跨越日期變更線時(shí)出錯(cuò)這是測地線計(jì)算的邊界情況。解決方案是換用局部投影或者把數(shù)據(jù)切分成不跨越這些特殊區(qū)域的部分分別處理。6.6 緩沖區(qū)結(jié)果有自相交或空洞原始數(shù)據(jù)有拓?fù)溴e(cuò)誤或者緩沖距離太大導(dǎo)致幾何自相交。先用“Fix geometries”工具修復(fù)原始數(shù)據(jù)再重新緩沖。如果緩沖距離超過地球周長的一半任何平面緩沖都會(huì)出問題這時(shí)候必須用測地線方法。問題現(xiàn)象可能原因排查方法解決方案緩沖區(qū)整體偏移坐標(biāo)系錯(cuò)誤檢查CRS和單位投影到正確的投影坐標(biāo)系面積偏小分段數(shù)太低對比理論面積提高Segments到50以上要素?cái)?shù)量不對Dissolve參數(shù)錯(cuò)誤檢查融合設(shè)置按需求選None/All/By field端點(diǎn)形狀異常End cap style錯(cuò)誤查看線端點(diǎn)改為Round幾何自相交原始數(shù)據(jù)拓?fù)溴e(cuò)誤運(yùn)行拓?fù)錂z查修復(fù)幾何后重新緩沖計(jì)算速度極慢頂點(diǎn)數(shù)過多或數(shù)據(jù)量大查看要素頂點(diǎn)數(shù)降低分段數(shù)、分塊處理、用PostGIS7. 幾個(gè)讓我印象深刻的實(shí)戰(zhàn)案例7.1 商業(yè)選址多環(huán)緩沖區(qū)的疊加分析前年做過一個(gè)連鎖便利店選址的項(xiàng)目。需求是在已有門店周圍找空白區(qū)域同時(shí)要考慮競爭對手的覆蓋范圍。我的做法是對自家門店做500米多環(huán)緩沖300、500、800米對競爭對手門店做500米單環(huán)緩沖然后把自家緩沖區(qū)減去競爭對手緩沖區(qū)得到“有效覆蓋空白區(qū)”。再疊加人口密度柵格找出人口密度高但覆蓋不足的區(qū)域。這個(gè)項(xiàng)目里緩沖區(qū)的融合策略很關(guān)鍵。自家門店的緩沖區(qū)用By field按門店編號融合這樣能知道每個(gè)門店的覆蓋范圍。競爭對手的緩沖區(qū)用All融合因?yàn)椴恍枰獏^(qū)分是哪家對手只需要知道“被覆蓋”這個(gè)事實(shí)。最后做擦除操作時(shí)用自家緩沖區(qū)圖層擦除競爭對手緩沖區(qū)圖層得到的結(jié)果就是每個(gè)自家門店的獨(dú)占覆蓋區(qū)。7.2 環(huán)境評估線緩沖區(qū)的端點(diǎn)處理另一個(gè)項(xiàng)目是做高速公路噪聲影響評估。需要對高速公路做200米緩沖區(qū)表示噪聲影響范圍。這里端點(diǎn)樣式就很重要了。高速公路的起點(diǎn)和終點(diǎn)如果直接截?cái)郌lat會(huì)低估端點(diǎn)附近的噪聲影響。用Round更合理因?yàn)樵肼晱亩它c(diǎn)也會(huì)向外傳播。但Round會(huì)在端點(diǎn)處生成半圓如果高速公路端點(diǎn)連接著其他道路這個(gè)半圓可能會(huì)覆蓋到不該覆蓋的區(qū)域。最后我的處理是主體用Round但在端點(diǎn)處手動(dòng)裁剪去掉不合理的部分。這個(gè)案例說明緩沖區(qū)參數(shù)沒有絕對的對錯(cuò)只有是否適合具體場景。工具是死的人是活的該手動(dòng)調(diào)整的時(shí)候就得手動(dòng)調(diào)整。7.3 應(yīng)急避難網(wǎng)絡(luò)緩沖區(qū)的必要性前面提到的急救站覆蓋分析補(bǔ)充一點(diǎn)細(xì)節(jié)。歐氏緩沖區(qū)和網(wǎng)絡(luò)緩沖區(qū)的差距主要來自河流、鐵路、封閉式小區(qū)這些“不可穿越”的障礙。做這類分析時(shí)如果只用歐氏緩沖區(qū)結(jié)果會(huì)過于樂觀。我的做法是先用歐氏緩沖區(qū)做初步篩選找出可能覆蓋不足的區(qū)域再用網(wǎng)絡(luò)緩沖區(qū)做精確分析。這樣既保證了效率又保證了精度。網(wǎng)絡(luò)緩沖區(qū)的實(shí)現(xiàn)我用的是QNEAT3插件的“Iso-Area”功能。輸入路網(wǎng)圖層、設(shè)施點(diǎn)圖層、時(shí)間或距離閾值輸出等時(shí)圈或等距圈。注意路網(wǎng)數(shù)據(jù)必須有正確的拓?fù)浣徊婵谝驍喾駝t路徑計(jì)算會(huì)出錯(cuò)。8. 寫在最后一些零散但有用的經(jīng)驗(yàn)緩沖區(qū)分析這個(gè)操作工具本身很簡單難的是對空間數(shù)據(jù)的理解和對自己分析目的的清晰認(rèn)知。我見過太多人把緩沖區(qū)當(dāng)成一個(gè)“點(diǎn)一下就行”的操作結(jié)果在后續(xù)分析里被各種問題折磨。如果你只能從這篇文章里記住一件事我希望是做緩沖區(qū)之前先問自己三個(gè)問題——數(shù)據(jù)在什么坐標(biāo)系緩沖距離的物理意義是什么融合策略是否符合分析目的這三個(gè)問題想清楚了八成的問題都不會(huì)發(fā)生。另外分享一個(gè)小技巧如果你經(jīng)常做緩沖區(qū)分析可以寫一個(gè)Python腳本或者QGIS模型把常用的參數(shù)組合固化下來。比如“投影到CGCS2000 3度帶 緩沖500米 分段50 全部融合”這一套流程做成一個(gè)一鍵操作。這樣既能保證每次操作的一致性又能節(jié)省大量時(shí)間。我自己維護(hù)了一套這樣的腳本覆蓋了選址分析、環(huán)境評估、服務(wù)區(qū)劃分等幾個(gè)常見場景用起來很順手。最后說一個(gè)容易被忽略的點(diǎn)緩沖區(qū)結(jié)果的屬性表。默認(rèn)情況下緩沖區(qū)結(jié)果會(huì)繼承原始要素的所有屬性。如果原始數(shù)據(jù)屬性很多緩沖區(qū)結(jié)果的屬性表會(huì)很大影響后續(xù)處理速度?？梢栽谏删彌_區(qū)時(shí)只保留必要的字段或者生成后用“Drop fields”工具清理。這個(gè)細(xì)節(jié)在數(shù)據(jù)量大的時(shí)候能省不少事。