[已解决]用zemax的ZPL语言测量激光光斑的均匀度，就像这样的探测器图片，可否指点一二，论...

泽深向山

我做了一个准直匀化的光学系统，想用zemax直接看看这个光斑的均匀度，我现在用的方法是，把这个图片输出来，用matlab写了一个程序，把它的平均值，最大值，最小值，均方根都求出来，然后代入到平均绝对误差的公式里面，来计算它的均匀度，但是老师说这样不能写进论文里面，因为数据变了，只能用作参考，所以想找找用ZPL编写宏的大神，能指点一二

最佳答案

月排行榜 / 总排行榜

帅三石

2022-6-7 22:20

这个是我之前在非序列里面编的宏，希望对你有所帮助。
temp=NSDD(1,0,0,0)
NSTR 1,0,0,0,0,0,0,0    # NSTR surf, source, split, scatter, usepol, ignore_err, rand_seed, save

EBAR=NSDD(1,5,-5,1)  #获取非零像素上的平均光照度
PRINT "光照度平均值=",EBAR

!DECLARE E,DOUBLE,1,320000
DECLARE A,DOUBLE,1,320000 #320000是探测器像素数

SUM=0
M=0
SUM1=0
FOR i=1,320000,1
  EI=NSDD(1,5,i,1)  # 获取单个像素的光照度
  HI=NSDD(1,5,i,0)  # 获取单个像素的能量

  IF(EI>0)
     AI=POWR((EI-EBAR),2)
    ! A(i)=AI
     SUM=SUM+AI
     M=M+1           # M是光线打到探测像素个数
     SUM1=SUM1+HI
  ENDIF

NEXT
PRINT "sum=",SUM
PRINT  "M=", M
PRINT "SUM1=",SUM1

R=1-(SQRT(SUM/M))/EBAR  # 求光斑的均匀性

跳转到最佳答案楼层

帅三石

这个是我之前在非序列里面编的宏，希望对你有所帮助。
temp=NSDD(1,0,0,0)
NSTR 1,0,0,0,0,0,0,0    # NSTR surf, source, split, scatter, usepol, ignore_err, rand_seed, save

EBAR=NSDD(1,5,-5,1)  #获取非零像素上的平均光照度
PRINT "光照度平均值=",EBAR

!DECLARE E,DOUBLE,1,320000
DECLARE A,DOUBLE,1,320000 #320000是探测器像素数

SUM=0
M=0
SUM1=0
FOR i=1,320000,1
  EI=NSDD(1,5,i,1)  # 获取单个像素的光照度
  HI=NSDD(1,5,i,0)  # 获取单个像素的能量

  IF(EI>0)
     AI=POWR((EI-EBAR),2)
    ! A(i)=AI
     SUM=SUM+AI
     M=M+1           # M是光线打到探测像素个数
     SUM1=SUM1+HI
  ENDIF

NEXT
PRINT "sum=",SUM
PRINT  "M=", M
PRINT "SUM1=",SUM1

R=1-(SQRT(SUM/M))/EBAR  # 求光斑的均匀性

泽深向山

帅三石 发表于 2022-6-7 22:20
这个是我之前在非序列里面编的宏，希望对你有所帮助。
temp=NSDD(1,0,0,0)
NSTR 1,0,0,0,0,0,0,0    # NS ...

感谢您的帮助！有个问题就是，我这个zemax仿真布置了四个探测器，然后怎么看这个宏运行出来的结果是哪张图？还有就是那个FOR的循环显示说语法有问题，我也不太懂这个zpl的语法，请问是否可以注释的更详细点

帅三石

泽深向山 发表于 2022-6-8 22:20
感谢您的帮助！有个问题就是，我这个zemax仿真布置了四个探测器，然后怎么看这个宏运行出来的结果是哪张 ...

关于你使用第几个探测器，请查看一下NSDD的用法，数字5代表的是探测器的序号，关于你说的FOR循环有问题，我这边能正常使用，可否把你的显示错误的截图发一下。（应该是提交时网卡了，请版主删掉重复的）

泽深向山

帅三石 发表于 2022-6-9 08:59
关于你使用第几个探测器，请查看一下NSDD的用法，数字5代表的是探测器的序号，关于你说的FOR循环有问题， ...

我探测器是500*500的，然后我就把320000改成了250000，然后把NSDD里面的5换成了相应的探测器，出现的问题就是这个

泽深向山

泽深向山 发表于 2022-6-9 10:07
我探测器是500*500的，然后我就把320000改成了250000，然后把NSDD里面的5换成了相应的探测器，出现的问 ...

请问NSDD的Def#是体现在哪里啊，NSDD(1,0,0,0),这四个哪个是？（实在对不起，我真的不太会）

帅三石

泽深向山 发表于 2022-6-9 10:14
请问NSDD的Def#是体现在哪里啊，NSDD(1,0,0,0),这四个哪个是？（实在对不起，我真的不太会）

这个是清除探测器数据，宏语言前两行，不用管。你检查一下宏语言标点符号是英文的吗，还有你前后变量一致吗，运行之前先进行一次光线追迹，我这边程序能正常运行

泽深向山

帅三石 发表于 2022-6-9 11:37
这个是清除探测器数据，宏语言前两行，不用管。你检查一下宏语言标点符号是英文的吗，还有你前后变量一致 ...

可以正常运行了，十分好用，感谢！！！

泽深向山

帅三石 发表于 2022-6-9 11:37
这个是清除探测器数据，宏语言前两行，不用管。你检查一下宏语言标点符号是英文的吗，还有你前后变量一致 ...

  想问一下是不是有一个print，没有打印出来？，平均光照度的那个

帅三石

泽深向山 发表于 2022-6-9 13:16
想问一下是不是有一个print，没有打印出来？，平均光照度的那个

对的，你这边好像没看到