FME中的栅格数据操作之十三——生成应用于网络制图平台的瓦片(三)
原文发布时间:2010-06-13
作者:毛毛虫
对于前面讲述的生成瓦片的处理,如果要处理的数据集不很大的时候,或许能顺利运行。但是通常要处理的数据会很大,这是会出现FME内存不足的情况。
我们这里再说一下前面我们生成瓦片的过程:先读入shp文件,然后选择一小部分需要生成瓦片的区域;之后栅格化这部分数据,然后生成瓦片输出。
这时存在下面一些问题:一是数据管理起来可能会很麻烦;二、产生的瓦片可能不完整。我们随机选择的区域(没有按Bing Maps网格的划分)可能会把一个瓦片分成几个部分,也就是这个瓦片一部分可能会在你选择的一个区域,而另一部分在另一个区域,这时候需要它们融合为一张。
我们改进下这个流程,即使用Bing Maps瓦片的边界作为边界来选择裁剪区域。同时,为了方便的管理数据,我们最好能采用一些方便的格式存储(最好是使用空间索引,如FFS或一些空间数据库的方式)。
这个改进的方法大致流程是:沿着Bing Maps瓦片的边界读入生成瓦片的区域,然后设置显示比例和颜色,使用转换器Clipper,用由Bing Maps瓦片生成的多边形裁剪设置好显示比例和颜色后的矢量数据,然后使用ImageRasterizer栅格化数据。
具体做法:创建两个工作空间,使用转换器WorkspaceRunner把两个工作空间连接起来。其中第一个工作空间用来生成裁剪区域,并运行第二个工作空间。第一个工作空间选定一个区域,当然用户可以在这里设置相关的参数(如显示级别)。再把这个区域发送到VirtualEarthTiler转换器中,生成需要Bing Maps的边界。然后从每个瓦片中提取最大/最小的X/Y的坐标值,作为运行第二个工作空间的参数。最用使用转换器WorkspaceRunner运行第二个工作空间。转换器WorkspaceRunner设置如下:
在FME WorkSpace参数中输入要运行的第二个工作空间,然后在参数列表中设置下第二个工作空间需要的参数值,就可以运行在这个工作空间中按设置的参数值运行第二个工作空间。
例如下图:
图中橘黄色区域是我们选择的区域,如果在5-7级进行切片的话,需要用到28个工作空间(用到的工作空间图中都用了Quadkey标识出来了)。对于这些瓦片的选择,我们可以使用Clipper转换器,在生成的5-7级的瓦片中选择落在提供的区域内部和跟该区域有交集的部分。然后提取每个选中瓦片的边界,即最大/最小的X/Y的坐标值传入到第二个工作空间。第二个工作空间的相关参数如下:
其中Minimum X、Minimum Y、Maximum X和Maximum Y就是第二个工作空间需要切片的范围。
第二个工作空间读入根据第一个工作空间设置的区域内的数据,对矢量数据进行显示级别设置和着色,裁剪数据最后进行栅格化。有时需要进一步切片数据。这时也是使用一个瓦片生成转换器VirtualEarthTiler和边界提取转换器BoundingBoxReplacer:
这是因为有一个工作空间数目与每个工作空间每秒切片情况的对比。
当工作空间太多时,电脑的运行开销也会激增,这样快速切片节省的开销无法抵消工作空间的消耗,也就达不到快速切片的目的,所有工作空间太多是不可取的。通过FME的相关测试,最佳的情况是每个workspace处理大概64-1000个瓦片的时候。
另一种途径是可以通过FME Server。今后我会写一个FME Server的相关专题再来讨论这个问题。