xiaoairen322的個人博客分享 http://www.ueservicedoffices.com/u/xiaoairen322

博文

GDAL空間數據轉換庫中的一些應用(ogr2ogr)

已有 760 次閱讀 2019-5-14 22:12 |個人分類:GMT|系統分類:科研筆記

  如果說小米是手機界的萬能之王,那么GDAL就是地學界的萬源之王


    GDAL(Geospatial Data Abstraction Library)是一個在X/MIT許可協議下的開源柵格空間數據轉換庫。它利用抽象數據模型來表達所支持的各種文件格式。它還有一系列命令行工具來進行數據轉換和處理。

OGR是GDAL項目的一個分支,功能與GDAL類似,只不過它提供對矢量數據的支持。

有很多著名的GIS類產品都使用了GDAL/OGR庫,包括ESRI的ARCGIS 9.3,Google Earth和跨平臺的GRASS GIS系統。利用GDAL/OGR庫,可以使基于Linux的地理空間數據管理系統提供對矢量和柵格文件數據的支持。

   最近由于想繪制一個區域的地形圖及行政區劃圖,現有的數據也是通過https://gadm.org/download_world.html和https://gadm.org/download_country_v3.html  還有基礎地理信息中心下發的中國矢量(shp)格式的圖,所以要用gmt繪制圖形必須得轉化為gmt的格所以gmt中文社區提供一種方法 GDAL (ogr2ogr)

  安裝GDAL(LIUNX )

  sudo apt-get install gdal-bin python-gdal

  安裝畢  測試 

  ogr2ogr --version

顯示版本即可,說明安裝成功 

  后面我們來說一下 我們轉換為gmt所能識別的格式

下面以下載所需要的shp 和 gpkg的數據為例 來講解ogr2ogr的用法 

Usage: ogr2ogr [--help-general] [-skipfailures] [-append] [-update]

               [-select field_list] [-where restricted_where]

               [-progress] [-sql <sql statement>] [-dialect dialect]

               [-preserve_fid] [-fid FID]

               [-spat xmin ymin xmax ymax] [-geomfield field]

               [-a_srs srs_def] [-t_srs srs_def] [-s_srs srs_def]

               [-f format_name] [-overwrite] [[-dsco NAME=VALUE] ...]

               dst_datasource_name src_datasource_name

               [-lco NAME=VALUE] [-nln name] [-nlt type] [-dim 2|3|layer_dim] [layer [layer ...]]


Advanced options :

               [-gt n]

               [-clipsrc [xmin ymin xmax ymax]|WKT|datasource|spat_extent]

               [-clipsrcsql sql_statement] [-clipsrclayer layer]

               [-clipsrcwhere expression]

               [-clipdst [xmin ymin xmax ymax]|WKT|datasource]

               [-clipdstsql sql_statement] [-clipdstlayer layer]

               [-clipdstwhere expression]

               [-wrapdateline][-datelineoffset val]

               [[-simplify tolerance] | [-segmentize max_dist]]

               [-addfields]

               [-relaxedFieldNameMatch]

               [-fieldTypeToString All|(type1[,type2]*)] [-unsetFieldWidth]

               [-fieldmap identity | index1[,index2]*]

               [-splitlistfields] [-maxsubfields val]

               [-explodecollections] [-zfield field_name]

               [-gcp pixel line easting northing [elevation]]* [-order n | -tps]


以gpkg為例  轉換為gmt 

 ogr2ogr -f GMT gadm3_CHN.gmt gadm36_CHN.gpkg


解釋一下 

 #ogr2ogr  為GDAL的一個轉化模塊

 #-f     為強制轉換選項

 #gadm36_CHN.gmt 為轉換后的文件

 #gadm36_CHN。gpkg 為待轉換的文件


轉換后的數據格式如下



shp格式n也一樣  一定要注意文件的順序 否則會報錯 ,有時候靈感來了 幾天未解決的幾秒鐘就解決了。很神奇


感謝gmt中文社區 



http://www.ueservicedoffices.com/blog-3387053-1179036.html

上一篇:關于GNSS工程控制網投影變形的一些事
下一篇:關于 envi 與sarscape的安裝一些

0

該博文允許注冊用戶評論 請點擊登錄 評論 (0 個評論)

數據加載中...
掃一掃,分享此博文

Archiver|手機版|科學網 ( 京ICP備14006957 )

GMT+8, 2019-8-19 12:29

Powered by ScienceNet.cn

Copyright © 2007- 中國科學報社

返回頂部
时时彩平台 乐山市  嘉黎县  时尚  巴东县  榆林市  安新县   鹤岗市  中方县  凯里市  灵石县   鄂伦春自治旗  竹溪县  阳西县  太康县  竹山县  宿迁市   长海县  天柱县  繁昌县  长阳