av一区二区在线观看_亚洲男人的天堂网站_日韩亚洲视频_在线成人免费_欧美日韩精品免费观看视频_久草视

您的位置:首頁技術文章
文章詳情頁

python 使用GDAL實現(xiàn)柵格tif轉矢量shp的方式小結

瀏覽:2日期:2022-07-28 13:05:19
前言

目前有一張tif格式的柵格影像,需要在web地圖上進行展示,使用動態(tài)切片WMS的方式,渲染速度比較慢,而且大的時候會出現(xiàn)模糊的問題。并且后面需要做多期影像的切換,渲染與加載效率也值得關注。

計劃是使用柵格轉矢量的方式,將柵格數(shù)據(jù)轉為矢量shp文件,然后進行矢量切片,使用Mapbox進行前端動態(tài)渲染。在網(wǎng)上查詢了很多資料,有人說使用d3-contour在node.js中生成或者使用rasterio在python中進行轉換,整體過程都比較麻煩,很不易實現(xiàn)。最終選定了使用GDAL進行柵格轉矢量的方法,代碼比較簡單。原始tif影像(12.8MB)如下:

python 使用GDAL實現(xiàn)柵格tif轉矢量shp的方式小結

核心函數(shù)

GDAL中柵格轉矢量的函數(shù)主要是以下兩個,二者的參數(shù)沒有任何區(qū)別,只是功能有區(qū)別:

FPolygonize(*args, **kwargs)

FPolygonize(Band srcBand, Band maskBand, Layer outLayer, int iPixValField, char options=None, GDALProgressFunc callback=0, void * callback_data=None) -> int

將每個像元轉成一個矩形。

Polygonize(*args, **kwargs) **

Polygonize(Band srcBand, Band maskBand, Layer outLayer, int iPixValField, char ** options=None, GDALProgressFunc callback=0, void * callback_data=None) -> int

將每個像元轉成一個矩形,然后將相似的像元進行合并。

轉換代碼

from osgeo import gdal, ogr, osrimport osimport datetimeimport numpy as nppath = 'Z_NAFP20210727.tif'if __name__ == ’__main__’: start_time = datetime.datetime.now() inraster = gdal.Open(path) # 讀取路徑中的柵格數(shù)據(jù) inband = inraster.GetRasterBand(1) # 這個波段就是最后想要轉為矢量的波段,如果是單波段數(shù)據(jù)的話那就都是1 prj = osr.SpatialReference() prj.ImportFromWkt(inraster.GetProjection()) # 讀取柵格數(shù)據(jù)的投影信息,用來為后面生成的矢量做準備 outshp = path[:-4] + '.shp' # 給后面生成的矢量準備一個輸出文件名,這里就是把原柵格的文件名后綴名改成shp了 drv = ogr.GetDriverByName('ESRI Shapefile') if os.path.exists(outshp): # 若文件已經存在,則刪除它繼續(xù)重新做一遍drv.DeleteDataSource(outshp) Polygon = drv.CreateDataSource(outshp) # 創(chuàng)建一個目標文件 Poly_layer = Polygon.CreateLayer(path[:-4], srs=prj, geom_type=ogr.wkbMultiPolygon) # 對shp文件創(chuàng)建一個圖層,定義為多個面類 newField = ogr.FieldDefn(’value’, ogr.OFTReal) # 給目標shp文件添加一個字段,用來存儲原始柵格的pixel value,浮點型, Poly_layer.CreateField(newField) gdal.Polygonize(inband, None, Poly_layer, 0) # 核心函數(shù),執(zhí)行的就是柵格轉矢量操作 # gdal.FPolygonize(inband, None, Poly_layer, 0) # 只轉矩形,不合并 Polygon.SyncToDisk() Polygon = None end_time = datetime.datetime.now() print('Succeeded at', end_time) print('Elapsed Time:', end_time - start_time) # 輸出程序運行所需時間轉換效果 使用FPolygonize

轉換之后的矢量數(shù)據(jù)有270MB,非常大,打開非常卡

python 使用GDAL實現(xiàn)柵格tif轉矢量shp的方式小結

使用Polygonize

合并之后的矢量數(shù)據(jù)有48MB,相對第一種方法數(shù)據(jù)量大大減少

python 使用GDAL實現(xiàn)柵格tif轉矢量shp的方式小結

到此這篇關于python 使用GDAL實現(xiàn)柵格tif轉矢量shp的文章就介紹到這了,更多相關python柵格tif轉矢量shp內容請搜索好吧啦網(wǎng)以前的文章或繼續(xù)瀏覽下面的相關文章希望大家以后多多支持好吧啦網(wǎng)!

標簽: Python 編程
相關文章:
主站蜘蛛池模板: 99re6热在线精品视频播放 | 99精品欧美一区二区三区 | 特黄小视频| 夜夜操天天干 | 国产成人精品一区二区三区在线 | 欧美日本在线观看 | 亚洲成人中文字幕 | 免费一区二区三区 | 亚洲精品在线免费观看视频 | 国产成人网 | 中文字幕日韩一区 | 蜜臀久久99精品久久久久久宅男 | 国产一区二区在线免费观看 | 毛色毛片免费看 | 国产精品96久久久久久 | 色在线免费 | 日韩久久久久 | 亚洲欧美在线观看视频 | 欧美激情精品久久久久 | 久久精品久久久久久 | 国产一区二区三区久久久久久久久 | 亚洲国产精品网站 | 久久久久国色av免费观看性色 | 黄色片a级 | 成人国产在线视频 | 免费久草 | 欧美中文一区 | 亚洲免费人成在线视频观看 | 看片地址 | 欧美日韩亚洲国产 | 国产97在线视频 | 欧美日韩综合视频 | 午夜精品一区二区三区在线播放 | 一区二区三区四区在线 | 久久一区视频 | 免费久久网 | 日韩欧美国产一区二区三区 | 日韩成人免费视频 | 久久国产综合 | 久久久久久国产精品三区 | 久久99精品久久久久久国产越南 |