IP Library Granted Patent US 12,731,273
Granted Patent B2
US 12,731,273 · App. 18/692,758 · Granted Sep 8, 2026

Systems and methods for image registration or alignment

Inventors: Olga Vorobyova (Pleasanton, CA); Brynn Claypoole (Pleasanton, CA); Dongyao Li (Pleasanton, CA); Neil Ira Weisenfeld (Pleasanton, CA); Didem Pelin Sarikaya (Pleasanton, CA); Peigeng Li (Pleasanton, CA); Guy Joseph (Pleasanton, CA); Eric Siegel (Pleasanton, CA); Naga Sudha Kodavatikanti (Cupertino, CA)
Assignee: 10x Genomics, Inc.
G06T7/337G06T7/13G06T2207/10024G06T2207/10056G06T2207/10064G06T2207/20036G06T2207/20101G06T2207/30024
View Patent ↗
Loading inventors, assignments & file history…
Monitor This Case
Get email alerts when status or documents change.
Order Certified Copies
Most orders are placed with the USPTO same day — all within 24 business hours.
Order via The Patent Place →
Pre-filled with this patent's details
Quick Facts
Patent No.
US 12,731,273
App. No.
18/692,758
Granted
Sep 8, 2026
Kind
B2
Abstract

Systems and methods for overlaying image data for a biological sample on spatial analyte data are provided. A first image of the sample on a first substrate and a second image of the sample on the first substrate overlayed on a second substrate are obtained. The second substrate includes spatial fiducials and capture spots. At least one of the first substrate and the second substrate is transparent. A registration for the first image and the second image is determined, using a first pattern of the sample in the first image and a second pattern of the sample in the second image. The registration is used to overlay the first image onto a spatial dataset including spatial analyte data for the capture spots from the sample. A frame of reference of the spatial dataset is known with respect to the second image, based on the spatial fiducials of the second image.

Claims (60)

1 . A method for overlaying image data for a biological sample onto spatial analyte data for a plurality of analytes of the biological sample, the method comprising:

using a computer system comprising one or more processing cores and a memory:

obtaining a first image of the biological sample on a first substrate;

receiving a second image of the biological sample on the first substrate overlayed on a second substrate, wherein the second substrate comprises one or more spatial fiducials and a set of capture spots, and wherein the set of capture spots comprises at least 1000 capture spots,

wherein at least one of the first substrate and the second substrate is transparent;

determining a registration for the first image and the second image, wherein the registration registers the first image and the second image to each other, using a first pattern of the biological sample in the first image and a second pattern of the biological sample in the second image; and

using the registration to overlay the first image onto a spatial dataset comprising spatial analyte data for the set of capture spots from the biological sample, wherein a frame of reference of the spatial dataset is known with respect to the second image, based on the one or more spatial fiducials of the second image,

thereby overlaying image data for the biological sample onto spatial analyte data for the plurality of analytes of the biological sample.

2 . The method of claim 1 , wherein the first image of the biological sample is a histological image of the biological sample.

3 . The method of claim 1 , wherein the determining the registration for the first image and the second image comprises:

identifying one or more landmarks, each respective landmark comprising a respective pair of reference positions including a first respective reference position for the first pattern of the biological sample in the first image and a corresponding second respective reference position for the second pattern of the biological sample in the second image, thereby identifying a first set of landmark coordinates for the first image and a second set of landmark coordinates for the second image;

using the first set of landmark coordinates for the first image and the second set of landmark coordinates for the second image to obtain a transformation between the first set of landmark coordinates and the second set of landmark coordinates with an alignment algorithm, thereby obtaining a first alignment of the first image with the second image.

4 . The method of claim 3 , wherein the identifying one or more landmarks is performed, on a display, via manual user selection of each respective pair of reference positions in the first image and the second image.

5 . The method of claim 3 , further comprising performing, on a display, via user interaction, a manual alignment of the first image and the second image.

6 . The method of claim 5 , wherein the manual alignment is selected from the group consisting of: a translation, a scaling, and a rotation.

7 . The method of claim 3 , further comprising refining the transformation using a refinement algorithm, thereby obtaining a second alignment of the first image with the second image.

8 . The method of claim 3 , wherein a respective similarity metric for a respective alignment of the first image with the second image is determined by evaluating an intensity of a first plurality of pixels in the first image compared to corresponding pixels in a second plurality of pixels in the second image that are overlayed onto each other in the respective alignment.

9 . The method of claim 3 , wherein the one or more landmarks comprises 3 landmarks.

10 . The method of claim 3 , wherein the alignment algorithm transforms the first set of landmark coordinates for the first image relative to the second set of landmark coordinates for the second image.

11 . The method of claim 3 , wherein the alignment algorithm transforms the second set of landmark coordinates for the second image relative to the first set of landmark coordinates for the first image.

12 . The method of claim 3 , further comprising displaying, on a display, an image overlay for the first image aligned with the second image, wherein the first image is displayed in color and the second image is displayed in grayscale.

13 . The method of claim 1 , wherein the determining the registration for the first and second image comprises:

obtaining a respective first area of the first pattern and a respective first centroid for the first area from an analysis of the first pattern;

obtaining a respective second area of the second pattern and a respective second centroid for the second area from an analysis of the second pattern;

scaling at least one of the first image and the second image based on a relative proportion of the first and second area;

centering the first image with the second image based on an alignment between the first and second centroid;

sampling, for each respective rotation angle in a plurality of rotation angles, a corresponding transformation between the first and second image, by a procedure that comprises:

rotating one of the first and second image relative to the other of the first and second image by the respective rotation angle to obtain a respective candidate rotational alignment of the first and second image;

refining the respective candidate rotational alignment to obtain a respective refined candidate rotational alignment; and

scoring the respective refined candidate rotational alignment of the first and second image by determining a similarity metric across a respective intensity of a first plurality of pixels in the first image compared to corresponding pixels in a second plurality of pixels in the second image that are overlayed onto each other by the respective candidate rotational alignment of the first and second image,

thereby obtaining a plurality of refined candidate rotational alignments and a corresponding plurality of similarity metrics; and

selecting the rotational alignment from among the plurality of refined candidate rotational alignments that yields the highest similarity metric in the corresponding plurality of similarity metrics as the registration for the first and second image.

14 . The method of claim 13 , wherein the plurality of rotation angles comprises 4 or more different rotation angles, wherein each rotation angle in the plurality of rotation angles is between 0 and 2π radians.

15 . The method of claim 13 , wherein the procedure further comprises:

rotating one of a mirror image of the first image and the second image relative to the other of the mirror image of the first image and the second image by the respective rotation angle to obtain a respective candidate rotational alignment of the mirror image of the first image and the second image; and

scoring the respective candidate rotational alignment of the mirror image of the first image and the second image by determining a similarity metric for the respective candidate rotational alignment of the mirror image of the first image and the second image.

16 . The method of claim 1 , wherein:

the first image is obtained without spatial fiducials,

the determining a registration for the first image and the second image comprises overlaying the one or more spatial fiducials of the second image onto the first image, and

the using the registration to overlay the first image onto the spatial dataset comprises associating the one or more spatial fiducials of the first image to the frame of reference of the spatial dataset, thereby allowing the first image to be in the same frame of reference as the spatial dataset.

17 . The method of claim 1 , wherein the first image and the second image have different image resolutions.

18 . The method of claim 1 , wherein the biological sample is prepared for imaging on the first substrate using a stain selected from the group consisting of: live/dead stain, trypan blue, periodic acid-Schiff reaction stain, Masson's trichrome, Alcian blue, van Gieson, reticulin, Azan, Giemsa, Toluidine blue, isamin blue, sudan black and osmium, acridine orange, Bismarck brown, carmine, Coomassie blue, cresyl violet, DAPI, eosin, ethidium bromide, acid fuchsine, hematoxylin, Hoechst stains, iodine, methyl green, methylene blue, neutral red, Nile blue, Nile red, osmium tetroxide, propidium iodide, rhodamine, safranin, or a combination thereof.

19 . The method of claim 1 , wherein the first image is obtained when the sample is not in contact with the second substrate and wherein the second image is obtained when the biological sample is in contact with the second substrate.

20 . A computer system comprising:

one or more processors;

memory; and

one or more programs, wherein the one or more programs are stored in the memory and configured to be executed by the one or more processors, the one or more programs including instructions for overlaying image data for a biological sample onto spatial analyte data for a plurality of analytes of the biological sample by a method comprising:

obtaining a first image of the biological sample on a first substrate;

receiving a second image of the biological sample on the first substrate overlayed on a second substrate, wherein the second substrate comprises one or more spatial fiducials and a set of capture spots, and wherein the set of capture spots comprises at least 1000 capture spots,

wherein at least one of the first substrate and the second substrate is transparent;

determining a registration for the first image and the second image, wherein the registration registers the first image and the second image to each other, using a first pattern of the biological sample in the first image and a second pattern of the biological sample in the second image; and

using the registration to overlay the first image onto a spatial dataset comprising spatial analyte data for the set of capture spots from the biological sample, wherein a frame of reference of the spatial dataset is known with respect to the second image, based on the one or more spatial fiducials of the second image,

thereby overlaying image data for the biological sample onto spatial analyte data for the plurality of analytes of the biological sample.

21 . A non-transitory computer readable storage medium storing one or more programs, the one or more programs comprising instructions, which when executed by an electronic device with one or more processors and a memory cause the electronic device to perform a method for overlaying image data for a biological sample onto spatial analyte data for a plurality of analytes of the biological sample, comprising:

obtaining a first image of the biological sample on a first substrate;

receiving a second image of the biological sample on the first substrate overlayed on a second substrate, wherein the second substrate comprises one or more spatial fiducials and a set of capture spots, and

wherein the set of capture spots comprises at least 1000 capture spots, wherein at least one of the first substrate and the second substrate is transparent;

determining a registration for the first image and the second image, wherein the registration registers the first image and the second image to each other, using a first pattern of the biological sample in the first image and a second pattern of the biological sample in the second image; and

using the registration to overlay the first image onto a spatial dataset comprising spatial analyte data for the set of capture spots from the biological sample, wherein a frame of reference of the spatial dataset is known with respect to the second image, based on the one or more spatial fiducials of the second image,

thereby overlaying image data for the biological sample onto spatial analyte data for the plurality of analytes of the biological sample.

Assignments (1)
ASSIGNMENT OF ASSIGNOR'S INTEREST Recorded Jul 16, 2024
From: VOROBYOVA, OLGA; CLAYPOOLE, BRYNN; LI, DONGYAO; WEISENFELD, NEIL IRA; LI, PEIGENG; JOSEPH, GUY; SARIKAYA, DIDEM PELIN; SIEGEL, ERIC; KODAVATIKANTI, NAGA SUDHA
To: 10X GENOMICS, INC.
Reel/Frame 068003/0138 →
Continuity (3)
Provisional Application 63350593 · Jun 9, 2022
Provisional Application 63245453 · Sep 17, 2021
Related Publication 20240378734A1 · Nov 14, 2024
References Cited (84)
US 7709198B2 · Luo et al. · 2010 [cited by applicant]
US 8604182B2 · Luo et al. · 2013 [cited by applicant]
US 8951726B2 · Luo et al. · 2015 [cited by applicant]
US 9012022B2 · George et al. · 2015 [cited by applicant]
US 9727810B2 · Fodor et al. · 2017 [cited by applicant]
US 9783841B2 · Nolan et al. · 2017 [cited by applicant]
US 10002316B2 · Fodor et al. · 2018 [cited by applicant]
US 10041949B2 · Bendall et al. · 2018 [cited by applicant]
US 10059990B2 · Boyden et al. · 2018 [cited by applicant]
US 10640816B2 · Beechem et al. · 2020 [cited by applicant]
US 10913397B2 · Pribisic · 2021 [cited by applicant]
US 11501440B2 · Weisenfeld et al. · 2022 [cited by applicant]
US 11514575B2 · Mellen et al. · 2022 [cited by applicant]
US 20110059865A1 · Smith et al. · 2011 [cited by applicant]
US 20160019695A1 · Chukka · 2016 [cited by examiner]
US 20180253871A1 · Bredno · 2018 [cited by examiner]
US 20200277663A1 · Ramachandran Iyer · 2020 [cited by applicant]
US 20210062272A1 · Williams et al. · 2021 [cited by applicant]
US 20210150707A1 · Weisenfeld et al. · 2021 [cited by applicant]
US 20210155982A1 · Yin · 2021 [cited by examiner]
US 20210183034A1 · Chang · 2021 [cited by examiner]
US 20220065787A1 · Kim · 2022 [cited by examiner]
WO WO2011094669 · 2011 [cited by applicant]
WO WO2011127099 · 2011 [cited by applicant]
WO WO2012140224 · 2012 [cited by applicant]
WO WO2014060483 · 2014 [cited by applicant]
WO WO2014163886A1 · 2014 [cited by applicant]
WO WO2014210225 · 2014 [cited by applicant]
WO WO2014210233 · 2014 [cited by applicant]
WO WO2015161173 · 2015 [cited by applicant]
WO WO2016007839A1 · 2016 [cited by applicant]
WO WO2016057552 · 2016 [cited by applicant]
WO WO2016162309 · 2016 [cited by applicant]
WO WO2016166128 · 2016 [cited by applicant]
WO WO2017027367A1 · 2017 [cited by applicant]
WO WO2017027368 · 2017 [cited by applicant]
WO WO2017144338A1 · 2017 [cited by applicant]
WO WO2017147483 · 2017 [cited by applicant]
WO WO2017222453A1 · 2017 [cited by applicant]
WO WO2018022809A1 · 2018 [cited by applicant]
WO WO2018045181A1 · 2018 [cited by applicant]
WO WO2018045186A1 · 2018 [cited by applicant]
WO WO2018057999A1 · 2018 [cited by applicant]
WO WO2018091676A1 · 2018 [cited by applicant]
WO WO2018107054A1 · 2018 [cited by applicant]
WO WO2018136856A1 · 2018 [cited by applicant]
WO WO2019068880A1 · 2019 [cited by applicant]
WO WO2019075091A1 · 2019 [cited by applicant]
WO WO2020123320A2 · 2020 [cited by applicant]
WO WO2020176788A1 · 2020 [cited by applicant]
WO WO2021092433A2 · 2021 [cited by applicant]
WO WO2022061150A2 · 2022 [cited by applicant]
10X, “Inside Visium spatial capture technology,” Oct. 14, 2019, product brochure, retrieved from the Internet: https://pages.10xgenomics.com/rs/446-PBO-704/images/10x_BR060_Inside_Visium_Spatial_Technology.pdf. [cited by applicant]
10X, “Spatially Resolved Transcriptomics: The Landscape of Tools for Spatial Gene Expression,” Sep. 9, 2020, webinar, retrieved from the Internet: https://www.10xgenomics.com/videos/80lda5wtby?autoplay=true. [cited by applicant]
10X, “Visium Spatial Gene Expression Solution,” May 29, 2020, product brochure, retrieved from the Internet: https://pages.10xgenomics.com/rs/446-PBO-704/images/10x_LIT054_Spatial-Gene-Expression_Flyer_Letter_Digital.pd… [cited by applicant]
Achim et al., 2015, “High-throughput spatial mapping of single-cell RNA-seq data to tissue of origin,” Nature Biotechnology 33: 503-509, doi:10.1038/nbt.3209. [cited by applicant]
Bolognesi et al., 2017, J. Histochem. Cytochem. 65(8): 431-444. [cited by applicant]
Boykov, et al., 2001, “Interactive graph cuts for optimal boundary and region segmentation of objects in N-D images,” Proc. IEEE Int. Conf. on Computer Vision, CD-ROM. [cited by applicant]
Canny, 1986, “A Computational Approach to Edge Detection,” IEEE Trans Pattern Anal Mach Intell. 8(6):679-98. [cited by applicant]
Carter et al., 2007, Applied Optics 46:421-427),. [cited by applicant]
Chen et al., Science 348(6233): aaa6090, 2015. [cited by applicant]
Chetverikov et al., 2002, “The Trimmed Iterative Closest Point Algorithm,” Object recognition supported by user interaction for service robots, Quebec City, Quebec, Canada, ISSN: 1051-4651. [cited by applicant]
Chetverikov et al., 2005, “Robust Euclidean alignment of 3D point sets; the trimmed iterative closest point algorithm,” Image and Vision Computing 23(3), pp. 299-309. [cited by applicant]
Chuang, et al., 2001, “A Bayesian approach to digital matting,” Proc. IEEE Conf. [cited by applicant]
Chui and Rangarajanb, 2003, “A new point matching algorithm for non-rigid registration,” Computer Vision and Image Understanding 89(2-3), pp. 114-141. [cited by applicant]
Gao et al., 2017, BMC Biology 15:50, doi:10.1186/s12915-017-0393-3. [cited by applicant]
Glass et al., 2009, J. Histochem. Cytochem. 57:899-905. [cited by applicant]
Greig et al., 1989, “Exact MAP estimation for binary images,” J. Roy. Stat. Soc. B. 51, 271-279. [cited by applicant]
Gupta et al., [cited by applicant]
Lee et al., [cited by applicant]
Lin et al., 2015, Nat Commun. 6:8390. [cited by applicant]
McCormick et al., “ITK: enabling reproducible research and open science,” Front Neuroinform 2014;8:13; doi:10.3389/fninf.2014.00013. [cited by applicant]
Myronenko, Andriy, et al. “Non-rigid point set registration: Coherent Point Drift”, Department of Computer Science and Electrical Engineering, OGI School of Science and Engineering, Beaverton, OR, pp. 1-8 (2006). [cited by applicant]
Myronenko, Andriy, et al. “Point Set Registration: Coherent Point Drift”, Department of Science and Engineering, School of Medicine, Oregon Health and Science University, Porland, OR, May 15, 2009, pp. 1-14. [cited by applicant]
Phogat, et al., 2012, “Different Image Registration Methods—An Overview,” Int J Sci Eng Res 5(12): 44-49. [cited by applicant]
Pirici et al., 2009, J. Histochem. Cytochem. 57:567-75. [cited by applicant]
Rodriques et al., [cited by applicant]
Rother et al., 2004, “‘GrabCut’—Interactive Foreground Extraction using Iterated Graph Cuts,” ACM Transactions on Graphics. 23(3):309-314, doi:10.1145/1186562.1015720. [cited by applicant]
Satija et al., 2015, “Spatial reconstruction of single-cell gene expression data,” Nature Biotechnology. 33, 495-502, doi:10.1038.nbt.3192. [cited by applicant]
Trejo et al., [cited by applicant]
Uchida, 2013, “Image processing and recognition for biological images,” Develop. Growth Differ., The Japanese Society of Developmental Biologists 55, 523-549. [cited by applicant]
Woo, J., “Multimodal Registration via Mutual Information Incorporating Geometric and Spatial Context,” IEEE Trans Image Process. Feb. 2015 ; 24(2): 757-769. [cited by applicant]
Yang, 2011, “The thin plate spline robust point matching (TPS-RPM) algorithm: A revisit,” Pattern Recognition Letters 32(7), pp. 910-918. [cited by applicant]
Yoo, et al., “Engineering and Algorithm Design for an Image Processing API: A Technical Report on ITK—The Insight Toolkit,” In Proc of Medicine Meets Virtual Reality, J Westwood, ed, IOS Press Amsterdam pp. 586-592 (200… [cited by applicant]