利用DotSpatial发布WMS, WFS服务

我们遇到的几个给政府部门做的GIS系统，一般都只要面子，只要好看，领导高兴得不得了，点点这里点点那里，哟，这按钮一点还会转，领导开心得跟朵花似的。。。要是搞个各种分析什么的全堆上来，他就嫌烦了。。。这不说是由于使用的功能比较简单，花钱买个ArcGIS太浪费了，要用个开源系统做一个，一来给客户节省开销，二来我们自己也能赚一点。然后就在挑开源GIS，主要到现在也就会写点C#的代码，Grass, QGis, GeoTools, GeoServer什么的就都没考虑，就在DotSpatial和SharpMap里选了DotSpatial。我们的要求不高，只要能完成常用的地图加载查询什么的就好了。首先写的WMS。

WMS

WMS服务必须实现的有 GetCapabilites, GetMap，可选的功能是GetFeatureInfo。（说明一下，在ArcGIS发布的ＷＭＳ服务里，还有GetStyles功能，但是在ＯＧＣ的1.3.0版本的ＷＭＳ说明文档里没有提到GetStyles，也就是说GetStyles并不是OGC标准功能）这里分条介绍一下：

    1. GetCapabilites

      其功能是获取WMS服务的功能列表，也就是描述一下我发布的这个WMS能实现哪几个功能。

      GetCapabilites请求有不同的版本，常见的是1.1.1和1.3.0的，这里我写的是1.3.0版本的，下面只针对1.3.0写一下必备的参数：       

                               

      是否必需

      描述

      REQUEST=GetCapabilities

      必需

            

      SERVICE=WMS

      必需

      大写的WMS

      VERSION=1.3.0

      可选

      这里我写的只支持1.3.0的，所以不设置的话就默认返回1.3.0的格式了

      FORMAT=MINE_type

      可选

      返回类型，我这里只返回text/tml格式

      UPDATESEQUENCE=string

      可选

      Sequence number or string for cache control，我这里忽略了

           
      例：http://sampleserver1.arcgisonline.com/ArcGIS/services/Specialty/ESRI_StatesCitiesRivers_USA/MapServer/WMSServer?version=1.3.0&request=GetMap&CRS=CRS:84&bbox=-178.217598,18.924782,-66.969271,71.406235&width=760&height=360&layers=0&styles=default&format=image/png

      在添加了一个一般处理程序wms.ashx之后，首先从HttpContext里获取各种参数，利用DotSpatial.Control.Map新建一个Map对象，接下来就是根据下面的返回格式进行解析拼凑了：

      xml太长了，直接附地址：http://sampleserver1.arcgisonline.com/ArcGIS/services/Specialty/ESRI_StatesCitiesRivers_USA/MapServer/WMSServer?version=1.3.0&request=GetCapabilities&service=WMS

      大体代码如下：
      string wmsNamespaceURI = "http://www.opengis.net/wms";
XmlDocument doc = new XmlDocument();
XmlNode rootNode = doc.CreateNode(capabilities.CreateNode(XmlNodeType.Element, "WMS_Capabilities", wmsNamespaceURI);
rootNode.Attributes.Append(doc.CreateAttribute(...));
...
doc.AppendChild(rootNode);

context.Response.Clear();
context.Response.ContentType = "text/xml";
XmlWriter writer = XmlWriter.Create(context.Response.OutputStream);
capabilities.WriteTo(writer);
writer.Close();
context.Response.End();

       

    2. GetMap

      GetMap方法主要是根据请求和地图来获取图片。

      这里依旧根据1.3.0版本简单介绍一下请求参数:

      是否必需

      描述

      REQUEST=GetMap

      必需

      VERSION=1.3.0

      必需

      LAYERS=layer_list

      必需

      用逗号分隔多个图层，请求全部置为空

      STYLES=style_list

      必需

      用逗号分隔多个图层，与LAYERS排序相同

      CRS=namespace:identifier

      必需

      如 CRS=EPSG:2436

      BBOX=minx,miny,maxx,maxy

      必需

      要请求的图上坐标范围

      WIDTH=output_width

      必需

      请求返回的图像的宽度

      HEIGHT=output_height

      必需

      FORMAT=output_format

      必需

      返回图像格式 image/png

      TRANSPARENT=TRUE|FALSE

      可选

      地图是否透明 true

      BGCOLOR=color_value

      可选

      背景颜色 White

      EXCEPTIONS=exception_format

      可选

      抛出异常默认为xml格式

      与GetCapabilities一样，也是取到参数之后，用参数初始化一个DotSpatial.Control.Map对象，把Map.Size设置成WIDTH和HEIGHT大小，这样导出的地图就正好符合输出尺寸了，设置好LAYERS和STYLES，缩放到BBOX，然后利用Map提供的SnapShot()方法获得图像并插入到Response中。后边代码：

      var img = map.SnapShot();

ImageEncodeInfo info;

ImageCodecInfo [] infos = ImageCodecInfo.GetImageEncoders();
 for (int i = 0; i < infos.Length; i++)
{
    if (infos[i].MimeType == pMimeType)
    {
        info = infos[i];
        break;
    }
}

byte[] buffer;
 using (var ms = new MemoryStream())
{
       img.Save(ms,  info, null);
       img.Dispose();
       buffer = ms.ToArray();
}
context.Response.Clear();
context.Response.ContentType = pMimeType;
context.Response.OutputStream.Write(buffer, 0, buffer.Length);
context.Response.End();

    3. GetFeatureInfo

      根据请求参数获取要素。

      还是1.3.0版本的参数设置：

      是否必需

      描述信息

      REQUEST=GetFeatureInfo

      必需

      VERSION=1.3.0

      必需

      #regionGetMap的设置地图部分

      LAYERS=layer_list

      必需

      用逗号分隔多个图层，请求全部置为空

      STYLES=style_list

      必需

      用逗号分隔多个图层，与LAYERS排序相同

      CRS=namespace:identifier

      必需

      如 CRS=EPSG:2436

      BBOX=minx,miny,maxx,maxy

      必需

      地图的当前范围

      WIDTH=output_width

      必需

      地图的当前范围

      HEIGHT=output_height

      必需

      请求返回的图像的宽度

      …其它可选参数

      #endregion

      I=pixel_column

      必需

      请求点在图片上的位置Ｘ

      J=pixel_row

      必需

      QUERY_LAYERS=layer_list

      必需

      查询的图层，逗号分隔，必须包含在LAYERS中

      INFO_FORMAT=output_format

      必需

      要素信息返回的格式Capabilities里<Request><FeatureInfo><Format>中的格式

      FEATURE_COUNT=number

      可选

      最大返回数量

      EXCEPTIONS=exception_format

      可选


      例：http://sampleserver1.arcgisonline.com/arcgis/services/Specialty/ESRI_StatesCitiesRivers_USA/MapServer/WMSServer?version=1.3.0&request=GetFeatureInfo&layers=0&styles=default&CRS=EPSG:4326&bbox=-125.192865,11.2289864971264,
      -66.105824,62.5056715028736&width=1044&height=906&format=text/html&I=500&J=400&query_layers=0&INFO_FORMAT=text/html

      GetFeatureInfo也与上面无异，首先设置好Map，然后通过DotSpatial.Control.IMapFeatureLayer.Select(Extent)工具选出要查找的要素，最后把要素的属性拼接出来。

      除了用Extent选择以外，还可以将得到的List<IFeature>用CQL_FILTER进行过滤。

      把Context.Response.ContentType设置成要返回的格式，Context.Response.Write(string)写入取得的字符串。

WFS

WFS(Web Feature Service)服务必须实现的功能有GetCapabilities 、 DescribeFeatureType 和 GetFeature，可选的功能有GetGmlObject 和 Tansaction。这里说明的版本为1.1.0

这里都是根据索引号为OGC 04-094 Version: 1.1.0的Web Feature Service Implementation Specification 里的说明。

    1. GetCapabilities

      这是所有OWS(OGC Web Service)都必须实现的功能，功能都是返回该服务的服务能力等元数据。

      下面是1.1.0版本的参数说明

       

      是否必需

      		  

      REQUEST=GetCapabilities

      必需

      		  


      还是拼xml，和WMS的GetCapabilities差不多了。请求如：http://…..?Request=GetCapabilities 

      <?xml version="1.0" encoding="UTF-8"?>
<wfs:WFS_Capabilities xmlns:ows="http://www.opengis.net/ows" version="1.1.0" xmlns:Basemap="http://sritserver-pc/ArcGIS/services/Basemap/MapServer/WFSServer" xmlns:xlink="http://www.w3.org/1999/xlink" xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" xmlns:gml="http://www.opengis.net/gml" xmlns:ogc="http://www.opengis.net/ogc" xmlns:wfs="http://www.opengis.net/wfs" xsi:schemaLocation="http://www.opengis.net/gml http://schemas.opengis.net/gml/3.1.1/base/gml.xsd http://www.opengis.net/ogc http://schemas.opengis.net/filter/1.1.0/filter.xsd http://www.opengis.net/ows http://schemas.opengis.net/ows/1.0.0/owsAll.xsd http://www.opengis.net/wfs http://schemas.opengis.net/wfs/1.1.0/wfs.xsd">
    <ows:ServiceIdentification>……
    <ows:ServiceProvider>……
    <ows:OperationsMetadata>……
    <wfs:FeatureTypeList>……
    <ogc:Filter_Capabilities>……
</wfs:WFS_Capabilities>

    2. DescribeFeatureType

      请求参数：

      <xsd:element name="DescribeFeatureType" type="wfs:DescribeFeatureTypeType"/>
<xsd:complexType name="DescribeFeatureTypeType">
    <xsd:complexContent>
        <xsd:extension base="wfs:BaseRequestType">
            <xsd:sequence>
                <xsd:element name="TypeName" type="xsd:QName"minOccurs="0" maxOccurs="unbounded"/>
            </xsd:sequence>
            <xsd:attribute name="outputFormat"type="xsd:string" use="optional"default="text/xml; subtype=gml/3.1.1"/>
        </xsd:extension>
    </xsd:complexContent>
</xsd:complexType>
       

      是否必需

      		  

      TypeName

      不必需

      如果不带TypeName参数则返回所有图层

      outputFormat

      不必需

      该属性用于设置返回的要素类型的描述语言，默认为text/xml; subtype=gml/3.1.1，如果设置错了不识别会自动跳过该属性


      返回格式大致如下：(http://…..?Request=DescribeFeatureType&TypeName=SQ).具体的格式可以查看OGC标准文档说明或arcgis发布的服务的返回。 

      <?xml version="1.0" encoding="UTF-8"?>
<xs:schema xmlns:Basemap="http://sritserver-pc/ArcGIS/services/Basemap/MapServer/WFSServer" xmlns:gml="http://www.opengis.net/gml" xmlns:xs="http://www.w3.org/2001/XMLSchema" elementFormDefault="qualified" targetNamespace="http://sritserverpc/ArcGIS/services/Basemap/MapServer/WFSServer"><xs:import schemaLocation="http://schemas.opengis.net/gml/3.1.1/base/gml.xsd" namespace="http://www.opengis.net/gml"/>
    <xs:element type="Basemap:SQType" name="SQ" substitutionGroup="gml:_Feature"/>
    <xs:complexType name="SQType">……
</xs:schema>

    3. GetFeature


      常用的几个参数如下：

       

      是否必需

      		  

      Request=GetFeature

      必需

      		  

      featureID=SQ.F1026__6

      		  

      直接取得SQ图层里featureID（不是ObjectID）为F…__6的要素，以F开头，最后6是ObjectID，中间的几个字母不知道是怎么对应的，不过可以根据图层的各信息来编码，只要不重复就可以

      TypeName=SQ

      如果featureID未指定，则必需，多个用逗号隔开

      		  

      BBOX=xmin,ymin,xmax,ymax

      非必需

      		  

      outputFormat

      非必需

      默认值text/gml; subtype=gml/3.1.1，还有text/xml; subtype=gml/2.1.2等

      maxFeatures

      非必需

      		  

      featureVersion=ALL

      非必需

      对于支持版本的系统返回指定版本的要素

      resultType=Results

      非必需

      还有Hits可选值，Hits时不返回所有结果，只返回结果的一些描述，如数量、结构等

      …

      		    

      其它参数…

      		    

      依旧是根据请求的参数 和 Map 来查找需要返回的要素，然后转换成 GML 的格式返回，这里列一下我用到的几个格式，不一定全：

      对于Shape，最顶层节点都是<Shape>

      点。 
      <gml:Point> 
    <gml:pos>x y<gml:pos> 
</gml:Point>

       
      线，多线 
      <gml:multiCurve>
    <gml:CurveMember><!--每一段对应一个CurveMember-->
        <gml:LineString>
            <gml:posList>x1 y1 x2 y2 x3 y3……</gml:posList>
        </gml:LineString>
    </gml:CurveMember>
    <gml:CurveMember>…    </gml:CurveMember>
</gml:multiCurve>
      面，多面 
      <gml:MultiSurface>
    <gml:surfaceMember><!--每一部分面对应一个surfaceMember-->
        <gml:Polygon>
            <gml:exterior><!--只有一个外环-->
                <gml:LinearRing>
                    <gml:posList>x1 y1 x2 y2 x3 y3......</gml:posList>
                </gml:LinearRing>
            </gml:exterior>
            <gml:interior><!--每一个内环对应一个interior-->
                <gml:LinearRing>
                    <gml:posList>x1 y1 x2 y2 x3 y3......</gml:posList>
                </gml:LinearRing>
            </gml:interior>
            <gml:interior>......</gml:interior>
        </gml:Polygon>
    </gml:surfaceMember>
    <gml:surfaceMember>......</gml:surfaceMember>
<gml:MultiSurface>


转载请注明出处

posted @ 2013-09-03 16:21  toffrey  阅读(3174)  评论(0)  编辑  收藏  举报