JPWO2020066348A1 - 判定方法 - Google Patents
判定方法 Download PDFInfo
- Publication number
- JPWO2020066348A1 JPWO2020066348A1 JP2020548139A JP2020548139A JPWO2020066348A1 JP WO2020066348 A1 JPWO2020066348 A1 JP WO2020066348A1 JP 2020548139 A JP2020548139 A JP 2020548139A JP 2020548139 A JP2020548139 A JP 2020548139A JP WO2020066348 A1 JPWO2020066348 A1 JP WO2020066348A1
- Authority
- JP
- Japan
- Prior art keywords
- phase difference
- sphere
- index value
- cells
- image
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 81
- 238000003384 imaging method Methods 0.000 claims abstract description 20
- 210000004027 cell Anatomy 0.000 claims description 70
- 239000002245 particle Substances 0.000 claims description 37
- 230000004083 survival effect Effects 0.000 claims description 24
- 210000000130 stem cell Anatomy 0.000 claims description 8
- 230000002093 peripheral effect Effects 0.000 claims description 2
- 230000003833 cell viability Effects 0.000 description 13
- 230000008569 process Effects 0.000 description 13
- 230000004069 differentiation Effects 0.000 description 12
- 238000010586 diagram Methods 0.000 description 11
- 230000006870 function Effects 0.000 description 11
- 238000009826 distribution Methods 0.000 description 10
- 238000001228 spectrum Methods 0.000 description 10
- 230000003287 optical effect Effects 0.000 description 7
- 238000003860 storage Methods 0.000 description 7
- 238000005311 autocorrelation function Methods 0.000 description 6
- 238000012545 processing Methods 0.000 description 6
- 238000005516 engineering process Methods 0.000 description 5
- 239000000523 sample Substances 0.000 description 5
- 230000035899 viability Effects 0.000 description 5
- 238000011156 evaluation Methods 0.000 description 4
- 210000002241 neurite Anatomy 0.000 description 4
- 210000001778 pluripotent stem cell Anatomy 0.000 description 4
- 238000004891 communication Methods 0.000 description 3
- 230000001066 destructive effect Effects 0.000 description 3
- 238000001093 holography Methods 0.000 description 3
- 239000013307 optical fiber Substances 0.000 description 3
- 108090000623 proteins and genes Proteins 0.000 description 3
- 238000012136 culture method Methods 0.000 description 2
- 238000012258 culturing Methods 0.000 description 2
- 230000013632 homeostatic process Effects 0.000 description 2
- 230000006698 induction Effects 0.000 description 2
- 230000001537 neural effect Effects 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 1
- 230000024245 cell differentiation Effects 0.000 description 1
- 230000001427 coherent effect Effects 0.000 description 1
- 230000000295 complement effect Effects 0.000 description 1
- 238000012888 cubic function Methods 0.000 description 1
- 230000006378 damage Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000005684 electric field Effects 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 239000007850 fluorescent dye Substances 0.000 description 1
- 210000001654 germ layer Anatomy 0.000 description 1
- 230000001939 inductive effect Effects 0.000 description 1
- 239000004973 liquid crystal related substance Substances 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 229910044991 metal oxide Inorganic materials 0.000 description 1
- 150000004706 metal oxides Chemical class 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 230000035699 permeability Effects 0.000 description 1
- 238000004886 process control Methods 0.000 description 1
- 238000012887 quadratic function Methods 0.000 description 1
- 239000013074 reference sample Substances 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
- 230000003595 spectral effect Effects 0.000 description 1
- 230000007480 spreading Effects 0.000 description 1
- 238000003892 spreading Methods 0.000 description 1
Images
Classifications
-
- C—CHEMISTRY; METALLURGY
- C12—BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
- C12Q—MEASURING OR TESTING PROCESSES INVOLVING ENZYMES, NUCLEIC ACIDS OR MICROORGANISMS; COMPOSITIONS OR TEST PAPERS THEREFOR; PROCESSES OF PREPARING SUCH COMPOSITIONS; CONDITION-RESPONSIVE CONTROL IN MICROBIOLOGICAL OR ENZYMOLOGICAL PROCESSES
- C12Q1/00—Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions
- C12Q1/02—Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions involving viable microorganisms
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
- G01N21/41—Refractivity; Phase-affecting properties, e.g. optical path length
- G01N21/45—Refractivity; Phase-affecting properties, e.g. optical path length using interferometric methods; using Schlieren methods
- G01N21/453—Holographic interferometry
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N33/00—Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
- G01N33/48—Biological material, e.g. blood, urine; Haemocytometers
- G01N33/50—Chemical analysis of biological material, e.g. blood, urine; Testing involving biospecific ligand binding methods; Immunological testing
- G01N33/5005—Chemical analysis of biological material, e.g. blood, urine; Testing involving biospecific ligand binding methods; Immunological testing involving human or animal cells
- G01N33/5008—Chemical analysis of biological material, e.g. blood, urine; Testing involving biospecific ligand binding methods; Immunological testing involving human or animal cells for testing or evaluating the effect of chemical or biological compounds, e.g. drugs, cosmetics
- G01N33/502—Chemical analysis of biological material, e.g. blood, urine; Testing involving biospecific ligand binding methods; Immunological testing involving human or animal cells for testing or evaluating the effect of chemical or biological compounds, e.g. drugs, cosmetics for testing non-proliferative effects
- G01N33/5032—Chemical analysis of biological material, e.g. blood, urine; Testing involving biospecific ligand binding methods; Immunological testing involving human or animal cells for testing or evaluating the effect of chemical or biological compounds, e.g. drugs, cosmetics for testing non-proliferative effects on intercellular interactions
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N33/00—Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
- G01N33/48—Biological material, e.g. blood, urine; Haemocytometers
- G01N33/50—Chemical analysis of biological material, e.g. blood, urine; Testing involving biospecific ligand binding methods; Immunological testing
- G01N33/5005—Chemical analysis of biological material, e.g. blood, urine; Testing involving biospecific ligand binding methods; Immunological testing involving human or animal cells
- G01N33/5008—Chemical analysis of biological material, e.g. blood, urine; Testing involving biospecific ligand binding methods; Immunological testing involving human or animal cells for testing or evaluating the effect of chemical or biological compounds, e.g. drugs, cosmetics
- G01N33/5044—Chemical analysis of biological material, e.g. blood, urine; Testing involving biospecific ligand binding methods; Immunological testing involving human or animal cells for testing or evaluating the effect of chemical or biological compounds, e.g. drugs, cosmetics involving specific cell types
- G01N33/5073—Stem cells
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/18—Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/0002—Inspection of images, e.g. flaw detection
- G06T7/0012—Biomedical image inspection
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/60—Analysis of geometric attributes
- G06T7/62—Analysis of geometric attributes of area, perimeter, diameter or volume
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V20/00—Scenes; Scene-specific elements
- G06V20/60—Type of objects
- G06V20/69—Microscopic objects, e.g. biological cells or cellular parts
- G06V20/695—Preprocessing, e.g. image segmentation
-
- C—CHEMISTRY; METALLURGY
- C12—BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
- C12M—APPARATUS FOR ENZYMOLOGY OR MICROBIOLOGY; APPARATUS FOR CULTURING MICROORGANISMS FOR PRODUCING BIOMASS, FOR GROWING CELLS OR FOR OBTAINING FERMENTATION OR METABOLIC PRODUCTS, i.e. BIOREACTORS OR FERMENTERS
- C12M1/00—Apparatus for enzymology or microbiology
- C12M1/34—Measuring or testing with condition measuring or sensing means, e.g. colony counters
-
- G—PHYSICS
- G03—PHOTOGRAPHY; CINEMATOGRAPHY; ANALOGOUS TECHNIQUES USING WAVES OTHER THAN OPTICAL WAVES; ELECTROGRAPHY; HOLOGRAPHY
- G03H—HOLOGRAPHIC PROCESSES OR APPARATUS
- G03H1/00—Holographic processes or apparatus using light, infrared or ultraviolet waves for obtaining holograms or for obtaining an image from them; Details peculiar thereto
- G03H1/04—Processes or apparatus for producing holograms
- G03H1/08—Synthesising holograms, i.e. holograms synthesized from objects or objects from holograms
- G03H1/0866—Digital holographic imaging, i.e. synthesizing holobjects from holograms
-
- G—PHYSICS
- G03—PHOTOGRAPHY; CINEMATOGRAPHY; ANALOGOUS TECHNIQUES USING WAVES OTHER THAN OPTICAL WAVES; ELECTROGRAPHY; HOLOGRAPHY
- G03H—HOLOGRAPHIC PROCESSES OR APPARATUS
- G03H1/00—Holographic processes or apparatus using light, infrared or ultraviolet waves for obtaining holograms or for obtaining an image from them; Details peculiar thereto
- G03H1/04—Processes or apparatus for producing holograms
- G03H1/08—Synthesising holograms, i.e. holograms synthesized from objects or objects from holograms
- G03H1/0866—Digital holographic imaging, i.e. synthesizing holobjects from holograms
- G03H2001/0883—Reconstruction aspect, e.g. numerical focusing
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10056—Microscopic image
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10064—Fluorescence image
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20048—Transform domain processing
- G06T2207/20056—Discrete and fast Fourier transform, [DFT, FFT]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20076—Probabilistic image processing
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30004—Biomedical image processing
- G06T2207/30024—Cell structures in vitro; Tissue sections in vitro
Landscapes
- Engineering & Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Chemical & Material Sciences (AREA)
- Immunology (AREA)
- General Health & Medical Sciences (AREA)
- Theoretical Computer Science (AREA)
- Biomedical Technology (AREA)
- Molecular Biology (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Pathology (AREA)
- Urology & Nephrology (AREA)
- Hematology (AREA)
- Data Mining & Analysis (AREA)
- Organic Chemistry (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Microbiology (AREA)
- Biotechnology (AREA)
- Cell Biology (AREA)
- Zoology (AREA)
- Wood Science & Technology (AREA)
- Mathematical Analysis (AREA)
- Proteomics, Peptides & Aminoacids (AREA)
- General Engineering & Computer Science (AREA)
- Mathematical Physics (AREA)
- Pure & Applied Mathematics (AREA)
- Computational Mathematics (AREA)
- Mathematical Optimization (AREA)
- Quality & Reliability (AREA)
- Radiology & Medical Imaging (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Medical Informatics (AREA)
- Food Science & Technology (AREA)
- Medicinal Chemistry (AREA)
- Toxicology (AREA)
- Tropical Medicine & Parasitology (AREA)
Abstract
Description
開示の技術は、複数の細胞の凝集体の状態を判定する判定方法に関する。
細胞の状態を評価または判定する技術として、例えば、以下の技術が知られている。特許文献1には、一細胞の表面の平坦度又は細胞集団の表面の平坦度を分化の度合いの指標として用いた多能性幹細胞の分化の度合いを判別する方法が記載されている。
特許文献2には、撮像画像における輝度に基づいて、分化した多能性幹細胞を含む分化コロニーと未分化多能性幹細胞のみを含む未分化コロニーと多層に重なった多能性幹細胞を含む多層コロニーとを識別する方法が記載されている。この方法では、輝度の第1の閾値より明るい輝度の領域を有するコロニーを分化コロニーであると判断する。また、第1の閾値に等しいかこれより暗い輝度の領域のみを有するコロニーを未分化コロニーであると判断する。また、第1の閾値に等しいかこれより暗く第2の閾値に等しいかこれより明るい輝度の領域のみを有するコロニーを未分化コロニーであると判断する。また、第2の閾値より暗い輝度の領域を有するコロニーを多層コロニーであると判断する。
特許文献3には、神経分化過程の細胞を撮像した撮像画像を入力する画像入力ステップと、神経分化過程の細胞に出現する神経突起を撮像画像に基づく原画像から抽出する神経突起抽出ステップと、抽出された神経突起の状態を判定する神経突起対応判定ステップと、を備えることを特徴とする細胞評価方法が記載されている。
特許文献4には、細胞状態を提示する方法であって、該細胞に由来する遺伝子から選択される少なくとも1つの遺伝子に関連する遺伝子状態を経時的にモニタすることより該細胞の経時プロファイルを得る工程と該経時プロファイルを提示する工程とを含む方法が記載されている。
細胞の大量生産が可能な培養手法として、細胞の凝集体であるスフェアを培地中に浮遊させた状態で培養する三次元培養法が知られている。三次元培養による細胞の製造工程においては、工程管理の容易性の観点から、細胞の品質をスフェアの状態のまま非破壊且つ簡便に評価する技術が求められる。しかしながら、現時点においては、三次元空間中にランダムに存在する様々なサイズのスフェアを評価する手法が確立されておらず、特にスフェアの内部における細胞の密度及び生存状況について直接観察することは困難である。このため、特許文献1−3にも記載されているように、従来の二次元培養の手法を適用した評価がなされているが、培養する細胞の数の増加に伴って、評価工数が増えるため、多くの人手及び多くの時間が必要となる。従来の二次元培養の手法を適用した評価では、スフェアを単一細胞に分解したり、特許文献4に記載されているように、蛍光色素を添加したりするといった細胞の破壊を伴う処理が必要とされる。
開示の技術は、三次元培養により形成される複数の細胞の凝集体の状態を非破壊且つ簡便に判定することを目的とする。
開示の技術に係る判定方法は、複数の細胞の凝集体を撮像したホログラムから凝集体の位相差画像を生成し、位相差画像を構成する複数の画素における位相差量の配列のランダム性を示す第1の指標値を導出し、第1の指標値に基づいて、凝集体を構成する細胞の状態を判定することを含む。開示の技術に係る判定方法によれば、三次元培養により形成される複数の細胞の凝集体の状態を非破壊且つ簡便に判定することが可能となる。
第1の指標値は、位相差画像における同一の位相差量の画素を結んだ等位相線によって囲まれた領域の形状の、円形からの乖離の程度に応じて定まる値であってもよい。具体的には、位相差画像を構成する複数の画素における所定範囲の位相差量の最小値をΦ0、最大値をΦNとし、上記所定範囲内における任意の位相Φにおける等位相線の周長をL(Φ)とし、周長L(Φ)の等位相線によって囲まれる領域の面積をA(Φ)とした場合、下記の(3)式によって定義される平均位相ゆらぎλを、第1の指標値として用いることができる。第1の指標値として平均位相ゆらぎλを用いることによって、位相差画像を構成する複数の画素における位相差量の配列のランダム性を定量化することができ、従って、凝集体を構成する細胞の状態を適確に判定することが可能となる。
また、位相差画像から凝集体の形状に依存する成分を除去する処理を施した形状成分除去画像に基づいて第1の指標値を導出してもよい。例えば、形状成分除去画像に基づいて導出された自己相関関数または二次元パワースペクトルを、第1の指標値として導出してもよい。
開示の技術に係る判定方法は、第1の指標値に基づいて、凝集体を構成する細胞の生存率に関する判定を行うことを含み得る。また、凝集体を構成する細胞が幹細胞である場合、開示の技術に係る判定方法は、第1の指標値に基づいて、凝集体を構成する幹細胞の未分化率に関する判定を行うことを含み得る。細胞の生存率または未分化率に関する判定を、第1の指標値に基づいて行うことによって、当該判定を非破壊且つ簡便に行うことが可能となる。
開示の技術に係る判定方法は、判定対象ロットに含まれる複数の凝集体について、第1の指標値と凝集体の粒径との相関性を示す第2の指標値を導出し、第2の指標値に基づいて、判定対象ロットについて判定を行うことを含み得る。これにより、判定対象ロットについての判定を、非破壊且つ簡便に行うことが可能となる。
開示の技術に係る判定方法は、第2の指標値に基づいて、判定対象ロットに含まれる細胞の生存率に関する判定を行うことを含み得る。また、凝集体を構成する細胞が幹細胞である場合、開示の技術に係る判定方法は、第2の指標値に基づいて、判定対象ロットに含まれる幹細胞の未分化率に関する判定を行うことを含み得る。判定対象ロットに含まれる細胞の生存率または未分化率に関する判定を、第2の指標値に基づいて行うことによって、当該判定を非破壊且つ簡便に行うことが可能となる。
開示の技術によれば、三次元培養により形成される複数の細胞の凝集体の状態を非破壊且つ簡便に判定することが可能となる。
以下、本発明の実施形態について図面を参照しつつ説明する。尚、各図面において、実質的に同一又は等価な構成要素又は部分には同一の参照符号を付している。
開示の技術の実施形態に係る判定方法は、複数の細胞の凝集体(スフェア)を撮像したホログラムから凝集体(スフェア)の位相差画像を生成し、位相差画像を構成する複数の画素における位相差量の配列のランダム性を示す第1の指標値を導出し、この第1の指標値に基づいて、凝集体(スフェア)を構成する細胞の状態を判定するという値である。この判定方法によれば、凝集体(スフェア)の状態の判定を非破壊的且つ簡便に行うことが可能である。
図1は、開示の技術の実施形態に係る判定方法の実施に用いる撮像システム1の構成の一例を示す図である。撮像システム1は、公知のデジタルホログラフィ技術を用いてスフェアのホログラムを取得するためのホログラム光学系10を含んで構成されている。
デジタルホログラフィ技術は、物体を透過または反射した物体光と、物体光に対してコヒーレントである参照光との干渉によって生じる像をイメージセンサで撮像し、撮像によって得られた画像について、光伝搬に基づく数値計算を実施することによって、物体からの光波の波面を復元する技術である。デジタルホログラフィ技術によれば、物体の位相分布を定量化し、また、焦点位置を機械的に移動させることなく、物体の三次元情報を取得することができる。
ホログラム光学系10は、レーザ光源11、ビームスプリッタ12、18、コリメートレンズ13、21、対物レンズ15、ダイクロイックミラー34、結像レンズ17及びCMOS(Complementary Metal Oxide Semiconductor)カメラ19を含んで構成されている。サンプルステージにセットされる試料14としてのスフェアは、コリメートレンズ13と対物レンズ15との間に配置される。
レーザ光源11には、例えば波長632.8nmのHeNeレーザを用いることが可能である。レーザ光源11から出射されたレーザ光は、ビームスプリッタ12により、2つのレーザ光に分割される。2つのレーザ光の一方は物体光とされ、他方は参照光とされる。物体光は、コリメートレンズ13によって平行光とされた後、サンプルステージにセットされた試料14であるスフェアに照射される。スフェアを透過した物体光による像は、対物レンズ15によって拡大される。対物レンズ15を透過した物体光は、ダイクロイックミラー34を透過して、結像レンズ17によって再び平行光とされた後、ビームスプリッタ18を介してCMOSカメラ19の撮像面に結像される。一方、参照光は、光ファイバ20によってコリメートレンズ21の手前まで導かれる。光ファイバ20から出射した参照光は、コリメートレンズ21によって平行光とされ、ビームスプリッタ18を介してCMOSカメラ19の撮像面に入射する。物体光と参照光との干渉によって生じるホログラムが、CMOSカメラ19によって記録される。なお、CMOSカメラ19の撮像面に入射する物体光と参照光の光軸方向が互いに異なったオフアキシャル光学系が構成されていてもよい。
本実施形態に係る撮像システム1によれば、スフェアを破壊することなく、またスフェアを構成する細胞にダメージを与えることなくスフェアの位相差画像を取得することができる。なお、上記した撮像システム1の構成は、一例に過ぎず、上記の構成に限定されるものではない。開示の技術に係る判定方法の実施には、デジタルホログラム技術を用いてホログラムを取得することができる、あらゆる撮像システムを利用することが可能である。
以下に、撮像システム1を用いて取得したスフェアのホログラムから、スフェアの位相差画像を取得する方法の一例について説明する。
はじめに、CMOSカメラ19によって取得した図2Aに例示するホログラムを、二次元フーリエ変換することにより、物体光のみの複素振幅成分を抜き出す。図2Bは、この処理によって得られるスフェアのフーリエ変換画像の一例である。
次に、例えば角スペクトル法を適用して任意の空間位置のスフェアの位相を示す画像を復元する。図2Cは、この処理によって得られるスフェアのアンラッピング前の位相差画像の一例である。この時点におけるスフェアの位相は、0〜2πの値に畳みこまれている。そこで、例えば、Unweighted Least Squares(重みなし最小2乗法)またはFlynn's Algorithm(フリンのアルゴリズム)などの位相接続(アンラッピング)手法を適用して2π以上の部分も接合していくことで、図2Dに例示するような最終的なスフェアの位相差画像を得ることができる。なお、アンラッピングの手法は数多く提案されており、位相不整合を生じない適切な手法を適宜選択すればよい。
図3は、位相差画像IPの概念を示す図である。図3の下段は、位相差画像IPの各画素kにおける位相差量を3次元表示した画像である。図3の上段は、位相差画像IPの各画素kにおける位相差量を平面上にグレースケールで示した図である。
ここで、位相差画像IPにおける位相差量θは、位相差画像IPの同一焦点面内に存在するバックグラウンド(スフェアの存在しない領域)の位相をθBとし、スフェアの存在する領域の位相をθSとした場合、下記の(1)式によって表わされる。また、本明細書中における「位相」という用語は、光を電磁波とみなした場合の電場振幅の位相であり、より一般的な意味で使用される。
また、位相差画像IPの各画素kにおける位相差量θkは、下記(2)式によって表わすことができる。但し、nkは位相差画像IPの各画像kに対応する部位におけるスフェアの屈折率であり、dkは位相差画像IPの各画素kに対応する部位におけるスフェアの厚さであり、λはホログラム光学系10における物体光の波長である。
スフェアの位相差画像は、スフェアを透過した物体光の光路長分布を示した画像である。スフェア内における光路長は、スフェアの屈折率とスフェアの厚さの積に相当することから、スフェアの位相差画像は、(2)式にも示されているように、スフェアの屈折率及び厚さ(形状)の情報を含んでいる。
スフェアに対して焦点が合っていない位相差画像からは、回折による広がりの影響によりスフェアの実態に合致した正確な情報が得られない。従って、CMOSカメラ19によって取得したホログラムから位相差画像を取得する際に、スフェアに焦点を合わせることが好ましい。ここで、「スフェアに焦点を合わせる」とは、球形状のスフェアの中央付近でスライスした位相差画像を得ることを意味する。スフェアに焦点が合った位相差画像を用いてスフェアの状態を判定することによって、より正確な判定結果を得ることができる。
位相差画像の焦点合わせは、人手によらず、自動化することが好ましい。焦点合わせを自動化することによって、作業者による任意性を排除し、更に、処理時間の短縮を図ることができる。本発明者らは、以下に説明する自動化可能な焦点合わせの手法を見出した。
図4の左のグラフは、スフェアの位相差画像における平面方向の位置と位相差量との関係の一例を示すグラフであり、実線がスフェアに焦点が合っている状態に対応し、点線がスフェアに焦点が合っていない状態に対応する。スフェアに焦点が合っている場合、位相差画像における特定の位置に急峻なピークが現れる。一方、スフェアに焦点が合っていない場合、焦点が合っている場合と比較してピークが低く且つなだらかになる。
図4の右のグラフは、スフェアの位相差画像における位相差量のヒストグラムの一例であり、実線がスフェアに焦点が合っている状態に対応し、点線がスフェアに焦点が合っていない状態に対応する。スフェアに焦点が合っている場合、カーブの半値幅w(位相差量のばらつき)は、相対的に大きくなり、スフェアに焦点が合っていない場合、カーブの半値幅w(位相差量のばらつき)は、相対的に小さくなる。
従って、互いに異なる焦点位置(スライス位置)毎にスフェアの位相差画像を取得し、取得した位相差画像の各々について、位相差量のヒストグラムにおけるカーブの半値幅w(位相差量のばらつき)を求め、求めた半値幅wのうち、最大の半値幅wを有する位相差画像を、スフェアに焦点が合った位相差画像として抽出することで焦点合わせを実現できる。
上記の焦点合わせは、コンピュータを用いて自動化することが可能である。図5は、上記の焦点合わせを自動で行うオートフォーカス処理を行うコンピュータ500のハードウェア構成の一例である。
コンピュータ500は、CPU(Central Processing Unit)501、一時記憶領域としての主記憶装置502、不揮発性の補助記憶装置503、CMOSカメラ19との間において通信を行うための通信I/F(InterFace)504、及び液晶ディスプレイ等の表示部505を含んで構成されている。CPU501、主記憶装置502、補助記憶装置503、通信I/F504、及び表示部505は、それぞれ、バス507に接続されている。補助記憶装置503には、上記のオートフォーカス処理の手順を記述したオートフォーカスプログラム506が格納されている。コンピュータ500は、CPU501がオートフォーカスプログラム506を実行することによって、オートフォーカス処理を行う。
図6は、コンピュータ500において実施されるオートフォーカス処理の流れの一例を示すフローチャートである。
ステップS1において、CPU501は、CMOSカメラ19からスフェアのホログラムを取得する。
ステップS2において、CPU501は、取得したホログラムから、焦点位置(スライス位置)が互いに異なる複数の位相差画像を生成する。
ステップS3において、CPU501は、焦点位置(スライス位置)毎の位相差画像の各々について、位相差量のばらつきを導出する。CPU501は、例えば、位相差画像における位相差量の最大値と最小値との差分を、当該位相差画像における位相差量のばらつきとして導出してもよい。
ステップS4において、CPU501は、焦点位置(スライス位置)が互いに異なる複数の位相差画像のうち、ステップS3にて導出した位相差量のばらつきが最大となる位相差画像を、焦点が合った位相差画像として抽出する。
図7は、スフェアの位相差画像における焦点位置(スライス位置)と位相差量のばらつきとの関係の一例を示したグラフである。図7には、焦点位置が、−400μm、−200μm、0μm、+200μm、及び+400μmに対応するスフェアの位相差画像が、グラフとともに例示されている。なお、図7では、位相差量のばらつきが最大となる焦点位置を0μmとしている。上記のオートフォーカス処理によれば、位相差量のばらつきが最大となる焦点位置0μmに対応する位相差画像が、焦点が合った位相差画像として抽出される。位相差量のばらつきが最大となる焦点位置0μmに対応する位相差画像において、スフェアの輪郭が最も鮮明となる。
上記のように、開示の技術の実施形態に係る判定方法においては、スフェアの位相差画像を構成する複数の画素における位相差量の配列のランダム性を示す第1の指標値を導出し、この第1の指標値に基づいて、スフェアを構成する細胞の状態を判定することを含む。
図8Aは、未分化状態を維持しているiPS細胞の、ロット内における割合(以下、未分化率という)が99%である培養ロットから抽出されたiPS細胞の凝集体であるスフェアの位相差画像の代表例を示す図である。図8Bは、図8Aに示す位相差画像における位相差量の分布を示すコンター図である。図9Aは、未分化率が87%である培養ロットから抽出されたiPS細胞の凝集体であるスフェアの位相差画像の代表例を示す図である。図9Bは、図9Aに示す位相差画像における位相差量の分布を示すコンター図である。
未分化率が相対的に高いロット内のスフェアは、未分化率が相対的に低い(すなわち、胚葉への分化が進行している)ロット内のスフェアと比較して、内部の均質性が高いと考えられる。従って、未分化率が相対的に高いロット内のスフェアの位相差画像においては、図8Bに示すように、同一の位相差量の画素を結んだ等位相線によって囲また領域の形状が円形に近い形状となり、等位相線が同心円状に分布する。一方、未分化率が相対的に低いロット内のスフェアの位相差画像においては、図9Bに示すように、位相差画像を構成する複数の画素における位相差量の配列がランダムとなり、等位相線によって囲まれた領域の形状がくずれる。すなわち、スフェアを構成する細胞の分化の進行状態は、スフェアの位相差画像を構成する複数の画素における位相差量の配列のランダム性に反映される。従って、スフェアの位相差画像を構成する複数の画素における位相差量の配列のランダム性を示す指標値によって、スフェアを構成する細胞の分化の進行状態を定量化することが可能である。
スフェアの位相差画像を構成する複数の画素における位相差量の配列のランダム性を示す指標値として、例えば、スフェアの位相差画像における同一の位相差量の画素を結んだ等位相線によって囲まれた領域の形状の、円形からの乖離の程度に応じて定まる指標値を用いることができる。
具体的には、スフェアの位相差画像を構成する複数の画素における位相差量の配列のランダム性を示す指標値として、下記の(3)式によって定義される平均位相ゆらぎλを用いることができる。(3)式において、Φ0は、位相差画像を構成する複数の画素における所定範囲の位相差量の最小値であり、ΦNは、上記所定範囲の位相差量の最大値である。L(Φ)は、上記所定範囲内における任意の位相Φにおける等位相線の周長であり、A(Φ)は、周長L(Φ)の等位相線によって囲まれる領域の面積である。
平均位相ゆらぎλは、等位相線によって囲まれる領域の形状が真円である場合に最小となり、等位相線によって囲まれる領域の形状が円形から離れる程大きくなる程大きくなる。すなわち、等位相線によって囲まれる領域の形状のランダム性が高くなる程、平均位相ゆらぎλは大きくなる。
図8A及び図8Bに示される、未分化率99%に対応するスフェアの位相差画像における平均位相ゆらぎλは13.6であり、図9A及び図9Bに示される、未分化率87%に対応するスフェアの位相差画像における平均位相ゆらぎλは18.9であった。未分化率が相対的に低いロット内のスフェアは、未分化状態を逸脱した複数の細胞を内部に含むことにより、スフェア内部の均質性が低下する。その結果、スフェアの位相差画像を構成する複数の画素における位相差量の配列のランダム性が増大し、平均位相ゆらぎλが大きくなる。
このように、スフェアの位相差画像を構成する複数の画素における位相差量の配列のランダム性を示す指標値として、平均位相ゆらぎλを用いることによって、スフェアを構成する細胞の分化の進行状態を、細胞を破壊することなく推定することが可能である。
平均位相ゆらぎλは、細胞の分化の進行状態の判定のみならず、スフェア内における細胞の生存率の判定に用いることもできる。図10Aは、ロット内における細胞の生存率が87.3%であるロットから抽出されたiPS細胞の凝集体であるスフェアの位相差画像の代表例を示す図である。図10Bは、図10Aに示す位相差画像における位相差量の分布を示すコンター図である。図11Aは、ロット内における細胞の生存率が59.2%である培養ロットから抽出されたiPS細胞の凝集体であるスフェアの位相差画像の代表例を示す図である。図11Bは、図11Aに示す位相差画像における位相差量の分布を示すコンター図である。
健全な細胞は、その恒常性から内部の屈折率は、媒質の屈折率とは異なる一定の値を維持するものと考えられる。一方、死細胞は、恒常性を喪失し、内部の屈折率が媒質の屈折率と略同じになるものと考えられる。従って、細胞の生存率が相対的に高い培養ロット内のスフェアは、細胞の生存率が相対的に低い培養ロット内のスフェアと比較して、内部の均質性が高いと考えられる。従って、細胞の生存率が相対的に高い培養ロット内のスフェアの位相差画像においては、図10Bに示すように、同一の位相差量の画素を結んだ等位相線によって囲まれた領域の形状が円形に近い形状となり、等位相線が同心円状に分布する。一方、細胞の生存率が相対的に低い培養ロット内のスフェアの位相差画像においては、図11Bに示すように、当該位相差画像を構成する複数の画素における位相差量の配列がランダムとなり、等位相線によって囲まれた領域の形状がくずれる。すなわち、スフェアを構成する細胞の生存率は、スフェアの位相差画像を構成する複数の画素における位相差量の配列のランダム性に反映される。従って、平均位相ゆらぎλを、スフェア内の細胞の生存率の判定に用いることが可能である。
図10A及び図10Bに示される、細胞の生存率87.3%に対応するスフェアの位相差画像における平均位相ゆらぎλは5.77であり、図11A及び図11Bに示される、細胞の生存率59.2%に対応するスフェアの位相差画像における平均位相ゆらぎλは11.08であった。細胞の生存率が相対的に低いロット内のスフェアは、多くの死細胞を内部に含むことによりスフェア内部の均質性が低下する。その結果、スフェアの位相差画像を構成する複数の画素における位相差量の配列のランダム性が増大し、平均位相ゆらぎλが大きくなる。
このように、スフェアの位相差画像を構成する複数の画素における位相差量の配列のランダム性を示す指標値として、平均位相ゆらぎλを用いることによって、スフェアを構成する細胞の生存率を、細胞を破壊することなく推定することが可能である。
開示の技術の実施形態に係る判定方法は、判定対象ロットに含まれる複数のスフェアについて、スフェアの位相差画像を構成する複数の画素における位相差量の配列のランダム性を示す上記の指標値(以下、第1の指標値という)とスフェアの粒径との相関性を示す第2の指標値を導出し、この第2の指標値に基づいて、判定対象ロットについて判定を行うことを含み得る。第1の指標値としては、平均位相ゆらぎλを用いることができる。
図12は、未分化率がそれぞれ、87%及び99%である2つのロットのそれぞれに含まれる複数のスフェアについて取得した、スフェア粒径と平均位相ゆらぎλとの相関性を示すグラフである。図12に示すように、スフェア粒径が大きくなる程、平均位相ゆらぎλが大きくなる傾向があることが判明した。平均位相ゆらぎλのスフェア粒径依存性は、細胞株、三次元培養プロセス及び分化誘導プロセスの特徴を表わしていると考えることができる。具体的には、ガス濃度及び誘導因子の透過性に起因して、粒径が小さいスフェアは分化の進行が均一となりやすい一方、粒径が大きいスフェアは分化の進行が不均一となりやすいこと、及び培養及びまたは分化誘導プロセスに応じてスフェア内の細胞密度のサイズ依存性が異なるといったことが、平均位相ゆらぎλがスフェア粒径依存性を持つ要因となると考えられる。
また、図12に示すように、未分化率87%のロットと、未分化率99%のロットとでは、スフェア粒径と平均位相ゆらぎλとの相関性に差異があることが判明した。すなわち、ロット間における未分化率の違いは、スフェア粒径と平均位相ゆらぎλとの相関性に反映される。図12に例示されるスフェア粒径と平均位相ゆらぎλとの相関性は、例えば、下記の(4)式に示される関数によってフィッティングすることが可能である。すなわち、スフェア粒径と平均位相ゆらぎλとの相関性は、(4)式の関数による近似式によって表わすことができる。(4)式においてXはスフェア粒径であり、Yは平均位相ゆらぎλであり、Aは定数である。(4)式における定数Aは、スフェア粒径と平均位相ゆらぎλとの相関性を示す(相関性を特徴付ける)指標値(第2の指標値)として用いることができる。
図12に示される未分化率87%のロットにおける スフェア粒径と平均位相ゆらぎλとの相関性は、(4)式の関数によるフィッティングにより、Y=e0.0122Xと表わすことができる。すなわち、未分化率87%のロットにおけるスフェア粒径と平均位相ゆらぎλとの相関性を特徴付ける指標値(第2の指標値)である定数Aは、0.0122である。また、図12に示される未分化率99%のロットにおけるスフェア粒径と平均位相ゆらぎλとの相関性は、(4)式の関数によるフィッティングにより、Y=e0.0107Xと表わすことができる。すなわち、未分化率99%のロットにおけるスフェア粒径と平均位相ゆらぎλとの相関性を特徴付ける指標値(第2の指標値)である定数Aは、0.0107である。
図13は、図12に示す、未分化率が互いに異なる2つのロットにおける、(4)式の定数Aと、未分化率との相関性を示すグラフである。図13に示すように、定数Aと未分化率との間には、一定の相関性がみられることから、例えば、定数Aに基づいて、判定対象ロット内における細胞の未分化率を推定することが可能である。このように、平均位相ゆらぎλ(第1の指標値)とスフェア粒径との相関性を示す指標値である定数A(第2の指標値)に基づいて、判定対象ロットについて未分化率を推定することが可能であり、また、推定された未分化率に基づいて、判定対象ロットについて、例えば良否判定を行うことが可能である。
図14は、ロット内における細胞の生存率がそれぞれ、59.2%、69.8%、81.5%、87.3%である4つのロットのそれぞれに含まれる複数のスフェアについて取得した、スフェア粒径と平均位相ゆらぎλとの相関性を示すグラフである。平均位相ゆらぎλが、スフェア粒径依存性を有する点は、上記したとおりである。図14に示すように、生存率が互いに異なるロット間において、スフェア粒径と平均位相ゆらぎλとの相関性に差異があることが判明した。すなわち、ロット間における細胞の生存率の違いは、スフェア粒径と平均位相ゆらぎλとの相関性に反映される。図14に示されるスフェア粒径と平均位相ゆらぎλとの相関性は、例えば、(4)式に示される関数によってフィッティングすることが可能である。すなわち、スフェア粒径と平均位相ゆらぎλとの相関性は、(4)式の関数による近似式によって表わすことができる。
図14に示される生存率59.2%のロットにおけるスフェア粒径と平均位相ゆらぎλとの相関性は、(4)式の関数によるフィッティングにより、Y=e0.0183Xと表わすことができる。すなわち、生存率59.2%のロットにおけるスフェア粒径と平均位相ゆらぎλとの相関性を特徴付ける指標値(第2の指標値)である定数Aは、0.0183である。また、図14に示される生存率69.8%のロットにおけるスフェア粒径と平均位相ゆらぎλとの相関性は、(4)式の関数によるフィッティングにより、Y=e0.0162Xと表わすことができる。すなわち、生存率69.8%のロットにおけるスフェア粒径と平均位相ゆらぎλとの相関性を特徴付ける指標値(第2の指標値)である定数Aは、0.0162である。図14に示される生存率81.5%のロットにおけるスフェア粒径と平均位相ゆらぎλとの相関性は、(4)式の関数によるフィッティングにより、Y=e0.0132Xと表わすことができる。すなわち、生存率81.5%のロットにおけるスフェア粒径と平均位相ゆらぎλとの相関性を特徴付ける指標値(第2の指標値)である定数Aは、0.0132である。また、図14に示される生存率87.3%のロットにおけるスフェア粒径と平均位相ゆらぎλとの相関性は、(4)式の関数によるフィッティングにより、Y=e0.0126Xと表わすことができる。すなわち、生存率87.3%のロットにおけるスフェア粒径と平均位相ゆらぎλとの相関性を特徴付ける指標値(第2の指標値)である定数Aは、0.0126である。
図15は、図14に示す、細胞の生存率が互いに異なる4つのロットにおける、(4)式の定数Aと、細胞の生存率との相関性を示すグラフである。図15に示すように、定数Aと細胞の生存率との間には、一定の相関性がみられることから、定数Aに基づいて、判定ロット内における細胞の生存率を推定することが可能である。このように、平均位相ゆらぎλ(第1の指標値)とスフェア粒径との相関性を示す指標値である定数A(第2の指標値)に基づいて、判定対象ロットについてロット内における細胞の生存率を推定することが可能であり、また、推定された細胞の生存率に基づいて、判定対象ロットについて、例えば良否判定を行うことが可能である。
以上の説明では、スフェアの位相差画像を構成する複数の画素における位相差量の配列のランダム性を示す指標値として、(3)式によって定義される平均位相ゆらぎλを用いる場合を例示したが、この態様に限定されない。
例えば、スフェアの位相差画像からスフェアの形状に依存する成分を除去する処理を施した形状成分除去画像に基づいて、スフェアの位相差画像を構成する複数の画素における位相差量の配列のランダム性を示す指標値を導出してもよい。ここで、スフェアの形状に依存する成分とは、スフェアの位相差画像において、スフェア表面の基本形状を構成する、画素サイズに対して緩やかな曲面であり、二次関数及び三次関数等の多項式曲面による近似によって導出される曲面成分である。この成分を、スフェアの位相差画像から除去することによって、スフェアの位相差画像における、位相差量のランダム性(空間ばらつき)のみを抽出した形状成分除去を得ることができる。図16は、スフェアの位相差画像(左図)からスフェアの形状に依存する成分を除去する処理を施すことにより、位相差量のランダム性(空間ばらつき)を抽出した形状成分除去画像(右図)の一例を示す図である。
例えば、形状成分除去画像(図16右図)について導出した自己相関関数を、スフェアの位相差画像を構成する複数の画素における位相差量の配列のランダム性を示す指標値として用いることが可能である。すなわち、自己相関関数によって、スフェアの位相差画像における位相差量のランダム性(空間ばらつき)を定量化することができる。従って、形状成分除去画像(図16右図)について導出した自己相関関数に基づいて、当該スフェアまたは当該スフェアを含む判定対象ロットについて、細胞の未分化率及び生存率を推定することが可能である。また、上記のようにして得られる自己相関関数を、例えば、細胞の生存率及び未分化率が既知である基準サンプルと比較することによって、当該スフェアまたは当該スフェアを含む判定対象ロットについて、例えば良否判定を行うことが可能である。
また、形状成分除去画像(図16右図)について導出した二次元パワースペクトルを、スフェアの位相差画像を構成する複数の画素における位相差量の配列のランダム性を示す指標値として用いることが可能である。位相差画像Φ(x,y)の二次元フーリエ変換スペクトルΦ(kx,ky)における二次元パワースペクトルPは、下記の(5)式によって表わされる。ここで、kx及びkyは、空間周波数である。
図17は、形状成分除去画像(図16右図)について導出した二次元パワースペクトル(ky=0とした)を示す図である。例えば、図17に例示されるパワースペクトルの傾き及び、関数フィテッィングにより、当該スフェアの位相差画像における位相差量のランダム性(空間ばらつき)を定量化することができる。そして、定量化された数値に基づいて、当該スフェアまたは当該スフェアを含む判定対象ロットについて、細胞の未分化率及び生存率を推定することが可能である。また、上記のようにして得られる二次元パワースペクトルを、例えば、細胞の生存率及び未分化率が既知である標準サンプルと比較することによって、当該スフェアまたは当該スフェアを含む判定対象ロットについて、例えば良否判定を行うことが可能である。なお、図17に示す例では、ky=0としたが、R=(kx2+ky2)1/2とし、Rについて図17に示すようなパワースペクトルをとってもよい。また、Wiener-Khintchineの定理より、パワースペクトルの逆フーリエ変換である自己相関関数についても同様のことができる。
1 撮像システム
10 ホログラム光学系
11 レーザ光源
12 ビームスプリッタ
13 コリメートレンズ
14 試料
15 対物レンズ
17 結像レンズ
18 ビームスプリッタ
19 CMOSカメラ
20 光ファイバ
21 コリメートレンズ
34 ダイクロイックミラー
500 コンピュータ
502 主記憶装置
503 補助記憶装置
504 通信インターフェース
505 表示部
506 オートフォーカスプログラム
507 バス
IP 位相差画像
θ 位相差量
θB バックグランドの位相
θS スフェアの存在する領域の位相
θk 1画素の位相差量
k 画素
w カーブの半値幅
λ 平均位相ゆらぎ
Claims (9)
-
複数の細胞の凝集体を撮像したホログラムから前記凝集体の位相差画像を生成し、
前記位相差画像を構成する複数の画素における位相差量の配列のランダム性を示す第1の指標値を導出し、
前記第1の指標値に基づいて、前記凝集体を構成する細胞の状態を判定する判定方法。
-
前記第1の指標値は、前記位相差画像における同一の位相差量の画素を結んだ等位相線によって囲まれた領域の形状の、円形からの乖離の程度に応じて定まる
請求項1に記載の判定方法。
-
前記位相差画像から前記凝集体の形状に依存する成分を除去する処理を施した形状成分除去画像に基づいて前記第1の指標値を導出する
請求項1に記載の判定方法。
-
前記第1の指標値に基づいて、前記凝集体を構成する細胞の生存率に関する判定を行う
請求項1から請求項4のいずれか1項に記載の判定方法。
-
前記凝集体を構成する細胞は幹細胞であり、
前記第1の指標値に基づいて、前記凝集体を構成する幹細胞の未分化率に関する判定を行う
請求項1から請求項5のいずれか1項に記載の判定方法。
-
判定対象ロットに含まれる複数の前記凝集体について、前記第1の指標値と前記凝集体の粒径との相関性を示す第2の指標値を導出し、
前記第2の指標値に基づいて、前記判定対象ロットについて判定を行う
請求項1から請求項6のいずれか1項に記載の判定方法。
-
前記第2の指標値に基づいて、前記判定対象ロットに含まれる細胞の生存率に関する判定を行う
請求項7に記載の判定方法。
-
前記凝集体を構成する細胞は幹細胞であり、
前記第2の指標値に基づいて、前記判定対象ロットに含まれる幹細胞の未分化率に関する判定を行う
請求項7または請求項8に記載の判定方法。
Applications Claiming Priority (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2018185583 | 2018-09-28 | ||
JP2018185583 | 2018-09-28 | ||
PCT/JP2019/031862 WO2020066348A1 (ja) | 2018-09-28 | 2019-08-13 | 判定方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
JPWO2020066348A1 true JPWO2020066348A1 (ja) | 2021-09-16 |
JP7123155B2 JP7123155B2 (ja) | 2022-08-22 |
Family
ID=69949349
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2020548139A Active JP7123155B2 (ja) | 2018-09-28 | 2019-08-13 | 判定方法 |
Country Status (4)
Country | Link |
---|---|
US (1) | US11869181B2 (ja) |
EP (1) | EP3859006A4 (ja) |
JP (1) | JP7123155B2 (ja) |
WO (1) | WO2020066348A1 (ja) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPWO2022224722A1 (ja) * | 2021-04-21 | 2022-10-27 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2007122655A2 (en) * | 2006-04-21 | 2007-11-01 | Consiglio Nazionale Delle Ricerche | Quantitative phase-contrast digital holography method for the numerical reconstruction of images, and relevant apparatus |
JP2009521216A (ja) * | 2005-12-22 | 2009-06-04 | フェイズ ホログラフィック イメージング ペーハーイー アーベー | 細胞サンプル分析のための方法と装置 |
JP2015192644A (ja) * | 2014-03-26 | 2015-11-05 | 株式会社Screenホールディングス | スフェロイドの評価方法およびスフェロイド評価装置 |
US20170261930A1 (en) * | 2011-07-19 | 2017-09-14 | Ovizio Imaging Systems NV/SA | Method and system for detecting and/or classifying cancerous cells in a cell sample |
WO2019176427A1 (ja) * | 2018-03-12 | 2019-09-19 | 富士フイルム株式会社 | 判定方法 |
Family Cites Families (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2004079007A2 (en) | 2003-03-04 | 2004-09-16 | National Institute Of Advanced Industrial Science And Technology | Time-lapse cell analysis method |
IL155435A0 (en) | 2003-04-14 | 2003-11-23 | Bromine Compounds Ltd | Solid biocide formulations |
US9008406B2 (en) | 2009-10-09 | 2015-04-14 | Kawasaki Jukogyo Kabushiki Kaisha | Method and apparatus for discriminating undifferentiated pluripotent stem cells, and automated culture method and system |
US9557549B2 (en) * | 2011-12-09 | 2017-01-31 | Massachusetts Institute Of Technology | Systems and methods for self-referenced quantitative phase microscopy |
JP6015113B2 (ja) | 2012-05-11 | 2016-10-26 | 株式会社ニコン | 細胞評価装置、細胞評価方法およびプログラム |
WO2014041935A1 (ja) | 2012-09-13 | 2014-03-20 | 浜松ホトニクス株式会社 | 多能性幹細胞の分化の度合いを判別する方法 |
US11016443B2 (en) * | 2017-10-20 | 2021-05-25 | Naeilhae, Co. Ltd. | Holographic reconstruction apparatus and method |
US11269294B2 (en) * | 2018-02-15 | 2022-03-08 | University Of Connecticut | Portable common path shearing interferometry-based holographic microscopy system with augmented reality visualization |
-
2019
- 2019-08-13 WO PCT/JP2019/031862 patent/WO2020066348A1/ja unknown
- 2019-08-13 JP JP2020548139A patent/JP7123155B2/ja active Active
- 2019-08-13 EP EP19865460.0A patent/EP3859006A4/en active Pending
-
2021
- 2021-02-23 US US17/182,569 patent/US11869181B2/en active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2009521216A (ja) * | 2005-12-22 | 2009-06-04 | フェイズ ホログラフィック イメージング ペーハーイー アーベー | 細胞サンプル分析のための方法と装置 |
US20100060897A1 (en) * | 2005-12-22 | 2010-03-11 | Phase Holographic Imaging Phi Ab | Method and apparatus for analysis of a sample of cells |
WO2007122655A2 (en) * | 2006-04-21 | 2007-11-01 | Consiglio Nazionale Delle Ricerche | Quantitative phase-contrast digital holography method for the numerical reconstruction of images, and relevant apparatus |
US20170261930A1 (en) * | 2011-07-19 | 2017-09-14 | Ovizio Imaging Systems NV/SA | Method and system for detecting and/or classifying cancerous cells in a cell sample |
JP2015192644A (ja) * | 2014-03-26 | 2015-11-05 | 株式会社Screenホールディングス | スフェロイドの評価方法およびスフェロイド評価装置 |
WO2019176427A1 (ja) * | 2018-03-12 | 2019-09-19 | 富士フイルム株式会社 | 判定方法 |
Non-Patent Citations (1)
Title |
---|
RECICLOUT, JPN6019044368, December 2009 (2009-12-01), pages 10 - 15, ISSN: 0004696971 * |
Also Published As
Publication number | Publication date |
---|---|
WO2020066348A1 (ja) | 2020-04-02 |
US11869181B2 (en) | 2024-01-09 |
EP3859006A1 (en) | 2021-08-04 |
EP3859006A4 (en) | 2021-11-24 |
US20210174507A1 (en) | 2021-06-10 |
JP7123155B2 (ja) | 2022-08-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP6995975B2 (ja) | 判定方法 | |
Luo et al. | Single-shot autofocusing of microscopy images using deep learning | |
JP6490219B2 (ja) | デジタルホログラフィにおけるオートフォーカスシステムおよびオートフォーカス方法 | |
Anand et al. | Automatic identification of malaria-infected RBC with digital holographic microscopy using correlation algorithms | |
US20190101482A1 (en) | Holographic method for characterising a particle in a sample | |
Trusiak et al. | DarkFocus: numerical autofocusing in digital in-line holographic microscopy using variance of computational dark-field gradient | |
JP2019133429A (ja) | 画像処理方法、画像判定方法、プログラムおよび記録媒体 | |
US20230085827A1 (en) | Single-shot autofocusing of microscopy images using deep learning | |
Anand et al. | Digital holographic microscopy for automated 3D cell identification: an overview | |
WO2018158901A1 (ja) | 細胞解析方法及び細胞解析装置 | |
JP7026853B2 (ja) | 画像処理装置、評価システム、画像処理プログラム及び画像処理方法 | |
US11869181B2 (en) | Determination method | |
Pitkäaho et al. | Partially coherent digital in-line holographic microscopy in characterization of a microscopic target | |
CN108593528B (zh) | 基于激光干涉的非球形粗糙粒子形状和尺寸测量方法 | |
JP7046211B2 (ja) | 判定方法 | |
Yao et al. | Automatic three-dimensional imaging for blastomere identification in early-stage embryos based on brightfield microscopy | |
JP7382289B2 (ja) | 画像処理方法、プログラムおよび記録媒体 | |
Castañeda et al. | Concurrent execution of phase compensation and automatic focusing procedures for telecentric off-axis digital holographic microscopy | |
Tai et al. | A forward reconstruction, holographic method to overcome the lens effect during 3D detection of semi-transparent, non-spherical particles | |
Zhao et al. | Adaptive spatial filtering based on fuzzy C-means and phase in off-axis digital holographic microscopy | |
Chang et al. | Combining optical diffraction tomography with imaging flow cytometry for characterizing morphology, hemoglobin content, and membrane deformability of live red blood cells | |
Martínez-Suárez et al. | Modified correlation method to autofocus particles in digital in-line holography | |
Kiss | Proper autofocus for better particle measurements | |
WO2024013202A1 (fr) | Procede et systeme de caracterisation de microorganismes par microscopie holographique numerique | |
Lu et al. | Comparison study of distinguishing cancerous and normal prostate epithelial cells by confocal and polarization diffraction imaging |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20210311 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20220201 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20220317 |
|
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: 20220719 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20220809 |
|
R150 | Certificate of patent or registration of utility model |
Ref document number: 7123155 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |