GDAL——命令使用专题——ogrinfo命令
前言
GDAL(Geospatial Data Abstraction Library)是一个在X/MIT许可协议下的开源栅格空间数据转换库。它利用抽象数据模型来表达所支持的各种文件格式。它还有一系列命令行工具来进行数据转换和处理。
OGR是GDAL项目的一个分支,功能与GDAL类似,只不过它提供对矢量数据的支持。
这一篇就讲讲常见的CDAL相关命令操作,使用好GDAL的相关命令,能够让你事半功倍。有人会说,现在都可视化操作了,还用啥命令。但是当你的机器性能比较低,或者说你在无界面华的操作系统下工作,那么命令就非常有用了,如:Linux等。
主要针对:ogrinfo命令、gdalinfo命令、gdalsrsinfo命令、gdallocationinfo命令、gdal_translate命令、gdaltransform命令 等。这一节讨论ogrinfo命令。
软件环境:在上次安装MapServer中,定义环境变量安装。详见《MapServer Configuring with IIS》。
GDAL通用、常规、多用途命令
在cmd中输入:ogrinfo --help-general 或者 gdalinfo --help-general
复制代码
Generic GDAL utility command options:
--version: report version of GDAL in use.
--license: report GDAL license info.
--formats: report all configured format drivers.
--format [format]: details of one format.
--optfile filename: expand an option file into the argument list.
--config key value: set system configuration option.
--debug [on/off/value]: set debug level.
--pause: wait for user input, time to attach debugger
--locale [locale]: install locale for debugging (i.e. en_US.UTF-8)
--help-general: report detailed help on general options.
复制代码
以上列表中的相关命令为通用命令。
通用命令的部分详细使用见ogrinfo命令中的相关部分。
ogrinfo命令
用途:列出OGR所支持的数据信息。
官网:https://gdal.org/ogrinfo.html
数据类型:OGR支持的数据类型
在cmd中输入:ogrinfo --formats
数据应用类型 数据类型 操作权限 说明
netCDF -raster,vector- (rw+s) Network Common Data Format
AmigoCloud -vector- (rw+) AmigoCloud
OCI -vector- (rw+) Oracle Spatial
PCIDSK -raster,vector- (rw+v) PCIDSK Database File
JP2OpenJPEG -raster,vector- (rwv) JPEG-2000 driver based on OpenJPEG library
PDF -raster,vector- (rw+vs) Geospatial PDF
MBTiles -raster,vector- (rw+v) MBTiles
DB2ODBC -raster,vector- (rw+) IBM DB2 Spatial Database
ESRI Shapefile -vector- (rw+v) ESRI Shapefile
MapInfo File -vector- (rw+v) MapInfo File
UK .NTF -vector- (rov) UK .NTF
OGR_SDTS -vector- (rov) SDTS
S57 -vector- (rw+v) IHO S-57 (ENC)
DGN -vector- (rw+v) Microstation DGN
OGR_VRT -vector- (rov) VRT - Virtual Datasource
REC -vector- (ro) EPIInfo .REC
Memory -vector- (rw+) Memory
BNA -vector- (rw+v) Atlas BNA
CSV -vector- (rw+v) Comma Separated Value (.csv)
NAS -vector- (rov) NAS-ALKIS
GML -vector- (rw+v) Geography Markup Language (GML)
GPX -vector- (rw+v) GPX
LIBKML -vector- (rw+v) Keyhole Markup Language (LIBKML)
KML -vector- (rw+v) Keyhole Markup Language (KML)
GeoJSON -vector- (rw+v) GeoJSON
ESRIJSON -vector- (rov) ESRIJSON
TopoJSON -vector- (rov) TopoJSON
Interlis 1 -vector- (rw+v) Interlis 1
Interlis 2 -vector- (rw+v) Interlis 2
OGR_GMT -vector- (rw+v) GMT ASCII Vectors (.gmt)
GPKG -raster,vector- (rw+vs) GeoPackage
SQLite -vector- (rw+v) SQLite / Spatialite
ODBC -vector- (rw+) ODBC
WAsP -vector- (rw+v) WAsP .map format
PGeo -vector- (ro) ESRI Personal GeoDatabase
MSSQLSpatial -vector- (rw+) Microsoft SQL Server Spatial Database
PostgreSQL -vector- (rw+) PostgreSQL/PostGIS
MySQL -vector- (rw+) MySQL
OpenFileGDB -vector- (rov) ESRI FileGDB
XPlane -vector- (rov) X-Plane/Flightgear aeronautical data
DXF -vector- (rw+v) AutoCAD DXF
CAD -raster,vector- (rovs) AutoCAD Driver
Geoconcept -vector- (rw+v) Geoconcept
GeoRSS -vector- (rw+v) GeoRSS
GPSTrackMaker -vector- (rw+v) GPSTrackMaker
VFK -vector- (ro) Czech Cadastral Exchange Data Format
PGDUMP -vector- (w+v) PostgreSQL SQL dump
OSM -vector- (rov) OpenStreetMap XML and PBF
GPSBabel -vector- (rw+) GPSBabel
SUA -vector- (rov) Tim Newport-Peace's Special Use Airspace Format
OpenAir -vector- (rov) OpenAir
OGR_PDS -vector- (rov) Planetary Data Systems TABLE
WFS -vector- (rov) OGC WFS (Web Feature Service)
WFS3 -vector- (ro) OGC WFS 3 client (Web Feature Service)
HTF -vector- (rov) Hydrographic Transfer Vector
AeronavFAA -vector- (rov) Aeronav FAA
Geomedia -vector- (ro) Geomedia .mdb
EDIGEO -vector- (rov) French EDIGEO exchange format
GFT -vector- (rw+) Google Fusion Tables
SVG -vector- (rov) Scalable Vector Graphics
CouchDB -vector- (rw+) CouchDB / GeoCouch
Cloudant -vector- (rw+) Cloudant / CouchDB
Idrisi -vector- (rov) Idrisi Vector (.vct)
ARCGEN -vector- (rov) Arc/Info Generate
SEGUKOOA -vector- (rov) SEG-P1 / UKOOA P1/90
SEGY -vector- (rov) SEG-Y
ODS -vector- (rw+v) Open Document/ LibreOffice / OpenOffice Spreadsheet
XLSX -vector- (rw+v) MS Office Open XML spreadsheet
ElasticSearch -vector- (rw+) Elastic Search
Walk -vector- (ro) Walk
Carto -vector- (rw+) Carto
SXF -vector- (rov) Storage and eXchange Format
Selafin -vector- (rw+v) Selafin
JML -vector- (rw+v) OpenJUMP JML
PLSCENES -raster,vector- (ro) Planet Labs Scenes API
CSW -vector- (ro) OGC CSW (Catalog - Service for the Web)
VDV -vector- (rw+v) VDV-451/VDV-452/INTREST Data Format
GMLAS -vector- (rwv) Geography Markup Language (GML) driven by application schemas
MVT -vector- (rw+v) Mapbox Vector Tiles
TIGER -vector- (rw+v) U.S. Census TIGER/Line
AVCBin -vector- (rov) Arc/Info Binary Coverage
AVCE00 -vector- (rov) Arc/Info E00 (ASCII) Coverage
HTTP -raster,vector- (ro) HTTP Fetching Wrapper
说明:
raster:表示栅格数据
vector:表示矢量数据
ro:表示只读('ro' is read-only driver)
rw:读、写以及复制(read or write (i.e. supports CreateCopy))
rw+:读、写以及更新(read, write and update (i.e. supports Create)read, write and update (i.e. supports Create))
v:虚拟读写(appended for formats supporting virtual IO (/vsimem, /vsigzip, /vsizip, etc))
s:数据子集(appended for formats supporting subdatasets)
用法:
复制代码
ogrinfo [--help-general] [-ro] [-q] [-where restricted_where|\@filename]
[-spat xmin ymin xmax ymax] [-geomfield field] [-fid fid]
[-sql statement|\@filename] [-dialect dialect] [-al] [-rl] [-so] [-fields={YES/NO}]
[-geom={YES/NO/SUMMARY/WKT/ISO_WKT}] [--formats] [[-oo NAME=VALUE] ...]
[-nomd] [-listmdd] [-mdd domain|`all`]*
[-nocount] [-noextent]
datasource_name [layer [layer ...]]
复制代码
ogrinfo --help-general
cmd输入:ogrinfo --help-general
用于获取CDAL相关命令的通用支持
复制代码
C:\WINDOWS\system32>ogrinfo --help-general #CDAL多用途命令,通用命令
Generic GDAL utility command options:
--version: report version of GDAL in use. #查看版本信息
--license: report GDAL license info. #查看许可信息
--formats: report all configured format drivers. #查看支持的驱动应用,详情可查看:https://gdal.org/ogr_formats.html
--format [format]: details of one format. #查看某一格式支持的详细信息 ,如:ogrinfo --format "ESRI Shapefile"。注意,在ogrinfo --formats列表中的文件名称,如果名称含空格,则必须用双引号查询 。
--optfile filename: expand an option file into the argument list. #执行操作文件中的命令,想见:--optfile
--config key value: set system configuration option. #全局变量配置
--debug [on/off/value]: set debug level. #设置代开,关闭DEBUG,或者DEBUG的级别
--pause: wait for user input, time to attach debugger #在进行调试的时候,等待用户输入
--locale [locale]: install locale for debugging (i.e. en_US.UTF-8) #调试本地化安装
--help-general: report detailed help on general options.
复制代码
ogrinfo --optfile操作
在你的shape文件所在文件夹下面新建“states_ugl.opt(名字可以任意取,后缀名也可以是txt)”文件。
states_ugl.opt文件中的内容为:-ro -al -so -rl
在cmd中输入:cd /d E:\SvnWorkspace\LY_WEB_GIS\branches\Documents\ms4w-mapserver-for-wimdows\release-1911-x64-gdal-2-3-3-mapserver-7-2-1\apps\Example1.1\data 定位到你的shape文件实际路径。
在cmd中输入:ogrinfo --optfile states_ugl.opt states_ugl.shp
通常这样做事为了执行一些比较连贯和复杂的操作。
--config
GDAL的运行时配置选项,与构建时配置脚本的选项不同。运行时配置选项适用于所有平台,并在运行时进行评估。它们可以通过编程、命令行开关或用户在环境中设置。
配置选项通常用于更改gdal和ogr驱动程序的默认行为,在某些情况下还用于更改gdal和ogr核心。它们本质上是用户可以设置的全局变量。
--config命令格式为: --config key value 来设置的。详见:https://trac.osgeo.org/gdal/wiki/ConfigOptions
下面以SHAPE_RESTORE_SHX来讲解如何使用。
首先打开你存放shape文件目录,把shx文件删除掉,如下图:
然后再cmd中输入:cd /d E:\SvnWorkspace\LY_WEB_GIS\branches\Documents\ms4w-mapserver-for-wimdows\release-1911-x64-gdal-2-3-3-mapserver-7-2-1\apps\Example1.1\data
在cmd中输入:ogrinfo --config SHAPE_RESTORE_SHX YES states_ugl.shp
文件创建成功。
ogrinfo --config SHAPE_RESTORE_SHX YES states_ugl.shp 表示:在打开shp文件时,如果没有shx文件,则尝试创建。
其他操作不一一讨论,明白格式以及在实际需求中会使用就好。
--debug
设置debug级别,值为on、off、0或1
如:ogrinfo db.sqlite --debug on
--pause
在进行调试的时候,等待用户输入。
--locale
调试本地化安装,如:en_US.UTF-8
以上是GDAL的一些通用命令使用方法,下面针对ogrinfo特有命令做相关介绍。
ogrinfo常用命令
-ro:read only,以只读形式打开文件
-al:all layer,获取所有layer的特征信息。如在cmd中输入:ogrinfo -ro -al states_ugl.shp,可看到下图:
其所有特征字段信息数据均会展现出来。
-rl:random layer,在GDAL2.2版本开始使用,表示启用随机读取模式,即按在数据集中找到的顺序迭代特性,而不是按每层迭代特性。对于像OSM、GMLAS等格式文件读取速度回快很多。
-so:summary only,仅显示摘要信息:显示部分字段类型信息,仅显示投影、架构、功能计数和范围等摘要信息。如在cmd中输入:ogrinfo -ro -al -so states_ugl.shp
-q:quite,显示主要信息,包括坐标系、层模式、范围和功能计数(但是个人觉得没这么多,只有:图层数量,名称,层模式,可能是我的shp文件问题)。在cmd中输入:ogrinfo -q states_ugl.shp
复制代码
#输入
ogrinfo -ro -q states_ugl.shp
#输出
1: states_ugl (Polygon)
复制代码
-where:类似sql语句查询,只对shape文件的属性进行匹配查询。如在cmd中输入:ogrinfo -ro -al -so -where CLASS='land' states_ugl.shp
或者通过导入文件的方式执行where查询(GDAL2.1版本以后才支持)
在数据源文件夹中创建一个where.txt文件,内容为CLASS='water',然后在cmd中运行:ogrinfo -ro -al -so -where @where.txt shates_ugl.shp
-sql:类似sql语句查询,同时支持单独指定sql文件(.txt、.sql都可以)。
如在cmd中输入:ogrinfo -ro -al -so -sql "SELECT CLASS,STATE FROM states_ugl" E:\SvnWorkspace\LY_WEB_GIS\branches\Documents\ms4w-mapserver-for-wimdows\release-1911-x64-gdal-2-3-3-mapserver-7-2-1\apps\tutorial\data
或者在cmd中输入:ogrinfo -ro -al -where @E:\SvnWorkspace\LY_WEB_GIS\branches\Documents\ms4w-mapserver-for-wimdows\release-1911-x64-gdal-2-3-3-mapserver-7-2-1\apps\Example1.1\data\where.txt -sql @E:\SvnWorkspace\LY_WEB_GIS\branches\Documents\ms4w-mapserver-for-wimdows\release-1911-x64-gdal-2-3-3-mapserver-7-2-1\apps\Example1.1\data\sql.txt E:\SvnWorkspace\LY_WEB_GIS\branches\Documents\ms4w-mapserver-for-wimdows\release-1911-x64-gdal-2-3-3-mapserver-7-2-1\apps\tutorial\data
这个语句有上面的基础很好理解。文件路径是绝对路径,除非你的cmd命令定位的数据源文件夹中存在sql过滤文件和where过滤文件。
-dialect:ogr版本sql语句。不大了解,标准sql语句也能使用,而且官网说性能地下。同时,sqlite支持的sql语句都能使用,所以使用通用sql语句就好,不做进一步的尝试与讲解。
-spat xmin ymin xmax ymax:根据数据范围过滤。
在cmd中输入:ogrinfo -ro -al -so states_ugl.shp 可以查看其数据范围为“ (-97.238976, 41.619778) - (-82.122902, 49.385620) ” ,上图有。
在cmd中输入:ogrinfo -ro -al -so -spat -90 42 -82 45 states_ugl.shp
-geomfield field:
-fid:OGR支持的一个特殊过滤字段,fid。可以在 -where 中使用,如: -where "fid in (200,201,202)"
在cmd中输入:ogrinfo -ro -al -where "fid in (200,201,202)" states_ugl.shp
-fields={YES/NO}:设置是否显示字段数据,如果是YES,显示字段数据信息,如果为NO,不显示。默认YES。
在cmd中输入:ogrinfo -ro -al -fields=NO states_ugl.shp
-geom={YES/NO/SUMMARY/WKT/ISO_WKT}:用于控制字段数据、几何图形的POLYGON数据显示以及输出(GDAL2.1版本的默认格式ISO_WKT)。
YES:默认值,全部显示,显示格式为OGC WKT格式输出。
NO:不显示POLYGON数据
如在cmd中输入:ogrinfo -ro -al -gemo=NO states_ugl.shp
SUMMARY:显示简要数据,但是多了点的个数
设置为WKT、ISO_WKT效果通YES。WKT处其实显示的是坐标信息。
-oo NAME=VALUE:数据集打开操作项。
有点类似 --config KEY VALUE。只是-00 NAME=VALUE 是 名称=值。表达方式不一样。
复制代码
#data.json的定义
{
"type": "FeatureCollection",
"features":
[
{
"type": "Feature",
"geometry": {
"type": "Point",
"coordinates": [ 2, 49 ]
},
"properties": {
"a_property": "foo",
"some_object": {
"a_property": 1,
"another_property": 2
}
}
}
]
}
复制代码
然后在cmd中输入:
ogrinfo data.json -al -oo FLATTEN_NESTED_ATTRIBUTES=yes
可以看到如下转换:
复制代码
INFO: Open of `data.json'
using driver `GeoJSON' successful.
Layer name: data
Geometry: Point
Feature Count: 1
Extent: (2.000000, 49.000000) - (2.000000, 49.000000)
Layer SRS WKT:
GEOGCS["WGS 84",
DATUM["WGS_1984",
SPHEROID["WGS 84",6378137,298.257223563,
AUTHORITY["EPSG","7030"]],
AUTHORITY["EPSG","6326"]],
PRIMEM["Greenwich",0,
AUTHORITY["EPSG","8901"]],
UNIT["degree",0.0174532925199433,
AUTHORITY["EPSG","9122"]],
AUTHORITY["EPSG","4326"]]
a_property: String (0.0)
some_object_a_property: Integer (0.0)
some_object_another_property: Integer (0.0)
OGRFeature(data):0
a_property (String) = foo
some_object_a_property (Integer) = 1
some_object_another_property (Integer) = 2
POINT (2 49)
复制代码
-nomd:阻止元数据打印。有些数据集可能包含许多元数据字符串。
-listmdd:列出所有的元数据。
-nomd 和 -listmdd 都是针对元数据的。
复制代码
#metadata 元数据
Metadata:
DBF_DATE_LAST_UPDATE=2002-03-07
#在states_ugl.shp 文件查看的时候,你们会发现有这了数据。
#-nomd不显示这条
#默认显示, -listmdd显示有所的。
复制代码
-mdd domain:从指定的域显示元数据。all可以显示所有域的元数据。
-nocount:没有Feature数量,即图形数量Feature Count: 204
-noextent:不输出空间范围。即没有 Extent: (-97.238976, 41.619778) - (-82.122902, 49.385620)
--formats:列出OGR支持的数据格式。上面讲过。
datasource_name:表示数据源、或者多个数据源文件所在目录、或者其他定义的全局虚拟名称。查看:https://gdal.org/ogr_formats.html 此网页,了解OGR支持的数据格式。
如:在cmd中知己输入(不定位到数据所在文件夹):ogrinfo -ro -al -so E:\SvnWorkspace\LY_WEB_GIS\branches\Documents\ms4w-mapserver-for-wimdows\release-1911-x64-gdal-2-3-3-mapserver-7-2-1\apps\tutorial\data
会罗列出当前文件下下面所有可以打开的数据文件的信息。
layer:获取请求图层数据信息。
如在cmd中输入:ogrinfo -ro -al -so E:\SvnWorkspace\LY_WEB_GIS\branches\Documents\ms4w-mapserver-for-wimdows\release-1911-x64-gdal-2-3-3-mapserver-7-2-1\apps\tutorial\data roads_ugl
则是显示图层名称为roads_ugl的数据信息
或者输入:ogrinfo -ro -al -so E:\SvnWorkspace\LY_WEB_GIS\branches\Documents\ms4w-mapserver-for-wimdows\release-1911-x64-gdal-2-3-3-mapserver-7-2-1\apps\tutorial\data roads_ugl states_ugl
其他命令官网
如果想学习而不等待博客更新,可以先上官网看看。
gdalinfo命令
官网:https://gdal.org/gdalinfo.html
dalsrs