CN105847633A - 图像处理设备及其图像处理方法 - Google Patents
图像处理设备及其图像处理方法 Download PDFInfo
- Publication number
- CN105847633A CN105847633A CN201610158846.7A CN201610158846A CN105847633A CN 105847633 A CN105847633 A CN 105847633A CN 201610158846 A CN201610158846 A CN 201610158846A CN 105847633 A CN105847633 A CN 105847633A
- Authority
- CN
- China
- Prior art keywords
- motion
- cell
- unit
- evaluation
- image processing
- 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
Classifications
-
- 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
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04N—PICTORIAL COMMUNICATION, e.g. TELEVISION
- H04N5/00—Details of television systems
- H04N5/14—Picture signal circuitry for video frequency region
- H04N5/144—Movement detection
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/20—Analysis of motion
- G06T7/215—Motion-based segmentation
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B11/00—Measuring arrangements characterised by the use of optical techniques
-
- 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/84—Systems specially adapted for particular applications
-
- 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/5014—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 toxicity
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01P—MEASURING LINEAR OR ANGULAR SPEED, ACCELERATION, DECELERATION, OR SHOCK; INDICATING PRESENCE, ABSENCE, OR DIRECTION, OF MOVEMENT
- G01P3/00—Measuring linear or angular speed; Measuring differences of linear or angular speeds
- G01P3/64—Devices characterised by the determination of the time taken to traverse a fixed distance
- G01P3/80—Devices characterised by the determination of the time taken to traverse a fixed distance using auto-correlation or cross-correlation detection means
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T1/00—General purpose image data processing
-
- 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
- G06T7/0014—Biomedical image inspection using an image reference approach
- G06T7/0016—Biomedical image inspection using an image reference approach involving temporal comparison
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/20—Analysis of motion
-
- 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)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Theoretical Computer Science (AREA)
- Computer Vision & Pattern Recognition (AREA)
- General Health & Medical Sciences (AREA)
- Multimedia (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Biomedical Technology (AREA)
- Immunology (AREA)
- Biochemistry (AREA)
- Analytical Chemistry (AREA)
- Biotechnology (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Radiology & Medical Imaging (AREA)
- Quality & Reliability (AREA)
- Medical Informatics (AREA)
- Urology & Nephrology (AREA)
- Organic Chemistry (AREA)
- Pathology (AREA)
- Toxicology (AREA)
- Zoology (AREA)
- Hematology (AREA)
- Microbiology (AREA)
- Molecular Biology (AREA)
- Wood Science & Technology (AREA)
- Medicinal Chemistry (AREA)
- Tropical Medicine & Parasitology (AREA)
- Sustainable Development (AREA)
- Signal Processing (AREA)
- Cell Biology (AREA)
- General Engineering & Computer Science (AREA)
- Genetics & Genomics (AREA)
- Food Science & Technology (AREA)
- Apparatus Associated With Microorganisms And Enzymes (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
- Measuring Or Testing Involving Enzymes Or Micro-Organisms (AREA)
Abstract
本技术涉及一种图像处理设备及其图像处理方法。所述图像处理设备包括:运动量生成单元,被配置为对于包含细胞的一个或多个图像,针对每个子区域生成运动量;以及评价单元,使用由所述运动量生成单元生成的运动量的一部分或全部评价所述细胞的运动传播状态。
Description
技术领域
本公开涉及显示控制设备和方法、图像处理设备和方法以及程序,尤其涉及能够容易并且非侵入性地观察脉动传播的显示控制设备和方法、图像处理设备和方法以及程序。
背景技术
在再生医学领域中,使用由培养细胞制造的经培养的细胞,已经实现了由于事故、疾病等原因丧失的身体细胞、组织、器官等的再生以及功能恢复。可作为这种经培养的细胞制造的细胞组织具有多样性。经培养的细胞的一个实例包括心肌细胞,且这些心肌细胞用于心脏治疗。经培养的心肌细胞本身具有与脉动对应的运动。因此,在经培养的心肌细胞的制造阶段,例如,需要在质量上评价运动是否有利。
在对经培养的心肌细胞进行这种质量评价方面,在当前情况下进行例如目视观察。而且,已经通过电极刺穿经培养的心肌细胞测量电位。然而,目视观察主要取决于观察者的主观视角,并且难以获得客观、精确的评价结果。而且,在测量电位的情况下,经培养的心肌细胞与电极进行接触,从而存在测量具有侵入性的问题。此外,根据电位测量可进行量化的信息限于脉动时间等。而且,待测量的物体限于放置在电极上。
因此,作为过去的技术,已知一种配置,其中,在通过拍摄心肌细胞所获得的成像屏幕内设置测量点,自动测量这些测量点的亮度,并且通过测量值测量心肌细胞的变形周期(例如,见专利文档1)。
顺便提一下,通过分析经培养的心肌细胞的相位差观察运动图像从而获得的不同区域内的脉动示出了通过与培养持续时间相关的方式进行的协同脉动。然而,该脉动示出了各种药物的配给所造成的波动。通过以某种方式检测这种波动,可提前在药物研制中评价药物毒性、影响等,并且近年来这已经受到关注。
例如,在过去已经有一种方法,其中,由设置在培养皿的底部上的电极检测细胞的外部场电位,并且由细胞的膜电位变化,获取细胞的脉动行为。而且,已经有一种方法,其中,将与钙连接并且发射光的荧光染料放入细胞内,并且检测根据细胞的激励(动作电位)波动的钙浓度,从而检测细胞的脉动节奏,并且评价细胞的信息传播图案。
引用列表
专利文献
专利文献1:日本专利申请公开号63-233392(图1)
发明内容
本发明要解决的问题
然而,在培养皿内设置电极并且检测电位的变化时,需要一种特别的培养皿。而且,脉动传播的检测取决于设置在培养皿上的电极的密度。因此,通过现有装置的密度难以检测复杂的传播图案。而且,在放入荧光染料的情况下,荧光染料较为昂贵,放入荧光染料的工作复杂并且费时,并且会造成变色。而且,在这些技术中,将电压施加给待观察的物体,并且放入荧光染料。因此,具有影响待观察对象的可能性。即,这些技术不会容易并且非侵入性地进行脉动传播观察。
鉴于以上内容,已经提出了本公开,并且本公开的目标在于,能够容易并且非侵入性地观察脉动传播。
问题的解决方案
根据本公开的一个方面,提供了一种图像处理设备,包括:运动量生成单元,被配置为对于包含细胞的一个或多个图像,针对每个子区域生成运动量;以及评价单元,被配置为使用由所述运动量生成单元生成的运动量的一部分或所有评价所述细胞的运动传播状态。
评价单元可将由运动量表示的细胞的运动传播状态与确定的细胞的运动为正常的情况的比较结果作为图像进行显示。
通过获得细胞图像的运动量与确定为细胞的运动正常的情况的差分,并评价差分,评价单元对细胞的运动传播状态的变化进行定量评价。
评价单元可对细胞图像内没有运动的静止区域的变化进行评价。
评价单元可对运动量大的区域的位移速度的变化进行评价。
评价单元可对运动量的重心的轨迹的变化进行评价。
评价单元可对有关细胞图像的运动的速度和方向的直方图进行评价。
评价单元可对在几个小时到几天的观察期间细胞的运动传播状态进行评价。
评价单元可评价细胞的运动传播状态,以便无论是否发生QT延长,都检测心肌细胞发生心律失常,QT延长为在心肌细胞的Q波和T波之间的时间延长。
而且,根据本公开的另一个方面,提供了一种图像处理设备的图像处理方法,所述方法包括以下步骤:由运动量生成单元生成包含细胞的一个或多个图像针对每个子区域生成运动量;以及由评价单元使用生成的运动量的一部分或全部评价细胞的运动传播状态。
而且,根据本公开的另一个方面,提供了一种程序,使计算机用作:运动量生成单元,被配置为对于包含细胞的一个或多个图像针对每个子区域生成运动量;以及评价单元,被配置为使用由所述运动量生成单元生成的运动量的一部分或全部评价所述细胞的运动传播状态。
在本公开的一个方面,通过包含细胞的一个或多个图像对细胞图像进行显示控制,使得针对每个子区域生成的运动量的一部分或全部与所述每个子区域相关联,并在所述细胞图像上叠加并显示。
在本公开另一个方面,相对于包含细胞的一个或多个图像,针对每个子区域生成运动量,并且使用生成的运动量的一部分或所有评价细胞的运动传播状态。
本发明的效应
根据本公开,可处理图像。尤其地,可容易并且非侵入性地观察脉动传播。
附图说明
图1为描述运动的协同性的示图;
图2为描述运动的协同性的示图;
图3为示出药物评价设备的主要配置实例的方框图;
图4为由给药造成的对脉动传播的影响的状态的可视化实例的示图;
图5为由给药造成的对脉动传播的影响的状态的另一个可视化实例的示图;
图6为由给药造成的对脉动传播的影响的状态的又一个可视化实例的示图;
图7为示出评价指标数据生成单元的主要配置实例的方框图;
图8为描述待评价物体的图像数据的配置实例的示图;
图9为示出运动检测单元的主要配置实例的方框图;
图10为示出帧图像数据的块划分实例的示图;
图11为描述运动检测数据的配置实例的示图;
图12为示出评价单元的主要配置实例的方框图;
图13为示出运动评价单元的主要配置实例的方框图;
图14为描述运动量重心评价的状态的实例的示图;
图15为描述相关直方图评价的状态的实例的示图;
图16为在图15之后描述相关直方图评价的状态的实例的示图;
图17为在图16之后描述相关直方图评价的状态的实例的示图;
图18为描述评价过程的流程的实例的流程图;
图19为描述评价指标数据生成过程的流程的实例的流程图;
图20为描述影响评价过程的流程的实例的流程图;
图21为描述运动评价过程的流程的实例的流程图;
图22为示出个人计算机的主要配置实例的方框图。
具体实施方式
在后文中,描述用于实现本公开的实施方式(在后文中称为实施方式)。按照以下顺序进行描述。
1、第一实施方式(药物评价设备)
2、第二实施方式(个人计算机)
<1、第一实施方式>
【经培养的心肌细胞】
在再生医学中,已经使用经培养的细胞(这些细胞为由从活体中收集的培养细胞制造的细胞组织)治疗人体的各种组织和器官。在图1A中所示的经培养的细胞1为通过培养和培育心肌细胞所获得的细胞。例如,经培养的心肌细胞为通过培养心肌细胞所获得的培养细胞,可用于进行心脏治疗等。而且,经培养的心肌细胞可用于药物研制过程中对心脏的毒性评价。
在活体中,心肌细胞连续地脉动,同时反复收缩和舒张。因此,在作为经培养的心肌细胞的培养细胞1中,在每个部分内的细胞与在图1B中的运动向量2一样在预定的方向移动,从而整个细胞反复收缩和舒张。实际上,心肌细胞具有自主脉动的部分以及根据周围的脉动而进行脉动的部分。即,可具有培养细胞1根据一个区域不自主脉动的情况。在这种情况下,通过使用电极,将电压从外部定期施加给培养细胞,培养细胞1可脉动。通过以这种方式从外部施加电压进行起搏的培养细胞1的脉动与自主脉动基本上相似。即,使用本技术,培养细胞1自主脉动的情况和培养细胞1通过施加电压进行脉动的情况观察起来相似。
图1C示出了将培养细胞1的观察区域划分成多个子区域(块),针对每个块观察运动量(运动向量),并观察其时移的情况。例如,图1C的曲线图4-1示出了块3-1的运动量的时移,曲线图4-2示出了块3-2的运动量的时移。
图2A的曲线图5-1至5-3示出了在曲线图4-1中所示的块3-1的细胞的运动量与在曲线图4-2中所示的块3-2的细胞的运动量之间的关系的时移。
在开始从活体中收集细胞时,如曲线图5-1中所示,块3-1中具有的细胞的运动量和块3-2中具有的细胞的运动量具有较差的相关性。然而,随着时间的流逝并且随着培养的继续,如图5-2中所示,双方逐步相关,并且随着时间的进一步流逝,如图5-3中所示,双方强烈相关。
即,与在图2B中所示的曲线图一样,在所培育的培养细胞1的多个位置之中的运动量的相关系数稳定地达到较大的值。即,每个区域内的细胞的运动的协同性变得更强。理想地,细胞的活动彼此相关,并且培养细胞1在总体上作为一个活体组织进行脉动。
所培育的培养细胞1由多个细胞构成,并且信号通过间隙连接在这些细胞之间传输。因此,在这些细胞之间,造成运动相关性,并且培养细胞1在总体上作为一个活体组织进行脉动。在这种所培育的培养细胞1中,在这些细胞之间传播脉动。
具有高相关性(协同性)并且传播脉动的培养细胞1j可用于评价药物。例如,为培养细胞1配给药物,观察给药前后的运动,从而可通过培养细胞1的运动的变化,评价药物的效应和毒性等。
【药物评价设备】
图3为示出药物评价设备的主要配置实例的方框图。
在图3中所示的药物评价设备100为通过观察心肌细胞110的运动对为培养的心肌细胞110配给的药物进行评价的设备。如图3中所示,药物评价设备100包括成像单元101、待评价物体的图像数据生成/记录单元102、评价指标数据生成单元103以及评价单元104。
成像单元101将作为待观察物体的培养的心肌细胞110成像。成像单元101可直接(不通过其他部件)将培养的心肌细胞110成像。或者可通过其他部件(例如,显微镜)对培养的心肌细胞110成像。
而且,培养的心肌细胞110可固定或者可不固定至成像单元101。由于药物评价设备100通常检测运动(位置的时间变化),所以,更有利地,培养的心肌细胞110固定至成像单元101。
在预定的周期内,成像单元101将培养的心肌细胞110成像。即,成像单元101获得作为目标的培养的心肌细胞110的运动图像。在给药之前和之后,成像单元101对培养的心肌细胞110成像。要注意的是,例如,在根据预定的条件给药之后,例如,在经过预定的时间之后,成像单元101可对培养的心肌细胞110成像几次。
成像单元101将通过成像获得的培养的心肌细胞110图像的图像信号111(运动图像)提供给待评价的对象的图像数据生成/记录单元102。
待评价的对象的图像数据生成/记录单元102基于从成像单元101提供的图像信号,生成待评价的对象的图像数据,并且例如在内部记录介质内储存和保持所生成的待评价的物体的图像数据。在此处生成的待评价的物体的图像数据变成由例如通过对经培养的心肌细胞110成像获得的图像信号生成的运动图像数据。
例如,待评价的对象的图像数据生成/记录单元102可从由成像单元101提供的多个帧图像中仅仅提取在该周期的一部分中的帧图像,并且可具有作为待评价的对象的图像数据的帧图像。可选地,例如,待评价的对象的图像数据生成/记录单元102可提取由成像单元101提供的每个帧图像的一部分的区域作为小型帧图像,并且可具有由作为待评价的物体的图像数据的小型帧图像构成的运动图像。
而且,例如,待评价的对象的图像数据生成/记录单元102可将任意的图像处理应用于由成像单元101提供的帧图像,并且可具有作为待评价的物体的图像数据的图像处理的结果。作为图像处理,放大、减小、旋转以及变形,可考虑校正亮度和颜色、提供清晰度、去除噪音、产生中间帧图像等。当然,可采用上述以外的任何图像处理。
待评价的对象的图像数据生成/记录单元102可在预定的时间或根据评价指标数据生成单元103的请求,将所储存的待评价的对象的图像数据112提供给评价指标数据生成单元103。
在所提供的待评价的对象的图像数据112的各帧图像内,评价指标数据生成单元103针对每个块对待观察的对象(经培养的心肌细胞110)进行运动检测,每个块为通过将待观察的物体(经培养的心肌细胞110)的图像的所有区域划分为多个区域所获得的子区域。
评价指标数据生成单元103将每个块的经检测的运动表示为运动向量,并且获得运动向量的大小(运动量)。要注意的是,该运动量为绝对值,因此在后文中,该运动量也称为运动量绝对值。
评价指标数据生成单元103将运动量绝对值作为评价指标数据113提供给评价单元104。
评价单元104将作为评价指标数据113提供的运动量绝对值映射到二维彩色图内,并且通过表示脉动传播状态,而且进一步通过在数量上评价其状态,评价药物。评价单元104输出二维彩色图和量化评价结果作为评价值114。
二维彩色图表示在彩色平面上的每个位置处的参数的大小。在这种情况下,该平面表示经培养的心肌细胞110的观察区域,并且该颜色表示运动量绝对值(该值的大小)。即,运动量绝对值的位置由在该平面上的位置表示,并且该值的大小由颜色表示。即,运动量绝对值在观察区域内的分布由颜色在该平面上的变化表示。
在图4到图6中示出了这种二维彩色图的实例。在图4中所示的八个彩色图分别表示运动量绝对值在观察区域内的分布,在彩色图上的位置表示在观察区域内的位置,并且更高的颜色密度表示更大的运动量绝对值。即,颜色密度较高的部分表示该部分强烈脉动(进行收缩或舒张)。要注意的是,为了进行描述,在彩色图上绘制箭头,并且并未在实际的彩色图上绘制这些箭头。
在八个彩色图中左边的四个彩色图表示在配给有机溶剂(控制)之后的脉动传播状态的实例。有机溶剂基本上不影响经培养的心肌细胞110的脉动。有机溶剂的实例包括二甲亚砜。这四个彩色图按照从上面开始的顺序分别表示紧接在给药之后(0ms)、在给药40ms之后、在给药80ms之后以及在给药120ms之后的脉动状态。
而且,在右边的四个彩色图表示在配给1-庚醇之后的脉动传播状态的实例。众所周知,1-庚醇阻碍了在细胞之间中继信号传输的间隙连接的功能。这四个彩色图按照从上面开始的顺序分别表示紧接在给药之后(0ms)、在给药80ms之后、在给药160ms之后以及在给药240ms之后的脉动状态。
在图4的左边的情况下(在仅仅配给有机溶剂的情况下),如箭头所示,随着时间的流逝,脉动(颜色密度较高的部分)使其面积扩大,同时在从右上侧到左下侧的方向大致线性传播,从而远离初始位置。这表示间隙连接正常运转并且顺利进行传播的状态。
相比之下,在图4的右边的情况下(在配给1-庚醇的实例的情况下),如箭头所示,随着时间的流逝,脉动(颜色密度较高的部分)线性传播,同时改变旋转的方向。即,脉动有时在返回初始位置(向后移动)的方向进行。此外,与左边相比,脉动需要更长的时间进行传播(即,传播的速度较慢)。这表示间隙连接的功能受碍并且未顺利进行传播的状态。
图5的彩色图比较在配给有机溶剂(控制)之后和在配给18-β-甘草次酸之后之间的脉动传播状态。众所周知,18-β-甘草次酸阻碍间隙连接,与1-庚醇一样。在图5的右边的这五个彩色图按照从上面开始的顺序分别表示紧接在给药之后(0ms)、在给药80ms之后、在给药120ms之后、在给药160ms之后以及在给药240ms之后的脉动状态。
在图5的左边的情况下(在仅仅配给有机溶剂的实例的情况下),如箭头所示,随着时间的流逝,脉动(颜色密度较高的部分)在从上侧到下侧的方向大致线性传播,从而远离初始位置,与图4的左侧的情况一样。这表示间隙连接正常运转并且顺利进行传播的状态。
相反,在图5的右边的情况下(在配给18-β-甘草次酸的实例的情况下),如箭头所示,继续在不同的方向脉动(颜色密度较高的部分),例如,从右侧到左侧、从下侧到上侧、从左侧到右侧等。而且,在传播的中间,可具有脉动传播临时停止(在待观察的物体区域内未造成脉动传播)的情况。此外,与左边相比,脉动需要更长的时间进行传播(即,传播的速度较慢)。这表示间隙连接的功能受阻并且未顺利进行传播的状态。
图6的彩色图比较在配给有机溶剂(控制)之后和在配给DL-索他洛尔之后的脉动传播状态。众所周知,DL-索他洛尔阻碍钾通道。图6的右边的这六个彩色图按照从上面开始的顺序分别表示紧接在给药之后(0ms)、在给药40ms之后、在给药80ms之后、在给药120ms之后、在给药160ms之后以及在给药200ms之后的脉动状态。
在图6的左边的情况下(在仅仅配给有机溶剂的实例的情况下),如箭头所示,随着时间的流逝,脉动(颜色密度较高的部分)在从左上侧到右下侧的方向大致线性传播,从而远离初始位置,与图4和图5的左侧的情况一样。这表示在每个细胞内的脉动稳定并且顺利进行传播的状态。
相比之下,在图6的右边的情况下(在配给DL-索他洛尔的实例的情况下),如箭头所示,脉动(颜色密度较高的部分)从左上侧传播到右下侧,但是该传播在观察区域内停止。然后,脉动开始从左下侧传播到右上侧。要注意的是,在观察周期内,在右侧内的一部分观察区域中未造成脉动。即,该脉动不在该区域的这个右侧部分内传播。此外,该传播速度比左边低。
在配给DL-索他洛尔时,由于在舒张过程中起作用的钾通道功能的变化,所以舒张过程变化。因此,每个细胞的脉动的波形变化,并且脉动并未与图4到图6的左侧的实例一样顺利地传播。
评价单元104通过这种二维彩色图表示各块的运动量绝对值。评价单元104在预定的帧间距内生成这种二维彩色图(可在每个帧内或者在多个帧内),如图4到6中所示。即,评价单元104通过二维彩色图的时间变化表示脉动传播。
评价单元104将二维彩色图作为图像向用户呈现。即,药物评价设备100可容易并且非侵入性地观察脉动传播状态,并且可呈现观察结果。用户可通过在二维彩色图的图像内显示的脉动传播的状态,容易并且非侵入性地评价所配给的药物的影响(效应、毒性等)。
要注意的是,评价单元104可使用二维彩色图量化评价脉动传播的状态。在图4的实例的情况下,在右侧图中所示的给药之后的传播状态(线性传播状态)与在左侧图中所示的给药之后的脉动传播状态(在观察区域内旋转的传播状态)不同。评价单元104比较例如在给药之前和之后的脉动传播的这种状态,并且可根据差分值的大小,确定在给药之后的脉动传播状态是否正常。
而且,评价单元104也可通过检测传播方向的变化或该方向(例如,反方向等)确定在给药之后的脉动传播状态是否正常。而且,评价单元104也可通过确定传播速度是否足够快来确定在给药之后的脉动传播状态是否正常。
评价单元104可输出这种评价的结果作为评价值114。即,药物评价设备100可容易并且非侵入性地观察脉动传播状态,并且可容易并且非侵入性地评价待评价的物体(给药造成的影响)。
要注意的是,药物评价设备100可将经培养的心肌细胞110以外的其他东西用作待观察的物体。例如,心肌细胞以外的细胞可用作待观察的对象,或者细胞以外的东西可用作待观察的对象。要注意的是,期望待观察的对象本身移动并且可通过评价该运动来评价为待观察的对象配给的药物。要注意的是,该运动可为自主(自发(主动))运动,或者可为从外部提供的电信号等的被动运动。
【评价指标数据生成单元】
图7为示出图1的评价指标数据生成单元103的主要配置实例的方框图。如图7中所示,评价指标数据生成单元103包括运动检测单元121、运动量绝对值计算单元122以及运动量绝对值储存单元123。
运动检测单元121从待评价的对象的图像数据生成/记录单元102输入所记录的待评价的对象的图像数据112,针对每个模块进行运动检测,并且将检测结果(运动向量)作为运动检测数据提供给运动量绝对值计算单元122。
运动量绝对值计算单元122计算运动量绝对值,该运动量绝对值为每个所提供的运动检测数据(运动向量)的大小。运动量绝对值计算单元122在运动量绝对值储存单元123内提供和储存所计算的运动量绝对值。
运动量绝对值储存单元123在预定的时间或根据评价单元104的请求,将所储存的运动量绝对值作为评价指标数据113提供给评价单元104。
要注意的是,运动检测单元121通过运动量绝对值储存单元123对待评价的对象的图像数据的各帧图像进行处理。
【待评价的对象的图像数据的结构】
图8示出了提供给评价指标数据生成单元103的待评价对象的图像数据112的结构实例。在具有预定长度的评价部分内进行成像(例如,T+1帧(T为任意自然数))。即,例如,提供给评价指标数据生成单元103的待评价对象的图像数据112由与评价部分对应的第一到第(T+1)帧图像数据132-1到132-(T+1)构成。
【运动检测单元的配置实例】
图9为示出运动检测单元121的主要配置实例的方框图。如图9中所示,运动检测单元121包括帧存储器141和运动向量计算单元142。帧存储器141将在每一个帧周期内按照顺序输入的帧图像数据132作为待评价对象的图像数据112进行保持。
运动向量计算单元142输入作为当前时间的待评价对象的图像数据112进行输入的帧图像数据和保持在帧存储器141内的前一个时间(在时间上的前一个时间)的帧图像数据。然后,运动向量计算单元142为每个块计算运动向量,该运动向量表示在这两个帧图像数据之间的运动。将所计算出的运动向量作为运动检测数据151提供给运动量绝对值计算单元122。
更详细地描述由图9的运动检测单元121执行的处理。运动向量计算单元142输入当前时间的帧图像数据132以及前一个时间(在时间上的前一个时间)的帧图像数据132。运动向量计算单元142将这些输入帧图像数据132分成M×N(M和N为任意自然数)个块161,如图10中所示,通过在帧图像之间进行块匹配等技术进行运动检测,并且为每个块161生成运动向量。各模块161由例如(16×16)个像素构成。
运动向量计算单元142使用第一到第(T+1)个帧图像数据132按照顺序执行运动检测过程。即,运动向量计算单元142使用(T+1)个帧图像数据生成(M×N×T)运动检测数据(运动向量)。运动向量计算单元142将通过这种方式计算出的运动向量作为运动检测数据提供给运动量绝对值计算单元122。
在完成使用第T到第(T+1)个帧图像数据132进行的最后的运动检测过程时,将由T个帧单元运动检测数据171-1到171-T构成的运动检测数据提供给运动量绝对值计算单元122,如图11中所示。
通过在每个帧周期内获得的当前时间的帧图像数据132和前一个时间(在时间上的前一个时间)的帧图像数据132上进行运动检测处理,获得帧单元运动检测数据171-1到171-T中的每一个。
例如,通过分别将第四个帧图像数据132-4和第三个帧图像数据132-3作为当前时间的和前一个时间的帧图像数据进行输入,进行运动检测,从而获得第三个帧单元运动检测数据171-3。
而且,通过(M×N)个块单元运动检测数据181形成帧单元运动检测数据171-1到171-T中的每一个。每个块单元运动检测数据181对应一个块161,并且用作表示针对相应块161检测的运动向量的数据。
通过这种方式,本实施方式的运动检测数据151具有对于每个帧单元运动检测数据171包括(M×N)个块单元运动检测数据181的结构。
【评价单元】
图12为示出评价单元104的主要配置实例的示图。如图12中所示,评价单元104包括运动量绝对值获取单元201、映射单元202、随时间变化分析单元203、运动评价单元204、显示单元205以及输出单元206。
运动量绝对值获取单元201从评价指标数据生成单元103的运动量绝对值储存单元123中获取作为评价指标数据113的期望的待评价对象的图像数据112(例如,由用户规定为待观察对象)的运动量绝对值。运动量绝对值获取单元201将所获取的运动量绝对值提供给映射单元202。
映射单元202根据块的坐标,将提供的运动量绝对值映射在平面上,并且生成如图4中所示的二维彩色图。映射单元202为多个帧(待评价的对象的图像数据112的所有或一部分帧)生成上述二维彩色图。
映射单元202将所生成的二维彩色图提供给显示单元205并且使显示单元205显示其图像,将二维彩色图提供给输出单元206并且使输出单元206将数据输出到药物评价设备100的外部(输出到其他设备等)。而且,映射单元202也将所生成的二维彩色图提供给随时间变化分析单元203。
要注意的是,映射单元202可仅仅生成能够表示运动量绝对值的分布的信息,并且所生成的信息不限于二维彩色图。例如,可采用三维或多维彩色图。而且,可采用灰度(密度或模式图案)图代替彩色图。而且,映射单元202可在弯曲的表面上映射运动量绝对值。
随时间变化分析单元203分析每个所提供的二维彩色图(运动量绝对值的分布)的时间变化(随着时间变化)的状态。例如,随时间变化分析单元203在二维彩色图内设置具有预定特征的区域(或点),作为关注区域(或关注点),并且获得其随着时间发生的变化。随时间变化分析单元203将随着时间变化(例如,关注区域的尺寸、形状、位置等随时间变化)的二维彩色图和分析结果提供给运动评价单元204。
运动评价单元204基于所提供的信息评价运动量分布。例如,运动评价单元204评价关注区域的运动等。例如,运动评价单元204确定在给药之前和之后,关注区域的路径是否发生变化,确定关注区域的前进方向是否发生极端变化,并且确定关注区域的前进速度是否发生变化。运动评价单元204将这种评价的结果的数据提供给显示单元205并且促使显示单元205显示数据,将该数据提供给输出单元206并且使输出单元206将数据输出到药物评价设备100的外部(输出给其他设备等)。
显示单元205包括任意的显示装置,并且使该显示装置显示由成像单元101获得的图像信号111(细胞图像)。而且,显示单元205产生由映射单元202提供的二维彩色图的图像,并且使该显示装置显示该图像。例如,显示单元205叠加或设置细胞图像(运动图像)的每帧图像和与帧图像对应的二维彩色图的图像,并且使显示装置显示该图像。而且,显示单元205产生由运动评价单元204提供的评价结果的图像,并且促使显示装置显示该图像。要注意的是,评价结果的图像可为静止图像或者可为运动图像。而且,显示单元205能够显示评价结果的图像以及细胞图像和二维彩色图中的一个或两个图像。在这种情况下,显示单元205可在细胞图像和二维彩色图中的一个或两个图像上叠加评价结果的图像,并且显示这些图像,或者可设置和显示这些图像。而且,在评价结果作为运动图像进行显示的情况下,显示单元205可显示评价结果的帧图像以及与细胞图像或二维彩色图的运动图像对应的帧图像。
输出单元206包括任意的输出接口,并且通过该输出接口将从映射单元202提供的二维彩色图的数据输出到位于药物评价设备100的外部的装置、网络等。而且,输出单元206通过该输出接口将从运动评价单元204提供的评价结果的数据输出给位于药物评价设备100的外部装置、网络等。
【运动评价单元】
图13为示出图12的运动评价单元204的主要配置实例的方框图。如图13中所示,运动评价单元204包括静止区域评价单元211、脉动传播速度评价单元212、运动量重心评价单元213以及相关直方图评价单元。
静止区域评价单元211评价未传播脉动的区域(即,未脉动的区域)的数量和面积或其时间变化。
在给药之后,具有在观察区域(静止区域)内具有未明确地移动(未脉动)的区域的情况。例如,在图6的情况中,在配给DL-索他洛尔之后,脉动未传播到观察区域的右侧部分。即,脉动停止。
DL-索他洛尔阻碍钾通道。在为经培养的心肌细胞配给DL-索他洛尔时,由于在舒张过程中起作用的钾通道功能发生变化,所以存在舒张过程发生变化的情况,从而脉动时间(例如,动作电位维持时间)可延长。而且,在某些情况下,脉动本身可停止。通过这种方式,由于药物的毒性,所以细胞的运动可停止。
静止区域(或其时间变化)的数量或面积可与配给的药物的毒性相关。例如,可具有以下情况:在给药之后,可评价静止区域的数量越大或者静止区域的面积越大,那么药物的毒性就越强。而且,通过从紧接在给药之后到静止区域增加时的时间、从静止区域增加时到静止区域减少时的时间等,可评价药物毒性的强度。
因此,静止区域评价单元211例如将静止区域用作关注区域,比较在给药前后在二维彩色图内的关注区域,并且通过其差别(变化量)的大小,评价静止区域的数量、尺寸、位置、形状等。当然,静止区域评价单元211可在给药后多个时刻内,在二维彩色图内比较关注区域。通过静止区域的这种评价,静止区域评价单元211评价配给的药物的毒性和影响。
脉动传播速度评价单元212在观察区域内评价脉动的传播速度或其时间变化。
如果由于配给的药物的毒性阻碍脉动传播,那么也可影响脉动的传播速度。例如,如果每个细胞的脉动时间延长,那么脉动传播速度可变化。即,脉动传播速度的变化可与配给的药物的毒性相关。例如,可具有以下情况:在给药之后,可评价,脉动传播速度的减小余地越大,那么药物的毒性就越强。而且,通过从紧接在给药之后到传播速度减小时的时间、从传播速度减小时到传播速度恢复时的时间等,可评价药物毒性的强度。
因此,脉动传播速度评价单元212例如将大幅脉动的部分用作关注区域,在给药之前和之后在二维彩色图内比较关注区域的运动(速度),并且通过关注区域的运动的速度的差值(在预定时间内的位移量),评价脉动传播速度的变化。当然,脉动传播速度评价单元212可在给药之后在多个时刻在二维彩色图内比较关注区域的运动。通过脉动传播速度的这种评价,脉动传播速度评价单元212评价所配给的药物的毒性和影响。
运动量重心评价单元213在观察区域内评价运动量的重心的位置和轨迹。
通常,在脉动传播时,在观察区域内的每个部分的运动改变,因此,运动量的重心的位置在观察区域内改变。即,可以运动量的重心位置的轨迹表示脉动传播。换言之,在由于所配给的药物的毒性,脉动传播的状态(路径、速度等)改变时,也影响运动量的重心位置改变的方式。
例如,如图4到图6所示,在给药之前在一个方向线性传播的脉动可传播以便进行旋转,或者可在给药之后在多个方向传播。例如,在给药造成产生静止区域时,脉动传播可缠绕在静止区域的外围周围(旋转)。而且,即使停止传播,未产生静止区域,在偏置在区域之间的细胞的脉动周期时,传播方向也可根据该偏置而改变(旋转)。
而且,可在静止区域内将脉动传播分成多个方向等。而且,在心肌细胞内,一部分细胞发挥起搏的作用,并且根据细胞的脉动,对其他细胞进行脉动传播。由于给药的影响,进行起搏的细胞可代替其他细胞。在某些情况下,进行起搏的细胞的数据可增大/减少。在这种情况下,脉动传播的路径大幅改变。
而且,可在传播的中间停止传播。例如,脉动传播在静止区域内停止。而且,传播速度可改变。例如,由于每个细胞的脉动时间延长,所以传播速度可改变。
由于如上所述的传播状态的变化,所以脉动可在彼此不同的时间从多个方向传播到相同的区域。在这种情况下,细胞的脉动混乱,并且可局部造成重新返回(这可造成心律失常)。
由于脉动传播的这种变化,所以运动量的重心的轨迹改变。因此,如图14A中所示,运动量重心评价单元213将在二维彩色图上的所有区域划分为垂直方向上N个并且水平方向上M个小区域,并且针对每个小区域(m、n)获得运动量v(m、n)。运动量重心评价单元213与以下公式(1)和(2)一样计算运动重心的坐标(Gx、Gy)。
【数学公式1】
【数学公式2】
在已经通过这种方式为每个二维彩色图获得运动量的重心的坐标时,如图14B中所示,运动量重心评价单元213在二维图上绘制每个二维彩色图(每时刻)的重心的位置,并且绘出变化(轨迹)。运动量重心评价单元213比较在给药之前和之后重心的这种轨迹,并且评价轨迹的变化。当然,运动量重心评价单元213可在给药之后在多个时刻内比较重心的轨迹。通过运动量的重心的轨迹的这种评价,运动量重心评价单元213评价所配给的药物的毒性和影响。
要注意的是,多个重心可用于评价轨迹。而且,可一起评价重心的运动量(脉动大小)。例如,运动量重心评价单元213可在三维图上绘制运动量的重心,该重心表示在XY坐标内的位置以及在Z坐标内的运动量。可使用任何方法计算重心的运动量。例如,可采用在重心位置获得的运动量作为重心的运动量,或者可采用重心附近的运动量的平均值作为运动量。可选地,可使用整个观察区域的运动量计算重心的运动量。
相关直方图评价单元214评价在观察区域内运动的分布(速度、方向等)。
心肌细胞在给药之前在正常的状态下按照预定的节奏反复进行脉动。因此,在观察区域内运动的分布(速度、方向等)变为大概恒定。在给药之后,在脉动传播的状态改变时,如上所述,由于所配给的药物的毒性,所以运动的分布也改变。因此,通过获得给药之前和之后运动的速度和方向的分布(直方图),并且评价分布的变化,相关直方图评价单元214评价所配给的药物的毒性和影响。
描述更具体的实例。例如,相关直方图评价单元214获得在相邻的小区域之间的脉动的相关系数d,如图15中所示。如图15的中间所示,该相关系数d为表示脉动的时间间隙的参数,并且如图15的右边所示,该间隙越大,那么值就越大。即,在这两个小区域的脉动时间改变半个周期的状态中,该相关系数d最大。在这两个小区域的脉动时间一致时,相关系数d的值为0。
在要处理的,相关直方图评价单元214针对待处理的小区域的左边、右边、顶部以及底部的每个相邻的小区域计算这种相关系数d,如图16A中所示(d0到d3)。相关直方图评价单元214与以下公式(3)和(4)一样使用相关系数d(d0到d3)计算水平方向上的运动量vx和垂直方向上的运动量vy。
【数学公式3】
【数学公式4】
相关直方图评价单元214由运动量vx和vy获得待处理的小区域(m、n)的运动量v(m、n)和其方向θ(m、n),如图16B中所示。相关直方图评价单元214通过这种方式为每个小区域获得运动量及其方向,并且生成如图17中所示的运动量(速度)和方向的直方图(频率分布)。
图17A示出了在给药之前的运动的速度和方向的直方图,并且图17B示出了在给药之后的运动的速度和方向的直方图。如图17A和17B中所示,在给药之前和之后的分布不同时,可确定给药造成影响。即,相关直方图评价单元214通过在给药之前和之后获得有关运动的速度和方向的分布(直方图),并且评价分布的变化,评价所配给的药物的毒性和影响。
通过相关直方图评价单元214的静止区域评价单元211分别将所获得的评价结果提供给评价结果数据生成单元215。评价结果数据生成单元215适当地总结评价结果,并且将这些结果作为评价结果数据提供给显示单元205和输出单元206(图12)。
要注意的是,运动评价单元204的上述配置为一个实例。运动评价单元204可具有任何配置,只要与待观察的对象的运动相关,就可评价任何指标。
如上所述,药物评价设备100可容易并且非侵入性地观察经培养的心肌细胞110的脉动传播状态,并且可容易并且非侵入性地评价待评价的对象(所配给的药物造成的影响)。
【评价过程的流程】
接下来,参照图18的流程图,描述由药物评价设备100执行的评价过程的流程的一个实例。
在开始评价过程时,在步骤S101中,药物评价设备100的成像单元101使作为待观察的对象的经培养的心肌细胞110成像。在步骤S102中,待评价的对象的图像数据生成/记录单元102由通过在步骤S101中成像从而获得的图像信号111生成待评价的对象的图像数据112。
在步骤S103中,评价指标数据生成单元103使用在步骤S102中生成的待评价的对象的图像数据112进行运动检测,计算运动量绝对值,并且生成评价指标数据113。在步骤S104中,评价单元104使用在步骤S103中生成的评价指标数据113生成二维彩色图,并且计算评价值114。
在步骤S105中,评价单元104输出在步骤S104中计算出的评家值114,并且终止评价过程。
【评价指标数据生成过程的流程】
接下来,参照图19的流程图,描述在图18的步骤S103中执行的评价指标数据生成过程的流程的一个实例。
在开始评价指标数据生成过程时,在步骤S121中,评价指标数据生成单元103的运动检测单元121针对每个块检测待评价的对象的运动,并且生成运动向量。在步骤S122中,运动量绝对值计算单元122计算在步骤S121中生成的每个块的运动向量的运动量绝对值。
在步骤S123中,运动量绝对值储存单元123储存在步骤S122中计算初的运动量绝对值。
在步骤S124中,运动检测单元121确定是否已经处理预定周期(评价部分)的数据。在确定在预定的评价部分内具有还未进行运动检测的帧图像时,运动检测单元121使过程返回到步骤S121中,并且对待处理的新的帧图像重复进行运动检测。
而且,在步骤S124中,确定假设在预定的评价部分内进行处理的所有帧图像已经进行运动检测,评价指标数据生成过程终止,并且该过程返回到图18并且继续进行步骤S104。
【相关评价过程的流程】
接下来,参照图20的流程图,描述在图18的步骤S104中执行的影响评价过程的流程的实例。
在开始影响评价过程时,在步骤S141中,运动量绝对值获取单元201从运动量绝对值储存单元123获取运动量绝对值。
在步骤S142中,映射单元202将在步骤S141中获取的运动量绝对值映射在平面上,并且产生二维彩色图。
在步骤S143中,在步骤S142中产生的运动量绝对值的每个二维彩色图中,随时间变化分析单元203分析二维彩色图的随时间变化。
在步骤S144中,运动评价单元204评价在步骤S143中分析的二维彩色图的时间变化(在二维彩色图之中的运动)。即,运动评价单元204评价给药的影响。运动评价单元204促使显示单元205显示评价值114,并且通过输出单元206将该值输出到药物评价设备100的外部。
在完成步骤S144中的过程时,运动评价单元204终止影响评价过程并且使该过程返回图18中。
【运动评价过程的流程】
接下来,参照图21的流程图,描述在图20的步骤S144中执行的运动评价过程的流程的实例。
在开始运动评价过程时,在步骤S161中,静止区域评价单元211评价静止区域。在步骤S162中,脉动传播速度评价单元212评价脉动传播速度。在步骤S163中,运动量重心评价单元213评价运动量的重心。在步骤S164中,相关直方图评价单元214通过直方图评价脉动传播的速度和方向。在步骤S165中,评价结果数据生成单元215生成评价结果数据,用于输出或者用于显示。
该过程在步骤S165完成时,在步骤S165中,评价结果数据生成单元215终止运动评价过程并且使该过程返回图20中。
如上所述,通过进行各种处理,药物评价设备100可容易并且非侵入性地观察待观察的对象,并且可容易并且非侵入性地评价给药对经培养的心肌细胞110的影响。
即,在本技术中,不需要特别的培养皿或荧光试剂。因此,可通过简单、非侵入性以及比较便宜的方式获取细胞脉动行为的变化,从而可容易并且精确地评价药物的毒性等。而且,有利于进行自动化。
要注意的是,药物通常具有两种毒性:在短期内显示影响的短期毒性,例如,在给药之后的大约几秒钟到几分钟,以及在经过长期之后显示影响的长期毒性,例如,在给药之后的几个小时到几天。在使用荧光试剂或电极时,影响待观察的对象。因此,不利于进行长期毒性观察。相反,本技术可非侵入性地观察待观察的对象。而且,在本技术中,可通过与短期毒性观察相似的方法,进行长期毒性观察。因此,本技术不仅有利于进行短期毒性观察,而且有利于进行长期毒性观察。
而且,本技术可通过彼此相似的方法进行短期毒性观察和长期毒性观察。因此,可进行长期毒性观察,例如,紧接在给药之后需要几天进行的观察。即,本技术也可观察和评价毒性随着时间的变化。
要注意的是,在作为待观察的对象的培养细胞通过密集的方式生长时,通过荧光试剂通常难以进行染色。然而,在本技术中,可非侵入性地观察待观察的对象,因此,不论培养细胞的生长程度如何,都可稳定地进行观察和评价。
顺便提一下,心肌细胞的脉动由收缩和舒张构成。通常,心肌的舒张与心电图中的T波对应,并且与心肌细胞膜的复极对应。这个T波延长通常称为QT延长,作为在Q波和T波之间的时间的延长,并且在显示该症状时,指出心律失常的可能性。例如,由于为经培养的心肌细胞配给药物而在阻碍进入钾通道中的离子的吸收和取出时,发生这种QT延长。例如,众所周知,DL-索他洛尔阻碍钾通道。即,在为经培养的心肌细胞配给DL-索他洛尔时,由于在舒张过程起作用的钾通道功能的变化,所以舒张过程改变。
然而,实际上,甚至在发生QT延长的情况下,在整个心肌细胞内发生大约均匀的QT延长时,脉动在这些细胞之中未产生大量间隙,因此,可具有未发生心律失常的情况。相反,如果如上所述,间隙连接受阻,那么即使未造成QT延长,也可产生心律失常。在本技术中,可容易并且非侵入性地观察脉冲的传播。因此,无论是否发生这种QT延长,都可观察到发生心律失常。
而且,在本技术中,由于观察区域可为相对窄的区域,例如,大约0.6平方毫米,所以通过少数细胞和少量试剂,可进行试验。而且,通过市场上可买到的高密度的培养板(1536孔板(每个孔的直径为1.7mm)或384孔板(每个孔的直径为3.6mm)),可充分地评价该试验,并且本技术有利于在药物的研制过程中进行初步筛选。要注意的是,在本技术中,无论观察的面积多大,都可通过相似的方法进行观察。因此,可容易地改变观察区。
而且,目前进行评价药物的毒性的方法,并且可能在未来提出一种新的评价方法或评价标准。由于可非侵入性地观察待观察对象的状态,所以本技术可应用于更多的不同评价方法或评价标准。
而且,本技术就可应用于评价任何物体的任何情况,只要通过观察经培养的心肌细胞110,可对该物体进行评价。例如,该物体可为气体、液体或固体。而且,在观察时,可为环境条件(例如,温度、湿度、大气压力、亮度、振动、磁场等)。
要注意的是,由于本技术可容易地观察脉动传播,所以本技术可应用于以下方法中:通过放入荧光染料并且检测根据细胞的激励(活动电位)进行波动的钙浓度,检测细胞的脉动节奏,并且评价细胞的信息传播模式。
<2、第二实施方式>
【个人计算机】
通过硬件或者通过软件,可执行上述一系列过程。在这种情况下,例如,可配置如图22中所示的个人计算机。
在图22中,根据储存在ROM(只读存储器)302内的程序或者从储存单元313装载到RAM(随机存取存储器)303上的程序,个人计算机300的CPU(中央处理单元)301执行各种处理。RAM 303适当地储存CPU 301执行各种过程所需要的数据等。
CPU 301、ROM 302以及RAM 303通过总线304彼此连接。输入/输出接口310也连接至总线304。
由键盘、鼠标等构成的输入单元311、由(由阴极射线管(CRT)、液晶显示器(LCD)等构成的)显示器和扬声器构成的输出单元312、由硬盘等配置的储存单元313、以及由调制解调器等配置的通信单元314连接至输入/输出接口310。通信单元314通过包括互联网的网络执行通信过程。
必要时,驱动器315还连接至输入/输出接口310,适当地安装可移动介质321(例如,磁盘、光盘、磁光盘或半导体存储器),并且必要时,在储存单元313上安装从其中读出的计算机程序。
在由软件执行上述一系列过程的情况下,通过网络或通过记录介质安装构成软件的程序。
例如,如图22中所示,记录介质与设备主体分开,记录介质不仅由可移动介质321配置,而且由ROM 302配置。分布该可移动介质,以将程序传送给用户,并且该可移动介质由在其上储存程序的磁盘(包括软盘)、光盘(包括光盘只读存储器(CD-ROM)和数字通用光盘(DVD))、磁光盘(包括微型光盘(MD))、半导体存储器等构成。要传送给用户的程序在组装到设备主体、包含在储存单元313内的硬盘等中的状态下记录在该ROM上。
要注意的是,由计算机执行的程序可为根据说明书中描述的顺序的时间序列处理的程序,或者可为并行或在必要时(例如,在调用时)处理的程序。
此外,在该说明书中,描述在记录介质上记录的程序的步骤不仅包括根据所描述的顺序在时间序列中执行的过程,而且包括不必按时间序列执行而是并行执行或单独执行的过程。
此外,在该说明书中,系统表示由多个装置配置的整个配置。
此外,在以上描述中,作为单个装置(或者作为处理单元)进行描述的配置可划分成多个装置(或处理单元)并且配置有这些装置或单元。与此相反,在以上描述中,作为多个装置(或者处理单元)进行描述的配置可集中并且配置为单个装置(或者处理单元)。此外,上述配置以外的配置也可加入这些装置(或者处理单元)的配置中。而且,在整个系统的配置或操作基本上相同时,某个装置(或者处理单元)的一部分配置可包含在另一个装置(或者另一个处理单元)的配置内。即,本技术不受到上述实施方式的限制,并且可在不背离本技术的主旨的范围内进行各种修改。
要注意的是,本技术可使用以下配置。
(1)一种显示控制设备,包括:
第一显示控制单元,被配置为通过包含细胞的一个或多个图像对细胞图像进行显示控制;以及
第二显示控制单元,被配置为以以下方式进行显示控制:将通过所述一个或多个图像针对每个子区域生成的运动量的一部分或全部与所述每个子区域相关联,并在所述细胞图像上叠加和显示。
(2)根据(1)所述的显示控制设备,其中
所述细胞为被动或主动移动的细胞,以及
所述第二显示控制单元显示以表现所述细胞的运动传播状态的方式显示所述运动量的一部分或全部。
(3)根据(2)所述的显示控制设备,其中,所述第二显示控制单元显示所述细胞的运动传播状态,作为在细胞图像内的运动量的分布随时间发生的变化。
(4)根据(3)所述的显示控制设备,其中,所述第二显示控制单元通过对多个时间中的每个时间在曲面或平面上显示表示所述运动量的分布的二维图,表现运动量的分布随时间发生的变化。
(5)根据(4)所述的显示控制设备,其中,所述第二显示控制单元将在细胞图像内的运动量的位置表现为在所述二维图上的位置,并且将运动量的绝对值表现为颜色、密度或图案。
(6)根据(1)到(5)中任一项所述的显示控制设备,其中
所述细胞图像为运动图像,以及
所述第二显示控制单元在所述细胞图像的各帧图像上叠加与所述帧图像对应的运动量的一部分或全部,并显示叠加的图像。
(7)根据(1)到(6)中任一项所述的显示控制设备,其中,所述第二显示控制单元显示运动量以及运动量的分布的评价结果。
(8)一种显示控制设备的显示控制方法,所述方法包括以下步骤:
由第一显示控制单元通过包含细胞的一个或多个图像对细胞图像进行显示控制;以及
由第二显示控制单元以下方式进行显示控制:将通过所述一个或多个图像针对每个子区域生成的运动量的一部分或全部与所述每个子区域相关联,并在所述细胞图像上叠加和显示。
(9)一种程序,其促使计算机用作:
第一显示控制单元,被配置为通过包含细胞的一个或多个图像对细胞图像进行显示控制;以及
第二显示控制单元,被配置为以以下方式进行显示控制:将通过所述一个或多个图像针对每个子区域生成的运动量的一部分或全部与所述每个子区域相关联,并在所述细胞图像上叠加和显示。
(10)一种图像处理设备,包括:
运动量生成单元,被配置为对于包括细胞的一个或多个图像,针对每个子区域生成运动量;以及
评价单元,被配置为使用由所述运动量生成单元生成的运动量的一部分或全部评价细胞的运动传播状态。
(11)根据(10)所述的图像处理设备,其中,所述评价单元将由所述运动量表示的所述细胞的运动传播状态与确定的所述细胞的运动为正常的情况的比较结果作为图像进行显示。
(12)根据(10)或(11)所述的图像处理设备,其中,通过获得所述细胞图像的运动量与确定为所述细胞的运动正常的情况的差分,并评价所述差分,所述评价单元对所述细胞的运动传播状态的变化进行定量评价。
(13)根据(12)所述的图像处理设备,其中,所述评价单元对所述细胞图像内没有运动的静止区域的变化进行评价。
(14)根据(12)或(13)所述的图像处理设备,其中,所述评价单元对所述运动量大的区域的移动速度的变化进行评价。
(15)根据(12)到(14)中任一项所述的图像处理设备,其中,所述评价单元对所述运动量的重心的轨迹的变化进评价。
(16)根据(12)到(15)中任一项所述的图像处理设备,其中,所述评价单元对有关所述细胞图像的运动的速度和方向的直方图进行评价。
(17)根据(10)到(16)中任一项所述的图像处理设备,其中,所述评价单元对在几个小时到几天的观察期间所述细胞的运动传播状态进行评价。
(18)根据(10)到(17)中任一项所述的图像处理设备,其中,所述评价单元评价所述细胞的运动传播状态,以便无论是否发生QT延长,都检测心肌细胞发生心律失常,所述QT延长为在心肌细胞的Q波和T波之间的时间延长。
(19)一种图像处理设备的图像处理方法,所述方法包括以下步骤:
由运动量生成单元对于包含细胞的一个或多个图像针对每个子区域生成运动量;以及
由评价单元使用生成的所述运动量的一部分或全部评价所述细胞的运动传播状态。
(20)一种程序,其促使计算机用作:
运动量生成单元,被配置为对于包含细胞的一个或多个图像针对每个子区域生成运动量;以及
评价单元,被配置为使用由所述运动量生成单元生成的所述运动量的一部分或全部评价所述细胞的运动传播状态。
参考符号列表
100 药物评价设备
101 成像单元
102 待评价对象的图像数据生成/记录单元
103 评价指标数据生成单元
104 评价单元
121 运动检测单元
122 运动量绝对值计算单元
123 运动量绝对值储存单元
141 帧存储器
142 运动向量计算单元
201 运动量绝对值获取单元
202 映射单元
203 随时间变化分析单元
204 运动评价单元
205 显示单元
206 输出单元
Claims (11)
1.一种图像处理设备,包括:
运动量生成单元,被配置为对于包含细胞的一个或多个图像,针对每个子区域生成运动量;以及
评价单元,使用由所述运动量生成单元生成的运动量的一部分或全部评价所述细胞的运动传播状态。
2.根据权利要求1所述的图像处理设备,其中,所述评价单元将由所述运动量表示的所述细胞的运动传播状态与确定为所述细胞的运动为正常的情况的比较结果作为图像进行显示。
3.根据权利要求1所述的图像处理设备,其中,通过获得所述细胞图像的运动量与确定为所述细胞的运动正常的情况的差分,并评价所述差分,所述评价单元对所述细胞的运动传播状态的变化进行定量评价。
4.根据权利要求3所述的图像处理设备,其中,所述评价单元对所述细胞图像内没有运动的静止区域的变化进行评价。
5.根据权利要求3所述的图像处理设备,其中,所述评价单元对所述运动量大的区域的移动速度的变化进行评价。
6.根据权利要求3所述的图像处理设备,其中,所述评价单元对所述运动量的重心的轨迹的变化进评价。
7.根据权利要求3所述的图像处理设备,其中,所述评价单元对有关所述细胞图像的运动的速度和方向的直方图进行评价。
8.根据权利要求1所述的图像处理设备,其中,所述评价单元对在几个小时到几天的观察期间所述细胞的运动传播状态进行评价。
9.根据权利要求1所述的图像处理设备,其中,所述评价单元评价所述细胞的运动传播状态,以便无论是否发生QT延长,都检测心肌细胞发生心律失常,所述QT延长为在心肌细胞的Q波和T波之间的时间延长。
10.一种图像处理设备的图像处理方法,所述方法包括以下步骤:
由运动量生成单元对于包含细胞的一个或多个图像针对每个子区域生成运动量;以及
由评价单元使用生成的所述运动量的一部分或全部评价所述细胞的运动传播状态。
11.根据权利要求10所述的图像处理方法,其中,
评价步骤包括通过所述评价单元将由所述运动量表示的所述细胞的运动传播状态与确定为所述细胞的运动为正常的情况的比较结果作为图像进行显示。
Applications Claiming Priority (5)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2011043330 | 2011-02-28 | ||
JP2011-043330 | 2011-02-28 | ||
JP2011-162640 | 2011-07-25 | ||
JP2011162640A JP6078943B2 (ja) | 2011-02-28 | 2011-07-25 | 画像処理装置および方法、並びに、プログラム |
CN2012800098935A CN103380209A (zh) | 2011-02-28 | 2012-02-28 | 显示控制设备和方法、图像处理设备和方法以及程序 |
Related Parent Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2012800098935A Division CN103380209A (zh) | 2011-02-28 | 2012-02-28 | 显示控制设备和方法、图像处理设备和方法以及程序 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105847633A true CN105847633A (zh) | 2016-08-10 |
CN105847633B CN105847633B (zh) | 2019-05-28 |
Family
ID=46757981
Family Applications (2)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2012800098935A Pending CN103380209A (zh) | 2011-02-28 | 2012-02-28 | 显示控制设备和方法、图像处理设备和方法以及程序 |
CN201610158846.7A Active CN105847633B (zh) | 2011-02-28 | 2012-02-28 | 图像处理设备及其图像处理方法 |
Family Applications Before (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2012800098935A Pending CN103380209A (zh) | 2011-02-28 | 2012-02-28 | 显示控制设备和方法、图像处理设备和方法以及程序 |
Country Status (6)
Country | Link |
---|---|
US (2) | US10311581B2 (zh) |
EP (1) | EP2682455B1 (zh) |
JP (1) | JP6078943B2 (zh) |
KR (1) | KR101900977B1 (zh) |
CN (2) | CN103380209A (zh) |
WO (1) | WO2012118049A1 (zh) |
Families Citing this family (23)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP5293587B2 (ja) | 2009-12-16 | 2013-09-18 | ソニー株式会社 | 表示制御装置、表示制御方法、プログラム |
JP5772817B2 (ja) | 2010-03-29 | 2015-09-02 | ソニー株式会社 | 画像処理装置および方法、並びに、プログラム |
JP6078943B2 (ja) | 2011-02-28 | 2017-02-15 | ソニー株式会社 | 画像処理装置および方法、並びに、プログラム |
JP5957948B2 (ja) | 2012-02-24 | 2016-07-27 | ソニー株式会社 | 表示制御装置、表示制御方法および記録媒体 |
JP6102166B2 (ja) * | 2012-10-10 | 2017-03-29 | 株式会社ニコン | 心筋細胞の運動検出方法、心筋細胞の培養方法、薬剤評価方法、画像処理プログラム及び画像処理装置 |
CN104884632B (zh) | 2012-12-27 | 2018-01-02 | 索尼公司 | 细胞分析系统及细胞分析方法 |
JP6494903B2 (ja) | 2013-02-14 | 2019-04-03 | ソニー株式会社 | 分析システム、分析プログラム及び分析方法 |
JP6217968B2 (ja) * | 2013-07-19 | 2017-10-25 | ソニー株式会社 | 画像処理装置および方法、並びにプログラム |
US10157459B2 (en) | 2013-07-19 | 2018-12-18 | Sony Corporation | Cell evaluation apparatus and method to search for an area to be observed in cardiomyocytes |
JP6454478B2 (ja) | 2013-09-12 | 2019-01-16 | オリンパス株式会社 | 心筋細胞への分化をモニタリングする方法 |
JP6942925B2 (ja) * | 2013-11-08 | 2021-09-29 | ソニーグループ株式会社 | 細胞分析システム、細胞分析プログラム及び細胞分析方法 |
JP6459458B2 (ja) * | 2014-03-03 | 2019-01-30 | ソニー株式会社 | 細胞評価装置 |
JP6274515B2 (ja) * | 2014-03-05 | 2018-02-07 | キヤノンメディカルシステムズ株式会社 | 医用画像処理装置 |
JP2016041023A (ja) * | 2014-08-14 | 2016-03-31 | ソニー株式会社 | 画像処理装置、画像処理方法及び画像処理プログラム。 |
JP2016140317A (ja) * | 2015-02-03 | 2016-08-08 | ソニー株式会社 | 情報処理装置、情報処理方法、プログラム及び観察システム |
JP6661165B2 (ja) * | 2015-03-27 | 2020-03-11 | 大日本印刷株式会社 | 画像解析装置、画像解析方法、及び画像解析プログラム |
GB201507454D0 (en) * | 2015-04-30 | 2015-06-17 | Phase Focus Ltd | Method and apparatus for determining temporal behaviour of an object |
JPWO2016189583A1 (ja) * | 2015-05-22 | 2018-03-15 | オリンパス株式会社 | 生体試料の解析方法 |
US10983039B2 (en) | 2016-01-27 | 2021-04-20 | Hitachi High-Tech Corporation | Observation device |
EP3611696A4 (en) * | 2017-02-16 | 2021-02-17 | Kyoto University | CELL EVALUATION PROCESS, CELL EVALUATION DEVICE AND CELL EVALUATION PROGRAM |
JP6424911B2 (ja) * | 2017-02-28 | 2018-11-21 | 株式会社ニコン | 心筋細胞の運動検出方法、心筋細胞の培養方法、薬剤評価方法、画像処理プログラム及び画像処理装置 |
JP6891299B2 (ja) * | 2018-01-04 | 2021-06-18 | オリンパス株式会社 | 細胞画像処理装置、細胞画像処理方法および細胞画像処理プログラム |
JP2019154816A (ja) * | 2018-03-13 | 2019-09-19 | ソニー・オリンパスメディカルソリューションズ株式会社 | 医療用画像処理装置、医療用観察装置、及び医療用観察装置の作動方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1728856A1 (en) * | 2005-06-01 | 2006-12-06 | Fraunhofer-Gesellschaft zur Förderung der angewandten Forschung e.V. | A method for optically determining the dynamic behaviour of contracting muscle cells |
JP2007121106A (ja) * | 2005-10-27 | 2007-05-17 | Foundation For Biomedical Research & Innovation | 心筋細胞のモニター方法およびモニター装置 |
US20080304732A1 (en) * | 2007-06-05 | 2008-12-11 | General Electric Company | Automatic characterization of cellular motion |
CN101404941A (zh) * | 2006-03-21 | 2009-04-08 | 皇家飞利浦电子股份有限公司 | 用于彩色组织多普勒成像的速度标尺的优化 |
CN101903532A (zh) * | 2008-03-24 | 2010-12-01 | 株式会社尼康 | 细胞观察的图像解析方法、图像处理程序和图像处理装置 |
Family Cites Families (25)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPS63233392A (ja) | 1987-03-20 | 1988-09-29 | Nippon Petrochem Co Ltd | 周期的に形状を変える物体の変形周期測定法 |
EP1152694B1 (en) * | 1999-12-07 | 2004-07-21 | Koninklijke Philips Electronics N.V. | Ultrasonic image processing method and system for displaying a composite image sequence of an artery segment |
US6555365B2 (en) * | 2000-12-07 | 2003-04-29 | Biocrystal, Ltd. | Microincubator comprising a cell culture apparatus and a transparent heater |
AU2002235939A1 (en) * | 2001-03-08 | 2002-09-19 | Universite Joseph Fourier | Quantitative analysis, visualization and movement correction in dynamic processes |
AU2003211104B2 (en) * | 2002-02-13 | 2009-01-29 | Reify Corporation | Method and apparatus for acquisition, compression, and characterization of spatiotemporal signals |
JP4412894B2 (ja) | 2002-11-28 | 2010-02-10 | 富士フイルム株式会社 | 画像抽出装置 |
US20050203708A1 (en) * | 2004-03-11 | 2005-09-15 | Srinka Ghosh | Method and system for microarray gradient detection and characterization |
JP2007222073A (ja) | 2006-02-23 | 2007-09-06 | Yamaguchi Univ | 画像処理により細胞運動特性を評価する方法、そのための画像処理装置及び画像処理プログラム |
CN101730736B (zh) | 2007-06-08 | 2014-01-15 | 国立大学法人东京医科齿科大学 | 心脏折返模型芯片以及利用心脏折返模型芯片的药剂评价装置和方法 |
CN101711276B (zh) | 2007-06-08 | 2014-05-07 | 国立大学法人东京医科齿科大学 | 模型细胞芯片,利用模型细胞芯片的药效评价装置以及药效评价方法 |
CN101945615B (zh) * | 2008-03-07 | 2014-01-29 | 株式会社日立医疗器械 | 超声波摄像装置 |
JP2009229274A (ja) | 2008-03-24 | 2009-10-08 | Nikon Corp | 細胞観察の画像解析方法、画像処理プログラム及び画像処理装置 |
JP2009229275A (ja) * | 2008-03-24 | 2009-10-08 | Nikon Corp | 細胞観察画像の画像解析方法、画像処理プログラム及び画像処理装置 |
JP2010026392A (ja) | 2008-07-23 | 2010-02-04 | Nikon Corp | 細胞観察画像の画像解析方法、画像処理プログラム及び画像処理装置 |
JP5438936B2 (ja) * | 2008-08-29 | 2014-03-12 | 株式会社東芝 | 超音波診断装置、画像処理装置及び画像処理プログラム |
JP5228826B2 (ja) * | 2008-11-17 | 2013-07-03 | 大日本印刷株式会社 | 細胞画像解析装置および細胞画像解析方法 |
JP5614928B2 (ja) * | 2008-12-05 | 2014-10-29 | 国立大学法人 東京医科歯科大学 | 心筋毒性検査装置、心筋毒性検査チップおよび心筋毒性検査方法 |
JP2010169823A (ja) | 2009-01-21 | 2010-08-05 | Nikon Corp | 細胞観察画像の画像処理方法、画像処理プログラム及び画像処理装置 |
JP5582520B2 (ja) | 2009-02-18 | 2014-09-03 | 国立大学法人名古屋大学 | 血流観測装置 |
JP5902092B2 (ja) * | 2009-10-19 | 2016-04-13 | セルラー ダイナミクス インターナショナル, インコーポレイテッド | 心筋細胞の生成 |
JP5293587B2 (ja) | 2009-12-16 | 2013-09-18 | ソニー株式会社 | 表示制御装置、表示制御方法、プログラム |
JP5772817B2 (ja) * | 2010-03-29 | 2015-09-02 | ソニー株式会社 | 画像処理装置および方法、並びに、プログラム |
JP6078943B2 (ja) | 2011-02-28 | 2017-02-15 | ソニー株式会社 | 画像処理装置および方法、並びに、プログラム |
JP5957948B2 (ja) | 2012-02-24 | 2016-07-27 | ソニー株式会社 | 表示制御装置、表示制御方法および記録媒体 |
WO2015041177A1 (ja) * | 2013-09-18 | 2015-03-26 | 株式会社ニコン | 画像解析装置、画像解析方法、画像解析プログラム、細胞の製造方法、細胞の培養方法、および細胞の製造装置 |
-
2011
- 2011-07-25 JP JP2011162640A patent/JP6078943B2/ja active Active
-
2012
- 2012-02-28 WO PCT/JP2012/054857 patent/WO2012118049A1/ja active Application Filing
- 2012-02-28 US US13/984,672 patent/US10311581B2/en active Active
- 2012-02-28 KR KR1020137021429A patent/KR101900977B1/ko active IP Right Grant
- 2012-02-28 EP EP12752715.8A patent/EP2682455B1/en active Active
- 2012-02-28 CN CN2012800098935A patent/CN103380209A/zh active Pending
- 2012-02-28 CN CN201610158846.7A patent/CN105847633B/zh active Active
-
2019
- 2019-04-22 US US16/391,157 patent/US10650534B2/en active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1728856A1 (en) * | 2005-06-01 | 2006-12-06 | Fraunhofer-Gesellschaft zur Förderung der angewandten Forschung e.V. | A method for optically determining the dynamic behaviour of contracting muscle cells |
JP2007121106A (ja) * | 2005-10-27 | 2007-05-17 | Foundation For Biomedical Research & Innovation | 心筋細胞のモニター方法およびモニター装置 |
CN101404941A (zh) * | 2006-03-21 | 2009-04-08 | 皇家飞利浦电子股份有限公司 | 用于彩色组织多普勒成像的速度标尺的优化 |
US20080304732A1 (en) * | 2007-06-05 | 2008-12-11 | General Electric Company | Automatic characterization of cellular motion |
CN101903532A (zh) * | 2008-03-24 | 2010-12-01 | 株式会社尼康 | 细胞观察的图像解析方法、图像处理程序和图像处理装置 |
Non-Patent Citations (3)
Title |
---|
A.KAMGOUE,J. OHAYON 等: "Quantification of Cardiomyocyte Contraction Based on Image Correlation Analysis", 《CYTOMERY PART A》 * |
KANG LI 等: "Cell Population Tracking and Lineage Construction with Spatiotemporal Context", 《MEDICAL IMAGE ANALYSIS》 * |
STACEY A.MASKARINEC 等: "Quantifying Cellular Traction Forces in Three Dimensions", 《PRO NALTL ACAD SCI USA》 * |
Also Published As
Publication number | Publication date |
---|---|
CN105847633B (zh) | 2019-05-28 |
US20190311483A1 (en) | 2019-10-10 |
JP2012194168A (ja) | 2012-10-11 |
JP6078943B2 (ja) | 2017-02-15 |
US10650534B2 (en) | 2020-05-12 |
KR101900977B1 (ko) | 2018-09-20 |
WO2012118049A1 (ja) | 2012-09-07 |
EP2682455A4 (en) | 2017-04-12 |
US20130321459A1 (en) | 2013-12-05 |
US10311581B2 (en) | 2019-06-04 |
CN103380209A (zh) | 2013-10-30 |
KR20140017552A (ko) | 2014-02-11 |
EP2682455A1 (en) | 2014-01-08 |
EP2682455B1 (en) | 2020-11-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105847633A (zh) | 图像处理设备及其图像处理方法 | |
CN102568001B (zh) | 图像处理设备和方法 | |
CN105378052B (zh) | 细胞评估装置和方法以及程序 | |
Baker et al. | Which way do I go? Neural activation in response to feedback and spatial processing in a virtual T-maze | |
JP5860868B2 (ja) | 生体の運動を決定するシステム及び方法 | |
CN105310703B (zh) | 一种基于用户瞳孔直径数据获取主观满意度的方法 | |
CN109938727A (zh) | 非人灵长类动物三维视觉刺激实验系统和方法 | |
US9972097B2 (en) | Image processing device, method, and medium for calculating propagation speed and direction of object pulsations | |
Göttlich et al. | Structural covariance of amygdala subregions is associated with trait aggression and endogenous testosterone in healthy individuals | |
EP0973431B1 (en) | Enhancing and extracting second order nonlinear response components | |
Wimmer et al. | Episodic memory retrieval is supported by rapid replay of episode content | |
Sebille et al. | Calcium sparks in muscle cells: interactive procedures for automatic detection and measurements on line-scan confocal images series | |
CN208709883U (zh) | 非人灵长类动物用的三维视觉刺激实验系统 | |
JP6436142B2 (ja) | 画像処理装置および方法、並びに、プログラム | |
Xue et al. | Inferring monosynaptic connections from paired dendritic spine Ca2+ imaging and large-scale recording of extracellular spiking | |
Engelmann et al. | Brain measurement and manipulation methods | |
Martinovic et al. | Division and spreading of attention across color | |
Xie | Development of an Artificial Intelligence-Based Video Monitoring System to Investigate Drug Effects on Daphnia Magna | |
Saberigarakani et al. | Light-field microscopy to study cardiac arrhythmias in zebrafish larvae | |
Anderson et al. | Chapter Five: Quantifying Two-Dimensional and Three-Dimensional Learning in Anatomy Using Electroencephalography | |
PAGE | 5. From Correlates to Theories of Consciousness | |
RU2106802C1 (ru) | Способ диагностики биообъектов | |
Holroyd et al. | Reward Processing when Evaluating Goals: Insight into Hierarchical Reinforcement Learning Chad C. Williams University of Victoria | |
JP2013053854A (ja) | 画像処理装置および方法、記録媒体、並びに、プログラム | |
Kaplan | Brain oscillations and novelty processing in human spatial memory |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |