从天地图下载瓦片构建ArcGIS Server底图服务

  • Post author:
  • Post category:其他


从天地图下载瓦片构建ArcGIS Server底图服务

我们可以通过调用天地图发布的API来将天地图作为自己的底图,但是这需要我们去学习如何天地图的API,并且要将天地图集成到自己的WebGIS系统中则比较繁琐,并且在公司不能使用外网的情况下,调用天地图的API将变为不可行。今天我们通过下载天地图的瓦片地图,然后按照ArcGIS的瓦片组织规则将天地图的瓦片构建为自己服务器上的底图服务。

步骤如下:

1、从天地图下载瓦片

我们看到的天地图其实是由一个个瓦片拼接而成的,而这些瓦片存储在天地图的服务上,并且每一个瓦片都有自己的URI,天地图上由四种类型的瓦片,分别是:vec(矢量底图瓦片)、cva(矢量注记瓦片)、img(影像底图瓦片)、cia(影像注记瓦片),他们对应的URI分别是:

矢量底图:

http://t0.tianditu.com/DataServer?T=vec_c&x=0&y=0&l=1

矢量注记:

http://t0.tianditu.com/DataServer?T=cva_c&x=0&y=0&l=1

影像底图:

http://t0.tianditu.com/DataServer?T=img_c&x=0&y=0&l=1

影像注记:

http://t0.tianditu.com/DataServer?T=cia_c&x=0&y=0&l=1

其中每个URI包含四个参数,分别是:

T=vec_c ; x = ? ; y = ? ; l = ?;

T代表底图类型,x代表瓦片的列号(起始位0),y代表瓦片的行号(起始位0),l代表底图的级别(起始位1,最大为20)。

天地图的切片规则是这样的,l=1时,整幅地图(全球地图)被切为两片,如图(l=1):

当l=2即以后,每个瓦片将被切位4片,如图(l=2):

2、通过坐标值计算在某一级别的行列号

通常,我们需要某一范围的地图,所以我们只需要下载一个范围内的地图,那如何通过一个范围的坐标系获得该范围的所有瓦片的行列好呢?

假设,X、Y为一个点的坐标值,而x、y为这个点在级别l时所处瓦片的行列号,通过下面一个公式来计算(X、Y)在级别l所处瓦片的行列号(x、y):

仅限在中国范围内

x = ceiling((180 + X)/s*256)-1;

y = ceiling((90-Y)/s*256)-1;

其中s代表在级别l时的分辨率,分辨率如下表:(分辨率(s)和级别(l)一一对应)

    this._listResolution.Add(0.70312500015485435);l=1
    this._listResolution.Add(0.35156250007742718);l=2
    this._listResolution.Add(0.17578125003871359);l=3
    this._listResolution.Add(0.0878906250193568); l=4
    this._listResolution.Add(0.0439453125096784); l=5
    this._listResolution.Add(0.0219726562548392); l=6
    this._listResolution.Add(0.0109863281274196); l=7
    this._listResolution.Add(0.0054931640637098); l=8
    this._listResolution.Add(0.0027465820318549957); l=9
    this._listResolution.Add(0.0013732910159274978);l=10
    this._listResolution.Add(0.00068664549607834132); l=11
    this._listResolution.Add(0.00034332275992416907); l=12
    this._listResolution.Add(0.00017166136807812298); l=13
    this._listResolution.Add(8.5830684039061379E-05); l=14
    this._listResolution.Add(4.2915342019530649E-05); l=15
    this._listResolution.Add(2.1457682893727977E-05); l=16
    this._listResolution.Add(1.0728841446864E-05); l=17
    this._listResolution.Add(5.3644207234319882E-06); l=18
    this._listResolution.Add(2.6822103617159941E-06); l=19
    this._listResolution.Add(1.341105180858E-06); l=20

公式讲解:

180 + X 代表从左至右的经度距离,s代表地图分辨率,就是每一个像素表示的距离(在这里单位为经纬度),而天地图的每一个切片为256*256像素的,所以每一个切片表示的距离为256*s,(180 + X)/s*256 则表示(X、Y)点所在切片的列号,可能是一个小数,所以则通过ceiling函数获得该数向上舍入(1.2则转化成2),由于行列号都是从0开始计数,故最后再减去1。如此便得到在 级别为 l 时,(X、Y)点的所在切片行列号为(x、y)。

我们只需要通过确定两个点的坐标就能确定一个范围,分别是点(Xmin,Ymin)和(Xmax,Ymax),则通过两个点获得行列的最大值和最小值,然后通过循环就能得到所有的瓦片,然后通过.NET中system.NET命名空间中的Webclient类来下载每一个切片,并将每一个切片保存在本机上,代码如下:


WebClient


client

= new

WebClient

();

client.

Headers

.

Add

(“user-agent”, “Mozilla/5.0 (Windows NT 6.1; Win64; x64; rv:21.0) Gecko/20130109 Firefox/21.0”);


string


address

=

string

.

Format

(strTileUrl, nCol, nRow, nlevel);


string


str2

=

Convert

.

ToString

((

long

) nCol, 0x10).

ToUpper

();


string


str3

= “00000000” + str2;

str3 = “C” + str3.

Substring

(str2.

Length

);


string


fileName

=

string

.

Format

(@”{0}\{1}.png”, strSavePath, str3);

client.

DownloadFile

(address, fileName);

3、如何重新组织下载的切片

下载的切片是散乱的,必须通过有效的组织才能作为ArcGIS Server的缓冲切片,则需要遵循ArcGIS Server对瓦片的组织方式:

ArcGIS Server中的缓存是以文件夹结合文件名来组织和读取瓦片,如图:

我们可以看到,文件夹的第一级别为 L00、L01。。。。。,代表底图的级别;

文件夹的第二级别为R000000A1、R000000A2.。。。。。,代表一个级别中,所有的行,R后面的值就是上面计算出来的y值的范围,其中后面是通过16进制的编码;

文件夹的第三级为C0000033A、C0000033B。。。。。,代表一行中所有的列,C后面的值就是上面计算出来的x值的范围,也是16进制的编码;

例如:瓦片

http://t0.tianditu.com/DataServer?T=cia_c&x=567&y=139&l=11

它的行列号和级别分别为139、567、11,则它存储的位置为【L11】—【R0000008B】—【L00000237】。

4、发布ArcGIS Server服务

使用(Xmin,Ymin)和(Xmax,Ymax)两个点做一个矩形,然后将该矩形发布为地图服务,并将地图切片,这样在arcgisservercache文件夹下则生成这个地图服务的缓存文件和配置文件,这时候只要把上面下载好的天地图瓦片文件夹复制到相应的文件夹下,然后修改配置文件:

只需要修改conf.xml和L00、L01等级别文件相互映射。将比例尺、分辨率和地图级别对应即可。

ArcGIS Server 是如何读取瓦片并显示在相应的位置的那?

就是通过conf.xml中的地图分辨率配合图片的文件位置和编号计算出瓦片的位置,比如位置为 【L00/R00000095/C00000355】的一个瓦片,然后通过L00从conf.xml中相应的solution,这个瓦片的位置为149行853列,瓦片为256*256像素,所以可以算出瓦片的经纬度,故显示在相应的位置。

上面这是我的猜想,然后通过修改L00中一个瓦片的名称,发现这个瓦片的位置也就显示在不同的地方了,故得出上面的猜想是正确的。

Conf.xml文件内同如下:

<LODInfos xsi:type=’typens:ArrayOfLODInfo’>

<LODInfo xsi:type=’typens:LODInfo’>

<LevelID>0</LevelID><Scale>225429.3893129775</Scale><Resolution>0.001373291015625</Resolution></LODInfo><LODInfo xsi:type=’typens:LODInfo’>

<LevelID>1</LevelID><Scale>112714.6946564888</Scale><Resolution>0.0006866455078125</Resolution></LODInfo><LODInfo xsi:type=’typens:LODInfo’>

<LevelID>2</LevelID><Scale>56357.34732824438</Scale><Resolution>3.4332275390625e-4</Resolution></LODInfo><LODInfo xsi:type=’typens:LODInfo’>

<LevelID>3</LevelID><Scale>28178.67366412219</Scale><Resolution>1.71661376953125e-4</Resolution></LODInfo><LODInfo xsi:type=’typens:LODInfo’>

<LevelID>4</LevelID><Scale>14089.3368320611</Scale><Resolution>8.58306884765625e-5</Resolution></LODInfo><LODInfo xsi:type=’typens:LODInfo’>

<LevelID>5</LevelID><Scale>7044.668416030548</Scale><Resolution>4.291534423828125e-5</Resolution></LODInfo><LODInfo xsi:type=’typens:LODInfo’>

<LevelID>6</LevelID><Scale>3522.334208015274</Scale><Resolution>2.145767211914063e-5</Resolution></LODInfo><LODInfo xsi:type=’typens:LODInfo’>

<LevelID>7</LevelID><Scale>1761.167104007637</Scale><Resolution>1.072883605957031e-5</Resolution></LODInfo><LODInfo xsi:type=’typens:LODInfo’>

<LevelID>8</LevelID><Scale>880.5835520038185</Scale><Resolution>5.364418029785156e-6</Resolution></LODInfo>

</LODInfos>

重启服务,完成配置过程,则将下载好的天地图瓦片发布为自己的地图服务:

转载注明出处  http://www.cnblogs.com/The-SEA/p/4364074.html

转载于:https://www.cnblogs.com/The-SEA/p/4364074.html