Login| Sign Up| Help| Contact|

Patent Searching and Data


Title:
METHOD FOR CREATING A COLORED LASER MARKING
Document Type and Number:
WIPO Patent Application WO/2021/259826
Kind Code:
A1
Abstract:
A method (1) for preparing a laser marking system (100) to create a colored laser mark on a specimen comprising the following steps: a) Providing a laser marking system (100) and a specimen (105) comprising a surface layer (105a), wherein the laser marking system comprises a preset number of laser parameters (12); b) Performing an exploration of a first gamut (2) specified by the laser marking system (100) and the specimen (105) comprising a surface layer (105a) including the following steps: aa) Creating (3) a design space (10) with a preset number of design points (11), wherein each design point (11) represents a combination of the preset number of laser parameters (12); bb) Performing (4) a marking of a sample on the specimen (105) for each design point (11); cc) Measuring (5) the sample using at least one detection device (106) and deter- mine for each design point a performance point (14), wherein the measured performance points (14) define a performance space (13); dd) Evaluating (6) the performance space (13) with regard to preset performance criteria using an evaluation device (107), wherein a Pareto front is determined comprising a subset of performance points; ee) Generating (7) an offspring design space (10a) with offspring design points (11a); ff) Creating (8) a first gamut (2) using the subset of performance points forming the Pareto front; wherein the steps bb) to dd) are iterated (9) for a preset iteration number, wherein in each iteration (9) the offspring design space (10a) of the previous iteration is used in step bb), wherein in each iteration the measured performance space is combined (15) with the performance space of the previous iteration (9) such that in step dd) the combined performance space (13a) is used.

Inventors:
SEIDEL HANS-PETER (DE)
BABAEI VAHID (DE)
CUCERCA SEBASTIAN (DE)
DIDYK PIOTR (CH)
Application Number:
PCT/EP2021/066753
Publication Date:
December 30, 2021
Filing Date:
June 21, 2021
Export Citation:
Click for automatic bibliography generation   Help
Assignee:
MAX PLANCK GESELLSCHAFT (DE)
International Classes:
B41M5/26; B41M5/34; G11B7/00; G11B7/005
Foreign References:
US5982924A1999-11-09
US20100149203A12010-06-17
US20090141295A12009-06-04
Other References:
DAVID PRICE ADAMSRYAN D MURPHYDAVID J SAIZDA HIRSCHFELDMA RODRIGUEZPAUL GA-BRIEL KOTULABH JARED: "Nanosecond pulsed laser irradiation of titanium: Oxide growth and effects on underlying metal", SURFACE AND COATINGS TECHNOLOGY, vol. 248, 2014, pages 38 - 45, XP028649424, DOI: 10.1016/j.surfcoat.2013.12.052
ARKADIUSZ J ANTONCZAKDARIUSZ KOCORIMACIEJ NOWAKPAWET KOZIOTKRZYSZTOF M ABRAMSKI: "Laser-induced colour marking - Sensitivity scaling for a stainless stee", APPLIED SURFACE SCIENCE, 2013
ARKADIUSZ J ANTONCZAKBOGUSZ STEPAKPAWET E KOZIOTKRZYSZTOF M ABRAMSKI: "The influence of process parameters on the laser-induced coloring of titanium", APPLIED PHYSICS A, 2014
THOMAS AUZINGERWOLFGANG HEIDRICHBERND BICKEL: "Computational design of nanostructural color for additive manufacturing", ACM TRANSACTIONS ON GRAPHICS (TOG, vol. 37, no. 4, 2018, pages 159
VAHID BABAEIROGER D HERSCH: "Juxtaposed color halftoning relying on discrete lines", IEEE TRANSACTIONS ON IMAGE PROCESSING, vol. 22, no. 2, 2012, pages 679 - 686, XP011492276, DOI: 10.1109/TIP.2012.2221727
VAHID BABAEIROGER D HERSCH: "Color Imaging XX: Displaying, Processing, Hardcopy, and Applications", vol. 9395, 2015, INTERNATIONAL SOCIETY FOR OPTICS AND PHOTONICS, article "Yule-Nielsen based multi-angle reflectance prediction of metallic halftones"
VAHID BABAEIROGER D HERSCH: "N-Ink Printer Characterization With Barycentric Subdivision", IEEE TRANSACTIONS ON IMAGE PROCESSING, vol. 25, no. 7, 2016, pages 3023 - 3031, XP011609870, DOI: 10.1109/TIP.2016.2560526
VAHID BABAEIROGER D HERSCH: "Color reproduction of metallic-ink images", JOURNAL OF IMAGING SCIENCE AND TECHNOLOGY, vol. 60, no. 3, 2016, pages 30503 - 1
MILTON BIRNBAUM: "Modulation of the reflectivity of semiconductors", JOURNAL OF APPLIED PHYSICS, vol. 36, no. 2, 1965, pages 657 - 658
KALYANMOY DEBAMRIT PRATAPSAMEER AGARWALTAMT MEYARIVAN: "A fast and elitist multiobjective genetic algorithm: NSGA-II", IEEE TRANSACTIONS ON EVOLUTIONARY COMPUTATION, vol. 6, no. 2, 2002, pages 182 - 197, XP011072884, DOI: 10.1109/4235.996017
A PEREZ DEL PINOJM FERNANDEZ-PRADASP SERRAJL MORENZA: "Coloring of titanium through laser oxidation: comparative study with anodizing", SURFACE AND COATINGS TECHNOLOGY, vol. 187, no. 1, 2004, pages 106 - 112, XP002518792
CARLOS M FONSECAPETER J FLEMING ET AL.: "Genetic Algorithms for Multiobjective Optimization: Formulation, Discussion and Generalization", ICGA, vol. 93, 1993, pages 416 - 423, XP055165790
MATHIEU HEBERT: "Color Imaging XIX: Displaying, Processing, Hardcopy, and Applications", vol. 9015, 2014, INTERNATIONAL SOCIETY FOR OPTICS AND PHOTONICS, article "Yule-Nielsen effect in halftone prints: graphical analysis method and improvement of the Yule-Nielsen transform"
ROGER D HERSCHPHILIPP DONZESYLVAIN CHOSSON: "ACM Transactions on Graphics (TOG", vol. 26, 2007, ACM, article "Color images visible under UV light", pages: 75
GUOWEI HONGM RONNIER LUOPETER A RHODES: "A study of digital camera colorimetric characterization based on polynomial modeling", COLOR RESEARCH & APPLICATION: ENDORSED BY INTER-SOCIETY COLOR COUNCIL, THE COLOUR GROUP (GREAT BRITAIN), CANADIAN SOCIETY FOR COLOR, COLOR SCIENCE ASSOCIATION OF JAPAN, DUTCH SOCIETY FOR THE STUDY OF COLOR, THE SWEDISH COLOUR CENTRE FOUNDATION, COLOU, vol. 26, no. 1, 2001, pages 76 - 84, XP009006927, DOI: 10.1002/1520-6378(200102)26:1<76::AID-COL8>3.0.CO;2-3
P LAAKSOS RUOTSALAINENH PANTSARR PENTTILA: "Relation of laser parameters in color marking of stainless steel", 12TH CONFERENCE ON LASER PROCESSING OF MATERIALS IN THE NORDIC COUNTRIES, 2009
C LANGLADEAB VANNESJM KRAFFTJR MARTIN: "Surface modification and tribologi-cal behaviour of titanium and titanium alloys after YAG-laser treatments", SURFACE AND COATINGS TECHNOLOGY, vol. 100, 1998, pages 383 - 387
SAMANTHA K LAWRENCEDAVID P ADAMSDAVID F BAHRNEVILLE R MOODY: "Mechanical and electromechanical behavior of oxide coatings grown on stainless steel 304L by nanosecond pulsed laser irradiation", SURFACE AND COATINGS TECHNOLOGY, vol. 235, 2013, pages 860 - 866, XP028749613, DOI: 10.1016/j.surfcoat.2013.09.013
KM TERCKAAJ ANTONCZAKB SZUBZDAMR WOJCIKBD STERPAKP SZYMCZYKM TRZCINSKIM OZIMEKKM ABRAMSKI: "Effects of laser-induced oxidation on the corrosion resistance of AISI 304 stainless steel. Journal of Laser Applications 28, 3 (2016), 032009. Achim Lewandowski, Marcus Ludl, Gerald Byrne, and Georg Dorffner. 2006", APPLYING THE YULE-NIELSEN EQUATION WITH NEGATIVE N. JOSA A, vol. 23, no. 8, 2006, pages 1827 - 1834
HUAGANG LIUWENXIONG LINMINGHUI HONG: "Surface coloring by laser irradiation of solid substrates", APL PHOTONICS, vol. 4, no. 5, 2019, pages 051101, XP012237815, DOI: 10.1063/1.5089778
LASZLO NANAIROBERT VAJTAITHOMAS F GEORGE.: "Laser-induced oxidation of metals: state of the art.", THIN SOLID FILMS, vol. 298, no. 1-2, 1997, pages 160 - 164
PETAR PJANICROGER D HERSCH: "Color and Imaging Conference", vol. 2013, 2013, SOCIETY FOR IMAGING SCIENCE AND TECHNOLOGY, article "Specular color imaging on a metallic substrate", pages: 61 - 68
JEAN-PIERRE REVEILLES: "Vision geometry", vol. 2573, 1995, INTERNATIONAL SOCIETY FOR OPTICS AND PHOTONICS, article "Combinatorial pieces in digital lines and planes", pages: 23 - 34
ROBERT ROLLESTONRAJA BALASUBRAMANIAN: "Color and Imaging Conferenc", vol. 1993, 1993, SOCIETY FOR IMAGING SCIENCE AND TECHNOLOGY, article "Accuracy of various types of Neugebauer model", pages: 32 - 37
JANOS SCHANDA: "Colorimetry: Understanding the CIE System", 2007, JOHN WILEY & SONS
ADRIANA SCHULZHARRISONWANGEITAN CRINSPUNJUSTIN SOLOMONWOJCIECH MATUSIK: "Interactive exploration of design trade-offs", ACM TRANSACTIONS ON GRAPHICS (TOG, vol. 37, no. 4, 2018, pages 131
CHRISTIAN SCHUMACHERBERND BICKELJAN RYSSTEVE MARSCHNERCHIARA DARAIOMARKUS GROSS: "Microstructures to control elasticity in 3D printing", ACM TRANSACTIONS ON GRAPHICS (TOG, vol. 34, no. 4, 2015, pages 136, XP055292197, DOI: 10.1145/2766926
GAURAV SHARMAWENCHENG WUEDUL N DALAL: "The CIEDE2000 color-difference formula: Implementation notes, supplementary test data, and mathematical observations", COLOR RESEARCH AND APPLICATION, vol. 30, no. 1, 2005, pages 21 - 30
PITCHAYA SITTHI-AMORNNICHOLAS MODLYWESTLEY WEIMERJASON LAWRENCE: "Genetic programming for shader simplification", ACM TRANSACTIONS ON GRAPHICS (TOG, vol. 30, no. 6, 2011, pages 1 - 12
ERIC J STOLLNITZVICTOR OSTROMOUKHOVDAVID H SALESIN: "Proceedings of the 25th annual conference on Computer graphics and interactive techniques", 1998, ACM, article "Reproducing color images using custom inks", pages: 267 - 274
VADIM VEIKOGALINA ODINTSOVAELENA GORBUNOVAEDUARD AGEEVALEXANDR SHIMKOYULIA KARLAGINAYAROSLAVA ANDREEVA: "Development of complete color palette based on spectrophotometric measurements of steel oxidation results for enhancement of color laser marking technology", MATERIALS & DESIGN, vol. 89, 2016, pages 684 - 688, XP029303792, DOI: 10.1016/j.matdes.2015.10.030
VP VEIKOAA SLOBODOVGV ODINTSOVA: "Availability of methods of chemical thermodynamics and kinetics for the analysis of chemical transformations on metal surfaces under pulsed laser action", LASER PHYSICS, vol. 23, no. 6, 2013, pages 066001
RUI WANGBOWEN YUJULIO MARCOTIANLEI HUDIEGO GUTIERREZHUJUN BAO: "Real-time rendering on a power budget", ACM TRANSACTIONS ON GRAPHICS (TOG, vol. 35, no. 4, 2016, pages 1 - 11, XP058275773, DOI: 10.1145/2897824.2925889
GUNTER WYSZECKIWALTER STANLEY STILES: "Color Science", vol. 8, 1982, WILEY NEW YORK
JAC YULEWJ NIELSEN: "The penetration of light into paper and its effect on halftone reproduction", PROC. TAGA, vol. 3, 1951, pages 65 - 76, XP000981221
BO ZHUMELINA SKOURASDESAI CHENWOJCIECH MATUSIK: "Two-scale topology optimization with microstructures", ACM TRANSACTIONS ON GRAPHICS (TOG, vol. 36, no. 5, 2017, pages 164
Attorney, Agent or Firm:
HANNKE BITTNER & PARTNER PATENT- UND RECHTSANWÄLTE MBB (DE)
Download PDF:
Claims:
Claims

1. A method (1) for preparing a laser marking system (100) to reproduce a laser-marked color image on a specimen comprising the following steps: a) Providing a laser marking system (100) and a specimen (105) comprising a sur face layer (105a), wherein the laser marking system comprises a preset number of laser parameters (12); b) Performing an exploration of a first gamut (2) specified by the laser marking sys tem (100) and the specimen (105) comprising a surface layer (105a) including the following steps: aa) Creating (3) a design space (10) with a preset number of design points (11), wherein each design point (11) comprises a combination of the preset num ber of laser parameters (12); bb) Performing (4) a marking of a sample on the specimen (105) for each design point (11); cc) Measuring (5) the sample using at least one detection device (106) and de termine for each design point a performance point (14), wherein the meas ured performance points (14) define a performance space (13); dd) Evaluating (6) the performance space (13) with regard to preset performance criteria using an evaluation device (107), wherein a Pareto front is deter mined comprising a subset of performance points; ee) Generating (7) an offspring design space (10a) with offspring design points (11a); ff) Creating (8) a first gamut (2) using the subset of performance points forming the Pareto front; wherein the steps bb) to dd) are iterated (9) for a preset iteration number, wherein in each iteration (9) the offspring design space (10a) of the previous iteration is used in step bb), wherein in each iteration the measured performance space is combined (15) with the performance space of the previous iteration (9) such that in step dd) the com bined performance space (13a) is used.

2. The method (1) according to claim 1 , wherein the laser system (100) comprises at least one pulsed laser (101) and at least one scanning device (103, 104), wherein by the scanning device (103, 104) a laser spot is movable relative to the specimen (105) or wherein by the scanning device (103, 104) the specimen (105) is movable relative to the laser spot.

3. The method (1) according to claim 2, wherein a design point (11, 11a) comprises at least one laser parameter (12) selected form: the frequency of the laser pulses, the power of a laser pulse, the width of a laser pulse, the speed of the laser beam relative to the specimen along a vector while marking, the line count, which defines the numbers of lines in a cluster representing the marked sample, the distance between the lines within a cluster representing the marked sample, the number of times a vector is marked, wherein a design point (11 , 11a) further comprises the parameter focal distance of the laser beam, type of me dium gas, ambient temperature.

4. The method (1) according to one of the claims 1 to 3, wherein the performance criteria in step dd) comprise at least one of: chromaticity, hue spread, resolution, performance space diversity, design space diversity, color repeat ability, color uniformity.

5. The method (1) according to one of the claims 1 to 3, wherein performance criteria in step dd) comprise at least one of: chromaticity, resolution, per formance space diversity, design space diversity, wherein the performance points (14) are projected in to a CIECH space , wherein a the CIECH space is split into a first number of circular sectors (15) forming a hue wheel (16), wherein the per formance points (14) within each sector (15) of the hue wheel (16) are evaluated re garding said performance criteria, wherein said evaluation is iterated for a preset iter ation number, wherein in each iteration the number of sectors (15) forming a hue wheel (16) is altered, wherein each performance point (14) is characterized by a fre quency vector, which represents the presence in a certain Pareto front.

6. The method (1) according to one of the previous claims, wherein an additional evaluation regarding the achromatic properties of the performance points (14) is performed by performing step b) using the performance criteria in step dd) lightness, resolution, performance space diversity, design space diversity.

7. The method (1) according to one of the previous claims, further comprising the step selecting a set of primary colors from the first gamut (2), wherein the selected primary colors form a second gamut.

8. The method (1) according to one of the previous claims, wherein the data relating to the design space (10, 10a, 10b) and the performance space (13, 13a) of the first gamut (2) and/or data related to the second gamut are stored in a da tabase (109).

9. A method (20) for reproducing a laser-marked color image on a specimen (105) com prising a surface layer (105a) comprising the following steps: a) Verifying (21) the database (109) regarding data related to the first gamut (2) and/or second gamut with regard to the type of the specimen (2) and the laser marking system (100), wherein said data is obtained by the method (1) for pre paring a laser marking system (100) according to one of the claims 1 to 8; b) Retrieving (22) data related to the first gamut (2) and/or second gamut from the database (109) or perform (23) the method (1) for preparing a laser marking sys tem (100) according to one of the claims 1 to 8; c) Providing (24) an input image (27) to be reproduced as laser marking on the specimen (105); d) Performing (25) a color management workflow (28, 28a) by which creates control data for the laser marking system derived from the input image (27); e) Perform (26) the marking according to the control data.

10. The method (20) according to claim 9, wherein the color management workflow (28) is a juxtaposed halftoning workflow (28a).

11. The method (20) according to one of claims 9 or 10, wherein the color management workflow (28) comprises the steps: aa) Applying (29) a forward color prediction model to construct a third gamut with regard to the second gamut and the use of juxtaposed halftoning; bb) Mapping (30) the input image (27) into the third gamut; cc) Perform (31) a color separation such that for each mapped color a corre sponding area coverage of each primary color is determined; dd) Binarize (32) the area-coverages using the juxtaposed halftoning method and create (33) raster halftone images; HA- 'WO 2021/259826 PARTNER PCT/EP2021/066753 )PT

- 39 - ee) Convert (34) the raster halftone images into vector data, wherein the control data comprise the said vector data.

12. The method (20) according to one of the previous claims, wherein the specimen (105) has a metallic surface layer (105a), wherein the laser marking is based on laser induced oxidation of the surface layer (105a) of the specimen (105) or laser induced structuring of the surface layer (105a) of the specimen (105) or the la ser induced generation of micro/nanoparticles on the surface layer (105a) of the spec imen (105).

Description:
Method for creating a colored laser marking

Description

The invention relates to a method for preparing a laser marking system to create a colored laser marking on a specimen and to a method for creating a colored laser marking on a spec imen comprising a surface. Creating visible patterns on surfaces using laser irradiation is a rapidly growing technology with many applications in object identification, customization, and authentication [Liu et al. 2019] Laser marking is an environmentally friendly, low maintenance process with no con sumables, dyes, or pigments. While mostly a monochromatic method, some materials exhibit a range of colors when treated with laser, as a result of complex physicochemical phenom- ena. Among such materials are stainless steel and titanium, some of the most important in dustrially metals. Despite the great potential, the industrial adoption of color laser marking is almost non-existent due to its challenging characterization. In the absence of such a charac terization, the relationship between the device’s design space (laser parameters) and perfor mance space (e.g. marked colors) is unknown. This relationship is too complex to capture with physics-based methods [Nanai et al.1997]. Instead, the current practice finds design pa rameters that lead to “interesting colors” by trial and error measurements. These primary col ors are then used to mark simple motives and logos. This brute-force color gamut exploration scales poorly with the laser marking high-dimensional design space, resulting in neglecting some design parameters. The coloration of different substrates using laser irradiation is an active field of research with a long history [Birnbaum 1965] There are many color formation mechanisms employing dif ferent laser sources and different materials; see [Liu et al. 2019] for a recent review. Surface oxidation is one of these mechanisms where the heat (generated by a laser) facilitates the reaction of materials with oxygen. Oxidation-induced colors are believed to stem from multi layer, heterogeneous mixture of structural colors (based on thin-film interference) [Del Pino et al. 2004] on one hand, and the traditional pigment-based color of oxides [Langlade et al. 1998] on the other hand. Despite a handful of initial efforts [Veiko et al. 2013], predicting the structure and composition of oxide layers is extremely difficult due to the complex thermody namics of the laser marking process [Nanai et al. 1997] Even with known material composi tion, predicting the surface color requires a challenging light-matter interaction model most likely based on an electromagnetic simulation [Auzinger et al. 2018]

For some popular metals, such as stainless steel and titanium, oxidation-based color laser marking has been extensively studied. This spans a range of laser-marked metals’ behav iors, from electromechanical [Lawrence et al. 2013] to corrosion resistance properties [teRcka et al. 2016] Further related to this invention is a class of studies focused on the ef fect of various process parameters on the marked colors [Laakso et al. 2009] Most of these works [Adams et al. 2014; Antohczak et al. 2013, 2014] rely on sampling and marking the process parameter space uniformly. As laser marking is time and material consuming, these methods cannot cope with the dimensions of the design space and end up ignoring a large portion of process parameters. It is worth noting that some empirical color discovery methods [Veiko et al. 2016] try to find different laser parameters that lead to the same color. But that color needs to be known beforehand. Moreover, these methods are restricted to interference- based colors within a limited range of laser energy and a limited number of parameters.

Formulating design problems based on multi-objective optimization and solving them by computing the Pareto front is known in the field of computer graphics. Notable examples are simplifying procedural shaders [Sitthi-Amorn et al. 2011] or minimizing power consumption in real-time rendering [Wang et al. 2016] In computational fabrication, exploring the perfor mance space of a process has attracted recent attention. With the advent of 3D additive technologies, these efforts are mainly focused on exploring the mechanical properties of 3D printed microstructures. As an example, Schumacher et al. [2015] precompute microstruc tures’ performance space defined by mechanical metamaterial families in order to accelerate their heterogeneous topology optimization. They first populate the performance space by perturbing the initial designs (in the design space) and then fill the unpopulated regions of the performance space through either interpolation or inverse optimization. For a similar pur pose, Zhu et al. [2017] combines a discrete, random perturbation of designs near the gamut boundaries with a continuous optimization that further expands the gamut by refining existing designs. In a more general-purpose method, Schulz et al. [2018] further emphasizes the im portance of exploring the performance gamut’s hypersurface (or Pareto front) instead of its hypervolume. A Pareto front captures a set of solutions in the performance space that are compromising different, potentially conflicting objectives. Although these methods serve as important sources of inspiration, it is not possible to rely on any of them as they depend on closed form, smooth characterization functions that connect the design and performance spaces. For example, the method of Schulz et al. [2018] requires a smooth (twice differentia ble) forward characterization of the process and works only with continuous design parame ters.

The object of the invention is to provide a method for preparing a laser marking system to create a colored laser marking on a specimen and to a method for creating a colored laser marking on a specimen comprising a surface such that the color laser marking is equipped with a high level of versatility.

The problem is solved by the method according to claim 1 and the method according to claim 10. The further dependent claims provide preferred embodiments.

According to the invention a method for preparing a laser marking system to reproduce a la ser-marked color image on a specimen comprises the following steps: a) Providing a laser marking system and a specimen comprising a surface layer, wherein the laser marking system comprises a preset number of laser parameters; b) Performing an exploration of a first gamut specified by the laser marking system and the specimen comprising a surface layer including the following steps: aa) Creating a design space with a preset number of design points, wherein each design point comprises a combination of the preset number of laser parame ters; bb) Performing a marking of a sample on the specimen for each design point; cc) Measuring the sample using at least one detection device and determine for each design point a performance point, wherein the measured performance points define a performance space; dd) Evaluating the performance space with regard to preset performance criteria using an evaluation device, wherein a Pareto front is determined comprising a subset of performance points; ee) Generating an offspring design space with offspring design points; ff) Creating a first gamut using the subset of performance points forming the Pa reto front; wherein the steps bb) to dd) are iterated for a preset iteration number, wherein in each itera tion the offspring design space of the previous iteration is used in step bb), wherein in each iteration the performance space and the measured performance space are combined such that in step dd) the combined performance space is used. Preferably in each iteration the de sign space is combined with the design space of the previous iteration such that in steps dd), ee) and ff) a combined design space is used. Preferably the design space is initially popu lated with randomly chosen design points. Preferably the population size is in the range of 50 to 500, more preferably in the range of 75 to 250. Preferably the evaluation device is inte grated in a control unit which controls the devices of the laser marking system. Preferably the method is executed by the control unit completely or in part. Preferably the evaluation device and/or the control unit are a computer, a processor unit, or a similar device.

The method provides a device characterization which is the prerequisite for any color repro duction system including laser marking. In the absence of an analytical function that maps laser marking parameters to marked colors, a data-driven method according to the invention is performed. The method provides a black-box model of the process ruling out a physics- based prediction of the laser-induced composition of oxides. This invention introduces the first systematic color discovery algorithm for laser marking systems. The present method is a non-exhaustive performance space exploration of a laser marking system. Preferably the surface layer is a metallic surface layer. However, the present method may also be applied on surface layers or specimen which are made of non-metallic materials. The present method may be applied for any kind of laser system and any kind of specimen.

Preferably step dd) involves a multi-objective optimization. Unlike a typical optimization, multi-objective optimization problems are evaluated based on multiple criteria. Very often, these criteria are in conflict. In the present case for example, some marked colors may be saturated but leave thick traces and lower the resolution. Hence, instead of a single optimal solution, there exists a set of optimal solutions, known as Pareto optimal solutions or Pareto set. The projected Pareto set into the performance space is called Pareto front. A Pareto front captures a set of solutions in the performance space that are compromising different, potentially conflicting objectives. A member of the Pareto front is not dominated by any other point in the performance space in all performance criteria. In other words, it is more perfor- mant than all other points in at least one criterion. Preferably a dense set of Pareto-optimal solutions to the color laser marking problem with the above objectives is uncovered. Preferably step dd) employs non-dominated sorting genetic algorithm (NSGAII) which is a sorting algorithm based on the performance point’s presence in multi-level Pareto fronts.

According to a preferred embodiment the laser system comprises at least one pulsed laser and at least one scanning device. Preferably by the scanning device a laser spot is movable relative to the specimen. Alternatively, or cumulatively by the scanning device the specimen is movable relative to the laser spot. The laser spot is preferably an area of the laser beam impinging on the specimen. With other words, it is conceivable that the specimen is in a fixed position and the laser beam is moved relative to the specimen. The scanning device could preferably be a galvanometric scanner, a movable mirror, or a similar device by which the di rection of the reflected laser beam can be controlled. Alternatively, it could be conceivable that the laser beam is not moved, and the specimen is moved relative to the laser beam. The scanning device could therefore preferably be a x-y-stage or a x-y-z-stage. It could be also conceivable that the laser beam as well as the specimen are movable by a scanning device. Further, the laser and the scanning device(s) are preferably controlled by the control unit.

The control unit may advantageously also be connected with the at least one detection de vice in order to receive the measured data which are then preferably used by the evaluation device. It is also conceivable that more than one laser beam is used for the color laser mark ing. Depending on the used method to create a color laser marking the type of pulsed laser may be chosen. The pulsed laser could therefore preferably be a nano- second laser a pico second-laser or a femto-second laser.

According to a preferred embodiment a design point comprises at least one laser parameter selected form: the frequency of the laser pulses, the power of a laser pulse, the width of a la ser pulse, the speed of the laser beam relative to the specimen along a vector while marking, the line count, which defines the numbers of lines in a cluster representing the marked sam ple, the distance between the lines within a cluster representing the marked sample, the number of times a vector is marked. Preferably the dimensionality of the design space is set by the number of laser parameter represented by a design point. Thus, a design space could be 7 dimensional incase all seven of the aforesaid laser parameters are comprised in the de sign space. The design space may however be adjusted according to the specific needs. It could therefore be any number and any combination of the aforementioned laser parameters. It is conceivable that the design space comprises further parameters which might influence the formation of colors. Such further parameters depend on the actual used method for creat ing a color on the specimen. Preferably a design point comprises the parameter focal dis tance of the laser beam, type of medium gas, the ambient temperature. The medium gas is the ambient gas surrounding the specimen during the laser marking. This medium gas could be for example air. A preferred method for color marking is laser induced oxidation of the sur face layer of the specimen. In such a method the type of the ambient medium gas is im portant in view of presence of oxygen and the amount of the present oxygen. The design space may therefore preferably comprise the relevant parameters which might influence the resulting color in the laser marking process.

According to a further preferred embodiment the performance criteria in step dd) comprise at least one of: chromaticity, hue spread, resolution, performance space diversity, design space diversity, color repeatability, color uniformity. Preferably the performance criteria in step dd) comprise all of the aforesaid criteria. According to a further preferred embodiment the criteria color repeatability, color uniformity are pruned. It is conceivable that further performance cri teria are considered. The type, the number and the combination of the used performance cri teria are preferably adjusted on the used specimen the laser system and further influencing factors.

The chromaticity may be preferably described by: where a * and b * are the color coordinates of the CIELAB color space [Wyszecki and Stiles 1982] Marked colors with large chroma produce more saturated color images. Preferably the hue spread (f H s) ensures the presence of high-chromaticity colors at all hue angles. The res olution may preferably be evaluated by measuring the thickness of a line marked by a set of given laser parameters. The resolution (fp;) may be preferably described by: where t is the line thickness. This criteria is preferably due to the preferred use of a line- based halftoning method. The performance space diversity (fp SD ) is preferably measured for each point in the performance space as the reciprocal of the distance to its closest two neighbors and is given by: i ¹ j, 0 £ i,j < |P| where p is the point to be scored and P its respective population in performance space. The design space diversity (f D s D ) is preferably measured analogously as the performance space diversity except in the design space.

According to a preferred embodiment the performance criteria in step dd) comprise at least one of: chromaticity, resolution, performance space diversity, design space diversity. Prefera bly performance criteria in step dd) comprise all of the aforesaid criteria. Preferably the hue spread is not performance criteria. Preferably wherein the performance points are projected in to a CIECH space. It is advantageous that the CIECH space is split into a first number of circular sectors. Preferably the circular sectors en bloc form a hue wheel. The performance points within each circular sector of the hue wheel are advantageously evaluated regarding said performance criteria. Advantageously said evaluation is iterated for a preset iteration number, wherein in each iteration the number of sectors forming a hue wheel is altered. Con sequently, for each iteration the area of the sectors is different. Preferably within each circu lar sector, the performance points are evaluated using a non-dominated sorting algorithm based on all said performance criteria except the hue spread. It is advantageous that this evaluation is repeated each time with a randomly chosen number of sectors, and with a ran dom angular offset. After each iteration, every performance point is assigned to a potentially different Pareto front. At the end of the preset number of iterations, every single performance point is characterized by its front frequency vector that represents the frequency of its pres ence in the first front, second front and so on. Preferably each performance point is charac terized by a frequency vector, which represents the presence in a certain Pareto front. By this procedure a gamut with a balanced hue spread is achieved. Thus, the performance crite rion hue spread may be evaluated effectively by evaluating all remaining performance criteria multiple times.

According to a further preferred embodiment an additional evaluation regarding the achro matic properties of the performance points is performed by performing step b) using the per formance criteria in step dd) lightness, resolution, performance space diversity, design space diversity. Thus, in the resulting first gamut both, the chromatic and the achromatic axes are covered.

According to a further preferred embodiment the method further comprises the step: select ing a set of primary colors from the first gamut, wherein the selected primary colors form a second gamut. The extraction of the primary colors is concerned with selecting a set of col ors that generates the maximum color gamut through the use of the preferred marking method for example a halftoning method. While, unlike printers, the number of primaries is not strictly limited, a smaller number of primaries lead to improved marking time as they cause fewer switching delays of the laser system. Not all colors in the explored gamut can be considered for primary extraction. Thus, before primary extraction, the gamut is preferably pruned by excluding colors that: 1) don’t satisfy a specified resolution requirement, 2) reveal a low repeatability, and 3) exhibit non-uniformity.

According to a further preferred embodiment the data relating to the design space and the performance space of the first gamut is stored in a database. This has the advantage that for a given laser system and a certain type of specimen the method for preparing the laser mark ing system preferably has to be done once. For the same type of specimen the relevant data regarding the first gamut may be retrieved from the database before further steps regarding the marking are implemented.

The object of the present invention is also solved by a method for creating a colored laser mark on a specimen comprising a surface.

The method for creating a colored laser mark on a specimen comprising a surface layer may comprise the single features or combinations of the features described above for the method for preparing a laser marking system to create a colored laser and vice versa. Further, the same advantages may apply for the method for creating a colored laser mark on a specimen comprising a surface layer as described above for the method for preparing a laser marking system to create a colored laser and vice versa.

The method for reproducing a laser-marked color image on a specimen comprising a surface layer comprises the following steps: a) Verifying the database regarding data related to the first gamut and/or second gamut with regard to the type of the specimen and the laser marking system, wherein said data is obtained by the method for preparing a laser marking sys tem according to one of the embodiments mentioned above; b) Retrieving data related to the first gamut and/or second gamut from the data base or perform the method for preparing a laser marking system according one of the embodiments mentioned above; c) Providing an input image to be reproduced as laser marking on the specimen; d) Performing a color management workflow by which creates control data for the laser marking system derived from the input image; e) Perform the marking according to the control data. In step a) the database is searched if for the specific type of the specimen data related to the first gamut and/or second gamut is present. Such data is obtained by performing the method for preparing a laser marking system according to one of the above-mentioned embodi ments. Such a verification is preferably done by a control unit. The user provides the control unit the specification of the specimen in particular the specification of the surface layer of the specimen. Preferably the surface layer is a metallic surface layer. However, the present method may also be applied on surface layers or specimen which are made of non-metallic materials. The present method may be applied for any kind of laser system and any kind of specimen.

According to a preferred embodiment the laser system comprises at least one pulsed laser and at least one scanning device. Preferably by the scanning device a laser spot is movable relative to the specimen. Alternatively, or cumulatively by the scanning device the specimen is movable relative to the laser spot. The laser spot is preferably an area of the laser beam impinging on the specimen. With other words, it is conceivable that the specimen is in a fixed position and the laser beam is moved relative to the specimen. The scanning device could preferably be a galvanometric scanner, a movable mirror, or a similar device by which the di rection of the reflected laser beam can be controlled. Alternatively, it could be conceivable that the laser beam is not moved, and the specimen is moved relative to the laser beam. The scanning device could there-fore preferably be a x-y-stage or a x-y-z-stage. It could be also conceivable that the laser beam as well as the specimen are movable by a scanning device. Further, the laser and the scanning device(s) are preferably controlled by the control unit.

The control unit may advantageously also be connected with the at least one detection de vice in order to receive the measured data which are then preferably used by the evaluation device. It is also conceivable that more than one laser beam is used for the color laser mark ing. Depending on the used method to create a color laser marking the type of pulsed laser may be chosen. The pulsed laser could therefore preferably be a nano- second laser a pico second-laser or a femto-second laser.

According to a further preferred embodiment the color management workflow is a halftoning workflow. By using multi-color halftoning through a color reproduction workflow the reproduc tion of arbitrary images is enabled and not only uniform colors. Further a preview of the im ages before the marking is enabled. Therefore, for the laser marking a halftoning technique is implemented. By using a halftoning technique the visual impression of a continuous tone image is reproduced by taking advantage of the low-pass filtering property of the human vis ual system. The halftoning method aims at creating bilevel images conveying the visual illu sion of a continuous tone image. Groups of colored and white pixels are printed with certain ratio and structure so that, when viewed by the eye, give the impression of continuous color. In color halftoning, a given number of color layers are halftoned separately. The final color halftone is the result of the color mixing of different halftone layers by overlaying them on top of each other. Existing color halftoning methods for printers are for example clustered dot and blue noise dithering. Herby a halftone layer is created for each color separately. The final color-halftone image is formed by the superposition of all the layers, wherein the screen dot layers partially overlap. In the laser marking an overlap of the colors and the superposition of the layers is not preferred. Preferably, the color management workflow is a juxtaposed half toning workflow. Accordingly, it is preferred that for the laser marking a juxtaposed halftoning technique is implemented. The preferred juxtaposed halftoning technique relies advanta geously on discrete line geometry, which provides subpixel precision for creating discrete thick lines. Preferably the continuous tone color image is converted into a set of binary im ages each corresponding to a primary color. The binary images are synthesized in the form of lines and places them next to each other without overlapping.

According to a further preferred embodiment the color management workflow comprises the steps aa) Applying a forward color prediction model to construct a third gamut with re gard to the second gamut and the use of juxtaposed halftoning; bb) Mapping the input image into the third gamut; cc) Perform a color separation such that for each mapped color a corresponding area-coverage of each primary color is determined; dd) Binarize the area-coverages using the juxtaposed halftoning method and cre ate raster halftone images; ee) Convert the raster halftone images into vector data, wherein the control data comprise the said vector data. Preferably the control data is sent to the laser.

In step aa) the third gamut is preferably created from the second gamut and/or the first gamut under the condition a halftoning method is used. Preferably the third gamut is gener ated by halftoning a set of primary colors. The forward color prediction model predicts the color of several thousands of halftones spanning the space of the relative area of primaries in each halftone, known as area coverages. The third gamut surface is then fitted to this volu metric point cloud and is later used in step bb) for gamut mapping. Preferably the Yule-Niel- sen (YN) prediction model to predict the multi-color, juxtaposed halftones of laser primary colors. In the gamut mapping the color space of the input image is translated to the color space of the third gamut. The color space is typically be displayed as a volume of achievable colors. In step cc) the color separation builds on the forward prediction model to compute the particular primaries and their area coverages that reproduce a given color (inside the color gamut).

The method for creating a colored laser mark on a specimen comprising a metallic surface may comprise the single features or combinations of the features described above for the method for preparing a laser marking system to create a colored laser and vice versa

According to a further preferred embodiment the laser marking is based on laser induced oxi dation of the surface layer of the specimen. Preferably this applies to the method for creating a colored laser mark on a specimen comprising a metallic surface and/or the method for pre paring a laser marking system to create a colored laser. In this approach the laser induces heating which leads to formation of a transparent or semitransparent oxide film on the sur face of the specimen. A with light illumination can be reflected from the top and bottom sur faces of the oxide film. A constructive interference of the reflected beams makes the surface appear a certain color, which is determined by film thickness, refractive index of the oxide, and the order of interference [Liu et al.2019].

Preferably, the laser marking is based on laser induced structuring of the surface layer of the specimen. Preferably this applies to the method for creating a colored laser mark on a speci men comprising a metallic surface and/or the method for preparing a laser marking system to create a colored laser. In this approach laser induced periodic surface structures (LIPSSs) are produced on the surface of the specimen by the laser system. The LIPSS act as a grat ing to give rise to iridescent colors due to the optical diffraction effect. The colors are not caused by pigments but originate from material surface micro/nanostructures, namely, struc tural colors.

Preferably, the laser marking is based on the laser induced generation of micro/nanoparticles on the surface layer of the specimen. Preferably this applies to the method for creating a col ored laser mark on a specimen comprising a metallic surface and/or the method for prepar ing a laser marking system to create a colored laser. In this approach surface structures are induced by the laser system. The surface structures which excite the surface colors are ran domly distributed without regularity, and the color does not vary with the viewing angle. Sur face plasmon resonance (SPR) effects arising from metallic nanostructures and nanoparti cles are the main causes for this type of coloring.

Preferably the laser marking is based on laser induced plasmonic colors on metals. Metal na noparticles exhibit scattering properties due to excited plasmons that depend on their shape, size, composition and the host medium. There are various techniques known which render plasmonic colors including laser interference lithography. By plasmonic colors precious met als such as gold and silver and also metals like copper and aluminum may be marked.

The method for creating a colored laser mark on a specimen comprising a surface layer and/or the method for preparing a laser marking system to create a colored laser marking may however also be applied in combination with other laser induced color marking methods. Other mechanism may enable marking on a wide range of metals and even non-metals. Due to the fact that the actual marking process is preferably treated as a black box and the method is a data driven method it is adaptable to various other laser induced color marking methods.

Preferably the specimen comprises at least a top surface layer made of a metal on which the laser marking is performed. It is also conceivable that the specimen is made completely of a metal. Advantageously the surface layer and/or the entire specimen is made of stainless- steel titanium or a similar metal. Preferably this applies to the method for creating a colored laser mark on a specimen comprising a metallic surface and/or the method for preparing a laser marking system to create a colored laser.

Further advantages, aims and properties of the present invention will be described by way of the appended drawings and the following description.

In the drawings:

Fig. 1 shows diagrams regarding the color response of the laser;

Fig. 2 shows a method for preparing a laser marking system (100) to create a col ored laser mark on a specimen;

Fig. 3 shows examples of different hue-wheel configurations used by the method;

Fig. 4 shows a method for creating a colored laser mark on a specimen;

Fig. 5 shows a schematic color reproduction workflow;

Fig. 6 shows a laser marking system; Fig. 7 shows a visualization of laser parameters;

Fig. 8 shows a color gamut evolution of a full exploration on AISI 304 stainless steel;

Fig. 9 shows a color gamut evolution of a random exploration on AISI 304 stainless steel;

Fig. 10 shows explored gamuts with different configurations on AISI 304 stainless steel;

Fig. 11a shows the average thicknesses over iterations with and without fp;

Fig. 11b shows the average reproduction error of the Yule-Nielsen model for different n-values;

Fig. 11c shows the extraction of chromatic primaries;

Fig. 12 shows multiple examples of original, gamut mapped and marked images;

Fig. 13 shows a comparison of the marking parameters from Antohczak et al. with the marking of the present laser marking system;

Fig. 14 shows marked images on AISI 304 and AISI 430;

Fig. 15 shows a color gamut evolution of a full exploration on AISI 430 stainless steel;

Fig. 16. shows a painting of Maria de’ Medici by Alessandro Allori, marked on AISI

304;

Fig. 17 shows laser marked images on stainless steel using the method of one em bodiment of this invention.

This invention provides means to equip color laser marking with the same level of versatility found in color printers. Assuming a blackbox model of the difficult device characterization, a measurement-based, data-driven performance space exploration is designed. Different per formance criteria are explored including the color gamut and marking resolution by consecu tive marking and measuring. For this, the process’s Pareto front is uncovered by formulating a multiobjective optimization problem and solving it using an evolutionary method augmented by a Monte-Carlo approach. The optimization explores the hidden corners of the 7- dimen sional design space in search of useful parameters that lead to a dense set of diverse, high- resolution colors. This invention goes far beyond the state of the art color image marking by introducing a complete color management workflow that takes an input image and laser- marks the closest approximation on metal surfaces. The color reproduction workflow adopts the principles of halftone-based color printing. It extracts a number of primary colors from the explored gamut and reproduces input colors by juxtaposing the extracted primaries next to each other in a controlled manner. The fabricated color images enjoy high resolution, intro duce no significant artifact, and demonstrate accurate color reproduction. The invention pro vides therefore a discovery method that automatically finds the desired design parameters of a black-box fabrication system and the first color-image reproduction workflow for laser mark ing on metals.

Device characterization is the prerequisite for any color reproduction system including laser marking. In the absence of an analytical function that maps laser marking parameters to marked colors, one must rely on data-driven methods. In a first attempt, one can sample the design space, mark and measure the sampled design points, and construct a look-up table. This exhaustive strategy is subject to the curse of dimensionality given the relatively large number of parameters involved in color laser marking. The fact that function evaluations re quire actual marking and measuring further slows down the process. Additionally, the non smooth color response to laser parameters renders interpolation schemes ineffective. This is shown in Figure 1 where color coordinates of marked patches may change abruptly in re sponse to marking parameters. It is contrasted with the smooth response of a typical printer to its control parameters. In Figure 1 the top left graph shows the color response of the laser vs. that of a typical printer. CIE L*, a* and b* values are plotted in red, green and blue re spectively. The laser-marked colors (on stainless steel AISI 304) are repeated three times. Their average (solid lines) and standard deviation (shaded region) are shown. The non smooth behavior of the laser-marked colors is apparent.

Figure 2 shows a method 1 for preparing a laser marking system 100 to reproduce a laser- marked color image on a specimen comprising the following steps: a) Providing a laser marking system 100 and a specimen 105 comprising a metallic sur face layer 105a, wherein the laser marking system comprises a preset number of laser parameters 12; b) Performing an exploration of a first gamut 2 specified by the laser marking system 100 and the specimen 105 comprising a metallic surface layer 105a including the following steps: aa) Creating 3 a design space 10 with a preset number of design points 11, wherein each design point 11 represents a combination of the preset number of laser pa rameters 12; bb) Performing 4 a marking of a sample on the specimen 105 for each design point 11; cc) Measuring 5 the sample using at least one detection device 106 and determine for each design point a performance point 14, wherein the measured performance points 14 define a performance space 13; dd) Evaluating 6 the performance space 13 with regard to preset performance criteria using an evaluation device, wherein a Pareto front is determined comprising a subset of performance points; ee) Generating 7 an offspring design space 10a with offspring design points 11a; ff) Creating 8 a first gamut 2 using the subset of performance points forming the Pa reto front; wherein the steps bb) to dd) are iterated 9 for a preset iteration number, wherein in each iter ation 9 the offspring design space 10a of the previous iteration is used in step bb), wherein in each iteration the measured performance space is combined 15 with the performance space of the previous iteration 9 such that in step dd) the combined performance space 13a is used. Preferably in each iteration the design space 10, 10a is combined with the design space 10, 10a of the previous iteration 9 such that in steps dd), ee) and ff) a combined de sign space 10b is used. Preferably the design space is initially populated with randomly cho sen design points. Preferably the population size is in the range of 50 to 500, more preferably in the range of 75 to 250. Preferably the evaluation device is a computer, a processor unit, or a similar device. Preferably the method 1 is executed completely or on part by a control unit. The control unit might be a processor, a computer or a similar device.

The laser system 100 is depicted in Figure 6 and comprises a preferably pulsed laser 101 and a scanning device 103, 104. According to one embodiment, the scanning device 103,

104 moves the laser spot relative to the specimen 105 or on the surface 105a of the speci men 105. The specimen 105 is therefore in a fixed position. The scanning device 103 could preferably be a galvanometric scanner, a movable mirror or a similar device by which the di rection of the reflected laser beam can be controlled. It is also conceivable that the specimen

105 is movable relative to the laser spot. The scanning device could therefore preferably be a x-y-stage or a x-y-z-stage. It could be also conceivable that the laser beam as well as the specimen are movable by a scanning device. Further, the laser 101 and the scanning de vice^) 103, 104 are controlled by the control unit 108. The control unit 108 may also be con nected with the at least one detection device 106 in order to receive the measured data which are then used by the evaluation device 107.

A design point 11 , 11a comprises at least one laser parameter 12 selected form: the fre quency of the laser pulses, the power of a laser pulse, the width of a laser pulse, the speed of the laser beam relative to the specimen along a vector while marking, the line count, which defines the numbers of lines in a cluster representing the marked sample, the distance be tween the lines within a cluster representing the marked sample, the number of times a vec tor is marked. A design point (11 , 11a) may further comprise the parameter focal distance of the laser beam, type of medium gas (in the present case air), ambient temperature. Gener ally, the design space might comprise relevant parameters which influence the formation of a color on a specific specimen. The performance criteria in step dd) comprise at least one of: chromaticity, hue spread, resolution, performance space diversity, design space diversity, color repeatability, color uniformity. In the present case the criteria color repeatability and color uniformity are, however, pruned.

The method 1 for preparing a laser marking system 100 provides a non-exhaustive perfor mance space 13 exploration of the laser marking system 100. Qualitatively speaking, the performance criteria favors diverse, saturated and high-resolution colors: the fundamental re quirements for color images. For solving this problem a multi-objective optimization is casted. Unlike a typical optimization, multi-objective optimization problems are evaluated based on multiple criteria. Very often, these criteria are in conflict. In the present case, for example, some marked colors may be saturated but leave thick traces and lower the resolution.

Hence, instead of a single optimal solution, there exists a set of optimal solutions, known as Pareto optimal solutions or Pareto set. The projected Pareto set into the performance space is called Pareto front. A member of the Pareto front is not dominated by any other point in the performance space in all criteria. In other words, it is more performant than all other points in at least one criterion. The goal is to uncover a dense set of Pareto-optimal solutions to the color laser marking problem with the above objectives. To this end, a multi-objective evolu tionary method is adopted, which is a successful tool for finding Pareto optimal solutions [Fonseca et al. 1993] The method, called non-dominated sorting genetic algorithm (NSGAII) [Deb et al. 2002] is well suited to our model-free characterization function, with both discrete and continuous parameters. At the heart of this method, is a evaluation or sorting method based on the members’ presence in multi-level Pareto fronts. The NSGA-II non-dominated sorting is insufficient for the present specific problem due to our hue diversity objective. Thus, a preferable Monte-Carlo approach is considered on top of the non-dominated evaluation and introduce a new evaluation method based on front frequencies. This method is called Monte-Carlo, multi-objective, genetic algorithm or MCMOGA for short.

In order to preferably adopt halftoning for color reproduction a set of primary colors that co vers both achromatic and chromatic axes is favorable. In a divide and conquer strategy, the- chromatic and achromatic (black and white) explorations may be separated, starting with ex plaining the former. High chromatic performance requires saturated colors corresponding to larger radii in the CIECH color space shown in Figure 3 [Schanda 2007] Furthermore, it re quires colors that span a range of different hues. Such colors mixed with black and white (through halftoning) generate a highly populated color gamut (2) that can be utilized for color image marking. By the method 1 for preparing a laser marking system 100 using a multi-ob jective optimization the performance criteria are maximized:

(1) Chromaticity, as marked colors with large chroma produce more saturated color im ages. f c ( a * , b * ) = sja * 2 + b *2 where a * and b * are the color coordinates of the CIELAB color space [Wyszecki and Stiles 1982]

(2) Hue spread (f H s) ensures the presence of high-chromaticity colors at all hue angles.

(3) Resolution; since we use a line-based halftoning, this criterion is evaluated by measur ing the thickness of a line marked by a set of given laser parameters where t is the line thickness.

(4) Performance space diversity (PSD); it is measured for each performance point (14) in the performance space (13) as the reciprocal of the distance to its closest two neigh bors i ¹ j, 0 £ i,j < |P| where p is the point to be scored and P its respective population in performance space.

(5) Design space diversity (f DSD ); it is measured as the performance space diversity except in the design space.

The performance criteria (1) to (4) are measured in the performance space while the perfor mance criterion (5) is measured in the design space. The performance criteria (1) to (3) are the qualities which are directly sought from laser marked images. The performance criterion (4) improves the convergence rate and criterion (5) helps avoiding local extrema by promot ing solo points in the performance space 14.

The method 1, navigates the laser’s design space 10, 10a, 10b in directions that lead to a dense Pareto set, i.e. , the set of designs (laser parameters) that improves the above perfor mance criteria. It is started with a random population in the design space 10, mark it and measure its performance, then iteratively evolve it into a larger population with as many Pa reto optimal solutions as possible. At each iteration, represented schematically in Figure 2, the Pareto set is promoted of its population to be passed along to the next iterations 9 using the genetic algorithm. The iterations are stopped when no significant improvement in the Pa reto front is observed any longer.

As in almost any genetic method, a fitness measure should be assigned to each member of the population. Fitter solutions are selected and used to create the next generation. Given the difficulty of assigning a single fitness value to multi-criterion objectives, the non-domi- nated evaluation method [Deb et al. 2002] evaluates the members of a population according to their presence in multi-level Pareto fronts. It starts with finding the first non-dominated front, i.e., all solutions or performance points 14 in a population that belong to the Pareto front. This is done by comparing each performance points’ 14 performance objective by ob jective to every other performance point 14 in the population. If a performance point 14 is more performant than all other performance points 14 in at least one criterion, it is labeled as a first-front performance point 14. The second non-dominated front is computed by temporar ily discarding the first front and repeating the above procedure. This procedure is continued until all members of the population are labeled with their respective fronts. This results in a number of disjoint subsets making up the whole population, each with its front label. Note that, in the spirit of the Pareto concept, members inside the same front are not sorted.

Figure 2 shows one iteration of the method 1. The method takes the starting population at it eration i (Pi) and generates an offspring generation Qi using the genetic method. Marking and measuring the design space (DS) yields the corresponding performance points 14 points in the performance space 13 (PS). Pi and Qi are combined into Ri and evaluated using the proposed method. Ri is added to the first gamut 2 and its fitter half of Pi+1 is passed as the starting generation to the next iteration 9.

Considering the hue spread objective the method is not favorable. The hue spread criterion helps the color gamut grow in all angular directions in a balanced manner. Without the hue spread criterion, the method may explore some specific hue angles more than others result ing in a non-uniform growth of the chromaticity gamut. For achieving a gamut 2 with bal anced hue spread a single solution cannot be evaluated but rather in combination with other solutions. This can quickly lead to nontrivial computation: for 10 angular samples in a popula tion of 200, (^) ³ 10 16 evaluations.

This combinatorial explosion can be avoided by resorting to a Monte-Carlo method. The per formance criteria in step dd) comprise at least one of: chromaticity, resolution, performance space diversity, design space diversity. Preferably all of said performance criteria are used. The performance points 14 are projected in to a CIECH space, wherein a the CIECH space is split into a first number of circular sectors 15 forming a hue wheel 16. Thus, the hue wheel 16 splits the CIECH space into a random number of circular sectors 15 (Figure 3). Within each sector 15, the performance points 14 are evaluated using the described non-dominated sorting algorithm based on all performance criteria except the hue spread. Said evaluation is iterated for a preset iteration number, wherein in each iteration the number of sectors 15 forming a hue wheel 16 is altered. Thus, the procedure is repeated each time with a ran domly chosen number of sectors 15, and with a random angular offset. After each turn of the hue wheel 16 or each iteration, every individual performance point 14 is assigned to a poten tially different front. At the end of this loop, every single performance point 14 is character ized by its front frequency vector that represents the frequency of its presence in the first front, second front and so on. The iteration is stopped when the change in front frequencies is below a certain threshold. The population of the performance space 13 is sorted based on the frequency of their "top" fronts where a single first front is worth more than any number of second fronts.

This procedure is schematically shown in Figure 3 in which multiple turns of the hue wheel 16, with different number of sectors 15 and angular offsets, ensure all performance points are evaluated in different configurations and are ranked in a proper way. In Figure 2, for ex ample, it is easy to see that some points may not be sorted properly using a single hue- wheel configuration. In Figure 3 Examples of different hue-wheel 16 configurations (left) used in the described method. First, the points within each circular sector of each hue wheel 16 are assigned to front labels (encircled numbers next to each point) using the conventional NDS method. Next, all front labels for each point are counted to form the front frequencies (right). As an example, the point shown with a star has been assigned two times to the first front and two times to the third front. Notice that, for the same point, only the first two hue- wheels 16 would not suffice as it would have been assigned only to the third front despite high potential for improving the gamut 2. Further, an additional evaluation regarding the achromatic properties of the performance points 14 is performed by performing step b) using the performance criteria in step dd) light ness, resolution, performance space diversity, design space diversity. Thus, in the chromatic exploration, the lightness values (CIE L*) are discarded. Two separate explorations are per formed for black and white colors on the lightness axis. For the black colors, minimize the chromaticity criteria is minimized, thereby encouraging low chromaticity colors. Additionally, Hue spread performance criteria is replaced with a lightness minimization. Exploring white colors is the same as the black colors except the lightness performance criteria is maxim ized.

After performing the method 1 for preparing a laser marking system 100 the performance space 13 of the laser marking system 100 is explored. Then the performance space 13 can be exploited for color image reproduction. This is done in this embodiment by adopt the prin ciples of halftone-based color printing for color laser marking. For this, a number of primary colors is found that meet the resolution requirement and produce the largest second color gamut. Afterwards a color management workflow is built that takes input color images and marks the closest approximation using the selected color primaries.

Thus, the method 1 may further comprise the step selecting a set of primary colors from the first gamut 2. The selected primary colors form then a second gamut. The primary extraction is concerned with selecting a set of colors that generates the maximum color gamut through halftoning. While, unlike printers, the number of primaries is not strictly limited, a smaller number of primaries lead to improved marking time as they cause fewer switching delays of the laser. Not all colors in the explored first gamut 2 can be considered for primary extraction. Thus, before primary extraction, the first gamut 2 is pruned by excluding colors that: 1) don’t satisfy the specified resolution requirement, 2) reveal low repeatability, and 3) exhibit non uniformity.

Similar to the gamut exploration described above the achromatic and the chromatic primaries are extracted separately. First the explored chromaticity gamut of the laser marking system composed of a discrete set of colors is explored. The convex hull of this set in the CIExy chromaticity space is found where x = X/(X + Y + Z), y = Y/(X + Y + Z) [Wyszecki and Stiles 1982] The reason for applying the convex hull in the CIExy is that, unlike CIEa*b* or CIECH, it is a linear space under halftoning. Colors inside the convex hull can be reproduced through halftoning with high accuracy. The colors in the convex set give the largest area and there fore the largest chromatic gamut. In order to reduce the number of primaries, those members of the convex set that don’t contribute to the gamut area significantly may further be ex cluded. The achromatic primary extraction selects the darkest and the brightest colors with negligible chromaticity from within the black and white explored gamuts, respectively.

The data relating to the design space 10, 10a, 10b and the performance space 13, 13a of the first gamut 2 and/or data related to the second gamut are stored in a database 109. The da tabase 109 may be connected to the control unit 108 and/or the evaluation unit 107.

The present invention comprises also a method 20 for creating a colored laser mark on a specimen 105 comprising a metallic surface 105a comprising the following steps: a) Verifying 21 the database 109 regarding data related to the first gamut 2 and/or second gamut with regard to the type of the specimen 2 and the laser marking system 100, wherein said data is obtained by the method 1 for preparing a laser marking system 100 according previous described embodiments; b) Retrieving 22 data related to the first gamut 2 and/or second gamut from the data base 109 or perform 23 the method 1 for preparing a laser marking system 100 accord ing to one of the previous described embodiments; c) Providing 24 an input image 27 to be reproduced as laser marking on the speci men 105; d) Performing 25 a color management workflow 28 which creates control data for the laser marking system derived from the input image 27; e) Perform 26 the marking according to the control data.

In step a) the database 109 is searched for the specific type of the specimen data related to the first gamut 2 and/or second gamut is present. Such data is obtained by performing the method 1 for preparing a laser marking system 100 according to one of the above-mentioned embodiments. Such a verification is preferably done by the control unit 108. The user pro vides the control unit 108 the specification of the specimen 105, in particular the specification of the surface layer 105a of the specimen 105. The retrieved data comprises the data regard ing the first and/or the second gamut which is matched to the used laser system and the spe cific specimen 105. The method is schematically shown in Figure 4. Further step d) is prefer ably performed by the control unit 108. In step e) the control unit 109 controls the laser 101 and the scanning device(s) 103, 104 accordingly.

Preferably the color management workflow 18 is a juxtaposed halftoning workflow. Accord ingly said method 20, wherein the color management workflow 28a) comprises the steps: aa) Applying 29 a forward color prediction model to construct a third gamut with regard to the second gamut and the use of juxtaposed halftoning; bb) Mapping 30 the input image 27 into the third gamut; cc) Perform 31 a color separation such that for each mapped color a correspond ing area coverage of each primary color is determined; dd) Binarize 32 the area-coverages using the juxtaposed halftoning method and create 33 raster halftone images; ee) Convert 34 the raster halftone images into vector data, wherein the control data comprise the said vector data.

The control data is then sent to the laser and step e) 28 may be performed.

A color management workflow ensures color reproducibility across different imaging devices. A real strength of the current method is to enable reproduction of arbitrary images and not only uniform colors using multi-color halftoning through a color reproduction workflow. It also enables a preview of the images before marking. The classic example is printing where the input images, from a camera for example, are printed as accurately as possible. Figure 5 sketches the color reproduction workflow for color laser marking. Given an input color in a given color space, e.g., sRGB, its reproducibility is ensured by mapping it into the color gamut of laser marking. The color separation computes the coverage of different laser pri mary colors which, when placed next to each other through halftoning, reproduce the input color. A typical printer’s color reproduction workflow generates different colors by spatial blending and superposition of multiple inks. Should such a workflow be imitated for the laser marking process, it needs to be ensured that both laser primary colors and their superposi tions are optimal. Exploring the design space for such an unlikely combination is significantly more difficult. Instead, different primary colors are placed strictly next to each other. This re sults in a considerably simpler exploration where only for a set of suitable primaries (and not their superpositions) is searched. In order to establish a color management workflow, juxta posed halftones of extracted primary colors are synthesized. The integral color of multi-pri mary halftones is predicted. This prediction model is numerically inverted in order to map the input colors into primary halftones.

In figure 5. The color reproduction workflow 28, 28a is depicted. An input image 27 is mapped to the second gamut of the laser marking system 100. In the color separation step 31, for each mapped color, the corresponding area coverages of each primary is computed (creating the gamut and color separation are built upon a color prediction model). The contin uous area-coverages are binarized 32 and placed next to each other using a juxtaposed half toning method 33. The raster halftone images are converted into vectors 34. Color halftoning converts a continuous tone color image into a set of binary images, each corresponding to one of the printer’s inks. The discrete-line juxtaposed halftoning [Babaei and Hersch 2012] synthesizes these binary images in the form of lines and places them next to each other without overlapping. In the original method designed for bitmap printers, using digital lines [Reveilles 1995] allows for subpixel thickness, low computational complexity, and, importantly to us, continuity. A continuous laser path ensures less switching delays, and therefore, faster marking with lower graininess caused by the two ends of each marked vec tor. As the original juxtaposed halftoning is designed for raster devices, the resulting raster images need to be transformed into vector representation suitable for our laser device. For this purpose, a naive line (a discrete line with unit thickness) as a mask and slide it on each halftone layer corresponding to each laser primary is used (Figure 5). This produces a list of vectors of different primaries which span the image plane and are sent to the laser device for marking.

The color prediction model has two roles in the color management workflow. First, it con structs the third color gamut generated by halftoning a set of primaries. It predicts the color of several thousands of halftones spanning the space of the relative area of primaries in each halftone, known as area coverages. The gamut surface is then fitted to this volumetric point cloud and is later used for gamut mapping 30. Second, the forward model is used in the color separation step 31 that computes the area coverages of the primaries for any input colors to be reproduced.

The Yule-Nielsen (YN) prediction model is used to predict the multi-color, juxtaposed half tones of laser primaries. The Yule-Nielsen equation [Yule and Nielsen 1951] predicts the Cl EX color coordinate (X t ) of a juxtaposed halftone as: where X is the Cl EX value of the i-th primary, and a, is its area coverage. The same equation applies for predicting Cl EX and CIEY color coordinates. The exponent n, called the Yule- Nielsen n-value is a tuning parameter.

Color separation builds on the forward prediction model to compute the particular primaries and their area coverages that reproduce a given color (inside the third color gamut). As the YN model is not analytically invertible, color separation is carried out by optimization: argminA Eoo ( Lab (Y/V(a)), c) a IMIi = 1 , a e [O, ΐR where c is the target color in the CIELAB color space and a is the optimization variable, i.e. the vector of area coverages of q primaries. As the CIEDE2000 color-difference formula [Sharma et al. 2005] is used for the distance metric, the modeled color using the YN model (YN(a)) should be converted to CIELAB from CIEXYZ (denoted by function Lab in the equa tion above. This equation searches for an area coverage vector that, after being marked, re sults in the minimum distance to the target color. As different primaries are juxtaposed, their relative area coverages should sum up to 1 and be non-negative.

The laser marking is based on laser induced oxidation of the surface layer (105a) of the specimen (105) or laser induced structuring of the surface layer (105a) of the specimen (105) or the laser induced generation of micro/nanoparticles on the surface layer (105a) of the specimen (105).

In the following different analyses and evaluations of both gamut exploration and image re production are presented. The laser marking system 100 is depicted in figure 6. In preferable experimental setup hardware may be used as described in the following. The laser marking device 101 comprises the main components in the form of a ytterbium fiber laser system (IPG Photonics YLPM- 1-4x200-20-20) and a galvanometric scanner (Scanlab InterScan III 10). The laser system (20 W, 1064 nm) generates a laser beam which is redirected by the scanning devices’s 103 Galvo Mirror system to any desired laser spot on the specimen. Equipped with an infrared F-Theta lens (f= 163 mm), the scanning device 103 is capable of imaging a planar field of 116 c 116 mm. An air filtering system blocks small particles from spreading in the room. In most of the experiments a 1 mm thick stainless steel type 1.4301 V2A (AISI 304) as specimen 105. Color laser marking is also possible on titanium.

In Figure 7 seven laser marking parameters are depicted including:

(1) Frequency: Defines the number of laser shots per second (1.6-1000 kHz, 100 Hz steps),

(2) Power: Adjusts the output power per shot (0-100%, 256 steps),

(3) Pulse width: Defines the duration of a single shot (4, 8, 14, 20, 30, 50, 100, 200 ns), and scanning parameters that forms a line cluster with properties:

(4) Speed: Defines the travel speed along a vector while marking (0-2000 mm/s, 1 mm/s steps),

(5) Line count: Defines the number of lines in a cluster (1-20 lines, 1 line steps),

(6) Hatching: Defines the distance between lines within a cluster (1-15 pm, 1 pm steps), (7) Pass count: Indicates the number of times a vector is marked (1-10 passes, 1 pass steps).

Figure 7 shows a visualization of laser (left) and scanning (right) parameters. The multipass, line cluster in the diagram on right forms the final color. It is worth noting that, due to tech nical limitations of the laser source, laser parameters cannot be used at arbitrary combina tions. Furthermore, it is not possible to vary laser parameters on the fly. For example, switch ing frequency and power takes 0.6 ms and 3 ms respectively; changing the pulse width takes about 2 seconds as it requires reestablishing the connection between the controller board and the laser. In our path planning, we therefore allow switching delays after changing these parameters ensuring the laser source can properly adapt to the new parameters.

For each point in the design space 10, 10a, 10b, its performance points 14 are measured in order to decide how to use that point in our exploration framework. All performance criteria, apart from the design space diversity, can be evaluated by measuring the thickness of a marked line cluster and the color of a marked patch. This is performed in two stages. First, for measuring the cluster’s thickness, a first detection device 106 in the form of a hand-held digital microscope (Reflecta DigiMicroscope USB 200) is used. In a second step, the thick ness of a given cluster is used to mark its corresponding patch by juxtaposing multiple clus ters within the desired area. Both hue and chromaticity, the pillars of the performance space 13 exploration, are computed from CIELAB, a perceptual color space. Therefore, a colorimet ric calibration [Hong et al. 2001] for measuring the color of marked patches is performed. The colorimetric calibration connects camera RGB signals to CIEXYZ coordinates through a form of regression. The CIEXYZ values then can be converted to the CIELAB coordinates using a set of well-known, analytical transformations [Wyszecki and Stiles 1982] For training the re gression, 121 printed color patches are used, with known spectra measured with an X-Rite i7 spectrophotometer, and obtain the ground-truth CIEXYZ values assuming D65 illumination. The same printed patches are captured with a second detection device 106 in form of Nikon D750 DSLR camera (with macro lens Tamron SP 90mm F/2.8 Di) obtaining raw RGB signals that have been corrected for spatial and temporal light fluctuations. It needs to be pointed out that this setup is a possible experimental setup. As already pointed out the at least one de tection device 106 could be connected to and controlled by the control unit 108. The above described measurements could then be performed automatically.

The colorimetric calibration shows high accuracy on a test set of 16 printed patches with an average DE 0 o = 2.26 and maximum 5.00. This calibration is therefore used to estimate the CIELAB color of marked patches. The structural nature of oxide colors causes a significant change in their appearance depending on the viewing and illumination geometry. It is ob served that laser-marked colors appear most saturated at specular and near-specular geom etries. Therefore, inspired by previous work on metallic prints [Pjanic and Hersch 2013], the color reproduction is confined to non-diffuse geometries. To this end, the stainless steel sub strate is illuminated with a large, diffuse area-light tilted approximately 45° from the sub strate’s normal and captured with the camera with a similar angle.

For evaluating the proposed gamut exploration algorithm, multiple runs are performed while discarding different objectives during different runs to show the objective’s effect on the ex ploration behavior. For a fair comparison, it is always started with the same randomly gener ated initial population. All generations have the same population size of 100. For the hue wheel, the random number of circular sections is limited between 4 and 72 while the random angular offset a is between 0 and 2TT. The stopping threshold for MCMOS is set to 0.001%.

In Figure 8 color gamut evolution of a full exploration (with fc, fHs, fR, fpsD, fosD) on AISI 304 stainless steel is shown. It demonstrates the evolution of the chromatic gamut, in hue- chroma polar diagram, when optimizing all performance criteria (referred to as the full explo ration). Overall, a decent evolution of colors with a symmetric, dense color gamut is seen. In terestingly, the purple to red regions are populated with a considerable delay, suggesting that some colors are more challenging to find than others.

In Figure 9 a color gamut evolution of a random exploration on AISI 304 stainless steel is shown. Compared to the full exploration (Figure 8), random marking does not lead to ade quate gamut growth. As the random exploration does not include the resolution objective, it is more illustrative to compare Figure 9d to Figure 10e as they both feature the same number of samples and none of them includes the resolution objective. The stagnant behavior of ran dom marking over time (Figure 9) and a lack of systematic resolution enhancement suggest that a very large number of samples is required to match the full gamut generated by the method of this invention.

In Figure 10 explored gamuts with different configurations on AISI 304 stainless steel are shown. In order to evaluate the effectiveness of the Monte-Carlo hue wheel method, two sim ilar explorations were run where the only difference is that the hue-spread objective f H s is en abled in one (Figure 10a) and disabled in the second (Figure 10d). It is observed that the MC approach promotes the hue diversity resulting in a symmetric color gamut. Ignoring the Monte- Carlo method introduces a bias toward areas with high chromaticity. In Figure 10d, for example, since the initial population (shown in Figure 9a) has a large number of chro matic yellow members, this area is emphasized during the exploration.

Marking high-quality images requires a set of diverse, saturated colors which are placed next to each other at a high spatial resolution. This criterion is defined by fp where design parame ters that mark thin line clusters encouraged. A comparison of two explorations with equal number of iterations, one with (Figure 10b) and another without (Figure 10e) the thickness minimization reveals that this objective slows down the color gamut growth and the overall gamut area by disfavoring saturated but thick colors. Crucially, however, it generates a denser gamut at lower thicknesses, visible when comparing gamuts that include only colors with small thicknesses (Figures 10c and 10f). A dense color gamut is very important during primary pruning (Section 4.1). Furthermore, Figure 11a shows the average thicknesses of the whole population at each iteration. The average thicknesses over iterations with fp is de picted in a continuous line. The dashed line represents the average thicknesses over itera- tions without f R . (only t < 80 pm were considered). Unlike the exploration without thickness objective, the full exploration shows a steady decrease in the marked line thicknesses.

The proposed color reproduction pipeline is evaluated by quantitative analysis and also a va riety of full-color marked images. After the pruning step a total of 6 primary colors is obtained including a black and white primary. Four chromatic primaries are shown in Figure 11c. Their parameters are reported in the following table:

In the pruning stage, the spatial and temporal repeatability is checked by marking the candi date primaries at four different locations of the substrate and compare their colors pairwise. Primaries having an average DEoo higher than 4 among all comparisons are discarded. Also for the progressive primary discarding, the number of primaries is reduced until the gamut area drops by more than 10%. For resolution pruning, the colors with thickness 40 ± 5 pm are kept. Additionally, colors with thickness around 20 pm are considered as juxtaposing two of them results in the target resolution.

For testing the accuracy of the Yule-Nielsen model, the primaries and also 92 test patches are marked with diverse area coverages of primaries. The resulting average DE 0 o error is 2.25 (Std=0.96, Min=0.50, Max=4.26) that demonstrates the high accuracy of the forward model. In Figure 11b the average reproduction error of the Yule-Nielsen model for different n- values is depicted. It is shown that the n-value equal to 1 works very well for the configura tion, reducing the present model to the widely known Neugebauer model [Rolleston and Bal- asubramanian 1993] There are a handful of physical and empirical interpretations of the Yule-Nielsen n-value in literature [Lewandowski et al. 2006] In the classic ink-on-paper prints, it accounts for the optical dot gain due to the lateral propagation of light inside the sub strate [Hebert 2014] Babaei and Hersch [2015] showed that this parameter is responsible for shadowing and masking in metallic-ink halftones. From the optimal n-value for our setup it can be inferred that, as expected, the subsurface scattering in metal is very negligible. Fur thermore, the marked primaries are very well leveled on the surface and cause no shadow ing or masking.

Figure 12 presents different results generated by the proposed image reproduction pipeline with the chosen primaries. Comparing the marked images with their gamut-mapped counter parts, it can be observed that the colors are reproduced faithfully. Furthermore, no significant artifacts are introduced in the laser-marked images. The gamut of the primary set allows marking diverse, vivid and relatively saturated colors as pointed by the image in the second row of Figure 12. Thanks to the high-resolution primaries, high spatial frequencies are pre served. This allows marking images with a vast level of details as shown in the bottom-right row of Figure 12. In the following, the question is studied of whether the marking parameters are transferable when using different marking settings or substrates. First a set of parameters reported in literature [Antohczak et al. 2013] are marked and show the results in Figure 13.

In Figure 13 the marking parameters from Antohczak et al. [2013], resulting in colors are shown in the middle row (as reported in the original paper). Same parameters marked on the same material (AISI 304) using the present device (bottom) lead to significant color differences (mean DEoo = 15.3) and a huge thickness variation. In the top the colors in the present gamut are shown closest to the reported colors in the middle. Despite using highly similar hardware and materials, the reported colors are not reproducible on the setup. Also, color thicknesses vary significantly making them unsuitable for halftoning and therefore image marking.

In the following table color differences are reported when marking on the present setup a general set of parameters in different circumstances:

The table shows the repeatability errors of color laser marking in form of DEoo mean (and standard deviation). The general set, consisting of 89 design points, is chosen to represent different colors in the explored gamut. A pure repeatability test is seen, on AISI 304 alloy, using the same marking settings leads to acceptable but not satisfactory accuracy. When the marking settings are changed by using a thicker substrate (2 mm), and therefore exiting the focal plane, the repeatability worsens. Finally, using a different alloy of stainless steel (AISI 430) results in the worst repeatability.

For comparison, in Table 1 , the result of the same experiments performed using the 6 extracted primaries are shown. Significantly higher accuracy in the pure repeatability experiment were observed as the primaries have been pruned against this circumstance. Interestingly, the primaries show acceptable repeatability when marked out of focus indicating that the extracted primaries are robust against some perturbations. However, the larger deviation when using a new substrate suggests that the primaries cannot be used for marking images on new substrates accurately. In Figure 14, an image on a new substrate (stainless steel AISI 430) is marked using primaries explored and extracted on the default substrate (AISI 304). While the image on the new substrate preserves the spatial details, there are signifi cant color shifts compared to the image marked on the default substrate. In order to show the present method is generalizable, a complete gamut exploration is performed, primary extrac tion and color reproduction on the new substrate. The results are shown in Figure 14 (bottom row). Figure 14 shows marked images on AISI 304 (top left) and AISI 430 (top right) with the same primaries explored and extracted on AISI 304 show significant color shifts. A newly ex plored and extracted set of primaries on AISI 430 shows good agreement between gamut- mapped (bottom left) and photograph (bottom right) of the same image marked on AISI 430. Furthermore, the full exploration on the new substrate, shown in Figure 15, leads to a differ ent gamut from the gamut obtained on the default substrate shown in Figure 8. Figure 15 shows a color gamut evolution of a full exploration (with fc, fas, fa, fasD, fosD) on AISI 430 stainless steel.

An iteration of the gamut exploration with a population size of 100 takes around 30 minutes. This includes marking the single clusters, measuring their thicknesses with a handheld mi croscope, marking the corresponding patches with proper distances of clusters, and finally capturing them with the colorimetric camera. The manual measurement of cluster thick nesses is the bottleneck as it takes approximately 20 minutes. Computing a new generation using the MCMOGA takes only a few seconds in Matlab. The marking time of an image is a function of the number of vectors (after halftone vectorization) and the marking speed of dif ferent primaries. A larger number of vectors causes more switching delays, making the mark ing time highly dependent on the image content in addition to its size. For example, the two marked images in the top row, and right side of the bottom row of Figure 12, despite a com parable image size (7 by 11 cm), required around 18 and 30 million vectors, and roughly 3 and 5.5 hours of marking time, respectively.

In this invention a computational framework is presented that enables a novel application of laser marking: color image reproduction. This method first characterizes the device using an evolutionary exploration of its performance space and then exploits that space for marking high-resolution, colorful images. A clear limitation of this method is the significant change of appearance from diffuse to non-diffuse configurations as shown in Figure 16, where a paint ing of Maria de’ Medici by Alessandro Allori, marked on AISI 304, and captured in non-diffuse (left) and diffuse (right) modes is depicted.

In Figure 17 laser marked images on stainless steel using method according to the present invention are shown (the plates are 13 c 13 cm). All the features disclosed in the application documents are claimed as being essential to the invention if, individually or in combination, they are novel over the prior art.

List of reference numerals

1 optical component

2 first gamut

3 step aa) of the method for preparing a laser marking system

4 step bb) of the method for preparing a laser marking system

5 step cc) of the method for preparing a laser marking system

6 step dd) of the method for preparing a laser marking system

7 step ee) of the method for preparing a laser marking system

8 step ff) of the method for preparing a laser marking system

9 iteration

10 design space 10a offspring design space 10b combined design space 11 design point 11a offspring design point 12 laser parameter

13 performance space 13a combined performance space

14 performance point

15 circular sectors

16 hue wheel 20 method for creating a colored laser mark on a specimen 21 step a) of the method for creating a colored laser mark on a specimen 22 step b) of the method for creating a colored laser mark on a specimen first al ternative

23 step b) of the method for creating a colored laser mark on a specimen second alternative

24 step c) of the method for creating a colored laser mark on a specimen

25 step d) of the method for creating a colored laser mark on a specimen

26 step e) of the method for creating a colored laser mark on a specimen 27 input image 28 color management workflow

28a juxtaposed halftoning workflow

29 step aa) of the color management workflow

30 step bb) of the color management workflow 31 step cc) of the color management workflow

32 step dd) of the color management workflow

33 step dd) of the color management workflow

34 step ee) of the color management workflow 100 laser marking system 101 laser 102 laser beam

103 scanning device

104 scanning device

105 specimen 105a surface layer

106 detection device

107 evaluation device

108 control unit 109 database

List of cited references

David Price Adams, Ryan D Murphy, David J Saiz, DA Hirschfeld, MA Rodriguez, Paul Ga briel Kotula, and BH Jared. 2014. Nanosecond pulsed laser irradiation of titanium: Oxide growth and effects on underlying metal. Surface and Coatings Technology 248 (2014), 38- 45.

Arkadiusz J Antohczak, Dariusz Kocoh, Maciej Nowak, Pawet Koziot, and Krzysztof M Abramski. 2013. Laser-induced colour marking - Sensitivity scaling for a stainless steel. In Applied Surface Science.

Arkadiusz J Antohczak, Bogusz St^pak, Pawet E Koziot, and Krzysztof M Abramski. 2014. The influence of process parameters on the laser-induced coloring of titanium. In Applied Physics A.

Thomas Auzinger, Wolfgang Heidrich, and Bernd Bickel. 2018. Computational design of nanostructural color for additive manufacturing. ACM Transactions on Graphics (TOG) 37, 4 (2018), 159.

Vahid Babaei and Roger D Hersch. 2012. Juxtaposed color halftoning relying on discrete lines. IEEE Transactions on Image Processing 22, 2 (2012), 679-686.

Vahid Babaei and Roger D Hersch. 2015. Yule-Nielsen based multi-angle reflectance predic tion of metallic halftones. In Color Imaging XX: Displaying, Processing, Hardcopy, and Appli cations, Vol. 9395. International Society for Optics and Photonics, 93950H.

Vahid Babaei and Roger D Hersch. 2016a. N-lnk Printer Characterization With Barycentric Subdivision. IEEE Transactions on Image Processing 25, 7 (2016), 3023-3031.

Vahid Babaei and Roger D Hersch. 2016b. Color reproduction of metallic-ink images. Journal of Imaging Science and Technology 60, 3 (2016), 30503-1.

Milton Birnbaum. 1965. Modulation of the reflectivity of semiconductors. Journal of Applied Physics 36, 2 (1965), 657-658.

Kalyanmoy Deb, Amrit Pratap, Sameer Agarwal, and TAMT Meyarivan. 2002. A fast and elit ist multiobjective genetic algorithm: NSGA-II. IEEE transactions on evolutionary computation 6, 2 (2002), 182-197.

A Perez Del Pino, JM Fernandez-Pradas, P Serra, and JL Morenza. 2004. Coloring of tita nium through laser oxidation: comparative study with anodizing. Surface and coatings tech nology 187, 1 (2004), 106-112.

Carlos M Fonseca, Peter J Fleming, et al. 1993. Genetic Algorithms for Multiobjective Optimi zation: Formulation, Discussion and Generalization. In lcga, Vol. 93. Citeseer, 416-423. Mathieu Hebert. 2014. Yule-Nielsen effect in halftone prints: graphical analysis method and improvement of the Yule-Nielsen transform. In Color Imaging XIX: Displaying, Processing, Hardcopy, and Applications, Vol. 9015. International Society for Optics and Photonics,

90150R.

Roger D Hersch, Philipp Donze, and Sylvain Chosson. 2007. Color images visible under UV light. In ACM Transactions on Graphics (TOG), Vol. 26. ACM, 75.

Guowei Hong, M Ronnier Luo, and Peter A Rhodes. 2001. A study of digital camera colori metric characterization based on polynomial modeling. Color Research & Application: En dorsed by Inter-Society Color Council, The Colour Group (Great Britain), Canadian Society for Color, Color Science Association of Japan, Dutch Society for the Study of Color, The Swedish Colour Centre Foundation, Colour Society of Australia, Centre Frangais de la Couleur 26, 1 (2001), 76-84.

P Laakso, S Ruotsalainen, H Pantsar, and R Penttila. 2009. Relation of laser parameters in color marking of stainless steel. In 12th Conference on Laser Processing of Materials in the Nordic Countries.

C Langlade, AB Vannes, JM Krafft, and JR Martin. 1998. Surface modification and tribologi cal behaviour of titanium and titanium alloys after YAG-laser treatments. Surface and Coat ings Technology 100 (1998), 383-387.

Samantha K Lawrence, David P Adams, David F Bahr, and Neville R Moody. 2013. Mechani cal and electromechanical behavior of oxide coatings grown on stainless steel 304L by nano second pulsed laser irradiation. Surface and Coatings Technology 235 (2013), 860-866.

KM teRcka, AJ Antonczak, B Szubzda, MR Wojcik, BD SteRpak, P Szymczyk, M Trzcinski,

M Ozimek, and KM Abramski. 2016. Effects of laser-induced oxidation on the corrosion re sistance of AISI 304 stainless steel. Journal of Laser Applications 28, 3 (2016), 032009.

Achim Lewandowski, Marcus Ludl, Gerald Byrne, and Georg Dorffner. 2006. Applying the Yule-Nielsen equation with negative n. JOSA A 23, 8 (2006), 1827-1834.

Huagang Liu, Wenxiong Lin, and Minghui Hong. 2019. Surface coloring by laser irradiation of solid substrates. APL Photonics 4, 5 (2019), 051101.

Laszlo Nanai, Robert Vajtai, and Thomas F George. 1997. Laser-induced oxidation of met als: state of the art. Thin Solid Films 298, 1-2 (1997), 160-164.

Petar Pjanic and Roger D Hersch. 2013. Specular color imaging on a metallic substrate. In Color and Imaging Conference, Vol. 2013. Society for Imaging Science and Technology, 61- 68.

Jean-Pierre Reveilles. 1995. Combinatorial pieces in digital lines and planes. In Vision geom etry IV, Vol. 2573. International Society for Optics and Photonics, 23-34. Robert Rolleston and Raja Balasubramanian. 1993. Accuracy of various types of Neugebauer model. In Color and Imaging Conference, Vol. 1993. Society for Imaging Sci ence and Technology, 32-37.

Janos Schanda. 2007. Colorimetry: Understanding the CIE System. John Wiley & Sons. Adriana Schulz, HarrisonWang, Eitan Crinspun, Justin Solomon, and Wojciech Matusik. 2018. Interactive exploration of design trade-offs. ACM Transactions on Graphics (TOG) 37, 4 (2018), 131.

Christian Schumacher, Bernd Bickel, Jan Rys, Steve Marschner, Chiara Daraio, and Markus Gross. 2015. Microstructures to control elasticity in 3D printing. ACM Transactions on Graphics (TOG) 34, 4 (2015), 136.

Gaurav Sharma, Wencheng Wu, and Edul N Dalai. 2005. The CIEDE2000 color-difference formula: Implementation notes, supplementary test data, and mathematical observations. Color research and application 30, 1 (2005), 21-30.

Pitchaya Sitthi-Amorn, Nicholas Modly, Westley Weimer, and Jason Lawrence. 2011. Ge netic programming for shader simplification. ACM Transactions on Graphics (TOG) 30, 6 (2011), 1-12.

Eric J Stollnitz, Victor Ostromoukhov, and David H Salesin. 1998. Reproducing color images using custom inks. In Proceedings of the 25th annual conference on Computer graphics and interactive techniques. ACM, 267-274.

Vadim Veiko, Galina Odintsova, Elena Gorbunova, Eduard Ageev, Alexandr Shimko, Yulia Karlagina, and Yaroslava Andreeva. 2016. Development of complete color palette based on spectrophotometric measurements of steel oxidation results for enhancement of color laser marking technology. Materials & Design 89 (2016), 684-688.

VP Veiko, AA Slobodov, and GV Odintsova. 2013. Availability of methods of chemical ther modynamics and kinetics for the analysis of chemical transformations on metal surfaces un der pulsed laser action. Laser Physics 23, 6 (2013), 066001.

Rui Wang, Bowen Yu, Julio Marco, Tianlei Hu, Diego Gutierrez, and Hujun Bao. 2016. Real time rendering on a power budget. ACM Transactions on Graphics (TOG) 35, 4 (2016), 1- 11.

Gunter Wyszecki and Walter Stanley Stiles. 1982. Color Science. Vol. 8. Wley New York.

JAC Yule and WJ Nielsen. 1951. The penetration of light into paper and its effect on halftone reproduction. In Proc. TAGA, Vol. 3. 65-76.

Bo Zhu, Melina Skouras, Desai Chen, and Wojciech Matusik. 2017. Two-scale topology opti mization with microstructures. ACM Transactions on Graphics (TOG) 36, 5 (2017), 164.