Python地圖柵格化實例

引言

shapefile是GIS中的一種非常重要的數(shù)據(jù)類型，由ESRI開發(fā)的空間數(shù)據(jù)開放格式，目前該數(shù)據(jù)格式已經(jīng)成為了GIS領域的開放標準。目前絕大多數(shù)開源以及收費的GIS軟件都支持該數(shù)據(jù)類型。事實上，shapefile文件指的一種文件存儲的方法，實際上該種文件是由多個文件組成的。組成一個shapefile有三種文件必不可少, '.shp','.shx','.dbf'文件。

成都創(chuàng)新互聯(lián)是網(wǎng)站建設專家,致力于互聯(lián)網(wǎng)品牌建設與網(wǎng)絡營銷，專業(yè)領域包括做網(wǎng)站、成都網(wǎng)站建設、電商網(wǎng)站制作開發(fā)、小程序開發(fā)、微信營銷、系統(tǒng)平臺開發(fā)，與其他網(wǎng)站設計及系統(tǒng)開發(fā)公司不同，我們的整合解決方案結合了恒基網(wǎng)絡品牌建設經(jīng)驗和互聯(lián)網(wǎng)整合營銷的理念，并將策略和執(zhí)行緊密結合，且不斷評估并優(yōu)化我們的方案，為客戶提供全方位的互聯(lián)網(wǎng)品牌整合方案！

geopandas對shapefile提供了很好的讀取與寫出支持。geopandas庫允許對幾何類型進行空間操作，其中DataFrame結構相當于GIS數(shù)據(jù)中的一張屬性表，使得可以直接操作矢量數(shù)據(jù)屬性表，使得python中操作地理數(shù)據(jù)更加方便。本實例通過geopandas實現(xiàn)對地理數(shù)據(jù)的操作。

開發(fā)準備

由于geopandas庫的安裝需要一些前提庫，因此需要先安裝一些庫

pip install pipwin
pipwin install gdal
pipwin install fiona
pip install geopandas

實測以上方法可以成功在windows下安裝（注：如果在Anaconda下安裝geopandas更為方便）

數(shù)據(jù)準備

下載地址

該數(shù)據(jù)是一段GPS掃描數(shù)據(jù)，包含經(jīng)緯度。

代碼實例

環(huán)境引入

import geopandas as gp
import matplotlib.pyplot as plt
from shapely import geometry
import math

GPS數(shù)據(jù)處理

lake_original_path = 'data.txt'
lake_original_data = ''
lake_points = []
# 讀取文件
with open(lake_original_path) as f:
    lake_original_data = f.read()
# 處理經(jīng)緯度坐標 并以Point的形式添加到list中
for xy in lake_original_data.split(';'):
    x, _, y = xy.partition(',')
    x = float(x.strip()) / 100
    y = float(y.strip()) / 100
    lake_points.append(geometry.Point(y, x))

創(chuàng)建要素

# 創(chuàng)建線狀要素
lake_line = geometry.LineString(lake_points)
# crs指定坐標系
lake_ = gp.GeoSeries(lake_line, crs='EPSG:4326')
# 保存shp文件
lake_.to_file("boundary.shp", driver='ESRI Shapefile', encoding='utf-8')
# 記錄邊界條件 用于構建柵格
x_min, y_min, x_max, y_max = lake_line.bounds[:4]
# 繪圖
lake_.plot()
plt.show()

構建柵格

# 柵格大小
GRID_WIDTH = 0.009 * 2 / 100
grid_rows_num = int(math.ceil((y_max - y_min) / float(GRID_WIDTH)))
grid_columns_num = int(math.ceil((x_max - x_min) / float(GRID_WIDTH)))
grids = []
for r in range(grid_rows_num):
    for c in range(grid_columns_num):
        grid_4coords = []
        # 左上角
        x_lt = x_min + c * GRID_WIDTH
        y_lt = y_max - r * GRID_WIDTH
        # 右上角
        x_rt = x_lt + GRID_WIDTH
        y_rt = y_lt
        # 左下角
        x_lb = x_lt
        y_lb = y_lt - GRID_WIDTH
        # 右下角
        x_rb = x_rt
        y_rb = y_lb
        # 兩個三角形拼接一個柵格
        grid_4coords.append(geometry.Point(x_lt,y_lt))
        grid_4coords.append(geometry.Point(x_rt,y_rt))
        grid_4coords.append(geometry.Point(x_rb,y_rb))
        grid_4coords.append(geometry.Point(x_lb,y_lb))
        grid_4coords.append(geometry.Point(x_lt,y_lt))
        # 創(chuàng)建一個網(wǎng)格
        grids.append(geometry.LineString(grid_4coords))
grid_ = gp.GeoSeries(grids)
grid_.to_file('E:\\just\\海韻湖智能技術實驗場\\data\\grids.shp',driver='ESRI Shapefile', encoding='utf-8')
grid_.plot()
plt.show()

要素疊加

# 要素疊加
elements = [lake_line]
elements += grids
elements_ = gp.GeoSeries(elements)
elements_.to_file('elements.shp', driver='ESRI Shapefile', encoding='utf-8')
elements_.plot()
plt.show()

參考鏈接

python-geopandas讀取、創(chuàng)建shapefile文件、geopandas學習教程

GeoPandas官方文檔

基于arcpy實現(xiàn)導出區(qū)域內網(wǎng)格中心坐標功能

gepandas安裝問題

當前名稱：Python地圖柵格化實例
文章來源：http://redsoil1982.com.cn/article48/dsogcep.html

成都網(wǎng)站建設公司_創(chuàng)新互聯(lián)，為您提供企業(yè)網(wǎng)站制作、外貿網(wǎng)站建設、網(wǎng)站收錄、響應式網(wǎng)站、網(wǎng)站策劃、微信小程序

廣告

聲明：本網(wǎng)站發(fā)布的內容（圖片、視頻和文字）以用戶投稿、用戶轉載內容為主，如果涉及侵權請盡快告知，我們將會在第一時間刪除。文章觀點不代表本網(wǎng)站立場，如需處理請聯(lián)系客服。電話：028-86922220；郵箱：631063699@qq.com。內容未經(jīng)允許不得轉載，或轉載時需注明來源： 創(chuàng)新互聯(lián)