张敏-基于FME的矢量地理底图数据分幅及栅格化.pdfVIP

基于FME 的矢量地理底图数据分幅及栅格化 武汉市测绘研究院 张敏 1 背景 在日常工作中，我们有时会遇到需要将大范围内多要素的地理底图数据根据特定的要求 分幅裁切的情况，并要求各要素保持原有的各种几何、属性信息。现有此类情况，需要将武 汉市全市域范围内MDB 格式的整张多要素地理底图数据裁切为1:2000 标准分幅的MDB 文 件，并配置线型、符号、注记等对其进行栅格化，输出为带坐标的图片文件，并制作瓦片数 据导入移动终端供外业调查使用。 2 数据概况 2.1 全市域多要素mdb 矢量数据 主要数据源为武汉市全市域范围的地理地图MDB 格式矢量数据，包括有：地名、城市 道路、公路、轨道交通、铁路、建筑、大型市政设施、水体、山体、绿地等类别，其中每个 大类中又包含若干个小要素类，比如城市道路分为城市主干路、城市其他道路等。这些矢量 数据中，都包含各自必须的属性信息，后面出图时可通过这些属性来生成文字标注。 图2.1 全市域地理地图要素 （局部） 2.2 全市域1：2000 标准分幅格网 为了将整张的全市域范围的地理底图进行分幅裁切，需要借助全市域的 1:2000 标准分 幅格网 shapefile 数据来完成。每个格网的属性信息中都包含有对应位置标准分幅文件的文 件名，即filename 字段的值。进行分幅文件输出的时候，需要将裁切后的要素按照对应的位 置关系拆分到对应的标准分幅2000 图中，并以规定的文件名来命名保存 （1:2000 比例尺地 形图数据文件名为8 位，前6 位分别取图幅坐标三位整公里数，X 坐标在前，Y 坐标在后。 第7、8 两位为0 ，如mdb）。 图2.2 1:2000 标准分幅格网 3 FME 实现 3.1 基本思路 3.1.1 分幅裁切矢量地理底图数据 运用 Clipper 转换器，利用标准分幅格网矢量数据，对整张的地理地图矢量文件进行裁 切，并保存为1:2000 标准分幅的多个 MDB 文件。这些分幅文件，后期将被用于数据处理、 栅格化。 3.1.2 配置要素符号样式 运用MapnikRasterizer 转换器将矢量数据输出成栅格图像。在输出之前，需要对各要素 的符号、线型、颜色、注记等进行配置，以使输出的图片能满足工作需要。需要注意的是， 对于需要生成注记的要素，其设置的属性字段名要与其矢量数据中完全一致，否则无法按要 求输出文字标注。 图3.1 栅格化样式参数配置 3.1.3 关联格网与分幅底图 由于裁切后的分幅MDB 底图文件中的要素类型、范围等会有所不同，而输出的图片尺 寸是固定统一的，就会导致相邻的图片中的要素无法正常接上。如下，右侧图幅中要素如果 未能布满整个图幅范围，那么图片输出的时候就会自动放大使其布满整个图片范围，但左图 中的要素都充满了幅图，所有导致两张图片中的要素缩放比例不一致： 图3.2 图幅未接边示例 为了解决这一问题，需要将每幅图对应的格网一起输出，这样每幅分幅底图中的要素范 围都一样，输出要素的缩放比例也一致，因此接边图幅中的要素都能正确接边。如何才能从 全市域范围的格网中选中与某图幅匹配的格网呢？可以根据分幅MDB 文件名称，批量利用 where 语句在格网中进行查找匹配。先根据每个读取进来的分幅文件名来创建一个新的属性 （假设为a ），使其属性值等于该文件名。新增一个FeatureReader 转换器来读取格网数据， 在其读取的where 条件语句中规定读取格网中filename 字段值与属性a 的值相等的格网（前 面提到过，格网数据的filename 字段值代表的就是与该格网对应的标准分幅图的文件名）， 这样就能将与图幅对应的格网一起栅格化输出。 图3.3 字段a 定义 图3.4 条件参数设置 3.1.4 循环执行栅格化 MapnikRasterizer 转换器每次只能输出一张图片，无法输出多张，所以需要添加Creator 、 WorkspaceRu