地面沉降光栅立体图研究制作

图1 光栅示意图

Fig.1 Sketch of lenticular

1.1 柱镜光栅图的立体成像原理

人的两眼在接收图像信息时,会根据物点光线的入射方向,判定该物点的方位及远近,立体视觉是由于双眼在观察物体时产生的双眼视差,经过大脑处理后而产生。光栅立体图像之所以能给人以立体感,就是基于其处理后的平面图像包含有符合人眼视觉特性的立体信息,同时又能被柱镜光栅精确地折射分像^[9](图2)。

图2

图2 光栅图像立体感形成原理

Fig.2 Principle sketch of stereo feeling of lenticular image

图2(a)为普通平面印刷产品,印刷面即为像平面,由于双眼观察a,b和c像点时视线交角均位于像平面上,人眼会判定平面印刷品上所有像点均在同一水平面上,没有立体感。而在图2(b)中,显示了在光栅底图上覆盖光栅板的情形。在光栅板的折射作用下,人的左眼看到像点a,b和c,右眼看到相应的立体像对同名点a',b'和c',由于双眼观看时视角不同而产生距离感: A浮在像平面上方,形成近景; C位于像平面下方,形成远景; B位于像平面上,形成中景。

设人的双眼距离为 $e$ ,眼睛到光栅立体图的观看距离为 $L$ ,立体像对中2个同名像点间的位错为 $Δ x$ ,那么视觉成像点与光栅立体图平面距离y为^[4]

(1)

y = \frac{L}{e \pm Δx} \overset{}{\sum_{}} Δx

式中: “ +”为成像在像平面上方; “ -”为成像在像平下方。

柱镜光栅成像的理论基础是利用光栅对底图进行折射和分像,再配合人的双眼视差形成立体感,因此有必要理解光栅底图的制作原理。传统光栅底图是由同一景物在不同视角下拍摄的图像按一定规则抽样、组合而成。理论上讲,只有2个或2个以上的互为立体像对的一组图像方可合成立体图像 ${^{[}}^{6}^{]}$ ,原理示意图如图3所示。

图3

图3 图像合成光栅底图原理示意图

Fig.3 Principle sketch of lenticular base image composed by 3 images

以3幅图像为例,图D,E,F表示同一景物在不同角度下互为立体像对的图像,将每一幅图像分割为3个条带,将切割的9个条带交织起来,这样每个光栅单元下都有3幅图像的一部分信息,整个光栅底图则包含了3幅图像的全部信息。光栅底图的制作原理决定了人眼直接观看光栅底图时只看到模糊不清的竖状条纹影像,只有在相应参数柱镜光栅的折射分像下,将抽样组合的信息还原为原始的立体像对,分别进入左右眼才能再现立体感。

由于柱镜光栅是由在平面上呈线性排列的多个圆柱透镜组成,这种结构使我们只能感觉到水平向的视差信息,这也决定了光栅底图的3个基本特点^{[11,12,13,14,15]}: ①立体图像必须是由对同一景物从不同视角获取的一组图像组成,由于光栅的单向立体特性,要求此组图像具有拍摄高度相同的特性; ②所合成的立体图像由平行于柱状光栅轴向的条纹组成,各个局部条带的交错放置情况如图3所示; ③为得到好的立体效果,在抽样组合时需首先实现主体景物成像在纸面上,亦即主体景物上点的位错(图2中的 $Δx$ )为0。

1.2 关键参数

光栅的特征是用特定的参数来描述的,主要包括: 光栅厚度、光栅线数、光栅视距、光栅透光率、光栅偏差值、光栅像素打印点和折射率等^{[6, 11-15]}。光栅厚度即为透镜的焦距。光栅线数(lpi)指在与柱镜延长向垂直的方向上,一英寸(约2.54 cm)距离内柱镜的条数。光栅线数与栅距成反比。光栅视距D指人与光栅图像之间可达到最佳立体感的距离,即最佳观察距离。D与光栅折射率n、光栅焦距 $f$ 、人眼瞳距e及栅距p,满足公式^[16]为

(2)

D = \frac{e}{2 \tan {\arcsin [n \sin (\arctan \frac{p}{2 f})]}}

由式(2)可以得出,在折射率和焦距一定,并假设人眼瞳距不变的情况下,不同线数的光栅材料,其对应的光栅视距是不同的: 栅距越小,光栅线数越密,光栅视距越小; 栅距越大,光栅线数越疏,光栅视距越大。因此近距离观察用的光栅板是薄且密的,远距离观察用的光栅板是厚且疏的^[2,16]。

2 数据源及其处理

2.1 数据源

地面沉降光栅立体图制作的源图是一幅包括地面沉降年均速率等值线并对等值面分级设色的数字平面图,名为“J50幅地表形变分布图”,具体信息如表1所示,地面沉降专题信息的提取和制图方法见参考文献[17,18,19]。

表1 地面沉降光栅立体图源数据信息

Tab.1 Summary information of the source data used to make the lenticular stereoscopic subsidence map

数据信息类型	数据信息内容
地表形变信息来源	ENVISAT ASAR时间序列数据
提取方法	相干目标短基线时间序列分析
基础地理底图	国家地理信息1:100万数字地形图
投影坐标	高斯-克吕格投影、西安80坐标系
数据周期	2008年1月—2010年10月
制图时间	2011年
制图单位	中国国土资源航空物探遥感中心

新窗口打开| 下载CSV

2.2 处理流程

柱镜光栅立体图的制作流程^{[3, 16]}主要分为4个步骤,即材质选择、制作光栅底图、印刷和印后处理(图4)。

图4

图4 光栅图制作流程

Fig.4 Flowchart of the production of a lenticular image

1)材质选择。根据材质的不同可以将光栅分为PP(聚丙烯)光栅、PS(聚苯乙烯)光栅、PVC(聚氯乙烯)光栅和PET(聚对苯二甲酸乙二醇酯)光栅。当前在制作较高精度的光栅立体图时,一般采用PET光栅板。这种材质环保无毒,透光度高,化学性能稳定,且具有耐酸、耐碱、不亲水、不亲油等特性。为了使地面沉降光栅立体图的三维效果尽可能细腻,使用高光栅线数的柱镜光栅,其光栅线数约为90.45 lpi,栅距约为0.28 mm。光栅材质为PET,折射率约为1.59,厚度约为0.60 mm,成图后适合手持近距离观看。

2)制作光栅底图。光栅底图是由图像序列抽样合成的,因此图像序列获取质量的高低直接影响到光栅底图的精度。图像序列的获取方法主要有4种: ①实景拍摄,在数字图像处理技术普及之前使用较多,是早期获取图像序列的主要手段,需要立体相机等特殊的设备,限制因素多; ②分层,原理是不断复制移动一张图像,通过图像的位移差来近似模拟角度视差,位移量越大,模拟的角度视差越大,景深的绝对值就越大,凹凸效果也就越明显; ③图像立体化,利用专业的图像处理软件根据物体的远近关系,勾画等距线进行分层切片,模拟生成多视角图像序列; ④仿实拍立体,通过现代计算机技术,在平面坐标要素的基础上,利用第三维度的等值线构建该维度上的数字模型,形成一种计算机虚拟立体显示系统,在系统中可以从任意角度来观测模型,获取图像多视角的信息。在光栅立体图的制作流程中,光栅底图的处理是一个关键环节。光栅立体图的推广在很大程度上得益于光栅底图合成技术的快速进步。随着数字图像处理技术的飞速发展,光栅立体印刷的图像处理也已进入了数字图像处理合成的阶段,利用计算机完全可以通过带有第三维信息的平面图像合成光栅底图,得到与实景拍摄相似的效果。因此,第③和④种方法,即图像立体化与仿实拍立体方法,成为了目前主流光栅底图的制作方法。

3)印刷和印后处理。立体印刷的主要方式有3种: ①普通胶印,首先印刷立体光栅底图,待纸张上的油墨干燥后,再在底图上覆盖一层与之参数相匹配的柱镜光栅板,然后进行光栅复合; ②丝网印刷,是将立体光栅板印刷在印好的立体图像表面,即把光栅看做图文进行印刷,印后需要模压成型; ③UV(紫外)印刷,直接在柱镜光栅板的背面(为一平面)上印刷光栅图像,配合紫外源在油墨干燥后覆以油墨保护层即可。

2.3 关键步骤

在制作光栅底图时,应用了仿实拍立体方法。首先在ArcGIS平台中,将J50幅地表形变分布图的地面沉降速率值作为模拟三维形态的依据,通过调整夸张因子,获得具有较好效果的三维场景。在此过程中,为了使观察者对不同因子下的三维场景具有更为形象的感受,将各种三维场景输出为红绿三维效果图,观察者佩戴红绿三维眼镜进行比较。以地面沉降显著的沧州市附近局部区域为例,源图和其对应的红绿三维效果如图5所示。

图5

图5 地面沉降显著区域源图和红绿三维图

Fig.5 Original map and red-green anaglyph of severe subsidence area

经比较确定满意的夸张因子和三维场景后,对场景从13个角度获取序列图像,相当于模仿了从13个角度进行实景拍摄的镜头,部分角度图像如图6所示。

图6

图6 不同观察角度的地面沉降场景序列

Fig.6 Series viewing images simulated for different angles of a subsidence scene

由图6可以发现,在固定地理位置的小窗口内,相同的地面沉降速率等值线在不同观察角度的镜头下的窗口内逐渐向右移,且形状也稍有一定程度的扭曲变化。

根据光栅底图分辨率等于光栅线数与镜头数的乘积,设置光栅底图分辨率约为1 175.82 dpi,按图3所示的原理将序列图像抽样合成为光栅底图,如图7所示。由图7可见该研究区光栅底图的合成效果,在地面沉降较为明显的地区是条带状模糊不清的影像; 在无形变地区,光栅底图与源图件一致。

图7

图7 地面沉降显著区域光栅底图

Fig.7 Lenticular base map of severe subsidence area

对光栅立体图的印刷及印后处理,采用了UV印刷和保护层处理技术,光栅底图和光栅板匹配程度高,立体效果明显。

受当前主流柱镜光栅材料特性所限,此次立体图的光栅线数仍不够大,对于地学要素的专业表达而言,清晰度与传统纸质图件相比尚有差距。此外,在本次立体图制作过程中还发现,目前柱镜光栅材料在立体图的幅面尺寸方面也有限制,这给大区域、大幅面地学现象的制图带来制约。相信随着光栅材料的发展,这2方面的缺点都将被克服。

3 结论与讨论

本文在阐述柱镜光栅立体图制作原理的基础上,结合现代数字图像处理技术,利用已有的“J50幅地表形变分布图”二维成果图,研究并制作了1:100万标准图幅的地面沉降光栅立体效果图。与已有的常规柱镜光栅立体产品相比,本次制作的立体效果图有3方面特点:

1)对地面沉降凹陷形态的表达形式做了新的探索。传统的人物、风景立体图多侧重表达在像平面之上的凸起形态,与地面沉降特征相比,图幅内场景立体特征更为明显; 当前常见的地形地貌立体图,所表达的起伏现象在变化剧烈程度上也比地面沉降更为明显,易于取得较好的立体视觉效果。本次对地面沉降凹陷形态的立体表达体现了新的特点。因此,为了理想表达地面沉降凹陷形态,在三维场景夸张程度调节、光栅底图远中近景设置等方面都做了成功探索。

2)不同图层之间的匹配精度高。以往三维光栅立体图像制作工艺多应用于摄影、广告等对图元间几何匹配精度没有严格要求的人文艺术领域,而本次光栅立体图包含了基础地理要素、地面沉降等值线和等值面等具有严格地理空间对应关系的不同图层,匹配精度明显高于以往传统产品。

3)立体图幅面大。本次光栅立体图成品尺寸约为700 mm×500 mm,在目前国内同类材质制品中较为突出。

研究表明,对大范围地面沉降的凹陷特征,可基于三维地理信息系统、数字图像处理、光栅立体成像与印刷技术,进行立体表达,图件新颖、直观、便于使用,有着较为广阔的应用前景。随着光栅材料的发展,制作幅面更大、分辨率更高的光栅立体图将成为可能。

本文将光栅立体成像技术引入地面沉降制图方向加以研究,尝试表现地面沉降的三维效果,取得了不错的效果,但也还存在着不足之处,主要表现在受光栅材质的制约,光栅线数仍不够多,清晰度仍不及传统二维图件,幅面尺寸也有所限制。但这一技术与地图制图行业的结合已经表现出了较大的优势和潜力。未来该技术在地图制图领域的发展和推广,以及在其他学科、行业的探索和研究仍需要做大量的工作。

志谢: 感谢北京万维时空科技有限公司邢耶在地面沉降专题图编制方面所做的工作。

参考文献

原文顺序

文献年度倒序

文中引用次数倒序

被引期刊影响因子

[1]

Roberts D

History of lenticular and related autostereoscopic methods[Z]

Hillsboro:Leap Technologies,LLC, 2003.

DOI:10.3969/j.issn.1674-5752.2009.05.001 URL [本文引用: 2]

[2]

史瑞芝

光栅立体印刷技术综述

[J]. 中国印刷与包装研究, 2009,1(5):1-9.

近十年来，立体印刷技术得到了迅猛发展，在立体印刷理论、设备、软件平台和制作工艺技术的研究、市场开拓等方面均取得了显著成果，在印前图像制作、数据输出、印刷成图工艺和光栅材料等方面均有所突破，并开始步入数据化、规范化生产的轨道。本文阐述了立体成像的视差理论及像对立体成像法和光栅立体成像法的原理，回顾了立体成像技术的发展历程，重点阐述了光栅立体印刷技术过程及关键技术，并对该技术未来的发展进行了展望。

Shi R

Summary of grating stereoscopic printing technology

[J]. China Printing and Packaging Study, 2009,1(5):1-9.

[3]

Johnson R

, Jacobsen G

Advances in lenticular lens arrays for visual display

[C]//Mouroulis P Z,Smith W J,Johnson R B.Current Developments in Lens Design and Optical Engineering VI.San Diego:SPIE, 2005.

[4]

Buchroithner M

, Knust

True-3D in cartography-current hard- and softcopy developments

[M] //Moore A,Drecki I.Geospatial Visualisation.Berlin Heidelberg:Springer, 2013: 41-65.

[5]

Moellering

Perspectives on 3D visualization of spatial geodata and future prospects

[M] // Buchroithner M.True-3D in Cartography.Berlin Heidelberg:Springer, 2012: 1-19.

[6]

史瑞芝, 安敬 .

一种基于柱镜光栅成像技术制作立体影像图的方法

[J].测绘通报, 2008(11):31-33.

URL [本文引用: 2]

论述光栅立体成像原理,以及采用此原理制作光栅立体影像图的技术方法。将一组有轻微角度差别拍摄的系列图像,按一定的规则抽样合成一幅图像,经柱镜光栅折射后进入人眼,从而产生立体感。在适合制作立体影像图的制图区域上空按照一定的要求拍摄影像,经过计算机对影像进行匹配与处理,制作成光栅图像。最后,将柱镜光栅与光栅图像精确对位并紧密压裱,形成立体光栅影像图。

Shi R

, An

A method of composite stereo photomap based on lenticular screen

[J].Bulletin of Surveying and Mapping,2008(11):31-33.

DOI:10.3969/j.issn.1007-791X.2010.03.016 URL [本文引用: 1]

[7]

彭爱华, 朱化凤 .

应用于立体印刷的矩阵式透镜光栅的原理

[J]. 燕山大学学报, 2010,34(3):262-265.

立体印刷中应用的透镜光栅为柱透镜光栅。本文探索了将矩阵式透镜光栅应用于立体印刷的原理,通过比较矩阵式透镜光栅和柱透镜光栅,指出矩阵式透镜光栅的优点。矩阵式透镜光栅制作的立体画比柱透镜光栅制作的立体画的观察角度更多,能使平面印刷图像取得更好的立体效果。本文给出了具体的制备参数。

Peng A

, Zhu H

Principium of matrix lens grating in three-dimensional printing

[J]. Journal of Yanshan University, 2010,34(3):262-265.

[8]

李鹏

立体印刷的成像原理

[J].印刷质量与标准化, 2007(3):8-11.

The imaging elements of 3D printing

[J].Printing Quality and Standardization,2007(3):8-11.

[9]

周立权

光栅立体印刷图像处理技术研究

[J]. 包装工程, 2010,31(3):107-110,136.

Zhou L

Research on image-processing technology of stereoscopic printing

[J]. Packaging Engineering, 2010,31(3):107-110,136.

[10]

董永贵, 沈立, 冯冠平 , 等.

一种基于柱透镜光栅的计算机辅助彩色立体图片合成方法

[J].光学技术, 1999(3):66-68.

DOI:10.3321/j.issn:1002-1582.1999.03.016 URL

对基于柱透镜光栅的立体照片的制作原理进行了理论分析，推导了了立体成像的光学原理，结合现代数字图像处理技术，对多幅图像的匹配及合成方法进行了分析，提出了和实现了一种计算机辅助彩色立体图片合成方法。

Dong Y

, Shen

, Feng G

, et al.

Computer-aided composite method of the lenticular screen covering color stereo-image

[J].Optical Technology,1999(3):66-68.

[11]

康立民

立体照片的拍摄及合成

[J]. 北京机械工业学院学报, 1996,11(2):30-42.

Kang L

Discusses shoot and composite of the stereogram

[J]. Journal of Beijing Institute of Machinery Industry, 1996,11(2):30-42.

[12]

大越孝敬. 三维成像技术[M]. 董太和译.北京: 机械工业出版社, 1982.

Takanori. Three-Dimensional Imaging Techniques[M]. Dong T H, transl.Beijing:China Machine Press, 1982.

[13]

侯卫兵

大幅面彩色喷涂系统及相关技术的研究

[D]. 北京:清华大学, 1997.

[14]

Ignat’ev N

Two modes of operation of a lens array for obtaining integral photography

[J]. Soviet Journal of Optical Technology, 1983,50(1):6-8.

[15]

Pollard S

, Mayhew J E

, Frisby J

Implementation details of the PMF stereo algorithm

[M] //Mayhew J E W, Frisby J P.3D Model Recognition from Stereoscopic Cues.Cambridge:MIT Press, 1991.

[16]

黄敏, 刘浩学 .

立体印刷制作工艺光栅板参数的选择

[J]. 北京印刷学院学报, 2005,13(1):9-12.

[本文引用: 3]

Huang

, Liu H

Selecting of the parameters of cylindrical lens sheet imaging in 3D printing

[J]. Journal of Beijing Institute of Graphic Communication, 2005,13(1):9-12.

[本文引用: 3]

[17]

范景辉, 李梅, 郭小方 , 等.

基于PSInSAR方法和ASAR数据监测天津地面沉降的试验研究

[J]. 国土资源遥感, 2007,19(4):23-27.doi: 10.6046/gtzyyg.2007.04.05.

URL Magsci [本文引用: 1]

概述了PSInSAR方法的技术流程，阐述了定标ENVISAT ASAR数据对PS候选点（Persistent Scatterer Cadidate，PSC）筛选的作用。研究表明：对ENVISAR ASAR 数据进行定标有助于筛选出更多的PS候选点；以定标得到的后向散射系数作为阈值，可以剔除散射信号统计特性较稳定但散射强度低的像素点，从而避免可能由这些点引入的相位误差。在初步实现PSInSAR方法的基础上，运用14景ENVISAT ASAR数据获得了天津地区的年均线性沉降速率，揭示的天津市地面沉降趋势与前人研究结果较为一致，获得的形变速率值的准确性尚待数据量的增多而进一步提高。

Fan J

, Li

, Guo X

, et al.

A preliminary study of the subsidence in Tianjin area using asar images based on psinsar technique

[J]. Remote Sensing for Land and Resources, 2007,19(4):23-27.doi: 10.6046/gtzyyg.2007.04.05.

Magsci [本文引用: 1]

[18]

王润生, 熊盛青, 聂洪峰 , 等.

遥感地质勘查技术与应用研究

[J]. 地质学报, 2011,85(11):1699-1743.

Magsci [本文引用: 1]

遥感的发展极大地拓宽了人类的视野和视觉能力，以其宏观性、综合性、多尺度、多层次的特点，已成为地质研究和地质勘查不可缺少的方面军，在地质调查、矿产勘查、地质环境评价、地质灾害监测和基础地质研究等方面都发挥了越来越大的作用。随着传感器分辨率（空间、光谱、时间、辐射）的不断提高，特别是高光谱和干涉雷达技术的发展，不仅极大地提高了遥感的观测尺度、对地物的分辨本领和识别的精细程度，而且使遥感地质发生了由宏观探测到微观探测，由定性解译到定量反演的质的飞跃，将遥感地质和应用都推向一个新的高度。本文首先简要介绍了遥感技术的发展及其意义，遥感地质找矿的原理、方法和遥感找矿模型，然后分别从岩矿波谱、遥感图像处理、多光谱蚀变异常提取、高光谱矿物识别与填图、InSAR地表形变调查与监测、遥感找矿模型、遥感地质灾害调查与监测7个方面介绍了遥感的基本原理、主要的技术方法、技术特点和技术关键、主要的应用领域，以及取得的一些研究成果和应用实例

Wang R

, Xiong S

, Nie H

, et al.

Remote sensing technology and its application in geological exploration

[J]. Acta Geologica Sinica, 2011,85(11):1699-1743.

Magsci [本文引用: 1]

[19]

范景辉, 燕云鹏, 葛大庆 , 等. 全国地表形变遥感地质(InSAR)调查技术指南[M]. 北京: 地质出版社, 2015.

Fan J

, Yan Y

, Ge D

, et al. Technical Guidance on Remote Sensing Geological Survey Based on InSAR for Ground Deformation in China[M]. Beijing: Geological Publishing House, 2015.