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