資訊內容
Python地理數據處理之分析使用GR進行矢量
wkz少兒編程網-Scratch_Python_教程_免費兒童編程學習平臺免費學習推薦：python視頻教程wkz少兒編程網-Scratch_Python_教程_免費兒童編程學習平臺
1、疊加分析wkz少兒編程網-Scratch_Python_教程_免費兒童編程學習平臺
??疊加分析操作：
??plot顏色：‘r’ 紅色， ‘g’ 綠色， ‘b’ 藍色， ‘c’ 青色， ‘y’ 黃色， ‘m’ 品紅， ‘k’ 黑色， ‘w’ 白色。wkz少兒編程網-Scratch_Python_教程_免費兒童編程學習平臺
??新奧爾良城市邊界、水體和濕地的簡單地圖：wkz少兒編程網-Scratch_Python_教程_免費兒童編程學習平臺
??1.新奧爾良城市沼澤區域分析：wkz少兒編程網-Scratch_Python_教程_免費兒童編程學習平臺
??2.計算城市的濕地面積：wkz少兒編程網-Scratch_Python_教程_免費兒童編程學習平臺
??注：通過空間過濾和屬性過濾，將不必要的要素過濾，這樣可以顯著減少處理時間。wkz少兒編程網-Scratch_Python_教程_免費兒童編程學習平臺
??3.兩圖層求交：wkz少兒編程網-Scratch_Python_教程_免費兒童編程學習平臺
# 將湖泊數據排除# 在內存中創建一個臨時圖層# 將圖層相交，將結果儲存在臨時圖層中water_lyr.SetAttributeFilter("Feature != 'Lake'")water_lyr.SetSpatialFilter(nola_geom)wetlands_area = 0for feat in water_lyr: intersect = feat.geometry().Intersection(nola_geom) # 求交 wetlands_area += intersect.GetArea()pcnt = wetlands_area / nola_geom.GetArea()print('{:.1%} of New Orleans is wetland'.format(pcnt))water_lyr.SetSpatialFilter(None)water_lyr.SetAttributeFilter("Feature != 'Lake'")memory_driver = ogr.GetDriverByName('Memory')temp_ds = memory_driver.CreateDataSource('temp')temp_lyr = temp_ds.CreateLayer('temp')nola_lyr.Intersection(water_lyr, temp_lyr)sql = 'SELECT SUM(OGR_GEOM_AREA) AS area FROM temp'lyr = temp_ds.ExecuteSQL(sql)pcnt = lyr.GetFeature(0).GetField('area') / nola_geom.GetArea()print('{:.1%} of New Orleans is wetland'.format(pcnt))28.7% of New Orleans is wetland2、鄰近分析（確定要素間的距離）wkz少兒編程網-Scratch_Python_教程_免費兒童編程學習平臺
??OGR包含兩個鄰近分析工具：測量幾何要素的距離；創建緩沖區。wkz少兒編程網-Scratch_Python_教程_免費兒童編程學習平臺
??1.確定美國有多少城市位于火山10英里（1英里=1609.3米）的范圍之內。確定火山附近城市數量的存在問題的方法：wkz少兒編程網-Scratch_Python_教程_免費兒童編程學習平臺
from osgeo import ogr shp_ds = ogr.Open(r'E:Google chromeDownloadgis with pythonosgeopy dataUS')volcano_lyr = shp_ds.GetLayer('us_volcanos_albers')cities_lyr = shp_ds.GetLayer('cities_albers')# 在內存中創建一個臨時層來存儲緩沖區memory_driver = ogr.GetDriverByName('memory')memory_ds = memory_driver.CreateDataSource('temp')buff_lyr = memory_ds.CreateLayer('buffer')buff_feat = ogr.Feature(buff_lyr.GetLayerDefn())# 緩緩沖每一個火山點，將結果添加到緩沖圖層中for volcano_feat in volcano_lyr: buff_geom = volcano_feat.geometry().Buffer(16000) tmp = buff_feat.SetGeometry(buff_geom) tmp = buff_lyr.CreateFeature(buff_feat)# 將城市圖層與火山緩沖區圖層相交result_lyr = memory_ds.CreateLayer('result')buff_lyr.Intersection(cities_lyr, result_lyr)print('Cities: {}'.format(result_lyr.GetFeatureCount())) Cities: 83??2.一個更好地確定火山附近城市數量方法：wkz少兒編程網-Scratch_Python_教程_免費兒童編程學習平臺
from osgeo import ogr shp_ds = ogr.Open(r'E:Google chromeDownloadgis with pythonosgeopy dataUS')volcano_lyr = shp_ds.GetLayer('us_volcanos_albers')cities_lyr = shp_ds.GetLayer('cities_albers')# 將緩沖區添加到一個復合多邊形，而不是一個臨時圖層multipoly = ogr.Geometry(ogr.wkbMultiPolygon)for volcano_feat in volcano_lyr: buff_geom = volcano_feat.geometry().Buffer(16000) multipoly.AddGeometry(buff_geom)# 將所有的緩沖區聯合在一起得到一個可以使用的多邊形作為空間過濾器cities_lyr.SetSpatialFilter(multipoly.UnionCascaded())print('Cities: {}'.format(cities_lyr.GetFeatureCount()))Cities: 78注：UnionCascaded()：有效地將所有的多邊形合并成一個復合多邊形
??第一個例子中，每當城市位于火山緩沖區內，就會復制到輸出結果中。說明一個城市位于多個16000米緩沖區內，將被列入不止一次。wkz少兒編程網-Scratch_Python_教程_免費兒童編程學習平臺
??3.計算特定的城市與火山的距離：wkz少兒編程網-Scratch_Python_教程_免費兒童編程學習平臺
import osfrom osgeo import ogrfrom ospybook.vectorplotter import VectorPlotter data_dir = r'E:Google chromeDownloadgis with pythonosgeopy data'shp_ds = ogr.Open(os.path.join(data_dir, 'US'))volcano_lyr = shp_ds.GetLayer('us_volcanos_albers')cities_lyr = shp_ds.GetLayer('cities_albers')# 西雅圖到雷尼爾山的距離volcano_lyr.SetAttributeFilter("NAME = 'Rainier'")feat = volcano_lyr.GetNextFeature()rainier = feat.geometry().Clone()cities_lyr.SetSpatialFilter(None)cities_lyr.SetAttributeFilter("NAME = 'Seattle'")feat = cities_lyr.GetNextFeature()seattle = feat.geometry().Clone()meters = round(rainier.Distance(seattle))miles = meters / 1600print('{} meters ({} miles)'.format(meters, miles))92656 meters (57.91 miles)??3. 用2.5D幾何對象，表示兩點之間的距離：wkz少兒編程網-Scratch_Python_教程_免費兒童編程學習平臺
# 2Dpt1_2d = ogr.Geometry(ogr.wkbPoint)pt1_2d.AddPoint(15, 15)pt2_2d = ogr.Geometry(ogr.wkbPoint)pt2_2d.AddPoint(15, 19)print(pt1_2d.Distance(pt2_2d))4.0# 2.5Dpt1_25d = ogr.Geometry(ogr.wkbPoint25D)pt1_25d.AddPoint(15, 15, 0)pt2_25d = ogr.Geometry(ogr.wkbPoint25D)pt2_25d.AddPoint(15, 19, 3)print(pt1_25d.Distance(pt2_25d))4.0??將高程Z值考慮進去，真正的距離是5。wkz少兒編程網-Scratch_Python_教程_免費兒童編程學習平臺
# 用2D計算面積ring = ogr.Geometry(ogr.wkbLinearRing)ring.AddPoint(10, 10)ring.AddPoint(10, 20)ring.AddPoint(20, 20)ring.AddPoint(20, 10)poly_2d = ogr.Geometry(ogr.wkbPolygon)poly_2d.AddGeometry(ring)poly_2d.CloseRings()print(poly_2d.GetArea())100.0# 用2.5D計算面積ring = ogr.Geometry(ogr.wkbLinearRing)ring.AddPoint(10, 10, 0)ring.AddPoint(10, 20, 0)ring.AddPoint(20, 20, 10)ring.AddPoint(20, 10, 10)poly_25d = ogr.Geometry(ogr.wkbPolygon25D)poly_25d.AddGeometry(ring)poly_25d.CloseRings()print(poly_25d.GetArea())100.0??2.5D的面積實際上是141。wkz少兒編程網-Scratch_Python_教程_免費兒童編程學習平臺
# 疊加操作同樣忽略了高程值Zprint(poly_2d.Contains(pt1_2d))print(poly_25d.Contains(pt1_2d))True True相關免費學習推薦：python教程（視頻）wkz少兒編程網-Scratch_Python_教程_免費兒童編程學習平臺
以上就是Python地理數據處理之分析使用GR進行矢量的詳細內容，更多請關注少兒編程網其它相關文章！wkz少兒編程網-Scratch_Python_教程_免費兒童編程學習平臺
- 上一篇
看看Python 高手都寫不出來的幾個錯誤
簡介對于剛入門的Pythonista在學習過程中運行代碼是或多或少會遇到一些錯誤，剛開始可能看起來比較費勁。隨著代碼量的積累，熟能生巧當遇到一些運行時錯誤時能夠很快的定位問題原題。下面整理了一些常見的17個錯誤，等你寫出的代碼不怎么出現這些錯誤的時候，你的Python功力就上一個臺階了。換一句話說，當你
- 下一篇
pycharm怎么設置背景顏色
簡介pycharm設置背景顏色的方法：1、打開pycharm，點擊【file】、【settings】；2、找到【editor】，選擇【colorscheme】；3、根據需要選擇合適的背景顏色即可。本文操作環境：windows10系統、pycharmcommunity2020.3、宏基S40-51電腦。具
相關資訊
