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