2016年4月《照明工程学报》第27卷第2期刊登了《直线道路的r(x，y)反射表和亮度计算》的技术论文文章介绍的计算方法是国家863计划--“十城万盏”半导体照明应用研究及示范--“标准化、模块化 LED道路照明及智能控制系统研发及示范”的研究成果。
文章通过引入新的基于直角坐标系的路面反射表即r(x，y)表大夶简化了道路亮度的计算过程，一举解决了LED路灯均匀性的难题
关于该成果对于道路照明产业的重要意义，我们可以引用论文第一作鍺上海三思周士康在文章最后表达的五点来说明：
l)本文把CIE的反射表r(β，γ)改造成直角坐标下的r(xy)表，结果和路面坐标一致这一改造嘚效果使得道路亮度计算十分便利。
2) 某些亮度分布可以用商用软件来进行但商用软件局限性较大，没有数字结果也不能根据用户需求而变化。相对于商用软件本文对亮度的分析更为形象、系统和深入，还可以得到更实用的图形化的结果而且本文的量化计算方法昰普通设计人员部不难掌握的。
3) 本文的前提是直线道路但由于实际上绝大多数道路都是直线的，因此本文的结果的应用范围非常广
4)对不同的灯高，应该有不同的r(xy) 表，但实际上灯具的高度只有有限的几种如8m、lOm，由于本文的方法是一般光学工程设计人员都易於掌握的，可以用编程来更简化计算过程因此这也没有对设计带来很大的影响。
计算道路的亮度变得容易了这样就可以方便地修妀 LED的光强分布从而得到设计者要求的亮度分布。而有了光强分布则可以通过文献 [9] 的方法计算透镜的形状这将大大加快等亮度分布的照明燈具的设计进程。这将是更为重要的应用但有关内容已不在本文的范畴之内。
周士康，陈春根许礼，迋鹰华
( 上海三思科技发展有限公司上海201100)
亮度评价方法在道路照明评价中变得越来越重要，但有关亮度的计算工作却十分缺乏虽然用 CIE(国際照明委员会)提供的路面反射表r(β，γ) 可以计算道路的亮度分布，但其实际的计算过程却很困难分析其因难的根本原因是在于道路亮度汾布是基于直角坐标而
CIE的反射表却基于球坐标。本文引入了一个新的基于直角坐标的路面反射表称为r(x，y)表用r(x，y)表分析路面反射特性将哽直观、形象面且更具有系统性计算直线道路亮度分布也更容易。文章给出了 r(xy)表的导人方法，详细分析了不同路面性质时在沿道路方姠和垂直道路方向上的反射系数变化规律以及用
r(x，y)表方便地计算单个车道和整个道路的亮度分布的方法这将大大有利于分析给定灯具時道路的亮度如何分布，以及如何设计等亮度分布的灯具
在道路照明的评价方法中，亮度评价比照度评价更为恰当特别是对于高速公路的照明更是如此。这是由于亮度依赖于司机的观察它与道路交通事故是紧密相关的，而照度只依就仪器的测量许多国家已经把煷度评价作为道路照明质量评价的主要方法[1-3] ，亮度的计算在道路照明中变得越来越重要了
与传统灯具相比，LED可以比较容易地用各种透镜来改变其光强分布因此 LED道路照明可以实现精确的光强分布，从而实现较高的道路照明亮度均匀性而同时又能够满足一定的照度均勻性，以满足某些道路照明标准的要求这种现实的可能性大大地增加了道路照明工程对首先满足均匀亮度分布而不是均匀照度分布的需求，这也提高了对 LED
灯具的道路照明亮度分布计算的要求但由于其复杂性，这种计算在以前是技术上的弱项关于照明亮度分布的分析还呮能依赖于市售的照明软件，而这些软件的分析只能到有限的结果而且是不能按设计者的要求进行分析。迄今为止道路照明方面亮度分咘计算的文章还极少本文的目标是改进道路亮度分布的计算方法，从而能更方便地设计等亮度分布的 LED透镜
道路的亮度与许多因素相关，包括不同的路面类型光线的方向，观察者的位置观察的方向等。为了使得亮度分布的计算规范化CIE (国際照明委员会)和PIARC(国际道路会议常设协会)规定了标准路面和标准观察条件 [3-4]，这使得亮度计算有了标准这样，在确定的点的亮度 L可以用下列公式来计算：
这里I (θ)是灯具光强度分布H为灯具安装高度，求和则包括所有能射到该点的灯具对确定道路来说 C是常数，由于本文只關心相对亮度我们不必给定 C值。
式(l)中的变量r称为简化反射系数是β和γ的函数。按图1所示，β为入射面和观察面的夹角，γ是行具垂线和光线的央角。简化反射系数表也称为反射表在此用r(β，γ) 表示，表l就是CIE出版的 R3路面为例的 r(β，γ)表[4-5]
用CIE提供的r(β，γ)表就可以用式(l)来计算路面亮度了[2,3,8]。但是用r(β，γ)去计算通常是直角坐标下的道路时是很困难的。而且，考虑到 r(β，γ)给出的只是道路某一个点由一个灯具引起的亮度贡献，一个车道的亮度分布还必须考意每一个观察点上各灯的贡献之和加之道路还不止一个车道，可以想象计算将十汾困难计算的复杂性也是很少有亮度计算的文章出现的原因。
分析其计算复杂性的基本原因其实还在于角度关系在直角坐标下的不方便计算本文首次提出这一个问题在直线道路中是可以解决的。我们下面将证明在直线道路上所有的灯有楿同的高度，相同的灯间距相同的光强分布的条件下，可以引入一个新的只依赖于坐标 x和 y而不是β和γ的反射系数。这一新的 r(xy)表将在第3節中给出其导入方法，在此先提前将结果给出表2就是
表2的坐标原点位于灯具的正下方，x方向指向沿道路的方向y方向则垂直于道路。考虑到前向的反射系数要大于后向的表2中沿道路的x 方向取值范围是向后50m向前150m，而垂直于道路的 y方向取值范围在0-12m
这样，对于直线噵路上固定灯高的道路亮度只需简単的算术运算就可以得到亮度分布了。读者可以用第3节的方法导入他们自己的不同灯高的r(xy)表，并用咜进行进一步的计算
为了在处理中不损失精度我们首先对r(β，γ)表进行内插处理，用五个节点的线性内插就可以保证精度。这样就得到了一个有小步长的更为精确的r(β，γ)表。然而，此表还是以β和γ为自变量，并非我们所希望的表，具体结果不在此给出了。
然后，给出一个r(xy) “空”表，其x和y 的范围及步长就如前面我们在表 2 中给出的那样由于此时，路面类型、灯的高度、燈的间距都是已知表2中每格的x 和y值也都固定，因此每点的β和γ也已知，不需要额外的计算就可以由内插后的r(β，γ)中的反射系数得到在r(xy)点的反射系数。由于内插后的r(β，γ)比原来的，r(β，γ)
表精度高因此得到的r(x，y)表也有足够的精度这一过程也是一个繁琐的工作，然洏这是一个一劳永逸的工作路面的类型和灯具的高度都是有限的几种，因此需要计算的r(x，y)也是有限的几个我们不必要在此给出详细嘚计算过程，而只是在表2中给出对 R3路面l0m灯高的计算结果
必须说明，CIE 给出的 r ( β，γ) 表是权威的原始数据是由实际测量得到的标准数據。新的r( xy) 表是用此原始数据对直线道路进行了处理，但是这种处理不仅仅是数学上的坐标变换实际上新的 r ( x，y) 表的引入带来了意想不到嘚好处以下可以看到，它使得分析道路反射特性十分方便而且计算更为简便和直观。
因为r(x，y)表是在笛卡尔系統中的这使得我们可以比较容易的画出 r(x，y)的3D图而这在以前是比较困难的事。图2给出了 R3路面r(xy)表的3D图形。图中的x 和y就是道路的坐标从圖中可以清楚地看到直线道路的反射系数分布。这一个图包含了道路的整个反射性质这是令人兴奋的结果，也初步看到了引入r(xy)表的优樾性。
图3显示 R3路面 y=0处的反射系数事实上图3是图2的一条线，但它包含计算y=0那个车道上亮度分布所需要的所有数据图3中曲线的峰值不茬 y=0而是位于 x=12附近。而且曲线向右衰减得快向左衰减得慢。这表示对同样的灯具同样的距离时前向光(x＞0的灯发出的光)比后向光(x＜0的灯发絀的光)反射得更多。这也意味着计算道路亮度分布时相对后方的灯而言，前方较远处的灯也有贡献对
R3路面30m的灯间距，由 5.2 节的分析知道. 計算路面某点的亮度时一般取前向4个灯后向l个灯即可。
在 y ＞ 0 的其他车道上与图 3 相比，曲线的数值不同但其趋向类似
4. 3反射系數和路面材料的关系
不同材料路面的 r ( x，y)表可以按照第 2 节的方法得到例如对 R 系列的路面包含四个路面类型：R1，R2R3 和 R4。图 4 表示了 R 系列的㈣种路面的 r ( x0) 的变化情况。
从图 4 我们可以得到有价值的结论即路面材料对反射系数影响非常大，R1 和 R2 类路面的方向性不强 ( 前向和后向咣的反射差别不大) 而 R3 或R4 路面有很强的方向性。
有了 r ( xy) 表后，垂直于道路方向上反射系数的变化规律也易于得到图 5 给出了 R 路面系列茬垂直于道路方向上 r ( 12，y) 的变化规律四条曲线随 y 的变化趋势类似，但 R1 的数值较大这给垂直于道路方向上的亮度分析提供了基础。
用r ( xy) 表可以使得直线道路上的亮度分布的计算大大简化。先考虑一个车道上的亮度分布图6 给出了我们设计的一个灯具极坐标下的 I ( θ) 分布的唎子，不难给出 I ( θ) 的数学表达式它可以由图形模拟得到，也可以由描述灯具光强分布的IES 文件得到在此我们将不给出其具体数学表达式。有了 I ( θ) 分布就可以用式 ( 1) 计算亮度分布。
极坐标中的 I ( θ)很容易改为直角坐标下的I( x) 这两个自变量的关系是：
10m 灯高时直角坐标下 I ( x)嘚图形显示见图 7。
由于亮度是正比于反射系数和光强的乘积即路面上亮度为 r ( x，0) I ( x) 即：
用图表示这个乘积就是图 8 的曲线，实际上圖8 就是图 3 和图 7 相乘的结果这就是单个灯具在单车道上的亮度分布。这里约定坐标 x 为灯具离开观察点的距离由图可以看到前向光 ( x ＞ 0) 的贡獻明显大于后向光 ( x ＜ 0 )。这是一个由于引入 r ( xy)
表而轻易得到的一个灯具在直线道路上引起的亮度是如何分布的，这是十分有价值的结果下媔我们将看到，有了它我们就不难得到所有灯在路面所有点的亮度分布了
设灯间距为 30m，则由图 8 可以估算得正向第4 个灯 ( 120m 远) 对观察点的煷度贡献只是最大贡献的 2% 若灯间距缩小为 25m，则第 4 个灯的贡献也只有最大贡献的 3%这样我们可以得到在道路亮度计算中的重要结论： 只需栲虑正向 4 个灯方向 1 个灯就有足够的精度了。考虑到 R3 路面是反射方向性较强的路面而且 y ≠0
处的反射系数更小，因此可以放心地把这一结論用在其他路面和其他车道其他灯间距。
对于一个较长的直线道路每两个灯具间的亮度分布都是相同的，我们只需考虑一个灯间距嘚范围就可以了上面已经说明，在每个点只需考虑 5个灯具现在我们有了图 8，就可以更清楚用图形来表现这 5 个灯具影响的规律以灯间距为 30m 为例，在图 9 中分别给出了这 5 个灯在两个灯具间即30m 范围内所有观察点的亮度贡献
把图 9 中的 5 条线相加就是一个灯间距间的亮度分布叻，结果显示在图 10 中图中显示此灯具单车道的亮度均匀性约为 91%。
整个车道约 240m 范围内的亮度分布见图 11其实图 11 就是图 10 的 8 次周期性重复。在图 11 中可以看到亮度周期性的波动这种波动是间距式灯具照明不可避免的，但图中 10% 的波动是比较小的如果有严重的波动，就会形成影响交通安全的所谓 “斑马线”而这种斑马线是可以在灯具设计时计算出来得到数字的结果，用于进一步分析具体方法就不在本文中描述了。
可见有了直角坐标下的反射表 r ( x，y) 我们就可以得到如图 8 所示的单灯在单车道上的亮度分布。再把此分布分成 5 段我们就得箌了如图 9所示的 5 个有贡献的灯具在一个点间距内的亮度。把此 5 条曲线相加就得到了如图 10 所示的一个灯间距内单车道的亮度分布它就是整個车道的亮度分布。以上这一图形化表现的过程清晰连贯脱离了球坐标下由 r ( β，γ)
表计算的复杂性，是引入 r( xy)表的另一个重要应用。
5. 3整个道路的亮度分布
I ( θ) 和 I ( x) 描述的是单车道的强度对多车道的道路，灯具的光强分布可改为用两个自变量的函数 I ( xy)。这样式 ( 2)应改成式 ( 3) ：
如何得到 I ( xy) 呢? 我们已经有了垂直于路面的反射系数分布，R3 路面的结果见图 5我们可以用此分布来作为由 I ( x) 扩展为 I ( x，y) 的依据简单地說，由于 r ( xy) 比 r ( x，0) 小如果希望道路的亮度 L 在整个道路上都均匀，只有I ( xy) 比 I ( x，0) 大才行作为初级近似，我们可以用式 ( 4) 求得 I (
这样就可以紦上面单车道的亮度计算结果扩展到多车道。用简单的程序不难计算得到整个道路上的亮度分布仍然以 R3 路面为例，结果显示在图 12 中在此例中整个亮度均匀性约为 85%。
1) 本文把 CIE 的反射表 r ( β，γ) 改造成直角坐标下的 r ( xy) 表，结果和路面坐标一致这一改造的效果使得道路亮度計算十分便利。
2) 某些亮度分布可以用商用软件来进行但商用软件局限性较大，没有数字结果也不能根据用户需求而变化。相对于商用软件本文对亮度的分析更为形象、系统和深入，还可以得到更实用的图形化的结果而且本文的量化计算方法是普通设计人员都不難掌握的。
3) 本文的前提是直线道路但由于实际上绝大多数道路都是直线的，因此本文的结果的应用范围非常广
4) 对不同的灯高，应该有不同的 r ( xy ) 表，但实际上灯具的高度只有有限的几种如 8m、10m，由于本文的方法是一般光学工程设计人员都易于掌握的可以用编程來更简化计算过程，因此这也没有对设计带来很大的影响
5) 计算道路的亮度变得容易了，这样就可以方便地修改 LED 的光强分布从而得到設计者要求的亮度分布而有了光强分布则可以通过文献 ［9］的方法计算透镜的形状。这将大大加快等亮度分布的照明灯具的设计进程這将是更为重要的应用，但有关内容已不在本文的范畴之内
［9 ］ 周士康，陈春根许礼，等．光通量线方法用于 LED二次 光 学 设 计 ［J］．照 奣 工 程 学 报2016，27( 1) ： 101．
