也不知道為啥,向量裁剪向量這麼普通的東西這麼難找,趕緊放出來讓大家用用
import os
import numpy as np
import geopandas as gpd
import warnings
warnings.filterwarnings(
'ignore'
,'geoseries.notna'
, userwarning)
shp_path =
'./osm/gis_osm_roads_free_1.shp'
#被裁剪的物件
mask =
'd:/toolset/test_data2/mask.shp'
#用於裁剪的範圍
shp_gpd = gpd.read_file(shp_path)
mask_gpd = gpd.read_file(mask)
print
("shp_gpd"
, shp_gpd.crs)
mask_clip = gpd.clip(shp_gpd, mask_gpd)
mask_clip = mask_clip[
~mask_clip.is_empty]
mask_clip.to_file(
"./osm/osm_road_anji.shp"
)#儲存結果
使用ogr裁剪向量資料
使用ogr裁剪向量資料 由來 近期有個需求,內容是這樣的 我們有兩個向量資料,現在要求以乙個向量檔案為底板,按字段對另乙個向量檔案進行分割,生成若干小的shpfile檔案 分析 經過分析之後,將步驟拆解如下 首先確保兩個shpfile投影座標系統一 如果出現不統一的情況,那麼用arcgis的工具pr...
狸貓的面試 專案描述 向量裁剪
專案描述 要求 在4s內,用單執行緒,完成對百萬級的向量圖形的裁剪 輸入 1.乙個多邊形 儲存在xml檔案中,順序儲存多邊形所有的端點 2.乙個向量圖形 line或者circle line儲存其兩個端點ab,circle儲存其圓心o和半徑r。實際有一百萬條直線和一百萬個圓 輸出 輸出裁剪後的圖形。例...
IDL實現向量 shp 裁剪柵格TASK(一
隨著envi idl版本的更新,idl對向量和柵格資料的處理也變得越來越簡單化。其提供了很多方便的介面,使得使用者呼叫和學習練習便捷成為了可能。最近接觸idl,發現好多網上的 都是延後的,新的介面 理解和編寫起來都比較方便,尤其是在做大量資料研究和應用時,使用批處理的方式顯得尤其重要。新的介面還在摸...