GeoServer 是 OpenGIS Web 服务器规范的 J2EE 实现的社区开源项目,利用 GeoServer 可以方便的发布地图数据,允许用户对特征数据进行更新、删除、插入操作,通过 GeoServer 可以比较容易的在用户之间迅速共享空间地理信息。本系列博文提供全面、完善的GeoServer部署解决方案,包括GeoServer环境搭建、地图数据处理、部署地图数据、发布地图服务等功能的详细介绍。文中内容来自本人工作中通过网络学习后总结而成,如有类同纯属巧合,同时欢迎广大网友前来交流。
系列目录导航:
GeoServer地图开发解决方案(一):环境搭建篇
GeoServer地图开发解决方案(二):地图数据处理篇
GeoServer地图开发解决方案(三):部署地图数据篇
GeoServer地图开发解决方案(四):发布Web地图服务(WMS)篇
GeoServer地图开发解决方案(五):基于Silverlight技术的地图客户端实现
我曾经写作过一篇关于微软Bing Maps的客户端实现的博文:《基于DeepZoom技术的Bing Maps客户端实现研究 》,详细介绍了如何使用Silverlight中的DeepZoom技术实现Bing Maps的客户端。本篇介绍的内容则为基于Web地图服务(Web Map Service,简称:WMS)的Silverlight地图客户端实现。
一、DeepZoom简介
DeepZoom技术以MultiScaleImage控件为核心,其内部有一个MultiScaleTileSource类型的源属性,主要用于设置MultiScaleImage控件所要呈现的数据源。基于Silverlight的Web GIS客户端实现也是通MultiScaleImage控件来实现,核心就在于通过MultiScaleTileSource属性针对不同的Web GIS地图瓦片数据(Image Tiles)提供商为MultiScaleImage控件实现一个数据源。因此本篇所需要做的工作就是针对WMS服务为MultiScaleImage控件实现一套加载数据源的算法。
二、WMS服务加载实现
实现WMS服务加载的算法其实非常简单,只需要了解WMS发布的方式、WMS地址的参数组成结构以及地图瓦片的投影原理就可以了,首先需要定义一个盒子对象作为访问WMS的边界参数对象。
public
class
BBox {
public
int
X {
get
;
set
; }
public
int
Y {
get
;
set
; }
public
int
Width {
get
;
set
; }
public
int
Height {
get
;
set
; }
public
BBox(
int
x,
int
y,
int
w,
int
h) {
this
.X
=
x;
this
.Y
=
y;
this
.Width
=
w;
this
.Height
=
h; } }
关于WMS服务加载的详细算法需要一些GIS理论基础才能够知道具体的实现原理,这里我就不逐一介绍,直接贴代码:
public
class
WMSTileSource : MultiScaleTileSource {
public
WMSTileSource() :
base
(
int
.MaxValue,
int
.MaxValue,
0x100
,
0x100
,
0
) { }
public
const
int
TILE_SIZE
=
256
;
///
<summary>
///
地球半径
///
</summary>
public
const
double
EARTH_RADIUS
=
6378137
;
///
<summary>
///
地球周长
///
</summary>
public
const
double
EARTH_CIRCUMFERENCE
=
EARTH_RADIUS
*
2
*
Math.PI;
public
const
double
HALF_EARTH_CIRCUMFERENCE
=
EARTH_CIRCUMFERENCE
/
2
;
///
<summary>
///
WMS服务地址
///
</summary>
private
const
string
TilePath
=
@"
http://localhost:8080/geoserver/wms?service=WMS&version=1.1.0&request=GetMap&layers=cq:CQ_County_region,cq:CQ_County_region_level&styles=&bbox={0},{1},{2},{3}&width=512&height=421&srs=EPSG:4326&&Format=image/png
"
;
public
string
GetQuadKey(
string
url) { var regex
=
new
Regex(
"
.*tiles/(.+)[.].*
"
); Match match
=
regex.Match(url);
return
match.Groups[
1
].ToString(); }
public
BBox QuadKeyToBBox(
string
quadKey,
int
x,
int
y,
int
zoomLevel) {
char
c
=
quadKey[
0
];
int
tileSize
=
2
<<
(
18
-
zoomLevel
-
1
);
if
(c
==
'
0
'
) { y
=
y
-
tileSize; }
else
if
(c
==
'
1
'
) { y
=
y
-
tileSize; x
=
x
+
tileSize; }
else
if
(c
==
'
3
'
) { x
=
x
+
tileSize; }
if
(quadKey.Length
>
1
) {
return
QuadKeyToBBox(quadKey.Substring(
1
), x, y, zoomLevel
+
1
); }
return
new
BBox(x, y, tileSize, tileSize); }
public
BBox QuadKeyToBBox(
string
quadKey) {
const
int
x
=
0
;
const
int
y
=
262144
;
return
QuadKeyToBBox(quadKey, x, y,
1
); }
public
double
XToLongitudeAtZoom(
int
x,
int
zoom) {
double
arc
=
EARTH_CIRCUMFERENCE
/
((
1
<<
zoom)
*
TILE_SIZE);
double
metersX
=
(x
*
arc)
-
HALF_EARTH_CIRCUMFERENCE;
double
result
=
RadToDeg(metersX
/
EARTH_RADIUS);
return
result; }
public
double
YToLatitudeAtZoom(
int
y,
int
zoom) {
double
arc
=
EARTH_CIRCUMFERENCE
/
((
1
<<
zoom)
*
TILE_SIZE);
double
metersY
=
HALF_EARTH_CIRCUMFERENCE
-
(y
*
arc);
double
a
=
Math.Exp(metersY
*
2
/
EARTH_RADIUS);
double
result
=
RadToDeg(Math.Asin((a
-
1
)
/
(a
+
1
)));
return
result; }
public
double
RadToDeg(
double
d) {
return
d
/
Math.PI
*
180.0
; }
private
static
string
TileXYToQuadKey(
int
tileX,
int
tileY,
int
levelOfDetail) { var quadKey
=
new
StringBuilder();
for
(
int
i
=
levelOfDetail; i
>
0
; i
--
) {
char
digit
=
'
0
'
;
int
mask
=
1
<<
(i
-
1
);
if
((tileX
&
mask)
!=
0
) { digit
++
; }
if
((tileY
&
mask)
!=
0
) { digit
++
; digit
++
; } quadKey.Append(digit); }
return
quadKey.ToString(); }
protected
override
void
GetTileLayers(
int
tileLevel,
int
tilePositionX,
int
tilePositionY, System.Collections.Generic.IList
<
object
>
tileImageLayerSources) {
int
zoom
=
tileLevel
-
8
;
if
(zoom
>
0
) {
string
quadKey
=
TileXYToQuadKey(tilePositionX, tilePositionY, zoom); BBox boundingBox
=
QuadKeyToBBox(quadKey);
double
lon
=
XToLongitudeAtZoom(boundingBox.X
*
TILE_SIZE,
18
);
double
lat
=
YToLatitudeAtZoom(boundingBox.Y
*
TILE_SIZE,
18
);
double
lon2
=
XToLongitudeAtZoom((boundingBox.X
+
boundingBox.Width)
*
TILE_SIZE,
18
);
double
lat2
=
YToLatitudeAtZoom((boundingBox.Y
-
boundingBox.Height)
*
TILE_SIZE,
18
);
string
wmsUrl
=
string
.Format(TilePath, lon, lat, lon2, lat2, TILE_SIZE); var veUri
=
new
Uri(wmsUrl); tileImageLayerSources.Add(veUri); } } }
前端通过一个按钮事件驱动触发加载WMS服务,按钮的XAML代码如下:
<
Button
Content
="WMS图层"
Height
="30"
Width
="80"
Name
="btnWms"
Click
="btnWms_Click"
/>
示例我就直接基于《基于DeepZoom技术的Bing Maps客户端实现研究 》一文中的示例扩展,对应的后台代码为如下代码块:
private
void
btnWms_Click(
object
sender, RoutedEventArgs e) { msi.Source
=
new
WMSTileSource(); }
版权说明
本文属原创文章,欢迎转载且注明文章出处,其版权归作者和博客园共有。为了保存作者的创作热情,请在转载后的明显位置标记本文出处。
作 者:Beniao
文章出处:http://beniao.cnblogs.com/ 或 http://www.cnblogs.com/