像差的变化数据提取
文章目录
- 像差的变化数据提取
- 前言
- 一、主要目标
- 二、主要代码
- 1.VBA代码
- 2.CODE V macro
- 输出结果
前言
在优化过程中，相差的观察非常重要，尤其是镜片变多以后哪一面的曲率，厚度以及非球面系数的变化对像差的影响需要总结，提炼以便后续更好的优化，此数据提出非常重要。
一、主要目标
将三级像差输出至相应的单元格。
二、主要代码
1.VBA代码
代码如下（示例）：
Sub getTHirdabrr()'定义接口
Dim Session As CVCommand'102是CODEV的版本，即10.2
Set Session = CreateObject("CodeV.Command.102")
'设置打开文件的初始位置
Session.SetStartingDirectory ("c:\CVUSER")
Session.StartCodeV
'载入(dbgauss.lens)，在命令处输入res cv_lens:dbgauss
result = Session.Command("res cv_lens:dbgauss;")'提取三级像差至 abr.txt文件
result = Session.Command("in cv_macro:forder;")
result = Session.Command("buf del ba;buf yes;in cv_macro:thirdabrget 1 0 1;buf mov b1;buf cop b0 I3..14 J2..7;buf exp abr.txt;go")'转化字符为数据
Dim s As Stringi = 1
Open "C:\CVUSER\abr.txt" For Input As #1
Do While Not EOF(1)
Line Input #1, sstr_txt = Split(s, vbTab)For j = 0 To UBound(str_txt)str_num = CDbl(str_txt(j))Worksheets("Sheet1").Range(Cells(i, j + 1), Cells(i, j + 1)) = str_numNext j
i = i + 1
Loop
Close #1'关闭CODEV实例
Session.StopCodeV
Set Session = NothingEnd Sub2.CODE V macro
需要自定义一个宏
!**************************************************************************
! Sequence: THIRDABRGET
!
! Purpose: Prints a table for surface range i to j of third order
! aberrations
!
! This macro can be called stand-alone, but is intended to be
! called by FORDER.
!
! Syntax: in THIRDABRGET [start_surf [end_surf [zoom_pos]]]
!
! Inputs: start_surf - first surface of surface range (if -1,
! computes sums at image surface)
! end_surf - ending surface of surface range (if -1,
! only computes contributions for one surface;
! if 0, computes to image surface)
! zoom_pos - zoom position (default 1)
! list_opd - list OPDs (default: N)
! list_pup - list pupil aberrations (default: N)
!
! Notes: Assumes FIFTHDEF has been run first.
! Examples of usage are as follows:
! in THIRDABRGET 1 3 ! print contributions for surfaces 1 to 3
! in THIRDABRGET 5 0 ! print contributions from surface 5 to image
! in THIRDABRGET 5 ! print contributions for surface 5 only
! in THIRDABRGET ! print totals only
! in THIRDABRGET 0 0 0 Y! print wave aberrations
! in THIRDABRGET 0 0 0 N Y! print transverse pupil aberrations
!
! Author: YS Date: 2024-02-08
!!**************************************************************************
! ARG0 "Prints out a table of third and fifth order aberrations. Output"
! ARG0 "can be surface contributions for specified surfaces, or can "
! ARG0 "include aberration sums at the image (default)."
!
! ARG1 NAME "Start surface"
! ARG1 TYPE lit
! ARG1 DEFAULT -1
! ARG1 HELP "Starting surface (negative 1 gives sums at image only)"
!
! ARG2 NAME "Ending surface"
! ARG2 TYPE lit
! ARG2 DEFAULT -1
! ARG2 HELP "Ending surface (negative = starting surface only; 0 = image surface)"
!
! ARG3 NAME "Zoom position"
! ARG3 TYPE lit
! ARG3 DEFAULT 1
! ARG3 HELP "Zoom position (default 1)"! ARG4 NAME "List wave aberrations"
! ARG4 TYPE STR
! ARG4 DEFAULT "N"
! ARG4 HELP "List wave aberrations"! ARG5 NAME "List pupil aberrations"
! ARG5 TYPE STR
! ARG5 DEFAULT "N"
! ARG5 HELP "List pupil aberrations"
!**************************************************************************rfd -1 -1 1 N Nlcl num ^junk ^surface ^surface1 ^surface2 ^z
lcl num ^dum ^ref_wllcl str ^format2 ^format3 ^format4 ^format5 ^format6
lcl str ^format2w ^format3w ^format4w ^format5w ^format6w
lcl str ^errorMessage
lcl str ^tempString1
lcl str ^tempString2
lcl str ^tempString3
lcl str ^aber_type ^list_opd ^list_pupver n
^format2 == "STO'4d.6d''4d.6d''4d.6d''4d.6d''4d.6d''4d.6d'"
^format3 == "'3d''4d.6d''4d.6d''4d.6d''4d.6d''4d.6d''4d.6d'"
^format4 == "SUM'4d.6d''4d.6d''4d.6d''4d.6d''4d.6d''4d.6d'"
^format5 == " '4d.6d''4d.6d''4d.6d''4d.6d''4d.6d''4d.6d'"
^format6 == " '4d.6d''4d.6d''4d.6d''4d.6d'"
^format2w == "STO'6d.4d''6d.4d''6d.4d''6d.4d''6d.4d'"
^format3w == "'3d''6d.4d''6d.4d''6d.4d''6d.4d''6d.4d'"
^format4w == "SUM'6d.4d''6d.4d''6d.4d''6d.4d''6d.4d'"
^format5w == " '6d.4d''6d.4d''6d.4d''6d.4d''6d.4d'"
^format6w == " '6d.4d''6d.4d''6d.4d''6d.4d''6d.4d'"^surface1 == str_to_num(rfstr(1))^surface2 == str_to_num(rfstr(2))^z == str_to_num(rfstr(3))
^z == roundf(absf(^z))^list_opd == upcase(rfstr(4))
if ^list_opd <> "Y"^list_opd == "N"
end if
^aber_type == upcase(rfstr(5))
if ^aber_type <> "Y"^aber_type == "IMAGE"
else^aber_type == "PUPIL"
end ifif (num z) > 1if (^z < 1) or ^z > (num z)^tempString1 == 'The specified zoom position, '^tempString2 == NUM_TO_STR(^z)^tempString3 == CONCAT(^tempString1,^tempString2)^tempString1 == ', is not a valid input.'^errorMessage == CONCAT(^tempString3,^tempString1)^dum == cverror(^errorMessage,0)rtnend if
end ifif ^surface1 = -1^surface1 == (FOC z^z)
elseif ^surface1 = 0^surface1 == 1else if ^surface1 > (FOC z^z)^tempString1 == 'The specified starting surface, '^tempString2 == NUM_TO_STR(^surface1)^tempString3 == CONCAT(^tempString1,^tempString2)^tempString1 == ', is not a valid input.'^errorMessage == CONCAT(^tempString3,^tempString1)^dum == cverror(^errorMessage,0)rtnend if
end ifif ^surface2 = -1^surface2 == ^surface1
else if ^surface2 = 0^surface2 == (FOC z^z)
elseif ^surface2 < ^surface1^surface2 == ^surface1else if ^surface2 > (FOC z^z)^tempString1 == 'The specified ending surface, '^tempString2 == NUM_TO_STR(^surface2)^tempString3 == CONCAT(^tempString1,^tempString2)^tempString1 == ', is not a valid input.'^errorMessage == CONCAT(^tempString3,^tempString1)^dum == cverror(^errorMessage,0)rtnend if
end if! When ^surface is 0, the variables are initialized, when ^surface is (NUM S)
! the final sums are converted to aberrations, otherwise the surface
! contributions are calculated.! This is to force a warning message for a zero last field
^dummy == (hcy si)for ^surface 0 ^surface2out n^junk == @FORDER(^surface,^z, ^aber_type) ! Evaluate aberrations (dummy variable)out! Output aberrationsif ^list_opd = "N"if ^surface = 0! Get unitsif (dim) = "I"^tempString1 == "inches"else if (dim) = "C"^tempString1 == "centimeters"else if (dim) = "M"^tempString1 == "millimeters"end ifif (afc z^z) <> 0if (adm) = 'RAD'^tempString1 == "radians"else if (adm) = 'MR'^tempString1 == "milliradians"else if (adm) = 'UR'^tempString1 == "microradians"else if (adm) = 'NR'^tempString1 == "nanoradians"else if (adm) = 'DEG'^tempString1 == "degrees"else if (adm) = 'MIN'^tempString1 == "minutes"else if (adm) = 'SEC'^tempString1 == "seconds"end ifend ifif ^aber_type = "IMAGE"^tempString2 == CONCAT(CONCAT("THIRD AND FIFTH ORDER IMAGE ABERRATIONS (", ^tempString1), ")")!wri!wri " " ^tempString2!wri!wri " SA3 TCO TAS SAG PTZ DST"!wri " SA5 TCO5 TAS5 SAG5 PTZ5 DST5"!wri " SA7 ECOM OTSA OSSA"else^tempString2 == CONCAT(CONCAT("THIRD AND FIFTH ORDER PUPIL ABERRATIONS (", ^tempString1), ")")!wri!wri " " ^tempString2!wri!wri " SA3P TCOP TASP SAGP PTZP DSTP"!wri " SA5P TCO5P TAS5P SAG5P PTZ5P DST5P"!wri " SA7P ECOMP OTSAP OSSAP"end if!wriels if ^surface <= ^surface2 AND ^surface >= ^surface1if ^surface = (STO z^z)if ^aber_type = "IMAGE"write q^format2 ^SA3 ^TCO ^TAS ^SAG ^PTZ ^DST!write q^format2 ^SA3 ^TCO ^TAS ^SAG ^PTZ ^DSTelse!write q^format2 ^SA3P ^TCOP ^TASP ^SAGP ^PTZP ^DSTPend ifels if ^surface = (FOC z^z)!writeif ^aber_type = "IMAGE"write q^format4 ^SA3 ^TCO ^TAS ^SAG ^PTZ ^DST!write q^format4 ^SA3 ^TCO ^TAS ^SAG ^PTZ ^DSTelse!write q^format4 ^SA3P ^TCOP ^TASP ^SAGP ^PTZP ^DSTPend ifelsif ^aber_type = "IMAGE"write q^format3 ^surface ^SA3 ^TCO ^TAS ^SAG ^PTZ ^DST!write q^format3 ^SA3 ^TCO ^TAS ^SAG ^PTZ ^DSTelse!write q^format3 ^surface ^SA3P ^TCOP ^TASP ^SAGP ^PTZP ^DSTPend ifend ifif ^aber_type = "IMAGE"!write q^format5 ^SA5 ^TCO5 ^TAS5 ^SAG5 ^PTZ5 ^DST5!write q^format6 ^SA7 ^ECOM ^OTSA ^OSSAelse!write q^format5 ^SA5P ^TCO5P ^TAS5P ^SAG5P ^PTZ5P ^DST5P!write q^format6 ^SA7P ^ECOMP ^OTSAP ^OSSAPend if!writeend ifelseif ^surface = 0^ref_wl == (ref)! wriif ^aber_type = "IMAGE"!wri q" '60c' '6d.4d' '4c'" " THIRD AND FIFTH ORDER IMAGE ABERRATIONS (Waves at " (wl w^ref_wl) "nm.)"!wri!wri " W040 W131 W220 W222 W311"!wri " W060 W151 W420 W422 W511"!wri " W080 W331 W333 W240 W242"else!wri q" '60c' '6d.4d' '4c'" " THIRD AND FIFTH ORDER PUPIL ABERRATIONS (Waves at " (wl w^ref_wl) "nm.)"! wri!wri " W040P W131P W220P W222P W311P"!wri " W060P W151P W420P W422P W511P"!wri " W080P W331P W333P W240P W242P"end if!wriels if ^surface <= ^surface2 AND ^surface >= ^surface1if ^surface = (STO z^z)if ^aber_type = "IMAGE"!write q^format2w ^W040 ^W131 ^W220 ^W222 ^W311else!write q^format2w ^W040P ^W131P ^W220P ^W222P ^W311Pend ifels if ^surface = (FOC z^z)writeif ^aber_type = "IMAGE"!write q^format4w ^W040 ^W131 ^W220 ^W222 ^W311else!write q^format4w ^W040P ^W131P ^W220P ^W222P ^W311Pend ifelsif ^aber_type = "IMAGE"!write q^format3w ^surface ^W040 ^W131 ^W220 ^W222 ^W311else!write q^format3w ^surface ^W040P ^W131P ^W220P ^W222P ^W311Pend ifend ifif ^aber_type = "IMAGE"!write q^format5w ^W060 ^W151 ^W420 ^W422 ^W511!write q^format6w ^W080 ^W331 ^W333 ^W240 ^W242else!write q^format5w ^W060P ^W151P ^W420P ^W422P ^W511P!write q^format6w ^W080P ^W331P ^W333P ^W240P ^W242Pend if!writeend ifend ifend for
ver y输出结果
code的输出
VBA输出
![BUUCTF-Real-[Tomcat]CVE-2017-12615](https://img-blog.csdnimg.cn/direct/5ae00e6044ea4453947915b6b513c052.png)
