TWI835323B - 用於多工螢光原位雜合影像之獲取及處理的系統及方法 - Google Patents

用於多工螢光原位雜合影像之獲取及處理的系統及方法 Download PDF

Info

Publication number
TWI835323B
TWI835323B TW111137276A TW111137276A TWI835323B TW I835323 B TWI835323 B TW I835323B TW 111137276 A TW111137276 A TW 111137276A TW 111137276 A TW111137276 A TW 111137276A TW I835323 B TWI835323 B TW I835323B
Authority
TW
Taiwan
Prior art keywords
image
images
sample
different
pixel
Prior art date
Application number
TW111137276A
Other languages
English (en)
Other versions
TW202307420A (zh
Inventor
張蘊青
謝丹
克洛伊 金
Original Assignee
美商應用材料股份有限公司
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by 美商應用材料股份有限公司 filed Critical 美商應用材料股份有限公司
Publication of TW202307420A publication Critical patent/TW202307420A/zh
Application granted granted Critical
Publication of TWI835323B publication Critical patent/TWI835323B/zh

Links

Classifications

    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16BBIOINFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR GENETIC OR PROTEIN-RELATED DATA PROCESSING IN COMPUTATIONAL MOLECULAR BIOLOGY
    • G16B40/00ICT specially adapted for biostatistics; ICT specially adapted for bioinformatics-related machine learning or data mining, e.g. knowledge discovery or pattern finding
    • G16B40/10Signal processing, e.g. from mass spectrometry [MS] or from PCR
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/62Systems in which the material investigated is excited whereby it emits light or causes a change in wavelength of the incident light
    • G01N21/63Systems in which the material investigated is excited whereby it emits light or causes a change in wavelength of the incident light optically excited
    • G01N21/64Fluorescence; Phosphorescence
    • G01N21/6428Measuring fluorescence of fluorescent products of reactions or of fluorochrome labelled reactive substances, e.g. measuring quenching effects, using measuring "optrodes"
    • CCHEMISTRY; METALLURGY
    • C12BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
    • C12QMEASURING 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/00Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions
    • C12Q1/68Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions involving nucleic acids
    • C12Q1/6813Hybridisation assays
    • C12Q1/6841In situ hybridisation
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/62Systems in which the material investigated is excited whereby it emits light or causes a change in wavelength of the incident light
    • G01N21/63Systems in which the material investigated is excited whereby it emits light or causes a change in wavelength of the incident light optically excited
    • G01N21/64Fluorescence; Phosphorescence
    • G01N21/645Specially adapted constructive features of fluorimeters
    • G01N21/6456Spatial resolved fluorescence measurements; Imaging
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/62Systems in which the material investigated is excited whereby it emits light or causes a change in wavelength of the incident light
    • G01N21/63Systems in which the material investigated is excited whereby it emits light or causes a change in wavelength of the incident light optically excited
    • G01N21/64Fluorescence; Phosphorescence
    • G01N21/645Specially adapted constructive features of fluorimeters
    • G01N21/6456Spatial resolved fluorescence measurements; Imaging
    • G01N21/6458Fluorescence microscopy
    • GPHYSICS
    • G02OPTICS
    • G02BOPTICAL ELEMENTS, SYSTEMS OR APPARATUS
    • G02B21/00Microscopes
    • G02B21/0004Microscopes specially adapted for specific applications
    • G02B21/002Scanning microscopes
    • G02B21/0024Confocal scanning microscopes (CSOMs) or confocal "macroscopes"; Accessories which are not restricted to use with CSOMs, e.g. sample holders
    • G02B21/0052Optical details of the image generation
    • G02B21/0076Optical details of the image generation arrangements using fluorescence or luminescence
    • GPHYSICS
    • G02OPTICS
    • G02BOPTICAL ELEMENTS, SYSTEMS OR APPARATUS
    • G02B21/00Microscopes
    • G02B21/16Microscopes adapted for ultraviolet illumination ; Fluorescence microscopes
    • GPHYSICS
    • G02OPTICS
    • G02BOPTICAL ELEMENTS, SYSTEMS OR APPARATUS
    • G02B21/00Microscopes
    • G02B21/36Microscopes arranged for photographic purposes or projection purposes or digital imaging or video purposes including associated control and data processing arrangements
    • G02B21/365Control or image processing arrangements for digital or video microscopes
    • GPHYSICS
    • G02OPTICS
    • G02BOPTICAL ELEMENTS, SYSTEMS OR APPARATUS
    • G02B21/00Microscopes
    • G02B21/36Microscopes arranged for photographic purposes or projection purposes or digital imaging or video purposes including associated control and data processing arrangements
    • G02B21/365Control or image processing arrangements for digital or video microscopes
    • G02B21/367Control or image processing arrangements for digital or video microscopes providing an output produced by processing a plurality of individual source images, e.g. image tiling, montage, composite images, depth sectioning, image comparison
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T3/00Geometric image transformations in the plane of the image
    • G06T3/40Scaling of whole images or parts thereof, e.g. expanding or contracting
    • G06T3/4038Image mosaicing, e.g. composing plane images from plane sub-images
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T5/00Image enhancement or restoration
    • G06T5/20Image enhancement or restoration using local operators
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T5/00Image enhancement or restoration
    • G06T5/50Image enhancement or restoration using two or more images, e.g. averaging or subtraction
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T5/00Image enhancement or restoration
    • G06T5/73Deblurring; Sharpening
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T5/00Image enhancement or restoration
    • G06T5/80Geometric correction
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/0002Inspection of images, e.g. flaw detection
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/0002Inspection of images, e.g. flaw detection
    • G06T7/0012Biomedical image inspection
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/30Determination of transform parameters for the alignment of images, i.e. image registration
    • G06T7/33Determination of transform parameters for the alignment of images, i.e. image registration using feature-based methods
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06VIMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V10/00Arrangements for image or video recognition or understanding
    • G06V10/20Image preprocessing
    • G06V10/30Noise filtering
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06VIMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V10/00Arrangements for image or video recognition or understanding
    • G06V10/40Extraction of image or video features
    • G06V10/50Extraction of image or video features by performing operations within image blocks; by using histograms, e.g. histogram of oriented gradients [HoG]; by summing image-intensity values; Projection analysis
    • G06V10/507Summing image-intensity values; Histogram projection analysis
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06VIMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V10/00Arrangements for image or video recognition or understanding
    • G06V10/70Arrangements for image or video recognition or understanding using pattern recognition or machine learning
    • G06V10/74Image or video pattern matching; Proximity measures in feature spaces
    • G06V10/75Organisation of the matching processes, e.g. simultaneous or sequential comparisons of image or video features; Coarse-fine approaches, e.g. multi-scale approaches; using context analysis; Selection of dictionaries
    • G06V10/758Involving statistics of pixels or of feature values, e.g. histogram matching
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06VIMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V20/00Scenes; Scene-specific elements
    • G06V20/60Type of objects
    • G06V20/69Microscopic objects, e.g. biological cells or cellular parts
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06VIMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
    • G06V20/00Scenes; Scene-specific elements
    • G06V20/60Type of objects
    • G06V20/69Microscopic objects, e.g. biological cells or cellular parts
    • G06V20/698Matching; Classification
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16BBIOINFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR GENETIC OR PROTEIN-RELATED DATA PROCESSING IN COMPUTATIONAL MOLECULAR BIOLOGY
    • G16B50/00ICT programming tools or database systems specially adapted for bioinformatics
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16BBIOINFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR GENETIC OR PROTEIN-RELATED DATA PROCESSING IN COMPUTATIONAL MOLECULAR BIOLOGY
    • G16B50/00ICT programming tools or database systems specially adapted for bioinformatics
    • G16B50/40Encryption of genetic data
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/62Systems in which the material investigated is excited whereby it emits light or causes a change in wavelength of the incident light
    • G01N21/63Systems in which the material investigated is excited whereby it emits light or causes a change in wavelength of the incident light optically excited
    • G01N21/64Fluorescence; Phosphorescence
    • G01N2021/6417Spectrofluorimetric devices
    • G01N2021/6419Excitation at two or more wavelengths
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/62Systems in which the material investigated is excited whereby it emits light or causes a change in wavelength of the incident light
    • G01N21/63Systems in which the material investigated is excited whereby it emits light or causes a change in wavelength of the incident light optically excited
    • G01N21/64Fluorescence; Phosphorescence
    • G01N2021/6417Spectrofluorimetric devices
    • G01N2021/6421Measuring at two or more wavelengths
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/62Systems in which the material investigated is excited whereby it emits light or causes a change in wavelength of the incident light
    • G01N21/63Systems in which the material investigated is excited whereby it emits light or causes a change in wavelength of the incident light optically excited
    • G01N21/64Fluorescence; Phosphorescence
    • G01N21/6428Measuring fluorescence of fluorescent products of reactions or of fluorochrome labelled reactive substances, e.g. measuring quenching effects, using measuring "optrodes"
    • G01N2021/6439Measuring fluorescence of fluorescent products of reactions or of fluorochrome labelled reactive substances, e.g. measuring quenching effects, using measuring "optrodes" with indicators, stains, dyes, tags, labels, marks
    • G01N2021/6441Measuring fluorescence of fluorescent products of reactions or of fluorochrome labelled reactive substances, e.g. measuring quenching effects, using measuring "optrodes" with indicators, stains, dyes, tags, labels, marks with two or more labels
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/75Systems in which material is subjected to a chemical reaction, the progress or the result of the reaction being investigated
    • G01N2021/755Comparing readings with/without reagents, or before/after reaction
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/01Arrangements or apparatus for facilitating the optical investigation
    • G01N21/03Cuvette constructions
    • G01N21/05Flow-through cuvettes
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2201/00Features of devices classified in G01N21/00
    • G01N2201/12Circuits of general importance; Signal processing
    • G01N2201/127Calibration; base line adjustment; drift compensation
    • G01N2201/12746Calibration values determination
    • G01N2201/12761Precalibration, e.g. for a given series of reagents
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/10Image acquisition modality
    • G06T2207/10056Microscopic image
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/10Image acquisition modality
    • G06T2207/10064Fluorescence image
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/30Subject of image; Context of image processing
    • G06T2207/30004Biomedical image processing
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/30Subject of image; Context of image processing
    • G06T2207/30168Image quality inspection
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2207/00Indexing scheme for image analysis or image enhancement
    • G06T2207/30Subject of image; Context of image processing
    • G06T2207/30204Marker

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Chemical & Material Sciences (AREA)
  • General Health & Medical Sciences (AREA)
  • Multimedia (AREA)
  • Analytical Chemistry (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Immunology (AREA)
  • Optics & Photonics (AREA)
  • Medical Informatics (AREA)
  • Biochemistry (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Molecular Biology (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • Pathology (AREA)
  • Biophysics (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Biotechnology (AREA)
  • Bioethics (AREA)
  • Databases & Information Systems (AREA)
  • Organic Chemistry (AREA)
  • Evolutionary Biology (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Artificial Intelligence (AREA)
  • Evolutionary Computation (AREA)
  • Software Systems (AREA)
  • Biomedical Technology (AREA)
  • Wood Science & Technology (AREA)
  • Genetics & Genomics (AREA)
  • Zoology (AREA)
  • Proteomics, Peptides & Aminoacids (AREA)
  • Quality & Reliability (AREA)
  • Chemical Kinetics & Catalysis (AREA)
  • Data Mining & Analysis (AREA)
  • Public Health (AREA)
  • Signal Processing (AREA)

Abstract

在一種螢光原位雜合成像系統中,作為巢狀迴圈地進行以下步驟:(1)一閥依序將一流通槽耦接至複數個不同試劑來源,以將樣本暴露至複數個不同試劑,(2)對於該複數個不同試劑的各試劑,一馬達依序將螢光顯微鏡相對於位於複數個不同視野處之樣本定位,(3)對於該複數個不同視野中的各視野,一可變頻激發光源依序發射複數個不同波長,(4)對於該複數個不同波長中的各波長,一致動器依序將該螢光顯微鏡相對於位於複數個不同垂直高度處之樣本定位,及(5)對於該複數個不同垂直高度中的各垂直高度,獲得一影像。

Description

用於多工螢光原位雜合影像之獲取及處理的系統及方法
本說明書有關多工螢光原位雜合法成像。
對於生物技術社群及製藥產業而言開發用於視覺化及量化一生物樣本(例如組織切除部、活組織切片、培養細胞)之內的多個生物分析物(例如DNA、RNA、及蛋白質)的方法有很大的興趣。科學家利用此類方法診斷/監測疾病、驗證生物標記、及調查治療方式。至今,例示性方法包括對一生物樣本藉功能性領域所標記之抗體或寡核苷酸(例如RNA或DNA)的多重成像。
多工螢光原位雜合法(mFISH)成像是空間轉錄體學中決定基因表現的強大技術。簡言之,樣本被曝露至多個寡核苷酸探針,該些探針目標為感興趣的RNA。這些探針具有不同標記方案,當補充的、經螢光標記的探針被導入該樣本時,該些不同標記方案將允許區分不同RNA物 種。接著藉由曝光至不同波長的激發光,獲取了連續數回合的螢光影像。對於各給定像素,來自不同波長之激發光的不同影像之螢光強度形成一信號序列。此序列接著與來自一碼簿的參考碼庫做比較,該碼簿將各代碼與一基因相關聯。最匹配的參考碼被用以識別出被表現於該影像中該像素處的一相關聯基因。
在一態樣中,一種螢光原位雜合成像系統包括:用以容納將暴露於一試劑中之螢光探針的一樣本的一流通槽、用以控制來自該流通槽之複數個試劑來源中之一者的流動的一閥、用以致使通過該流通槽的流體流動的一泵、包括一可變頻激發光源和經定位以接收從該樣本螢光發射之光的一攝影機的一螢光顯微鏡、用以致使該流通槽及該螢光顯微鏡之間的相對垂直運動的一致動器、用以致使該流通槽及該螢光顯微鏡之間的相對橫向運動的一馬達、及一控制系統。該控制系統經配置以作為巢狀迴圈地進行下列步驟:致使該閥依序將該流通槽耦接至複數個不同試劑來源,以將該樣本暴露至複數個不同試劑,對於該複數個不同試劑中的各試劑,致使該馬達依序將該螢光顯微鏡相對於位於複數個不同視野處之樣本定位,對於該複數個不同視野中的各視野,致使該可變頻激發光源依序發射複數個不同波長,對於該複數個不同波長中的各波長,致使該致動器依序將該螢光顯微鏡相對於位於複數個不同垂直高 度處之樣本定位,對於該複數個不同垂直高度中的各垂直高度,獲得位於該相應垂直高度處的一影像,該相應垂直高度涵蓋該樣本的該相應視野,該樣本讓該相應試劑的相應螢光探針被激發於該相應波長。
可在實施中包括以下特徵之一或更多者。該致動器可為壓電致動器。一資料處理系統可經配置以決定對於個影像的一影像品質值,及若影像品質值不符合一品質閾值則篩除(screen)影像。該影像品質值可為清晰度品質值或亮度品質值。該攝影機可為CMOS攝影機。
另一態樣中,一種螢光原位雜合成像及分析系統包括:用以容納將暴露於一試劑中之螢光探針的一樣本的一流通槽、用以依序於複數個不同成像參數組合獲得收集該樣本的複數個影像的一螢光顯微鏡、及一資料處理系統。該資料處理系統包括一線上預處理系統及一離線處理系統,該線上預處理系統經配置以隨著該等影像被收集而從該螢光顯微鏡依序接收該等影像,並進行即時(on-the-fly)影像預處理以移除該等影像的實驗性假影及提供RNA影像斑點銳化,該離線處理系統經配置以在收集該複數個影像之後,對具有一相同視野之影像進行對位,並將該複數個影像中的強度值解碼以辨識出所表現的基因。
可在實施中包括以下特徵之一或更多者。該資料處理系統可經配置以決定一清晰度品質值或一亮度品質值。該清晰度品質值的決定可藉由整合整個該影像之影像 亮度的梯度的絕對值。若該清晰度品質值在一最小清晰度值之上,則可決定該影像品質值符合閾值。該亮度品質值的決定可藉由整合整個影像的影像亮度。若該亮度品質值在一最小亮度值之上,則可決定該影像品質值符合閾值。
另一態樣中,一種螢光原位雜合成像系統包括:用以容納將暴露至一試劑中之螢光探針的一樣本的一流通槽、一螢光顯微鏡、及一控制系統。該螢光顯微鏡包括照亮該樣本的一可變頻激發光源、在一濾光輪上的複數個發射帶通濾光片、用以旋轉該濾光輪的一致動器、及經定位以接收從該樣本螢光發射且已通過該濾光輪上的一濾光片的光的一攝影機。該控制系統經配置以致使該可變頻激發光源發射具有一所選波長的一光束,致使該致動器旋轉該濾光輪以定位將一所選濾光片定位在該樣本與該攝影機之間的一光路徑中,從該攝影機獲得一影像,及協調該可變頻激發光源及濾光輪,使得當獲得該影像時,該所選濾光片具有關聯於該等螢光探針被該所選波長激發時的發射的一發射帶通。
另一態樣中,一種電腦程式產品具有指令以致使一或更多個處理器對複數個影像的各影像應用一或更多個濾光片並在該複數個影像的各影像上進行解卷積,以產生複數個經過濾的影像,對該複數個經過濾的影像進行對位以產生複數個經對位及過濾的影像,對該複數個經對位及過濾的影像之至少兩個像素的各像素藉由以下步驟來解碼該像素:從一碼簿中的複數個碼字識別出一碼字,該碼字 提供對於該像素而言該複數個經對位及過濾的影像之資料數值的一最佳匹配,以及,對至少兩個像素之各像素,決定關聯於該碼字的一基因,並儲存有關該像素處表現了該基因的一指示。
可在實施中包括以下特徵之一或更多者。可有指令用以決定該閾值,藉由決定用於該影像之亮度值直方圖,並決定在該直方圖中一分位數極限處的亮度值。該分位數極限可為99.5%或更多,例如99.9%以上。用以應用一或更多個濾光片的指令可包括:用以應用一低通濾光片的指令。用以進行解卷積的指令可包括用以藉一2D點擴散函數來解卷積的指令。
另一態樣中,一種電腦程式產品具有指令以接收一樣本的複數個初始影像,該複數個初始影像共享一共同視野。對於該複數個初始影像的各初始影像,該初始影像經分割成為具相同尺寸的複數個次影像,從該複數個次影像計算出一合併影像,該合併影像的尺寸等於一個次影像的尺寸,藉而產生複數個縮小的影像,使該複數個影像的各影像具有一對應縮小影像,對於該複數個縮小影像的至少一個縮小影像,決定一轉換來對該複數個縮小影像進行對位,對於該至少一個縮小的影像的各者,對該對應初始影像應用針對該縮小的影像所決定的該轉換,來產生一經轉換影像,該經轉換影像的尺寸等於該初始影像的尺寸,並藉而產生複數個經轉換影像,及對於該複數個經轉換影像之至少兩個像素的各像素,藉由以下步驟來解碼該像 素:從一碼簿中的複數個碼字識別出一碼字,該碼字提供對於該像素而言該複數個經轉換影像中之資料數值的一最佳匹配。對於至少兩個像素之各像素,決定關聯於該碼字的一基因,並儲存有關該像素處表現了該基因的一指示。
可在實施中包括以下特徵之一或更多者。用以將該初始影像分割成複數個次影像的指令可包括指令以將該初始影像分割成一矩形網格,其具有複數個行(column)及複數個列(row)。該矩形網格可有相等數目的行與列。初始影像可被分割成剛好N個影像,其中N=Z2而Z是一偶數整數。N可為2。次影像可具有不超過1024x1024像素的尺寸。
另一態樣中,一種電腦程式產品具有指令以接收代表一樣本之一共同視野的複數個影像,對該複數個影像進行對位以產生複數個經對位影像,對於該經對位影像之複數個像素的各像素藉以下步驟來解碼該像素:從一碼簿中的複數個碼字識別出一碼字,該碼字提供對於該像素而言該複數個經對位影像中之資料數值的一最佳匹配,以及對於經識別為對於一或更多個像素之一最佳匹配的各碼字以及對於該一或更多個像素之各像素,決定是否針對該像素的一影像字的一位元比符合用於該碼字的一閾值。該影像字是形成自針對該像素的該複數個經對位影像中的該等資料數值。對於經決定為符合該閾值的至少一個像素,決定關聯於該碼字的一基因並儲存有關該基因被表現於該像素處的一指示。其位元比不符合該閾值的像素被篩除。
可在實施中包括以下特徵之一或更多者。決定該閾值的步驟可包括決定在來自該樣本之該複數個影像的位元值直方圖中位於一第一分位數極限處的一值。該分位數極限可具有一預設值。決定該閾值之後,可藉由決定在來自該樣本之該複數個影像的位元值直方圖中位於不同的第二分位數極限處的一值來決定一更新閾值。
另一態樣中,一種電腦程式產品具有用以進行以下步驟的指令:接收一樣本的複數個影像,該複數個影像代表一共同視野,進行該複數個影像的對位以產生複數個經對位影像,對於該經對位影像之複數個像素的各像素來解碼該像素,其乃藉由從一碼簿中的複數個碼字識別出一碼字,該碼字提供對於該像素而言該複數個經對位影像中之資料數值的一最佳匹配,及對於經識別為對於一或更多個像素之一最佳匹配的各碼字以及對於該一或更多個像素之各像素,決定是否針對該像素的一影像字的一位元亮度符合用於該碼字的一閾值。該影像字是形成自針對該像素的該複數個經對位影像中的該等資料數值。對於經決定為符合該閾值的各像素,決定關聯於該碼字的一基因並儲存有關該像素處表現了該基因的一指示。位元比不符合該閾值的像素被篩除。
可在實施中包括以下特徵之一或更多者。決定該閾值的步驟可包括決定在來自該樣本之該複數個影像的位元值直方圖中位於一第一分位數極限處的一值。決定該閾值之後,可藉由決定在來自該樣本之該複數個影像的位元 值直方圖中位於不同的第二分位數極限處的一值來決定一更新閾值。
另一態樣中,一種電腦程式產品具有指令以接收一樣本的複數個初始影像。該複數個初始影像包括複數個不同視野之各視野的多重影像。指令經儲存以進行一影像處理演算法,在該演算法中包括對該視野之多重影像進行影像過濾或對該視野之多重影像進行對位中之一或更多者。對於各視野,決定一或更多個初始參數值,其中不同視野具有不同初始參數值。使用用於不同視野之不同初始解碼參數值,來對各視野進行該影像處理演算法。決定一或更多個接續的解碼參數值,並利用用於該不同視野的相同接續參數值,對各視野進行該影像處理演算法及一解碼演算法,來產生一組空間轉錄體學資料(spatial transcriptomic data)。
可在實施中包括以下特徵之一或更多者。對於一視野之複數個像素的各像素,該影像處理演算法及解碼演算法可藉由從一碼簿之複數個碼字識別出提供對複數個影像中之該像素的資料數值的最佳匹配的一碼字來解碼該像素,以產生複數個最匹配的碼字。對於各視野,可決定在各像素字與對該像素字而言之最佳匹配碼字之間的適配度,藉而提供對於各視野的複數個適配度。決定該適配度的操作可包括決定該像素字與該最佳匹配碼字之間的一距離。對於各視野,可從用於該視野之該複數個適配度決定一平均適配度,藉而提供複數個平均適配度。
另一態樣中,一種電腦程式產品具有指令以接收一樣本的複數個初始影像,該複數個初始影像包括用於複數個不同視野之各視野的多重影像。對於各視野,基於用於該視野的該多重影像產生一組空間轉錄體學資料。各組空間轉錄體學資料包括一或更多個資料組對,該一或更多個資料組對的各資料組對包括一位置與針對該位置的一經表現基因。決定來自鄰接視野之多重影像的重疊區域。從對於各視野的該組空間轉錄體學資料,產生用於涵蓋該複數個不同視野之一合併影像的一組合併空間轉錄體學資料。在產生該組合併空間轉錄體學資料之中,來自不同視野的至少一個經重複計數之資料組對被篩除,其代表在該樣本內之相同位置處的相同經表現基因。
可在實施中包括以下特徵之一或更多者。決定重疊區域的步驟可包括,對於各列,在來自第一FOV之一第一影像與來自一第二鄰接FOV之一第二影像之間進行水平位移估計。決定重疊區域的步驟可包括,對於各列,在來自第一FOV之該第一影像與來自該第二鄰接FOV之該第二影像之間進行垂直位移估計。該水平位移估計及該垂直位移估計可各包括該第一影像與該第二影像之間的交互相關性。
另一態樣中,一種電腦程式產品具有指令以接收一樣本的複數個影像,該複數個影像包括針對複數個不同視野之各視野的多重影像。在該複數個影像的各影像上,在該影像上進行影像品質檢核。進行該影像品質檢核之 後,從該複數個影像之各影像移除假影(artifact)。在移除假影之後,銳化在該複數個影像的各影像中的影像斑點。在銳化該等影像斑點之後,對於各多重影像,進行該多重影像的對位以產生複數個經對位影像。在進行對位之後,對於各多重影像,並對於一相應多重影像之該經對位影像的複數個像素的各像素,藉由以下步驟來解碼該像素:從一碼簿中的複數個碼字識別出一碼字,該碼字提供對於該像素而言該相應多重經對位影像中之資料數值的一最佳匹配。
本案之實施的優點能包括(但不限於)以下之一或更多者。
多工螢光原位雜合法(mFISH)影像的獲取與為決定基因表現的接續處理能以適合於工業採用的方式經自動化。能以統一的影像品質來處理需要多個視野(FOV)的大量樣本大小。能改善多樣本中或多模式研究中的行程對行程(run-to-run)一致性。要將不同時間所捕捉的相同FOV之螢光影像對位,能不使用分離的基準(例如不需要基準微粒(bead))、並使用比一些習用技術更少的資源且更準確地進行。
本案之一或更多個實施例的細節經闡述在隨附圖式及以下說明中。其他特徵、態樣、及優點將從說明、圖式及申請專利範圍中顯而易見。
100:mFISH成像及影像處理設備
10:樣本
110:流通槽
112:雜合試劑來源
112a:試劑
112b:清洗流體
114:閥
116:泵
118:載物台
118a:直線致動器
118b:Z軸致動器
119:化學廢料管理次系統
120:螢光顯微鏡
122:激發光源
122a:雷射模組
124:顯微鏡主體
128:濾光輪
128a:濾光片
128b:致動器
130:激發光
132:螢光
134:攝影機
136:物鏡
138:多頻帶分光鏡
140:控制系統
142:電腦
144:第一微控制器
146:第二微控制器
150:資料處理系統
152:資料處理裝置
154:儲存裝置
156:遠端存儲
202:步驟
202a:步驟
202b:步驟
203:步驟
204:步驟
206:步驟
208:步驟
210:步驟
212:步驟
214:步驟
252:步驟
254:步驟
302:影像
304a,304b,304c,304d:影像區塊
306:影像
308:堆疊
310:影像
312:堆疊
400:對數直方圖
402:第一峰值
404:第二峰值
410:閾值
602:直方圖
604:尾端
610:折線
620:直方圖
702:像素字之矩陣
704,724:列
706,726:行
710:像素字
722:碼簿/參考碼字之矩陣
730:參考碼字
802,812:直方圖
804,814:閾值
第1圖是用於多工螢光原位雜合法成像之設備的示意圖。
第2圖是資料處理方法的流程圖。
第3圖示意地圖示利用交互資訊演算法的影像對位程序。
第4圖是一影像中之強度的對數直方圖。
第5A圖是一所收集影像。
第5B圖圖示從第5A圖之影像所生成的一遮罩。
第6A圖是來自一經過濾且解卷積的影像的強度對數直方圖。
第6B圖是標準化之後的強度對數直方圖。
第7圖圖示一解碼方法。
第8A圖是對於一樣本碼字所識別的像素的位元比之對數直方圖。
第8B圖是對於一樣本碼字所識別的像素的亮度之對數直方圖。
第9圖是拼接(stitching)方法的流程圖。
第10圖是資料處理方法之另一實施的流程圖,在該資料處理方法中影像處理經即時進行。
不同圖式中的相似參考元件符號及標記指示出相似的元件。
雖然多工螢光原位雜合法(mFISH)成像是種強大的技術,目前為止其主要被使用在學術情境中。產業情境中將此種技術用於空間轉錄體學的運用(例如被進行研究中之生技公司利用針對空間轉錄體學之mFISH來開發意圖產品化的藥品),造成各式各樣的議題。例如:
●系統應被完全或接近完全地自動化,使得熟悉生物學但不熟悉影像處理或軟體程式化技術的使用者(例如研究科學家)僅需鍵入有關所調查之生物態樣的參數,而不需要做出有關影像處理的決策或具有程式專業。此種自動化事實上並非小事,而包含各種陷阱。
●當前的最新技術採取來自通用光學顯微鏡的樣本影像,其具有適合小樣本的小影像尺寸。在產業情境中運用的樣本可能更大且可需要多個視野以完整地成像。此種多個FOV能被人工地拼接,但半自動的做法遭受可能的問題像是不良的影像對位、過度計數經表現基因、或在處理拼接的影像的不同區域上的不一致性所致的缺少可信賴度。
●此程序需要足夠高的產出量以能在工業上可施行。然而,當前最新技術花費可觀的成像時間並在資料分析上需要許多人工步驟。對於大規模生命科學及臨床研究需要成像產出量及品質管控。此外,資料處理的可規模化能力(scalability)受限於運算時間。
●因為各FOV經個別處理,能發現跨於一感興趣區域的顯著變異。額外地,來自不同儀器的影像在被發送至軟體管線(pipeline)之前需要經良好地校準及轉換。習知地,影像處理參數需要被人工調整以容許系統或樣本變異。
●當前的解決方案強烈假設完美的影像,此在實驗條件與校準方法之間創造了差距。
以下討論的一整體架構及各種個別組件及步驟能解決上述議題之一或更多個。
一般,此整體架構及個別組件及步驟致能一個涉及最少人工步驟的端到端mFISH成像及解碼程序。該系統整合完全自動化硬體架構及即時影像處理。該些方法包括回饋控制、影像校準、及廣域最佳化。該端到端程序時間能經減少,並能增加每樣本的產出量。該設備能簡單的被研究科學家在工業情境中使用。控制度量能在下游分析之前確保資料品質;早期偵測減少了樣本的浪費。演算法能致能對大成像區域的可規模化能力,使其有統一的及一致的空間基因表現。能改善對大規模研究而言的資料準確度及可重複性;結果能輕易被解讀及在軟體套件中進一步分析。
參看第1圖,一多工螢光原位雜合法(mFISH)成像及影像處理設備100包括一流通槽110以固定一樣本10、一螢光顯微鏡120以獲得樣本10的影像、及一控制系統140以控制mFISH成像及影像處理設備100之各種組件 的操作。控制系統140能包括執行控制軟體的電腦142(例如具有記憶體、處理器、等等)。
螢光顯微鏡120包括一激發光源122,該激發光源能產生多個不同波長的激發光130。尤其,激發光源122能於不同時間產生具不同波長的窄頻寬光束。例如,能由一多波長連續波雷射系統來提供激發光源122,例如能被獨立地啟用以產生具不同波長之雷射光束的多個雷射模組122a。來自雷射模組122a的輸出能被多路傳輸至一共同光束路徑中。
螢光顯微鏡120包括一顯微鏡主體124,後者包括各種光學組件以將來自光源122的激發光導向流通槽110。例如,來自光源122的激發光能被耦合至一多模式光纖中,被一組透鏡重新聚焦及放大,接著被一核心成像組件(像是高數值孔徑(NA)物鏡136)導入樣本中。當需要切換激發通道時,能停用多個雷射模組122a中之一者且啟用另一雷射模組122a,而該些裝置之間的同步是由一或更多個微控制器144、146完成。
物鏡136(或整個顯微鏡主體124)能被安裝在耦接至一Z軸致動器的可垂直移動架座上。對Z位置的調整(例如藉由控制該Z軸致動器的一微控制器146)能致能焦點位置的精細調校。替代地,或者額外地,流通槽110(或在流通槽110中支撐樣本的載物台118,例如軸向壓電載物台)可以被Z軸致動器118b垂直移動。此一壓電載物台能允許精確的且迅速的多平面影像獲取。
將要成像的樣本10經定位在流通槽110中。流通槽110能為一腔室,其具有面積為大約2cm乘以2cm的截面區域(平行於物體或顯微鏡的影像平面)。樣本10能被支撐在流通槽內的一載物台118上,且該載物台(或整個流通槽)能為可橫向移動的,例如被一對直線致動器118a移動以允許XY運動。此允許對不同橫向地偏移的視野(FOV)中之樣本10的影像獲取。替代地,顯微鏡主體124可以被承載在一可橫向移動的載物台上。
對流通槽110的一入口經連接至一組雜合試劑來源112。能由控制系統140控制一閥114以在來源之間切換,以選擇哪種試劑112a被供給至流通槽110。各試劑包括一組不同的一或更多個寡核苷酸探針。各探針標的為感興趣的不同RNA序列,且具有不同的一組一或更多個螢光材料(例如磷),該一或更多個螢光材料被不同波長組合所激發。除了試劑112a以外,能有清洗流體112b(例如DI水)的來源。
流通槽110的出口經連接至一泵116(例如蠕動泵),其也被控制系統140控制以控制流體(例如試劑或清洗流體)通過流通槽110的流動。來自流通槽110的使用過的溶液能由泵116傳遞至一化學廢料管理次系統119。
操作上,控制系統140致使光源122發射激發光130,該激發光致使樣本10中螢光材料的螢光,例如結合至樣本中RNA的探針的螢光以及被該激發光之波長激發的探針的螢光。發射螢光132、還有向後傳播激發光(例如 從樣本、載物台、等等散射的激發光)被顯微鏡主體124的物鏡136收集。
所收集的光能被顯微鏡主體124中的多頻帶分光鏡138所過濾,以將發射螢光與向後傳播照明光分離,而發射螢光被傳遞至攝影機134。多頻帶分光鏡138能包括用於在各式各樣激發波長之下預期來自該等探針的各發射波長的一通帶。使用單一多頻帶分光鏡(相較於多個分光鏡或一個可移動分光鏡)能提供改善的系統穩定性。
攝影機134能為一高解析度(例如2048x2048像素)CMOS(例如科學CMOS)攝影機,且能經安裝於物鏡的即刻影像平面。可以有其他攝影機類型(例如CCD)。當被信號(例如來自微控制器者)觸發時,來自攝影機的影像資料能被捕捉,例如被發送至資料處理系統150。因此,攝影機134能收集來自該樣本的影像序列。
為進一步移除殘留的激發光並將激發通道之間的串音最小化,各雷射發射波長能與一對應帶通發射濾光片128a配對。各濾光片128a能有10~50nm(例如14~32nm)的波長。在一些實施方式中,濾光片比起產生自激發光的探針之螢光材料的頻寬要窄,例如若探針的螢光材料有拖曳的光譜輪廓。
濾光片經安裝在一高速濾光輪128上,該濾光輪可被一致動器127b旋轉。濾光輪128能被安裝於無限空間處以最小化成像路徑中的光學像差。在通過濾光輪128的發射濾光片之後,純淨的螢光信號能被鏡筒透鏡重新聚焦 並被攝影機134捕捉。分光鏡138能經定位在物鏡136與濾光輪128之間的光路徑中。
為促進該系統的高速、同步的操作,控制系統140能包括兩個微控制器144、146,該些微控制器經採用以協調的方式發送觸發信號(例如TTL信號)至螢光顯微鏡120的組件。第一微控制器144直接被電腦142運行,並觸發濾光輪128的致動器128b以切換發射濾光片128a於不同顏色通道。第一微控制器144也觸發第二微控制器146,該第二微控制器發送數位信號至光源122以控制何種波長的光被傳遞至樣本10。例如,第二微控制器146能發送開/關信號至光源122的個別雷射模組,以控制哪個雷射模組是啟用的,而因此控制何種波長的光被用於該激發光。在完成到新的激發通道的切換後,第二微控制器146控制用於壓電載物台的馬達以選擇成像高度。最後,第二微控制器146發送觸發信號至攝影機134以用於影像獲取。
電腦142與設備100之裝置組件之間的通訊被控制軟體所協調。此控制軟體能整合全部裝置組件的驅動器到一單一架構中,而因此能允許使用者將該成像系統當作單一儀器來操作(而不是需要分別控制許多裝置)。
該控制軟體支援顯微鏡的互動式操作以及成像結果的立即視覺化。額外地,該控制軟體能提供一程式介面,其允許使用者設計及自動化他們的成像工作流程。在腳本語言中能指定一組預設的工作流程描述檔(script)。
在一些實施方式中,控制系統140經配置(即被該控制軟體及/或該工作流程描述檔所配置)以在依以下順序(從最內迴圈到最外迴圈)的迴圈中獲得螢光影像(也簡單稱為「經收集影像」或僅是「影像」):z軸、顏色通道、橫向位置、及試劑。
這些迴圈可被下列虛擬碼表示:
Figure 111137276-A0305-02-0021-1
對於z軸迴圈,控制系統140致使載物台118逐步經過多個垂直位置。由於載物台118的垂直位置是由壓電致動器所控制,要調整位置所需的時間少而此迴圈中各步驟非常快。
首先,樣本可能足夠地厚(例如數微米)以使整個樣本可預期有多個影像平面。例如,可能出現多個細胞層,或甚至在一細胞內可有基因表現上的垂直變異。此外,對於薄的樣本,不一定可事先得知焦平面的垂直位置(例如熱漂移所致)。還有,樣本10可在流通槽110內垂直漂移。成像於多個Z軸位置處能確保涵蓋到厚的樣本中的大部分細胞,並能協助識別出薄的樣本中的最佳聚焦位置。
對於顏色通道迴圈,控制系統140致使光源122逐步經過不同波長的激發光。例如,啟用雷射模組中之一者,停用其他的雷射模組,且旋轉發射濾光輪128以將適當的濾光片帶到樣本10與攝影機134之間的光的光學路徑中。
對於橫向位置,控制系統140使光源122逐步經過不同橫向位置以獲得該樣本的不同視野(FOV)。例如,於該迴圈的各步,能驅動支撐載物台118的直線致動器以橫向地位移該載物台。在一些實施方式中,步數及橫向運動經選擇使得累加的FOV涵蓋了整個樣本10。在一些實施方式中,橫向運動經選擇使得FOV局部重疊。
對於試劑,控制系統140使設備100逐步經過多個不同的可用試劑。例如,於該迴圈的各步驟,控制系統140能控制閥114去將流通槽110連接至清洗流體112b,使泵116抽吸該清洗流體通過該細胞達一第一時間段以清洗目前的試劑,接著控制閥114去將流通槽110連接至不同的新試劑,並接著抽吸該新試劑通過該細胞達一第二時間 段長,該第二時間段足以讓新試劑中的探針結合至適當的RNA序列。由於要清洗流通槽及讓新試劑中的探針結合需要一些時間,要調整試劑所需的時間相較於調整橫向位置、顏色通道、或z軸而言是最長的。
作為結果,對於z軸、顏色通道(激發波長)、橫向FOV、及試劑之可能值的各個組合獲得一個螢光影像。因為最內迴圈具有最快的調整時間,而逐漸外圍的迴圈具有逐漸較慢的調整時間,此種配置方式提供為獲得用於這些參數之數值的組合的影像最省時的技術。
資料處理系統150經用以處理影像並決定基因表現,以產生空間轉錄體學資料。最低限度下,資料處理系統150包括一資料處理裝置152(例如一或更多個處理器,其由儲存在一電腦可讀取媒體上的軟體控制)及一區域儲存裝置154(例如非揮發性電腦可讀取媒體,其接收由攝影機134獲取的影像)。例如,資料處理裝置152能為安裝有GPU處理器或FPGA板的工作站。資料處理系統150也能藉由網路連接至遠端存儲156,例如藉由網際網路連接至雲端存儲。
在一些實施方式中,資料處理系統150隨著影像被接收而進行即時影像處理。特定言之,在資料獲取在進行中的同時,資料處理裝置152能進行影像預處理步驟,像是過濾及解卷積,上述步驟能在儲存裝置154中的影像資料上進行但不需要整個資料集合。由於過濾及解卷積是 資料處理管線中的主要瓶頸,隨著影像獲取發生來預處理能顯著地縮短離線處理時間而因此改善產出量。
第2圖圖示一資料處理方法的流程圖,在該資料處理方法中處理的進行是在已獲取全部影像之後。流程開始於該系統接收原始影像檔及支援檔案(步驟202)。尤其,該資料處理系統能從攝影機接收整個原始影像的集合,例如對於z軸、顏色通道(激發波長)、橫向FOV、及試劑之可能數值的各個組合有一影像。
額外地,該資料處理系統能接收一參考表現檔(expression file)(例如FPKM(每百萬對應次讀取每個序列之千鹼基對的片段)檔案)、資料綱目(data schema)、一或更多個染色影像(例如DAPI影像)。參考表現檔能被用以在傳統序列結果及mFISH結果之間交叉檢核。
從攝影機接收的影像檔能可選地包括元資料,其為該影像所被拍攝時的硬體參數值(像是載物台位置、像素大小、激發通道、及其他)。資料綱目提供基於硬體參數來排序該些影像的的規則,使得影像被以適當順序放入一或更多個影像堆疊中。若沒有包括元資料,資料綱目能將該些影像的順序關聯於用於產生該些影像之z軸、顏色通道、橫向FOV及試劑的數值。
染色影像將被覆蓋上轉錄體學資訊而呈現給使用者。
為了篩去品質不足的影像,經收集的影像能在更密集的處理之前受到一或更多個品質度量(步驟203)。僅有符合該些品質度量的影像被通過以用於進一步處理。此能顯著地減少該資料處理系統上的處理負荷。
為了偵測到聚焦的失敗,能對各個經收集影像決定一清晰度品質值。作為一例,清晰度品質的計算能藉由將整個影像之強度中的梯度加總:
Figure 111137276-A0305-02-0025-2
若清晰度品質值在一閾值之下,則該影像能被識別為失焦。該閾值可能是預設的經驗決定值,或者該閾值可能是基於針對具相同FOV之影像的清晰度值所計算。例如,能計算清晰度品質值的標準差,該閾值能基於該標準差來設定。
如上述,焦平面的z軸位置不一定為確切已知。失焦的影像能被篩去。在一些實施方式中,在各次z軸迴圈中,僅保留單一個影像,例如具有最佳清晰度品質的影像。在一些實施中,來自單一個z軸迴圈的多個影像可能符合清晰度品質,而該多個影像被組合(例如融接(blend))以形成單一影像。
為了偵測感興趣的區域,能對各個經收集影像決定一亮度品質值。亮度品質能被用以決定是否該影像中出現任何細胞。例如,該影像中全部像素的強度值能被加總並與一閾值做比較。若該總數少於該閾值,則此可能表示該影像中基本上沒有東西(即該影像中沒有細胞),且沒 有感興趣的資訊而該影像不需要被處理。做為另一例,能計算該影像中強度值的標準差並與一閾值做比較。若該標準差在該閾值之下,此表示該影像大致均勻而缺少表示出被激發之探針的亮斑點。再一次,此一經收集影像沒有感興趣的資訊而該影像不需要被處理。
以上兩個品質檢核為例示性的而可能使用其他的品質度量。例如,能使用測量斑點品質的品質度量。
另一品質檢核能偵測雜合間位移(inter-hybridization shift)。若在成像行程之間發生顯著的樣本漂移,則後續的對位可能是不可靠的。樣本漂移能被各種影像處理演算法(例如相位相關性)所偵測。尤其,經收集影像能與用於相同FOV但用於先前試劑的一影像作比較。若該些影像沒有足夠相關,例如不符合閾值相關性,則該經收集影像沒有感興趣的資訊而該影像不需要被處理。
接下來,各影像經處理以移除實驗性假影(步驟204)。由於各RNA分子將與不同激發通道的探針雜合多次,跨於多個通道、多回合影像堆疊的嚴格對齊對於顯現跨於整個FOV的RNA識別是有益的。移除實驗性假影的操作能包括場域攤平(field flattening)及/或色差校正。在一些實施方式中,場域攤平的進行在色差校正之前。
激發光中之照明輪廓線(illumination profile)的不完美(例如高斯或其他不均勻分佈,而非跨於FOV的均勻分佈)可能導致跨於該FOV的不均勻螢光信 號位準,而此種空間不均勻性往往在激發通道之間有所差異。此種不均勻能從照明場域傳遞至跨於FOV的螢光信號位準,而因此使空間轉錄體學分析結果有偏差。為了消除此種假影,利用一螢光標的對各顏色通道測量照明輪廓線(例如測試影像),該螢光標的在該FOV內具有均勻的螢光團密度。這些經測量的照明輪廓線接著被使用作為成像資料上的通道特定的標準化因數。例如,用於各顏色通道(即所要測量的各發射波長)的照明輪廓線能被儲存成一標準化影像,且在用於一顏色通道的一經收集影像中各像素的強度值被用於該相關聯顏色通道之該標準化影像之相應像素的強度值所除。
色差校正造成跨於不同發射通道改變的光學扭曲。換言之,由於探針被針對不同波長量測,在成像路徑中的色差(若有的話)可顯著地使下游分析結果降級。為了相對於此種可能的硬體瑕疵而維持分析結果的強大性,能在資料處理管線中提供一種基於軟體的色差校正解決方案。尤其,利用具有已知圖案(例如矩形圖案)之螢光斑點的標的來針對各顏色通道產生一測試影像。對於各顏色通道計算一轉換,該轉換將測試影像中的經測量斑點帶回成為與該已知圖案對齊。這些轉換被儲存,而當收集一樣本之影像時的正常操作中,選擇一適當轉換並基於對經收集影像的顏色通道對各經收集影像應用該轉換。
各影像經過處理以提供RNA影像斑點銳化(步驟206)。RNA影像斑點銳化能包括應用濾光片來移除細胞的背景及/或藉由點擴散函數來解卷積以銳化RNA斑點。
為了將RNA斑點區別於相對明亮的背景,對該影像應用一低通濾光片,例如應用至經場域攤平且色彩校正的影像以移除圍繞RNA斑點的細胞的背景。該經過濾的影像經進一步藉由一2-D點擴散函數(PSF)解卷積以銳化該RNA斑點,再與具有一半像素寬度的一2-D高斯核心卷積以稍微使斑點圓滑。
具有相同FOV的影像經對位以將其中的特徵(例如細胞或細胞胞器)對齊(步驟208)。為了正確地辨識影像序列中的RNA種類,在不同回合之影像中的特徵被對齊(例如)到次像素精確度。然而,由於mFISH樣本是以水相經成像且被一馬達載物台四處移動,經過數小時之久的成像程序下的樣本漂移及載物台漂移能轉變成為影像特徵的位移,若該位移維持位被解決則可能破壞轉錄體學分析。換言之,即使假設螢光顯微鏡對流通槽或支撐件有精確的可重複對齊,樣本在稍後的影像中可能不再於相同位置中,而此可能在解碼中導入錯誤或就是使無法解碼。
要將影像對位的一種習用技術是在載玻片上的載體材料內置放基準標記,例如螢光微粒。一般,樣本及基準標記微粒將大致一致地移動。能基於其尺寸及形狀來在影像中辨識這些微粒。對微粒的位置做比較允許將兩個影像對位,例如仿射轉換的計算。
然而,可能在不使用螢光微粒或類似基準標記、卻僅使用影像內的資料之下來對位影像。尤其,高強度區域應跨於相同FOV的多個影像大概位於相同位置。能被利用於影像間對位的技術包括相位相關性演算法及交互資訊(MI)演算法。MI演算法一般比相位相關性演算法更精確及強大。因此,能將基於交互資訊(MI)的漂移校正方法利用於對齊。然而,精確度的代價可能是高運算時間。
第3圖圖示用於利用交互資訊演算法來對齊具有相同FOV之經收集影像的程序,該程序保留MI對位演算法的高精確度同時減少運算負荷。具有相同FOV之經收集影像的集合能被認為是一原始影像堆疊;該堆疊中影像的順序由以上討論之資料綱目所設定。影像堆疊中的各經收集影像302可為(例如)2048x2048像素。一般而言,在具有此種大小的影像上進行交互資訊演算法可能需要過量的運算能力。經收集影像302經平均地切割成四個影像區塊304a、304b、304c、304d。假設每個經收集影像302是2048x2048像素,則該些影像區塊為1024x1024像素。各影像區塊304a、304b、304c、304d乘上一個相位因數φ,經相位偏移的影像區塊被重疊及加總以產生一複合影像,而該些複合影像的實部及虛部被加在一起以產生一加總影像306。此加總影像306具有原始的影像302的¼大小,例如1024x1024像素。
此程序針對原始影像堆疊中的各原始影像302重複,以產生具有縮小尺寸之影像306的堆疊308。交互資訊 演算法接著被利用堆疊308進行,以產生針對堆疊308中各縮小影像306的轉換。該轉換能包括平移及/或縮放(scaling)轉換。為了對齊原始堆疊中的影像,針對各縮小影像所決定的轉換被套用至原始堆疊中的相應經收集影像。此應造成具有原始尺寸(例如2048x2048像素)之經收集影像310的一堆疊312,其中用於相同FOV的影像310對齊。
此種切片及重疊的策略保留影像中的大部分特徵,而因此仍應具有高精確度卻減少了用於位移計算的暫時性影像大小,因此具有減小的運算負荷。儘管第3圖及以上討論聚焦在分割成為對應於原始影像之右上、右下、左上、及左下象限之四個區塊,但能夠使用另外的區塊個數,例如影像能被分割成為其他形狀(例如矩形而非正方形),只要讓區塊具有相同的維度。例如,原始影像可以被分割成256x2048像素的8個影像區塊。一般,影像能被分割成具有複數個行及複數個列的矩形網格,例如有相等的行數及列數。該影像能被分割成為N個影像,其中N=Z2且Z是偶數整數(對第3圖中的範例N=2)。
回到第2圖,由於步驟206的過濾及解卷積程序經獨立地實施在個別圖框(即位元)上,因此它們能進行於步驟208的對位之前。此外,由於用於對位的影像堆疊是在背景收集及銳化以顯現RNA的空間特徵,因此能改善對位精確度。為了驗證後面這個優點,來自低聚合(confluency)細胞樣本、高聚合細胞樣本、及腦組織 樣本的三個資料集合針對藉由對位優先對比於過濾優先的做法的對位精確度來做比較。對位品質的估算乃藉由跨於全部層將經對位影像堆疊加總,接著計算經加總影像的最大值及最小值。過濾優先的做法一般而言產生經攤平影像的較大的最大值及較小的最小值,此表示出對位品質的改良。
對位之後能進行對位品質檢核。若經正確地對位,各影像中的亮點應重疊使得總亮度增加。參考未經對位之影像的堆疊能被加總,類似地經對位影像310的堆疊312(見第3圖)能經加總,例如:
Figure 111137276-A0305-02-0031-3
IU(x,y)的最大值能與IR(x,y)的最大值比較。若max(IR(x,y))>max(IU(x,y))則能視為對位成功。相反地,若max(IR(x,y))
Figure 111137276-A0305-02-0031-7
max(IU(x,y)),則資料處理系統能指示出對位失敗並發出一錯誤碼。此FOV可以從接續的分析排除,及/或可能需要取得新影像。
可選地,在對位之後能針對各經收集影像計算一遮罩。簡言之,各像素的強度值與一閾值做比較。若強度值高於該閾值則遮罩中的一相應像素被設為1,而若強度值低於該閾值則該相應像素被設為0。該閾值能為實驗決定的預先決定值,或者能從影像中的強度值計算出。一般,該 遮罩能對應於該樣本內細胞的位置;細胞之間的空間不應發螢光而應具有低強度。
為產生遮罩,對於各視野,從該視野的影像序列產生一經壓縮圖框P(x,y)。經壓縮圖框P(x,y)的產生能藉由在該FOV的影像上進行乘積(product)運算,例如:
Figure 111137276-A0305-02-0032-4
其中i代表經對位影像310的堆疊312中的層編號(見第3圖)。能基於log P(x,y)的直方圖來產生閾值。
例如,第4圖圖示一經收集影像中經壓縮圖框值(P(x,y))的對數直方圖400。一般,該對數直方圖將包括強度中的一初始急遽增加(其可能是第一局部最大值的前緣)及強度中的一最終急遽減少(可能是第二局部最大值的後緣)。該直方圖能經歷峰值及/或邊緣查找(edge-finding)演算法以辨識出最靠近直方圖下端的一第一峰值402(或者若沒有明顯的峰值則是強度中急遽增加的一邊緣),以及最靠近直方圖上端的第二峰值404(或者若沒有明顯的峰值則是強度中急遽減少的一邊緣)。能基於該兩個經識別出的點402、404的強度值來設定一閾值410。例如,閾值410能為強度值的平均。替代峰值查找法,該直方圖能經受大津(Otsu)閾值方法來決定閾值410。此種技巧可能是有用的,因為直方圖不會永遠展現兩個清楚的峰值。
能從遮罩移除椒鹽效應(Salt and pepper effects)。在連續的「開啟(on)」像素區域外側的單 一的「開啟」像素,或者在連續的「開啟」像素區域內側的單一的「關閉(off)」像素,能被轉換成「關閉」或「開啟」以匹配其周圍區域。一般,在低於一閾值大小(例如50個像素)的鄰接像素群集能被指示成錯誤且被改變以匹配周圍區域。藉由使用遮罩,能保留記憶體空間且減少計算工作負擔。
第5A圖中圖示一經收集影像,而第5B圖圖示自第5A圖之影響產生的一遮罩。
回到第2圖,在將一FOV的影像對位之後,能進行空間轉錄體學分析(步驟210)。
為了使直方圖更一致,能在空間轉錄體學分析之前進行FOV標準化。在一些實施方式中,FOV標準化發生在對位之後。替代地,FOV標準化能發生在對位之前。FOV標準化可以被認為是過濾之一部分。即使經過以上討論的過濾及解卷積,在經收集影像中的一問題在於可能在最亮的像素之中有顯著的亮度差異,例如超過數個數量級。例如,參看第6A圖,在一經收集影像中的強度值之直方圖602有很長的尾端604。此能在解碼流程導入不確定性。FOV標準化能夠抵銷此種錯誤來源。
首先,影像中全部的強度值被相對於該影像中之最大強度值標準化。例如,決定最大強度值,及用該最大值去除全部的強度值,使得強度值在0及IMAX(例如1)之間變化。
接著,分析該影像中的強度值以決定包括最高強度值的一上分位數,例如99%以上之分位數(即前1%)、99.5%以上之分位數、或99.9%以上之分位數、或99.95%以上之分位數(即前0.05%)。此分位數的截止在第6A圖中以折線610表示。位於此分位數極限的強度值(例如在第6A圖之範例中於99.95%分位數處的強度值log(-2))能經決定及儲存。具有強度值在該上分位數內的全部像素被重設成具有最大強度值(例如1)。接著其餘像素的強度值被依比例調整以運行到相同的最大值。為完成此,沒有在該上分位數之中的像素的強度值除以用於該分位數極限的經儲存強度值。第6B圖圖示產生自FOV標準化的強度值直方圖620。
將參看第7圖解釋解碼。然而,暫時回到第3圖,用於一特定FOV的對齊影像能被視為包括多個影像層310的堆疊312,其各影像層為X乘以Y像素(例如2048x2048像素)。影像層的數目B依照顏色通道的數目(例如激發波長的數目)及雜合的數目(例如試劑的數目)之組合而定,例如B=N_雜合*N_通道。
轉向第7圖,在標準化之後,此影像堆疊能被當作一個2-D像素字矩陣702來運算。矩陣702能有P個列704,其中P=X*Y,以及B個行706,其中B是對於一給定FOV的堆疊中的影像個數,例如N_雜合*N_通道。各列704對應於像素中之一像素(跨於該堆疊之多個影像的相同像素),來自列704的數值提供一像素字710。各行706 提供字710中的數值之一數值,即來自該像素之影像層的強度值。如上述,該些數值能經標準化,例如在0與IMAX之間變化。不同強度值在第7圖中被表示成相應細胞的不同深淺程度。
若全部像素通過到解碼步驟,則全部P個字將被如下地處理。然而,細胞邊界外側的像素能被2-D遮罩篩去(見以上第4B圖)而不經處理。因此,能在以下分析中顯著減少計算負荷。
資料處理系統150儲存一碼簿722,該碼簿被用以解碼影像資料來辨識於特定像素處表現的基因。碼簿722包括多個參考碼字,各參考碼字關聯於一特定基因。如第7圖中顯示,碼簿722能被表示成具有G個列724及B個行726的一2D矩陣,其中G是碼字的個數,例如基因的個數(儘管可以由多個碼字代表相同的基因)。各列724對應於參考碼字730之一者,而各行706提供參考碼字730中該些數值之一,如由先前校準及測試已知基因所建立。對於各行,參考碼730中的數值可能是二元的,即「開啟」或「關閉」。例如,各值能為0抑或IMAX(例如1)。第7圖中將開啟及關閉值由相應細胞的淺色或深色陰影表示。
對於各個將要解碼的像素,可計算在像素字710與各參考碼字730之間的一距離d(p,i)。例如,像素字710與參考碼字730之間的距離能被計算成歐幾里德距離,例如該像素字中各數值與該參考碼字中的相應數值之間的差的平方的總和。此計算能被表示成:
Figure 111137276-A0305-02-0036-5
其中Ip,x是來自像素字之矩陣702的數值,而Ci,x是來自參考碼字之矩陣722的數值。能使用其他度量替代歐幾里德距離,例如差的絕對值的總和、餘弦角、相關性、等等。
一旦針對給定像素計算了用於各碼字的距離值,則決定最小的距離值,提供最小距離值的該碼字被選作最匹配的碼字。換句話說,資料處理設備決定min(d(p,1),d(p,2),…d(p,B)),並決定數值b作為提供該最小值的i的值(在1與B之間)。決定對應於該最匹配的碼字的基因(例如從將碼字關聯於基因的一查找表),而該像素被標記(tagged)為表現該基因。
回到第2圖,資料處理設備能濾除錯誤的標註(callout)。濾除錯誤標註的一種技巧是在指示出有一基因的表現之處的距離值d(p,b)大於一閾值(例如若d(p,b)>D1MAX)處去除標記。
還有另一個用於濾除錯誤標註的技巧是拒絕所計算位元比落在一閾值之下的碼字。位元比經計算成來自應當為開啟的層(從碼字決定)的影像字的強度值之平均,除以來自應當為關閉的層(仍從碼字決定)的影像字的強度值之平均。例如,針對第7圖中所示的碼字730,對應於行3~行6及行12~行14的位元為開啟的,而對應於行1~行2及行7~行11的位元是關閉的。因此,假設此碼字被選為用於像素p的最佳匹配,位元比BR能經計算成
Figure 111137276-A0305-02-0037-6
位元比BR經與一閾值THBR作比較。在一些實施方式中,閾值THBR是從先前的量測實證地決定。然而,在一些實施方式中,能基於從樣本獲得的量測結果針對一特定碼字自動地計算閾值THBR。第8A圖圖示針對被一特定碼字標記之全部像素的位元比強度值之對數直方圖802。基於此直方圖802計算一閾值804。例如,能利用一預設分位數(例如50%)來計算一初始閾值804。然而,在以下討論的最佳化期間能調整該分位數。
還有另一個用於濾除錯誤標註的技巧是拒絕所計算位元亮度落在一閾值之下的碼字。位元亮度計算成來自應當為開啟的層(從碼字決定)的影像字的強度值之平均。例如,針對第7圖中所示的碼字730,對應於行3~行6及行12~行14的位元為開啟的。因此,假設此碼字被選為用於像素p的最佳匹配,位元亮度BB能被計算成BB=[I(p,3)+I(p,4)+I(p,5)+I(p,6)+I(p,12)+I(p,13)+I(p,14)]/7
位元亮度BB經與一閾值THBB作比較。在一些實施方式中,閾值THBB是從先前的量測實證地決定。然而,在一些實施方式中,能基於從樣本獲得的量測結果針對一特定碼字自動地計算閾值THBB。第8B圖圖示針對被一特定碼字標記之全部像素的位元亮度值之對數直方圖812。基於此直方圖812計算出閾值814。例如,一預設分位數(例 如80%或90%)能被用以計算初始閾值814。然而,在以下討論的最佳化期間能調整該分位數。
回到第2圖,資料處理設備現能進行最佳化及再解碼(步驟212)。最佳化能包括對該些解碼參數的基於機器學習之最佳化,隨後藉更新後的解碼參數回到步驟210以更新空間轉錄體學分析結果。此循環能重複直到該些解碼參數穩定化為止。
對解碼參數的最佳化將使用績效函數(merit function),例如FPKM/TPM相關性、空間相關性、或信心比例(confidence ratio)。能在績效函數中作為變數包括的參數包括經用以移除細胞的背景的濾光片的形狀(例如頻率範圍的開始及結束、等等)、經用以銳化RNA斑點之點擴散函數的數據孔徑值、在FOV之標準化中使用的分位數界限Q、位元比閾值THBR、位元亮度閾值THBB(或用以決定位元比閾值THBR和位元亮度閾值THBB的分位數)、及/或一像素字能被認為匹配一碼字的最大距離D1max
此績效函數可有效地為不連續函數,習知的梯度跟隨(gradient following)演算法不足以辨識理想參數值。能使用機器學習模型來在參數值上收斂。
接著,資料處理設備能進行跨於全部FOV之參數值的統一(unification)。由於各FOV經個別處理,各場域可能經歷了不同標準化、設閾值、及過濾設定。作為結果,高對比的影像可能產生的直方圖會有變異致使在安 靜區域中的錯誤標註。統一的結果是全部FOV使用相同的參數值。此能顯著地移除來自安靜區域中之背景雜訊的標註,並能在大樣本範圍中提供清楚且未偏差的空間圖案。
要選擇將被跨於全部FOV使用的一參數值,有各種作法是可能的。一個選項是簡單地挑一個預先決定的FOV(例如第一個量測的FOV或接近樣本中心的FOV),並使用用於該預先決定的FOV的參數值。另一選項是對跨於多個FOV之參數的值取平均,接著利用該平均值。另一個選項是決定哪個FOV造成像素字與經標記的碼字之間的最佳適配。例如,在經標記碼字與用於該些碼字之像素字之間有最小的平均距離d(p,b1)的FOV能被決定而接著被選擇。
替代地,能在空間轉錄體學分析之前進行統一,例如在步驟208與210之間。在此情況中,將針對在步驟210之前所選擇參數(例如步驟204或206中使用的參數)進行統一,但不會針對步驟210及212中的參數進行。
資料處理設備現能進行拼接及分區(segmentation)(步驟214)。拼接將多個FOV合併成單一影像。能利用各種技術進行拼接。參看第9圖,一種做法是,對於在一起將形成樣本之合併影像的FOV的各列以及對於該列內的各FOV,決定各FOV的一水平位移。一旦計算了水平位移,對於FOV的各列計算一垂直位移。能基於交互相關性(例如相位相關性)來計算該水平及垂直位移。有了用於各FOV的水平及垂直位移,能生成單一合 併影像,並能基於該水平及垂直位移將基因座標轉移至合併影像。
能加入有關該合併螢光影像中之特定座標處(如從FOV中之座標及該FOV的水平及垂直位移所決定者)表現了一基因的一指示,例如當作元資料。此指示能被稱為「標註(callout)」。
染色影像(例如DAPI影像)能利用針對第9圖討論之技術來拼接在一起以產生合併染色影像。在一些實施方式中,不需要從所收集螢光影像產生合併螢光影像;一旦決定了用於各FOV的水平及垂直位移,就能計算在合併染色影像內的基因座標。染色影像能對經收集螢光影像經對位。能加入有關該合併染色影像中之特定座標處(如從FOV中之座標及該FOV的水平及垂直位移所決定者)表現了一基因的一指示,例如當作元資料,以提供標註。
拼接的影像中仍有一潛在問題。尤其,重疊區域中可能有一些基因被重複計數。為了移除重複計數,能計算在經標記為表現一基因之各像素與經標記為表現相同基因的其他附近像素之間的一距離(例如歐幾里德距離)。若該距離在一閾值之下,其中一個標註能被移除。若一個像素群集被標記為表現一基因,能使用更複雜的技術。
將合併影像(例如染色細胞的影像)分割成為對應於細胞之區域能利用各種已知技術來進行。分區常見在對影像的拼接之後進行,但能發生在對合併影像加標註之前或之後。
具有指示出基因表現位置之標註的分區影像現能被儲存並呈現給使用者(例如在視覺顯示器上)以供分析。
儘管以上討論假設各FOV中使用單一個z軸影像,但這並不是必要的。來自不同z軸位置的影像能被分別地處理;效果上不同z軸位置提供一組新的FOV。
如上述,一個問題在於資料處理(包括影像處理步驟)是運算上高需求且可能需要大量時間的。一個議題是,習知技術中全部的影像是在影像處理起始前獲得的。然而,隨著接收影像的同時進行即時影像處理是可能的。
第10圖圖示一種經修改資料處理方法,其中影像處理經即時進行。未被特定討論的步驟能以類似第2圖之方法的方式來進行。一開始,資料處理系統能接收支援檔案(步驟202a),不過此步驟可以延後到影像獲取之後。對某些檔案而言(例如資料綱目)可能需要被載入資料處理系統以用於將進行的影像處理。
該設備開始獲取新影像的程序(步驟202b),例如逐步藉由迴圈來調整z軸、顏色通道、橫向位置、及試劑,並獲取用於這些參數之各可能數值組合的影像。
隨著各影像的獲取,進行品質管控(QC)測試以決定是否該影像具有足夠的品質以被預處理(步驟203),若是,則進行對該經獲取影像的即時預處理(步驟252)。在此前後文中,「即時」指示出對一影像的處理是平行於對接續影像之獲取來進行的。即時預處理能包括移除實驗性假影(步驟204)及RNA影像斑點銳化(步驟206)。 可注意到由於這些影像預處理步驟是針對各影像個別地進行,且不依賴相同FOV中的其他影像或其他FOV的影像,這種方式中移動處理步驟不會不利地影響處理品質。一般,大量運算負荷能被移至當改變試劑的時候。因此,從載入樣本至流通槽中到接收具有標註之分區影像的整體時間經顯著地縮短,此能顯著地改善產出量及加速研究。
在一些實施方式中,影像經「即時」預處理,也就是對一影像的處理是平行於該設備對接續的影像的收集來進行的,且完成於資料處理系統處接收到接續的影像之前。
本說明書在相關於系統及電腦程式組件上利用「經配置」一詞。對於具有將被配置以進行特定操作或動作之一或更多個電腦的一系統而言,表示該系統已在其上安裝軟體、韌體、硬體、或以上之組合,該軟體、韌體、硬體、或其組合在運作中致使該系統進行該等操作或動作。對於將被配置以進行特定操作或動作之一或更多個電腦程式而言,其表示該一或更多個程式包括指令,當該等指令被資料處理設備執行時致使該設備進行該等操作或動作。
本說明書中所述的所請標的實施例與功能性操作能經實施在數位電子電路系統中、在有形地實現的電腦軟體或韌體中、在電腦硬體(包括本說明書中所揭露的結構與其結構均等物)中、或在以上之一或更多者之組合中。本說明書中所述所請標的實施例能經實施成一或更多個電 腦程式,即具電腦程式指令的一或更多個模組,該等指令經編碼在一有形的非暫時性儲存媒體上已供由資料處理設備執行(或控制該資料處理設備的操作)。電腦儲存媒體可能是機器可讀取儲存裝置、機器可讀取儲存基板、隨機或序列存取記憶體裝置、或以上之一或更多者的組合。替代地或額外地,該等程式指令能經編碼在人工產生的傳播信號上,例如機器產生的電性、光學、或電磁信號,該信號經產生以將資訊編碼以用於傳輸至適當的接收器設備,供資料處理設備執行。
「資料處理設備」一詞指的是資料處理硬體,並涵蓋用於處理資料的全部種類的設備、裝置、及機器,包括(例如)可程式化處理器、電腦、或多個處理器或電腦。該設備也可能是(或進一步包括)特殊用途邏輯電路系統,例如FPGA(現場可程式化閘陣列)或ASIC(應用特定積體電路)。除了硬體,該設備能可選擇地包括創造用於電腦程式之一執行環境的代碼,例如構成韌體、協定疊、資料庫管理系統、作業系統、或以上的一或更多者之組合的代碼。
電腦程式(也可稱為或描述為程式、軟體、軟體應用、應用、模組、軟體模組、描述檔(script)、或代碼)能以任何形式的程式語言寫成,包括編譯語言或直譯語言,或宣告式或程序式語言;且電腦程式能以任何形式部署,包括做為一獨立程式或作為模組、組件、次常式、或其他適合用於計算環境中的單元。程式可(但不一定)對 應於檔案系統中的一檔案。程式能被儲存在保存其他程式或資料的一檔案的一部分中,例如在一標記語言文件中所儲存的一或更多個描述檔中,在專屬於有問題之程式的單一檔案中,或在多個協作檔案中(例如儲存一或更多個模組、次程式、或部分代碼的檔案)。能將電腦程式部署以在一電腦或多個電腦上執行,該多個電腦位在一個地點處或跨於多個地點分散並由一資料通訊網路互連。
本說明書中所述程序及邏輯流程能由執行一或更多個電腦程式的一或更多個可程式化電腦所執行,以藉由在輸入資料上操作並產生輸出來執行功能。該等程序及邏輯流程也能由特殊用途邏輯電路系統進行,例如FPGA或ASIC,或藉由特殊用途邏輯電路系統與一或更多個經程式化之電腦的組合。
適合用於一電腦程式之執行的電腦能為基於一般用途微處理器或特殊用途微處理器或兩者的,或者是基於任何其他種類的中央處理單元。一般,中央處理單元將從唯讀記憶體或隨機存取記憶體或兩者來接收指令及資料。電腦的基本元件是用於執行指令的中央處理單元與用於儲存指令及資料的一或更多個記憶體裝置。中央處理單元及記憶體能補充有(或併入)特殊用途邏輯電路系統。一般,電腦也將包括用於儲存資料的一或更多個大量儲存裝置、或與該一或更多個大量儲存裝置操作耦合以自其接收資料或傳輸資料至該大量儲存裝置,該大量儲存裝置例如磁 性、磁光學碟片、光學碟片。然而,電腦不必然具有此類裝置。
適合用於儲存電腦程式指令及資料的電腦可讀取媒體包括所有形式的非揮發性記憶體、媒體和記憶體裝置,包括例如半導體記憶體裝置(例如EPROM、EEPROM、及快閃記憶體裝置);磁碟(例如內部硬碟或可移除式碟片);磁性光碟;及CD-ROM和DVD-ROM碟片。
為供與使用者互動之用,本說明書中所述標的之實施例能經實施在一電腦上,該電腦具有用於對使用者顯示資訊之顯示裝置(例如CRT(映像管)或LCD(液晶顯示器)監視器)、及使用者能用以向電腦提供輸入的鍵盤和指向裝置(例如滑鼠或軌跡球)。也能使用其他種類的裝置來提供與使用者的互動;例如,提供給使用者的回饋能為任意形式的感覺回饋,例如視覺回饋、聽覺回饋、或觸覺回饋;而來自使用者的輸入能以任何形式接收,包括聲音、口語、或觸覺輸入。此外,電腦能藉由發送文件至使用者使用的裝置或從該裝置接收文件來與該使用者互動;例如,藉由回應於從使用者之裝置上的網頁瀏覽器接收的請求,而向該網頁瀏覽器發送網頁。另外,電腦能藉由向一個人裝置(例如正執行傳訊應用的智慧型手機)發送文字訊息或其他形式的訊息來與使用者互動,並從該使用者接收回覆性訊息作為回報。
用於實施機器學習模型的資料處理設備也能包括(例如)特殊用途硬體加速器單元,其用於處理機器學習訓練或生產(即推論)工作量上共同的且計算密集之部分。
能利用一機器學習架構來實施及部署機器學習模型,該機器學習架構例如TensorFlow架構、Microsoft Cognitive Toolkit架構、Apache Singa架構、或Apache MXNet架構。
本說明書中所述標的之實施例能經實施在包括一後端組件(例如資料伺服器)、或包括中介軟體組件(例如應用伺服器)、或包括前端組件(例如客戶端電腦,該客戶端電腦具有圖形化使用者介面、網頁瀏覽器、或app,使用者藉由該圖形化使用者介面、該網頁瀏覽緝或該app能與本說明書中所述標的之實施互動)、或此類後端、中介軟體、或前端組件中一或更多者之任意組合之計算系統中。該系統的組件能由任何形式的數位資料通訊或其介質所互連,例如通訊網路。通訊網路的範例包括區域網路(LAN)及廣域網路(WAN),例如網際網路。
儘管本說明書包含許多特定的實施方式細節,這些不應被解讀為對任何發明之範疇或對所請標的之範疇的限制,卻應被認為是對特定於某些發明之某些實施例之特徵的描述。本說明書中在不同實施例之情境中說明的特定特徵也能結合在單一實施例中而實施。相反地,在單一實施例的情境中說明的各不同特徵也能被獨立地實施在多個實施例中或在任何適當的次組合中。此外,雖然以上可能 將特徵描述為以特定組合來作動且甚至一開始經聲明為如此,但來自一個經聲明之組合的一或更多個特徵在某些情況中能從該組合中去除,且該經聲明之組合可針對一個次組合或是一個次組合的變化型。
類似地,儘管在圖式中描繪並記載於申請專利範圍中的操作有特定順序,此不應被理解為要達成所欲效果,就必須要求此類操作必須以所圖示之特定順序或循序地進行,或者全部所描繪之操作必須被進行。在某些情況中,多工及平行處理可能是有益的。此外,上述實施例中各不同系統模組及組件的區分不應被理解為在全部實施例中均必須有此種區分,並應理解所述程式組件及系統一般能被一起整合在單一軟體產品中或包裝成為多個軟體產品。
以上已說明所請標的之特定實施例。其他實施例在以下申請專利範圍之範疇內。例如,申請專利範圍中所記動作能以不同順序進行並仍達成所欲結果。作為一例,隨附圖中描繪的程序不必然需要有所圖式的特定順序(或循序)以達成所欲結果。在一些情況中,多工及平行處理可為有益的。
100:mFISH成像及影像處理設備
10:樣本
110:流通槽
112:雜合試劑來源
112a:試劑
112b:清洗流體
114:閥
116:泵
118:載物台
118a:直線致動器
118b:Z軸致動器
119:化學廢料管理次系統
120:螢光顯微鏡
122:激發光源
122a:雷射模組
124:顯微鏡主體
128:濾光輪
128a:濾光片
128b:致動器
130:激發光
132:螢光
134:攝影機
136:物鏡
138:多頻帶分光鏡
140:控制系統
142:電腦
144:第一微控制器
146:第二微控制器
150:資料處理系統
152:資料處理裝置
154:儲存裝置
156:遠端存儲

Claims (5)

  1. 一種螢光原位雜合成像及分析系統,包含:一流通槽,用以容納一樣本,該樣本將暴露於一試劑中之螢光探針;一螢光顯微鏡,用以依序於複數個不同成像參數組合獲得收集該樣本的複數個影像;及一資料處理系統,包括:一線上預處理系統,該線上預處理系統經配置以隨著該等影像被收集而從該螢光顯微鏡依序接收該等影像,並進行即時影像預處理以移除該等影像的實驗性假影及提供RNA影像斑點銳化,及一離線處理系統,該離線處理系統經配置以在收集該複數個影像之後,對具有一相同視野之影像進行對位,並將該複數個影像中的強度值解碼以辨識出所表現的基因。
  2. 如請求項1所述之系統,其中該資料處理系統經配置以決定用於各影像的一影像品質值,且若該影像品質值不符合一品質閾值則篩除該影像以不經預處理。
  3. 如請求項1所述之系統,其中該資料處理系統經配置以決定用於各影像的複數個影像品質值,且若該等影像品質值之任意者不符合一相關聯品質閾值則篩除該影像以不經預處理。
  4. 如請求項1所述之系統,其中預處理以移除 實驗性假影包括以下操作:場域攤平及/或色差校正。
  5. 如請求項1所述之系統,其中斑點銳化包括以下操作:應用一低通濾光片及/或以一2-D點擴散函數進行解卷積。
TW111137276A 2019-12-17 2020-12-17 用於多工螢光原位雜合影像之獲取及處理的系統及方法 TWI835323B (zh)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
US201962949391P 2019-12-17 2019-12-17
US62/949,391 2019-12-17

Publications (2)

Publication Number Publication Date
TW202307420A TW202307420A (zh) 2023-02-16
TWI835323B true TWI835323B (zh) 2024-03-11

Family

ID=76316250

Family Applications (2)

Application Number Title Priority Date Filing Date
TW109144629A TWI782368B (zh) 2019-12-17 2020-12-17 用於多工螢光原位雜合影像之獲取及處理的系統及方法
TW111137276A TWI835323B (zh) 2019-12-17 2020-12-17 用於多工螢光原位雜合影像之獲取及處理的系統及方法

Family Applications Before (1)

Application Number Title Priority Date Filing Date
TW109144629A TWI782368B (zh) 2019-12-17 2020-12-17 用於多工螢光原位雜合影像之獲取及處理的系統及方法

Country Status (7)

Country Link
US (5) US11630067B2 (zh)
EP (1) EP4078151A4 (zh)
JP (1) JP7420950B2 (zh)
KR (1) KR20220110589A (zh)
CN (1) CN114868010A (zh)
TW (2) TWI782368B (zh)
WO (1) WO2021127019A1 (zh)

Families Citing this family (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9928403B2 (en) * 2016-02-09 2018-03-27 Molecular Devices, Llc System and method for image analysis of multi-dimensional data
EP3644044B1 (en) * 2018-10-24 2020-12-23 Leica Biosystems Imaging, Inc. Camera exposure control when acquiring fluorescence in situ hybridization images
JP7280107B2 (ja) * 2019-05-10 2023-05-23 株式会社エビデント 画像処理方法、プログラム、画像処理装置、画像処理システム、及び、顕微鏡システム
WO2021127019A1 (en) 2019-12-17 2021-06-24 Applied Materials, Inc. System and method for acquisition and processing of multiplexed fluorescence in-situ hybridization images
US20220083832A1 (en) 2020-09-16 2022-03-17 10X Genomics, Inc. Methods and systems for barcode-assisted image registration and alignment
JP7554083B2 (ja) * 2020-09-18 2024-09-19 株式会社Screenホールディングス 画像処理方法、画像処理装置および撮像装置
DE102021115588A1 (de) * 2021-06-16 2022-12-22 Carl Zeiss Meditec Ag Multi-task-lernen von weisslichtaufnahmen für ein operationsmikroskop
DE112021008418T5 (de) * 2021-11-02 2024-08-22 Leica Microsystems Cms Gmbh Verfahren zur Bereitstellung von Positionsinformationen zum Auffinden einer Zielposition in einer mikroskopischen Probe, Verfahren zur Untersuchung und/oder Verarbeitung einer solchen Zielposition und Mittel zur Durchführung dieser Verfahren
WO2023108139A2 (en) 2021-12-10 2023-06-15 10X Genomics, Inc. Multi-resolution in situ decoding
WO2023172915A1 (en) 2022-03-08 2023-09-14 10X Genomics, Inc. In situ code design methods for minimizing optical crowding
WO2023192302A1 (en) 2022-03-29 2023-10-05 10X Genomics, Inc. Spectral unmixing combined with decoding for super-multiplexed in situ analysis
KR102483745B1 (ko) * 2022-04-06 2023-01-04 주식회사 포트래이 공간전사체정보 분석장치 및 이를 이용한 분석방법
CN115602245B (zh) * 2022-09-09 2023-10-03 郑州思昆生物工程有限公司 一种荧光图筛选的方法、装置、设备以及存储介质
CN116798523B (zh) * 2023-06-01 2024-07-30 南京金域医学检验所有限公司 一种抗中性粒细胞胞浆抗体的图谱识别判定系统

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20050244863A1 (en) * 2002-09-19 2005-11-03 The Chancellor, Master And Scholars Of The University Of Oxford Molecular arrays and single molecule detection
CN1754095A (zh) * 2002-12-23 2006-03-29 爱慕斯塔尔图像与模拟:战略、分析与实现公司 用于生物芯片的芯片读出器及其相关方法
CN101553727A (zh) * 2006-08-04 2009-10-07 伊康尼西斯公司 用于显微系统的图象处理方法
CN102719357A (zh) * 2012-05-31 2012-10-10 博奥生物有限公司 一种实时监控微阵列芯片分析过程的杂交系统
TWI654309B (zh) * 2016-11-08 2019-03-21 台達電子國際(新加坡)私人有限公司 多通道螢光偵測裝置

Family Cites Families (52)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP3224860B2 (ja) 1992-06-19 2001-11-05 ポーラ化成工業株式会社 顕微鏡画像評価システム
US6007994A (en) * 1995-12-22 1999-12-28 Yale University Multiparametric fluorescence in situ hybridization
US5827660A (en) * 1996-08-09 1998-10-27 University Of Massachusetts Caged fluorochrome-labeled probes and subtraction of autoflourescence background
US6404906B2 (en) 1997-03-03 2002-06-11 Bacus Research Laboratories,Inc. Method and apparatus for acquiring and reconstructing magnified specimen images from a computer-controlled microscope
US5880473A (en) * 1997-07-28 1999-03-09 Applied Imaging, Inc. Multifluor-fluorescence in-situ hybridization (M-FISH) imaging techniques using multiple multiband filters with image registration
US5834203A (en) * 1997-08-25 1998-11-10 Applied Spectral Imaging Method for classification of pixels into groups according to their spectra using a plurality of wide band filters and hardwire therefore
US6900435B1 (en) * 1999-02-14 2005-05-31 Nanoptics, Inc. Deconvolving far-field images using scanned probe data
US6834122B2 (en) * 2000-01-22 2004-12-21 Kairos Scientific, Inc. Visualization and processing of multidimensional data using prefiltering and sorting criteria
US6995841B2 (en) * 2001-08-28 2006-02-07 Rice University Pulsed-multiline excitation for color-blind fluorescence detection
JP2003088399A (ja) 2001-09-14 2003-03-25 Fuji Photo Film Co Ltd ハイブリダイゼーション方法および装置ならびに生化学解析用ユニット
JP3902939B2 (ja) 2001-10-26 2007-04-11 株式会社日立ハイテクノロジーズ 標本中の微小領域測定装置及び方法
JP4497923B2 (ja) 2001-12-05 2010-07-07 ザ ジェイ. デビッド グラッドストーン インスティテューツ ロボット顕微鏡検査システム
KR100437253B1 (ko) 2002-05-13 2004-06-23 학교법인 한양학원 마이크로어레이 모사 이미지 생성 시스템 및 그 방법
CA2492071A1 (en) 2002-06-18 2003-12-24 Lifespan Biosciences, Inc. Computerized image capture of structures of interest within a tissue sample
JP4111099B2 (ja) 2003-08-21 2008-07-02 ソニー株式会社 生化学分析方法及び装置
US20050265588A1 (en) * 2004-02-03 2005-12-01 Bioimagene, Inc. Method and system for digital image based flourescent in situ hybridization (FISH) analysis
US20080232658A1 (en) * 2005-01-11 2008-09-25 Kiminobu Sugaya Interactive Multiple Gene Expression Map System
US20070031043A1 (en) 2005-08-02 2007-02-08 Perz Cynthia B System for and method of intelligently directed segmentation analysis for automated microscope systems
SG10201405158QA (en) 2006-02-24 2014-10-30 Callida Genomics Inc High throughput genome sequencing on dna arrays
US7826977B2 (en) 2006-08-25 2010-11-02 The Trustees Of Columbia University In The City Of New York Systems and methods for high-speed image scanning
US20120082978A1 (en) * 2006-09-15 2012-04-05 Linda Pilarski Cell Analysis On Microfluidic Chips
US8244021B2 (en) * 2006-12-20 2012-08-14 Ventana Medical Systems, Inc. Quantitative, multispectral image analysis of tissue specimens stained with quantum dots
CN101910399B (zh) * 2007-10-30 2015-11-25 考利达基因组股份有限公司 用于核酸高通量测序的装置
US20100034444A1 (en) * 2008-08-07 2010-02-11 Helicos Biosciences Corporation Image analysis
US8488863B2 (en) * 2008-11-06 2013-07-16 Los Alamos National Security, Llc Combinational pixel-by-pixel and object-level classifying, segmenting, and agglomerating in performing quantitative image analysis that distinguishes between healthy non-cancerous and cancerous cell nuclei and delineates nuclear, cytoplasm, and stromal material objects from stained biological tissue materials
US20100330555A1 (en) * 2009-06-25 2010-12-30 General Electric Company Accuracy fluorescence in-situ hybridization assay of samples with apoptotic cells
US10353190B2 (en) * 2009-12-30 2019-07-16 Koninklijke Philips N.V. Sensor for microscopy
JP5497587B2 (ja) 2010-03-23 2014-05-21 株式会社日立ハイテクノロジーズ マイクロ流路チップ及びマイクロアレイチップ
ES2490965T3 (es) 2011-04-29 2014-09-04 F. Hoffmann-La Roche Ag Procedimiento para la operación de una celda de trabajo para muestras automatizadas
EP2831591B1 (en) * 2012-03-30 2019-05-08 GE Healthcare Bio-Sciences Corp. Immunoflu0rescence and fluorescent-based nucleic acid analysis on a single sample
AU2013295679A1 (en) 2012-07-25 2015-01-29 Theranos, Inc. Image analysis and measurement of biological samples
JP5663541B2 (ja) 2012-09-19 2015-02-04 株式会社日立ハイテクノロジーズ 反応容器,並列処理装置、及びシーケンサ
US9064304B2 (en) 2013-03-18 2015-06-23 General Electric Company Image quality assessment of microscopy images
US10168259B2 (en) * 2014-06-12 2019-01-01 University Of Notre Dame Microfluidic devices, systems, and methods for imaging tissue samples
US10240146B2 (en) * 2014-07-30 2019-03-26 President And Fellows Of Harvard College Probe library construction
EP3227830B8 (en) * 2014-12-03 2020-07-01 Ventana Medical Systems, Inc. Methods, systems, and apparatuses for quantitative analysis of heterogeneous biomarker distribution
US10078778B2 (en) 2015-01-15 2018-09-18 Massachusetts Institute Of Technology Systems, methods, and apparatus for in vitro single-cell identification and recovery
WO2016120440A1 (en) * 2015-01-29 2016-08-04 Ventana Medical Systems, Inc. Dot detection, color classification of dots and counting of color classified dots
WO2016120442A1 (en) * 2015-01-30 2016-08-04 Ventana Medical Systems, Inc. Foreground segmentation and nucleus ranking for scoring dual ish images
EP3252452A1 (en) * 2016-05-25 2017-12-06 The Board of Trustees of the Leland Stanford Junior University Method for imaging and analysis of a biological specimen
US10203491B2 (en) 2016-08-01 2019-02-12 Verily Life Sciences Llc Pathology data capture
EP4339956A3 (en) 2016-11-08 2024-06-19 President and Fellows of Harvard College Multiplexed imaging using merfish, expansion microscopy, and related technologies
JP6968177B2 (ja) * 2016-12-22 2021-11-17 ベンタナ メディカル システムズ, インコーポレイテッド 一次染色および免疫組織化学画像に基づくコンピュータ採点
US10573003B2 (en) 2017-02-13 2020-02-25 Amit Sethi Systems and methods for computational pathology using points-of-interest
CN110476101B (zh) 2017-06-13 2022-02-18 谷歌有限责任公司 用于病理学的增强现实显微镜
US10417816B2 (en) * 2017-06-16 2019-09-17 Nauto, Inc. System and method for digital environment reconstruction
US11222415B2 (en) * 2018-04-26 2022-01-11 The Regents Of The University Of California Systems and methods for deep learning microscopy
WO2019219651A1 (en) * 2018-05-15 2019-11-21 Ventana Medical Systems, Inc. Quantitation of signal in stain aggregates
EP3830742A1 (en) * 2018-07-27 2021-06-09 Ventana Medical Systems, Inc. Systems for automated in situ hybridization analysis
CN112823376A (zh) * 2018-10-15 2021-05-18 文塔纳医疗系统公司 图像增强以实现改善的核检测和分割
EP3924722A1 (en) * 2019-02-13 2021-12-22 Ventana Medical Systems, Inc. Systems and methods for computing the contributions of autofluorescence in multichannel images
WO2021127019A1 (en) 2019-12-17 2021-06-24 Applied Materials, Inc. System and method for acquisition and processing of multiplexed fluorescence in-situ hybridization images

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20050244863A1 (en) * 2002-09-19 2005-11-03 The Chancellor, Master And Scholars Of The University Of Oxford Molecular arrays and single molecule detection
CN1754095A (zh) * 2002-12-23 2006-03-29 爱慕斯塔尔图像与模拟:战略、分析与实现公司 用于生物芯片的芯片读出器及其相关方法
CN101553727A (zh) * 2006-08-04 2009-10-07 伊康尼西斯公司 用于显微系统的图象处理方法
CN102719357A (zh) * 2012-05-31 2012-10-10 博奥生物有限公司 一种实时监控微阵列芯片分析过程的杂交系统
TWI654309B (zh) * 2016-11-08 2019-03-21 台達電子國際(新加坡)私人有限公司 多通道螢光偵測裝置

Also Published As

Publication number Publication date
US20210181111A1 (en) 2021-06-17
US20210182531A1 (en) 2021-06-17
EP4078151A1 (en) 2022-10-26
US11783916B2 (en) 2023-10-10
US20210199584A1 (en) 2021-07-01
EP4078151A4 (en) 2024-01-03
TW202307420A (zh) 2023-02-16
US11624708B2 (en) 2023-04-11
TWI782368B (zh) 2022-11-01
CN114868010A (zh) 2022-08-05
JP2023507595A (ja) 2023-02-24
US20230253073A1 (en) 2023-08-10
US20210183034A1 (en) 2021-06-17
JP7420950B2 (ja) 2024-01-23
WO2021127019A1 (en) 2021-06-24
TW202138564A (zh) 2021-10-16
US11630067B2 (en) 2023-04-18
KR20220110589A (ko) 2022-08-08

Similar Documents

Publication Publication Date Title
TWI835323B (zh) 用於多工螢光原位雜合影像之獲取及處理的系統及方法
US10115013B2 (en) Method and apparatus for single-particle localization using wavelet analysis
KR20220015368A (ko) 조직학적 이미지 내의 종양 및 수술 후 종양 마진 평가에 대한 컴퓨터 지원 검토
CA2658827C (en) System and method for co-registering multi-channel images of a tissue micro array
KR20220012214A (ko) 디지털 병리학을 위한 인공 지능 처리 시스템 및 자동화된 사전-진단 워크플로우
JP6594294B2 (ja) 顕微鏡画像の画像品質評価
JP6053327B2 (ja) 顕微鏡システム、標本画像生成方法及びプログラム
JP5996334B2 (ja) 顕微鏡システム、標本画像生成方法及びプログラム
CN105190686A (zh) 多重获取载玻片成像中的引用
JP7214756B2 (ja) ステイン集合体における信号の定量
Tang Computational optimisation of optical projection tomography for 3D image analysis
US20220064718A1 (en) Alignment beads for mfish
US20230230234A1 (en) Cell body segmentation using machine learning
Selinummi On Algorithms for Two and Three Dimensional High Throughput Light Microscopy