GDAL WKT 学习 理解 使用

WKT

WKT表达几何类型、坐标信息、参考系统、投影定义字符串的起始表示格式通常是由类型关键字和一对方括号构成的，可以方便地被计算机程序解析和生成，也可以方便地被人类读写和理解。具体来说：

• 几何类型通常以一个大写字母作为关键字，放在方括号内，如“POINT”、“LINESTRING”、“POLYGON”等表示点、线、面。
• 参考系统通常以关键字“PROJCS”（投影坐标系统）或“GEOGCS”（地理坐标系统）作为起始关键字，放在方括号内。其后紧跟着参考系的详细信息，包括基准面、椭球体、基准点等。

地理坐标系统

下面是一张tif的投影定义字符串，我们来逐个理解它

GEOGCS["NAD83",

DATUM["North_American_Datum_1983",SPHEROID["GRS 1980",6378137,298.257222101004,AUTHORITY["EPSG","7019"]],AUTHORITY["EPSG","6269"]],

PRIMEM["Greenwich",0],

UNIT["degree",0.0174532925199433,AUTHORITY["EPSG","9122"]],

AXIS["Latitude",NORTH],

AXIS["Longitude",EAST],

AUTHORITY["EPSG","4269"]]

• GEOGCS：表示这是一个地理坐标系统。
• "NAD83"：坐标系统的名称为NAD83。
• DATUM：定义了该坐标系统使用的基准面和基准点。这里使用的基准面是北美1983大地基准（North_American_Datum_1983），而基准点则没有给出具体的描述。
• SPHEROID：定义了椭球体的参数，包括椭球体的长半轴、扁率等。这里使用的是GRS 1980椭球体，其长半轴为6378137米，扁率为298.257222101004。AUTHORITY["EPSG","7019"]指定了此椭球体的EPSG代码为7019。
• AUTHORITY：指定了此坐标系统的EPSG代码，其中“6269”是指定的北美基准面 1983的代码。
• PRIMEM：定义了使用的基准经线。这里使用的是格林威治子午线，即经度为0度的经线。
• UNIT：定义了坐标系统的单位。这里使用的是度（degree），其换算系数为0.0174532925199433弧度。AUTHORITY["EPSG","9122"]指定了此单位的EPSG代码为9122。
• AXIS：定义了坐标轴的方向。这里使用的是北纬和东经方向。
• AUTHORITY：最后一个AUTHORITY指定了该坐标系统的EPSG代码为4269。

如何获得坐标系信息

proj提供了工具projinfo来查询某个坐标系的WKT和PROJ格式：http://t.csdn.cn/WZpSM

EPSG网站提供了一个全球性的地理空间参考系统数据库，包括椭球体、投影、地理坐标系、转换等内容：https://epsg.org/home.html

点

用 POINT 关键字表示，后跟坐标序列，例如：POINT (30 10)。

线

用 LINESTRING 关键字表示，后跟坐标序列，例如：LINESTRING (30 10, 10 30, 40 40)。

面

用 POLYGON 关键字表示，后跟坐标序列，表示多边形的外环。如果有多个内环，每个内环都用一对圆括号括起来。

例如：POLYGON ((30 10, 40 40, 20 40, 10 20, 30 10), (20 30, 35 35, 30 20, 20 30))。

除了这些基本几何类型，WKT 还可以表示多点、多线、多面、集合等更复杂的几何类型 这里就不多赘述了