JP2021043828A - Analysis method of metal solidified crystal structure, and analyzing device of metal solidified crystal structure - Google Patents
Analysis method of metal solidified crystal structure, and analyzing device of metal solidified crystal structure Download PDFInfo
- Publication number
- JP2021043828A JP2021043828A JP2019166769A JP2019166769A JP2021043828A JP 2021043828 A JP2021043828 A JP 2021043828A JP 2019166769 A JP2019166769 A JP 2019166769A JP 2019166769 A JP2019166769 A JP 2019166769A JP 2021043828 A JP2021043828 A JP 2021043828A
- Authority
- JP
- Japan
- Prior art keywords
- crystal
- crystal structure
- observation image
- average shape
- region
- 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.)
- Pending
Links
- 239000013078 crystal Substances 0.000 title claims abstract description 298
- 238000004458 analytical method Methods 0.000 title claims abstract description 99
- 239000002184 metal Substances 0.000 title claims abstract description 80
- 238000000034 method Methods 0.000 claims abstract description 87
- 238000004364 calculation method Methods 0.000 claims abstract description 55
- 238000009826 distribution Methods 0.000 claims abstract description 53
- 230000007246 mechanism Effects 0.000 claims description 57
- 239000007787 solid Substances 0.000 claims description 7
- 230000008569 process Effects 0.000 description 19
- 230000005484 gravity Effects 0.000 description 10
- 238000013500 data storage Methods 0.000 description 9
- 230000007704 transition Effects 0.000 description 8
- 238000007405 data analysis Methods 0.000 description 7
- 238000010586 diagram Methods 0.000 description 7
- 239000007769 metal material Substances 0.000 description 7
- 230000008901 benefit Effects 0.000 description 6
- 230000008859 change Effects 0.000 description 6
- 238000004519 manufacturing process Methods 0.000 description 6
- 238000006243 chemical reaction Methods 0.000 description 4
- 238000007711 solidification Methods 0.000 description 4
- 230000008023 solidification Effects 0.000 description 4
- 230000015572 biosynthetic process Effects 0.000 description 2
- 239000002131 composite material Substances 0.000 description 2
- 238000005260 corrosion Methods 0.000 description 2
- 230000007797 corrosion Effects 0.000 description 2
- 230000002950 deficient Effects 0.000 description 2
- 238000001887 electron backscatter diffraction Methods 0.000 description 2
- 238000011156 evaluation Methods 0.000 description 2
- 238000010191 image analysis Methods 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- 238000004445 quantitative analysis Methods 0.000 description 2
- 238000003786 synthesis reaction Methods 0.000 description 2
- 230000002194 synthesizing effect Effects 0.000 description 2
- 238000000342 Monte Carlo simulation Methods 0.000 description 1
- 238000002003 electron diffraction Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 230000003595 spectral effect Effects 0.000 description 1
- 238000001308 synthesis method Methods 0.000 description 1
Images
Landscapes
- Image Analysis (AREA)
Abstract
Description
本発明は、金属凝固体の結晶組織を解析する技術に関し、特に、金属凝固体の柱状晶−等軸晶遷移による柱状晶領域と等軸晶領域とを定量的に解析する方法、および該方法を実施するための金属凝固体結晶組織の解析装置に関するものである。 The present invention relates to a technique for analyzing the crystal structure of a solidified metal body, and in particular, a method for quantitatively analyzing a columnar crystal region and an equiaxed crystal region due to a columnar crystal-equal axis crystal transition of a solid metal solid body, and the method It relates to an analyzer for analyzing a solidified metal crystal structure for carrying out the above.
種々の機械装置で使用される金属材料は、当該機械装置の稼働環境で要求される様々な性状(例えば、機械的特性、熱的特性、耐食性)を満たす必要がある。一般的に、耐食性や熱的特性は金属材料の組成との関連性が強く、機械的特性は金属材料の微細組織との関連性が強い。また、金属材料の微細組織は製造方法との関連性が極めて強い。 The metal material used in various mechanical devices needs to satisfy various properties (for example, mechanical properties, thermal properties, corrosion resistance) required in the operating environment of the mechanical device. In general, corrosion resistance and thermal properties are strongly related to the composition of the metal material, and mechanical properties are strongly related to the microstructure of the metal material. Moreover, the fine structure of the metal material is extremely closely related to the manufacturing method.
すなわち、望ましい機械的特性を有する金属材料を得るためには微細組織の制御が必要不可欠であり、望ましい微細組織を得るためには製造方法との関連性を明確化した上で、製造方法を制御する必要がある。 That is, control of the microstructure is indispensable for obtaining a metal material having desirable mechanical properties, and control of the manufacturing method is performed after clarifying the relationship with the manufacturing method in order to obtain the desired microstructure. There is a need to.
金属材料の微細組織を評価する方法は、現在までに数多く提案されているが、一つの評価項目を単純数値化したものが多い(例えば、平均結晶粒径、所定の相の占有率など)。従来の単純数値化した一項目指標は、微細組織の観察領域内全体がほぼ均等である場合には有効である。しかしながら、微細組織の観察領域内に大きな分布や偏りがある場合には、単純数値化した一項目指標で分布や偏りがある状態を解析するのが困難であった。 Many methods for evaluating the microstructure of metallic materials have been proposed so far, but many of them are simple numerical values of one evaluation item (for example, average crystal grain size, occupancy of a predetermined phase, etc.). The conventional simple numerical one-item index is effective when the entire observation area of the microstructure is almost even. However, when there is a large distribution or bias in the observation area of the microstructure, it is difficult to analyze the state of the distribution or bias using a simple numerical one-item index.
一方、近年は画像解析の技術が飛躍的に向上しており、微細組織観察の領域全体を二次元的に解析評価できるようになってきた。例えば、特許文献1(特開2004-192055)には、広視野をカバーする仮想環境の合成画像を、透視投影画像を元画像として投影系の内部パラメータにより設定される変換表に基づいて画素変換することにより合成する手順からなる広視野分割合成方法であって、前記変換表を設定する少なくとも一つの内部パラメータに対し、オンライン変更指示がなされたとき第2 の画素変換手順によって簡易変換表を作成し広視野画像を合成する手順を特徴とする広視野分割合成方法、が開示されている。 On the other hand, in recent years, the technique of image analysis has been dramatically improved, and it has become possible to analyze and evaluate the entire area of microstructure observation in two dimensions. For example, in Patent Document 1 (Japanese Patent Laid-Open No. 2004-192055), a composite image of a virtual environment covering a wide field of view is converted into pixels based on a conversion table set by internal parameters of the projection system using a perspective projection image as an original image. This is a wide-field division synthesis method consisting of a procedure for synthesizing by performing the above, and when an online change instruction is given to at least one internal parameter for setting the conversion table, a simple conversion table is created by the second pixel conversion procedure. A wide-field split-composite method, which comprises a procedure for synthesizing a wide-field image, is disclosed.
金属材料は、多くの場合で溶融凝固によって成形された凝固体をベースとして利用される。金属凝固体で望ましい機械的特性を達成するためには、前述したように、金属凝固体における結晶組織を制御する必要があり、金属凝固体の結晶組織を制御するために、当該結晶組織と製造方法との関連性を明確化した上で、製造方法を制御する必要がある。すなわち、製造方法を制御する前に、金属凝固体の結晶組織を明確化する必要がある。 Metallic materials are often used as a base for solidified bodies formed by melt solidification. As described above, in order to achieve the desired mechanical properties of the solidified metal, it is necessary to control the crystal structure of the solidified metal, and in order to control the crystal structure of the solidified metal, the crystal structure and production thereof. It is necessary to control the manufacturing method after clarifying the relationship with the method. That is, it is necessary to clarify the crystal structure of the solidified metal before controlling the production method.
ここで、金属凝固体では、その凝固過程において、しばしば結晶組織内に柱状晶−等軸晶遷移が起こることが知られている。柱状晶は、等軸晶に比して種々の性状に異方性を有することから、金属凝固体全体としての性状を把握するためには、柱状晶領域と等軸晶領域とを定量的に解析する必要がある。 Here, in a solidified metal, it is known that a columnar-equal axis transition often occurs in the crystal structure during the solidification process. Since columnar crystals have anisotropy in various properties compared to equiaxed crystals, the columnar crystal region and equiaxed crystal region are quantitatively determined in order to grasp the properties of the entire metal solidified body. Need to be analyzed.
金属凝固体の結晶組織における柱状晶領域/等軸晶領域の定量的な解析を考えると、従来の単純数値化した一項目指標で全体を把握することは困難である。また、本発明者等の検討によると、特許文献1の画像解析技術でも、金属凝固体の結晶組織における柱状晶領域/等軸晶領域の定量的な解析は困難であると考えられた。 Considering the quantitative analysis of the columnar crystal region / equiaxed crystal region in the crystal structure of the solidified metal, it is difficult to grasp the whole by the conventional simple numerical one-item index. Further, according to the study by the present inventors, it is considered that even with the image analysis technique of Patent Document 1, it is difficult to quantitatively analyze the columnar crystal region / equiaxed crystal region in the crystal structure of the solidified metal.
したがって、本発明の目的は、金属凝固体の結晶組織における柱状晶領域と等軸晶領域とを定量的に解析する方法、および該方法を実施するための金属凝固体結晶組織の解析装置を提供することにある。 Therefore, an object of the present invention is to provide a method for quantitatively analyzing a columnar crystal region and an equiaxed crystal region in the crystal structure of a metal solidified body, and an apparatus for analyzing the metal solidified body crystal structure for carrying out the method. To do.
(I)本発明の一態様は、金属凝固体の結晶組織を解析する方法であって、
前記結晶組織の観察画像を取得する結晶組織画像取得ステップと、
前記観察画像の結晶粒を識別する結晶粒識別ステップと、
前記結晶粒の形状異方性量を算出する結晶粒形状異方性量算出ステップと、
前記観察画像の中で解析単位領域を設定する解析単位領域設定ステップと、
前記解析単位領域の中で前記結晶粒の部分面積を算出する結晶粒部分面積算出ステップと、
前記形状異方性量と前記部分面積とから前記解析単位領域の平均形状異方性量を算出する平均形状異方性量算出ステップと、
前記平均形状異方性量の分布を作成する平均形状異方性量分布作成ステップと、
前記平均形状異方性量の分布から柱状晶領域と等軸晶領域とを判別するしきい値を設定して柱状晶領域/等軸晶領域を判定する柱状晶領域/等軸晶領域判定ステップと、
を有することを特徴とする金属凝固体結晶組織の解析方法、を提供するものである。
(I) One aspect of the present invention is a method for analyzing the crystal structure of a solidified metal body.
A crystal structure image acquisition step for acquiring an observation image of the crystal structure, and
A crystal grain identification step for identifying crystal grains in the observation image, and
A step of calculating the amount of crystal grain shape anisotropy for calculating the amount of crystal grain shape anisotropy, and
An analysis unit area setting step for setting an analysis unit area in the observation image, and
A crystal grain partial area calculation step for calculating the partial area of the crystal grain in the analysis unit region, and
An average shape anisotropic amount calculation step for calculating the average shape anisotropic amount of the analysis unit region from the shape anisotropic amount and the partial area, and
The step of creating the average shape anisotropic amount distribution and the step of creating the average shape anisotropic amount distribution,
A columnar crystal region / equiaxed crystal region determination step for determining a columnar crystal region / equiaxed crystal region by setting a threshold value for discriminating between a columnar crystal region and an equiaxed crystal region from the distribution of the average shape anisotropic amount.
The present invention provides a method for analyzing a solidified metal crystal structure, which is characterized by having.
(II)本発明の他の一態様は、上記の金属凝固体結晶組織の解析方法を実施する装置であって、
前記金属凝固体結晶組織の前記観察画像を取得する結晶組織観察画像取得機構と、
前記観察画像の前記結晶粒の前記形状異方性量を算出する結晶粒形状異方性量算出機構と、
前記解析単位領域に基づいて前記結晶粒の前記部分面積を算出する結晶粒部分面積算出機構と、
前記形状異方性量と前記部分面積とから前記平均形状異方性量を算出する平均形状異方性量算出機構と、
前記平均形状異方性量の分布を作成する平均形状異方性量分布作成機構と、
前記平均形状異方性量を二分類する平均形状異方性量二値化機構と、
を有することを特徴とする金属凝固体結晶組織の解析装置、を提供するものである。
(II) Another aspect of the present invention is an apparatus for carrying out the above-mentioned method for analyzing a solidified metal crystal structure.
A crystal structure observation image acquisition mechanism for acquiring the observation image of the metal solidified crystal structure,
A crystal grain shape anisotropy calculation mechanism for calculating the shape anisotropy of the crystal grains in the observation image,
A crystal grain partial area calculation mechanism that calculates the partial area of the crystal grain based on the analysis unit region,
An average shape anisotropy calculation mechanism that calculates the average shape anisotropy from the shape anisotropy and the partial area,
An average shape anisotropic amount distribution creation mechanism that creates a distribution of the average shape anisotropic amount, and
An average shape anisotropy binarization mechanism that classifies the average shape anisotropy into two, and
The present invention provides an analyzer for analyzing a solidified metal crystal structure, which is characterized by having.
本発明によれば、金属凝固体の結晶組織における柱状晶領域と等軸晶領域とを定量的に解析する方法、および該方法を実施するための金属凝固体結晶組織の解析装置を提供することができる。 According to the present invention, a method for quantitatively analyzing a columnar crystal region and an equiaxed crystal region in the crystal structure of a metal solidified body, and an apparatus for analyzing the metal solidified body crystal structure for carrying out the method are provided. Can be done.
以下、本発明の実施形態について図面を参照しながら具体的に説明する。本発明は、ここで取り挙げた実施形態に限定されるものではなく、その発明の技術的思想を逸脱しない範囲で公知技術と適宜組み合わせたり公知技術に基づいて改良したりすることが可能である。 Hereinafter, embodiments of the present invention will be specifically described with reference to the drawings. The present invention is not limited to the embodiments taken up here, and can be appropriately combined with a known technique or improved based on the known technique without departing from the technical idea of the invention. ..
[第1実施形態]
図1は、本発明に係る金属凝固体結晶組織の解析方法におけるプロセスの一例を示すフロー図である。図1に示したように、本発明の金属凝固体結晶組織の解析方法は、大別して、結晶粒特徴量解析プロセスと、特徴量分布解析プロセスと、柱状晶−等軸晶遷移解析プロセスとからなる。
[First Embodiment]
FIG. 1 is a flow chart showing an example of a process in the method for analyzing a solidified metal crystal structure according to the present invention. As shown in FIG. 1, the analysis method of the metal solidified crystal structure of the present invention is roughly classified into a crystal grain feature amount analysis process, a feature amount distribution analysis process, and a columnar crystal-equal axis crystal transition analysis process. Become.
より具体的には、結晶粒特徴量解析プロセスは、結晶組織画像取得ステップS101と、結晶粒識別ステップS102と、結晶粒形状異方性量算出ステップS103とを有する。特徴量分布解析プロセスは、解析単位領域設定ステップS201と、結晶粒部分面積算出ステップS202と、平均形状異方性量算出ステップS203と、平均形状異方性量分布作成ステップS204とを有する。柱状晶−等軸晶遷移解析プロセスは、柱状晶領域/等軸晶領域判定ステップS301と、柱状晶領域/等軸晶領域識別出力ステップS302とを有する。以下、各ステップをより詳細に説明する。 More specifically, the crystal grain feature amount analysis process includes a crystal structure image acquisition step S101, a crystal grain identification step S102, and a crystal grain shape anisotropic amount calculation step S103. The feature amount distribution analysis process includes an analysis unit region setting step S201, a crystal grain partial area calculation step S202, an average shape anisotropic amount calculation step S203, and an average shape anisotropic amount distribution creation step S204. The columnar crystal-equal axis crystal transition analysis process includes a columnar crystal region / equiaxed crystal region determination step S301 and a columnar crystal region / equiaxed crystal region identification output step S302. Hereinafter, each step will be described in more detail.
(結晶組織画像取得ステップ)
本ステップS101は、金属凝固体の結晶組織の観察画像を取得するステップである。前述したように、本発明の目的の一つは、金属凝固体の結晶組織における柱状晶領域と等軸晶領域とを定量的に解析することであることから、観察する結晶組織は、金属凝固体中の柱状晶−等軸晶遷移を観察し易いように、金属凝固体の深さ方向を含む断面を選択することが好ましい。
(Crystal structure image acquisition step)
This step S101 is a step of acquiring an observation image of the crystal structure of the solidified metal. As described above, one of the objects of the present invention is to quantitatively analyze the columnar crystal region and the equiaxed crystal region in the crystal structure of the solidified metal, so that the crystal structure to be observed is metal solidified. It is preferable to select a cross section including the depth direction of the solidified metal body so that the transition between columnar crystals and equiaxed crystals in the body can be easily observed.
結晶組織の画像取得方法としては、各結晶粒の形状を明確に観察し撮影できるかぎり特段の限定はなく、従前の観察方法/画像取得方法を利用できる。例えば、走査型電子顕微鏡や偏光顕微鏡を用いて観察した画像を好適に利用できる。 The image acquisition method of the crystal structure is not particularly limited as long as the shape of each crystal grain can be clearly observed and photographed, and the conventional observation method / image acquisition method can be used. For example, an image observed using a scanning electron microscope or a polarizing microscope can be preferably used.
(結晶粒識別ステップ)
本ステップS102は、取得した観察画像の各結晶粒を識別するステップである。各結晶粒の識別方法も特段の限定はなく、従前の識別方法/画像処理方法を利用できる。例えば、走査型電子顕微鏡を用いた後方散乱電子回折(SEM/EBSD)法、偏光顕微鏡観察画像の各ピクセルのRGB値を基にして各結晶粒を識別する方法、材料組織シミュレーション(例えば、フェーズフィールド法、セルオートマトン法、モンテカルロ法)を利用して各結晶粒を識別する方法を好適に利用できる。
(Crystal grain identification step)
This step S102 is a step of identifying each crystal grain of the acquired observation image. The method for identifying each crystal grain is not particularly limited, and the conventional identification method / image processing method can be used. For example, backscattered electron diffraction (SEM / EBSD) method using a scanning electron microscope, a method of identifying each grain based on the RGB value of each pixel of a polarizing microscope observation image, material structure simulation (for example, phase field). A method for identifying each crystal grain by using a method (method, cell automaton method, Monte Carlo method) can be preferably used.
一例を説明すると、SEM/EBSD法では、取得した観察画像の隣接ピクセル間で結晶方位の角度差が所定の範囲内(例えば15°以内)の場合に同一結晶粒と判定することで、各結晶粒を識別することができる。図2は、観察画像においてピクセル毎の結晶方位判定により各結晶粒を識別した例を示す模式図である。図2においては、観察画像内に5個の結晶粒が存在すると判定される。このとき、各結晶粒を識別するためのインデックス(例えば、結晶粒G1〜結晶粒G5)を付与する。 To explain an example, in the SEM / EBSD method, when the angular difference in crystal orientation between adjacent pixels of the acquired observation image is within a predetermined range (for example, within 15 °), it is determined that each crystal is the same crystal grain. Grains can be identified. FIG. 2 is a schematic diagram showing an example in which each crystal grain is identified by determining the crystal orientation for each pixel in the observation image. In FIG. 2, it is determined that five crystal grains are present in the observation image. At this time, an index for identifying each crystal grain (for example, crystal grain G 1 to crystal grain G 5 ) is given.
(結晶粒形状異方性量算出ステップ)
本ステップS103は、識別された各結晶粒の形状異方性量を算出するステップである。図3は、識別された結晶粒において、該結晶粒の外形を楕円近似して形状異方性量を算出するプロセスの例を示す模式図である。
(Crystal grain shape anisotropic amount calculation step)
This step S103 is a step of calculating the amount of shape anisotropy of each of the identified crystal grains. FIG. 3 is a schematic diagram showing an example of a process of calculating the amount of shape anisotropy by elliptical approximation of the outer shape of the identified crystal grain.
図3に示したように、観察画像内に結晶粒Gi(iは自然数)が識別された場合、観察画像に対してx-y直交座標系を当てはめると、結晶粒Giの外形は、結晶粒Giを定義するピクセルの集合の外周上の各頂点(座標(xj, yj)とする)の集合で表される。 As shown in FIG. 3, when crystal grains G i (i is a natural number) are identified in the observation image, when the xy Cartesian coordinate system is applied to the observation image, the outer shape of the crystal grains G i is the crystal grains. It is represented by a set of each vertex (coordinates (x j , y j )) on the outer circumference of the set of pixels that define G i.
つぎに、当該結晶粒Giの外形を、次の式(1)で規定される楕円に近似する。なお、式(1)中の係数Ai、Bi、Ci、Di、Eiは、式(2)で規定される誤差Error(i)が最小になるように選定する。 Next, the outer shape of the crystal grain G i is approximated to an ellipse defined by the following equation (1). The coefficients A i , B i , C i , D i , and E i in Eq. (1) are selected so that the error Error (i) specified in Eq. (2) is minimized.
つぎに、式(1)で規定される近似楕円から長軸径aiと短軸径biとを求め、それらの比を結晶粒Giの形状異方性量piとして式(3)で定義する。式(3)の定義から、形状異方性量piの値が小さい方が、形状異方性が大きいことを意味する。 Next, the semimajor axis a i and the minor axis diameter b i are obtained from the approximate ellipse defined by the equation (1), and their ratio is defined by the equation (3) as the shape anisotropic amount p i of the crystal grain G i. To do. From the definition of equation (3), the smaller the value of the shape anisotropy amount p i , the larger the shape anisotropy.
上記プロセスを観察画像内で識別された各結晶粒Giに対して行う。 The above process is performed for each grain G i identified in the observation image.
なお、上記では、結晶粒Giの形状異方性量piを近似楕円の長軸径aiと短軸径biとの比から求めたが、本発明はそれに限定されるものではない。結晶粒Giの形状異方性量piを適切に定義できるかぎり他の手段でもよい。例えば、結晶粒Giの外形に対する外接楕円または内接楕円を定め、その長軸径短軸径比から形状異方性量piを求めてもよい。また、結晶粒Giの外形に対する外接矩形または内接矩形を定め、その長辺短辺比から形状異方性量piを求めてもよい。 In the above, it was determined from the ratio of the crystal grains G major axis diameter of the shape anisotropy amount p i the approximation ellipses i a i and minor axis b i, the present invention is not limited thereto. Shape anisotropy amount p i of grain G i may be any other means as long as can properly define. For example, an circumscribed ellipse or an inscribed ellipse with respect to the outer shape of the crystal grain G i may be determined, and the shape anisotropic amount p i may be obtained from the semi-major axis-minor axis diameter ratio thereof. Further, the circumscribed rectangle or the inscribed rectangle with respect to the outer shape of the crystal grain G i may be determined, and the shape anisotropic amount p i may be obtained from the long side short side ratio thereof.
(解析単位領域設定ステップ)
本ステップS201は、観察画像内で解析単位領域を設定するステップである。解析単位領域Rn(nは自然数)とは、観察画像全体の結晶組織を解析するための基本単位であり、観察画像全体を網羅するように2以上の領域に分割したものである。また、解析単位領域Rnの最小面積は、観察画像内の最小結晶粒以上の面積となるように設定するのが好ましい。
(Analysis unit area setting step)
This step S201 is a step of setting the analysis unit area in the observed image. The analysis unit region R n (n is a natural number) is a basic unit for analyzing the crystal structure of the entire observation image, and is divided into two or more regions so as to cover the entire observation image. Further, it is preferable that the minimum area of the analysis unit region R n is set to be an area equal to or larger than the minimum crystal grain in the observed image.
なお、分割する方向としては、結晶組織における柱状晶−等軸晶遷移を解析し易くするために、少なくとも金属凝固体の深さ方向を分割するように解析単位領域Rnを設定することが好ましい。解析単位領域Rnの形状に特段の限定はないが、次工程以降の計算容易性の観点からは同一形状に分割することが好ましく、形状を数式で定義し易い矩形状がより好ましい。 As the division direction, it is preferable to set the analysis unit region R n so as to divide at least the depth direction of the solidified metal body in order to facilitate the analysis of the columnar crystal-isoaxial crystal transition in the crystal structure. .. The shape of the analysis unit region R n is not particularly limited, but from the viewpoint of ease of calculation in the next step and thereafter, it is preferable to divide the shape into the same shape, and a rectangular shape in which the shape can be easily defined by a mathematical formula is more preferable.
本ステップS201は、次ステップS202の前であればよく、ステップS103の後である必然性はない。例えば、ステップS101やステップS102の直後に行ってもよい。 This step S201 may be before the next step S202, and does not necessarily have to be after step S103. For example, it may be performed immediately after step S101 or step S102.
(結晶粒部分面積算出ステップ)
本ステップS202は、各解析単位領域Rn内で各結晶粒Giの部分面積を算出するステップである。図4は、取得した観察画像に対して、解析単位領域設定ステップおよび結晶粒部分面積算出ステップを行う際のプロセス例を示す模式図である。図4では、金属凝固体の表面に平行方向に9 mmで金属凝固体の深さ方向に3 mmの観察画像を用い、当該観察画像の左下角を原点とするx-y直交座標系を当てはめた。
(Crystal grain partial area calculation step)
This step S202 is a step of calculating the partial area of each crystal grain G i within each analysis unit region R n. FIG. 4 is a schematic diagram showing a process example when performing the analysis unit region setting step and the crystal grain partial area calculation step on the acquired observation image. In FIG. 4, an observation image of 9 mm parallel to the surface of the solid metal body and 3 mm in the depth direction of the solid metal body was used, and an xy Cartesian coordinate system with the lower left corner of the observation image as the origin was applied.
また、図4では、解析単位領域設定ステップS201として、幅(x軸方向)が観察画像の幅と同じであり高さ(y軸方向)が観察画像の高さを10分割した矩形状の解析単位領域Rn(n=1〜10)を設定した例を示した。結晶粒部分面積算出ステップS202においては、結晶粒形状異方性量算出ステップS103で導出した結晶粒Giの近似楕円(式(1)参照)が、解析単位領域設定ステップS201で設定した解析単位領域Rnの中で占める部分面積Ai-nを算出する。 Further, in FIG. 4, as the analysis unit area setting step S201, a rectangular analysis in which the width (x-axis direction) is the same as the width of the observed image and the height (y-axis direction) is the height of the observed image divided into 10. An example in which the unit area R n (n = 1 to 10) is set is shown. In the crystal grain partial area calculation step S202, the approximate ellipse of the crystal grain G i derived in the crystal grain shape anisotropic amount calculation step S103 (see equation (1)) is the analysis unit region R set in the analysis unit region setting step S201. Calculate the partial area A in that occupies in n.
上記プロセスを各解析単位領域Rnおよび各結晶粒Giに対して行う。 The above process is performed for each analysis unit region R n and each grain G i .
(平均形状異方性量算出ステップ)
本ステップS203は、解析単位領域Rn内で各結晶粒Giの形状異方性量piと部分面積Ai-nとから形状異方性量の面積荷重平均(解析単位領域Rnの平均形状異方性量pn(ave))を算出するステップである。平均形状異方性量pn(ave)の算出を数式で表すと、次の式(4)のようになる。
(Step for calculating average shape anisotropic amount)
This step S203, the average shape anisotropy amount p n of the analysis unit regions R n within an area weighted average of the shape anisotropy amount from the shape anisotropy amount p i and the partial area A in each grain G i (analysis unit region R n (ave)) ) is a step to calculate. The calculation of the average shape anisotropic amount p n (ave) can be expressed by a mathematical formula as shown in the following equation (4).
当該プロセスを各解析単位領域Rnに対して行う。 The process is performed for each analysis unit region R n .
(平均形状異方性量分布作成ステップ)
本ステップS204は、前ステップS203で算出した解析単位領域Rnの平均形状異方性量pn(ave)の分布を作成するステップである。図5は、図4の観察画像に対して第1実施形態の平均形状異方性量分布作成ステップを行った結果例を示す平均形状異方性量の分布グラフである。ここで、平均形状異方性量pn(ave)は、解析単位領域Rn全体の特徴量であることから、解析単位領域Rnの重心位置に代表点があると考えることができる。よって、図5において、縦軸のプロット位置(y座標)は各解析単位領域Rnの高さの中央位置とした。
(Step for creating average shape anisotropic amount distribution)
This step S204 is a step of creating a distribution of the average shape anisotropy amount p n (ave) of the analysis unit region R n calculated in the previous step S203. FIG. 5 is a distribution graph of the average shape anisotropic amount showing an example of the result of performing the average shape anisotropic amount distribution creation step of the first embodiment on the observation image of FIG. Here, the average shape anisotropy weight p n (ave) can be considered because it is a feature of the entire analysis unit region R n, and there is a representative point in the center of gravity of the analysis unit region R n. Therefore, in FIG. 5, the plot position (y coordinate) on the vertical axis is set to the center position of the height of each analysis unit region R n.
図5に示したように、「y=0.15〜2.25」においては「pn(ave)<0.25」であるが、「y=2.55〜2.85」では「pn(ave)>0.4」になっており、「y=2.25〜2.55」の間で急激に変化していることが判る。 As shown in FIG. 5, "p n (ave) <0.25" in "y = 0.15 to 2.25", but "p n (ave) >0.4" in "y = 2.55 to 2.85". It can be seen that there is a rapid change between "y = 2.25 and 2.55".
(柱状晶領域/等軸晶領域判定ステップ)
本ステップS301は、平均形状異方性量の分布から柱状晶領域と等軸晶領域とを判別するしきい値を設定して、該しきい値に基づいて柱状晶領域/等軸晶領域を定量的に判定するステップである。図5の例では、前述したように、「y=2.25〜2.55」の間で平均形状異方性量pn(ave)に急激な変化が生じている(隣接する平均形状異方性量pn(ave)の間の差が最も大きい)。そこで、例えば「y=2.25」と「y=2.55」との中点である「y=2.4」にしきい値を設定することができる。
(Columnar crystal region / equiaxed crystal region determination step)
In this step S301, a threshold value for discriminating between the columnar crystal region and the equiaxed crystal region is set from the distribution of the average shape anisotropic amount, and the columnar crystal region / equiaxed crystal region is quantitatively determined based on the threshold value. It is a step to judge. In the example of FIG. 5, as described above, the average shape anisotropy amount p n (ave) suddenly changes between “y = 2.25 to 2.55” (adjacent average shape anisotropy amount p n (ave)). The difference between them is the largest). Therefore, for example, a threshold value can be set at "y = 2.4", which is the midpoint between "y = 2.25" and "y = 2.55".
図5の平均形状異方性量の分布グラフと図4の観察画像とを比較すると、「y=2.4」のしきい値は、柱状晶領域と等軸晶領域との境界を良く表していると言える。言い換えると、ステップS204で平均形状異方性量pn(ave)の分布を作成し、ステップS301で柱状晶領域/等軸晶領域の判別しきい値を設定することによって、観察画像全体で柱状晶領域/等軸晶領域を定量的に判定できることが判る。 Comparing the distribution graph of the average shape anisotropic amount in FIG. 5 with the observation image in FIG. 4, it can be said that the threshold value of “y = 2.4” well represents the boundary between the columnar crystal region and the equiaxed crystal region. .. In other words, by creating the distribution of the average shape anisotropic amount p n (ave) in step S204 and setting the discrimination threshold value of the columnar crystal region / equiaxed crystal region in step S301, the columnar crystal region in the entire observation image. / It can be seen that the equiaxed crystal region can be quantitatively determined.
(柱状晶領域/等軸晶領域出力ステップ)
本ステップS302は、前ステップS301で得られた柱状晶領域/等軸晶領域の判定結果を出力するステップである。本発明において、本ステップS302で行う出力手段に特段の限定はなく、利用者が望む出力形式/出力方法(例えば、ディスプレー画面表示、紙面印刷、データ出力)を適宜利用できる。
(Columnar crystal region / equiaxed crystal region output step)
This step S302 is a step of outputting the determination result of the columnar crystal region / equiaxed crystal region obtained in the previous step S301. In the present invention, the output means performed in step S302 is not particularly limited, and the output format / output method desired by the user (for example, display screen display, paper printing, data output) can be appropriately used.
第1実施形態の金属凝固体結晶組織の解析方法は、解析単位領域設定ステップS201において、金属凝固体の深さ方向のみを分割するように解析単位領域Rnを設定したことから(図4では、金属凝固体の深さ方向が高さ方向となる観察画像に対して、観察画像の幅と同じ幅で観察画像の高さ方向を分割した高さを有する解析単位領域Rnを設定したことから)、金属凝固体の深さ方向の全体傾向を把握する(一次元的に定量判定する)のに適している。また、解析単位領域Rnの分割総数が比較的少ないことから、計算回数を少なくすることができる(計算時間を短くできる/短時間で計算結果が得られる)という利点もある。 In the analysis method of the metal solidified body crystal structure of the first embodiment, since the analysis unit region R n was set so as to divide only the depth direction of the metal solidified body in the analysis unit region setting step S201 (in FIG. 4). , The analysis unit region R n having the height obtained by dividing the height direction of the observation image by the same width as the width of the observation image was set for the observation image in which the depth direction of the metal solidified body is the height direction. It is suitable for grasping the overall tendency of the solidified metal in the depth direction (one-dimensional quantitative determination). In addition, since the total number of divisions of the analysis unit region R n is relatively small, there is an advantage that the number of calculations can be reduced (calculation time can be shortened / calculation results can be obtained in a short time).
[第2実施形態]
第2実施形態の金属凝固体結晶組織の解析方法は、解析単位領域設定ステップS201において、解析単位領域を第1実施形態よりも細かく/小さく設定して、観察画像全体を二次元的に定量判定する点で異なる。以下、第1実施形態と異なる点を重点的に説明する。
[Second Embodiment]
In the method for analyzing the metal solidified crystal structure of the second embodiment, in the analysis unit region setting step S201, the analysis unit region is set finer / smaller than that of the first embodiment, and the entire observed image is two-dimensionally quantitatively determined. It differs in that it does. Hereinafter, the points different from the first embodiment will be mainly described.
まず、結晶組織画像取得ステップS101〜結晶粒形状異方性量算出ステップS103は、第1実施形態と同じである。 First, the crystal structure image acquisition step S101 to the crystal grain shape anisotropic amount calculation step S103 are the same as those in the first embodiment.
一方、解析単位領域設定ステップS201において、上述したように、解析単位領域を第1実施形態よりも細かく設定する。例えば、図4の観察画像に対して、y軸方向で10分割しx軸方向で30分割した矩形状の解析単位領域Rn,m(n=1〜10,m=1〜30)を設定する。 On the other hand, in the analysis unit area setting step S201, as described above, the analysis unit area is set more finely than in the first embodiment. For example, for the observation image of FIG. 4, a rectangular analysis unit region R n, m (n = 1 to 10, m = 1 to 30) divided into 10 in the y-axis direction and 30 in the x-axis direction is set. To do.
つぎに、第1実施形態と同様にして、結晶粒部分面積算出ステップS202において、結晶粒形状異方性量算出ステップS103で導出した結晶粒Giの近似楕円(式(1)参照)が、解析単位領域設定ステップS201で設定した解析単位領域Rn,mの中で占める部分面積Ai-n,mを算出する。 Next, in the same manner as in the first embodiment, in the crystal grain partial area calculation step S202, the approximate ellipse of the crystal grain G i derived in the crystal grain shape anisotropic amount calculation step S103 (see equation (1)) is the analysis unit. The partial area A in, m occupied in the analysis unit area R n, m set in the area setting step S201 is calculated.
つぎに、平均形状異方性量分布作成ステップS204を行う。解析単位領域Rn,mの平均形状異方性量pn,m(ave)の算出としては、式(4)における「Ai-n」を「Ai-n,m」と読み替えればよい。また、第2実施形態においては、解析単位領域Rn,mを比較的細かく設定して観察画像全体を二次元的に評価することから、算出される平均形状異方性量pn,m(ave)に対して色および/または濃淡を当てはめて表記すると見やすい。 Next, step S204 for creating the average shape anisotropic amount distribution is performed. To calculate the average shape anisotropic amount p n, m (ave) in the analysis unit region R n, m , “A in ” in Eq. (4) may be read as “A in, m”. Further, in the second embodiment, the average shape anisotropic amount p n, m (ave) calculated from the two-dimensional evaluation of the entire observed image by setting the analysis unit regions R n, m relatively finely. It is easy to see if the color and / or shade is applied to the notation.
図6は、図4の観察画像に対して第2実施形態の平均形状異方性量分布作成ステップを行った結果例を示す平均形状異方性量の分布マップである。図6に示したように、平均形状異方性量分布マップは、図5の平均形状異方性量分布グラフに比して、観察画像全体の様子をより理解し易い利点がある。例えば、図6において、上から2段目までの方が3段目以降よりも濃色であり、平均形状異方性量pn,m(ave)が大きい(すなわち、形状異方性が小さい)ことが読み取れる。 FIG. 6 is a distribution map of the average shape anisotropic amount showing an example of the result of performing the average shape anisotropic amount distribution creation step of the second embodiment on the observation image of FIG. As shown in FIG. 6, the average shape anisotropic amount distribution map has an advantage that it is easier to understand the state of the entire observed image as compared with the average shape anisotropic amount distribution graph of FIG. For example, in FIG. 6, the color from the top to the second stage is darker than that of the third and subsequent stages, and the average shape anisotropy amount p n, m (ave) is large (that is, the shape anisotropy is small). Can be read.
なお、平均形状異方性量分布マップの左上角の「×」は欠損領域を示す。図4から判るように、当該領域は、観察画像の中で結晶粒Giが存在しない領域であることから、領域内で結晶粒Giの部分面積Ai-n,mがゼロとなり、式(4)で平均形状異方性量pn,m(ave)を算出できない。言い換えると、本発明では、式(4)で平均形状異方性量を算出できない領域を欠損領域として取り扱うものとする。 The “x” in the upper left corner of the average shape anisotropic amount distribution map indicates a defective region. As can be seen from FIG. 4, since this region is a region in which the crystal grains G i do not exist in the observation image, the partial areas A in and m of the crystal grains G i become zero in the region, and the equation (4) ) Cannot calculate the average shape anisotropy amount p n, m (ave). In other words, in the present invention, a region in which the average shape anisotropic amount cannot be calculated by the formula (4) is treated as a defective region.
図6の平均形状異方性量分布マップは、図5の平均形状異方性量分布グラフに比して、観察画像全体の様子をより理解し易い利点があるが、柱状晶領域/等軸晶領域の判別が曖昧になり易い(観測者の主観に依存する)という点に改善の余地がある。そこで、第2実施形態の柱状晶領域/等軸晶領域判定ステップS301では、k平均法によって前ステップS204で算出した平均形状異方性量pn,m(ave)を2つに分類することで、柱状晶領域/等軸晶領域のしきい値を設定する。 The average shape anisotropic amount distribution map of FIG. 6 has an advantage that the state of the entire observed image is easier to understand than the average shape anisotropic amount distribution graph of FIG. 5, but the columnar crystal region / equiaxed crystal region can be discriminated. There is room for improvement in that it tends to be ambiguous (depending on the subjectivity of the observer). Therefore, in the columnar crystal region / equiaxed crystal region determination step S301 of the second embodiment, the average shape anisotropic amount p n, m (ave) calculated in the previous step S204 by the k-means method is classified into two. Set the threshold value for the columnar / equiaxed region.
k平均法とは、多数のデータをk個のクラスタに分類する手法であり、一般的には次のような流れで行われる。(a)各データに対してランダムにクラスタを割り振る。(b)割り振ったデータを基にして各クラスタの重心を計算する。クラスタ重心の計算は、しばしば割り当てられたデータの各要素の算術平均が使用される。(c)各データと各クラスタ重心との距離を求め、各データを最も近い重心のクラスタに割り当て直す。(d)上記の処理で全てのデータクラスタ割り当てが変化しなかった場合、あるいは変化量が事前に設定した一定のしきい値を下回った場合に、収束したと判断して処理を終了する。そうでない場合は新しく割り振られたクラスタから重心を再計算して上記の処理を繰り返す。 The k-means method is a method of classifying a large number of data into k clusters, and is generally performed in the following flow. (A) Randomly allocate clusters for each data. (B) Calculate the center of gravity of each cluster based on the allocated data. The calculation of the center of gravity of the cluster often uses the arithmetic mean of each element of the assigned data. (C) Obtain the distance between each data and the center of gravity of each cluster, and reassign each data to the cluster with the closest center of gravity. (D) When all the data cluster allocations do not change in the above process, or when the amount of change falls below a predetermined threshold value, it is determined that the data clusters have converged and the process ends. If not, the center of gravity is recalculated from the newly allocated cluster and the above process is repeated.
本実施形態では、算出した各平均形状異方性量pn,m(ave)を2個のクラスタ(ここでは、柱状晶領域、等軸晶領域)に分類することになる。k平均法の上記(a)に関しては、例えば、コンピータの乱数生成を利用することができる。上記(d)のしきい値に関しては、例えば、1×10-3〜1×10-6を適宜利用すればよい。また、上記(a)〜(d)の一連を複数回行って、(d)の変化量が最も小さい値を示した結果を採用してもよい。クラスタ重心の値が大きい方が、より小さい形状異方性を意味することから等軸晶領域と判定され、クラスタ重心の値が小さい方が、より大きい形状異方性を意味することから柱状晶領域と判定される。 In the present embodiment, each calculated average shape anisotropic amount p n, m (ave) is classified into two clusters (here, a columnar crystal region and an equiaxed crystal region). Regarding the above (a) of the k-means method, for example, random number generation of a computer can be used. Regarding the threshold value of (d) above, for example, 1 × 10 -3 to 1 × 10 -6 may be appropriately used. Further, the series of the above steps (a) to (d) may be performed a plurality of times, and the result showing the value in which the amount of change in (d) is the smallest may be adopted. A larger value of the cluster center of gravity means a smaller shape anisotropy, so it is determined to be an equiaxed crystal region, and a smaller value of the cluster center of gravity means a larger shape anisotropy. It is determined to be an area.
図7は、算出した平均形状異方性量pn,m(ave)をk平均法によって二値化した結果を示す柱状晶領域/等軸晶領域マップである。k平均法の条件としては、(d)のしきい値を1×10-5とし、(a)〜(d)の一連を5回行って(d)の変化量が最も小さい値を示した結果を採用した。 FIG. 7 is a columnar crystal region / equiaxed crystal region map showing the result of binarizing the calculated average shape anisotropic amount p n, m (ave) by the k-means method. As the conditions of the k-means method, the threshold value of (d) was set to 1 × 10 -5, and the series of (a) to (d) was repeated 5 times to show the value in which the amount of change in (d) was the smallest. The result was adopted.
図7に示したように、柱状晶領域/等軸晶領域マップは、観察画像全体の様子を理解できる利点に加えて、図6の平均形状異方性量分布マップに比して、柱状晶領域/等軸晶領域をより明確に判別できる利点がある。 As shown in FIG. 7, the columnar crystal region / equiaxed crystal region map has the advantage that the state of the entire observed image can be understood, and in addition, the columnar crystal region / equiaxed crystal region map is compared with the average shape anisotropic amount distribution map of FIG. There is an advantage that the equiaxed crystal region can be discriminated more clearly.
なお、上記の説明では、k平均法によって算出した各平均形状異方性量pn,m(ave)を二分類したが、各平均形状異方性量pn,m(ave)を適切に二分類できるかぎり、該k平均法に限定されるものではない。例えば、スペクトラル クラスタリング法、ウォード階層的クラスタリング法などを適宜利用することができる。 In the above explanation, each average shape anisotropic amount p n, m (ave) calculated by the k-means method is classified into two, but as long as each average shape anisotropic amount p n, m (ave) can be appropriately classified into two categories. , The k-means method is not limited. For example, a spectral clustering method, a ward hierarchical clustering method, or the like can be appropriately used.
柱状晶領域/等軸晶領域判定ステップS301の後、第1実施形態と同様に、柱状晶領域/等軸晶領域出力ステップS302として、得られた柱状晶領域/等軸晶領域の判定結果を利用者が望む出力形式/出力方法で出力すればよい。 After the columnar crystal region / equiaxed crystal region determination step S301, as in the first embodiment, as the columnar crystal region / equiaxed crystal region output step S302, the obtained columnar crystal region / equiaxed crystal region determination result is obtained. The output format / output method desired by the user may be used for output.
[第3実施形態]
第3実施形態の金属凝固体結晶組織の解析方法は、平均形状異方性量算出ステップS203および平均形状異方性量分布作成ステップS204において、第1実施形態や第2実施形態のように解析単位領域の平均形状異方性量をスカラーとして算出・表示するのではなく、ベクトルとして算出・表示する点で異なる。以下、第2実施形態と異なる点を重点的に説明する。
[Third Embodiment]
The method for analyzing the metal solidified crystal structure of the third embodiment is the average of the analysis unit regions as in the first embodiment and the second embodiment in the average shape anisotropic amount calculation step S203 and the average shape anisotropic amount distribution creation step S204. The difference is that the amount of shape anisotropy is not calculated and displayed as a scalar, but is calculated and displayed as a vector. Hereinafter, the points different from the second embodiment will be mainly described.
まず、結晶組織画像取得ステップS101〜結晶粒部分面積算出ステップS202は、第2実施形態と同様に行えばよい。 First, the crystal structure image acquisition step S101 to the crystal grain partial area calculation step S202 may be performed in the same manner as in the second embodiment.
一方、平均形状異方性量算出ステップS203においては、上述したように、解析単位領域Rn,mの平均形状異方性量をスカラーとして算出するのではなく、ベクトル(平均形状異方性ベクトルpvec n,m(ave))として算出する。ベクトルを定義するためには、大きさと方向とを規定する必要がある。 On the other hand, in the average shape anisotropy calculation step S203, as described above, the average shape anisotropy of the analysis unit regions R n and m is not calculated as a scalar, but a vector (average shape anisotropy vector p vec n, Calculate as m (ave)). In order to define a vector, it is necessary to specify the magnitude and direction.
平均形状異方性ベクトルpvec n,m(ave)の算出方法としては、平均形状異方性量pn,m(ave)の算出方法と同様の考え方で、解析単位領域Rn,mを構成する各結晶粒Giの近似楕円の長軸方向を面積荷重合成する(各結晶粒Giの部分面積Ai-n,mの比率の長さを有するベクトルを合成する)方法でベクトルの方向を導出することができる。平均形状異方性ベクトルpvec n,m(ave)の大きさ(ノルム)としては、形状異方性が大きい方が大きな値を有することが望ましいので、例えば、第2実施形態の平均形状異方性量pn,m(ave)を逆数化した1/pn,m(ave)(解析単位領域Rn,mの平均形状アスペクト比に相当)を採用することができる。 As the calculation method of the average shape anisotropy vector p vec n, m (ave) , the analysis unit region R n, m is constructed in the same way as the calculation method of the average shape anisotropy amount p n, m (ave). deriving a direction of the vector in each crystal grain G i approximation surface load synthesize major axis of the ellipse of (partial area a in each grain G i, to synthesize a vector having a length ratio of m) mETHOD be able to. As the magnitude (norm) of the average shape anisotropy vector p vec n, m (ave) , it is desirable that the larger the shape anisotropy, the larger the value. It is possible to adopt 1 / p n, m (ave) (corresponding to the average shape aspect ratio of the analysis unit region R n, m ), which is the reciprocal of the square quantities p n, m (ave).
また、ベクトル合成に際しては、ベクトル方向の正負を定める必要があるが、例えば、金属凝固体が凝固する際の抜熱方向(図4においてはy軸の負方向)をベクトルの正方向としたり、金属凝固体の凝固方向(図4においてはy軸の正方向)をベクトルの正方向としたりして揃えればよい。なお、金属凝固体が凝固する際の抜熱方向と凝固方向とは、通常逆方向になる。 Further, in vector synthesis, it is necessary to determine the positive and negative directions of the vector. For example, the heat removal direction (negative direction of the y-axis in FIG. 4) when the solidified metal solidifies is set to the positive direction of the vector. The solidification direction of the solidified metal (the positive direction of the y-axis in FIG. 4) may be aligned with the positive direction of the vector. The heat removal direction and the solidification direction when the metal solidified body is solidified are usually opposite to each other.
つぎに、平均形状異方性量分布作成ステップS204を行う。図8は、図4の観察画像に対して第3実施形態の平均形状異方性量分布作成ステップを行った結果例を示す平均形状異方性ベクトルの分布マップである。本例では、前ステップS203の段階で、解析単位領域Rn,mを、図6と同じくy軸方向で10分割しx軸方向で30分割した矩形状の領域とした。平均形状異方性ベクトルpvec n,m(ave)の方向を算出するためのベクトル合成においては、金属凝固体が凝固する際の抜熱方向(y軸の負方向)をベクトルの正方向にして揃えた。 Next, step S204 for creating the average shape anisotropic amount distribution is performed. FIG. 8 is a distribution map of the average shape anisotropy vector showing an example of the result of performing the average shape anisotropy distribution creation step of the third embodiment on the observation image of FIG. In this example, at the stage of the previous step S203, the analysis unit areas R n and m are divided into 10 in the y-axis direction and 30 in the x-axis direction as in FIG. In vector synthesis for calculating the direction of the average shape anisotropy vector p vec n, m (ave) , the heat removal direction (negative direction of the y-axis) when the metal solidified body solidifies is set to the positive direction of the vector. I prepared it.
また、第1実施形態と同様の考え方により、平均形状異方性ベクトルpvec n,m(ave)は、解析単位領域Rn,m全体の特徴量であることから、解析単位領域Rn,mの重心位置に代表点があると考えることができる。よって、図7において、各ベクトルの始点は、各解析単位領域の重心位置とした。 Further, based on the same concept as in the first embodiment, since the average shape anisotropy vector p vec n, m (ave) is a feature quantity of the entire analysis unit region R n, m , the analysis unit region R n, It can be considered that there is a representative point at the position of the center of gravity of m. Therefore, in FIG. 7, the starting point of each vector is the position of the center of gravity of each analysis unit region.
図8に示したように、平均形状異方性ベクトル分布マップは、図6の平均形状異方性量分布マップの情報に加えて、結晶粒の異方性の様子(金属凝固体が凝固する際の抜熱方向)をより理解し易い利点がある。例えば、図8において、上から2段目までの方が3段目以降よりもベクトル長さが短いことから、形状異方性が小さいことが読み取れる。また、ベクトルの向きおよび集中具合から、金属凝固体が凝固する際の熱の流れを読み取ることができる。 As shown in FIG. 8, in the average shape anisotropy vector distribution map, in addition to the information of the average shape anisotropy distribution map of FIG. 6, the state of anisotropy of the crystal grains (when the metal solidified body solidifies). There is an advantage that it is easier to understand the heat removal direction). For example, in FIG. 8, since the vector length from the top to the second stage is shorter than that from the third stage onward, it can be read that the shape anisotropy is small. In addition, the heat flow when the metal solidified body solidifies can be read from the direction and concentration of the vector.
ただし、図8の平均形状異方性ベクトル分布マップは、図6と同様の点(柱状晶領域/等軸晶領域の判別が観測者の主観に依存するという点)に改善の余地を有する。そこで、本実施形態においても、平均形状異方性量分布作成ステップS204のつぎに、第2実施形態と同様の柱状晶領域/等軸晶領域判定ステップS301を行うことは好ましい。それにより、図7と図8とを合体したような柱状晶領域/等軸晶領域マップが得られる(図示は省略する)。 However, the average shape anisotropy vector distribution map of FIG. 8 has room for improvement in the same points as in FIG. 6 (the point that the discrimination between the columnar crystal region / equiaxed crystal region depends on the subjectivity of the observer). Therefore, also in this embodiment, it is preferable to perform the same columnar crystal region / equiaxed crystal region determination step S301 as in the second embodiment after the average shape anisotropic amount distribution creation step S204. As a result, a columnar crystal region / equiaxed crystal region map as if FIG. 7 and FIG. 8 are combined can be obtained (not shown).
柱状晶領域/等軸晶領域判定ステップS301の後、第2実施形態と同様に、柱状晶領域/等軸晶領域出力ステップS302として、得られた柱状晶領域/等軸晶領域マップを利用者が望む出力形式/出力方法で出力すればよい。 After the columnar crystal region / equiaxed crystal region determination step S301, the user uses the obtained columnar crystal region / equiaxed crystal region map as the columnar crystal region / equiaxed crystal region output step S302 as in the second embodiment. You can output in the output format / output method you desire.
なお、上記では、第2実施形態をベースにして平均形状異方性量算出ステップS203〜柱状晶領域/等軸晶領域判定ステップS301を説明したが、本実施形態は、それに限定されるものではなく、第1実施形態をベースにしてステップS203〜ステップS301を行ってもよい。 In the above description, the average shape anisotropic amount calculation step S203 to the columnar crystal region / equiaxed crystal region determination step S301 have been described based on the second embodiment, but the present embodiment is not limited thereto. Steps S203 to S301 may be performed based on the first embodiment.
[第4実施形態]
本実施形態では、第1〜第3実施形態の方法を実施するための金属凝固体結晶組織の解析装置について説明する。図9は、本発明に係る金属凝固体結晶組織の解析装置の概略構成を示す模式図である。
[Fourth Embodiment]
In this embodiment, an analyzer for analyzing a solidified metal crystal structure for carrying out the methods of the first to third embodiments will be described. FIG. 9 is a schematic view showing a schematic configuration of an analyzer for analyzing a solidified metal crystal structure according to the present invention.
金属凝固体結晶組織の解析装置100は、金属凝固体の柱状晶−等軸晶遷移による柱状晶領域と等軸晶領域とを定量的に解析する装置であって、大きく分けて、金属凝固体の結晶組織の観察画像を取得する観察画像データ取得部10と、観察画像データから結晶粒の特徴量解析、結晶組織の特徴量分布解析、および柱状晶−等軸晶遷移の定量的解析を行う観察画像データ解析部20と、観察画像データおよび各種解析データを記憶するデータ記憶部30と、解析条件の入力および解析結果の出力を行う入出力部40と、から構成される。以下、各部の内部構成をより具体的に説明する。
The metal solidified body crystal
画像データ取得部10は、金属凝固体結晶組織の観察画像を取り込む結晶組織観察画像取得機構11を有し、結晶組織画像取得ステップS101を担う。結晶組織観察画像取得機構11は、結晶組織観察自体を行う顕微鏡類(顕微鏡、該顕微鏡の制御装置、顕微鏡観察データの解析装置など)を含んでいてもよいし、別装置の顕微鏡類で観察した結晶組織の画像データを取り込む機構であってもよい。 The image data acquisition unit 10 has a crystal structure observation image acquisition mechanism 11 that captures an observation image of the metal solidified crystal structure, and is responsible for the crystal structure image acquisition step S101. The crystal structure observation image acquisition mechanism 11 may include microscopes (microscope, control device for the microscope, analysis device for microscope observation data, etc.) for observing the crystal structure itself, or observed with a microscope of another device. It may be a mechanism for capturing image data of the crystal structure.
取り込んだ観察画像データは、データ記憶部30のデータ記憶機構31において記憶される。データ記憶機構31は、必要なデータが記憶できるかぎり特段の限定はなく、従前のデータ記憶装置(例えば、ランダムアクセスメモリ(RAM)、ハードディスク(HD)、ソリッドステートドライブ(SSD)など)を適宜利用できる。
The captured observation image data is stored in the
観察画像データ解析部20は、取得した観察画像の各結晶粒を識別する結晶粒識別機構21と、識別された各結晶粒の近似形状を算出し該近似形状から形状異方性量を算出する結晶粒形状異方性量算出機構22と、別途設定される解析単位領域に基づいて各結晶粒の部分面積を算出する結晶粒部分面積算出機構23と、各結晶粒の形状異方性量と部分面積とから各解析単位領域の平均形状異方性量を算出する平均形状異方性量算出機構24と、平均形状異方性量の分布を作成する平均形状異方性量分布作成機構25と、各平均形状異方性量を二値化する平均形状異方性量二値化機構26と、各結晶粒の近似形状の長軸方向と部分面積とから各解析単位領域の平均形状異方性ベクトルを算出する平均形状異方性ベクトル算出機構27と、を有する。各機構内で計算された結果は、データ記憶機構31において記憶される。
The observation image
観察画像データ解析部20は、基本的に結晶粒識別ステップS102、結晶粒形状異方性量算出ステップS103、結晶粒部分面積算出ステップS202、平均形状異方性量算出ステップS203、平均形状異方性量分布作成ステップS204、および柱状晶領域/等軸晶領域判定ステップS301を担う。
The observation image
必要とされる解析が適切に行われるかぎり、観察画像データ解析部20の各機構の装置構成に特段の限定はなく、従前の解析装置(例えば、コンピータ)を適宜利用できる。
As long as the required analysis is properly performed, the device configuration of each mechanism of the observation image
なお、画像データ取得部10の結晶組織観察画像取得機構11自体が観察画像の各結晶粒を識別する機能を有する場合、観察画像データ解析部20の結晶粒識別機構21は必須の構成ではない。また、第3実施形態に示したようなベクトル解析を不要とする装置構成の場合、観察画像データ解析部20の平均形状異方性ベクトル算出機構27も必須の構成ではない。
When the crystal structure observation image acquisition mechanism 11 of the image data acquisition unit 10 itself has a function of identifying each crystal grain of the observation image, the crystal
入出力部40は、解析条件(例えば、解析単位領域Rnのサイズ、柱状晶領域/等軸晶領域の識別しきい値の条件)の入力を行う入力機構41と、解析結果の出力を行う出力機構42とを有し、解析単位領域設定ステップS201、および柱状晶領域/等軸晶領域出力ステップS302を担う。入力された各種解析条件、および柱状晶領域/等軸晶領域判定の結果は、データ記憶機構31において記憶される。
The input /
必要とされる入力および望まれる出力ができるかぎり、入力機構41および出力機構42の装置構成に特段の限定はなく、従前の入出力装置(例えば、キーボード、ディスプレー、プリンタ)を適宜利用できる。
As long as the required input and the desired output can be obtained, the device configuration of the
上述した実施形態は、本発明の理解を助けるために説明したものであり、本発明は、記載した具体的な構成のみに限定されるものではない。例えば、実施形態の構成の一部を当業者の技術常識の構成に置き換えることが可能であり、実施形態の構成に当業者の技術常識の構成を加えることも可能である。また、実施形態同士を適宜組み合わせてもよい。すなわち、本発明は、本明細書の実施形態の構成の一部について、発明の技術的思想を逸脱しない範囲で、削除・他の構成による置換・他の構成の追加をすることが可能である。 The above-described embodiments have been described for the purpose of assisting the understanding of the present invention, and the present invention is not limited to the specific configurations described. For example, it is possible to replace a part of the configuration of the embodiment with the configuration of the common general technical knowledge of those skilled in the art, and it is also possible to add the configuration of the common general technical knowledge of the person skilled in the art to the configuration of the embodiment. In addition, the embodiments may be combined as appropriate. That is, the present invention can delete, replace, or add another configuration to a part of the configuration of the embodiment of the present specification without departing from the technical idea of the invention. ..
100・・・金属凝固体結晶組織の解析装置、
10・・・観察画像データ取得部、11・・・結晶組織観察画像取得機構、
20・・・観察画像データ解析部、21・・・結晶粒識別機構、
22・・・結晶粒形状異方性量算出機構、23・・・結晶粒部分面積算出機構、
24・・・平均形状異方性量算出機構、25・・・平均形状異方性量分布作成機構、
26・・・平均形状異方性量二値化機構、27・・・平均形状異方性ベクトル算出機構、
30・・・データ記憶部、31・・・データ記憶機構、
40・・・入出力部、41・・・入力機構、42・・・出力機構。
100 ・ ・ ・ Analysis device for metal solidified crystal structure,
10 ... Observation image data acquisition unit, 11 ... Crystal structure observation image acquisition mechanism,
20 ... Observation image data analysis unit, 21 ... Crystal grain identification mechanism,
22 ... Crystal grain shape anisotropic amount calculation mechanism, 23 ... Crystal grain partial area calculation mechanism,
24 ... Average shape anisotropic amount calculation mechanism, 25 ... Average shape anisotropic amount distribution creation mechanism,
26 ... Average shape anisotropy binarization mechanism, 27 ... Average shape anisotropy vector calculation mechanism,
30 ... Data storage unit, 31 ... Data storage mechanism,
40: Input / output unit, 41: Input mechanism, 42: Output mechanism.
Claims (11)
前記結晶組織の観察画像を取得する結晶組織画像取得ステップと、
前記観察画像の結晶粒を識別する結晶粒識別ステップと、
前記結晶粒の形状異方性量を算出する結晶粒形状異方性量算出ステップと、
前記観察画像の中で解析単位領域を設定する解析単位領域設定ステップと、
前記解析単位領域の中で前記結晶粒の部分面積を算出する結晶粒部分面積算出ステップと、
前記形状異方性量と前記部分面積とから前記解析単位領域の平均形状異方性量を算出する平均形状異方性量算出ステップと、
前記平均形状異方性量の分布を作成する平均形状異方性量分布作成ステップと、
前記平均形状異方性量の分布から柱状晶領域と等軸晶領域とを判別するしきい値を設定して柱状晶領域/等軸晶領域を判定する柱状晶領域/等軸晶領域判定ステップと、
を有することを特徴とする金属凝固体結晶組織の解析方法。 This is a method for analyzing the crystal structure of a solid metal body.
A crystal structure image acquisition step for acquiring an observation image of the crystal structure, and
A crystal grain identification step for identifying crystal grains in the observation image, and
A step of calculating the amount of crystal grain shape anisotropy for calculating the amount of crystal grain shape anisotropy,
An analysis unit area setting step for setting an analysis unit area in the observation image, and
A crystal grain partial area calculation step for calculating the partial area of the crystal grain in the analysis unit region, and
An average shape anisotropic amount calculation step for calculating the average shape anisotropic amount of the analysis unit region from the shape anisotropic amount and the partial area, and
The step of creating the average shape anisotropic amount distribution and the step of creating the average shape anisotropic amount distribution,
A columnar crystal region / equiaxed crystal region determination step for determining a columnar crystal region / equiaxed crystal region by setting a threshold value for discriminating between a columnar crystal region and an equiaxed crystal region from the distribution of the average shape anisotropic amount.
A method for analyzing a solidified metal crystal structure.
前記結晶粒形状異方性量算出ステップは、該結晶粒の外形を楕円に近似して、該楕円の短軸径と長軸径との比から前記形状異方性量を算出するステップであることを特徴とする金属凝固体結晶組織の解析方法。 In the method for analyzing a solidified metal crystal structure according to claim 1,
The crystal grain shape anisotropic amount calculation step is characterized in that the outer shape of the crystal grain is approximated to an ellipse and the shape anisotropic amount is calculated from the ratio of the minor axis diameter to the major axis diameter of the ellipse. A method for analyzing the crystal structure of a solidified metal.
前記観察画像は、前記金属凝固体の深さ方向を含む断面の観察画像であり、
前記深さ方向をx-y直交座標系のy軸方向に当てはめた場合、前記解析単位領域設定ステップは、前記解析単位領域として、前記観察画像のy軸方向長さを複数個に分割した長さと前記観察画像のx軸方向長さとを有する矩形領域を設定するステップであることを特徴とする金属凝固体結晶組織の解析方法。 In the method for analyzing a solidified metal crystal structure according to claim 2.
The observation image is an observation image of a cross section including the depth direction of the solidified metal body.
When the depth direction is applied to the y-axis direction of the xy Cartesian coordinate system, the analysis unit region setting step is the length obtained by dividing the y-axis length of the observation image into a plurality of the analysis unit regions and the above. A method for analyzing a solid metal crystal structure, which is a step of setting a rectangular region having a length in the x-axis direction of an observation image.
前記観察画像は、前記金属凝固体の深さ方向を含む断面の観察画像であり、
前記深さ方向をx-y直交座標系のy軸方向に当てはめた場合、前記解析単位領域設定ステップは、前記解析単位領域として、前記観察画像のy軸方向長さを複数個に分割した長さと前記観察画像のx軸方向長さを複数個に分割した長さとを有する矩形領域を設定するステップであることを特徴とする金属凝固体結晶組織の解析方法。 In the method for analyzing a solidified metal crystal structure according to claim 2.
The observation image is an observation image of a cross section including the depth direction of the solidified metal body.
When the depth direction is applied to the y-axis direction of the xy Cartesian coordinate system, the analysis unit region setting step is the length obtained by dividing the length of the observation image in the y-axis direction into a plurality of the analysis unit regions and the above. A method for analyzing a metal solidified crystal structure, which is a step of setting a rectangular region having a length obtained by dividing the length in the x-axis direction of an observation image into a plurality of parts.
前記柱状晶領域/等軸晶領域解析ステップにおける前記しきい値の設定をk平均法によって行うことを特徴とする金属凝固体結晶組織の解析方法。 In the method for analyzing a solidified metal crystal structure according to claim 4,
A method for analyzing a metal solidified crystal structure, which comprises setting the threshold value in the columnar crystal region / equiaxed crystal region analysis step by the k-means method.
前記平均形状異方性量算出ステップにおける前記平均形状異方性量がベクトルであり、
前記ベクトルは、前記形状異方性量の逆数と前記部分面積とから求められる面積荷重平均の大きさと、前記結晶粒の外形を近似した前記楕円の長軸方向と前記部分面積とから求められる面積荷重平均の方向とを有することを特徴とする金属凝固体結晶組織の解析方法。 In the method for analyzing a solidified metal crystal structure according to claim 4 or 5.
The average shape anisotropic amount in the average shape anisotropic amount calculation step is a vector.
The vector is an area load average obtained from the reciprocal of the shape anisotropic amount and the area load average obtained from the partial area, and the major axis direction of the ellipse and the partial area that approximate the outer shape of the crystal grains. A method for analyzing a metal solidified crystal structure, which comprises the direction of.
前記金属凝固体結晶組織の前記観察画像を取得する結晶組織観察画像取得機構と、
前記観察画像の前記結晶粒の前記形状異方性量を算出する結晶粒形状異方性量算出機構と、
前記解析単位領域に基づいて前記結晶粒の前記部分面積を算出する結晶粒部分面積算出機構と、
前記形状異方性量と前記部分面積とから前記平均形状異方性量を算出する平均形状異方性量算出機構と、
前記平均形状異方性量の分布を作成する平均形状異方性量分布作成機構と、
前記平均形状異方性量を二分類する平均形状異方性量二値化機構と、
を有することを特徴とする金属凝固体結晶組織の解析装置。 An apparatus for carrying out the method for analyzing a solidified metal crystal structure according to any one of claims 1 to 5.
A crystal structure observation image acquisition mechanism for acquiring the observation image of the metal solidified crystal structure,
A crystal grain shape anisotropy calculation mechanism for calculating the shape anisotropy of the crystal grains in the observation image,
A crystal grain partial area calculation mechanism that calculates the partial area of the crystal grain based on the analysis unit region,
An average shape anisotropy calculation mechanism that calculates the average shape anisotropy from the shape anisotropy and the partial area,
An average shape anisotropic amount distribution creation mechanism that creates a distribution of the average shape anisotropic amount, and
An average shape anisotropy binarization mechanism that classifies the average shape anisotropy into two, and
An analyzer for analyzing a solidified metal crystal structure.
前記観察画像の前記結晶粒を識別する結晶粒識別機構を更に有することを特徴とする金属凝固体結晶組織の解析装置。 In the analysis apparatus for the crystal structure of a solidified metal according to claim 7.
An analysis device for a metal solidified crystal structure, further comprising a crystal grain identification mechanism for identifying the crystal grains in the observation image.
前記金属凝固体結晶組織の前記観察画像を取得する結晶組織観察画像取得機構と、
前記観察画像の前記結晶粒の前記形状異方性量を算出する結晶粒形状異方性量算出機構と、
前記解析単位領域に基づいて前記結晶粒の前記部分面積を算出する結晶粒部分面積算出機構と、
前記形状異方性量と前記部分面積とから前記平均形状異方性量を算出する平均形状異方性量算出機構と、
前記平均形状異方性量の分布を作成する平均形状異方性量分布作成機構と、
前記平均形状異方性量を二分類する平均形状異方性量二値化機構と、
前記結晶粒の外形を近似した前記楕円の長軸方向と前記部分面積と前記平均形状異方性量とから平均形状異方性ベクトルを算出する平均形状異方性ベクトル算出機構と、
を有することを特徴とする金属凝固体結晶組織の解析装置。 An apparatus for carrying out the method for analyzing a solidified metal crystal structure according to claim 6.
A crystal structure observation image acquisition mechanism for acquiring the observation image of the metal solidified crystal structure,
A crystal grain shape anisotropy calculation mechanism for calculating the shape anisotropy of the crystal grains in the observation image,
A crystal grain partial area calculation mechanism that calculates the partial area of the crystal grain based on the analysis unit region,
An average shape anisotropy calculation mechanism that calculates the average shape anisotropy from the shape anisotropy and the partial area,
An average shape anisotropic amount distribution creation mechanism that creates a distribution of the average shape anisotropic amount, and
An average shape anisotropy binarization mechanism that classifies the average shape anisotropy into two, and
An average shape anisotropy vector calculation mechanism that calculates an average shape anisotropy vector from the long axis direction of the ellipse that approximates the outer shape of the crystal grain, the partial area, and the average shape anisotropy amount.
An analyzer for analyzing a solidified metal crystal structure.
前記観察画像の前記結晶粒を識別する結晶粒識別機構を更に有することを特徴とする金属凝固体結晶組織の解析装置。 In the analysis apparatus for the crystal structure of a solidified metal according to claim 9,
An analysis device for a metal solidified crystal structure, further comprising a crystal grain identification mechanism for identifying the crystal grains in the observation image.
前記解析単位領域を設定入力する入力機構と、
解析結果の出力を行う出力機構とを更に有することを特徴とする金属凝固体結晶組織の解析装置。 The analyzer for analyzing a solidified metal crystal structure according to any one of claims 7 to 10.
An input mechanism for setting and inputting the analysis unit area, and
An analysis device for a solidified metal crystal structure, which further has an output mechanism for outputting analysis results.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2019166769A JP2021043828A (en) | 2019-09-13 | 2019-09-13 | Analysis method of metal solidified crystal structure, and analyzing device of metal solidified crystal structure |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2019166769A JP2021043828A (en) | 2019-09-13 | 2019-09-13 | Analysis method of metal solidified crystal structure, and analyzing device of metal solidified crystal structure |
Publications (1)
Publication Number | Publication Date |
---|---|
JP2021043828A true JP2021043828A (en) | 2021-03-18 |
Family
ID=74863109
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2019166769A Pending JP2021043828A (en) | 2019-09-13 | 2019-09-13 | Analysis method of metal solidified crystal structure, and analyzing device of metal solidified crystal structure |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP2021043828A (en) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113894270A (en) * | 2021-09-26 | 2022-01-07 | 中国科学院金属研究所 | A method for predicting microporosity of superalloy equiaxed blades or structural castings |
-
2019
- 2019-09-13 JP JP2019166769A patent/JP2021043828A/en active Pending
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113894270A (en) * | 2021-09-26 | 2022-01-07 | 中国科学院金属研究所 | A method for predicting microporosity of superalloy equiaxed blades or structural castings |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111524137B (en) | Cell identification counting method and device based on image identification and computer equipment | |
CN115038965B (en) | Metallographic structure phase classification method, metallographic structure phase classification device, metallographic structure phase learning method, metallographic structure phase learning device, metallographic structure phase prediction method and metallographic material material property prediction device | |
US6835931B2 (en) | Chemical prefiltering for phase differentiation via simultaneous energy dispersive spectrometry and electron backscatter diffraction | |
JP7418639B2 (en) | Particle analysis data generation method, particle analysis data generation program, and particle analysis data generation device | |
CN115398228B (en) | Method and device for determining shooting conditions of metal structure, shooting method and device, phase classification method and device, material property prediction method and device | |
JP7547954B2 (en) | How to generate a quality attribute map | |
CN112435214A (en) | Pollen detection method and device based on prior frame linear scaling and electronic equipment | |
EP3770853B1 (en) | Image processing method, computer program, and recording medium | |
US9557299B2 (en) | Adaptive data collection for local states of a material | |
Wei et al. | 3D nanostructural characterisation of grain boundaries in atom probe data utilising machine learning methods | |
JP2021043828A (en) | Analysis method of metal solidified crystal structure, and analyzing device of metal solidified crystal structure | |
Mike et al. | Nonlandmark classification in paleobiology: computational geometry as a tool for species discrimination | |
Ghosh et al. | Developing virtual microstructures and statistically equivalent representative volume elements for polycrystalline materials | |
CN105527461A (en) | Material structure quantitative analysis method based on transmission electron microscope HAADF image | |
Calvat et al. | Learning Metal Microstructural Heterogeneity through Spatial Mapping of Diffraction Latent Space Features | |
JP2018125019A (en) | Image processing apparatus and image processing method | |
Boyadjian et al. | Detecting defects in materials using deep convolutional neural networks | |
CN114646653A (en) | In-situ electron diffraction data processing and analysis method and application based on machine learning | |
JP6329651B1 (en) | Image processing apparatus and image processing method | |
Matsumoto et al. | CT image-based 3D inflorescence estimation of Chrysanthemum seticuspe | |
JP7076052B1 (en) | Porosity measuring device and porosity measuring method | |
EP4538703A1 (en) | Device, method, program, and system for calculating alloy composition | |
EP3588435B1 (en) | Image processing method, computer program and recording medium | |
JP2019219310A (en) | Device and method for evaluating coating film | |
JP2004279292A (en) | Electromagnetic field analyzer, electromagnetic field analyzing method, and recording medium |