一种采用直角坐标下的路面反射表计算道路亮度分布的方法
一种采用直角坐标下的路面反射表计算道路亮度分布的方法
【技术领域】
[0001] 本发明设及道路照明领域，特别是设及一种采用直角坐标下的路面反射表计算道 路亮度分布的方法。
【背景技术】
[0002] 亮度评价方法在道路照明评价中变得越来越重要，但有关亮度的计算工作却十分 缺乏。虽然用CIE(国际照明委员会)提供的路面反射表Ηβ，丫）可W计算道路的亮度分布， 但其实际的计算过程却很困难。
[0003] 在道路照明的评价方法中，亮度评价比照度评价更为恰当，特别是对于高速公路 的照明更是如此。运是由于亮度依赖于司机的观察，它与道路交通事故是紧密相关的，而照 度只依赖仪器的测量。许多国家已经把亮度评价作为道路照明质量评价的主要方法，亮度 的计算在道路照明中变得越来越重要了。
[0004] 与传统灯具相比，L抓可W比较容易地用各种透镜来改变其光强分布，因此L抓道 路照明可W实现精确的光强分布，从而实现较高的道路照明亮度均匀性，而同时又能够满 足一定的照度均匀性，W满足某些道路照明标准的要求。运种现实的可能性大大地增加了 道路照明工程对首先满足均匀亮度分布而不是均匀照度分布的需求，运也提高了对Lm)灯 具的道路照明亮度分布计算的要求。但由于其复杂性，运种计算在W前是技术上的弱项，关 于照明亮度分布的分析还只能依赖与市售的照明软件，而运些软件的分析只能到有限的结 果，而且是不能按设计者的要求进行分析。迄今为止道路照明方面亮度分布计算的文章还 极少。
[0005] 道路的亮度与许多因素相关，包括不同的路面类型，光线的方向，观察者的位置， 观察的方向等。为了使得亮度分布的计算规范化，CIE(国际照明委员会)和PIARC(国际道路 会议常设协会)规定了标准路面和标准观察条件，运使得亮度计算有了标准。
[0006] 如图1所示，β为入射面和观察面的夹角，丫是灯具垂线和光线的夹角。简化反射系 数表也称为反射表，在此用Γ(β，丫）表示，下列的表1就是CIE出版的R3路面为例的Γ(β，丫） 表。
[0007] 表1CIE出版的R3路面的Γ(β，丫）表 [000引
[0009]
[0010] 用CIE提供的Γ(β，丫）表就可w用来计算路面亮度了。但是用Γ(β，丫）去计算通常 是直角坐标下的道路时是很困难的。而且，考虑到r(e，丫）给出的只是道路上某一个点由一 个灯具引起的亮度贡献，一个车道的亮度分布还必须考虑许每一个观察点上各灯的贡献之 和，加之道路还不止一个车道，可W想象计算将十分困难。
【发明内容】
[0011] 鉴于W上所述现有技术的缺点，本发明的目的在于提供一种直角坐标下的路面反 射表计算道路亮度分布的方法，用于解决现有技术中由于道路亮度分布是基于直角坐标而 CIE的反射表却基于球坐标，从而计算十分困难的问题。
[0012] 本发明提供的一种采用直角坐标下的路面反射表计算道路亮度分布的方法，包括 步骤:S1、将球坐标下的路面反射表Γ(β，丫）转换为直角坐标下的路面反射表r(x，y);S2、根 据所述直角坐标下的路面反射表Hx，y)分析路面特性;S3、结合所述路面特性根据所述直 角坐标下的路面反射表r(x，y)计算道路亮度分布。
[0013] 于本发明的一实施方式中，所述步骤S1包括步骤:S11、对所述球坐标下的路面反 射表r(0，丫）进行内插处理；S12、计算所述直角坐标下的路面反射表;r(x，y)中相应数据； S13、将道路中每一车道中有亮度贡献的灯具的所述球坐标下的路面反射表r(i3, 丫）的反射 系数相加得到对应点新的所述直角坐标下的路面反射表r(x，y)的反射系数。
[0014] 于本发明的一实施方式中，所述有亮度贡献的灯具的贡献值至少为2% W上;所述 相应数据包括路面类型、灯具的高度和灯具的间距。
[0015] 于本发明的一实施方式中，所述步骤S2还包括:不同的所述道路类型对应不同的 所述球坐标下的路面反射表r(e，丫），并根据不同的所述球坐标下的路面反射表Γ(β，丫）转 换为对应的所述直角坐标下的路面反射表r (X，y)。
[0016] 于本发明的一实施方式中，所述步骤S2还包括步骤:根据所述直角坐标下的路面 反射表r(x，y)画出路面的所述直角坐标下的路面反射表r(x，y)的3D曲线图。
[0017] 于本发明的一实施方式中，所述步骤S2还包括步骤:取所述3D曲线图中对应车道 的2D曲线图，并根据所述2D曲线图分析对应车道在车道延伸方向前后的灯光亮度衰减情 况。
[0018] 于本发明的一实施方式中，所述步骤S3还包括步骤:S31、计算单个灯具对道路亮 度分布的贡献;S32、计算单个车道的道路亮度分布;S33、计算整个道路的道路亮度分布。
[0019] 于本发明的一实施方式中，所述步骤S31还包括:根据公式：8:::燃邮音，将灯的极 坐标下的强度分布Ι(θ)转换为灯的直角坐标下的强度分布I(x);其中θ为光线与灯具法线 反向的夹角，Η为灯具安装高度。
[0020] 于本发明的一实施方式中，所述步骤S32还包括:单车道的亮度正比于反射系数和 光强的乘积;所述单车道的亮度分布的计算公式为:L(x)= Σ0 · (r(x，0) · Ι(χ))/Η2。
[0021] 于本发明的一实施方式中，所述步骤S33还包括:所述整个道路的亮度分布的计算 公式为：L(x，y)=IC· (;r(x，y) · I(x，y))/H2;其中，作为初级近似计算， 乃―= 户，C为常数。
[0022] 如上所述，本发明具有W下有益效果:方便地计算单个车道和整个道路的亮度分 布，大大有利于分析给定灯具时道路的亮度如何分布，W及如何设计等亮度分布的灯具等。
【附图说明】
[0023] 图1显示为路面反射的角度关系示意图。
[0024] 图2显示为本发明采用直角坐标下的路面反射表计算道路亮度分布的方法一实施 例的流程方框图。
[0025] 图3显示为本发明采用直角坐标下的路面反射表计算道路亮度分布的方法一实施 例的路面r(x，y)表的3D曲线图。
[0026] 图4显示为图3中y = 0时的2D曲线图。
[0027] 图5显示为本发明采用直角坐标下的路面反射表计算道路亮度分布的方法一实施 例的R系列沿路方向的反射系数示意图。
[0028] 图6显示为本发明采用直角坐标路面下的反射表计算道路亮度分布的方法一实施 例的R系列垂直道路方向的反射系数示意图。
[0029] 图7显示为本发明用直角坐标路面下的反射表计算道路亮度分布的方法一实施例 极坐标下的强度分布图。
[0030] 图8显示为本发明用直角坐标下的路面反射表计算道路亮度分布的方法一实施例 转换后直角坐标下的强度分布图。
[0031] 图9显示为本发明用直角坐标下的路面反射表计算道路亮度分布的方法一实施例 单车道上单个灯具贡献的亮度分布图。
[0032] 图10显示为本发明用直角坐标下的路面反射表计算道路亮度分布的方法一实施 例单车道上单个灯间距内的由五个灯具引起的各自的的亮度贡献图。
[0033] 图11显示为本发明用直角坐标下的路面反射表计算道路亮度分布的方法一实施 例单车道上单个灯间距内的亮度分布图。
[0034] 图12显示为本发明用直角坐标下的路面反射表计算道路亮度分布的方法一实施 例整个单车道的亮度分布图。
[0035] 图13显示为本发明用直角坐标下的路面反射表计算道路亮度分布的方法一实施 例整个车道的3D亮度分布图。
【具体实施方式】
[0036] W下通过特定的具体实例说明本发明的实施方式，本领域技术人 员可由本说明书 所掲露的内容轻易地了解本发明的其他优点与功效。本发明还可w通过另外不同的具体实 施方式加 W实施或应用，本说明书中的各项细节也可W基于不同观点与应用，在没有背离 本发明的精神下进行各种修饰或改变。
[0037] 请参阅图2至图13。需要说明的是，本实施例中所提供的图示仅W示意方式说明本 发明的基本构想，遂图式中仅显示与本发明中有关的组件而非按照实际实施时的组件数 目、形状及尺寸绘制，其实际实施时各组件的型态、数量及比例可为一种随意的改变，且其 组件布局型态也可能更为复杂。
[0038] 如图2所示，本发明提供了一种采用直角坐标路面反射表计算道路亮度分布的方 法，包括步骤:S1、将球坐标下的路面反射表Γ(β，丫）转换为直角坐标下的路面反射表Hx， y);S2、根据所述直角坐标路面反射表r(x，y)分析路面特性;S3、结合所述路面特性根据所 述直角坐标下的路面反射表r (X，y)计算道路亮度分布。坐标原点位于灯具的正下方，X方向 指向沿道路的方向，y方向则垂直于道路。
[0039] 进一步地，所述步骤S1包括步骤:S11、对所述球坐标下的路面反射表Γ(β，丫）进行 内插处理;S12、设定所述直角坐标下的路面反射表r(x，y)中相应数据；S13、将道路中每一 车道中有贡献的所述球坐标路面反射表Γ(β，γ )的反射系数相加得到对应点新的所述直角 坐标下的路面反射表r(x，y)的反射系数。在本发明一优选实施例中，所述相应数据包括路 面类型、灯的高度和灯的间距。
[0040] 为了在处理中不损失精度，我们首先对Ηβ，丫）表进行内插处理，本实施例中，用 五个节点的线性内插就可W保证精度。运样就得到了一个步长更小的更为精确的Ηβ，丫） 表，具体结果不在此给出。然而，此表还是we和丫为自变量，并非我们所希望的表。然后， 给出一个自变量为直角坐标的X和y的r(x，y)表，其X和y的范围及步长就如表2所示。表2是 路面类型为R3和10米灯高为例的r(x，y)表。
[0041 ]表2R3路面H=10m的r(x，y)表
[0042]
[0043] 进一步地，所述步骤S2还包括：不同的所述道路类型对应不同的所述球坐标下的 路面反射表Γ(β，丫），并根据不同的所述球坐标下的路面反射表Γ(β，丫）转换为对应的所述 直角坐标下的路面反射表r(x，y)。由于此时，当路面类型、灯的高度、灯的间距都是已知，表 中每格的X和y值也都固定，因此每点的0和丫也已知，只要对每一车道中对有贡献的灯具在 该(X，y)点的r (β，丫）中的反射系数相加就是被观察点新的（X，y)点的反射系数。路面的类 型和灯具的高度都是有限的几种，因此，需要计算的r(x，y)也是有限的几个。
[0044] 进一步地，如图3所示，所述步骤S2还包括步骤:根据所述直角坐标下的路面反射 表r(x，y)画出路面的所述直角坐标下的路面反射表r(x，y)的3D曲线图。
[0045] 本实施例中，直角坐标下的路面反射表r(x，y)是在笛卡尔系统中的，运使得我们 可W比较容易的画出直角坐标下的路面反射表Hx，y)的3D图，而运在W前是比较困难的 事。图3给出了 R3路面直角坐标下的路面反射表r(x，y)的3D图形。图中的X和y就是道路的坐 标，从图中可W清楚地看到直线道路的反射系数分布。运一个图包含了道路的整个反射性 质，可W初步体现直角坐标下的路面反射表r(x，y)的优越性。
[0046] 进一步地，如图4所示，所述步骤S2还包括步骤:取所述3D曲线图中对应车道的2D 曲线图，并根据所述2D曲线图分析对应车道在车道延伸方向前后的灯光衰减情况。图4显示 为图3中y = 0时的2D曲线图。由此可见，图4包含了计算y = 0那个车道上亮度分布所需要的 所有数据。图4中曲线的峰值不在x = 0而是位于x=12附近。而且曲线向右衰减得快，向左衰 减得慢。运表示对同样的灯具同样的距离时，前向光(x>〇的灯发出的光)比后向光(x<〇的 灯发出的光)反射得更多。运也意味着计算道路亮度分布时，相等后方的灯而言，前方较远 处的灯也有贡献。对R3路面30m的灯间距，由分析知道，计算路面某点的亮度时，一般取前向 4个灯后向1个灯即可。在y>0的其他车道上，与图4相比，曲线的数值不同但其趋向类似。
[0047] 上述实施例均WR系列路面的R3路面为例。在本发明中，不同材料路面的直角坐标 下的路面反射表r(x，y)可W按照上述的方法得到。例如对R系列的路面包含四个路面类型： R1，R2，R3和R4。图5表示了 R系列的四种路面在r (X，0)时的变化情况。如图5所示，不需数字 分析从图5我们就可W得到有价值的结论，即路面材料对反射系数影响非常大，R1和R2类路 面的方向性不强（前向和后向光的反射差别不大），而R3或R4路面有很强的方向性。有了直 角坐标下的路面反射表r(x，y)后，垂直于道路方向上反射系数的变化规律也易于得到。如 图6所示，图6给出了 R路面系列在垂直于道路方向上r(12，y)的变化规律。四条曲线随y的变 化趋势类似，但R1的数值较大，运给垂直于道路方向上的亮度分析提供了基础。
[004引进一步地，所述步骤S3还包括步骤:S31、计算单个灯的道路亮度分布;S32、计算单 个车道道路亮度分布;S33、计算整个道路道路亮度分布。
[0049] 其中，所述步骤S31还包括:根据公式：e = 将灯的极坐标下的强度分布I (9)转换为灯的直角坐标下的强度分布I(x);其中Θ为光线与灯具法线的夹角，Η为灯的高 度。
[0050] 用直角坐标下的路面反射表r(x，y)可W使得直线道路上的亮度分布的计算大大 简化。如图7所示，图7给出了我们设计的一个灯具极坐标下的Ι(θ)分布的例子，不难给出I (9)的数学表达式，它可W由图形模拟得到，也可W由描述灯具光强分布的IES文件得到，在 此我们将不给出其具体数学表达式。有了 Ι(θ)分布，就可W用公式0 = 十算亮度分 布。计算后的光强分布I(x)如图8所示。
[0051] 进一步地，所述步骤S32还包括:单车道的亮度正比于反射系数和光强的乘积，路 面的高度y = 〇;所述单车道的亮度分布的计算公式为:L(x)= ΣΟ · (;r(x，0) · Ι(χ))/Η2ηΗ 为灯具安装高度，求和则包括所有能射到该点的灯具，对确定道路来说C是常数，由于本文 只关屯、相对亮度，我们不必给定C值。
[0052] 由于亮度是正比于反射系数和光强的乘积，即路面上亮度为r(x，0)I(x)，如图9所 示。运就是单个灯具在单车道上的亮度分布。运里约定坐标X为灯具离开观察点的距离。由 图9可W看到前向光(x>0)的贡献明显大于后向光(x<0)。运是一个由于引入直角坐标下 的路面反射表r(x，y)而轻易知道的一个灯具在直线道路上引起的亮度是如何分布的，运是 十分有价值的结果，下面我们将看到，有了它我们就不难得到所有灯在路面所有点的亮度 分布了。
[0053] 设灯间距为30m，则由图9可W估算得正向第4个灯（120m远)对观察点的亮度贡献 只是最大贡献的2%，若灯间距缩小为25m，则第4个灯的贡献也只有最大贡献的3%。运样我 们可W得到在道路亮度计算中的重要结论:只需考虑正向4个灯方向1个灯就有足够的精度 了。考虑到R3路面是反射方向性较强的路面，而且y^O处的反射系数更小，因此，可W放屯、 地把运一结论用在其他路面和其他车道其他灯间距。
[0054] 对于一个较长的直线道路，每两个灯具间的亮度分布都是相同的，我们只需考虑 一个灯间距的范围就可W了，上面已经说明，在每个点只需考虑5个灯具。现在我们根据图 9,就可W更清楚用图形来表现运5个灯具影响的规律。W灯间距为30m为例，在图10中分别 给出 了运5个灯在两个灯具间即30m范围内所有观察点的亮度贡献。把图10中的5条线相加 就是一个灯间距间的亮度分布了，结果显示在图11中。图11中显示此灯具单车道的亮度均 匀性约为91 %。
[0055] 整个车道约240m范围内的亮度分布见图12,其实图12就是图11的8次周期性重复。 在图12中可W看到亮度周期性的波动，运种波动是间距式灯具照明不可避免的，但图中 10%的波动是比较小的。如果有严重的波动，就会形成影响交通安全的所谓"斑马线"，而运 种斑马线是可W在灯具设计时计算出来得到数字的结果，用于进一步分析，具体方法就不 在本发明中描述了。由此可见，有了直角坐标下的反射表r(x，y)，我们就可W得到如图9所 示的单灯在单车道上的亮度分布。再把此分布分成5段，我们就得到了如图10所示的5个有 贡献的灯具在一个点间距内的亮度。把此5条曲线向加就得到了如图11所示的一个灯间距 内单车道的亮度分布，它就是整个车道的亮度分布。W上运一图形化表现的过程清晰连贯。
[0056] Ι(θ)和Ι(χ)描述的是单车道的强度，对多车道的道路，灯具的光强分布可改为用 两个自变量的函数I(x，y)。运样公式应改为:L(x，y)= Σ0 · (r(x，y) · I(x，y))/H2。如何得 到I(x，y)呢?我们已经有了垂直于路面的反射系数分布，R3路面的结果见图6。我们可W用 此分布来作为由Ι(χ)扩展为I(x，y)的依据。简单地说，由于r(x，y)比r(x，0)小，如果希望道 路的亮度L在整个道路上都均匀，只有I(x，y)比l(x，0)大才行。作为初级近似，我们可W用
下式求得I(x，y): 运样，就可W把上面单车道的亮度计算结果扩展到 0 多车道。用简单的程序不难得计算到整个道路上的亮度分布。仍然WR3路面为例，结果显示 在图13中，在此例中整个亮度均匀性约为85%。
[0057] 综上，本发明把CIE的反射表；Γ(β，丫）改造成直角坐标下的;r(x，y)表，结果由于和 路面坐标一致，运一改造的效果使得道路亮度计算十分便利。
[0058] 现有技术中某些亮度分布可W用商用软件来进行，但商用软件局限性较大，没有 数字结果，也不能根据用户需求而变化。相对于商用软件，本文对亮度的分析更为形象、系 统和深入，还可W得到更实用的图形化的结果，而且本文的量化的计算方法是普通设计人 员都不难掌握的。
[0059] 本文的前提是直线道路，但由于实际上绝大多数道路都是直线的，因此本文的结 果的应用范围的非常广的。
[0060] 对不同的灯高，应该有不同的r(x，y)表，但由于实际上的灯具的高度只有有限的 几种，如8m，10m，由于本文的方法是一般光学工程设计人员都易于掌握的，可W用编程来更 简化计算过程，因此运也没有对设计带来很大的影响。
[0061] 由于计算道路的亮度变得容易了，运样就可W方便地修改LED的光强分布从而得 到设计者要求的亮度分布。而有了光强分布则可W计算透镜的形状。运将大大加快等亮度 分布的照明灯具的设计进程。运将是更为重要的应用，但有关内容已不在本文的范畴之内。
[0062] 上述实施例仅例示性说明本发明的原理及其功效，而非用于限制本发明。任何熟 悉此技术的人±皆可在不违背本发明的精神及范畴下，对上述实施例进行修饰或改变。因 此，举凡所属技术领域中具有通常知识者在未脱离本发明所掲示的精神与技术思想下所完 成的一切等效修饰或改变，仍应由本发明的权利要求所涵盖。
【主权项】
1. 一种采用直角坐标下的路面反射表计算道路亮度分布的方法，其特征在于，包括步 骤： 51、 将球坐标下的路面反射表Γ(β，γ )转换为直角坐标下的路面反射表r(X，y); 52、 根据所述直角坐标下的路面反射表r (X，y)分析路面特性； 53、 结合所述路面特性根据所述直角坐标下的路面反射表r(x，y)计算道路亮度分布。2. 根据权利要求1所述的采用直角坐标下的路面反射表计算道路亮度分布的方法，其 特征在于，所述步骤Sl包括步骤： 511、 对所述球坐标下的路面反射表r(i3, γ )进行内插处理； 512、 计算所述直角坐标下的路面反射表r (X，y)中相应数据； 513、 将道路中每一车道中有亮度贡献的灯具的所述球坐标下的路面反射表Γ(β，γ )的 反射系数相加得到对应点新的所述直角坐标下的路面反射表r(x，y)的反射系数。3. 根据权利要求2所述的采用直角坐标下的路面反射表计算道路亮度分布的方法，其 特征在于，所述有亮度贡献的灯具的贡献值至少为2%以上;所述相应数据包括路面类型、 灯具的高度和灯具的间距。4. 根据权利要求3所述的采用直角坐标下的路面反射表计算道路亮度分布的方法，其 特征在于，所述步骤S2还包括： 不同的所述道路类型对应不同的所述球坐标下的路面反射表r(i3, γ )，并根据不同的 所述球坐标下的路面反射表r(i3, γ )转换为对应的所述直角坐标下的路面反射表r(x，y)。5. 根据权利要求1所述的采用直角坐标下的路面反射表计算道路亮度分布的方法，其 特征在于，所述步骤S2还包括步骤： 根据所述直角坐标下的路面反射表r(x，y)画出路面的所述直角坐标下的路面反射表r (x，y)的3D曲线图。6. 根据权利要求5所述的采用直角坐标路面下的反射表计算道路亮度分布的方法，其 特征在于，所述步骤S2还包括步骤： 取所述3D曲线图中对应车道的2D曲线图，并根据所述2D曲线图分析对应车道在车道延 伸方向前后的灯光亮度衰减情况。7. 根据权利要求1所述的采用直角坐标下的路面反射表计算道路亮度分布的方法，其 特征在于，所述步骤S3还包括步骤： 531、 计算单个灯具对道路亮度分布的贡献； 532、 计算单个车道的道路亮度分布； 533、 计算整个道路的道路亮度分布。8. 根据权利要求7所述的采用直角坐标下的路面反射表计算道路亮度分布的方法，其 特征在于，所述步骤S31还包括： 根据公式：将灯的极坐标下的强度分布Ι(θ)转换为灯的直角坐标下的强度 分布Ι(χ); 其中Θ为光线与灯具法线反向的夹角，H为灯具安装高度。9. 根据权利要求8所述的采用直角坐标下的路面反射表r(x，y)计算道路亮度分布的方 法，其特征在于，所述步骤S32还包括： 单车道的亮度正比于反射系数和光强的乘积； 所述单车道的亮度分布的计算公式为:L(x)= SC · (r(x，0) · I(x))/H2。10.根据权利要求9所述的采用直角坐标下的路面反射表计算道路亮度分布的方法，其 特征在于，所述步骤S33还包括： 所述整个道路的亮度分布的计算公式为:L(x，y)= SC · (r(x，y) · I(x，y))/H2; 其中，作为初级近似计算，C为常数。
【专利摘要】本发明提供了一种采用直角坐标下的路面反射表计算道路亮度分布的方法，包括步骤：S1、将球坐标下的路面反射表r(β，γ)转换为直角坐标下的路面反射表r(x，y)；S2、根据所述直角坐标下的路面反射表r(x，y)分析路面特性；S3、结合所述路面特性根据所述直角坐标下的路面反射表r(x，y)计算道路亮度分布。本发明提供的一种采用直角坐标下的路面反射表计算道路亮度分布的方法，方便地计算单个车道和整个道路的亮度分布，大大有利于分析给定灯具时道路的亮度分布，以及设计等亮度分布的灯具等。
【IPC分类】G06F17/50
【公开号】CN105488249
【申请号】CN201510799893
【发明人】周士康, 陈春根, 许礼, 王鹰华
【申请人】上海三思电子工程有限公司, 上海三思科技发展有限公司, 嘉善三思光电技术有限公司, 三思光电科技（上海）有限公司
【公开日】2016年4月13日
【申请日】2015年11月19日
【技术领域】
[0001] 本发明设及道路照明领域，特别是设及一种采用直角坐标下的路面反射表计算道 路亮度分布的方法。
【背景技术】
[0002] 亮度评价方法在道路照明评价中变得越来越重要，但有关亮度的计算工作却十分 缺乏。虽然用CIE(国际照明委员会)提供的路面反射表Ηβ，丫）可W计算道路的亮度分布， 但其实际的计算过程却很困难。
[0003] 在道路照明的评价方法中，亮度评价比照度评价更为恰当，特别是对于高速公路 的照明更是如此。运是由于亮度依赖于司机的观察，它与道路交通事故是紧密相关的，而照 度只依赖仪器的测量。许多国家已经把亮度评价作为道路照明质量评价的主要方法，亮度 的计算在道路照明中变得越来越重要了。
[0004] 与传统灯具相比，L抓可W比较容易地用各种透镜来改变其光强分布，因此L抓道 路照明可W实现精确的光强分布，从而实现较高的道路照明亮度均匀性，而同时又能够满 足一定的照度均匀性，W满足某些道路照明标准的要求。运种现实的可能性大大地增加了 道路照明工程对首先满足均匀亮度分布而不是均匀照度分布的需求，运也提高了对Lm)灯 具的道路照明亮度分布计算的要求。但由于其复杂性，运种计算在W前是技术上的弱项，关 于照明亮度分布的分析还只能依赖与市售的照明软件，而运些软件的分析只能到有限的结 果，而且是不能按设计者的要求进行分析。迄今为止道路照明方面亮度分布计算的文章还 极少。
[0005] 道路的亮度与许多因素相关，包括不同的路面类型，光线的方向，观察者的位置， 观察的方向等。为了使得亮度分布的计算规范化，CIE(国际照明委员会)和PIARC(国际道路 会议常设协会)规定了标准路面和标准观察条件，运使得亮度计算有了标准。
[0006] 如图1所示，β为入射面和观察面的夹角，丫是灯具垂线和光线的夹角。简化反射系 数表也称为反射表，在此用Γ(β，丫）表示，下列的表1就是CIE出版的R3路面为例的Γ(β，丫） 表。
[0007] 表1CIE出版的R3路面的Γ(β，丫）表 [000引
[0009]
[0010] 用CIE提供的Γ(β，丫）表就可w用来计算路面亮度了。但是用Γ(β，丫）去计算通常 是直角坐标下的道路时是很困难的。而且，考虑到r(e，丫）给出的只是道路上某一个点由一 个灯具引起的亮度贡献，一个车道的亮度分布还必须考虑许每一个观察点上各灯的贡献之 和，加之道路还不止一个车道，可W想象计算将十分困难。
【发明内容】
[0011] 鉴于W上所述现有技术的缺点，本发明的目的在于提供一种直角坐标下的路面反 射表计算道路亮度分布的方法，用于解决现有技术中由于道路亮度分布是基于直角坐标而 CIE的反射表却基于球坐标，从而计算十分困难的问题。
[0012] 本发明提供的一种采用直角坐标下的路面反射表计算道路亮度分布的方法，包括 步骤:S1、将球坐标下的路面反射表Γ(β，丫）转换为直角坐标下的路面反射表r(x，y);S2、根 据所述直角坐标下的路面反射表Hx，y)分析路面特性;S3、结合所述路面特性根据所述直 角坐标下的路面反射表r(x，y)计算道路亮度分布。
[0013] 于本发明的一实施方式中，所述步骤S1包括步骤:S11、对所述球坐标下的路面反 射表r(0，丫）进行内插处理；S12、计算所述直角坐标下的路面反射表;r(x，y)中相应数据； S13、将道路中每一车道中有亮度贡献的灯具的所述球坐标下的路面反射表r(i3, 丫）的反射 系数相加得到对应点新的所述直角坐标下的路面反射表r(x，y)的反射系数。
[0014] 于本发明的一实施方式中，所述有亮度贡献的灯具的贡献值至少为2% W上;所述 相应数据包括路面类型、灯具的高度和灯具的间距。
[0015] 于本发明的一实施方式中，所述步骤S2还包括:不同的所述道路类型对应不同的 所述球坐标下的路面反射表r(e，丫），并根据不同的所述球坐标下的路面反射表Γ(β，丫）转 换为对应的所述直角坐标下的路面反射表r (X，y)。
[0016] 于本发明的一实施方式中，所述步骤S2还包括步骤:根据所述直角坐标下的路面 反射表r(x，y)画出路面的所述直角坐标下的路面反射表r(x，y)的3D曲线图。
[0017] 于本发明的一实施方式中，所述步骤S2还包括步骤:取所述3D曲线图中对应车道 的2D曲线图，并根据所述2D曲线图分析对应车道在车道延伸方向前后的灯光亮度衰减情 况。
[0018] 于本发明的一实施方式中，所述步骤S3还包括步骤:S31、计算单个灯具对道路亮 度分布的贡献;S32、计算单个车道的道路亮度分布;S33、计算整个道路的道路亮度分布。
[0019] 于本发明的一实施方式中，所述步骤S31还包括:根据公式：8:::燃邮音，将灯的极 坐标下的强度分布Ι(θ)转换为灯的直角坐标下的强度分布I(x);其中θ为光线与灯具法线 反向的夹角，Η为灯具安装高度。
[0020] 于本发明的一实施方式中，所述步骤S32还包括:单车道的亮度正比于反射系数和 光强的乘积;所述单车道的亮度分布的计算公式为:L(x)= Σ0 · (r(x，0) · Ι(χ))/Η2。
[0021] 于本发明的一实施方式中，所述步骤S33还包括:所述整个道路的亮度分布的计算 公式为：L(x，y)=IC· (;r(x，y) · I(x，y))/H2;其中，作为初级近似计算， 乃―= 户，C为常数。
[0022] 如上所述，本发明具有W下有益效果:方便地计算单个车道和整个道路的亮度分 布，大大有利于分析给定灯具时道路的亮度如何分布，W及如何设计等亮度分布的灯具等。
【附图说明】
[0023] 图1显示为路面反射的角度关系示意图。
[0024] 图2显示为本发明采用直角坐标下的路面反射表计算道路亮度分布的方法一实施 例的流程方框图。
[0025] 图3显示为本发明采用直角坐标下的路面反射表计算道路亮度分布的方法一实施 例的路面r(x，y)表的3D曲线图。
[0026] 图4显示为图3中y = 0时的2D曲线图。
[0027] 图5显示为本发明采用直角坐标下的路面反射表计算道路亮度分布的方法一实施 例的R系列沿路方向的反射系数示意图。
[0028] 图6显示为本发明采用直角坐标路面下的反射表计算道路亮度分布的方法一实施 例的R系列垂直道路方向的反射系数示意图。
[0029] 图7显示为本发明用直角坐标路面下的反射表计算道路亮度分布的方法一实施例 极坐标下的强度分布图。
[0030] 图8显示为本发明用直角坐标下的路面反射表计算道路亮度分布的方法一实施例 转换后直角坐标下的强度分布图。
[0031] 图9显示为本发明用直角坐标下的路面反射表计算道路亮度分布的方法一实施例 单车道上单个灯具贡献的亮度分布图。
[0032] 图10显示为本发明用直角坐标下的路面反射表计算道路亮度分布的方法一实施 例单车道上单个灯间距内的由五个灯具引起的各自的的亮度贡献图。
[0033] 图11显示为本发明用直角坐标下的路面反射表计算道路亮度分布的方法一实施 例单车道上单个灯间距内的亮度分布图。
[0034] 图12显示为本发明用直角坐标下的路面反射表计算道路亮度分布的方法一实施 例整个单车道的亮度分布图。
[0035] 图13显示为本发明用直角坐标下的路面反射表计算道路亮度分布的方法一实施 例整个车道的3D亮度分布图。
【具体实施方式】
[0036] W下通过特定的具体实例说明本发明的实施方式，本领域技术人 员可由本说明书 所掲露的内容轻易地了解本发明的其他优点与功效。本发明还可w通过另外不同的具体实 施方式加 W实施或应用，本说明书中的各项细节也可W基于不同观点与应用，在没有背离 本发明的精神下进行各种修饰或改变。
[0037] 请参阅图2至图13。需要说明的是，本实施例中所提供的图示仅W示意方式说明本 发明的基本构想，遂图式中仅显示与本发明中有关的组件而非按照实际实施时的组件数 目、形状及尺寸绘制，其实际实施时各组件的型态、数量及比例可为一种随意的改变，且其 组件布局型态也可能更为复杂。
[0038] 如图2所示，本发明提供了一种采用直角坐标路面反射表计算道路亮度分布的方 法，包括步骤:S1、将球坐标下的路面反射表Γ(β，丫）转换为直角坐标下的路面反射表Hx， y);S2、根据所述直角坐标路面反射表r(x，y)分析路面特性;S3、结合所述路面特性根据所 述直角坐标下的路面反射表r (X，y)计算道路亮度分布。坐标原点位于灯具的正下方，X方向 指向沿道路的方向，y方向则垂直于道路。
[0039] 进一步地，所述步骤S1包括步骤:S11、对所述球坐标下的路面反射表Γ(β，丫）进行 内插处理;S12、设定所述直角坐标下的路面反射表r(x，y)中相应数据；S13、将道路中每一 车道中有贡献的所述球坐标路面反射表Γ(β，γ )的反射系数相加得到对应点新的所述直角 坐标下的路面反射表r(x，y)的反射系数。在本发明一优选实施例中，所述相应数据包括路 面类型、灯的高度和灯的间距。
[0040] 为了在处理中不损失精度，我们首先对Ηβ，丫）表进行内插处理，本实施例中，用 五个节点的线性内插就可W保证精度。运样就得到了一个步长更小的更为精确的Ηβ，丫） 表，具体结果不在此给出。然而，此表还是we和丫为自变量，并非我们所希望的表。然后， 给出一个自变量为直角坐标的X和y的r(x，y)表，其X和y的范围及步长就如表2所示。表2是 路面类型为R3和10米灯高为例的r(x，y)表。
[0041 ]表2R3路面H=10m的r(x，y)表
[0042]
[0043] 进一步地，所述步骤S2还包括：不同的所述道路类型对应不同的所述球坐标下的 路面反射表Γ(β，丫），并根据不同的所述球坐标下的路面反射表Γ(β，丫）转换为对应的所述 直角坐标下的路面反射表r(x，y)。由于此时，当路面类型、灯的高度、灯的间距都是已知，表 中每格的X和y值也都固定，因此每点的0和丫也已知，只要对每一车道中对有贡献的灯具在 该(X，y)点的r (β，丫）中的反射系数相加就是被观察点新的（X，y)点的反射系数。路面的类 型和灯具的高度都是有限的几种，因此，需要计算的r(x，y)也是有限的几个。
[0044] 进一步地，如图3所示，所述步骤S2还包括步骤:根据所述直角坐标下的路面反射 表r(x，y)画出路面的所述直角坐标下的路面反射表r(x，y)的3D曲线图。
[0045] 本实施例中，直角坐标下的路面反射表r(x，y)是在笛卡尔系统中的，运使得我们 可W比较容易的画出直角坐标下的路面反射表Hx，y)的3D图，而运在W前是比较困难的 事。图3给出了 R3路面直角坐标下的路面反射表r(x，y)的3D图形。图中的X和y就是道路的坐 标，从图中可W清楚地看到直线道路的反射系数分布。运一个图包含了道路的整个反射性 质，可W初步体现直角坐标下的路面反射表r(x，y)的优越性。
[0046] 进一步地，如图4所示，所述步骤S2还包括步骤:取所述3D曲线图中对应车道的2D 曲线图，并根据所述2D曲线图分析对应车道在车道延伸方向前后的灯光衰减情况。图4显示 为图3中y = 0时的2D曲线图。由此可见，图4包含了计算y = 0那个车道上亮度分布所需要的 所有数据。图4中曲线的峰值不在x = 0而是位于x=12附近。而且曲线向右衰减得快，向左衰 减得慢。运表示对同样的灯具同样的距离时，前向光(x>〇的灯发出的光)比后向光(x<〇的 灯发出的光)反射得更多。运也意味着计算道路亮度分布时，相等后方的灯而言，前方较远 处的灯也有贡献。对R3路面30m的灯间距，由分析知道，计算路面某点的亮度时，一般取前向 4个灯后向1个灯即可。在y>0的其他车道上，与图4相比，曲线的数值不同但其趋向类似。
[0047] 上述实施例均WR系列路面的R3路面为例。在本发明中，不同材料路面的直角坐标 下的路面反射表r(x，y)可W按照上述的方法得到。例如对R系列的路面包含四个路面类型： R1，R2，R3和R4。图5表示了 R系列的四种路面在r (X，0)时的变化情况。如图5所示，不需数字 分析从图5我们就可W得到有价值的结论，即路面材料对反射系数影响非常大，R1和R2类路 面的方向性不强（前向和后向光的反射差别不大），而R3或R4路面有很强的方向性。有了直 角坐标下的路面反射表r(x，y)后，垂直于道路方向上反射系数的变化规律也易于得到。如 图6所示，图6给出了 R路面系列在垂直于道路方向上r(12，y)的变化规律。四条曲线随y的变 化趋势类似，但R1的数值较大，运给垂直于道路方向上的亮度分析提供了基础。
[004引进一步地，所述步骤S3还包括步骤:S31、计算单个灯的道路亮度分布;S32、计算单 个车道道路亮度分布;S33、计算整个道路道路亮度分布。
[0049] 其中，所述步骤S31还包括:根据公式：e = 将灯的极坐标下的强度分布I (9)转换为灯的直角坐标下的强度分布I(x);其中Θ为光线与灯具法线的夹角，Η为灯的高 度。
[0050] 用直角坐标下的路面反射表r(x，y)可W使得直线道路上的亮度分布的计算大大 简化。如图7所示，图7给出了我们设计的一个灯具极坐标下的Ι(θ)分布的例子，不难给出I (9)的数学表达式，它可W由图形模拟得到，也可W由描述灯具光强分布的IES文件得到，在 此我们将不给出其具体数学表达式。有了 Ι(θ)分布，就可W用公式0 = 十算亮度分 布。计算后的光强分布I(x)如图8所示。
[0051] 进一步地，所述步骤S32还包括:单车道的亮度正比于反射系数和光强的乘积，路 面的高度y = 〇;所述单车道的亮度分布的计算公式为:L(x)= ΣΟ · (;r(x，0) · Ι(χ))/Η2ηΗ 为灯具安装高度，求和则包括所有能射到该点的灯具，对确定道路来说C是常数，由于本文 只关屯、相对亮度，我们不必给定C值。
[0052] 由于亮度是正比于反射系数和光强的乘积，即路面上亮度为r(x，0)I(x)，如图9所 示。运就是单个灯具在单车道上的亮度分布。运里约定坐标X为灯具离开观察点的距离。由 图9可W看到前向光(x>0)的贡献明显大于后向光(x<0)。运是一个由于引入直角坐标下 的路面反射表r(x，y)而轻易知道的一个灯具在直线道路上引起的亮度是如何分布的，运是 十分有价值的结果，下面我们将看到，有了它我们就不难得到所有灯在路面所有点的亮度 分布了。
[0053] 设灯间距为30m，则由图9可W估算得正向第4个灯（120m远)对观察点的亮度贡献 只是最大贡献的2%，若灯间距缩小为25m，则第4个灯的贡献也只有最大贡献的3%。运样我 们可W得到在道路亮度计算中的重要结论:只需考虑正向4个灯方向1个灯就有足够的精度 了。考虑到R3路面是反射方向性较强的路面，而且y^O处的反射系数更小，因此，可W放屯、 地把运一结论用在其他路面和其他车道其他灯间距。
[0054] 对于一个较长的直线道路，每两个灯具间的亮度分布都是相同的，我们只需考虑 一个灯间距的范围就可W了，上面已经说明，在每个点只需考虑5个灯具。现在我们根据图 9,就可W更清楚用图形来表现运5个灯具影响的规律。W灯间距为30m为例，在图10中分别 给出 了运5个灯在两个灯具间即30m范围内所有观察点的亮度贡献。把图10中的5条线相加 就是一个灯间距间的亮度分布了，结果显示在图11中。图11中显示此灯具单车道的亮度均 匀性约为91 %。
[0055] 整个车道约240m范围内的亮度分布见图12,其实图12就是图11的8次周期性重复。 在图12中可W看到亮度周期性的波动，运种波动是间距式灯具照明不可避免的，但图中 10%的波动是比较小的。如果有严重的波动，就会形成影响交通安全的所谓"斑马线"，而运 种斑马线是可W在灯具设计时计算出来得到数字的结果，用于进一步分析，具体方法就不 在本发明中描述了。由此可见，有了直角坐标下的反射表r(x，y)，我们就可W得到如图9所 示的单灯在单车道上的亮度分布。再把此分布分成5段，我们就得到了如图10所示的5个有 贡献的灯具在一个点间距内的亮度。把此5条曲线向加就得到了如图11所示的一个灯间距 内单车道的亮度分布，它就是整个车道的亮度分布。W上运一图形化表现的过程清晰连贯。
[0056] Ι(θ)和Ι(χ)描述的是单车道的强度，对多车道的道路，灯具的光强分布可改为用 两个自变量的函数I(x，y)。运样公式应改为:L(x，y)= Σ0 · (r(x，y) · I(x，y))/H2。如何得 到I(x，y)呢?我们已经有了垂直于路面的反射系数分布，R3路面的结果见图6。我们可W用 此分布来作为由Ι(χ)扩展为I(x，y)的依据。简单地说，由于r(x，y)比r(x，0)小，如果希望道 路的亮度L在整个道路上都均匀，只有I(x，y)比l(x，0)大才行。作为初级近似，我们可W用
下式求得I(x，y): 运样，就可W把上面单车道的亮度计算结果扩展到 0 多车道。用简单的程序不难得计算到整个道路上的亮度分布。仍然WR3路面为例，结果显示 在图13中，在此例中整个亮度均匀性约为85%。
[0057] 综上，本发明把CIE的反射表；Γ(β，丫）改造成直角坐标下的;r(x，y)表，结果由于和 路面坐标一致，运一改造的效果使得道路亮度计算十分便利。
[0058] 现有技术中某些亮度分布可W用商用软件来进行，但商用软件局限性较大，没有 数字结果，也不能根据用户需求而变化。相对于商用软件，本文对亮度的分析更为形象、系 统和深入，还可W得到更实用的图形化的结果，而且本文的量化的计算方法是普通设计人 员都不难掌握的。
[0059] 本文的前提是直线道路，但由于实际上绝大多数道路都是直线的，因此本文的结 果的应用范围的非常广的。
[0060] 对不同的灯高，应该有不同的r(x，y)表，但由于实际上的灯具的高度只有有限的 几种，如8m，10m，由于本文的方法是一般光学工程设计人员都易于掌握的，可W用编程来更 简化计算过程，因此运也没有对设计带来很大的影响。
[0061] 由于计算道路的亮度变得容易了，运样就可W方便地修改LED的光强分布从而得 到设计者要求的亮度分布。而有了光强分布则可W计算透镜的形状。运将大大加快等亮度 分布的照明灯具的设计进程。运将是更为重要的应用，但有关内容已不在本文的范畴之内。
[0062] 上述实施例仅例示性说明本发明的原理及其功效，而非用于限制本发明。任何熟 悉此技术的人±皆可在不违背本发明的精神及范畴下，对上述实施例进行修饰或改变。因 此，举凡所属技术领域中具有通常知识者在未脱离本发明所掲示的精神与技术思想下所完 成的一切等效修饰或改变，仍应由本发明的权利要求所涵盖。
【主权项】
1. 一种采用直角坐标下的路面反射表计算道路亮度分布的方法，其特征在于，包括步 骤： 51、 将球坐标下的路面反射表Γ(β，γ )转换为直角坐标下的路面反射表r(X，y); 52、 根据所述直角坐标下的路面反射表r (X，y)分析路面特性； 53、 结合所述路面特性根据所述直角坐标下的路面反射表r(x，y)计算道路亮度分布。2. 根据权利要求1所述的采用直角坐标下的路面反射表计算道路亮度分布的方法，其 特征在于，所述步骤Sl包括步骤： 511、 对所述球坐标下的路面反射表r(i3, γ )进行内插处理； 512、 计算所述直角坐标下的路面反射表r (X，y)中相应数据； 513、 将道路中每一车道中有亮度贡献的灯具的所述球坐标下的路面反射表Γ(β，γ )的 反射系数相加得到对应点新的所述直角坐标下的路面反射表r(x，y)的反射系数。3. 根据权利要求2所述的采用直角坐标下的路面反射表计算道路亮度分布的方法，其 特征在于，所述有亮度贡献的灯具的贡献值至少为2%以上;所述相应数据包括路面类型、 灯具的高度和灯具的间距。4. 根据权利要求3所述的采用直角坐标下的路面反射表计算道路亮度分布的方法，其 特征在于，所述步骤S2还包括： 不同的所述道路类型对应不同的所述球坐标下的路面反射表r(i3, γ )，并根据不同的 所述球坐标下的路面反射表r(i3, γ )转换为对应的所述直角坐标下的路面反射表r(x，y)。5. 根据权利要求1所述的采用直角坐标下的路面反射表计算道路亮度分布的方法，其 特征在于，所述步骤S2还包括步骤： 根据所述直角坐标下的路面反射表r(x，y)画出路面的所述直角坐标下的路面反射表r (x，y)的3D曲线图。6. 根据权利要求5所述的采用直角坐标路面下的反射表计算道路亮度分布的方法，其 特征在于，所述步骤S2还包括步骤： 取所述3D曲线图中对应车道的2D曲线图，并根据所述2D曲线图分析对应车道在车道延 伸方向前后的灯光亮度衰减情况。7. 根据权利要求1所述的采用直角坐标下的路面反射表计算道路亮度分布的方法，其 特征在于，所述步骤S3还包括步骤： 531、 计算单个灯具对道路亮度分布的贡献； 532、 计算单个车道的道路亮度分布； 533、 计算整个道路的道路亮度分布。8. 根据权利要求7所述的采用直角坐标下的路面反射表计算道路亮度分布的方法，其 特征在于，所述步骤S31还包括： 根据公式：将灯的极坐标下的强度分布Ι(θ)转换为灯的直角坐标下的强度 分布Ι(χ); 其中Θ为光线与灯具法线反向的夹角，H为灯具安装高度。9. 根据权利要求8所述的采用直角坐标下的路面反射表r(x，y)计算道路亮度分布的方 法，其特征在于，所述步骤S32还包括： 单车道的亮度正比于反射系数和光强的乘积； 所述单车道的亮度分布的计算公式为:L(x)= SC · (r(x，0) · I(x))/H2。10.根据权利要求9所述的采用直角坐标下的路面反射表计算道路亮度分布的方法，其 特征在于，所述步骤S33还包括： 所述整个道路的亮度分布的计算公式为:L(x，y)= SC · (r(x，y) · I(x，y))/H2; 其中，作为初级近似计算，C为常数。
【专利摘要】本发明提供了一种采用直角坐标下的路面反射表计算道路亮度分布的方法，包括步骤：S1、将球坐标下的路面反射表r(β，γ)转换为直角坐标下的路面反射表r(x，y)；S2、根据所述直角坐标下的路面反射表r(x，y)分析路面特性；S3、结合所述路面特性根据所述直角坐标下的路面反射表r(x，y)计算道路亮度分布。本发明提供的一种采用直角坐标下的路面反射表计算道路亮度分布的方法，方便地计算单个车道和整个道路的亮度分布，大大有利于分析给定灯具时道路的亮度分布，以及设计等亮度分布的灯具等。
【IPC分类】G06F17/50
【公开号】CN105488249
【申请号】CN201510799893
【发明人】周士康, 陈春根, 许礼, 王鹰华
【申请人】上海三思电子工程有限公司, 上海三思科技发展有限公司, 嘉善三思光电技术有限公司, 三思光电科技（上海）有限公司
【公开日】2016年4月13日
【申请日】2015年11月19日
最新回复(0)