建正方形網(wǎng)格：從坐標(biāo)系選擇到自動(dòng)化腳本全解析)
1. 項(xiàng)目緣起為什么我們需要在QGIS中創(chuàng)建正方形網(wǎng)格在地理信息處理的實(shí)際工作中我們常常會(huì)遇到一個(gè)看似簡(jiǎn)單卻至關(guān)重要的需求如何將一個(gè)不規(guī)則的研究區(qū)域高效、精確地分割成一個(gè)個(gè)規(guī)整的單元無(wú)論是進(jìn)行城市地塊的抽樣調(diào)查、生態(tài)學(xué)中的樣方分析、還是區(qū)域性的資源評(píng)估將連續(xù)的空間離散化為規(guī)則的網(wǎng)格都是后續(xù)統(tǒng)計(jì)分析、模型構(gòu)建和可視化表達(dá)的基礎(chǔ)。手動(dòng)繪制這些網(wǎng)格不僅效率低下而且難以保證幾何精度和拓?fù)湟恢滦?。QGIS作為一款功能強(qiáng)大的開(kāi)源GIS桌面軟件其“創(chuàng)建網(wǎng)格”工具正是為此而生。它允許用戶基于指定的范圍、形狀如矩形、菱形、六邊形和尺寸批量生成覆蓋目標(biāo)區(qū)域的規(guī)則幾何圖形。其中正方形網(wǎng)格因其計(jì)算簡(jiǎn)便、易于理解和空間索引高效成為最常用的網(wǎng)格類(lèi)型。然而從“知道有這個(gè)工具”到“能穩(wěn)定產(chǎn)出符合預(yù)期、坐標(biāo)正確的網(wǎng)格”中間往往隔著對(duì)坐標(biāo)系、參數(shù)設(shè)置和輸出結(jié)果處理的深刻理解。很多新手甚至是有一定經(jīng)驗(yàn)的用戶都可能在“創(chuàng)建網(wǎng)格”這一步就踩坑比如生成的網(wǎng)格大小不對(duì)、覆蓋范圍有偏差或者后續(xù)疊加數(shù)據(jù)時(shí)出現(xiàn)錯(cuò)位。本文將從一個(gè)資深GIS從業(yè)者的視角手把手帶你深入QGIS的“創(chuàng)建網(wǎng)格”功能。我們不僅會(huì)完成“點(diǎn)擊按鈕生成網(wǎng)格”這個(gè)動(dòng)作更會(huì)拆解其背后的每一個(gè)參數(shù)邏輯特別是坐標(biāo)系如EPSG:3857的選擇與影響并分享在實(shí)際項(xiàng)目中如何避免常見(jiàn)陷阱確保你得到的每一個(gè)正方形網(wǎng)格都精準(zhǔn)、可用。2. 核心工具解析“創(chuàng)建網(wǎng)格”的功能界面與參數(shù)全解“創(chuàng)建網(wǎng)格”工具位于QGIS的“矢量”菜單下的“研究工具”子菜單中。打開(kāi)后你會(huì)看到一個(gè)參數(shù)對(duì)話框。這個(gè)對(duì)話框里的每一個(gè)選項(xiàng)都直接決定了最終輸出網(wǎng)格的形態(tài)和空間屬性。我們來(lái)逐一拆解并解釋其背后的“為什么”。2.1 網(wǎng)格類(lèi)型與范圍定義首先在“網(wǎng)格類(lèi)型”下拉菜單中選擇“矩形多邊形”。這里有一個(gè)關(guān)鍵點(diǎn)雖然我們最終想要的是正方形但QGIS是通過(guò)指定“寬度”和“高度”來(lái)定義矩形的。當(dāng)我們將寬度和高度設(shè)置為相同的數(shù)值時(shí)得到的就是正方形。選擇“多邊形”而非“線”是因?yàn)槎噙呅尉W(wǎng)格才能作為一個(gè)獨(dú)立的面要素進(jìn)行后續(xù)的面積計(jì)算、空間連接等操作。接下來(lái)是“網(wǎng)格范圍”。這是第一個(gè)容易出錯(cuò)的環(huán)節(jié)。你有四種方式定義范圍畫(huà)布范圍直接使用當(dāng)前地圖視圖顯示的范圍。這非常方便但風(fēng)險(xiǎn)在于如果你縮放或平移了地圖范圍就變了可重復(fù)性差。圖層范圍選擇某個(gè)已有矢量圖層如一個(gè)表示研究區(qū)的面圖層工具會(huì)計(jì)算該圖層的整體外包矩形即能完全包含該圖層所有要素的最小矩形。這是最常用且推薦的方式因?yàn)樗_保了網(wǎng)格能完全覆蓋你的目標(biāo)區(qū)域。書(shū)簽范圍如果你之前保存了地圖書(shū)簽可以調(diào)用。手動(dòng)輸入坐標(biāo)直接輸入左下角和右上角的X、Y坐標(biāo)。這種方式最精確但需要你事先知道精確的范圍坐標(biāo)。注意無(wú)論選擇哪種方式這里定義的范圍都是一個(gè)矩形。如果你的研究區(qū)形狀極不規(guī)則比如一個(gè)彎曲的河流流域使用其外包矩形創(chuàng)建網(wǎng)格會(huì)在研究區(qū)外圍產(chǎn)生大量“無(wú)效”的網(wǎng)格單元。后續(xù)通常需要用一個(gè)“按位置選擇”或“裁剪”工具將落在實(shí)際研究區(qū)之外的網(wǎng)格剔除。2.2 網(wǎng)格間距理解“寬度”與“高度”的單位這是核心參數(shù)。在“水平間距”和“垂直間距”中分別輸入你想要的網(wǎng)格單元的寬度和高度。例如輸入“100”期望得到邊長(zhǎng)為100單位的正方形。這里的坑在于“單位”。這個(gè)“100”單位是什么米公里度這完全取決于你為輸出網(wǎng)格所選擇的坐標(biāo)系CRS。如果輸出CRS是地理坐標(biāo)系如EPSG:4326, WGS84單位是十進(jìn)制度。在地球表面1度經(jīng)度的實(shí)際長(zhǎng)度隨緯度變化在赤道最長(zhǎng)向兩極遞減至01度緯度的長(zhǎng)度相對(duì)恒定約111公里。此時(shí)設(shè)置“100”度將得到一個(gè)在現(xiàn)實(shí)世界中巨大且扭曲的“正方形”這通常不是我們想要的。如果輸出CRS是投影坐標(biāo)系如EPSG:3857 Web Mercator, 或EPSG:32649 UTM Zone 49N單位通常是米。此時(shí)設(shè)置“100”就意味著邊長(zhǎng)100米的正方形。這才是進(jìn)行面積計(jì)算和空間分析的理想選擇。因此在輸入間距值之前你必須明確你希望網(wǎng)格在現(xiàn)實(shí)世界中的實(shí)際尺寸例如100米 x 100米并確保為輸出圖層選擇了合適的、以米為單位的投影坐標(biāo)系。工具界面中的“輸出CRS”選項(xiàng)就是用來(lái)定義這一點(diǎn)的。2.3 輸出CRS坐標(biāo)系選擇的深遠(yuǎn)影響“輸出CRS”是決定網(wǎng)格空間正確性的靈魂參數(shù)。它有兩個(gè)主要選項(xiàng)項(xiàng)目CRS使用當(dāng)前QGIS項(xiàng)目設(shè)置的坐標(biāo)系。指定CRS手動(dòng)選擇一個(gè)坐標(biāo)系。強(qiáng)烈建議選擇“指定CRS”并選擇一個(gè)適合你研究區(qū)域的投影坐標(biāo)系。原因如下保證尺寸準(zhǔn)確如上所述投影坐標(biāo)系使用線性單位米、英尺能保證網(wǎng)格邊長(zhǎng)在地圖上的測(cè)量值與實(shí)際地面距離對(duì)應(yīng)。避免幾何變形正方形網(wǎng)格在投影坐標(biāo)系下才能保持較好的形狀。在地理坐標(biāo)系下由于將球面展開(kāi)為平面帶來(lái)的變形網(wǎng)格會(huì)變成梯形或其他形狀。與底圖/業(yè)務(wù)數(shù)據(jù)對(duì)齊如果你的在線底圖如OSM、谷歌地圖或業(yè)務(wù)數(shù)據(jù)使用的是Web MercatorEPSG:3857那么將網(wǎng)格也創(chuàng)建為相同的坐標(biāo)系可以確保完美的疊加沒(méi)有偏移。關(guān)于EPSG:3857這是一個(gè)全球性的投影坐標(biāo)系被幾乎所有網(wǎng)絡(luò)地圖服務(wù)使用。它的單位是“米”但在高緯度地區(qū)面積和形狀變形非常嚴(yán)重。如果你的研究區(qū)域范圍不大比如一個(gè)城市且主要為了與在線地圖疊加顯示使用3857是方便的選擇。但如果需要進(jìn)行精確的面積計(jì)算或大范圍的空間分析應(yīng)優(yōu)先使用更優(yōu)的局部投影坐標(biāo)系如UTM通用橫軸墨卡托分區(qū)。你可以通過(guò)QGIS的“按圖層范圍自動(dòng)選擇”功能或查詢當(dāng)?shù)氐臏y(cè)繪標(biāo)準(zhǔn)來(lái)確定最適合的UTM帶號(hào)如EPSG:32650 for UTM Zone 50N。2.4 其他參數(shù)與輸出水平/垂直格網(wǎng)數(shù)這是一個(gè)與“間距”互斥的參數(shù)。如果你知道范圍的總寬度是1000米并想分成10列那么可以設(shè)置“水平格網(wǎng)數(shù)”為10工具會(huì)自動(dòng)計(jì)算間距為100米。我個(gè)人更傾向于直接使用“間距”因?yàn)槲覍?duì)實(shí)地尺寸更敏感。網(wǎng)格對(duì)齊點(diǎn)可以設(shè)置網(wǎng)格的起始偏移點(diǎn)。通常不需要調(diào)整除非你有特殊的對(duì)齊要求如讓網(wǎng)格的角點(diǎn)與某個(gè)已知坐標(biāo)點(diǎn)對(duì)齊。輸出指定生成網(wǎng)格的保存路徑和格式如GeoPackage、Shapefile。點(diǎn)擊“運(yùn)行”一個(gè)全新的、覆蓋指定范圍的規(guī)則網(wǎng)格圖層就會(huì)加載到你的地圖中。3. 實(shí)戰(zhàn)演練從零生成一個(gè)用于抽樣分析的正方形網(wǎng)格讓我們通過(guò)一個(gè)完整的模擬案例將上述理論付諸實(shí)踐。假設(shè)我們有一個(gè)“某市森林公園邊界.shp”的面數(shù)據(jù)我們需要在其范圍內(nèi)創(chuàng)建邊長(zhǎng)為500米的正方形網(wǎng)格用于規(guī)劃生態(tài)監(jiān)測(cè)樣地。3.1 數(shù)據(jù)準(zhǔn)備與坐標(biāo)系檢查首先將“森林公園邊界.shp”加載到QGIS中。右鍵點(diǎn)擊圖層選擇“屬性”切換到“信息”選項(xiàng)卡。查看“坐標(biāo)系”信息。假設(shè)它顯示為“EPSG:4547”這是一個(gè)適用于該地區(qū)的投影坐標(biāo)系單位是米。這很好我們的數(shù)據(jù)本身就在一個(gè)合適的投影坐標(biāo)系中。如果數(shù)據(jù)是地理坐標(biāo)系如EPSG:4326你需要先考慮是否要將其重新投影到一個(gè)投影坐標(biāo)系。你可以在導(dǎo)入時(shí)使用“圖層CRS”進(jìn)行設(shè)置或者后續(xù)使用“導(dǎo)出”或“重投影圖層”工具來(lái)處理?；驹瓌t是用于創(chuàng)建網(wǎng)格的“范圍圖層”和“輸出網(wǎng)格”最好在同一個(gè)、以米為單位的投影坐標(biāo)系下以避免不必要的轉(zhuǎn)換誤差。3.2 執(zhí)行“創(chuàng)建網(wǎng)格”工具打開(kāi)“矢量” - “研究工具” - “創(chuàng)建網(wǎng)格”。網(wǎng)格類(lèi)型選擇“矩形多邊形”。網(wǎng)格范圍選擇“圖層范圍”并在下拉菜單中選擇“森林公園邊界”。水平間距輸入500。垂直間距輸入500。輸出CRS選擇“指定CRS”。點(diǎn)擊右側(cè)的“選擇CRS”按鈕地球圖標(biāo)。在過(guò)濾器中輸入“4547”選擇它。關(guān)鍵點(diǎn)這里我們選擇與源數(shù)據(jù)相同的CRS確保網(wǎng)格與原始邊界完美套合。輸出點(diǎn)擊“...”按鈕將結(jié)果保存為GeoPackage格式命名為“森林公園_500米網(wǎng)格.gpkg”。點(diǎn)擊“運(yùn)行”。瞬間一個(gè)整齊的網(wǎng)格圖層就覆蓋在了森林公園之上。3.3 結(jié)果驗(yàn)證與后處理生成網(wǎng)格后不能直接使用必須進(jìn)行驗(yàn)證和清理。驗(yàn)證1尺寸檢查。使用“測(cè)量工具”設(shè)置為“米”為單位隨機(jī)測(cè)量幾個(gè)網(wǎng)格的邊長(zhǎng)。確認(rèn)它們確實(shí)是500米左右由于范圍除間距可能不是整數(shù)邊緣網(wǎng)格可能略小。驗(yàn)證2范圍檢查。觀察網(wǎng)格是否完全覆蓋了森林公園。由于我們用的是圖層的“外包矩形”網(wǎng)格范圍肯定會(huì)比實(shí)際公園邊界大四周會(huì)多出很多覆蓋了公園外部區(qū)域的網(wǎng)格。后處理裁剪網(wǎng)格。我們需要只保留落在森林公園內(nèi)部的網(wǎng)格。使用“矢量” - “地理處理工具” - “裁剪”工具?！拜斎雸D層”選擇我們剛生成的“森林公園_500米網(wǎng)格”?！隘B加圖層”選擇“森林公園邊界”。運(yùn)行。得到的新圖層“裁剪后的網(wǎng)格”就只包含完全位于公園內(nèi)部的完整或部分網(wǎng)格了。后處理添加唯一標(biāo)識(shí)。為了方便管理為每個(gè)網(wǎng)格賦予一個(gè)ID。打開(kāi)“裁剪后的網(wǎng)格”圖層的屬性表。點(diǎn)擊“打開(kāi)字段計(jì)算器”創(chuàng)建一個(gè)新字段例如叫g(shù)rid_id。在表達(dá)式框中輸入row_number這將為每一行即每個(gè)網(wǎng)格生成一個(gè)從1開(kāi)始的序列號(hào)。點(diǎn)擊“確定”。現(xiàn)在每個(gè)網(wǎng)格都有了自己的編號(hào)。至此一套可用于布設(shè)監(jiān)測(cè)樣地、進(jìn)行分區(qū)統(tǒng)計(jì)的標(biāo)準(zhǔn)化正方形網(wǎng)格就制作完成了。你可以基于grid_id來(lái)隨機(jī)抽樣或者將其他數(shù)據(jù)如植被類(lèi)型、動(dòng)物觀測(cè)點(diǎn)與這些網(wǎng)格進(jìn)行空間連接計(jì)算每個(gè)網(wǎng)格內(nèi)的統(tǒng)計(jì)值。4. 深度避坑坐標(biāo)系3857 vs. 局部投影帶來(lái)的隱秘陷阱“創(chuàng)建網(wǎng)格”操作中最頑固、最不易察覺(jué)的問(wèn)題幾乎都與坐標(biāo)系有關(guān)。我們重點(diǎn)剖析EPSG:3857這個(gè)“熟悉的陌生人”。4.1 EPSG:3857的便利性與局限性便利性EPSG:3857Web Mercator是網(wǎng)絡(luò)地圖的“世界語(yǔ)”。幾乎所有的在線瓦片地圖Google Maps, OpenStreetMap, 天地圖等都使用它。在QGIS中加載XYZ Tiles連接這些在線地圖時(shí)項(xiàng)目CRS會(huì)自動(dòng)切換到3857。如果你在此狀態(tài)下直接用畫(huà)布范圍創(chuàng)建網(wǎng)格并選擇“項(xiàng)目CRS”作為輸出那么網(wǎng)格會(huì)與底圖完美貼合視覺(jué)上非常整齊。局限性即大坑嚴(yán)重的面積變形Web Mercator投影在高緯度地區(qū)會(huì)極度夸大面積。例如在赤道附近一個(gè)3857坐標(biāo)系下邊長(zhǎng)為100米的網(wǎng)格其實(shí)際地面面積接近1公頃。但在北緯60度地區(qū)同一個(gè)“100米”網(wǎng)格所代表的實(shí)際地面面積會(huì)小得多。這意味著如果你在3857下創(chuàng)建了“500米網(wǎng)格”并用于面積相關(guān)的分析如計(jì)算網(wǎng)格內(nèi)的森林面積占比結(jié)果將是完全錯(cuò)誤的。非等角投影雖然它被宣傳為“保角”但實(shí)際上在非標(biāo)準(zhǔn)緯線上并不完全保角網(wǎng)格的形狀也會(huì)發(fā)生微小變形并非嚴(yán)格的正方形。4.2 一個(gè)真實(shí)的對(duì)比實(shí)驗(yàn)假設(shè)我們?cè)诒本┑貐^(qū)約北緯40度創(chuàng)建一個(gè)邊長(zhǎng)為1000米的網(wǎng)格。方案A錯(cuò)誤示范項(xiàng)目CRS為EPSG:3857加載OSM底圖以當(dāng)前視圖范圍創(chuàng)建網(wǎng)格輸出CRS使用“項(xiàng)目CRS”即3857。方案B正確做法將項(xiàng)目CRS設(shè)置為更適合該區(qū)域的投影如EPSG:32650UTM Zone 50N。重新加載數(shù)據(jù)或?qū)?shù)據(jù)重投影然后以圖層范圍創(chuàng)建網(wǎng)格輸出CRS指定為EPSG:32650。對(duì)比結(jié)果視覺(jué)對(duì)齊方案A的網(wǎng)格與OSM底圖嚴(yán)絲合縫。方案B的網(wǎng)格在OSM底圖3857上會(huì)有可見(jiàn)的旋轉(zhuǎn)和偏移因?yàn)榈讏D被動(dòng)態(tài)投影到了3857顯示而網(wǎng)格是32650的。這常常讓新手困惑認(rèn)為B方案“錯(cuò)了”。實(shí)際上視覺(jué)上的錯(cuò)位是投影顯示造成的B方案的數(shù)據(jù)在它自己的坐標(biāo)系下才是正確的。面積真相使用QGIS的“字段計(jì)算器”和$area函數(shù)計(jì)算兩個(gè)方案下網(wǎng)格的面積。方案A3857網(wǎng)格$area返回的值接近 1,000,000 平方米即1平方公里但這只是地圖上的“投影面積”。方案B32650網(wǎng)格$area返回的值也接近 1,000,000 平方米。關(guān)鍵來(lái)了使用“測(cè)量工具”設(shè)置為“橢球體”模式如WGS84分別測(cè)量?jī)蓚€(gè)網(wǎng)格的實(shí)際地面距離。你會(huì)發(fā)現(xiàn)方案B的網(wǎng)格邊長(zhǎng)更接近真實(shí)的1000米而方案A的網(wǎng)格實(shí)際地面邊長(zhǎng)要小于1000米。其“投影面積”1平方公里是虛高的。實(shí)操心得永遠(yuǎn)不要為了與在線底圖視覺(jué)對(duì)齊的便利而犧牲空間數(shù)據(jù)的幾何精度。正確的 workflow 是1) 為你的分析確定一個(gè)合適的局部投影坐標(biāo)系CRS。2) 將所有業(yè)務(wù)數(shù)據(jù)包括用于定義范圍的圖層都轉(zhuǎn)換或創(chuàng)建到這個(gè)CRS下。3) 在此CRS下創(chuàng)建網(wǎng)格、進(jìn)行分析計(jì)算。4) 如果需要在線底圖作為參考可以添加X(jué)YZ TilesQGIS會(huì)動(dòng)態(tài)將其投影到你的項(xiàng)目CRS中顯示此時(shí)底圖可能會(huì)有變形和偏移但你的業(yè)務(wù)數(shù)據(jù)之間的空間關(guān)系是精確的。記住底圖是“參考背景”你的分析數(shù)據(jù)才是“主角”。4.3 處理“創(chuàng)建網(wǎng)格”時(shí)常見(jiàn)的坐標(biāo)系報(bào)錯(cuò)有時(shí)在運(yùn)行工具時(shí)會(huì)遇到“圖層范圍在輸出CRS中無(wú)效”或類(lèi)似錯(cuò)誤。這通常是因?yàn)椤胺秶鷪D層”的坐標(biāo)系與“輸出CRS”不兼容或者范圍坐標(biāo)值超出了目標(biāo)CRS的有效域例如試圖在UTM投影中創(chuàng)建一個(gè)經(jīng)度跨度過(guò)大的網(wǎng)格。解決方法統(tǒng)一坐標(biāo)系確保范圍圖層和輸出CRS一致或至少是能夠相互轉(zhuǎn)換的。檢查范圍值在范圍圖層的屬性中查看其實(shí)際坐標(biāo)范圍判斷是否合理。使用項(xiàng)目CRS一個(gè)穩(wěn)妥的方法是先將項(xiàng)目和所有相關(guān)圖層的CRS都設(shè)置為你的目標(biāo)分析CRS然后使用“畫(huà)布范圍”或“圖層范圍”來(lái)創(chuàng)建網(wǎng)格。5. 進(jìn)階應(yīng)用與腳本化讓網(wǎng)格創(chuàng)建融入自動(dòng)化工作流對(duì)于需要頻繁創(chuàng)建不同尺寸網(wǎng)格或作為復(fù)雜模型一部分的用戶手動(dòng)點(diǎn)擊GUI效率太低。QGIS的強(qiáng)大之處在于其處理框架和PyQGIS允許我們將這個(gè)過(guò)程腳本化。5.1 使用處理模型設(shè)計(jì)器構(gòu)建可復(fù)用流程我們可以將“創(chuàng)建網(wǎng)格”和后續(xù)的“裁剪”、“添加ID”等步驟打包成一個(gè)自定義模型。打開(kāi)“處理” - “圖形模型設(shè)計(jì)器”。從右側(cè)算法列表中將“創(chuàng)建網(wǎng)格”、“裁剪”、“字段計(jì)算器”拖入畫(huà)布。進(jìn)行連接將“范圍圖層”作為“創(chuàng)建網(wǎng)格”和“裁剪”的輸入將“創(chuàng)建網(wǎng)格”的輸出連接到“裁剪”的輸入圖層將“裁剪”的輸出連接到“字段計(jì)算器”的輸入。為每個(gè)算法設(shè)置參數(shù)。對(duì)于“創(chuàng)建網(wǎng)格”你可以將“間距”和“輸出CRS”設(shè)置為模型參數(shù)這樣每次運(yùn)行模型時(shí)都可以動(dòng)態(tài)輸入。保存模型。以后你只需要運(yùn)行這個(gè)模型輸入范圍圖層、間距和CRS就能一鍵得到裁剪好并帶ID的最終網(wǎng)格。5.2 使用PyQGIS腳本批量生成以下是一個(gè)示例PyQGIS腳本演示了如何在Python控制臺(tái)或外部腳本中自動(dòng)化完成網(wǎng)格創(chuàng)建。這個(gè)腳本假設(shè)你已經(jīng)有一個(gè)名為study_area的面圖層矢量層并希望以指定邊長(zhǎng)和CRS創(chuàng)建網(wǎng)格。# 導(dǎo)入必要的模塊 from qgis.core import (QgsProject, QgsVectorLayer, QgsGeometry, QgsRectangle, QgsCoordinateReferenceSystem, QgsFields, QgsField, QgsFeature, QgsProcessingFeedback, QgsProcessingContext) from qgis.analysis import QgsNativeAlgorithms import processing from processing.core.Processing import Processing # 初始化Processing框架 Processing.initialize() # 1. 定義參數(shù) input_layer_name study_area # 你的范圍圖層名 grid_size 500 # 網(wǎng)格邊長(zhǎng)單位取決于輸出CRS output_crs_code EPSG:32650 # 輸出坐標(biāo)系例如UTM 50N output_path /path/to/your/output/grid_500m.gpkg # 輸出文件路徑 # 2. 獲取輸入圖層 project QgsProject.instance() input_layer project.mapLayersByName(input_layer_name)[0] if not input_layer: print(f錯(cuò)誤未找到名為 {input_layer_name} 的圖層。) exit() # 3. 獲取輸入圖層的范圍在其自身坐標(biāo)系下 extent input_layer.extent() xmin, ymin, xmax, ymax extent.xMinimum(), extent.yMinimum(), extent.xMaximum(), extent.yMaximum() # 4. 設(shè)置算法參數(shù) params { TYPE: 2, # 2 代表矩形多邊形 EXTENT: f{xmin},{xmax},{ymin},{ymax}, HSPACING: grid_size, VSPACING: grid_size, HOVERLAP: 0, VOVERLAP: 0, CRS: output_crs_code, OUTPUT: memory: # 先輸出到臨時(shí)內(nèi)存層 } # 5. 執(zhí)行“創(chuàng)建網(wǎng)格”算法 feedback QgsProcessingFeedback() context QgsProcessingContext() result processing.run(native:creategrid, params, feedbackfeedback, contextcontext) grid_layer result[OUTPUT] # 6. 可選裁剪網(wǎng)格到輸入圖層范圍 if True: # 如果你想裁剪將此設(shè)為T(mén)rue clip_params { INPUT: grid_layer, OVERLAY: input_layer, OUTPUT: memory: } result_clip processing.run(native:clip, clip_params, feedbackfeedback, contextcontext) final_layer result_clip[OUTPUT] else: final_layer grid_layer # 7. 為最終圖層添加ID字段并賦值 final_layer.startEditing() provider final_layer.dataProvider() if final_layer.fields().indexFromName(grid_id) -1: provider.addAttributes([QgsField(grid_id, QVariant.Int)]) final_layer.updateFields() idx final_layer.fields().indexFromName(grid_id) for i, feat in enumerate(final_layer.getFeatures()): final_layer.changeAttributeValue(feat.id(), idx, i1) final_layer.commitChanges() # 8. 將結(jié)果保存到文件并加載到項(xiàng)目 options QgsVectorFileWriter.SaveVectorOptions() options.driverName GPKG options.fileEncoding UTF-8 error QgsVectorFileWriter.writeAsVectorFormatV2(final_layer, output_path, QgsCoordinateTransformContext(), options) if error[0] QgsVectorFileWriter.NoError: # 加載保存的圖層 vlayer QgsVectorLayer(output_path, Final_Grid_500m, ogr) if vlayer.isValid(): project.addMapLayer(vlayer) print(f成功網(wǎng)格已保存并加載{output_path}) else: print(加載生成的圖層失敗。) else: print(f保存文件時(shí)出錯(cuò){error[1]})這個(gè)腳本展示了完整的流程獲取范圍、創(chuàng)建網(wǎng)格、裁剪、添加屬性。你可以將其保存為.py文件在QGIS的Python控制臺(tái)中運(yùn)行或通過(guò)“處理”-“腳本”-“創(chuàng)建新腳本”將其集成到處理工具箱中實(shí)現(xiàn)一鍵運(yùn)行。5.3 處理大型區(qū)域與性能優(yōu)化當(dāng)研究區(qū)域非常大如一個(gè)省或國(guó)家創(chuàng)建小尺寸網(wǎng)格如100米會(huì)產(chǎn)生海量要素?cái)?shù)百萬(wàn)甚至上千萬(wàn)個(gè)多邊形可能導(dǎo)致QGIS卡頓甚至崩潰。優(yōu)化策略分塊處理將大區(qū)域劃分為多個(gè)子區(qū)域如按市、縣分別創(chuàng)建網(wǎng)格后再合并。可以使用“分割矢量圖層”工具按屬性分割范圍圖層。使用數(shù)據(jù)庫(kù)將輸出直接保存到PostGIS等空間數(shù)據(jù)庫(kù)中而不是Shapefile或GeoPackage文件。數(shù)據(jù)庫(kù)在處理海量要素時(shí)性能更優(yōu)并支持空間索引。調(diào)整網(wǎng)格尺寸在滿足分析精度的前提下適當(dāng)增大網(wǎng)格邊長(zhǎng)能指數(shù)級(jí)減少要素?cái)?shù)量。在創(chuàng)建后簡(jiǎn)化屬性創(chuàng)建網(wǎng)格時(shí)默認(rèn)會(huì)帶有id、left、top、right、bottom等字段。如果不需要可以在創(chuàng)建后刪除這些字段減少文件大小。通過(guò)理解工具原理、規(guī)避常見(jiàn)陷阱、并掌握自動(dòng)化方法你就能將QGIS的“創(chuàng)建網(wǎng)格”功能從一個(gè)簡(jiǎn)單的圖形生成器轉(zhuǎn)變?yōu)橹文憧臻g分析工作的強(qiáng)大、可靠且高效的核心組件。記住精準(zhǔn)的網(wǎng)格是空間量化和模型化的基石多花幾分鐘在坐標(biāo)系和參數(shù)設(shè)置上能為后續(xù)所有分析省去數(shù)小時(shí)的糾錯(cuò)時(shí)間。 QGIS創(chuàng)建正方形網(wǎng)格：從坐標(biāo)系選擇到自動(dòng)化腳本全解析
1. 項(xiàng)目緣起為什么我們需要在QGIS中創(chuàng)建正方形網(wǎng)格在地理信息處理的實(shí)際工作中我們常常會(huì)遇到一個(gè)看似簡(jiǎn)單卻至關(guān)重要的需求如何將一個(gè)不規(guī)則的研究區(qū)域高效、精確地分割成一個(gè)個(gè)規(guī)整的單元無(wú)論是進(jìn)行城市地塊的抽樣調(diào)查、生態(tài)學(xué)中的樣方分析、還是區(qū)域性的資源評(píng)估將連續(xù)的空間離散化為規(guī)則的網(wǎng)格都是后續(xù)統(tǒng)計(jì)分析、模型構(gòu)建和可視化表達(dá)的基礎(chǔ)。手動(dòng)繪制這些網(wǎng)格不僅效率低下而且難以保證幾何精度和拓?fù)湟恢滦?。QGIS作為一款功能強(qiáng)大的開(kāi)源GIS桌面軟件其“創(chuàng)建網(wǎng)格”工具正是為此而生。它允許用戶基于指定的范圍、形狀如矩形、菱形、六邊形和尺寸批量生成覆蓋目標(biāo)區(qū)域的規(guī)則幾何圖形。其中正方形網(wǎng)格因其計(jì)算簡(jiǎn)便、易于理解和空間索引高效成為最常用的網(wǎng)格類(lèi)型。然而從“知道有這個(gè)工具”到“能穩(wěn)定產(chǎn)出符合預(yù)期、坐標(biāo)正確的網(wǎng)格”中間往往隔著對(duì)坐標(biāo)系、參數(shù)設(shè)置和輸出結(jié)果處理的深刻理解。很多新手甚至是有一定經(jīng)驗(yàn)的用戶都可能在“創(chuàng)建網(wǎng)格”這一步就踩坑比如生成的網(wǎng)格大小不對(duì)、覆蓋范圍有偏差或者后續(xù)疊加數(shù)據(jù)時(shí)出現(xiàn)錯(cuò)位。本文將從一個(gè)資深GIS從業(yè)者的視角手把手帶你深入QGIS的“創(chuàng)建網(wǎng)格”功能。我們不僅會(huì)完成“點(diǎn)擊按鈕生成網(wǎng)格”這個(gè)動(dòng)作更會(huì)拆解其背后的每一個(gè)參數(shù)邏輯特別是坐標(biāo)系如EPSG:3857的選擇與影響并分享在實(shí)際項(xiàng)目中如何避免常見(jiàn)陷阱確保你得到的每一個(gè)正方形網(wǎng)格都精準(zhǔn)、可用。2. 核心工具解析“創(chuàng)建網(wǎng)格”的功能界面與參數(shù)全解“創(chuàng)建網(wǎng)格”工具位于QGIS的“矢量”菜單下的“研究工具”子菜單中。打開(kāi)后你會(huì)看到一個(gè)參數(shù)對(duì)話框。這個(gè)對(duì)話框里的每一個(gè)選項(xiàng)都直接決定了最終輸出網(wǎng)格的形態(tài)和空間屬性。我們來(lái)逐一拆解并解釋其背后的“為什么”。2.1 網(wǎng)格類(lèi)型與范圍定義首先在“網(wǎng)格類(lèi)型”下拉菜單中選擇“矩形多邊形”。這里有一個(gè)關(guān)鍵點(diǎn)雖然我們最終想要的是正方形但QGIS是通過(guò)指定“寬度”和“高度”來(lái)定義矩形的。當(dāng)我們將寬度和高度設(shè)置為相同的數(shù)值時(shí)得到的就是正方形。選擇“多邊形”而非“線”是因?yàn)槎噙呅尉W(wǎng)格才能作為一個(gè)獨(dú)立的面要素進(jìn)行后續(xù)的面積計(jì)算、空間連接等操作。接下來(lái)是“網(wǎng)格范圍”。這是第一個(gè)容易出錯(cuò)的環(huán)節(jié)。你有四種方式定義范圍畫(huà)布范圍直接使用當(dāng)前地圖視圖顯示的范圍。這非常方便但風(fēng)險(xiǎn)在于如果你縮放或平移了地圖范圍就變了可重復(fù)性差。圖層范圍選擇某個(gè)已有矢量圖層如一個(gè)表示研究區(qū)的面圖層工具會(huì)計(jì)算該圖層的整體外包矩形即能完全包含該圖層所有要素的最小矩形。這是最常用且推薦的方式因?yàn)樗_保了網(wǎng)格能完全覆蓋你的目標(biāo)區(qū)域。書(shū)簽范圍如果你之前保存了地圖書(shū)簽可以調(diào)用。手動(dòng)輸入坐標(biāo)直接輸入左下角和右上角的X、Y坐標(biāo)。這種方式最精確但需要你事先知道精確的范圍坐標(biāo)。注意無(wú)論選擇哪種方式這里定義的范圍都是一個(gè)矩形。如果你的研究區(qū)形狀極不規(guī)則比如一個(gè)彎曲的河流流域使用其外包矩形創(chuàng)建網(wǎng)格會(huì)在研究區(qū)外圍產(chǎn)生大量“無(wú)效”的網(wǎng)格單元。后續(xù)通常需要用一個(gè)“按位置選擇”或“裁剪”工具將落在實(shí)際研究區(qū)之外的網(wǎng)格剔除。2.2 網(wǎng)格間距理解“寬度”與“高度”的單位這是核心參數(shù)。在“水平間距”和“垂直間距”中分別輸入你想要的網(wǎng)格單元的寬度和高度。例如輸入“100”期望得到邊長(zhǎng)為100單位的正方形。這里的坑在于“單位”。這個(gè)“100”單位是什么米公里度這完全取決于你為輸出網(wǎng)格所選擇的坐標(biāo)系CRS。如果輸出CRS是地理坐標(biāo)系如EPSG:4326, WGS84單位是十進(jìn)制度。在地球表面1度經(jīng)度的實(shí)際長(zhǎng)度隨緯度變化在赤道最長(zhǎng)向兩極遞減至01度緯度的長(zhǎng)度相對(duì)恒定約111公里。此時(shí)設(shè)置“100”度將得到一個(gè)在現(xiàn)實(shí)世界中巨大且扭曲的“正方形”這通常不是我們想要的。如果輸出CRS是投影坐標(biāo)系如EPSG:3857 Web Mercator, 或EPSG:32649 UTM Zone 49N單位通常是米。此時(shí)設(shè)置“100”就意味著邊長(zhǎng)100米的正方形。這才是進(jìn)行面積計(jì)算和空間分析的理想選擇。因此在輸入間距值之前你必須明確你希望網(wǎng)格在現(xiàn)實(shí)世界中的實(shí)際尺寸例如100米 x 100米并確保為輸出圖層選擇了合適的、以米為單位的投影坐標(biāo)系。工具界面中的“輸出CRS”選項(xiàng)就是用來(lái)定義這一點(diǎn)的。2.3 輸出CRS坐標(biāo)系選擇的深遠(yuǎn)影響“輸出CRS”是決定網(wǎng)格空間正確性的靈魂參數(shù)。它有兩個(gè)主要選項(xiàng)項(xiàng)目CRS使用當(dāng)前QGIS項(xiàng)目設(shè)置的坐標(biāo)系。指定CRS手動(dòng)選擇一個(gè)坐標(biāo)系。強(qiáng)烈建議選擇“指定CRS”并選擇一個(gè)適合你研究區(qū)域的投影坐標(biāo)系。原因如下保證尺寸準(zhǔn)確如上所述投影坐標(biāo)系使用線性單位米、英尺能保證網(wǎng)格邊長(zhǎng)在地圖上的測(cè)量值與實(shí)際地面距離對(duì)應(yīng)。避免幾何變形正方形網(wǎng)格在投影坐標(biāo)系下才能保持較好的形狀。在地理坐標(biāo)系下由于將球面展開(kāi)為平面帶來(lái)的變形網(wǎng)格會(huì)變成梯形或其他形狀。與底圖/業(yè)務(wù)數(shù)據(jù)對(duì)齊如果你的在線底圖如OSM、谷歌地圖或業(yè)務(wù)數(shù)據(jù)使用的是Web MercatorEPSG:3857那么將網(wǎng)格也創(chuàng)建為相同的坐標(biāo)系可以確保完美的疊加沒(méi)有偏移。關(guān)于EPSG:3857這是一個(gè)全球性的投影坐標(biāo)系被幾乎所有網(wǎng)絡(luò)地圖服務(wù)使用。它的單位是“米”但在高緯度地區(qū)面積和形狀變形非常嚴(yán)重。如果你的研究區(qū)域范圍不大比如一個(gè)城市且主要為了與在線地圖疊加顯示使用3857是方便的選擇。但如果需要進(jìn)行精確的面積計(jì)算或大范圍的空間分析應(yīng)優(yōu)先使用更優(yōu)的局部投影坐標(biāo)系如UTM通用橫軸墨卡托分區(qū)。你可以通過(guò)QGIS的“按圖層范圍自動(dòng)選擇”功能或查詢當(dāng)?shù)氐臏y(cè)繪標(biāo)準(zhǔn)來(lái)確定最適合的UTM帶號(hào)如EPSG:32650 for UTM Zone 50N。2.4 其他參數(shù)與輸出水平/垂直格網(wǎng)數(shù)這是一個(gè)與“間距”互斥的參數(shù)。如果你知道范圍的總寬度是1000米并想分成10列那么可以設(shè)置“水平格網(wǎng)數(shù)”為10工具會(huì)自動(dòng)計(jì)算間距為100米。我個(gè)人更傾向于直接使用“間距”因?yàn)槲覍?duì)實(shí)地尺寸更敏感。網(wǎng)格對(duì)齊點(diǎn)可以設(shè)置網(wǎng)格的起始偏移點(diǎn)。通常不需要調(diào)整除非你有特殊的對(duì)齊要求如讓網(wǎng)格的角點(diǎn)與某個(gè)已知坐標(biāo)點(diǎn)對(duì)齊。輸出指定生成網(wǎng)格的保存路徑和格式如GeoPackage、Shapefile。點(diǎn)擊“運(yùn)行”一個(gè)全新的、覆蓋指定范圍的規(guī)則網(wǎng)格圖層就會(huì)加載到你的地圖中。3. 實(shí)戰(zhàn)演練從零生成一個(gè)用于抽樣分析的正方形網(wǎng)格讓我們通過(guò)一個(gè)完整的模擬案例將上述理論付諸實(shí)踐。假設(shè)我們有一個(gè)“某市森林公園邊界.shp”的面數(shù)據(jù)我們需要在其范圍內(nèi)創(chuàng)建邊長(zhǎng)為500米的正方形網(wǎng)格用于規(guī)劃生態(tài)監(jiān)測(cè)樣地。3.1 數(shù)據(jù)準(zhǔn)備與坐標(biāo)系檢查首先將“森林公園邊界.shp”加載到QGIS中。右鍵點(diǎn)擊圖層選擇“屬性”切換到“信息”選項(xiàng)卡。查看“坐標(biāo)系”信息。假設(shè)它顯示為“EPSG:4547”這是一個(gè)適用于該地區(qū)的投影坐標(biāo)系單位是米。這很好我們的數(shù)據(jù)本身就在一個(gè)合適的投影坐標(biāo)系中。如果數(shù)據(jù)是地理坐標(biāo)系如EPSG:4326你需要先考慮是否要將其重新投影到一個(gè)投影坐標(biāo)系。你可以在導(dǎo)入時(shí)使用“圖層CRS”進(jìn)行設(shè)置或者后續(xù)使用“導(dǎo)出”或“重投影圖層”工具來(lái)處理?；驹瓌t是用于創(chuàng)建網(wǎng)格的“范圍圖層”和“輸出網(wǎng)格”最好在同一個(gè)、以米為單位的投影坐標(biāo)系下以避免不必要的轉(zhuǎn)換誤差。3.2 執(zhí)行“創(chuàng)建網(wǎng)格”工具打開(kāi)“矢量” - “研究工具” - “創(chuàng)建網(wǎng)格”。網(wǎng)格類(lèi)型選擇“矩形多邊形”。網(wǎng)格范圍選擇“圖層范圍”并在下拉菜單中選擇“森林公園邊界”。水平間距輸入500。垂直間距輸入500。輸出CRS選擇“指定CRS”。點(diǎn)擊右側(cè)的“選擇CRS”按鈕地球圖標(biāo)。在過(guò)濾器中輸入“4547”選擇它。關(guān)鍵點(diǎn)這里我們選擇與源數(shù)據(jù)相同的CRS確保網(wǎng)格與原始邊界完美套合。輸出點(diǎn)擊“...”按鈕將結(jié)果保存為GeoPackage格式命名為“森林公園_500米網(wǎng)格.gpkg”。點(diǎn)擊“運(yùn)行”。瞬間一個(gè)整齊的網(wǎng)格圖層就覆蓋在了森林公園之上。3.3 結(jié)果驗(yàn)證與后處理生成網(wǎng)格后不能直接使用必須進(jìn)行驗(yàn)證和清理。驗(yàn)證1尺寸檢查。使用“測(cè)量工具”設(shè)置為“米”為單位隨機(jī)測(cè)量幾個(gè)網(wǎng)格的邊長(zhǎng)。確認(rèn)它們確實(shí)是500米左右由于范圍除間距可能不是整數(shù)邊緣網(wǎng)格可能略小。驗(yàn)證2范圍檢查。觀察網(wǎng)格是否完全覆蓋了森林公園。由于我們用的是圖層的“外包矩形”網(wǎng)格范圍肯定會(huì)比實(shí)際公園邊界大四周會(huì)多出很多覆蓋了公園外部區(qū)域的網(wǎng)格。后處理裁剪網(wǎng)格。我們需要只保留落在森林公園內(nèi)部的網(wǎng)格。使用“矢量” - “地理處理工具” - “裁剪”工具?！拜斎雸D層”選擇我們剛生成的“森林公園_500米網(wǎng)格”?！隘B加圖層”選擇“森林公園邊界”。運(yùn)行。得到的新圖層“裁剪后的網(wǎng)格”就只包含完全位于公園內(nèi)部的完整或部分網(wǎng)格了。后處理添加唯一標(biāo)識(shí)。為了方便管理為每個(gè)網(wǎng)格賦予一個(gè)ID。打開(kāi)“裁剪后的網(wǎng)格”圖層的屬性表。點(diǎn)擊“打開(kāi)字段計(jì)算器”創(chuàng)建一個(gè)新字段例如叫g(shù)rid_id。在表達(dá)式框中輸入row_number這將為每一行即每個(gè)網(wǎng)格生成一個(gè)從1開(kāi)始的序列號(hào)。點(diǎn)擊“確定”。現(xiàn)在每個(gè)網(wǎng)格都有了自己的編號(hào)。至此一套可用于布設(shè)監(jiān)測(cè)樣地、進(jìn)行分區(qū)統(tǒng)計(jì)的標(biāo)準(zhǔn)化正方形網(wǎng)格就制作完成了。你可以基于grid_id來(lái)隨機(jī)抽樣或者將其他數(shù)據(jù)如植被類(lèi)型、動(dòng)物觀測(cè)點(diǎn)與這些網(wǎng)格進(jìn)行空間連接計(jì)算每個(gè)網(wǎng)格內(nèi)的統(tǒng)計(jì)值。4. 深度避坑坐標(biāo)系3857 vs. 局部投影帶來(lái)的隱秘陷阱“創(chuàng)建網(wǎng)格”操作中最頑固、最不易察覺(jué)的問(wèn)題幾乎都與坐標(biāo)系有關(guān)。我們重點(diǎn)剖析EPSG:3857這個(gè)“熟悉的陌生人”。4.1 EPSG:3857的便利性與局限性便利性EPSG:3857Web Mercator是網(wǎng)絡(luò)地圖的“世界語(yǔ)”。幾乎所有的在線瓦片地圖Google Maps, OpenStreetMap, 天地圖等都使用它。在QGIS中加載XYZ Tiles連接這些在線地圖時(shí)項(xiàng)目CRS會(huì)自動(dòng)切換到3857。如果你在此狀態(tài)下直接用畫(huà)布范圍創(chuàng)建網(wǎng)格并選擇“項(xiàng)目CRS”作為輸出那么網(wǎng)格會(huì)與底圖完美貼合視覺(jué)上非常整齊。局限性即大坑嚴(yán)重的面積變形Web Mercator投影在高緯度地區(qū)會(huì)極度夸大面積。例如在赤道附近一個(gè)3857坐標(biāo)系下邊長(zhǎng)為100米的網(wǎng)格其實(shí)際地面面積接近1公頃。但在北緯60度地區(qū)同一個(gè)“100米”網(wǎng)格所代表的實(shí)際地面面積會(huì)小得多。這意味著如果你在3857下創(chuàng)建了“500米網(wǎng)格”并用于面積相關(guān)的分析如計(jì)算網(wǎng)格內(nèi)的森林面積占比結(jié)果將是完全錯(cuò)誤的。非等角投影雖然它被宣傳為“保角”但實(shí)際上在非標(biāo)準(zhǔn)緯線上并不完全保角網(wǎng)格的形狀也會(huì)發(fā)生微小變形并非嚴(yán)格的正方形。4.2 一個(gè)真實(shí)的對(duì)比實(shí)驗(yàn)假設(shè)我們?cè)诒本┑貐^(qū)約北緯40度創(chuàng)建一個(gè)邊長(zhǎng)為1000米的網(wǎng)格。方案A錯(cuò)誤示范項(xiàng)目CRS為EPSG:3857加載OSM底圖以當(dāng)前視圖范圍創(chuàng)建網(wǎng)格輸出CRS使用“項(xiàng)目CRS”即3857。方案B正確做法將項(xiàng)目CRS設(shè)置為更適合該區(qū)域的投影如EPSG:32650UTM Zone 50N。重新加載數(shù)據(jù)或?qū)?shù)據(jù)重投影然后以圖層范圍創(chuàng)建網(wǎng)格輸出CRS指定為EPSG:32650。對(duì)比結(jié)果視覺(jué)對(duì)齊方案A的網(wǎng)格與OSM底圖嚴(yán)絲合縫。方案B的網(wǎng)格在OSM底圖3857上會(huì)有可見(jiàn)的旋轉(zhuǎn)和偏移因?yàn)榈讏D被動(dòng)態(tài)投影到了3857顯示而網(wǎng)格是32650的。這常常讓新手困惑認(rèn)為B方案“錯(cuò)了”。實(shí)際上視覺(jué)上的錯(cuò)位是投影顯示造成的B方案的數(shù)據(jù)在它自己的坐標(biāo)系下才是正確的。面積真相使用QGIS的“字段計(jì)算器”和$area函數(shù)計(jì)算兩個(gè)方案下網(wǎng)格的面積。方案A3857網(wǎng)格$area返回的值接近 1,000,000 平方米即1平方公里但這只是地圖上的“投影面積”。方案B32650網(wǎng)格$area返回的值也接近 1,000,000 平方米。關(guān)鍵來(lái)了使用“測(cè)量工具”設(shè)置為“橢球體”模式如WGS84分別測(cè)量?jī)蓚€(gè)網(wǎng)格的實(shí)際地面距離。你會(huì)發(fā)現(xiàn)方案B的網(wǎng)格邊長(zhǎng)更接近真實(shí)的1000米而方案A的網(wǎng)格實(shí)際地面邊長(zhǎng)要小于1000米。其“投影面積”1平方公里是虛高的。實(shí)操心得永遠(yuǎn)不要為了與在線底圖視覺(jué)對(duì)齊的便利而犧牲空間數(shù)據(jù)的幾何精度。正確的 workflow 是1) 為你的分析確定一個(gè)合適的局部投影坐標(biāo)系CRS。2) 將所有業(yè)務(wù)數(shù)據(jù)包括用于定義范圍的圖層都轉(zhuǎn)換或創(chuàng)建到這個(gè)CRS下。3) 在此CRS下創(chuàng)建網(wǎng)格、進(jìn)行分析計(jì)算。4) 如果需要在線底圖作為參考可以添加X(jué)YZ TilesQGIS會(huì)動(dòng)態(tài)將其投影到你的項(xiàng)目CRS中顯示此時(shí)底圖可能會(huì)有變形和偏移但你的業(yè)務(wù)數(shù)據(jù)之間的空間關(guān)系是精確的。記住底圖是“參考背景”你的分析數(shù)據(jù)才是“主角”。4.3 處理“創(chuàng)建網(wǎng)格”時(shí)常見(jiàn)的坐標(biāo)系報(bào)錯(cuò)有時(shí)在運(yùn)行工具時(shí)會(huì)遇到“圖層范圍在輸出CRS中無(wú)效”或類(lèi)似錯(cuò)誤。這通常是因?yàn)椤胺秶鷪D層”的坐標(biāo)系與“輸出CRS”不兼容或者范圍坐標(biāo)值超出了目標(biāo)CRS的有效域例如試圖在UTM投影中創(chuàng)建一個(gè)經(jīng)度跨度過(guò)大的網(wǎng)格。解決方法統(tǒng)一坐標(biāo)系確保范圍圖層和輸出CRS一致或至少是能夠相互轉(zhuǎn)換的。檢查范圍值在范圍圖層的屬性中查看其實(shí)際坐標(biāo)范圍判斷是否合理。使用項(xiàng)目CRS一個(gè)穩(wěn)妥的方法是先將項(xiàng)目和所有相關(guān)圖層的CRS都設(shè)置為你的目標(biāo)分析CRS然后使用“畫(huà)布范圍”或“圖層范圍”來(lái)創(chuàng)建網(wǎng)格。5. 進(jìn)階應(yīng)用與腳本化讓網(wǎng)格創(chuàng)建融入自動(dòng)化工作流對(duì)于需要頻繁創(chuàng)建不同尺寸網(wǎng)格或作為復(fù)雜模型一部分的用戶手動(dòng)點(diǎn)擊GUI效率太低。QGIS的強(qiáng)大之處在于其處理框架和PyQGIS允許我們將這個(gè)過(guò)程腳本化。5.1 使用處理模型設(shè)計(jì)器構(gòu)建可復(fù)用流程我們可以將“創(chuàng)建網(wǎng)格”和后續(xù)的“裁剪”、“添加ID”等步驟打包成一個(gè)自定義模型。打開(kāi)“處理” - “圖形模型設(shè)計(jì)器”。從右側(cè)算法列表中將“創(chuàng)建網(wǎng)格”、“裁剪”、“字段計(jì)算器”拖入畫(huà)布。進(jìn)行連接將“范圍圖層”作為“創(chuàng)建網(wǎng)格”和“裁剪”的輸入將“創(chuàng)建網(wǎng)格”的輸出連接到“裁剪”的輸入圖層將“裁剪”的輸出連接到“字段計(jì)算器”的輸入。為每個(gè)算法設(shè)置參數(shù)。對(duì)于“創(chuàng)建網(wǎng)格”你可以將“間距”和“輸出CRS”設(shè)置為模型參數(shù)這樣每次運(yùn)行模型時(shí)都可以動(dòng)態(tài)輸入。保存模型。以后你只需要運(yùn)行這個(gè)模型輸入范圍圖層、間距和CRS就能一鍵得到裁剪好并帶ID的最終網(wǎng)格。5.2 使用PyQGIS腳本批量生成以下是一個(gè)示例PyQGIS腳本演示了如何在Python控制臺(tái)或外部腳本中自動(dòng)化完成網(wǎng)格創(chuàng)建。這個(gè)腳本假設(shè)你已經(jīng)有一個(gè)名為study_area的面圖層矢量層并希望以指定邊長(zhǎng)和CRS創(chuàng)建網(wǎng)格。# 導(dǎo)入必要的模塊 from qgis.core import (QgsProject, QgsVectorLayer, QgsGeometry, QgsRectangle, QgsCoordinateReferenceSystem, QgsFields, QgsField, QgsFeature, QgsProcessingFeedback, QgsProcessingContext) from qgis.analysis import QgsNativeAlgorithms import processing from processing.core.Processing import Processing # 初始化Processing框架 Processing.initialize() # 1. 定義參數(shù) input_layer_name study_area # 你的范圍圖層名 grid_size 500 # 網(wǎng)格邊長(zhǎng)單位取決于輸出CRS output_crs_code EPSG:32650 # 輸出坐標(biāo)系例如UTM 50N output_path /path/to/your/output/grid_500m.gpkg # 輸出文件路徑 # 2. 獲取輸入圖層 project QgsProject.instance() input_layer project.mapLayersByName(input_layer_name)[0] if not input_layer: print(f錯(cuò)誤未找到名為 {input_layer_name} 的圖層。) exit() # 3. 獲取輸入圖層的范圍在其自身坐標(biāo)系下 extent input_layer.extent() xmin, ymin, xmax, ymax extent.xMinimum(), extent.yMinimum(), extent.xMaximum(), extent.yMaximum() # 4. 設(shè)置算法參數(shù) params { TYPE: 2, # 2 代表矩形多邊形 EXTENT: f{xmin},{xmax},{ymin},{ymax}, HSPACING: grid_size, VSPACING: grid_size, HOVERLAP: 0, VOVERLAP: 0, CRS: output_crs_code, OUTPUT: memory: # 先輸出到臨時(shí)內(nèi)存層 } # 5. 執(zhí)行“創(chuàng)建網(wǎng)格”算法 feedback QgsProcessingFeedback() context QgsProcessingContext() result processing.run(native:creategrid, params, feedbackfeedback, contextcontext) grid_layer result[OUTPUT] # 6. 可選裁剪網(wǎng)格到輸入圖層范圍 if True: # 如果你想裁剪將此設(shè)為T(mén)rue clip_params { INPUT: grid_layer, OVERLAY: input_layer, OUTPUT: memory: } result_clip processing.run(native:clip, clip_params, feedbackfeedback, contextcontext) final_layer result_clip[OUTPUT] else: final_layer grid_layer # 7. 為最終圖層添加ID字段并賦值 final_layer.startEditing() provider final_layer.dataProvider() if final_layer.fields().indexFromName(grid_id) -1: provider.addAttributes([QgsField(grid_id, QVariant.Int)]) final_layer.updateFields() idx final_layer.fields().indexFromName(grid_id) for i, feat in enumerate(final_layer.getFeatures()): final_layer.changeAttributeValue(feat.id(), idx, i1) final_layer.commitChanges() # 8. 將結(jié)果保存到文件并加載到項(xiàng)目 options QgsVectorFileWriter.SaveVectorOptions() options.driverName GPKG options.fileEncoding UTF-8 error QgsVectorFileWriter.writeAsVectorFormatV2(final_layer, output_path, QgsCoordinateTransformContext(), options) if error[0] QgsVectorFileWriter.NoError: # 加載保存的圖層 vlayer QgsVectorLayer(output_path, Final_Grid_500m, ogr) if vlayer.isValid(): project.addMapLayer(vlayer) print(f成功網(wǎng)格已保存并加載{output_path}) else: print(加載生成的圖層失敗。) else: print(f保存文件時(shí)出錯(cuò){error[1]})這個(gè)腳本展示了完整的流程獲取范圍、創(chuàng)建網(wǎng)格、裁剪、添加屬性。你可以將其保存為.py文件在QGIS的Python控制臺(tái)中運(yùn)行或通過(guò)“處理”-“腳本”-“創(chuàng)建新腳本”將其集成到處理工具箱中實(shí)現(xiàn)一鍵運(yùn)行。5.3 處理大型區(qū)域與性能優(yōu)化當(dāng)研究區(qū)域非常大如一個(gè)省或國(guó)家創(chuàng)建小尺寸網(wǎng)格如100米會(huì)產(chǎn)生海量要素?cái)?shù)百萬(wàn)甚至上千萬(wàn)個(gè)多邊形可能導(dǎo)致QGIS卡頓甚至崩潰。優(yōu)化策略分塊處理將大區(qū)域劃分為多個(gè)子區(qū)域如按市、縣分別創(chuàng)建網(wǎng)格后再合并。可以使用“分割矢量圖層”工具按屬性分割范圍圖層。使用數(shù)據(jù)庫(kù)將輸出直接保存到PostGIS等空間數(shù)據(jù)庫(kù)中而不是Shapefile或GeoPackage文件。數(shù)據(jù)庫(kù)在處理海量要素時(shí)性能更優(yōu)并支持空間索引。調(diào)整網(wǎng)格尺寸在滿足分析精度的前提下適當(dāng)增大網(wǎng)格邊長(zhǎng)能指數(shù)級(jí)減少要素?cái)?shù)量。在創(chuàng)建后簡(jiǎn)化屬性創(chuàng)建網(wǎng)格時(shí)默認(rèn)會(huì)帶有id、left、top、right、bottom等字段。如果不需要可以在創(chuàng)建后刪除這些字段減少文件大小。通過(guò)理解工具原理、規(guī)避常見(jiàn)陷阱、并掌握自動(dòng)化方法你就能將QGIS的“創(chuàng)建網(wǎng)格”功能從一個(gè)簡(jiǎn)單的圖形生成器轉(zhuǎn)變?yōu)橹文憧臻g分析工作的強(qiáng)大、可靠且高效的核心組件。記住精準(zhǔn)的網(wǎng)格是空間量化和模型化的基石多花幾分鐘在坐標(biāo)系和參數(shù)設(shè)置上能為后續(xù)所有分析省去數(shù)小時(shí)的糾錯(cuò)時(shí)間。