accosin.m

来自「zemax与matlab调用程序」· M 代码 · 共 869 行 · 第 1/4 页

function Fixups = ReadNewConfigData(fID, Config, ACCOSSurface, InputFixups)
global ACCOSMaterials ZEMAXMaterials
ZEMAXSurface = zFindLabel(ACCOSSurface);
Fixups = InputFixups;
disp('Processing New Config Data.');
while 1
  CfgLine = fgetl(fID);
  if ~ischar(CfgLine), return; end; % (Unexpected) End of file
  CfgLine = upper(CfgLine); % Convert to upper case
  [Directive, Rest] = strtok(CfgLine, ' ,'); % Take the directive from the start of the line
  switch Directive
      case 'CFG', [Token, Rest] = strtok(Rest, ' ,'); Config = str2num(Token);
                  Fixups.NumberConfigs = max([Fixups.NumberConfigs, Config]);
                  [Token, Rest] = strtok(Rest, ' ,'); if (Token), ACCOSSurface = str2num(Token); ZEMAXSurface = zFindLabel(ACCOSSurface); end;
      case 'SAY', ApertureData = zGetSystemAper; EPD = ApertureData(3);
                  Fixups.ConfigOperands = Fixups.ConfigOperands + 1;
                  Fixups.ConfigData{Fixups.ConfigOperands} = sprintf('APER 0 1 %11.9f 0', EPD); 
                  [Token, Rest] = strtok(Rest, ' ,'); EPD = 2 * str2double(Token);
                  Fixups.ConfigOperands = Fixups.ConfigOperands + 1;
                  Fixups.ConfigData{Fixups.ConfigOperands} = sprintf('APER 0 %i %11.9f 0',Config, EPD);                   
      case 'SURF',[Token, Rest] = strtok(Rest, ' ,'); ACCOSSurface = str2num(Token); ZEMAXSurface = zFindLabel(ACCOSSurface);
      case 'SCY', [Token, Rest] = strtok(Rest, ' ,');
                  if (strcmp(Token, 'FANG'))
                      Angle = zGetField(2); Angle = Angle(2);
                      Fixups.ConfigOperands = Fixups.ConfigOperands + 1;
                      Fixups.ConfigData{Fixups.ConfigOperands} = sprintf('YFIE 2 1 %11.9f 0', Angle); 
                      [Token, Rest] = strtok(Rest, ' ,');
                      Angle = str2double(Token);
                      Fixups.ConfigOperands = Fixups.ConfigOperands + 1;
                      Fixups.ConfigData{Fixups.ConfigOperands} = sprintf('YFIE 2 %i %11.9f 0', Config, Angle); 
                  end;  
      case 'WV', [Token, Rest] = strtok(Rest, ' ,'); WaveNumber = 0;
                 while (Token)
                     WaveNumber = WaveNumber + 1;
                     WV = zGetWave(WaveNumber); WV = WV(1);
                     Fixups.ConfigOperands = Fixups.ConfigOperands + 1;
                     Fixups.ConfigData{Fixups.ConfigOperands} = sprintf('WAVE %i 1 %11.9f 0', WaveNumber, WV);
                     WV = str2double(Token);
                     Fixups.ConfigOperands = Fixups.ConfigOperands + 1;
                     Fixups.ConfigData{Fixups.ConfigOperands} = sprintf('WAVE %i %i %11.9f 0', WaveNumber, Config, WV);
                     [Token, Rest] = strtok(Rest, ' ,');
                 end;
     case 'SPTWT',[Token, Rest] = strtok(Rest, ' ,'); WaveNumber = 0;
                 while (Token)
                     WaveNumber = WaveNumber + 1;
                     Weight = zGetWave(WaveNumber); Weight = Weight(2);
                     Fixups.ConfigOperands = Fixups.ConfigOperands + 1;
                     Fixups.ConfigData{Fixups.ConfigOperands} = sprintf('WLWT %i 1 %11.9f 0', WaveNumber, Weight);
                     Weight = str2double(Token);
                     Fixups.ConfigOperands = Fixups.ConfigOperands + 1;
                     Fixups.ConfigData{Fixups.ConfigOperands} = sprintf('WLWT %i %i %11.9f 0', WaveNumber, Config, Weight);
                     [Token, Rest] = strtok(Rest, ' ,');
                 end;
      case 'ASTOP',[Token, Rest] = strtok(Rest, ' ,');
                 ApertureData = zGetSystemAper; StopSurf = ApertureData(2);
                 Fixups.ConfigOperands = Fixups.ConfigOperands + 1;
                 Fixups.ConfigData{Fixups.ConfigOperands} = sprintf('STPS 0 1 %i 0', StopSurf);
                 Fixups.ConfigOperands = Fixups.ConfigOperands + 1;
                 Fixups.ConfigData{Fixups.ConfigOperands} = sprintf('STPS 0 %i %i 0', Config, ZEMAXSurface); 
      case 'CLAP',[Token, Rest] = strtok(Rest, ' ,');
                 Apmx = zGetAperture(ZEMAXSurface); Apmx = Apmx(3);
                 Fixups.ConfigOperands = Fixups.ConfigOperands + 1;
                 Fixups.ConfigData{Fixups.ConfigOperands} = sprintf('APMX %i 1 %11.9f 0', ZEMAXSurface, Apmx);
                 Apmx = str2double(Token);
                 Fixups.ConfigOperands = Fixups.ConfigOperands + 1;
                 Fixups.ConfigData{Fixups.ConfigOperands} = sprintf('APMX %i %i %11.9f 0', ZEMAXSurface, Config, Apmx);
      case 'CLAPD', % Handle CLAPD by setting aperture large 
                 Apmx = zGetAperture(ZEMAXSurface); Apmx = Apmx(3);
                 Fixups.ConfigOperands = Fixups.ConfigOperands + 1;
                 Fixups.ConfigData{Fixups.ConfigOperands} = sprintf('APMX %i 1 %11.9f 0', ZEMAXSurface, Apmx);
                 Apmx = 10000.0;
                 Fixups.ConfigOperands = Fixups.ConfigOperands + 1;
                 Fixups.ConfigData{Fixups.ConfigOperands} = sprintf('APMX %i %i %11.9f 0', ZEMAXSurface, Config, Apmx);
      case 'CV', Curvature = zGetSurfaceData(ZEMAXSurface, 2);
                 Fixups.ConfigOperands = Fixups.ConfigOperands + 1;
                 Fixups.ConfigData{Fixups.ConfigOperands} = sprintf('CRVT %i 1 %11.9f 0', ZEMAXSurface, Curvature);
                 [Token, Rest] = strtok(Rest, ' ,'); Curvature = str2double(Token);
                 Fixups.ConfigOperands = Fixups.ConfigOperands + 1;
                 Fixups.ConfigData{Fixups.ConfigOperands} = sprintf('CRVT %i %i %11.9f 0', ZEMAXSurface, Config, Curvature);
      case 'CC', Conic = zGetSurfaceData(ZEMAXSurface, 6);
                 Fixups.ConfigOperands = Fixups.ConfigOperands + 1;
                 Fixups.ConfigData{Fixups.ConfigOperands} = sprintf('CONN %i 1 %11.9f 0', ZEMAXSurface, Conic);          
                 [Token, Rest] = strtok(Rest, ' ,'); Conic = str2double(Token);
                 Fixups.ConfigOperands = Fixups.ConfigOperands + 1;
                 Fixups.ConfigData{Fixups.ConfigOperands} = sprintf('CONN %i %i %11.9f 0', ZEMAXSurface, Config, Conic);
      case 'AD', AD = zGetSurfaceParameter(ZEMAXSurface, 2);
                 Fixups.ConfigOperands = Fixups.ConfigOperands + 1;
                 Fixups.ConfigData{Fixups.ConfigOperands} = sprintf('PAR2 %i 1 %11.9f 0', ZEMAXSurface, AD);                 
                 [Token, Rest] = strtok(Rest, ' ,'); AD = str2double(Token);
                 Fixups.ConfigOperands = Fixups.ConfigOperands + 1;
                 Fixups.ConfigData{Fixups.ConfigOperands} = sprintf('PAR2 %i %i %11.9f 0', ZEMAXSurface, Config, AD);
      case 'AE', AE = zGetSurfaceParameter(ZEMAXSurface, 3);
                 Fixups.ConfigOperands = Fixups.ConfigOperands + 1;
                 Fixups.ConfigData{Fixups.ConfigOperands} = sprintf('PAR3 %i 1 %11.9f 0', ZEMAXSurface, AE);          
                 [Token, Rest] = strtok(Rest, ' ,'); AE = str2double(Token);
                 Fixups.ConfigOperands = Fixups.ConfigOperands + 1;
                 Fixups.ConfigData{Fixups.ConfigOperands} = sprintf('PAR3 %i %i %11.9f 0', ZEMAXSurface, Config, AE);
      case 'AF', AF = zGetSurfaceParameter(ZEMAXSurface, 4);
                 Fixups.ConfigOperands = Fixups.ConfigOperands + 1;
                 Fixups.ConfigData{Fixups.ConfigOperands} = sprintf('PAR4 %i 1 %11.9f 0', ZEMAXSurface, AF);          
                 [Token, Rest] = strtok(Rest, ' ,'); AF = str2double(Token);
                 Fixups.ConfigOperands = Fixups.ConfigOperands + 1;
                 Fixups.ConfigData{Fixups.ConfigOperands} = sprintf('PAR4 %i %i %11.9f 0', ZEMAXSurface, Config, AF);
      case 'AG', AG = zGetSurfaceParameter(ZEMAXSurface, 2);
                 Fixups.ConfigOperands = Fixups.ConfigOperands + 1;
                 Fixups.ConfigData{Fixups.ConfigOperands} = sprintf('PAR5 %i 1 %11.9f 0', ZEMAXSurface, AG);          
                 [Token, Rest] = strtok(Rest, ' ,'); AG = str2double(Token);
                 Fixups.ConfigOperands = Fixups.ConfigOperands + 1;
                 Fixups.ConfigData{Fixups.ConfigOperands} = sprintf('PAR5 %i %i %11.9f 0', ZEMAXSurface, Config, AG);
      case 'TH', Thickness = zGetSurfaceData(ZEMAXSurface, 3);
                 Fixups.ConfigOperands = Fixups.ConfigOperands + 1;
                 Fixups.ConfigData{Fixups.ConfigOperands} = sprintf('THIC %i 1 %11.9f 0', ZEMAXSurface, Thickness);          
                 [Token, Rest] = strtok(Rest, ' ,'); Thickness = str2double(Token);
                 Fixups.ConfigOperands = Fixups.ConfigOperands + 1;
                 Fixups.ConfigData{Fixups.ConfigOperands} = sprintf('THIC %i %i %11.9f 0', ZEMAXSurface, Config, Thickness); 
      case 'RD', Curvature = zGetSurfaceData(ZEMAXSurface, 2);
                 Fixups.ConfigOperands = Fixups.ConfigOperands + 1;
                 Fixups.ConfigData{Fixups.ConfigOperands} = sprintf('CRVT %i 1 %11.9f 0', ZEMAXSurface, Curvature);          
                 [Token, Rest] = strtok(Rest, ' ,'); Radius = str2double(Token);
                 Fixups.ConfigOperands = Fixups.ConfigOperands + 1;
                 if (Radius ~= 0)
                   Fixups.ConfigData{Fixups.ConfigOperands} = sprintf('CRVT %i %i %11.9f 0', ZEMAXSurface, Config, 1/Radius); 
                 else
                   Fixups.ConfigData{Fixups.ConfigOperands} = sprintf('CRVT %i %i 0.0 0', ZEMAXSurface, Config); 
                 end;
      case {'SCHOTT', 'HOYA', 'OHARA', 'MATL', 'GLASS'},
                 Glass = zGetSurfaceData(ZEMAXSurface, 4);
                 Fixups.ConfigOperands = Fixups.ConfigOperands + 1;
                 Fixups.ConfigData{Fixups.ConfigOperands} = sprintf('GLSS %i 1 %s', ZEMAXSurface, Glass);                
                 Token = strtok(Rest, ' ,'); 
                 Found = find(strcmp(ACCOSMaterials, Token)); % Look for ZEMAX equivalent
                 if (Found) Glass = ZEMAXMaterials{Found(1)}; else Glass = Token; end; % Take a flyer if not found
                 Fixups.ConfigOperands = Fixups.ConfigOperands + 1;
                 Fixups.ConfigData{Fixups.ConfigOperands} = sprintf('GLSS %i %i %s', ZEMAXSurface, Config, Glass); 
      case 'REFL',
                 Glass = zGetSurfaceData(ZEMAXSurface, 4);
                 Fixups.ConfigOperands = Fixups.ConfigOperands + 1;
                 Fixups.ConfigData{Fixups.ConfigOperands} = sprintf('GLSS %i 1 %s', ZEMAXSurface, Glass); 
                 Fixups.ConfigOperands = Fixups.ConfigOperands + 1;
                 Fixups.ConfigData{Fixups.ConfigOperands} = sprintf('GLSS %i %i MIRROR', ZEMAXSurface, Config); 
      case 'AIR', 
                 Glass = zGetSurfaceData(ZEMAXSurface, 4);
                 Fixups.ConfigOperands = Fixups.ConfigOperands + 1;
                 Fixups.ConfigData{Fixups.ConfigOperands} = sprintf('GLSS %i 1 %s', ZEMAXSurface, Glass); 
                 Fixups.ConfigOperands = Fixups.ConfigOperands + 1;
                 Fixups.ConfigData{Fixups.ConfigOperands} = sprintf('GLSS %i %i', ZEMAXSurface, Config); 
      case 'EOS', return;
      otherwise, disp(['Config Directive "' Directive '" Ignored.']);
  end;
end;


function DoFixups(LensFileName, Fixups)
% In this function a number of touchups are performed directly on the .zmx file created with
% zSaveFile above. This is mandatory because certain things cannot currently be done through the
% ZEMAX DDE server. There are also some bugs requiring workarounds eg. model glasses.
% disp(['Doing fixups on ' LensFileName]);
fin = fopen([LensFileName '.tmp'], 'rt'); % This is the file needing fixups
fout = fopen(LensFileName, 'wt'); % This is name of the touched up file.
while 1
    LenLine = fgetl(fin);
    if ~ischar(LenLine), break, end; % End of file
    LenLine = upper(LenLine); % Convert to upper case
    [Directive, Rest] = strtok(LenLine, ' ,'); % Take the directive from the start of the line
    if (Directive)
      switch Directive
        case 'NAME', fprintf(fout, 'NAME %s\n', Fixups.Title);    % Put in the real title
        case 'NOTE', NoteNumber = str2num(strtok(Rest, ' '));
                     if (NoteNumber <= size(Fixups.Comments,2))
                         fprintf(fout, 'NOTE %1i %s\n', NoteNumber, Fixups.Comments{NoteNumber});
                     else
                         fprintf(fout, '%s\n', LenLine);
                     end;
        case 'SURF', SurfNumber = str2num(strtok(Rest, ' ')); fprintf(fout, '%s\n', LenLine);
                     if (SurfNumber > 0)&(SurfNumber <= size(Fixups.TiltDecBefore,2)),
                       if (Fixups.TiltDecBefore{SurfNumber}),
                         fprintf(fout, '  %s\n', Fixups.TiltDecBefore{SurfNumber});
                         %fprintf(fout, '  %s\n', Fixups.TiltDecAfter{SurfNumber});                         
                       end;
                     end;
        case 'GLAS', if (SurfNumber > 0)&(SurfNumber <= size(Fixups.ModelGlasses,2)),  % Here we patch the model glasses just a little
                       if (Fixups.ModelGlasses{SurfNumber}),
                           fprintf(fout, '  %s\n', Fixups.ModelGlasses{SurfNumber});
                       else
                           fprintf(fout, '%s\n', LenLine);
                       end;
                     else
                       fprintf(fout, '%s\n', LenLine);
                     end;
        case 'MNUM', if (Fixups.NumberConfigs > 1),
                       fprintf(fout, 'MNUM %i\n', Fixups.NumberConfigs);
                       if (Fixups.ConfigOperands > 0)
                         Fixups.ConfigData = sort(Fixups.ConfigData); % Sort the data before writing into the ZEMAX file
                         fprintf(fout, '%s\n', Fixups.ConfigData{1});
                         if (Fixups.ConfigOperands > 1),
                           for i = 2:Fixups.ConfigOperands,
                             if (~strcmp(Fixups.ConfigData{i}, Fixups.ConfigData{i-1})), % Avoid writing duplicate lines
                               fprintf(fout, '%s\n', Fixups.ConfigData{i});
                             end;
                           end;
                         end;
                       end;
                     end;
        case 'MOFF', % Write nothing
        otherwise,   fprintf(fout, '%s\n', LenLine);
      end;
    end;
end;
fclose(fin);
fclose(fout);
delete([LensFileName '.tmp']);