RU2749893C1 - Способ совмещения изображений флуоресценции, прибор и система секвенирования генов и запоминающая среда - Google Patents
Способ совмещения изображений флуоресценции, прибор и система секвенирования генов и запоминающая среда Download PDFInfo
- Publication number
- RU2749893C1 RU2749893C1 RU2020134588A RU2020134588A RU2749893C1 RU 2749893 C1 RU2749893 C1 RU 2749893C1 RU 2020134588 A RU2020134588 A RU 2020134588A RU 2020134588 A RU2020134588 A RU 2020134588A RU 2749893 C1 RU2749893 C1 RU 2749893C1
- Authority
- RU
- Russia
- Prior art keywords
- path
- pixel
- sums
- obtaining
- lines
- Prior art date
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/30—Determination of transform parameters for the alignment of images, i.e. image registration
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/30—Determination of transform parameters for the alignment of images, i.e. image registration
- G06T7/33—Determination of transform parameters for the alignment of images, i.e. image registration using feature-based methods
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/30—Determination of transform parameters for the alignment of images, i.e. image registration
- G06T7/33—Determination of transform parameters for the alignment of images, i.e. image registration using feature-based methods
- G06T7/337—Determination of transform parameters for the alignment of images, i.e. image registration using feature-based methods involving reference images or patches
-
- C—CHEMISTRY; METALLURGY
- C12—BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
- C12Q—MEASURING OR TESTING PROCESSES INVOLVING ENZYMES, NUCLEIC ACIDS OR MICROORGANISMS; COMPOSITIONS OR TEST PAPERS THEREFOR; PROCESSES OF PREPARING SUCH COMPOSITIONS; CONDITION-RESPONSIVE CONTROL IN MICROBIOLOGICAL OR ENZYMOLOGICAL PROCESSES
- C12Q1/00—Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions
- C12Q1/68—Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions involving nucleic acids
- C12Q1/6869—Methods for sequencing
-
- C—CHEMISTRY; METALLURGY
- C40—COMBINATORIAL TECHNOLOGY
- C40B—COMBINATORIAL CHEMISTRY; LIBRARIES, e.g. CHEMICAL LIBRARIES
- C40B20/00—Methods specially adapted for identifying library members
- C40B20/02—Identifying library members by their fixed physical location on a support or substrate
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/0002—Inspection of images, e.g. flaw detection
- G06T7/0012—Biomedical image inspection
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
- G06T7/11—Region-based segmentation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/70—Determining position or orientation of objects or cameras
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16B—BIOINFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR GENETIC OR PROTEIN-RELATED DATA PROCESSING IN COMPUTATIONAL MOLECULAR BIOLOGY
- G16B25/00—ICT specially adapted for hybridisation; ICT specially adapted for gene or protein expression
-
- C—CHEMISTRY; METALLURGY
- C12—BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
- C12Q—MEASURING OR TESTING PROCESSES INVOLVING ENZYMES, NUCLEIC ACIDS OR MICROORGANISMS; COMPOSITIONS OR TEST PAPERS THEREFOR; PROCESSES OF PREPARING SUCH COMPOSITIONS; CONDITION-RESPONSIVE CONTROL IN MICROBIOLOGICAL OR ENZYMOLOGICAL PROCESSES
- C12Q2563/00—Nucleic acid detection characterized by the use of physical, structural and functional properties
- C12Q2563/107—Nucleic acid detection characterized by the use of physical, structural and functional properties fluorescence
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10056—Microscopic image
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10064—Fluorescence image
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/20—Special algorithmic details
- G06T2207/20021—Dividing image into blocks, subimages or windows
-
- 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/30072—Microarray; Biochip, DNA array; Well plate
-
- 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/30241—Trajectory
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Vision & Pattern Recognition (AREA)
- General Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Health & Medical Sciences (AREA)
- Organic Chemistry (AREA)
- General Health & Medical Sciences (AREA)
- Molecular Biology (AREA)
- Proteomics, Peptides & Aminoacids (AREA)
- Medical Informatics (AREA)
- Wood Science & Technology (AREA)
- Biotechnology (AREA)
- Zoology (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Genetics & Genomics (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- Biophysics (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Radiology & Medical Imaging (AREA)
- Quality & Reliability (AREA)
- Immunology (AREA)
- Chemical Kinetics & Catalysis (AREA)
- General Chemical & Material Sciences (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Medicinal Chemistry (AREA)
- Microbiology (AREA)
- Evolutionary Biology (AREA)
- Bioinformatics & Computational Biology (AREA)
- General Engineering & Computer Science (AREA)
- Investigating, Analyzing Materials By Fluorescence Or Luminescence (AREA)
- Image Analysis (AREA)
Abstract
Изобретение относится к области вычислительной техники для совмещения флуоресцентных изображений. Технический результат заключается в повышении оптимизации работы по позиционированию и совмещению флуоресцентных групп на флуоресцентном изображении. Способ совмещения изображений флуоресценции включает получение изображения флуоресценции биочипа; выбор внутренней локальной области изображения флуоресценции; получение сумм значений пикселей во внутренней локальной области изображения флуоресценции; выбор множества первых шаблонных линий согласно шаблонному параметру; выполнение коррекции пиксельного уровня на локальной области линии пути; получение других перекрестков пути на биочипе согласно перекресткам пути пиксельного уровня и выполнение коррекции пиксельного уровня на других перекрестках пути; корректирование положения линии пути пиксельного уровня; получение положений субпиксельного уровня всех сайтов, равномерно распределенных на биочипе. 4 н. и 6 з.п. ф-лы, 13 ил.
Description
ОБЛАСТЬ ТЕХНИКИ, К КОТОРОЙ ОТНОСИТСЯ ИЗОБРЕТЕНИЕ
[0001] Объект настоящего изобретения относится к секвенированию генов, и более конкретно к способу совмещения флуоресцентных изображений, прибору для секвенирования генов, системе секвенирования генов и запоминающей среде.
ПРЕДПОСЫЛКИ ИЗОБРЕТЕНИЯ
[0002] В этом разделе раскрываются предпосылки или контекст для формулы изобретения и подробного описания настоящего изобретения. Представленное здесь описание нельзя рассматривать как известный уровень техники только потому, что оно включено в данный раздел.
[0003] Секвенирование генов относится к анализу нуклеотидной последовательности специфичных фрагментов ДНК, представляющей собой расположение оснований аденина (A), тимина (T), цитозина (C) и гуанина (G). Каждое из указанных выше четырех оснований содержит разную флуоресцентную группу, и разные флуоресцентные группы после возбуждения испускают флуоресценцию на разных длинах волн (цвета). Одним из самых часто используемых способов секвенирования в настоящее время является идентификация типа основания, подлежащего синтезу, путем идентификации длины волны флуоресценции. Таким образом, можно считывать нуклеотидную последовательность. В технологии секвенирования второго поколения используется система формирования изображений высокого разрешения, полученных при помощи микроскопа, для захвата изображений флуоресцентных молекул ДНК-наношариков (то есть DNB, «ДНК-наношарики») на биочипе (чип для секвенирования генов). Изображения флуоресцентных молекул отправляют на программное обеспечение для распознавания оснований. Путем декодирования сигналов изображений можно получить нуклеотидную последовательность. Во время самого процесса секвенирования, если имеется несколько изображений одной сцены, изображения одной сцены необходимо выравнивать друг с другом посредством способа позиционирования и совмещения. Затем посредством алгоритмов извлекают сигналы пятна, после чего анализируют и обрабатывают информацию о яркости, чтобы получить нуклеотидную последовательность. С развитием технологии секвенирования второго поколения приборы для секвенирования часто включают программное обеспечение для обработки и анализа данных секвенирования в реальном времени. Большинство приборов для секвенирования оснащены алгоритмами совмещения и позиционирования.
[0004] Большинство существующих технологий совмещения сопоставляют сходства на основании присущих признаков изображений флуоресценции и выполняют извлечение и совмещение для разных целевых признаков. Однако сигналы флуоресцентных молекул в изображениях высокого разрешения, полученных при помощи микроскопа, представляют собой сигналы точечных источников света. В большинстве случаев применяют способ местного центра тяжести, который извлекает значения пикселей, находящихся в окрестности каждого сайта, для вычисления центра тяжести. Однако имеется большая плотность флуоресцентных сигналов (базовых сигналов) технологии секвенирования второго поколения, и некоторые целевые сайты могут не быть достаточно яркими, что вызывает проблемы для существующих алгоритмов позиционирования.
СУЩНОСТЬ ИЗОБРЕТЕНИЯ
[0005] Таким образом, для оптимизации работы по позиционированию и совмещению флуоресцентных групп на флуоресцентном изображении необходимы способ совмещения изображений флуоресценции, прибор для секвенирования генов, система секвенирования генов и запоминающая среда.
[0006] Один аспект настоящего изобретения предусматривает способ совмещения изображений флуоресценции, применяемый к биочипу. Пиксельное расстояние между линиями пути на биочипе устанавливают в соответствии с шаблонным параметром. Способ совмещения изображений флуоресценции включает получение по меньшей мере одного изображения флуоресценции биочипа. Выбирают внутреннюю локальную область изображения флуоресценции. Получают суммы значений пикселей во внутренней локальной области изображения флуоресценции соответственно вдоль первого направления и второго направления, при этом первое направление перпендикулярно второму направлению. Множество первых шаблонных линий выбирают в соответствии с шаблонным параметром, суммы значений пикселей соответственно вдоль первого направления и второго направления обходят множеством первых шаблонных линий, чтобы найти минимальное общее значение сумм значений пикселей, соответствующих множеству первых шаблонных линий, положение минимального общего значения соответствует линии пути. Коррекцию пиксельного уровня выполняют на локальной области линии пути, и перекресток пути линий пути после коррекции пиксельного уровня представляет собой перекресток пути пиксельного уровня. Другие перекрестки пути на биочипе получают согласно перекресткам пути пиксельного уровня, и коррекцию пиксельного уровня выполняют на других перекрестках пути. Положение линии пути пиксельного уровня корректируют посредством способа центра тяжести, чтобы получить положение субпиксельного уровня линии пути. Положения субпиксельного уровня всех сайтов, равномерно распределенных на биочипе, получают посредством способа разделения сетками равного размера.
[0007] Кроме того, «получение сумм значений пикселей во внутренней локальной области изображения флуоресценции соответственно вдоль первого направления и второго направления» включает выбор множества вторых шаблонных линий. Множество вторых шаблонных линий перемещают во внутренней локальной области флуоресцентного изображения соответственно вдоль первого направления и второго направления. Вычисляют сумму значений пикселей по серой шкале во внутренней локальной области флуоресцентного изображения, покрываемых множеством вторых шаблонных линий, при этом сумма значений пикселей по серой шкале представляет собой сумму уровней яркости пикселей, покрываемых множеством вторых шаблонных линий.
[0008] Кроме того, «выбор множества первых шаблонных линий согласно шаблонному параметру, обход сумм значений пикселей соответственно вдоль первого направления и второго направления множеством первых шаблонных линий, чтобы найти минимальное общее значение сумм значений пикселей, соответствующих множеству первых шаблонных линий» включает выбор множества первых шаблонных линий согласно шаблонному параметру. Вычисляют общее значение сумм значений пикселей соответственно вдоль первого направления и второго направления, соответствующих множеству первых шаблонных линий. Получают минимальное общее значение среди общего значения сумм значений пикселей.
[0009] Кроме того, «выполнение коррекции пиксельного уровня на локальной области линии пути» включает получение сумм значений пикселей локальной области линии пути вдоль первого направления и второго направления соответственно. Выбирают множество третьих шаблонных линий, разнесенных друг от друга на заданное расстояние, чтобы найти суммы значений пикселей локальной области линии пути посредством обхода. Получают минимальное общее значение сумм значений пикселей, соответствующих множеству третьих шаблонных линий, разнесенных друг от друга на заданное расстояние. Получают положение пиксельного уровня линии пути согласно положению, соответствующему минимальному общему значению.
[0010] Кроме того, «получение положения пиксельного уровня линии пути согласно положению, соответствующему минимальному общему значению», включает получение положения пиксельного уровня впадины W-образной линии согласно положению минимального общего значения, при этом суммы значений пикселей локальной области линии пути включают W-образную линию. Положение пиксельного уровня линии пути получают согласно положению пиксельного уровня впадины.
[0011] Кроме того, «корректирование положения линии пути пиксельного уровня посредством способа центра тяжести» включает получение локальной области линии пути пиксельного уровня. Получают центр тяжести локальной области линии пути пиксельного уровня. Положение субпиксельного уровня линии пути получают в соответствии с центром тяжести.
[0012] Кроме того, «получение положений субпиксельного уровня всех сайтов, равномерно распределенных на биочипе, посредством способа разделения сетками равного размера» включает получение блочной области, образованной перекрестками пути двух соседних линий пути субпиксельного уровня вдоль первого направления и второго направления, при этом сайты расположены на блочной области согласно заданному правилу. Получают положение субпиксельного уровня всех сайтов на блочной области посредством способа разделения сетками равного размера.
[0013] Другой аспект настоящего изобретения предусматривает систему секвенирования генов, применяемую к биочипу. Пиксельное расстояние между линиями пути на биочипе устанавливают в соответствии с шаблонным параметром. Система секвенирования генов содержит модуль получения изображений, приспособленный для получения по меньшей мере одного изображения флуоресценции биочипа. Модуль выбора области приспособлен для выбора внутренней локальной области изображения флуоресценции. Модуль получения суммы значений пикселей приспособлен для получения сумм значений пикселей во внутренней локальной области изображения флуоресценции соответственно вдоль первого направления и второго направления, причем первое направление перпендикулярно второму направлению. Модуль нахождения минимального общего значения приспособлен для выбора множества первых шаблонных линий согласно шаблонному параметру, обхода сумм значений пикселей соответственно вдоль первого направления и второго направления множеством первых шаблонных линий, чтобы найти минимальное общее значение сумм значений пикселей, соответствующих множеству первых шаблонных линий, причем положение минимального общего значения соответствует линии пути. Модуль коррекции пиксельного уровня приспособлен для выполнения коррекции пиксельного уровня на локальной области линии пути, причем перекресток пути линий пути после коррекции пиксельного уровня представляет собой перекресток пути пиксельного уровня. Модуль получения других перекрестков пути приспособлен для получения других перекрестков пути на биочипе согласно перекресткам пути пиксельного уровня и выполнения коррекции пиксельного уровня на других перекрестках пути. Модуль коррекции центра тяжести приспособлен для коррекции положения линии пути пиксельного уровня посредством способа центра тяжести, чтобы получить положение субпиксельного уровня линии пути. Модуль получения сайта субпиксельного уровня приспособлен для получения положений субпиксельного уровня всех сайтов, равномерно распределенных на биочипе, посредством способа разделения сетками равного размера.
[0014] Другой аспект настоящего изобретения предусматривает прибор для секвенирования генов, содержащий процессор и запоминающее устройство, в котором хранится одна или более компьютерных программ, которые при исполнении процессором вызывают выполнение процессором указанного выше способа совмещения изображений флуоресценции.
[0015] Другой аспект настоящего изобретения предусматривает энергонезависимую запоминающую среду, в которой хранятся компьютерные команды. Когда компьютерные команды исполняются процессором, процессор приспособлен для выполнения указанного выше способа совмещения изображений флуоресценции.
КРАТКОЕ ОПИСАНИЕ ГРАФИЧЕСКИХ МАТЕРИАЛОВ
[0016] На фиг. 1 представлена блок-схема варианта осуществления способа совмещения флуоресцентных изображений согласно настоящему изобретению.
[0017] На фиг. 2 представлена структурная схема варианта осуществления прибора для секвенирования генов согласно настоящему изобретению.
[0018] На фиг. 3 представлена примерная структурная схема прибора секвенирования генов, представленной на фиг. 2.
[0019] На фиг. 4A представлена схема варианта осуществления части изображения флуоресценции согласно настоящему изобретению.
[0020] На фиг. 4B представлена увеличенная схема части, представленной на фиг. 4A.
[0021] На фиг. 5A представлена схема варианта осуществления внутренней локальной области изображения флуоресценции согласно настоящему изобретению.
[0022] На фиг. 5B представлен график, на котором показаны суммы значений пикселей вдоль вертикального направления внутренней локальной области, представленной на фиг. 5A.
[0023] На фиг. 6A представлен график, на котором показаны три первых шаблонных линии, осуществляющих поиск сумм значений пикселей, представленных на фиг. 5B.
[0024] На фиг. 6B представлен другой график, на котором показаны три первых шаблонных линии, осуществляющих поиск сумм значений пикселей, представленных на фиг. 5B.
[0025] На фиг. 6C представлен еще один график, на котором показаны три первых шаблонных линии, осуществляющих поиск сумм значений пикселей, представленных на фиг. 5B.
[0026] На фиг. 7 представлен график, на котором показано общее значение сумм значений пикселей, соответствующих трем первым шаблонным линиям.
[0027] На фиг. 8A представлена схема части линии пути согласно настоящему изобретению.
[0028] На фиг. 8B представлен график, на котором показаны суммы значений пикселей части линии пути вдоль вертикального направления.
[0029] На фиг. 8C представлен увеличенный график выбранной области, представленной на фиг. 8B.
[0030] На фиг. 9 представлен график, на котором показано общее значение, полученное суммами значений пикселей, представленными на фиг. 8B.
[0031] На фиг. 10 представлена схема, на которой показаны другие перекрестки пути, полученные посредством перекрестков пути пиксельного уровня согласно настоящему изобретению.
[0032] На фиг. 11A представлена схема части линии пути пиксельного уровня.
[0033] На фиг. 11B представлен график, на котором показаны суммы значений пикселей части линии пути пиксельного уровня вдоль горизонтального направления.
[0034] На фиг. 12A представлена схема положений линии пути пиксельного уровня.
[0035] На фиг. 12B представлена схема положений линии пути субпиксельного уровня.
[0036] На фиг. 13 представлена схема, на которой показаны субпиксельные положения сайтов, полученные посредством способа разделения сетками равного размера.
[0037] Символьное описание основных компонентов:
1 Прибор для секвенирования генов
10 Запоминающее устройство
20 Дисплей
30 Процессор
11 Модуль получения изображений
12 Модуль выбора области
13 Модуль получения сумм значений пикселей
14 Модуль нахождения минимального общего значения
15 Модуль коррекции пиксельного уровня
16 Модуль получения других перекрестков пути
17 Модуль коррекции центра тяжести
18 Модуль получения сайта субпиксельного уровня
[0038] Теперь будут описаны варианты реализации настоящей технологии только в качестве примера со ссылкой на прилагаемые фигуры.
ПОДРОБНОЕ ОПИСАНИЕ
[0039] Для того, чтобы иметь возможность понять объект изобретения, признаки и преимущества вариантов осуществления настоящего изобретения, варианты реализации настоящего изобретения теперь будут описаны только с помощью вариантов осуществления со ссылкой на графические материалы. Следует отметить, что непротиворечивые детали и признаки в вариантах осуществления настоящего изобретения могут быть объединены друг с другом.
[0040] В следующем описании особые детали объясняются, чтобы сделать понятными варианты осуществления настоящего изобретения. Описанные варианты осуществления являются не всеми вариантами осуществления настоящего изобретения, а лишь их частью. На основании вариантов осуществления настоящего изобретения другие варианты осуществления, полученные специалистом в данной области техники без творческих усилий, должны попадать в рамки объема настоящего изобретения.
[0041] Если не указано иное, все технические термины, используемые в настоящем документе, имеют такое же значение, которое обычно понимает специалист в данной области техники. Технические термины, используемые в данном документе, не следует рассматривать как ограничивающие объем вариантов осуществления.
[0042] На фиг. 1 представлена блок-схема варианта осуществления способа совмещения флуоресцентных изображений согласно настоящему изобретению. Обращаясь к фиг. 1, способ совмещения изображений флуоресценции может включать следующие этапы.
[0043] S101: получают по меньшей мере одно изображение флуоресценции биочипа.
[0044] В одном варианте осуществления биочип может представлять собой чип генного секвенирования. Изображение флуоресценции может представлять собой изображение сигналов флуоресценции, захваченное во время процесса секвенирования. Во время процесса секвенирования камера микроскопа может быть использована для захвата изображения флуоресцентных сигналов биочипа. Поле обзора камеры микроскопа является малым, по существу 768,6 мкм×648 мкм. Камера микроскопа может захватывать сотни полей обзора (FOV) для биочипа. Область между двумя соседними горизонтальными линиями пути и двумя соседними вертикальными линиями пути в каждом поле обзора определяют как блок. Блок может быть внутренним блоком или внешним блоком. Каждый блок биочипа содержит некоторое количество равномерно распределенных сайтов. Сайты могут поглощать молекулы ДНК-наношариков (DNB). Молекулы ДНК-наношариков могут представлять собой продукты амплификации, содержащие фрагменты ДНК. Молекула ДНК-наношариков содержит флуоресцентные группы, когда синтезируют основание. Флуоресцентные группы при возбуждении могут испускать флуоресцентные сигналы. Часть сайтов расположены согласно заданному правилу, формируя первый набор линий пути, параллельных друг другу, вдоль первого направления и второй набор линий пути, параллельных друг другу, вдоль второго направления. Первое направление может быть горизонтальным направлением, а второе направление может быть вертикальным направлением. Пересечения между первым набором линий пути и вторым набором линий пути называют «перекрестками пути».
[0045] S102: выбирают внутреннюю локальную область изображения флуоресценции.
[0046] В одном варианте осуществления флуоресцентные группы могут быть зафиксированы на биочипе согласно заданному правилу. Благодаря особым конструкции и процессу некоторые положения на биочипе не содержат сайтов, то есть в этих положениях не сформированы флуоресцентные группы. Когда флуоресцентные группы являются яркими в случайных положениях больше чем 25 % (четыре основания, включая аденин (A), тимин (T), цитозин (C) и гуанин (G) сбалансированы), проступают несветящиеся граничные линии. Каждая проступившая граничная линия может состоять из трех положений флуоресцентных групп. Флуоресцентные группы в центральном ряду являются яркими, а флуоресцентные группы по обеим сторонам от центрального ряда не являются яркими. Флуоресцентные группы в центральном ряду формируют линию пути, а флуоресцентные группы по обеим сторонам от центрального ряда формируют темную линию. Таким образом, проступившая граничная линия может содержать линию пути и темные линии по обеим сторонам от линии пути. Искажение формирования изображения можно игнорировать в области темных линий.
[0047] В одном варианте осуществления область, составляющая 80 % ширины вдоль первого направления и 10 % длины вдоль второго направления изображения флуоресценции, может быть выбрана как внутренняя локальная область изображения флуоресценции. Внутренняя локальная область изображения флуоресценции может содержать по меньшей мере одну линию пути вдоль каждого из первого направления и второго направления.
[0048] S103: получают суммы значений пикселей во внутренней локальной области изображения флуоресценции вдоль первого направления и второго направления. Первое направление перпендикулярно второму направлению.
[0049] В этом варианте осуществления получение сумм значений пикселей во внутренней локальной области изображения флуоресценции вдоль первого направления и второго направления может включать этапы выбора некоторого количества вторых шаблонных линий, и такое количество может быть также равно только единице. Вторые шаблонные линии перемещают во внутренней локальной области флуоресцентного изображения соответственно вдоль первого направления и второго направления. Вычисляют сумму значений пикселей по серой шкале во внутренней локальной области флуоресцентного изображения, покрываемых вторыми шаблонными линиями. То есть сумма значений пикселей по серой шкале представляет собой сумму уровней яркости пикселей, покрытых вторыми шаблонными линиями. После перемещения вторых шаблонных линий на внутренней локальной области изображения флуоресценции соответственно вдоль первого направления и второго направления можно получить суммы значений пикселей во внутренней локальной области изображения флуоресценции соответственно вдоль первого направления и второго направления. Граничная линия соответствует положению наименьшего значения из сумм значений пикселей. В настоящем изобретении для простоты сумма значений пикселей по серой шкале называется «суммой значений пикселей».
[0050] S104: в соответствии с шаблонным параметром выбирают некоторое количество первых шаблонных линий. Первые шаблонные линии обходят суммы значений пикселей соответственно вдоль первого направления и второго направления, чтобы найти минимальное общее значение сумм значений пикселей, соответствующих первым шаблонным линиям. Положение минимального общего значения соответствует положению линии пути.
[0051] В этом варианте осуществления первые шаблонные линии выбирают в соответствии с шаблонным параметром, и может существовать три первых шаблонных линии. Шаблонный параметр требует, чтобы пиксельное расстояние между двумя первыми шаблонными линиями было постоянным. Пиксельное расстояние между каждыми двумя первыми шаблонными линиями может быть одинаковым или разным. Поиск сумм значений пикселей внутренней локальной области изображения флуоресценции вдоль первого направления и второго направления последовательно осуществляют посредством выбранных первых шаблонных линий, чтобы найти общее значение сумм значений пикселей, соответствующих первым шаблонным линиям. Можно понять, что когда все первые шаблонные линии находятся возле положения минимальных сумм значений пикселей, общее значение сумм пикселей, соответствующих первым шаблонным линиям, является минимальным. Положение минимального общего значения соответствует линии пути. Положение линии пути, полученное в этот момент, является приблизительным положением линии пути. Шаблонный параметр в настоящем изобретении обозначает параметр для проектирования биочипа.
[0052] S105: на локальной области каждой линии пути выполняют коррекцию пиксельного уровня. Перекресток пути линий пути после коррекции пиксельного уровня представляет собой перекресток пути пиксельного уровня.
[0053] В этом варианте осуществления суммы значений пикселей локальной области линии пути вдоль первого направления и второго направления получают соответственно. Приводя в качестве примера вертикальное направление, выбирают некоторое количество третьих шаблонных линий, разнесенных друг от друга на заданное расстояние, чтобы находить суммы значений пикселей локальной области линии пути посредством пересечения. Суммы значений пикселей локальной области линии пути включают линию, которая имеет форму «W». Две впадины на W-образной линии соответствуют темным линиям по обеим сторонам от линии пути, и сумма значений пикселей темных линий имеет наименьшее значение. Средняя вершина в W-образной линии соответствует линии пути, и сумма значений пикселей линии пути имеет наибольшее значение. Получают минимальное общее значение сумм значений пикселей, соответствующих третьим шаблонным линиям, разнесенным друг от друга на заданное расстояние. Минимальное общее значение суммы значений пикселей, соответствующих третьим шаблонным линиям, соответствует одной впадине W-образной линии. Поскольку пиксельное расстояние между впадиной и вершиной является постоянным, положение вершины можно получить согласно положению впадины. Понятно, что положение вершины соответствует положению линии пути, так что положение пиксельного уровня линии пути можно получить согласно положению, соответствующему минимальному общему значению. Перекресток пути линии пути для коррекции пиксельного уровня представляет собой перекресток пути пиксельного уровня. Следует понимать, что перекресток пути является виртуальной точкой. То есть фактический сайт не обязательно установлен в положении этой точки, и этот сайт не обязательно излучает свет. Первая шаблонная линия, вторая шаблонная линия и третья шаблонная линия на самом деле не существуют, но представляют собой виртуальные линии для облегчения описания настоящего изобретения.
[0054] S106: другие перекрестки пути на биочипе получают согласно перекресткам пути пиксельного уровня, и коррекцию пиксельного уровня выполняют на других перекрестках пути.
[0055] В этом варианте осуществления некоторые из сайтов расположены на биочипе согласно заданному правилу, формируя первый набор линий пути, параллельных друг другу, вдоль первого направления и второй набор линий пути, параллельных друг другу, вдоль второго направления. Можно понять, что расположение между первым набором линий пути и вторым набором линий пути является регулярным, и расположение перекрестков пути также является регулярным. Когда положение пиксельного уровня перекрестков пути известно, приблизительное положение других перекрестков пути на биочипе можно получить согласно правилам. Затем выполняют коррекцию пиксельного уровня на других перекрестках пути, чтобы получить положение пиксельного уровня других перекрестков пути.
[0056] S107: положение линии пути пиксельного уровня корректируют посредством способа центра тяжести, чтобы получить положение субпиксельного уровня линии пути.
[0057] В этом варианте осуществления, приводя в качестве примера вертикальное направление, получают локальную область линии пути пиксельного уровня. Область шириной 3 пикселя и длиной 50 пикселей можно выбрать как локальную область линии пути пиксельного уровня. Получают центр тяжести локальной области линии пути пиксельного уровня. Линия пути, проходящая через центр тяжести вдоль вертикального направления, представляет собой линию пути субпиксельного уровня, то есть положение субпиксельного уровня линии пути получают в соответствии с центром тяжести. Аналогично получают положение субпиксельного уровня линии пути вдоль горизонтального направления, и перекресток пути линий пути субпиксельного уровня представляет собой перекресток пути субпиксельного уровня.
[0058] S108: положения субпиксельного уровня сайтов, равномерно распределенных на биочипе, получают посредством способа разделения сетками равного размера.
[0059] В этом варианте осуществления получают блочную область, образованную перекрестками пути двух соседних линий пути субпиксельного уровня вдоль первого направления и второго направления. Сайты расположены на блочной области согласно заданному правилу. Положение субпиксельного уровня всех сайтов на блочной области может быть получено посредством способа разделения сетками равного размера.
[0060] В способе совмещения изображений флуоресценции настоящего изобретения получают по меньшей мере одно изображение флуоресценции биочипа. Выбирают внутреннюю локальную область изображения флуоресценции. Получают суммы значений пикселей во внутренней локальной области изображения флуоресценции вдоль первого направления и второго направления. В соответствии с шаблонным параметром выбирают некоторое количество первых шаблонных линий. Первые шаблонные линии обходят суммы значений пикселей вдоль первого направления и второго направления, чтобы найти минимальное общее значение сумм значений пикселей, соответствующих первым шаблонным линиям. Положение минимального общего значения соответствует положению линии пути. На локальной области каждой линии пути выполняют коррекцию пиксельного уровня. Перекресток пути линий пути после коррекции пиксельного уровня представляет собой перекресток пути пиксельного уровня. Другие перекрестки пути на биочипе получают согласно перекресткам пути пиксельного уровня, и коррекцию пиксельного уровня выполняют на других перекрестках пути. Положение линии пути пиксельного уровня корректируют посредством способа центра тяжести, чтобы получить положение субпиксельного уровня линии пути. Положения субпиксельного уровня сайтов, равномерно распределенных на биочипе, получают посредством способа разделения сетками равного размера. С помощью варианта осуществления настоящего изобретения операцию позиционирования и совмещения флуоресцентной группы на флуоресцентном изображении можно оптимизировать.
[0061] Приведенное выше описание объясняет способ, представленный в настоящем изобретении. Далее объясняется вариант осуществления прибора или устройства для секвенирования генов согласно настоящему изобретению.
[0062] В варианте осуществления настоящего изобретения предусмотрен прибор для секвенирования генов, который содержит запоминающее устройство, процессор и компьютерную программу, хранящуюся в запоминающем устройстве и исполняемую процессором. Процессор исполняет программу так, чтобы выполнять этапы способа совмещения изображений флуоресценции. Следует отметить, что прибор для секвенирования генов может содержать платформу чипа, оптическую систему и систему жидкостного контура. Платформа чипа может поддерживать биочип. Оптическая система может захватывать изображения флуоресценции. Система жидкостного контура может выполнять биохимические реакции с использованием заданных реагентов.
[0063] На фиг. 2 представлена схема варианта осуществления прибора для секвенирования генов согласно настоящему изобретению. Обращаясь к фиг. 2, прибор 1 для секвенирования генов содержит запоминающее устройство 10 для хранения системы 100 секвенирования генов. Система 100 секвенирования генов может получать по меньшей мере одно изображение флуоресценции биочипа. Выбирают внутреннюю локальную область изображения флуоресценции. Получают суммы значений пикселей во внутренней локальной области изображения флуоресценции вдоль первого направления и второго направления. В соответствии с шаблонным параметром выбирают некоторое количество первых шаблонных линий. Первые шаблонные линии обходят суммы значений пикселей вдоль первого направления и второго направления соответственно, чтобы найти минимальное общее значение сумм значений пикселей, соответствующих первым шаблонным линиям. Положение минимального общего значения соответствует линии пути. На локальной области каждой линии пути выполняют коррекцию пиксельного уровня. Перекресток пути линий пути после коррекции пиксельного уровня представляет собой перекресток пути пиксельного уровня. Другие перекрестки пути на биочипе получают согласно перекресткам пути пиксельного уровня, и коррекцию пиксельного уровня выполняют на других перекрестках пути. Положение линии пути пиксельного уровня корректируют посредством способа центра тяжести, чтобы получить положение субпиксельного уровня линии пути. Положения субпиксельного уровня сайтов, равномерно распределенных на биочипе, получают посредством способа разделения сетками равного размера. Позиционирование и совмещение флуоресцентных групп на флуоресцентных изображениях в настоящем изобретении являются оптимизированными.
[0064] В этом варианте осуществления прибор 1 для секвенирования генов может также содержать экран 20 дисплея и процессор 30. Запоминающее устройство 10 и экран 20 дисплея могут быть электрически соединены с процессором 30.
[0065] Запоминающее устройство 10 может представлять собой устройство хранения любого из разных типов для хранения различных данных. Например, запоминающее устройство 10 может представлять собой запоминающее устройство или внутренний накопитель прибора 1 для секвенирования генов. Запоминающее устройство 10 также может представлять собой карту памяти, которая может быть подключена к прибору 1 для секвенирования генов, такую как флеш-память, карта SM (Smart Media Card) или карта SD (Secure Digital Card, Secure digital card). В дополнение, запоминающее устройство 10 может содержать высокоскоростное оперативное запоминающее устройство и может также содержать энергонезависимое запоминающее устройство, такое как жесткий диск, запоминающее устройство, подключаемый жесткий диск; карту памяти с микросхемой (Smart Media Card, SMC), карту Secure Digital (SD), флеш-карту, по меньшей мере одно устройство хранения на магнитном диске, устройство флеш-памяти или другое энергозависимое твердотельное устройство хранения. Запоминающее устройство 10 может хранить данные различных типов, например, различные типы приложений, установленных на приборе 1 для секвенирования генов, и данные для настройки или данные, полученные из указанного выше способа совмещения изображений флуоресценции.
[0066] Экран 20 дисплея установлен в приборе 1 для секвенирования генов для отображения информации.
[0067] Процессор 30 может исполнять способ совмещения изображений флуоресценции и различное программное обеспечение, установленное на приборе 1 для секвенирования генов, такое как операционная система и прикладное программное обеспечение/программное обеспечение для отображения. Процессор 30 включает, но без ограничения, центральное процессорное устройство (Central Processing Unit, CPU), микроконтроллерное устройство (Micro Controller Unit, MCU) и другие устройства для интерпретации компьютерных команд и обработки данных в компьютерном программном обеспечении.
[0068] Система 100 секвенирования генов может содержать один или более модулей. Один или более модулей хранятся в запоминающем устройстве 10 прибора 1 для секвенирования генов и могут исполняться одним или более процессорами (то есть процессором 30 в этом варианте осуществления) для реализации варианта осуществления настоящего изобретения. Например, обращаясь к фиг. 3, система 100 секвенирования генов может содержать модуль 11 получения изображений, модуль 12 выбора области, модуль 13 получения суммы значений пикселей, модуль 14 нахождения минимального общего значения, модуль 15 коррекции пиксельного уровня, модуль 16 получения других перекрестков пути, модуль 17 коррекции центра тяжести и модуль 18 получения сайта субпиксельного уровня. Модули, упоминаемые в этом варианте осуществления настоящего изобретения, могут представлять собой программные сегменты, которые выполняют конкретную функцию. Модуль является больше подходящим для описания исполнения программного обеспечения в процессоре, чем программа.
[0069] Понятно, что на основе каждого варианта осуществления вышеописанного способа совмещения флуоресцентных изображений прибор 1 для секвенирования генов может содержать некоторые или все из модулей, показанных на фиг. 3. Функции каждого модуля являются следующими. Следует отметить, что те же объекты и связанные пояснения, приведенные в указанных выше вариантах осуществления способа совмещения изображений флуоресценции, также могут быть применены к следующим представлениям функций каждого модуля. Они не повторяются здесь в целях краткости и во избежание повторения.
[0070] Модуль 11 получения изображений приспособлен для получения по меньшей мере одного изображения флуоресценции биочипа.
[0071] Модуль 12 выбора области приспособлен для выбора внутренней локальной области изображения флуоресценции.
[0072] Модуль 13 получения суммы значений пикселей приспособлен для получения сумм значений пикселей во внутренней локальной области изображения флуоресценции вдоль первого направления и второго направления. Первое направление перпендикулярно второму направлению.
[0073] Модуль 14 нахождения минимального общего значения приспособлен для выбора некоторого количества первых шаблонных линий в соответствии с шаблонным параметром. Первые шаблонные линии обходят суммы значений пикселей вдоль первого направления и второго направления соответственно, чтобы найти минимальное общее значение сумм значений пикселей, соответствующих первым шаблонным линиям. Положение минимального общего значения соответствует положению линии пути.
[0074] Модуль 15 коррекции пиксельного уровня приспособлен для выполнения коррекции пиксельного уровня на локальной области каждой линии пути. Перекресток пути линий пути после коррекции пиксельного уровня представляет собой перекресток пути пиксельного уровня.
[0075] Модуль 16 получения других перекрестков пути приспособлен для получения других перекрестков пути на биочипе согласно перекресткам пути пиксельного уровня и выполнения коррекции пиксельного уровня на других перекрестках пути.
[0076] Модуль 17 коррекции центра тяжести приспособлен для коррекции положения линии пути пиксельного уровня посредством способа центра тяжести, чтобы получать положение субпиксельного уровня линии пути.
[0077] Модуль 18 получения сайта субпиксельного уровня приспособлен для получения положений субпиксельного уровня сайтов, равномерно распределенных на биочипе, посредством способа разделения сетками равного размера.
[0078] Настоящее изобретение дополнительно предусматривает вариант осуществления энергонезависимой машиночитаемой запоминающей среды, на которой хранится компьютерная программа. Этапы способа совмещения изображений флуоресценции в любом описанном выше варианте осуществления выполняются, когда компьютерная программа исполняется процессором.
[0079] Если система секвенирования генов, прибор для секвенирования генов и модули/блоки, встроенные в компьютерное оборудование, реализованы в форме функциональных блоков программного обеспечения и продаются или используются как независимые продукты, они могут храниться в машиночитаемой запоминающей среде. На основании этого часть или все этапы вышеописанного способа также могут быть выполнены посредством компьютерной программы под командованием соответствующего аппаратного обеспечения. Компьютерная программа может храниться в машиночитаемой запоминающей среде. Этапы указанного выше способа могут выполняться, когда программа исполняется процессором. Компьютерная программа содержит компьютерный программный код, который имеет форму исходного кода, объектного кода, исполняемого файла или некоторые промежуточные формы. Машиночитаемая запоминающая среда может включать любые объект или устройство, способные содержать компьютерный программный код, среду записи, U-диск, переносной жесткий диск, магнитный диск, оптический диск, компьютерное запоминающее устройство, постоянное запоминающее устройство (ROM), оперативное запоминающее устройство (RAM), электрический несущий сигнал, телекоммуникационный сигнал и носитель для распространения программного обеспечения.
[0080] Процессор может представлять собой центральное процессорное устройство (Central Processing Unit, CPU) или другой процессор общего назначения, процессор цифровой обработки сигналов (Digital Signal Processor, DSP), специализированную интегральную схему (ASIC), готовую программируемую вентильную матрицу (FPGA) или другое программируемое логическое устройство, логический элемент на дискретных компонентах, транзисторное логическое устройство или дискретный аппаратный компонент. Процессор общего назначения может представлять собой микропроцессор или любой традиционный процессор. Процессор является центром управления системы секвенирования генов и прибора для секвенирования генов. Различные интерфейсы и линии используются для соединения вместе разных частей всей системы секвенирования генов/прибора для секвенирования генов.
[0081] В запоминающем устройстве могут храниться компьютерная программа и/или модули. Процессор может запускать или исполнять компьютерную программу и/или модули, хранящиеся в запоминающем устройстве, и также вызывать данные, хранящиеся в запоминающем устройстве, чтобы выполнять функции системы секвенирования генов/прибора для секвенирования генов. Запоминающее устройство может содержать область хранения программ и область хранения данных. В области хранения программ может храниться операционная система, прикладная программа, требующаяся для по меньшей мере одной функции (такой как функция воспроизведения звука, функция отображения изображений и т. д.). В дополнение, запоминающее устройство может содержать высокоскоростное оперативное запоминающее устройство и может также содержать энергонезависимое запоминающее устройство, такое как жесткий диск, запоминающее устройство, подключаемый жесткий диск; карту памяти с микросхемой (SMC), карту Secure Digital (SD), флеш-карту, по меньшей мере одно устройство хранения на магнитном диске, устройство флеш-памяти или другое энергозависимое твердотельное устройство хранения.
[0082] Обращаясь к фиг. 4A и фиг. 4B, на фиг. 4A представлена схема варианта осуществления части изображения флуоресценции согласно настоящему изобретению. На фиг. 4B представлена увеличенная схема части, представленной на фиг. 4A. Флуоресцентные группы распределены на биочипе согласно заданному правилу. Флуоресцентные группы формируют в некоторых положениях на биочипе посредством особых схемы и процесса. Как показано на фиг. 4A, когда флуоресцентные группы являются яркими в случайных положениях больше чем 25 % (четыре основания сбалансированы), проступают несветящиеся граничные линии. Обращаясь к фиг. 4A, две граничные линии расположены вдоль горизонтального направления (на положение граничной линии указывает черная стрелка), и три граничные линии расположены вдоль вертикального направления. Граничные линии вдоль горизонтального и вертикального направлений увеличены, как показано на фиг. 4B. Приводя в качестве примера вертикальное направление, каждая проступившая граничная линия может состоять из трех положений флуоресцентных групп. Флуоресцентные группы в центральном ряду являются яркими, а флуоресцентные группы по обеим сторонам от центрального ряда не являются яркими. Флуоресцентные группы в центральном ряду формируют линию пути, а флуоресцентные группы по обеим сторонам от центрального ряда формируют темную линию. Понятно, что проступившая граничная линия может содержать линию пути и темные линии по обеим сторонам от линии пути. Ссылка на биохимический матричный чип, раскрытый в поданной согласно PCT патентной заявке PCT/US2011/050047, обеспечит лучшее понимание расположения сайтов и линий пути на биочипе.
[0083] Обращаясь к фиг. 5A и фиг. 5B, на фиг. 5A представлена схема варианта осуществления внутренней локальной области изображения флуоресценции согласно настоящему изобретению. На фиг. 5B представлен график, на котором показаны суммы значений пикселей вдоль вертикального направления внутренней локальной области, представленной на фиг. 5A. Как показано на фиг. 5A, область, составляющая 80 % ширины вдоль первого направления и 10 % длины вдоль второго направления изображения флуоресценции, выбрана как внутренняя локальная область изображения флуоресценции (то есть часть, очерченная белой прямоугольной рамкой). Во внутренней локальной области изображения флуоресценции включены четыре граничные линии вдоль вертикального направления и две граничные линии вдоль горизонтального направления. Приводя в качестве примера вертикальное направление, вычисляют суммы значений пикселей внутренней локальной области изображения флуоресценции вдоль вертикального направления. Как показано на фиг. 5B, горизонтальная ось, представленная на фиг. 5B, представляет координаты пикселей во внутренней локальной области флуоресцентного изображения на биочипе, а вертикальная ось представляет суммы значений пикселей по серой шкале. Четыре минимальных суммы значений пикселей содержатся на фиг. 5B и соответствуют четырем граничным линиям вдоль вертикального направления во внутренней локальной области изображения флуоресценции. Положение, отмеченное черным эллиптическим значком, представляет собой вторую граничную линию вдоль вертикального направления во внутренней частичной области изображения флуоресценции (в направлении слева направо), как показано на фиг. 5A.
[0084] Обращаясь к фиг. 6A, фиг. 6B и фиг. 6C, на фиг. 6A представлен график, на котором показаны три первых шаблонных линии, осуществляющих поиск сумм значений пикселей, представленных на фиг. 5B. На фиг. 6B представлен другой график, на котором показаны три первых шаблонных линии, осуществляющих поиск сумм значений пикселей, представленных на фиг. 5B. На фиг. 6C представлен еще один график, на котором показаны три первых шаблонных линии, осуществляющих поиск сумм значений пикселей, представленных на фиг. 5B. Три первых шаблонных линии последовательно производят поиск сумм значений пикселей. Расстояния между каждыми двумя первыми шаблонными линиями являются постоянными, но отличающимися друг от друга. Например, расстояние между первой шаблонной линией и второй шаблонной линией меньше, чем расстояние между второй первой шаблонной линией и третьей первой шаблонной линией. Как показано на фиг. 6A и фиг. 6C, три первых шаблонных линии в один момент времени не находятся в положениях минимальных сумм значений пикселей. Как показано на фиг. 6B, все три первых шаблонных линии находятся в положениях минимальных сумм значений пикселей (в направлении слева направо). В это время общее значение сумм значений пикселей трех первых шаблонных линий является минимальным. Положение минимального общего значения суммы значений пикселей является приблизительным положением линии пути. В отношении поиска сумм значений пикселей на фиг. 5B посредством трех первых шаблонных линий можно понять, что три виртуальных первых шаблонных линии с фиксированным промежутком между ними могут перемещаться слева направо так, чтобы проходить кривую сумм значений пикселей, представленную на фиг. 5B, при этом суммы значений пикселей считываются (то есть их обходят) тремя первыми шаблонными линиями в положениях, которые они проходят. Затем вычисляют общее значение трех сумм значений пикселей, и это общее значение добавляют на график, представленный на фиг. 7. Каждый раз, когда получают общее значение, это общее значение добавляют к вертикальной ординате, соответствующей положению левой первой шаблонной линии.
[0085] Обращаясь к фиг. 7, показано общее значение сумм значений пикселей, соответствующих трем первым шаблонным линиям. Как показано на фиг. 7, черный круг соответствует положению, в котором общее значение сумм значений пикселей трех первых шаблонных линий является минимальным. Положение минимального общего значения соответствует приблизительному положению линии пути.
[0086] Обращаясь к фиг. 8A, 8B и 8C, на фиг. 8A представлена схема части линии пути согласно настоящему изобретению. На фиг. 8B представлен график, на котором показаны суммы значений пикселей части линии пути вдоль вертикального направления. На фиг. 8C представлен увеличенный график выбранной области, представленной на фиг. 8B. Как показано на фиг. 8A, локальная область линии пути очерчена белой прямоугольной рамкой, и ее обходят для нахождения сумм значений пикселей локальной области линии пути вдоль одного направления (то есть горизонтального или вертикального направления). Приводя в качестве примера вертикальное направление, как показано на фиг. 8B, W-образная линия (отмеченная черным эллипсом) содержится в суммах значений пикселей локальной области линии пути вдоль вертикального направления. При увеличении W-образной линии в W-образной линии содержатся две впадины и одна вершина, как показано на фиг. 8C. Понятно, что две впадины соответствуют темным линиям по обеим сторонам от линии пути на биочипе, и эта сумма значений пикселей темных линий является наименьшей. Вершина на W-образной линии соответствует линии пути, и эта сумма значений пикселей линии пути является наибольшей.
[0087] Обращаясь к фиг. 9, показано общее значение, полученное суммами значений пикселей, представленными на фиг. 8B. Для сумм значений пикселей локальной области линии пути вдоль вертикального направления, как показано на фиг. 8B, две третьих шаблонных линии с заданным расстоянием, равным четырем, выбраны для обхода сумм значений пикселей локальной области линии пути. Расстояние, равное четырем, выбрано согласно шаблонному параметру биочипа. Получают суммы значений пикселей, соответствующих двум третьим шаблонным линиям, с заданным расстоянием, равным четырем. Затем получают минимальную сумму значений пикселей. Как показано на фиг. 9, минимальная сумма значений пикселей находится между (40, 60) на горизонтальной оси. Положение минимальной суммы значений пикселей соответствует положению пиксельного уровня одной впадины W-образной линии. Положение пиксельного уровня линии пути можно получить в соответствии с положением пиксельного уровня впадины.
[0088] Обращаясь к фиг. 10, показаны другие перекрестки пути, полученные посредством перекрестков пути пиксельного уровня согласно настоящему изобретению. Как показано на фиг. 10, светло-серые точки обозначают полученные перекрестки пути пиксельного уровня. Расположение первого набора линий пути и второго набора линий пути на биочипе является регулярным. Расположение перекрестков пути линий пути также является регулярным. Когда получено положение пиксельного уровня перекрестков пути, положение других перекрестков пути на биочипе можно получить согласно соответствующим правилам (темно-серые точки обозначают другие перекрестки пути). Полученные положения других перекрестков пути являются примерными положениями, и необходимо выполнить коррекцию пиксельного уровня на других перекрестках пути, чтобы получить положения пиксельного уровня других перекрестков пути.
[0089] Обращаясь к фиг. 11A и 11B, на фиг. 11A представлена схема части линии пути пиксельного уровня. На фиг. 11B представлен график, на котором показаны суммы значений пикселей части линии пути пиксельного уровня вдоль горизонтального направления. Как показано на фиг. 11A, область шириной 3 пикселя и длиной 50 пикселей можно выбрать как локальную область линии пути пиксельного уровня (как очерчено белой линией). Получают центр тяжести локальной области линии пути пиксельного уровня. Линия пути, проходящая через центр тяжести вдоль вертикального направления, является линией пути субпиксельного уровня. Как показано на фиг. 11B, горизонтальная ось представляет координаты пикселей в локальной области линии пути пиксельного уровня вдоль горизонтального направления, а вертикальная ось представляет суммы значений пикселей по серой шкале. Интервал координат горизонтальной оси составляет (0, 2). Сумма значений пикселей является максимальной, когда координата пикселя равна 1, и сумма значений пикселей равна нулю, когда координата пикселя равна 2. Когда координата пикселя равна 0, сумма значений пикселей находится в промежутке (1000, 1200). Используя локальную область линии пути пиксельного уровня, показанную на фиг. 11B, можно определить приблизительное положение центра тяжести на локальной области линии пути пиксельного уровня. Например, центр тяжести находится в самом левом положении локальной области линии пути пиксельного уровня.
[0090] Обращаясь к фиг. 12A и 12B, на фиг. 12A представлена схема положений линии пути пиксельного уровня. На фиг. 12B представлена схема положений линии пути субпиксельного уровня. Как показано на фиг. 12A, положение пиксельного уровня определенной линии пути очерчено черным квадратом. Когда пиксель подразделяют (положение линии пути пиксельного уровня корректируют посредством способа центра тяжести, чтобы получить положение субпиксельного уровня линии пути), получают положение субпиксельного уровня, как показано на фиг. 12B. Светло-серая точка обозначает положение центра тяжести на локальной области линии пути. Понятно, что прямая линия, проходящая через черную точку, является линией пути, и линия пути является линией пути субпиксельного уровня.
[0091] На фиг. 13 представлена схема, на которой показаны субпиксельные положения сайтов, полученные посредством способа разделения сетками равного размера. Обращаясь к фиг. 13, черный круг обозначает перекресток пути линии пути субпиксельного уровня вдоль первого направления и линии пути субпиксельного уровня вдоль второго направления. Положение перекрестка пути является положением субпиксельного уровня. Блочная область между перекрестками пути двух соседних линий пути субпиксельного уровня вдоль первого направления и второго направления. Сайты расположены на блочной области согласно заданному правилу. Субпиксельное положение сайта на блочной области можно получить посредством способа разделения сетками равного размера.
[0092] Следует понимать, что в различных вариантах осуществления настоящего изобретения терминал и способ могут быть реализованы различными путями. Например, система, описанная выше, является лишь иллюстративной. Разделение модулей основано только на логической функции, и другие способы описания модулей также могут быть включены в фактическую реализацию.
[0093] Для специалистов в данной области техники очевидно, что варианты осуществления настоящего изобретения не ограничены приведенными выше деталями и также могут быть реализованы в других особых формах без отступления от идеи или основных характерных признаков вариантов осуществления настоящего изобретения. Поэтому независимо от точки зрения варианты осуществления следует рассматривать как иллюстративные и неограничивающие. Объем вариантов осуществления настоящего изобретения определен прилагаемой формулой изобретения, а не приведенным выше описанием. Поэтому все изменения в пределах смысла и объема эквивалентных элементов формулы изобретения включены в варианты осуществления настоящего изобретения. Никакие ссылочные позиции в формуле изобретения не следует рассматривать как ограничивающие формулу изобретения. Множество блоков, модулей или устройств, заявленных в системе, устройстве или терминале согласно формуле изобретения, также могут быть реализованы посредством того же блока, модуля или устройства посредством программного или аппаратного обеспечения.
[0094] Хотя информация и преимущества настоящих вариантов изобретения были изложены в приведенном выше описании вместе с деталями о структурах и функциях настоящих вариантов осуществления, описание является лишь иллюстративным. В детали могут быть внесены изменения, особенно что касается формы, размера и расположения частей, в рамках принципов настоящих приведенных для примера вариантов осуществления в полной мере, обозначенной общепринятым значением терминов, в которых выражена прилагаемая формула изобретения.
Claims (45)
1. Способ совмещения изображений флуоресценции, применяемый к биочипу, где пиксельное расстояние между линиями пути на биочипе устанавливают в соответствии с шаблонным параметром, причем способ совмещения изображений флуоресценции включает:
получение по меньшей мере одного изображения флуоресценции биочипа;
выбор внутренней локальной области изображения флуоресценции;
получение сумм значений пикселей во внутренней локальной области изображения флуоресценции соответственно вдоль первого направления и второго направления, причем первое направление перпендикулярно второму направлению;
выбор множества первых шаблонных линий согласно шаблонному параметру, обход сумм значений пикселей соответственно вдоль первого направления и второго направления множеством первых шаблонных линий, чтобы найти минимальное общее значение сумм значений пикселей, соответствующих множеству первых шаблонных линий, причем положение минимального общего значения соответствует положению линии пути;
выполнение коррекции пиксельного уровня на локальной области линии пути, причем перекресток пути линий пути после коррекции пиксельного уровня представляет собой перекресток пути пиксельного уровня;
получение других перекрестков пути на биочипе согласно перекресткам пути пиксельного уровня и выполнение коррекции пиксельного уровня на других перекрестках пути;
корректирование положения линии пути пиксельного уровня посредством способа центра тяжести, чтобы получить положение субпиксельного уровня линии пути; и
получение положений субпиксельного уровня всех сайтов, равномерно распределенных на биочипе, посредством способа разделения сетками равного размера.
2. Способ совмещения изображений флуоресценции по п. 1, отличающийся тем, что получение сумм значений пикселей во внутренней локальной области изображения флуоресценции соответственно вдоль первого направления и второго направления включает:
выбор множества вторых шаблонных линий;
перемещение множества вторых шаблонных линий во внутренней локальной области флуоресцентного изображения соответственно вдоль первого направления и второго направления; и
вычисление суммы значений пикселей по серой шкале во внутренней локальной области флуоресцентного изображения, покрываемых множеством вторых шаблонных линий, при этом сумма значений пикселей по серой шкале представляет собой сумму уровней яркости пикселей, покрываемых множеством вторых шаблонных линий.
3. Способ совмещения изображений флуоресценции по п. 2, отличающийся тем, что выбор множества первых шаблонных линий согласно шаблонному параметру, обход сумм значений пикселей соответственно вдоль первого направления и второго направления множеством первых шаблонных линий, чтобы найти минимальное общее значение сумм значений пикселей, соответствующих множеству первых шаблонных линий, включает:
выбор множества первых шаблонных линий согласно шаблонному параметру;
вычисление общего значения сумм значений пикселей соответственно вдоль первого направления и второго направления, соответствующих множеству первых шаблонных линий; и
получение минимального общего значения среди общего значения сумм значений пикселей.
4. Способ совмещения изображений флуоресценции по п. 3, отличающийся тем, что выполнение коррекции пиксельного уровня на локальной области линии пути включает:
получение сумм значений пикселей локальной области линии пути вдоль первого направления и второго направления соответственно;
выбор множества третьих шаблонных линий, разнесенных друг от друга на заданное расстояние, чтобы найти суммы значений пикселей локальной области линии пути путем обхода;
получение минимального общего значения сумм значений пикселей, соответствующих множеству третьих шаблонных линий, разнесенных друг от друга на заданное расстояние; и
получение положения пиксельного уровня линии пути согласно положению, соответствующему минимальному общему значению.
5. Способ совмещения изображений флуоресценции по п. 4, отличающийся тем, что получение положения пиксельного уровня линии пути согласно положению, соответствующему минимальному общему значению, включает:
получение положения пиксельного уровня впадины W-образной линии согласно положению минимального общего значения, при этом суммы значений пикселей локальной области линии пути включают W-образную линию; и
получение положения пиксельного уровня линии пути согласно положению пиксельного уровня впадины.
6. Способ совмещения изображений флуоресценции по п. 5, отличающийся тем, что корректирование положения линии пути пиксельного уровня посредством способа центра тяжести включает:
получение локальной области линии пути пиксельного уровня;
получение центра тяжести локальной области линии пути пиксельного уровня; и
получение положения субпиксельного уровня линии пути в соответствии с центром тяжести.
7. Способ совмещения изображений флуоресценции по п. 6, отличающийся тем, что получение положений субпиксельного уровня всех сайтов, равномерно распределенных на биочипе, посредством способа разделения сетками равного размера включает:
получение блочной области, образованной перекрестками пути двух соседних линий пути субпиксельного уровня вдоль первого направления и второго направления, при этом сайты расположены на блочной области согласно заданному правилу; и
получение положения субпиксельного уровня всех сайтов на блочной области посредством способа разделения сетками равного размера.
8. Система секвенирования генов, применяемая к биочипу, где пиксельное расстояние между линиями пути на биочипе устанавливается в соответствии с шаблонным параметром, причем система секвенирования генов содержит:
модуль получения изображений, приспособленный для получения по меньшей мере одного изображения флуоресценции биочипа;
модуль выбора области, приспособленный для выбора внутренней локальной области изображения флуоресценции;
модуль получения суммы значений пикселей, приспособленный для получения сумм значений пикселей во внутренней локальной области изображения флуоресценции соответственно вдоль первого направления и второго направления, причем первое направление перпендикулярно второму направлению;
модуль нахождения минимального общего значения, приспособленный для выбора множества первых шаблонных линий согласно шаблонному параметру, обхода сумм значений пикселей соответственно вдоль первого направления и второго направления множеством первых шаблонных линий, чтобы найти минимальное общее значение сумм значений пикселей, соответствующих множеству первых шаблонных линий, причем положение минимального общего значения соответствует линии пути;
модуль коррекции пиксельного уровня, приспособленный для выполнения коррекции пиксельного уровня на локальной области линии пути, причем перекресток пути линий пути после коррекции пиксельного уровня представляет собой перекресток пути пиксельного уровня;
модуль получения других перекрестков пути, приспособленный для получения других перекрестков пути на биочипе согласно перекресткам пути пиксельного уровня и выполнения коррекции пиксельного уровня на других перекрестках пути;
модуль коррекции центра тяжести, приспособленный для коррекции положения линии пути пиксельного уровня посредством способа центра тяжести, чтобы получить положение субпиксельного уровня линии пути; и
модуль получения сайта субпиксельного уровня, приспособленный для получения положений субпиксельного уровня всех сайтов, равномерно распределенных на биочипе, посредством способа разделения сетками равного размера.
9. Прибор для секвенирования генов, содержащий:
процессор; и
запоминающее устройство, в котором хранится одна или более компьютерных программ, которые при исполнении процессором вызывают выполнение процессором способа совмещения изображений флуоресценции по любому из пп. 1-7.
10. Энергонезависимая запоминающая среда, в которой хранятся компьютерные команды, при этом, когда компьютерные команды исполняются процессором, процессор приспособлен для выполнения способа совмещения изображений флуоресценции по любому из пп. 1-7.
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
PCT/CN2018/082578 WO2019196019A1 (zh) | 2018-04-10 | 2018-04-10 | 荧光图像配准方法、基因测序仪及系统、存储介质 |
Publications (1)
Publication Number | Publication Date |
---|---|
RU2749893C1 true RU2749893C1 (ru) | 2021-06-18 |
Family
ID=68162773
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
RU2020134588A RU2749893C1 (ru) | 2018-04-10 | 2018-04-10 | Способ совмещения изображений флуоресценции, прибор и система секвенирования генов и запоминающая среда |
Country Status (11)
Country | Link |
---|---|
US (2) | US11682125B2 (ru) |
EP (1) | EP3779868B1 (ru) |
JP (1) | JP7072081B2 (ru) |
KR (1) | KR102504987B1 (ru) |
CN (1) | CN111971711B (ru) |
AU (1) | AU2018418609B2 (ru) |
BR (1) | BR112020020949A2 (ru) |
CA (1) | CA3096723C (ru) |
RU (1) | RU2749893C1 (ru) |
SG (1) | SG11202010040YA (ru) |
WO (1) | WO2019196019A1 (ru) |
Families Citing this family (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20230092006A1 (en) * | 2020-02-12 | 2023-03-23 | Mgi Tech Co., Ltd. | Optical imaging system and biochemical substance detection system using same |
CN113255649B (zh) * | 2021-06-21 | 2023-09-19 | 北博(厦门)智能科技有限公司 | 一种基于图像识别的图像分割框选方法及终端 |
WO2024001051A1 (zh) * | 2022-06-29 | 2024-01-04 | 深圳华大生命科学研究院 | 空间组学单细胞数据获取方法及装置和电子设备 |
CN115331735B (zh) * | 2022-10-11 | 2023-03-17 | 青岛百创智能制造技术有限公司 | 芯片解码方法及装置 |
WO2024140739A1 (zh) * | 2022-12-30 | 2024-07-04 | 广东润鹏生物技术有限公司 | 测序图像的光斑融合方法、装置、存储介质及计算机设备 |
WO2024202904A1 (ja) * | 2023-03-31 | 2024-10-03 | ソニーグループ株式会社 | データ処理方法、塩基配列データ生成システム、プログラム、及び核酸特定方法 |
CN117237441B (zh) * | 2023-11-10 | 2024-01-30 | 湖南科天健光电技术有限公司 | 亚像素定位方法、系统、电子设备和介质 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20150080231A1 (en) * | 2010-04-30 | 2015-03-19 | Complete Genomics, Inc. | Method and system for accurate alignment and registration of array for dna sequencing |
CN105303187A (zh) * | 2015-12-10 | 2016-02-03 | 北京中科紫鑫科技有限责任公司 | 一种dna测序的图像识别方法及装置 |
RU2015143532A (ru) * | 2013-03-13 | 2017-04-19 | Дженвив Байосайенсиз, Инк. | Системы и способы для обнаружения клеток с применением сконструированных трансдукционных частиц |
CN107730541A (zh) * | 2016-08-12 | 2018-02-23 | 广州康昕瑞基因健康科技有限公司 | 图像配准方法和系统及图像拍摄对位方法和系统 |
Family Cites Families (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR100463336B1 (ko) * | 2001-10-11 | 2004-12-23 | (주)가이아진 | 바이오칩 이미지 분석 시스템 및 그 방법 |
AU2003298655A1 (en) * | 2002-11-15 | 2004-06-15 | Bioarray Solutions, Ltd. | Analysis, secure access to, and transmission of array images |
US20050281462A1 (en) * | 2004-06-16 | 2005-12-22 | Jayati Ghosh | System and method of automated processing of multiple microarray images |
US9671344B2 (en) * | 2010-08-31 | 2017-06-06 | Complete Genomics, Inc. | High-density biochemical array chips with asynchronous tracks for alignment correction by moiré averaging |
JP5822664B2 (ja) * | 2011-11-11 | 2015-11-24 | 株式会社Pfu | 画像処理装置、直線検出方法及びコンピュータプログラム |
CN102821238B (zh) * | 2012-03-19 | 2015-07-22 | 北京泰邦天地科技有限公司 | 宽视场超高分辨率成像系统 |
CN103150550B (zh) * | 2013-02-05 | 2015-10-28 | 长安大学 | 一种基于运动轨迹分析的道路行人事件检测方法 |
AU2014296253A1 (en) * | 2013-07-30 | 2016-02-04 | President And Fellows Of Harvard College | Quantitative DNA-based imaging and super-resolution imaging |
US9536334B2 (en) * | 2013-10-18 | 2017-01-03 | Schlumberger Technology Corporation | Image strip display |
CN105427328B (zh) | 2015-12-10 | 2017-05-31 | 北京中科紫鑫科技有限责任公司 | 一种dna测序的互信息图像配准方法及装置 |
CN105550990B (zh) | 2015-12-10 | 2017-07-28 | 北京中科紫鑫科技有限责任公司 | 一种基于傅里叶配准的dna图像处理方法及装置 |
CN105427327B (zh) * | 2015-12-10 | 2017-05-31 | 北京中科紫鑫科技有限责任公司 | 一种dna测序的图像配准方法及装置 |
JP6687039B2 (ja) * | 2016-02-05 | 2020-04-22 | 株式会社リコー | 物体検出装置、機器制御システム、撮像装置、物体検出方法、及びプログラム |
-
2018
- 2018-04-10 KR KR1020207032386A patent/KR102504987B1/ko active IP Right Grant
- 2018-04-10 EP EP18914625.1A patent/EP3779868B1/en active Active
- 2018-04-10 BR BR112020020949-0A patent/BR112020020949A2/pt unknown
- 2018-04-10 AU AU2018418609A patent/AU2018418609B2/en active Active
- 2018-04-10 WO PCT/CN2018/082578 patent/WO2019196019A1/zh active Application Filing
- 2018-04-10 CN CN201880005564.0A patent/CN111971711B/zh active Active
- 2018-04-10 RU RU2020134588A patent/RU2749893C1/ru active
- 2018-04-10 JP JP2020555468A patent/JP7072081B2/ja active Active
- 2018-04-10 SG SG11202010040YA patent/SG11202010040YA/en unknown
- 2018-04-10 CA CA3096723A patent/CA3096723C/en active Active
- 2018-04-10 US US17/043,226 patent/US11682125B2/en active Active
-
2023
- 2023-05-24 US US18/201,309 patent/US12039744B2/en active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20150080231A1 (en) * | 2010-04-30 | 2015-03-19 | Complete Genomics, Inc. | Method and system for accurate alignment and registration of array for dna sequencing |
RU2015143532A (ru) * | 2013-03-13 | 2017-04-19 | Дженвив Байосайенсиз, Инк. | Системы и способы для обнаружения клеток с применением сконструированных трансдукционных частиц |
CN105303187A (zh) * | 2015-12-10 | 2016-02-03 | 北京中科紫鑫科技有限责任公司 | 一种dna测序的图像识别方法及装置 |
CN107730541A (zh) * | 2016-08-12 | 2018-02-23 | 广州康昕瑞基因健康科技有限公司 | 图像配准方法和系统及图像拍摄对位方法和系统 |
Also Published As
Publication number | Publication date |
---|---|
CN111971711A (zh) | 2020-11-20 |
WO2019196019A1 (zh) | 2019-10-17 |
EP3779868A4 (en) | 2021-11-17 |
EP3779868B1 (en) | 2024-05-08 |
CA3096723C (en) | 2023-10-24 |
BR112020020949A2 (pt) | 2021-03-02 |
JP2021517701A (ja) | 2021-07-26 |
KR102504987B1 (ko) | 2023-02-28 |
CN111971711B (zh) | 2024-08-02 |
AU2018418609A1 (en) | 2020-11-12 |
US11682125B2 (en) | 2023-06-20 |
SG11202010040YA (en) | 2020-11-27 |
US20230298188A1 (en) | 2023-09-21 |
EP3779868A1 (en) | 2021-02-17 |
US20210019895A1 (en) | 2021-01-21 |
JP7072081B2 (ja) | 2022-05-19 |
US12039744B2 (en) | 2024-07-16 |
AU2018418609B2 (en) | 2021-10-07 |
KR20200142040A (ko) | 2020-12-21 |
CA3096723A1 (en) | 2019-10-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
RU2749893C1 (ru) | Способ совмещения изображений флуоресценции, прибор и система секвенирования генов и запоминающая среда | |
US11972583B2 (en) | Fluorescence image registration method, gene sequencing instrument, and storage medium | |
US11164316B2 (en) | Image processing systems and methods for displaying multiple images of a biological specimen | |
US11960073B2 (en) | System and method for template-based image analysis | |
US7356184B2 (en) | Analysis of array/assay data utilizing transformation of scatter plot to image plot | |
CN114868010A (zh) | 用于多重荧光原位杂化图像的获取及处理的系统及方法 | |
Jung et al. | An automatic block and spot indexing with k-nearest neighbors graph for microarray image analysis | |
CN110111842A (zh) | 图像清晰度分析及对焦方法、测序仪、系统与存储介质 | |
KR20110111470A (ko) | 현미경 검사 | |
WO2024000288A1 (zh) | 图像拼接方法、基因测序系统及相应的基因测序仪 | |
Proksch et al. | Toward quantitative super-resolution microscopy: molecular maps with statistical guarantees | |
CN118154647A (zh) | 图像配准方法、基因测序仪、电子设备 | |
CN114693717A (zh) | 图像处理方法、装置、设备及计算机可读存储介质 | |
Wilson | Quantitative analysis of microscopy |