JP2017101971A - Element identification device, element identification program, and element identification method - Google Patents

Element identification device, element identification program, and element identification method Download PDF

Info

Publication number
JP2017101971A
JP2017101971A JP2015233915A JP2015233915A JP2017101971A JP 2017101971 A JP2017101971 A JP 2017101971A JP 2015233915 A JP2015233915 A JP 2015233915A JP 2015233915 A JP2015233915 A JP 2015233915A JP 2017101971 A JP2017101971 A JP 2017101971A
Authority
JP
Japan
Prior art keywords
groups
element identification
atomic
electron microscope
positions
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
JP2015233915A
Other languages
Japanese (ja)
Other versions
JP6551198B2 (en
Inventor
貴司 山▲崎▼
Takashi Yamazaki
貴司 山▲崎▼
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Fujitsu Ltd
Original Assignee
Fujitsu Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Fujitsu Ltd filed Critical Fujitsu Ltd
Priority to JP2015233915A priority Critical patent/JP6551198B2/en
Publication of JP2017101971A publication Critical patent/JP2017101971A/en
Application granted granted Critical
Publication of JP6551198B2 publication Critical patent/JP6551198B2/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Abstract

PROBLEM TO BE SOLVED: To provide an element identification device capable of identifying an element in a short time.SOLUTION: The element identification device includes: a calculation part that makes association an atom position in an image of scanning transmission electron microscope with a plurality of groups 52 and 54 based on a signal strength corresponding to the atom in the scanning transmission electron microscope image, calculates shortest routes 63 and 65 each of which passes all the atom positions in each of the plurality of groups, and replaces a part of the atom positions among the plurality of groups on the basis of the calculation result; and an element identification part that, after the calculation part has replaced a part of the atom position, identifies the atom position in the plurality of groups as the respective elements.SELECTED DRAWING: Figure 14

Description

本発明は、元素識別装置、元素識別プログラムおよび元素識別方法に関し、例えば、走査透過型電子顕微鏡像から元素を識別する元素識別装置、元素識別プログラムおよび元素識別方法に関する。   The present invention relates to an element identification device, an element identification program, and an element identification method. For example, the present invention relates to an element identification device, an element identification program, and an element identification method for identifying an element from a scanning transmission electron microscope image.

透過型電子顕微鏡(TEM:Transmission Electron Microscope)法は、例えば材料の原子レベルでの構造解析や組成解析を行なうときに用いられる。透過型電子顕微鏡法に、電子高エネルギー損失分光(EELS:Electron Energy- Loss Spectroscopy)法またはエネルギー分散X線分光(EDX:Energy Dispersive X-ray Spectroscopy)法を併用することで、原子レベルの元素識別が可能となっている(例えば、特許文献1および2)。   A transmission electron microscope (TEM) method is used, for example, when performing structural analysis or composition analysis at the atomic level of a material. Elemental element identification at the atomic level by using electron electron-loss spectroscopy (EELS) or energy dispersive X-ray spectroscopy (EDX) in combination with transmission electron microscopy (For example, Patent Documents 1 and 2).

特開2001−124709号公報JP 2001-124709 A 特開2001−153820号公報JP 2001-153820 A

TEM法とEELS法またはEDX法とを併用する方法では、原子ごとに分析を行なう。このため、画像撮影時間が長くなる。これにより、例えば試料のドリフトまたは電子線照射によるダメージが生じる。   In the method using both the TEM method and the EELS method or EDX method, analysis is performed for each atom. For this reason, the image capturing time becomes longer. As a result, for example, sample drift or damage due to electron beam irradiation occurs.

本元素識別装置、元素識別プログラムおよび元素識別方法は、短時間で元素識別を可能とすることを目的とする。   An object of the present element identification device, element identification program, and element identification method is to enable element identification in a short time.

走査透過型電子顕微鏡像における原子に対応した信号強度に基づき前記走査透過型電子顕微鏡像における原子位置をそれぞれ複数のグループに対応させ、前記複数のグループそれぞれ内において、原子位置を全て通る最短ルートを算出し、算出結果に基づき前記複数のグループ間で原子位置の一部を入れ換える演算部と、前記演算部が前記原子位置の一部を入れ換えた後、前記複数のグループ内の原子位置をそれぞれ異なる元素に識別する元素識別部と、を具備することを特徴とする元素識別装置を用いる。   Based on the signal intensity corresponding to the atoms in the scanning transmission electron microscope image, the atomic positions in the scanning transmission electron microscope image are made to correspond to a plurality of groups, respectively, and within each of the plurality of groups, the shortest route that passes through all the atom positions The calculation unit that calculates and replaces a part of the atomic positions between the plurality of groups based on the calculation result, and the atomic position in the plurality of groups is different after the calculation unit replaces a part of the atomic positions. An element identification device including an element identification unit for identifying an element is used.

コンピュータに、走査透過型電子顕微鏡像における原子に対応した信号強度に基づき前記走査透過型電子顕微鏡像における原子位置をそれぞれ複数のグループに対応させ、前記複数のグループそれぞれ内において、原子位置を全て通る最短ルートを算出し、算出結果に基づき前記複数のグループ間で原子位置の一部を入れ換え、前記原子位置の一部を入れ換えた後、前記複数のグループ内の原子位置をそれぞれ異なる元素に識別させることを特徴とする元素識別プログラムを用いる。   Based on the signal intensity corresponding to the atoms in the scanning transmission electron microscope image, the computer associates the atomic positions in the scanning transmission electron microscope image with a plurality of groups, and passes all the atomic positions in each of the plurality of groups. The shortest route is calculated, and based on the calculation result, a part of the atomic positions are replaced between the plurality of groups, and after a part of the atomic positions are replaced, the atomic positions in the plurality of groups are identified by different elements. An element identification program is used.

コンピュータに実行させる元素識別方法であって、走査透過型電子顕微鏡像における原子に対応した信号強度に基づき前記走査透過型電子顕微鏡像における原子位置をそれぞれ複数のグループに対応させ、前記複数のグループそれぞれ内において、原子位置を全て通る最短ルートを算出し、算出結果に基づき前記複数のグループ間で原子位置の一部を入れ換え、前記原子位置の一部を入れ換えた後、前記複数のグループ内の原子位置をそれぞれ異なる元素に識別させることを特徴とする元素識別方法を用いる。   An element identification method to be executed by a computer, the atomic positions in the scanning transmission electron microscope image corresponding to a plurality of groups based on the signal intensity corresponding to the atoms in the scanning transmission electron microscope image, and each of the plurality of groups And calculating a shortest route that passes through all the atom positions, replacing a part of the atom positions between the plurality of groups based on a calculation result, and replacing a part of the atom positions, and then replacing the atoms in the plurality of groups. An element identification method is used in which different elements are identified by different positions.

本元素識別装置、元素識別プログラムおよび元素識別方法によれば、短時間で元素識別を可能とすることができる。   According to this element identification device, element identification program, and element identification method, element identification can be performed in a short time.

図1は、実施例1における元素識別装置が用いられるシステムを示すブロック図である。FIG. 1 is a block diagram illustrating a system in which the element identification device according to the first embodiment is used. 図2は、実施例1において用いる電子顕微鏡の模式図である。FIG. 2 is a schematic diagram of an electron microscope used in the first embodiment. 図3は、実施例1に係る元素識別装置として機能するコンピュータのブロック図である。FIG. 3 is a block diagram of a computer that functions as an element identification device according to the first embodiment. 図4は、実施例1に係る元素識別装置の機能ブロック図である。FIG. 4 is a functional block diagram of the element identification device according to the first embodiment. 図5は、実施例1におけるコンピュータが行なう処理を示すフローチャートである。FIG. 5 is a flowchart illustrating processing performed by the computer according to the first embodiment. 図6は、実施例1における元素識別処理を示すフローチャートである。FIG. 6 is a flowchart showing the element identification process in the first embodiment. 図7は、電子顕微鏡像の例を示す図である。FIG. 7 is a diagram illustrating an example of an electron microscope image. 図8は、図7の信号強度のヒストグラムである。FIG. 8 is a signal intensity histogram of FIG. 図9は、図7の領域Aのピーク抽出像を示す図である。FIG. 9 is a diagram showing a peak extraction image of region A in FIG. 図10は、抽出したピーク強度のヒストグラムである。FIG. 10 is a histogram of the extracted peak intensity. 図11は、図10のヒストグラムにおいて識別困難な範囲を示す図である。FIG. 11 is a diagram showing a range that is difficult to identify in the histogram of FIG. 図12(a)から図12(c)は、図7の電子顕微鏡像に対応する原子位置を示す図である。FIG. 12A to FIG. 12C are diagrams showing atomic positions corresponding to the electron microscope image of FIG. 図13(a)から図13(d)は、各グループに含まれる原子位置を示す図(その1)である。FIG. 13A to FIG. 13D are diagrams (part 1) showing atomic positions included in each group. 図14(a)から図14(c)は、各グループに含まれる原子位置を示す図(その2)である。FIG. 14A to FIG. 14C are diagrams (part 2) showing the positions of atoms included in each group. 図15(a)および図15(b)は、元素識別処理における電子顕微鏡像に対応する原子位置を示す図(その1)である。FIG. 15A and FIG. 15B are diagrams (part 1) illustrating atomic positions corresponding to an electron microscope image in element identification processing. 図16(a)および図16(b)は、元素識別処理における電子顕微鏡像に対応する原子位置を示す図(その2)である。FIG. 16A and FIG. 16B are diagrams (part 2) illustrating atomic positions corresponding to an electron microscope image in element identification processing.

図面を参照し、実施例について説明する。   Embodiments will be described with reference to the drawings.

図1は、実施例1における元素識別装置が用いられるシステムを示すブロック図である。図1に示すように、システムは、元素識別装置30、電子顕微鏡32、像入力装置34および画像出力装置36を備えている。元素識別装置30は、例えばコンピュータである。電子顕微鏡32は、走査透過型電子顕微鏡である。像入力装置34は、例えばコンピュータ等である。画像出力装置36は、例えば液晶表示装置等の表示装置、プリンタ等の印刷装置、またはこれらを有するコンピュータである。   FIG. 1 is a block diagram illustrating a system in which the element identification device according to the first embodiment is used. As shown in FIG. 1, the system includes an element identification device 30, an electron microscope 32, an image input device 34, and an image output device 36. The element identification device 30 is, for example, a computer. The electron microscope 32 is a scanning transmission electron microscope. The image input device 34 is, for example, a computer. The image output device 36 is, for example, a display device such as a liquid crystal display device, a printing device such as a printer, or a computer having these.

電子顕微鏡32は、原子レベルの電子顕微鏡像を撮影し、元素識別装置30に電子顕微鏡像を出力する。像入力装置34は、電子顕微鏡32が過去に撮影した電子顕微鏡像を元素識別装置30に出力する。電子顕微鏡32および像入力装置34の少なくとも一方が設けられていればよい。元素識別装置30は、電子顕微鏡像内の元素を識別し、元素識別画像を画像出力装置36に出力する。画像出力装置36は、元素識別画像を表示する。元素識別装置30と、電子顕微鏡32、像入力装置34および画像出力装置36とは、例えばインターネット等で接続されていてもよい。元素識別装置30と、電子顕微鏡32、像入力装置34および画像出力装置36と、の画像等の入出力は可搬型記憶媒体を用いて行なってもよい。   The electron microscope 32 captures an atomic level electron microscope image and outputs the electron microscope image to the element identification device 30. The image input device 34 outputs an electron microscope image taken by the electron microscope 32 in the past to the element identification device 30. It is sufficient that at least one of the electron microscope 32 and the image input device 34 is provided. The element identification device 30 identifies an element in the electron microscope image and outputs an element identification image to the image output device 36. The image output device 36 displays an element identification image. The element identification device 30, the electron microscope 32, the image input device 34, and the image output device 36 may be connected, for example, via the Internet. Input / output of images and the like between the element identification device 30, the electron microscope 32, the image input device 34, and the image output device 36 may be performed using a portable storage medium.

図2は、実施例1において用いる電子顕微鏡の模式図である。図2に示すように、電子顕微鏡32内に試料40が設けられている。試料40の上面に電子線41が照射される。電子線41は、試料40の上面において原子レベルまで収束されている。電子線41は、試料40の上面を走査される。試料40を透過した電子は、高角度に散乱される電子43と低角度に散乱される電子45とに分類される。高角度に散乱される電子43は、試料40内で熱散漫散乱された電子である。低角度に散乱される電子45は、試料40内で電子に弾性散乱された電子である。環状の暗視野検出器42は、電子43を検出する。円状の明視野検出器44は電子45を検出する。   FIG. 2 is a schematic diagram of an electron microscope used in the first embodiment. As shown in FIG. 2, a sample 40 is provided in the electron microscope 32. An electron beam 41 is irradiated on the upper surface of the sample 40. The electron beam 41 is converged to the atomic level on the upper surface of the sample 40. The electron beam 41 is scanned on the upper surface of the sample 40. The electrons transmitted through the sample 40 are classified into electrons 43 scattered at a high angle and electrons 45 scattered at a low angle. The electrons 43 scattered at a high angle are electrons diffused and diffused in the sample 40. The electrons 45 scattered at a low angle are electrons that are elastically scattered by the electrons in the sample 40. The annular dark field detector 42 detects the electrons 43. A circular bright field detector 44 detects electrons 45.

走査された電子線41の位置に応じ暗視野検出器42および明視野検出器44が出力する信号強度を像として形成する。暗視野検出器42が出力した信号強度の像が暗視野像であり、HAADF(High-angle Annular Dark Field) STEM像と呼ばれる。明視野検出器44が出力した信号強度の像が明視野像であり、HABF(High-angle Bright Field) STEM像と呼ばれる。HAADF STEM像の明暗を反転するとほぼHABF STEM像となる。HAADF STEM像における原子位置の信号強度は主に原子番号に由来する。すなわち、信号強度は原子番号が大きいほど大きくなる。実施例1では、HAADF STEM像を用いる。   The signal intensity output from the dark field detector 42 and the bright field detector 44 is formed as an image according to the position of the scanned electron beam 41. The image of the signal intensity output from the dark field detector 42 is a dark field image and is called a HAADF (High-angle Annular Dark Field) STEM image. An image of the signal intensity output from the bright field detector 44 is a bright field image and is called a HABF (High-angle Bright Field) STEM image. When the brightness of the HAADF STEM image is reversed, a HABF STEM image is obtained. The signal intensity at the atomic position in the HAADF STEM image is mainly derived from the atomic number. That is, the signal intensity increases as the atomic number increases. In Example 1, a HAADF STEM image is used.

図3は、実施例1に係る元素識別装置として機能するコンピュータのブロック図である。コンピュータ10は、CPU(Central Processing Unit)11、表示装置12、入力装置13、出力装置14、記憶装置15、記憶媒体用ドライブ16、通信インターフェース17および内部バス18を備えている。表示装置12は、例えば液晶パネル等の表示パネルを含み、コマンドまたはデータ等を表示する。入力装置13は、例えばキーボード、マウスおよびタッチパネル等であり、コマンドまたはデータ等を入力する。出力装置14は、例えばプリンタであり、コマンドまたはデータ等を出力する。記憶装置15は、例えばRAM(Random Access Memory)等の揮発性メモリ、またはフラッシュメモリもしくはハードディスク等の不揮発性メモリであり、プログラム、処理中または処理後のデータを記憶する。記憶媒体用ドライブ16は、記憶媒体19に格納されたプログラムをインストールする際に用いる。また、記憶媒体用ドライブ16は、電子顕微鏡像を記憶媒体19から取得する、または元素識別画像を記憶媒体19に記憶させる。通信インターフェース17は、電子顕微鏡32または像入力装置34から電子顕微鏡像データを取得する。また、通信インターフェース17は、画像出力装置36に元素識別画像を出力する。内部バス18は、コンピュータ10内の各装置を接続する。   FIG. 3 is a block diagram of a computer that functions as an element identification device according to the first embodiment. The computer 10 includes a central processing unit (CPU) 11, a display device 12, an input device 13, an output device 14, a storage device 15, a storage medium drive 16, a communication interface 17, and an internal bus 18. The display device 12 includes a display panel such as a liquid crystal panel, for example, and displays commands or data. The input device 13 is a keyboard, a mouse, a touch panel, or the like, for example, and inputs commands or data. The output device 14 is a printer, for example, and outputs a command or data. The storage device 15 is, for example, a volatile memory such as a RAM (Random Access Memory), or a nonvolatile memory such as a flash memory or a hard disk, and stores a program, data during or after processing. The storage medium drive 16 is used when a program stored in the storage medium 19 is installed. The storage medium drive 16 acquires an electron microscope image from the storage medium 19 or stores an element identification image in the storage medium 19. The communication interface 17 acquires electron microscope image data from the electron microscope 32 or the image input device 34. Further, the communication interface 17 outputs an element identification image to the image output device 36. The internal bus 18 connects each device in the computer 10.

プログラムを格納するコンピュータ10が読み取り可能な記憶媒体19として可搬型記憶媒体を用いることができる。可搬型記憶媒体としては、例えば、CD−ROM(Compact Disc Read Only Memory)ディスク、DVD(Digital Video Disc)ディスク、ブルーレイディスクまたはUSB(Universal Serial Bus)メモリ等を用いることができる。記憶媒体19として、フラッシュメモリまたはHDD(Hard Disk Drive)等を用いてもよい。   A portable storage medium can be used as the storage medium 19 that can be read by the computer 10 storing the program. As the portable storage medium, for example, a CD-ROM (Compact Disc Read Only Memory) disc, a DVD (Digital Video Disc) disc, a Blu-ray disc, or a USB (Universal Serial Bus) memory can be used. As the storage medium 19, a flash memory or an HDD (Hard Disk Drive) may be used.

図4は、実施例1に係る元素識別装置の機能ブロック図である。コンピュータ10は、図3に図示した各部材とソフトウエアとの協働により図4のように機能する。像取得部20は、電子顕微鏡32または像入力装置34から電子顕微鏡像を取得する。演算部22は、電子顕微鏡像における原子に対応した信号強度に基づき、原子位置のグループ分けを行なう。詳細は後述する。元素識別部24は、分けられた複数のグループ内の原子位置をそれぞれ異なる元素に識別する。詳細は後述する。画像生成部26は、元素の識別結果に基づき元素識別画像を生成し、出力する。   FIG. 4 is a functional block diagram of the element identification device according to the first embodiment. The computer 10 functions as shown in FIG. 4 in cooperation with each member shown in FIG. 3 and software. The image acquisition unit 20 acquires an electron microscope image from the electron microscope 32 or the image input device 34. The computing unit 22 groups atomic positions based on the signal intensity corresponding to the atoms in the electron microscope image. Details will be described later. The element identifying unit 24 identifies atomic positions in the divided groups as different elements. Details will be described later. The image generation unit 26 generates and outputs an element identification image based on the element identification result.

図5は、実施例1におけるコンピュータが行なう処理を示すフローチャートである。図6は、実施例1における元素識別処理を示すフローチャートである。図5に示すように、像取得部20は、電子顕微鏡32または像入力装置34から電子顕微鏡像を取得する(ステップS10)。例えば、CPU11は、通信インターフェース17を介し電子顕微鏡像を取得する。演算部22は、電子顕微鏡像からピーク抽出像を演算する(ステップS12)。演算部22は、ピーク強度のヒストグラムを演算する(ステップS14)。   FIG. 5 is a flowchart illustrating processing performed by the computer according to the first embodiment. FIG. 6 is a flowchart showing the element identification process in the first embodiment. As shown in FIG. 5, the image acquisition unit 20 acquires an electron microscope image from the electron microscope 32 or the image input device 34 (step S10). For example, the CPU 11 acquires an electron microscope image via the communication interface 17. The computing unit 22 computes a peak extraction image from the electron microscope image (step S12). The calculation unit 22 calculates a histogram of peak intensity (step S14).

ステップS10からS14の処理について、図7から図10を用い説明する。図7は、電子顕微鏡像の例を示す図であり、SrTiO上にLaCoOを成長させた試料の断面の暗視野画像(すなわちHAADF STEM像)である。図7に示すように、ドメイン56と58の境界55を挟んで、左側がSrTiO(チタン酸ストロンチウム)のドメイン56であり、右側がLaCoO(コバルト酸タンタル)のドメイン58である。原子像の明るさ(信号強度)は、元素の原子番号に依存する。 The processing in steps S10 to S14 will be described with reference to FIGS. FIG. 7 is a diagram showing an example of an electron microscope image, which is a dark field image (ie, HAADF STEM image) of a cross section of a sample obtained by growing LaCoO 3 on SrTiO 3 . As shown in FIG. 7, the left side is the domain 56 of SrTiO 3 (strontium titanate) and the right side is the domain 58 of LaCoO 3 (tantalum cobaltate) across the boundary 55 between the domains 56 and 58. The brightness (signal intensity) of the atomic image depends on the atomic number of the element.

図8は、図7の信号強度のヒストグラムであり、信号強度(任意座標)に対する電子顕微鏡像内のピクセル(画素)数を示す。図8に示すように、ヒストグラムには2つのピークが存在する。図8の中には、Ti(チタン)、Sr(ストロンチウム)、Co(コバルト)およびLa(ランタン)の4つの元素の信号が含まれる。しかし、図8のヒストグラムから各元素に分離することは難しい。これは、図7には、原子の中心位置以外にも信号が存在するためである。   FIG. 8 is a histogram of the signal intensity of FIG. 7 and shows the number of pixels (pixels) in the electron microscope image with respect to the signal intensity (arbitrary coordinates). As shown in FIG. 8, there are two peaks in the histogram. In FIG. 8, signals of four elements of Ti (titanium), Sr (strontium), Co (cobalt), and La (lanthanum) are included. However, it is difficult to separate each element from the histogram of FIG. This is because there are signals in FIG. 7 other than the center position of the atoms.

そこで、ステップS14において、図7の電子顕微鏡像から各原子の信号強度がピークとなるピクセルを抽出してピーク抽出画像とする。図9は、図7の領域Aのピーク抽出像を示す図である。図9に示すように、各原子(TiおよびSr)におけるピークの信号のピクセルのみが白くなっている。すなわち、各原子(TiおよびSr)におけるピークの信号のピクセルのみが抽出されている。電子顕微鏡像内における抽出されたピクセルの位置を原子位置とし、抽出されたピクセルの信号強度を信号強度がピークとなる強度(すなわちピーク強度)とする。実施例1では、各原子において、最も信号強度の大きいピクセルを抽出している。例えば、各原子において信号強度をガウス分布となると仮定してフィッティングすることで、原子位置とピーク強度を演算してもよい。   Therefore, in step S14, a pixel having a peak signal intensity of each atom is extracted from the electron microscope image of FIG. FIG. 9 is a diagram showing a peak extraction image of region A in FIG. As shown in FIG. 9, only the pixel of the peak signal in each atom (Ti and Sr) is white. That is, only the pixel of the peak signal in each atom (Ti and Sr) is extracted. The position of the extracted pixel in the electron microscope image is an atomic position, and the signal intensity of the extracted pixel is the intensity at which the signal intensity reaches a peak (that is, peak intensity). In the first embodiment, the pixel having the highest signal intensity is extracted for each atom. For example, the atomic position and peak intensity may be calculated by fitting assuming that the signal intensity is Gaussian in each atom.

図10は、抽出したピーク強度のヒストグラムであり、ピーク強度(任意座標)に対するピーク抽出像のピクセル数を示す。図10に示すように、ヒストグラムには3つのピークが存在する。ピーク強度の小さいピークはTiに対応し、ピーク強度の大きいピークはLaに相当する。しかし、SrとCoのピーク強度のピークは重なっており分離できていない。   FIG. 10 is a histogram of the extracted peak intensity, and shows the number of pixels of the peak extracted image with respect to the peak intensity (arbitrary coordinates). As shown in FIG. 10, there are three peaks in the histogram. A peak with a low peak intensity corresponds to Ti, and a peak with a high peak intensity corresponds to La. However, the peak intensity peaks of Sr and Co overlap and cannot be separated.

図5に戻り、元素識別部24は、ピーク強度のヒストグラムに基づき、各原子の元素を識別する(ステップS16)。例えば、図10において、3つのピークを各元素と識別する。例えば、最もピーク強度が小さいピークを、原子番号の最も小さいTiと識別する。最もピーク強度が大きいピークを、最も原子番号の大きなLaと識別する。これにより、図10のグループ51に含まれる原子位置をTi、グループ53に含まれる原子位置をLaと識別する。   Returning to FIG. 5, the element identifying unit 24 identifies the element of each atom based on the peak intensity histogram (step S <b> 16). For example, in FIG. 10, three peaks are identified from each element. For example, the peak having the smallest peak intensity is identified as Ti having the smallest atomic number. The peak having the highest peak intensity is identified as La having the largest atomic number. Thereby, the atomic position included in the group 51 of FIG. 10 is identified as Ti, and the atomic position included in the group 53 is identified as La.

図5に戻り、演算部22は、識別困難なピークがあるか判定する(ステップS18)。例えば、図10のように、演算部22は、構成する元素数とピーク強度のヒストグラムにおけるピーク数が異なる場合、識別困難と判定する。NoのときステップS24に進む。Yesのとき、演算部22は、識別困難な範囲を設定する(ステップS20)。図11は、図10のヒストグラムにおいて識別困難な範囲を示す図である。図11に示すように、ヒストグラムは図10と同じである。範囲50は、SrとCoのピーク強度が重なっており、ヒストグラムから元素の識別が困難な範囲である。そこで、演算部22は、範囲50を識別困難範囲と設定する。   Returning to FIG. 5, the calculation unit 22 determines whether there is a peak that is difficult to identify (step S <b> 18). For example, as illustrated in FIG. 10, when the number of constituent elements and the number of peaks in the peak intensity histogram are different, the calculation unit 22 determines that identification is difficult. When No, the process proceeds to step S24. When the result is Yes, the calculation unit 22 sets a range that is difficult to identify (step S20). FIG. 11 is a diagram showing a range that is difficult to identify in the histogram of FIG. As shown in FIG. 11, the histogram is the same as FIG. The range 50 is a range in which the peak intensities of Sr and Co overlap and it is difficult to identify elements from the histogram. Therefore, the calculation unit 22 sets the range 50 as a difficult-to-identify range.

図5に戻り、演算部22は、元素識別処理を行なう(ステップS22)。図6に示すように、演算部22は、識別困難範囲を複数のグループに分割する(ステップS30)。図11のように、範囲50を互いに重ならず、かつ隣接する2つのグループ52と54に分割する。   Returning to FIG. 5, the calculation unit 22 performs element identification processing (step S <b> 22). As shown in FIG. 6, the calculating part 22 divides | segments the identification difficulty range into a some group (step S30). As shown in FIG. 11, the range 50 is divided into two adjacent groups 52 and 54 that do not overlap each other.

図12(a)から図12(c)は、図7の電子顕微鏡像に対応する原子位置を示す図であり、それぞれ、図11に示す範囲50、グループ52および54に含まれるピーク強度を有するピクセルの原子位置を示す図である。図12(a)に示すように、範囲50に含まれるピーク強度を有するピクセルの原子位置は全体に分布しており、SrとCoが識別されていない。図12(b)に示すように、グループ52に含まれるピーク強度の小さい原子位置は主にドメイン56に位置するが、ドメイン58に位置するものもある。図12(c)に示すように、グループ54に含まれるピーク強度の大きい原子位置はドメイン58に位置するが、ドメイン58に原子位置が表示されていない領域もある。このように、範囲50をグループ52と54に分割しただけでは、元素の識別ができていないことがわかる。   12 (a) to 12 (c) are diagrams showing atomic positions corresponding to the electron microscope image of FIG. 7, and have peak intensities included in the range 50 and the groups 52 and 54 shown in FIG. 11, respectively. It is a figure which shows the atomic position of a pixel. As shown in FIG. 12A, the atomic positions of the pixels having the peak intensity included in the range 50 are distributed throughout, and Sr and Co are not identified. As shown in FIG. 12 (b), atomic positions with small peak intensities included in the group 52 are mainly located in the domain 56, but some are located in the domain 58. As shown in FIG. 12 (c), the atomic position having a high peak intensity included in the group 54 is located in the domain 58, but there is also a region where the atomic position is not displayed in the domain 58. Thus, it can be seen that the elements cannot be identified by simply dividing the range 50 into the groups 52 and 54.

図6に戻り、演算部22は、各グループ52、54において、原子位置をすべて通る一筆書きの最短距離のルートを演算する(ステップS32)。この問題は巡回セールスマン問題と呼ばれ、さまざまな解法が知られている。実施例1ではやきなまし(シュミレーテッドアニーリング)法を用いる。演算部22は、ルート上において隣接する原子位置の間隔が最近接原子間距離より長い箇所が存在するか判定する(ステップS34)。Noのとき終了する。Yesのとき、演算部22は、最近接原子間距離より長い原子位置を他グループに移動させる(ステップS36)。   Returning to FIG. 6, the computing unit 22 computes the shortest distance route of one-stroke writing passing through all the atom positions in each of the groups 52 and 54 (step S <b> 32). This problem is called the traveling salesman problem, and various solutions are known. In the first embodiment, a annealed (simulated annealing) method is used. The computing unit 22 determines whether there is a place where the interval between adjacent atomic positions on the route is longer than the nearest interatomic distance (step S34). If no, the process ends. When the result is Yes, the calculation unit 22 moves an atomic position longer than the closest interatomic distance to another group (step S36).

演算部22は、ルート上において隣接する原子位置間距離が一定の距離より長い経路で挟まれた原子位置がないか判定する(ステップS38)。一定の距離は、例えば最近接原子間距離の2倍とする。NoのときステップS42に進む。Yesのとき、演算部22は、この経路に挟まれた原子位置をすべて他のグループに移動する(ステップS40)。演算部22は、2つ前の原子位置の配置と同じか判定する(ステップS42)。Yesのとき終了する。NoのときステップS32に戻る。   The computing unit 22 determines whether there is an atomic position sandwiched by a route having a distance between adjacent atomic positions on the route that is longer than a certain distance (step S38). The fixed distance is, for example, twice the closest interatomic distance. When No, the process proceeds to step S42. When the result is Yes, the calculation unit 22 moves all the atomic positions sandwiched by this path to another group (step S40). The calculation unit 22 determines whether the arrangement is the same as the previous atomic position (step S42). When yes, the process ends. If no, the process returns to step S32.

まず、ステップS38およびS40を除くフローについて、説明する。図13(a)から図14(c)は、各グループに含まれる原子位置を示す図である。図13(a)は、例えば図11の範囲50に含まれる原子位置を模式的に示している。図13(b)から図14(c)は、グループ52および54に含まれる原子位置を模式的に示している。ドメイン56および58は未知の場合もあるが、元素62はドメイン56に含まれ、元素64はドメイン58に含まれているとする。境界55はドメイン56と58の境界である。   First, the flow excluding steps S38 and S40 will be described. FIG. 13A to FIG. 14C are diagrams showing atomic positions included in each group. FIG. 13A schematically shows, for example, atomic positions included in the range 50 of FIG. FIGS. 13B to 14C schematically show the positions of atoms included in the groups 52 and 54. Although the domains 56 and 58 may be unknown, the element 62 is included in the domain 56 and the element 64 is included in the domain 58. The boundary 55 is a boundary between the domains 56 and 58.

図13(a)に示すように、ドメイン56に元素62が6個配列され、ドメイン58に元素64が6個配列されている。範囲50に含まれる原子位置では、元素62および64が分類されていないことがわかる。図13(b)に示すように、ステップS30において範囲50を2つのグループ52および54に分割する。原子位置64aはグループ52に含まれており、ドメイン56および58と完全には分類されていない。図13(c)に示すように、ステップS32において各グループ52および54において巡回セールスマン問題を解き、原子位置をすべて通る一筆書きで最短距離のルート63および65を演算する。   As shown in FIG. 13A, six elements 62 are arranged in the domain 56, and six elements 64 are arranged in the domain 58. It can be seen that at the atomic positions included in the range 50, the elements 62 and 64 are not classified. As shown in FIG. 13B, the range 50 is divided into two groups 52 and 54 in step S30. Atomic position 64 a is included in group 52 and is not completely classified as domains 56 and 58. As shown in FIG. 13C, in step S32, the traveling salesman problem is solved in each of the groups 52 and 54, and the shortest distance routes 63 and 65 are calculated by a single stroke passing through all the atomic positions.

図13(d)に示すように、ステップS34においてルート63および65上に隣接する原子位置間距離が最近接原子間距離より長い経路63aおよび65aが存在する。経路63aおよび65aは、同じドメイン56および58内の隣接する原子位置間を接続する経路としては不自然である。この場合、経路63aおよび65aが接続する原子位置64aおよび64b(丸で囲んだ原子位置)は本来のグループに属していない可能性がある。原子位置64aのように、1つの原子位置が経路63aに挟まれているとき、挟まれた原子位置64aを抽出する。経路65aのように、経路65aが単独で存在するとき、ルート65の走査方向により経路65aの両端のいずれか一方を原子位置64bとして抽出する。   As shown in FIG. 13D, there are paths 63a and 65a in which the distance between adjacent atomic positions on the routes 63 and 65 is longer than the closest interatomic distance in step S34. The paths 63a and 65a are unnatural as paths connecting adjacent atomic positions in the same domains 56 and 58. In this case, the atom positions 64a and 64b (circled atom positions) to which the paths 63a and 65a are connected may not belong to the original group. When one atomic position is sandwiched between the paths 63a as in the atomic position 64a, the sandwiched atomic position 64a is extracted. When the path 65a exists alone as in the path 65a, either one of both ends of the path 65a is extracted as the atomic position 64b depending on the scanning direction of the route 65.

図14(a)に示すように、ステップS36において、原子位置64aをグループ54に移動し、原子位置64bをグループ52に移動する。原子位置64bは、グループ52に移動しなくてもよい原子位置である。しかし、図13(d)において、グループ54に原子位置64aが抜けていたため、経路65aが抽出され、結果的に原子位置64bがグループ52に移動している。図14(b)に示すように、ステップS32に戻り、原子位置をすべて通る一筆書きで最短距離のルート63および65を演算する。ステップS34においてルート63上に隣接する原子位置間距離が最近接原子間距離より長い経路63aが存在する。グループ54では長い経路は存在しない。図14(c)において、長い経路63aに接続する原子位置64bをグループ54に移動する。これにより、ドメイン56および58にそれぞれ元素62および64を分類できる。ステップS42において、この後、同じ演算を繰り返しても同じ結果となるため、演算を終了する。このように、ステップS30からS36を繰り返すことにより、元素の識別が可能となる。   As shown in FIG. 14A, the atomic position 64a is moved to the group 54 and the atomic position 64b is moved to the group 52 in step S36. The atomic position 64 b is an atomic position that does not have to move to the group 52. However, in FIG. 13D, the atomic position 64 a is missing from the group 54, so the path 65 a is extracted, and as a result, the atomic position 64 b has moved to the group 52. As shown in FIG. 14B, the process returns to step S32, and the shortest distance routes 63 and 65 are calculated by a single stroke passing through all the atom positions. In step S <b> 34, a path 63 a in which the distance between adjacent atomic positions on the route 63 is longer than the closest interatomic distance exists. In group 54, there is no long path. In FIG. 14C, the atomic position 64b connected to the long path 63a is moved to the group 54. Thereby, the elements 62 and 64 can be classified into the domains 56 and 58, respectively. In step S42, since the same result is obtained even if the same calculation is repeated thereafter, the calculation is terminated. In this way, the elements can be identified by repeating steps S30 to S36.

図12(b)および図12(c)に対し、ステップS32からS42を繰り返した演算例を説明する。図15(a)から図16(b)は、元素識別処理における電子顕微鏡像に対応する原子位置を示す図である。各図において左側がグループ52、右側がグループ54に対応する。各グループ52および54内の左側がドメイン56に相当し、右側がドメイン58に相当する。実際の演算ではドメイン56および58は未知であることもあるが、わかりやすくするためドメイン56および58を用いて説明する。   An example of calculation in which steps S32 to S42 are repeated with respect to FIGS. 12B and 12C will be described. FIG. 15A to FIG. 16B are diagrams showing atomic positions corresponding to an electron microscope image in element identification processing. In each figure, the left side corresponds to the group 52 and the right side corresponds to the group 54. The left side in each group 52 and 54 corresponds to the domain 56, and the right side corresponds to the domain 58. In an actual calculation, the domains 56 and 58 may be unknown, but the domains 56 and 58 will be used for easy understanding.

図15(a)は、図12(a)および図12(b)の原子位置について巡回セールスマン問題を解いた結果である。グループ52および54についてルート63および65が設定されている。グループ52では、ドメイン58内にルート63上の原子位置が多く残っている。グループ54では、ドメイン56内にルート65上の原子位置が残っている。図15(b)は、ステップS32からS42を10回繰り返した結果である。グループ54では、ドメイン56内にルート65上の原子位置は境界55付近以外存在しない。しかし、グループ52ではドメイン58内にルート63上の原子位置が多数残っている。   FIG. 15A shows the result of solving the traveling salesman problem with respect to the atomic positions in FIGS. 12A and 12B. Routes 63 and 65 are set for the groups 52 and 54. In the group 52, many atomic positions on the route 63 remain in the domain 58. In the group 54, the atomic position on the route 65 remains in the domain 56. FIG. 15B shows the result of repeating steps S32 to S42 10 times. In the group 54, there is no atomic position on the route 65 in the domain 56 except near the boundary 55. However, in the group 52, many atomic positions on the route 63 remain in the domain 58.

図16(a)は、ステップS32からS42を50回繰り返した結果である。グループ52では、ドメイン58内にルート63上の原子位置のブロック67が存在する。ブロック67のルート63上の両端は経路63bを介しドメイン56に接続されている。経路63bは、ルート63上の隣接する原子位置間距離が最近接原子間距離の2倍より長い。このような経路63bに挟まれたブロック67は、ドメイン56と同じドメインに属するには不自然である。ステップS38において演算部22はYesと判定する。ステップS40において、演算部22は、経路63bで挟まれたブロック67をグループ54に移動する。   FIG. 16A shows the result of repeating steps S32 to S42 50 times. In the group 52, a block 67 of an atomic position on the route 63 exists in the domain 58. Both ends of the block 67 on the route 63 are connected to the domain 56 via a path 63b. In the path 63b, the distance between adjacent atomic positions on the route 63 is longer than twice the distance between the nearest atoms. Such a block 67 sandwiched between the paths 63 b is unnatural to belong to the same domain as the domain 56. In step S38, the calculation unit 22 determines Yes. In step S40, the calculation unit 22 moves the block 67 sandwiched by the path 63b to the group 54.

図16(b)に示すように、ブロック67をグル−プ54に移動すると、グループ52および54はドメイン56および58の原子位置に分類される。この後、ステップS32からS42を繰り返しても結果は変わらず収束する。よって、ステップS42において演算部22はYesと判定する。収束した後も2つの配置を繰り返すことがある。このため、ステップS42では、ステップS32からS42のループにおける2回前の原子位置の配置から変化していないかを判定することにより、収束したか否かを判定している。   As shown in FIG. 16 (b), when the block 67 is moved to the group 54, the groups 52 and 54 are classified into the atomic positions of the domains 56 and 58. After this, even if steps S32 to S42 are repeated, the result does not change and converges. Therefore, in step S42, the calculation unit 22 determines Yes. Even after convergence, the two arrangements may be repeated. For this reason, in step S42, it is determined whether or not it has converged by determining whether or not it has changed from the arrangement of atomic positions two times before in the loop of steps S32 to S42.

図6に戻り、ステップS42において演算部22がYesと判定したとき、終了し、図5のステップS24に進む。元素識別部24は、グループ52および54ごとに元素を識別し、元素識別画像を生成する。その後、元素識別画像を画像出力装置36に出力する(ステップS24)。例えば、グループ52内の原子位置をSrと識別し、グループ54内の原子位置をCoと判定する。ステップS16において識別したTiおよびLaの原子位置を合わせ、電子顕微鏡に対応する元素識別画像を生成する。例えば、元素識別画像は電子顕微鏡像の原子像に元素の種類に応じ色をつけた画像である。また、元素識別画像は電子顕微鏡像内の特定の元素のみの原子像を表示した画像でもよい。例えば元素識別画像として、それぞれTi、Sr、CaおよびLaに対応する原子像を表示した4つの画像を生成してもよい。また、元素識別部24は、元素識別画像を生成せず、元素識別結果を数値データとして生成してもよい。その後終了する。   Returning to FIG. 6, when the calculation unit 22 determines Yes in step S <b> 42, the process ends, and the process proceeds to step S <b> 24 in FIG. 5. The element identification unit 24 identifies the element for each of the groups 52 and 54 and generates an element identification image. Thereafter, the element identification image is output to the image output device 36 (step S24). For example, the atomic position in the group 52 is identified as Sr, and the atomic position in the group 54 is determined as Co. The atomic positions of Ti and La identified in step S16 are aligned, and an element identification image corresponding to the electron microscope is generated. For example, the element identification image is an image in which an atomic image of an electron microscope image is colored according to the type of element. The element identification image may be an image displaying an atomic image of only a specific element in the electron microscope image. For example, four images displaying atomic images respectively corresponding to Ti, Sr, Ca, and La may be generated as element identification images. Moreover, the element identification part 24 may produce | generate an element identification result as numerical data, without producing | generating an element identification image. Then exit.

例えば、実施例1の元素識別プログラムおよび元素識別方法を市販のパーソナルコンピュータで実行した場合、1辺が512ピクセルの電子顕微鏡像から原子位置が1000個以上の元素を数秒程度で識別できる。一方、EELS法を用いて同程度の分解能で元素を識別する場合、1ピクセル当たりのEELS分析時間が約0.1秒とすると、1辺が512ピクセルの元素識別の時間は約7時間となる。このように、実施例1では、元素識別を極めて高速で行なうことが可能となる。   For example, when the element identification program and the element identification method of Example 1 are executed by a commercially available personal computer, an element having an atomic position of 1000 or more can be identified in about several seconds from an electron microscope image having one side of 512 pixels. On the other hand, when elements are identified with the same resolution using the EELS method, if the EELS analysis time per pixel is about 0.1 seconds, the element identification time for 512 pixels per side is about 7 hours. . Thus, in Example 1, element identification can be performed at an extremely high speed.

以上、図7のように元素の種類およびドメインが既知の電子顕微鏡像を例に説明した。電子顕微鏡像の元素組成およびドメインの情報がない場合(例えば試料の結晶情報が不明な場合)であっても元素識別が可能である。例えば、ユーザが結晶情報を開示せずにインターネット等を介し、電子顕微鏡像を元素識別装置30に送信する。元素識別装置30は自動的に元素を識別し、結果をインターネット等を介しユーザの画像出力装置36に送信する。このように、ユーザは結晶情報を一切開示しなくとも、元素識別画像を得ることができる。実施例1を用い、このようなクラウドサービス等を行なうこともできる。   As described above, the electron microscope image in which the element type and domain are known as shown in FIG. 7 has been described as an example. Even when there is no element composition and domain information in the electron microscope image (for example, when the crystal information of the sample is unknown), element identification is possible. For example, the user transmits an electron microscope image to the element identification device 30 via the Internet or the like without disclosing crystal information. The element identification device 30 automatically identifies the element and transmits the result to the image output device 36 of the user via the Internet or the like. Thus, the user can obtain an element identification image without disclosing any crystal information. Such a cloud service or the like can also be performed using the first embodiment.

実施例1によれば、図6のステップS30のように、演算部22は、走査透過型電子顕微鏡像における原子に対応した信号強度に基づき像における原子位置をそれぞれ複数のグループ52および54に対応させる。ステップS32のように、複数のグループ52および54それぞれ内において、原子位置を全て通る最短ルート63および65を算出する。ステップS34およびS36のように、演算部22は、算出結果に基づき複数のグループ52および54間で原子位置の一部を入れ換える。図5のステップS22およびS24のように、元素識別部24は、入れ換え後の複数のグループ52および54内の原子位置をそれぞれ異なる元素に識別する。これより、TEM法とEELS法またはEDX法とを併用する方法に比べ、元素識別時間を短縮できる。これにより、例えば試料のドリフトまたは電子線照射によるダメージを抑制できる。また、EELS法またはEDX法のような高価な装置を用いなくてもよくなる。   According to the first embodiment, as shown in step S30 of FIG. 6, the calculation unit 22 corresponds the atomic positions in the image to the plurality of groups 52 and 54 based on the signal intensity corresponding to the atoms in the scanning transmission electron microscope image, respectively. Let As in step S32, the shortest routes 63 and 65 passing through all the atom positions are calculated in each of the plurality of groups 52 and 54. As in steps S34 and S36, the computing unit 22 interchanges some of the atomic positions between the plurality of groups 52 and 54 based on the calculation result. As in steps S22 and S24 of FIG. 5, the element identifying unit 24 identifies the atomic positions in the plurality of groups 52 and 54 after replacement as different elements. Thereby, element identification time can be shortened compared with the method using TEM method, EELS method, or EDX method together. Thereby, the damage by the drift of a sample or electron beam irradiation can be suppressed, for example. Further, it is not necessary to use an expensive apparatus such as the EELS method or the EDX method.

実施例1では、電子顕微鏡像としてHAADF STEM像を説明したが、電子顕微鏡像は、HABF STEM像でもよい。電子顕微鏡像としてHABF STEM像を用いる場合、ステップS12の前に、信号強度に−1を乗ずる等の反転処理を行なう。電子顕微鏡像は、信号強度が原子番号に依存していればよい。   In Example 1, although the HAADF STEM image was demonstrated as an electron microscopic image, a HABF STEM image may be sufficient as an electron microscopic image. When the HABF STEM image is used as the electron microscope image, inversion processing such as multiplying the signal intensity by −1 is performed before step S12. In the electron microscope image, the signal intensity may be dependent on the atomic number.

図5のステップS12のように、演算部22は、原子に対応した信号強度がピークとなる強度に基づき電子顕微鏡像における原子位置をそれぞれ複数のグループ52および54に対応させる。これにより、図10のように、信号強度による元素の識別がより容易になる。図8を用い元素識別が可能な場合、ピーク強度を抽出せず、図8の信号強度から元素識別処理を行なってもよい。   As in step S <b> 12 of FIG. 5, the calculation unit 22 associates the atomic positions in the electron microscope image with the plurality of groups 52 and 54 based on the intensity at which the signal intensity corresponding to the atoms peaks. Thereby, as shown in FIG. 10, the identification of the element by the signal intensity becomes easier. When element identification is possible using FIG. 8, the element identification processing may be performed from the signal intensity of FIG. 8 without extracting the peak intensity.

図6のステップS30および図11のように、演算部22は、元素の識別が困難な信号強度の範囲50を複数の範囲に分割する。演算部22は、信号強度が複数の範囲内の原子位置をそれぞれ複数のグループ52および54に対応させる。これにより、元素識別が困難な範囲内の元素を識別できる。   As shown in step S30 of FIG. 6 and FIG. 11, the calculation unit 22 divides the signal intensity range 50 in which element identification is difficult into a plurality of ranges. The computing unit 22 associates atomic positions having signal intensities within a plurality of ranges with a plurality of groups 52 and 54, respectively. Thereby, the element in the range where element identification is difficult can be identified.

図6のステップS30において分割するグループの数は任意であるが、複数のグループの数は、元素の識別が困難な信号強度の範囲50において識別が困難な元素の数と同じであることが好ましい。これにより、グループを元素に1対1で対応させることができる。   Although the number of groups to be divided in step S30 in FIG. 6 is arbitrary, the number of the plurality of groups is preferably the same as the number of elements that are difficult to identify in the signal intensity range 50 where it is difficult to identify the elements. . Thereby, a group can be made to correspond to an element on a one-to-one basis.

ステップS34およびS36のように、演算部22は、最短ルート上において隣接する原子位置の距離が第1距離(例えば最近接原子間距離)より長い原子位置を他のグループに移動する。最短ルートにおいて隣接する原子位置の距離が最近接原子間距離より長い場合、原子位置の属するグループが正しくない可能性がある。そこで、この原子位置を他のグループに移動する。これにより、元素の識別が容易となる。第1距離は最近接原子間距離以上であればよい。第1距離が最近接原子間距離の2倍以上のとき、原子位置が属するグループが正しくないときであっても、原子位置を他のグループに移動しないことが起こりえる。よって、第1距離は最近接原子間距離以上かつ最近接原子間距離の2倍より短い一定の距離であることが好ましい。また、最近接原子間距離は、結晶方位によって異なってもよい。例えば、図14(b)のグループ54内のルート65では、斜め方向の最近接原子間距離は横方向および縦方向の最近接原子間距離より大きくてもよい。   As in steps S34 and S36, the calculation unit 22 moves an atomic position where the distance between adjacent atomic positions on the shortest route is longer than a first distance (for example, the closest interatomic distance) to another group. When the distance between adjacent atom positions in the shortest route is longer than the distance between the nearest atoms, the group to which the atom position belongs may not be correct. Therefore, this atomic position is moved to another group. This facilitates element identification. The first distance may be equal to or greater than the nearest interatomic distance. When the first distance is more than twice the nearest interatomic distance, even if the group to which the atomic position belongs is not correct, the atomic position may not be moved to another group. Therefore, the first distance is preferably a constant distance that is equal to or greater than the nearest interatomic distance and shorter than twice the nearest interatomic distance. The closest interatomic distance may vary depending on the crystal orientation. For example, in the route 65 in the group 54 in FIG. 14B, the closest interatomic distance in the oblique direction may be larger than the closest interatomic distance in the horizontal direction and the vertical direction.

図16(a)において、ステップS38およびS40を行なわない場合、図16(b)に収束するまで、ステップS32からS42を繰り返すことになる。そこで、ステップS38およびS40のように、演算部22は、最短ルート上において隣接する原子位置の距離が第2距離(例えば最近接原子間距離の2倍)より長い経路で挟まれた原子位置を他のグループに移動する。図16(a)のグループ52のように、隣接する原子位置の距離が最近接原子間距離の2倍より長い経路63bで挟まれた原子位置のブロック67は、属するグループが正しくない可能性が高い。そこで、ブロック67ごとグループを移動する。これにより、ステップS38およびS40を行なわない場合に比べ、演算時間を短くできる。第2距離は最近接原子間距離の2倍以上の一定の距離であればよい。   In FIG. 16A, when steps S38 and S40 are not performed, steps S32 to S42 are repeated until convergence to FIG. 16B. Therefore, as in steps S38 and S40, the calculation unit 22 determines an atomic position sandwiched by a path in which the distance between adjacent atomic positions on the shortest route is longer than a second distance (for example, twice the distance between nearest atoms). Move to another group. As in the group 52 in FIG. 16A, the block 67 at the atomic position sandwiched by the path 63b in which the distance between adjacent atomic positions is longer than twice the distance between the nearest atoms may belong to an incorrect group. high. Therefore, the group is moved together with the block 67. Thereby, the calculation time can be shortened compared with the case where steps S38 and S40 are not performed. The second distance may be a certain distance that is at least twice the distance between the nearest atoms.

元素識別部24は、演算部22がグループ52および54それぞれ内において、ステップS32からS42を繰り返し、ステップS42のように原子位置の一部を入れ換えが実行されなくなったときに、元素に識別する。これにより、演算の収束後に元素の識別を行なうことができる。   In the groups 52 and 54, the element identification unit 24 repeats steps S32 to S42 in each of the groups 52 and 54, and identifies elements as elements when replacement of part of atomic positions is no longer performed as in step S42. Thereby, the element can be identified after the convergence of the calculation.

図7のように、電子顕微鏡像は複数のドメイン56および58を含む。このような場合、元素を識別することにより、図16(b)のように、ドメイン56および58の境界55を決定できる。よって、ドメイン56および58を抽出することができる。   As shown in FIG. 7, the electron microscope image includes a plurality of domains 56 and 58. In such a case, by identifying the element, the boundary 55 between the domains 56 and 58 can be determined as shown in FIG. Therefore, the domains 56 and 58 can be extracted.

図5のステップS16のように、元素識別部24は、元素の識別が困難な信号強度の範囲50以外の信号強度から別の元素を識別する。これにより、別の元素を簡単に識別できる。   As in step S <b> 16 of FIG. 5, the element identification unit 24 identifies another element from a signal intensity other than the signal intensity range 50 in which it is difficult to identify the element. Thereby, another element can be easily identified.

図5のステップS24のように、画像生成部26は、元素識別部24が識別した元素に基づき、電子顕微鏡像に対応し元素を識別した画像を生成する。これにより、元素識別結果を視覚的に表現できる。   As in step S <b> 24 of FIG. 5, the image generation unit 26 generates an image that identifies the element corresponding to the electron microscope image based on the element identified by the element identification unit 24. Thereby, an element identification result can be expressed visually.

以上、本発明の実施例について詳述したが、本発明は係る特定の実施例に限定されるものではなく、特許請求の範囲に記載された本発明の要旨の範囲内において、種々の変形・変更が可能である。   Although the embodiments of the present invention have been described in detail above, the present invention is not limited to such specific embodiments, and various modifications and changes can be made within the scope of the gist of the present invention described in the claims. It can be changed.

なお、以上の説明に関して更に以下の付記を開示する。
(付記1)走査透過型電子顕微鏡像における原子に対応した信号強度に基づき前記走査透過型電子顕微鏡像における原子位置をそれぞれ複数のグループに対応させ、前記複数のグループそれぞれ内において、原子位置を全て通る最短ルートを算出し、算出結果に基づき前記複数のグループ間で原子位置の一部を入れ換える演算部と、前記演算部が前記原子位置の一部を入れ換えた後、前記複数のグループ内の原子位置をそれぞれ異なる元素に識別する元素識別部と、を具備することを特徴とする元素識別装置。
(付記2)前記演算部は、前記原子に対応した信号強度がピークとなる強度に基づき前記走査透過型電子顕微鏡像における原子位置をそれぞれ複数のグループに対応させることを特徴とする付記1に記載の元素識別装置。
(付記3)前記演算部は、元素の識別が困難な信号強度の範囲を複数の範囲に分割し、信号強度が前記複数の範囲内の原子位置をそれぞれ前記複数のグループに対応させることを特徴とする付記1または2記載の元素識別装置。
(付記4)前記複数のグループの数は、前記元素の識別が困難な信号強度の範囲において識別が困難な元素の数と同じであることを特徴とする付記3記載の元素識別装置。
(付記5)前記演算部は、前記最短ルート上において隣接する原子位置の距離が最近接原子間距離以上の第1距離より長い原子位置を他のグループに移動することを特徴とする付記1から4のいずれか一項記載の元素識別装置。
(付記6)前記演算部は、前記最短ルート上において隣接する原子位置の距離が最近接原子間距離の2倍以上の第2距離より長い経路で挟まれた原子位置を他のグループに移動することを特徴とする付記5記載の元素識別装置。
(付記7)前記元素識別部は、前記演算部が前記複数のグループそれぞれ内において、前記最短ルートを算出し、前記複数のグループ間で原子位置の一部を入れ換える演算を繰り返し、前記原子位置の一部の入れ換えが実行されなくなったときに、前記複数のグループ内の原子位置をそれぞれ異なる元素に識別することを特徴とする付記1から6のいずれか一項記載の元素識別装置。
(付記8)前記走査透過型電子顕微鏡像は複数のドメインを含むことを特徴とする付記1から7のいずれか一項記載の元素識別装置。
(付記9)前記元素識別部は、前記元素の識別が困難な信号強度の範囲以外の信号強度から別の元素を識別することを特徴とする付記3または4記載の元素識別装置。
(付記10)前記元素識別部が識別した元素に基づき、前記走査透過型電子顕微鏡像に対応し元素を識別した画像を生成する画像生成部を具備することを特徴とする付記1から9のいずれか一項記載の元素識別装置。
(付記11)コンピュータに、走査透過型電子顕微鏡像における原子に対応した信号強度に基づき前記走査透過型電子顕微鏡像における原子位置をそれぞれ複数のグループに対応させ、前記複数のグループそれぞれ内において、原子位置を全て通る最短ルートを算出し、算出結果に基づき前記複数のグループ間で原子位置の一部を入れ換え、前記原子位置の一部を入れ換えた後、前記複数のグループ内の原子位置をそれぞれ異なる元素に識別させることを特徴とする元素識別プログラム。
(付記12)コンピュータに実行させる元素識別方法であって、走査透過型電子顕微鏡像における原子に対応した信号強度に基づき前記走査透過型電子顕微鏡像における原子位置をそれぞれ複数のグループに対応させ、前記複数のグループそれぞれ内において、原子位置を全て通る最短ルートを算出し、算出結果に基づき前記複数のグループ間で原子位置の一部を入れ換え、前記原子位置の一部を入れ換えた後、前記複数のグループ内の原子位置をそれぞれ異なる元素に識別させることを特徴とする元素識別方法。
In addition, the following additional notes are disclosed regarding the above description.
(Appendix 1) Based on the signal intensity corresponding to the atoms in the scanning transmission electron microscope image, the atomic positions in the scanning transmission electron microscope image are respectively associated with a plurality of groups, and all the atomic positions in each of the plurality of groups An arithmetic unit that calculates a shortest route to pass, and replaces part of the atomic positions between the plurality of groups based on the calculation result; and after the arithmetic unit replaces part of the atomic positions, the atoms in the plurality of groups An element identification device comprising: an element identification unit for identifying each position by a different element.
(Supplementary note 2) The supplementary note 1, wherein the calculation unit associates atomic positions in the scanning transmission electron microscope image with a plurality of groups based on the intensity at which the signal intensity corresponding to the atoms reaches a peak. Element identification device.
(Additional remark 3) The said calculating part divides | segments the range of the signal intensity | strength where it is difficult to identify an element into a some range, and signal intensity makes each atomic position in the said some range correspond to the said some group, respectively. The element identification device according to Supplementary Note 1 or 2.
(Supplementary note 4) The element identification device according to supplementary note 3, wherein the number of the plurality of groups is the same as the number of elements that are difficult to identify in a signal intensity range in which the elements are difficult to identify.
(Additional remark 5) The said calculating part moves the atomic position where the distance of the adjacent atomic position on the said shortest route is longer than 1st distance more than the nearest interatomic distance to another group from the additional remark 1 characterized by the above-mentioned. 5. The element identification device according to any one of 4.
(Additional remark 6) The said calculating part moves the atomic position on which the distance of the adjacent atomic position on the said shortest route is pinched | interposed by the path | route longer than 2nd distance more than twice the distance between nearest atoms to another group The element identification device according to supplementary note 5, wherein
(Additional remark 7) The said element identification part repeats the calculation which the said calculating part calculates the said shortest route in each of these groups, and replaces a part of atomic position between these groups, The element identification device according to any one of appendices 1 to 6, wherein when a part of the replacement is not performed, the atomic positions in the plurality of groups are identified by different elements.
(Supplementary note 8) The element identification device according to any one of supplementary notes 1 to 7, wherein the scanning transmission electron microscope image includes a plurality of domains.
(Supplementary note 9) The element identification device according to supplementary note 3 or 4, wherein the element identification unit identifies another element from a signal intensity outside a signal intensity range in which it is difficult to identify the element.
(Supplementary note 10) Any one of Supplementary notes 1 to 9, further comprising an image generation unit that generates an image that identifies the element corresponding to the scanning transmission electron microscope image based on the element identified by the element identification unit. An element identification device according to claim 1.
(Supplementary Note 11) A computer is caused to correspond to each atom position in the scanning transmission electron microscope image based on the signal intensity corresponding to the atom in the scanning transmission electron microscope image, and in each of the plurality of groups, an atom The shortest route that passes through all the positions is calculated, and based on the calculation result, a part of the atomic positions are replaced between the plurality of groups, and after a part of the atomic positions are replaced, the atomic positions in the plurality of groups are different from each other. An element identification program characterized by allowing an element to be identified.
(Supplementary note 12) An element identification method to be executed by a computer, wherein the atomic positions in the scanning transmission electron microscope image are respectively associated with a plurality of groups based on the signal intensity corresponding to the atoms in the scanning transmission electron microscope image, Within each of the plurality of groups, the shortest route that passes through all the atomic positions is calculated, and based on the calculation result, a part of the atomic positions is replaced between the plurality of groups, and after a part of the atomic positions are replaced, A method for identifying an element, characterized in that each atom position in a group is identified by a different element.

10 コンピュータ
20 像取得部
22 演算部
24 元素識別部
26 画像生成部
30 元素識別装置
32 電子顕微鏡
34 像入力装置
36 画像出力装置
50 識別困難範囲
52、54 グループ
55 境界
56、58 ドメイン
62、64 元素
64a、64b 原子位置
63、65 ルート
63a、63b、65a 経路
DESCRIPTION OF SYMBOLS 10 Computer 20 Image acquisition part 22 Operation part 24 Element identification part 26 Image generation part 30 Element identification apparatus 32 Electron microscope 34 Image input apparatus 36 Image output apparatus 50 Difficult-of-identification range 52, 54 Group 55 Boundary 56, 58 Domain 62, 64 Element 64a, 64b Atomic position 63, 65 Route 63a, 63b, 65a Route

Claims (10)

走査透過型電子顕微鏡像における原子に対応した信号強度に基づき前記走査透過型電子顕微鏡像における原子位置をそれぞれ複数のグループに対応させ、前記複数のグループそれぞれ内において、原子位置を全て通る最短ルートを算出し、算出結果に基づき前記複数のグループ間で原子位置の一部を入れ換える演算部と、
前記演算部が前記原子位置の一部を入れ換えた後、前記複数のグループ内の原子位置をそれぞれ異なる元素に識別する元素識別部と、
を具備することを特徴とする元素識別装置。
Based on the signal intensity corresponding to the atoms in the scanning transmission electron microscope image, the atomic positions in the scanning transmission electron microscope image are made to correspond to a plurality of groups, respectively, and within each of the plurality of groups, the shortest route that passes through all the atom positions A calculation unit that calculates and replaces part of atomic positions between the plurality of groups based on the calculation result;
After the arithmetic unit replaces a part of the atomic positions, an element identifying unit that identifies the atomic positions in the plurality of groups as different elements, and
An element identification device comprising:
前記演算部は、前記原子に対応した信号強度がピークとなる強度に基づき前記走査透過型電子顕微鏡像における原子位置をそれぞれ複数のグループに対応させることを特徴とする請求項1に記載の元素識別装置。   2. The element identification according to claim 1, wherein the calculation unit associates each atomic position in the scanning transmission electron microscope image with a plurality of groups based on an intensity at which a signal intensity corresponding to the atom reaches a peak. apparatus. 前記演算部は、元素の識別が困難な信号強度の範囲を複数の範囲に分割し、信号強度が前記複数の範囲内の原子位置をそれぞれ前記複数のグループに対応させることを特徴とする請求項1または2記載の元素識別装置。   The calculation unit divides a range of signal intensity where it is difficult to identify an element into a plurality of ranges, and the signal intensity corresponds to atomic positions in the plurality of ranges respectively corresponding to the plurality of groups. 3. The element identification device according to 1 or 2. 前記複数のグループの数は、前記元素の識別が困難な信号強度の範囲において識別が困難な元素の数と同じであることを特徴とする請求項3記載の元素識別装置。   4. The element identification device according to claim 3, wherein the number of the plurality of groups is the same as the number of elements that are difficult to identify in a signal intensity range in which the elements are difficult to identify. 前記演算部は、前記最短ルート上において隣接する原子位置の距離が最近接原子間距離以上の第1距離より長い原子位置を他のグループに移動することを特徴とする請求項1から4のいずれか一項記載の元素識別装置。   5. The calculation unit according to claim 1, wherein the calculation unit moves an atomic position having a distance between adjacent atomic positions on the shortest route that is longer than a first distance that is equal to or greater than the nearest interatomic distance to another group. An element identification device according to claim 1. 前記演算部は、前記最短ルート上において隣接する原子位置の距離が最近接原子間距離の2倍以上の第2距離より長い経路で挟まれた原子位置を他のグループに移動することを特徴とする請求項5記載の元素識別装置。   The arithmetic unit moves an atomic position sandwiched by a path longer than a second distance in which the distance between adjacent atomic positions on the shortest route is twice or more the closest interatomic distance to another group, The element identification device according to claim 5. 前記元素識別部は、前記演算部が前記複数のグループそれぞれ内において、前記最短ルートを算出し、前記複数のグループ間で原子位置の一部を入れ換える演算を繰り返し、前記原子位置の一部の入れ換えが実行されなくなったときに、前記複数のグループ内の原子位置をそれぞれ異なる元素に識別することを特徴とする請求項1から6のいずれか一項記載の元素識別装置。   The element identification unit repeats an operation in which the calculation unit calculates the shortest route in each of the plurality of groups, and replaces part of the atomic positions between the plurality of groups, thereby replacing part of the atomic positions. 7. The element identification device according to claim 1, wherein when the operation is stopped, the atomic positions in the plurality of groups are identified as different elements. 前記走査透過型電子顕微鏡像は複数のドメインを含むことを特徴とする請求項1から7のいずれか一項記載の元素識別装置。   The element identification apparatus according to claim 1, wherein the scanning transmission electron microscope image includes a plurality of domains. コンピュータに、
走査透過型電子顕微鏡像における原子に対応した信号強度に基づき前記走査透過型電子顕微鏡像における原子位置をそれぞれ複数のグループに対応させ、
前記複数のグループそれぞれ内において、原子位置を全て通る最短ルートを算出し、
算出結果に基づき前記複数のグループ間で原子位置の一部を入れ換え、
前記原子位置の一部を入れ換えた後、前記複数のグループ内の原子位置をそれぞれ異なる元素に識別させることを特徴とする元素識別プログラム。
On the computer,
Based on the signal intensity corresponding to the atoms in the scanning transmission electron microscope image, the atomic positions in the scanning transmission electron microscope image correspond to a plurality of groups,
Within each of the plurality of groups, calculate the shortest route through all the atom positions,
Based on the calculation results, replace some of the atomic positions between the plurality of groups,
An element identification program, wherein after replacing a part of the atomic positions, the atomic positions in the plurality of groups are identified by different elements.
コンピュータに実行させる元素識別方法であって、
走査透過型電子顕微鏡像における原子に対応した信号強度に基づき前記走査透過型電子顕微鏡像における原子位置をそれぞれ複数のグループに対応させ、
前記複数のグループそれぞれ内において、原子位置を全て通る最短ルートを算出し、
算出結果に基づき前記複数のグループ間で原子位置の一部を入れ換え、
前記原子位置の一部を入れ換えた後、前記複数のグループ内の原子位置をそれぞれ異なる元素に識別させることを特徴とする元素識別方法。
An element identification method executed by a computer,
Based on the signal intensity corresponding to the atoms in the scanning transmission electron microscope image, the atomic positions in the scanning transmission electron microscope image correspond to a plurality of groups,
Within each of the plurality of groups, calculate the shortest route through all the atom positions,
Based on the calculation results, replace some of the atomic positions between the plurality of groups,
After exchanging a part of the atomic positions, the element identifying method includes identifying different atomic positions in the plurality of groups.
JP2015233915A 2015-11-30 2015-11-30 Element identification device, element identification program and element identification method Active JP6551198B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2015233915A JP6551198B2 (en) 2015-11-30 2015-11-30 Element identification device, element identification program and element identification method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2015233915A JP6551198B2 (en) 2015-11-30 2015-11-30 Element identification device, element identification program and element identification method

Publications (2)

Publication Number Publication Date
JP2017101971A true JP2017101971A (en) 2017-06-08
JP6551198B2 JP6551198B2 (en) 2019-07-31

Family

ID=59017907

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2015233915A Active JP6551198B2 (en) 2015-11-30 2015-11-30 Element identification device, element identification program and element identification method

Country Status (1)

Country Link
JP (1) JP6551198B2 (en)

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2013011474A (en) * 2011-06-28 2013-01-17 Akita Univ EVALUATION METHOD OF FINE TISSUE STRUCTURE OF Mg-Li ALLOY
JP2013025509A (en) * 2011-07-20 2013-02-04 Fujitsu Ltd Program, method, and device for determining traveling route
KR101290869B1 (en) * 2012-04-26 2013-07-29 세종대학교산학협력단 Apparatus for monitoring electronic structures in atom and method for thereof
JP2015056253A (en) * 2013-09-11 2015-03-23 富士通株式会社 Sample analyzing program, sample analyzing apparatus, and sample analyzing method
JP2015206102A (en) * 2014-04-23 2015-11-19 株式会社ノリタケカンパニーリミテド Platinum hollow nanoparticle, catalyst carrying the particle, and method for producing the catalyst

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2013011474A (en) * 2011-06-28 2013-01-17 Akita Univ EVALUATION METHOD OF FINE TISSUE STRUCTURE OF Mg-Li ALLOY
JP2013025509A (en) * 2011-07-20 2013-02-04 Fujitsu Ltd Program, method, and device for determining traveling route
KR101290869B1 (en) * 2012-04-26 2013-07-29 세종대학교산학협력단 Apparatus for monitoring electronic structures in atom and method for thereof
JP2015056253A (en) * 2013-09-11 2015-03-23 富士通株式会社 Sample analyzing program, sample analyzing apparatus, and sample analyzing method
JP2015206102A (en) * 2014-04-23 2015-11-19 株式会社ノリタケカンパニーリミテド Platinum hollow nanoparticle, catalyst carrying the particle, and method for producing the catalyst

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
DILLON ET AL.: "Direct Observation of Multilayer Adsorption on Alumina Grain Boundaries", J. AM. CERAM. SOC., vol. 90, 3, JPN6019020506, 2007, pages 996 - 998, ISSN: 0004048047 *
ITAKURA ET AL.: "Atomic-Resolution X-ray Energy-Dispersive Spectroscopy Chemical Mapping of Substitutional Dy Atoms i", JAPANESE JOURNAL OF APPLIED PHYSICS, vol. 52, JPN6019020507, 2013, pages 050201, ISSN: 0004048048 *
柴田直哉他: "埋もれた界面中の原子1個の走査型透過電子顕微鏡直接観察", J. VAC. SOC. JPN, vol. 54, 4, JPN6019020505, 2011, pages 270 - 274, ISSN: 0004048046 *

Also Published As

Publication number Publication date
JP6551198B2 (en) 2019-07-31

Similar Documents

Publication Publication Date Title
TWI647442B (en) Defect classification device and defect classification method
Durmaz et al. A deep learning approach for complex microstructure inference
JP6731370B2 (en) Image processing system and computer program for performing image processing
KR20160078955A (en) Method and system for classifying and identifying individual cells in a microscopy image
US20150362446A1 (en) Phase Analyzer, Phase Analysis Method, and Surface Analyzer
CN108762740B (en) Page data generation method and device and electronic equipment
US20180342087A1 (en) Image Processing Apparatus, Surface Analyzer, and Image Processing Method
US20230237777A1 (en) Information processing apparatus, learning apparatus, image recognition apparatus, information processing method, learning method, image recognition method, and non-transitory-computer-readable storage medium
WO2019155850A1 (en) Autofocus method for scanning electron microscope
JP2021184766A (en) Cell observation device and cell observation method
US11508047B2 (en) Charged particle microscope device and wide-field image generation method
US10055638B2 (en) Image processing device, imaging device, microscope device, image processing method, and image processing program
JP6551198B2 (en) Element identification device, element identification program and element identification method
JP2012142143A (en) Charged particle beam device and observation image generating method
JP2014240780A (en) Sample structure analysis method, transmission electron microscope, and program
JP2020071619A (en) Optical-character-recognition-result correction assisting apparatus and correction assisting program
WO2020158261A1 (en) Image-matching determination method, image-matching determination device, and computer-readable recording medium on which program for executing image-matching determination method on computer is recorded
CN107020845A (en) The image processing apparatus and image processing method of time needed for estimation print processing
KR20200117880A (en) Using Images from Secondary Microscope Detectors to Automatically Generate Labeled Images from Primary Microscope Detectors
KR20220040466A (en) charged particle beam device
JP6595856B2 (en) Charged particle apparatus and measurement method
JP2013114893A (en) Scanning electron microscope, and automatic focusing method for the same
WO2021166142A1 (en) Pattern matching device, pattern measurement system, and non-transitory computer-readable medium
JP2020003837A (en) Identification apparatus and identification method
JP6490852B1 (en) Unread determination threshold setting method and unread determination threshold setting apparatus

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20180810

TRDD Decision of grant or rejection written
A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

Effective date: 20190604

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20190531

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20190617

R150 Certificate of patent or registration of utility model

Ref document number: 6551198

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150