网络地图服务（WMS）详解
根据具体实例详细论述了GIS中网络地图服务（WMS）的内容。
阅读本文之前可参考前文：《地图服务器GeoServer的安装与配置》与《GeoServer发布地图服务（WMS、WFS）》。
1.概述
经过前文的介绍，相信我们对WMS/WFS服务已经有了一个非常直观的认识，最起码我们知道了地图服务的数据从何而来，又是如何发布出去的。更进一步的说，OGC给WMS/WFS分别定义了非常详尽的标准规范，其主要内容是通过HTTP（Hypertext Transfer Protocol，超文本传输协议 ）实现的各种操作。在这里，我们要介绍的WMS支持标准HTTP的Get和Post请求，不过通常情况下基于GET方式（将参数以键值对的形式放置在URL中）进行服务请求较为方便。
HTTP协议用于客户端和服务器之间的通信，为了满足各种各样的通信需求，HTTP协议规定了多种不同的请求方法，其中Get和Post是最基础的两种。Get用于请求已被URI识别的资源，Post则用来传输实体的主体，两者的功能比较相近，但其实有着不小的区别。不过在本文中，我们将会主要看到这个区别：Get请求是将参数拼接到URL上进行参数传递的，而POST是将参数写入到请求正文中传递的；因此，Get请求有URL长度的限制，通常用于查询；Post请求则没有长度限制，通常用于增添、删除和修改。
根据《OpenGIS Web Map Server Implementation Specification》的1.3.0版本，WMS主要支持如下几种操作，如下表1所示所示：
【表1 WMS支持的操作】
| 操作 | 必要性 | 描述 |
|---|---|---|
| GetCapabilities | 必要 | 查询服务的元数据，包括支持的操作、参数、以及图层列表等 |
| GetMap | 必要 | 检查询指定区域和内容的地图图像 |
| GetFeatureInfo | 可选 | 查询地图上像素位置的基础数据，包括几何和属性值 |
2.GetCapabilities
GetCapabilities操作的目的是获取服务的元数据，主要包括WMS服务器提供的操作、请求参数等信息。GetCapabilities请求的参数如下表2所示：
【表2 WMS GetCapabilities请求参数】
| 请求参数 | 必要性 | 描述 |
|---|---|---|
| VERSION=version | 可选 | 请求的版本：值是1.1.0，1.1.1，1.3.0三者之一。 |
| SERVICE=WMS | 必要 | 服务的类型。 |
| REQUEST=GetCapabilities | 必要 | 请求操作的名称。 |
| FORMAT=MIME_type | 可选 | 服务元数据的格式。 |
| UPDATESEQUENCE=string | 可选 | 用于缓存控制的序列号或字符串。 |
根据表2所示的参数，我们在浏览器中输入如下地址：
http://localhost:8080/geoserver/ows?service=WMS&version=1.3.0&request=GetCapabilities
浏览器响应的结果如下图所示。可以看到返回的结果其实是一个XML文件，这个文件中描述的就是WMS服务具体的元数据信息。
3.GetMap
GetMap操作用于请求服务器生成地图并返回，这个操作也是WMS服务最主要的功能。不难想象，为了保证正确的输出结果，我们必须配置足够多的参数，比如地理范围，空间参考、输出宽高以及格式等。具体的请求参数如下表3所示：
【表3 WMS GetMap请求参数】
| 请求参数 | 必要性 | 描述 |
|---|---|---|
| VERSION=1.3.0 | 必要 | 请求的版本。 |
| REQUEST=GetMap | 必要 | 请求操作的名称。 |
| LAYERS=layer_list | 必要 | 一个或多个地图图层列表，使用逗号分隔。 |
| STYLES=style_list | 必要 | 渲染样式列表，使用逗号分隔。 |
| CRS=namespace:identifier | 必要 | 坐标参考系统。 |
| BBOX=minx,miny,maxx,maxy | 必要 | 边界框包围盒（从左下到右上，使用CRS的单位）。 |
| WIDTH=output_width | 必要 | 地图图像的宽度。 |
| HEIGHT=output_height | 必要 | 地图图像的高度。 |
| FORMAT=output_format | 必要 | 地图的输出格式。 |
| TRANSPARENT=TRUE|FALSE | 可选 | 地图背景透明度（默认不透明）。 |
| BGCOLOR=color_value | 可选 | 背景颜色的十六进制红绿蓝颜色值（默认白色）。 |
| EXCEPTIONS=exception_format | 可选 | WMS报告异常的格式（默认为XML）。 |
| TIME=time | 可选 | 图层需要的时间值。 |
| ELEVATION=elevation | 可选 | 图层需要的高程值。 |
| Other sample dimension(s) | 可选 | 其他维度的值，备用。 |
这些参数具体的参数值可以参考GetCapabilities操作返回的XML文件。例如，我们填充好所有的必要参数，在浏览器中输入如下地址来访问前文中配置好的矢量数据集：
http://localhost:8080/geoserver/wms?
VERSION=1.3.0
&REQUEST=GetMap
&LAYERS=test:multipolygons
&STYLES
&CRS=EPSG:4326
&BBOX=38.8920794,-77.0495419,38.9096339,-77.0199063
&WIDTH=768
&HEIGHT=454
&FORMAT=image/png
这时，浏览器就会返回一个PNG格式的图像地图文件，如下图所示。GetMap响应的地图图像大小、格式等并不是固定的，具体取决于请求的参数。如果我们对矢量和栅格足够熟悉的话，我们能大概猜出WMS服务返回地图文件大概的处理算法流程：
- 由于数据源是矢量数据，所以我们关键是矢量栅格化的算法。
- 根据传递的参数坐标参考系统CRS、地理范围BBox、宽WIDTH以及高HEIGHT足够帮助我们确定一个栅格数据集。
- 通过矢量栅格化算法对栅格数据集填充像素值，这中间涉及到图像的重采样操作。
- 最后将栅格数据集的文件流传输到浏览器端。
如果我们想访问前文中配置好的栅格数据集，URL请求参数大致相同，如下所示：
http://localhost:8080/geoserver/wms?
VERSION=1.3.0
&REQUEST=GetMap
&LAYERS=test:japan_tmo_2011057_geo
&STYLES
&CRS=EPSG:4326
&BBOX=23.034667983398435,120.15380860839844,52.990722641601565,153.88916014160156
&WIDTH=768
&HEIGHT=682
&FORMAT=image/png
此时，浏览器也会返回一张PNG格式的图像地图，如下图所示：
有一点要注意的是在WMS另一个广泛使用的版本1.1.1中，边界框包围盒的参数BBOX的x、y坐标顺序是与1.3.0版本是相反的，WMS版本1.1.1的请求URL如下所示：
http://localhost:8080/geoserver/wms?
VERSION=1.1.1
&REQUEST=GetMap
&LAYERS=test:multipolygons
&STYLES
&CRS=EPSG:4326
&BBOX=-77.0495419,38.8920794,-77.0199063,38.9096339
&WIDTH=768
&HEIGHT=454
&FORMAT=image/png
这个问题与GDAL2.0升级3.0版本时地理坐标系相关参数x和y需要进行调换有点类似，经过查阅相关资料才知道，原因是要保证x和y的顺序要与坐标参考系定义的轴顺序保持一致。但是在笔者看来，这种设计多少有点破坏向下兼容的原则。
4.GetFeatureInfo
GetFeatureInfo是可选操作，用于查询地图上给定位置的空间和属性信息。这个操作有点类似于鼠标在数字地图上拾取一点，高亮显示该点位置的要素并显示该要素关联的属性。不难理解，由于WMS使用的HTTP协议是无状态的，所以GetFeatureInfo请求参数包含了GetMap请求的大部分参数（除VERSION和REQUEST之外的所有参数）。毕竟，要获取地图某一点的要素信息，就先得把地图给获取了。GetFeatureInfo操作的标准请求参数如下表4所示：
【表4 WMS GetFeatureInfo请求参数】
| 请求参数 | 必要性 | 描述 |
|---|---|---|
| VERSION=1.3.0 | 必要 | 请求的版本。 |
| REQUEST=GetFeatureInfo | 必要 | 请求操作的名称。 |
| map request part | 必要 | GetMap请求的大部分参数。 |
| QUERY_LAYERS=layer_list | 必要 | 要查询的层级列表，以逗号分隔。 |
| INFO_FORMAT=output_format | 必要 | 特征信息返回格式（MIME类型）。 |
| FEATURE_COUNT=number | 可选 | 返回的特征数目 (默认为1)。 |
| I=pixel_column | 必要 | x方向像素坐标。 |
| J=pixel_row | 必要 | y方向像素坐标。 |
| EXCEPTIONS=exception_format | 可选 | WMS 报告异常的格式（默认为XML）。 |
我们按照表4所示的参数来查询前文中配置好的基于矢量数据集的地图，具体在某个像素位置的要素信息，可以在浏览器中输入如下地址：
http://localhost:8080/geoserver/wms?
VERSION=1.3.0
&REQUEST=GetFeatureInfo
&LAYERS=test:multipolygons
&STYLES
&CRS=EPSG:4326
&BBOX=38.8920794,-77.0495419,38.9096339,-77.0199063
&WIDTH=768
&HEIGHT=454
&FORMAT=image/png
&QUERY_LAYERS=test:multipolygons
&INFO_FORMAT=application/json
&I=384
&J=227
此时，浏览器会返回一个JSON，描述了该像素位置的几何信息，其关联的属性信息，以及一些元数据信息。如下图所示：
同样的，如果使用类似的参数来查询前文中配置好的基于栅格数据集的地图，在浏览器中输入如下地址：
http://localhost:8080/geoserver/wms?
VERSION=1.3.0
&REQUEST=GetFeatureInfo
&LAYERS=test:japan_tmo_2011057_geo
&STYLES
&CRS=EPSG:4326
&BBOX=23.034667983398435,120.15380860839844,52.990722641601565,153.88916014160156
&WIDTH=768
&HEIGHT=682
&FORMAT=image/png
&QUERY_LAYERS=test:japan_tmo_2011057_geo
&INFO_FORMAT=application/json
&I=384
&J=341
此时，浏览器会返回一个JSON，不过内容就简单多了，返回的是该地图在指定像素位置的包含RGB三个波段的像素值。如下图所示。这其实也不难理解，如果把栅格数据当作是点要素的集合，那么点的像素值自然就是其关联的属性数据了。
