利用DotSpatial发布WMS, WFS服务

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

WMS

WMS服务必须实现的有 GetCapabilites, GetMap,可选的功能是GetFeatureInfo。(说明一下,在ArcGIS发布的WMS服务里,还有GetStyles功能,但是在OGC的1.3.0版本的WMS说明文档里没有提到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

      必需

      请求点在图片上的位置X

      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 yuantf 阅读(...) 评论(...) 编辑 收藏