DE102015110316B4 - Vorrichtungen, Zytometer, Verfahren und Computerprogramm zum Bereitstellen von Information über zumindest eine Sequenz - Google Patents

Vorrichtungen, Zytometer, Verfahren und Computerprogramm zum Bereitstellen von Information über zumindest eine Sequenz Download PDF

Info

Publication number
DE102015110316B4
DE102015110316B4 DE102015110316.8A DE102015110316A DE102015110316B4 DE 102015110316 B4 DE102015110316 B4 DE 102015110316B4 DE 102015110316 A DE102015110316 A DE 102015110316A DE 102015110316 B4 DE102015110316 B4 DE 102015110316B4
Authority
DE
Germany
Prior art keywords
sequence
maximum
signal states
designed
signal
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
DE102015110316.8A
Other languages
English (en)
Other versions
DE102015110316A1 (de
Inventor
Michael Bassler
Stephan Quint
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Fraunhofer Gesellschaft zur Forderung der Angewandten Forschung eV
Original Assignee
Fraunhofer Gesellschaft zur Forderung der Angewandten Forschung eV
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Fraunhofer Gesellschaft zur Forderung der Angewandten Forschung eV filed Critical Fraunhofer Gesellschaft zur Forderung der Angewandten Forschung eV
Priority to DE102015110316.8A priority Critical patent/DE102015110316B4/de
Priority to CN201680036244.2A priority patent/CN107771278B/zh
Priority to PCT/EP2016/064492 priority patent/WO2016207262A1/de
Priority to US15/738,162 priority patent/US10551295B2/en
Publication of DE102015110316A1 publication Critical patent/DE102015110316A1/de
Application granted granted Critical
Publication of DE102015110316B4 publication Critical patent/DE102015110316B4/de
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N15/00Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
    • G01N15/10Investigating individual particles
    • G01N15/14Optical investigation techniques, e.g. flow cytometry
    • G01N15/1456Optical investigation techniques, e.g. flow cytometry without spatial resolution of the texture or inner structure of the particle, e.g. processing of pulse signals
    • G01N15/1459Optical investigation techniques, e.g. flow cytometry without spatial resolution of the texture or inner structure of the particle, e.g. processing of pulse signals the analysis being performed on a sample stream
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N15/00Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
    • G01N15/10Investigating individual particles
    • G01N15/14Optical investigation techniques, e.g. flow cytometry
    • G01N15/1434Optical arrangements
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/62Systems in which the material investigated is excited whereby it emits light or causes a change in wavelength of the incident light
    • G01N21/63Systems in which the material investigated is excited whereby it emits light or causes a change in wavelength of the incident light optically excited
    • G01N21/64Fluorescence; Phosphorescence
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • G06F17/15Correlation function computation including computation of convolution operations
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N15/00Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
    • G01N15/10Investigating individual particles
    • G01N15/14Optical investigation techniques, e.g. flow cytometry
    • G01N15/1434Optical arrangements
    • G01N2015/1447Spatial selection

Landscapes

  • Physics & Mathematics (AREA)
  • Chemical & Material Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Dispersion Chemistry (AREA)
  • Data Mining & Analysis (AREA)
  • Computational Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Computing Systems (AREA)
  • Algebra (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Apparatus Associated With Microorganisms And Enzymes (AREA)
  • Investigating, Analyzing Materials By Fluorescence Or Luminescence (AREA)
  • Investigating Or Analysing Biological Materials (AREA)

Abstract

Vorrichtung (20) zum Bereitstellen von Information über zumindest eine Sequenz, wobei die zumindest eine Sequenz zeitlich aufeinanderfolgende Signalzustände beschreibt, die Vorrichtung (20) umfassendeine Schnittstelle (22), die ausgebildet ist, um Information über eine Anzahl der Signalzustände zu erhalten; undein Berechnungsmodul (24), welches ausgebildet ist, um:eine Mehrzahl von möglichen Sequenzen zu generieren, basierend auf der Information über die Anzahl der Signalzustände,für zumindest eine Teilmenge der möglichen Sequenzen Korrelationsfunktionen zwischen einer Sequenz und zumindest einer zeitlichen Skalierung der Sequenz zu berechnen, wobei eine Korrelationsfunktion ein Hauptmaximum und ein oder mehrere Nebenmaxima umfasst,die zumindest eine Sequenz basierend auf den Korrelationsfunktionen zu bestimmen, wobei die Folge der Signalzustände innerhalb der zumindest einen Sequenz derart gewählt ist, dass ein Nebenmaximum in einer Korrelationsfunktion der Sequenz mit der zumindest einen zeitlichen Skalierung der Sequenz gegenüber einem durch unterschiedliche Anordnungen der Signalzustände in der Sequenz maximal in einer Korrelationsfunktion erzielbaren Nebenmaximum reduziert ist, unddie Information über die zumindest eine Sequenz basierend auf der zumindest einen Sequenz zu bestimmen und über die Schnittstelle (22) bereitzustellen.

Description

  • Technisches Gebiet
  • Ausführungsbeispiele befassen sich mit einer Vorrichtung, einem Verfahren und einem Computerprogramm zum Bereitstellen von Information über zumindest eine Sequenz, mit einer Vorrichtung, einem Verfahren und einem Computerprogramm für ein Zytometer zur Bereitstellung von Information über ein oder mehrere Zellen in einem Medium in einem Kanal, und mit einem Zytometer.
  • Hintergrund
  • In der medizinischen und biologischen Forschung und in der Analyse und Erkennung von Krankheiten basieren viele Erkenntnisse auf einem Analysieren von Zellen. Für die Analyse großer Mengen von Zellen kommen oft Durchflusszytometer zum Einsatz. In einem Durchflusszytometer werden Zellen in einer Lösung mit einer hohen Geschwindigkeit durch einen Kanal geleitet. Diese Zellen senden ein optisches Signal aus, beispielsweise angeregt durch eine Lichtquelle, etwa einen Laser. Dieses optische Signal wird von dem Durchflusszytometer erfasst und ermöglicht eine Bestimmung der Eigenschaften der Zellen in der Lösung, etwa die Anzahl von bestimmten Zelltypen, oder ihre Größe und andere Eigenschaften.
  • Eine Weiterentwicklung der Durchflusszytometrie basiert auf dem Prinzip der Räumlich Modulierten Frequenz (RMF). In konventionellen Flusszytometern wird der Kanal, in dem die Lösung mit den Zellen am Sensor vorbeigeführt wird, so verengt, dass die Zellen den Sensor einzeln passieren. Dabei sind hochpräzise Optiken für den Detektor und den Laser für die Anregung und Erfassung der einzelnen Zellen nötig, die einen komplexen Aufbau und großen Platzbedarf des Flusszytometers bedingen. Durchflusszytometer, die auf dem Prinzip der Räumlich Modulierten Frequenz basieren, können jedoch den Durchfluss mehrerer Zellen gleichzeitig unterstützen, beispielsweise bei einer geringeren Komplexität des optischen Aufbaus. Die Zellen in Lösung im Kanal werden gleichzeitig angeregt, beispielsweise von einem Laser, und das emittierte Licht wird von einem räumlichen Filter, etwa einer Filtermaske, beeinflusst. Der räumliche Filter sorgt dafür, dass der Sensor für jede angeregte Zelle ein Signal basierend auf dem Filter erfasst. Durch eine Analyse des Signals, welches von den Signalen anderer Zellen überlagert sein kann, kann das Zytometer eine Anzahl, Geschwindigkeit und andere Eigenschaften der Zellen in der Lösung bestimmen. Die Qualität der Erkennung der Zellen und Zelleigenschaften basiert dabei auf dem Filter, welcher sich beispielsweise durch eine Sequenz zeitlich aufeinanderfolgende Signalzustände beschreiben lässt.
  • Weitere Informationen können beispielsweise in folgenden Dokumenten gefunden werden:
    • P. Kiesel, M. Bassler, M. Beck and N. Johnson, Appl. Phys. Lett., 2009, 94, p. 41107;
    • N. Levanon and E. Mozeson, Radar Signals, Wiley, New York, 2004;
    • R. Turyn and S. Storer, On binary sequences, Proceedings of the American Mathematical Society 12, 394 (1961);
    • S. Mertens, Exhaustive search for low-autocorrelation binary sequences, Journal of Physics A: Mathematical and General 29, L473 (1996);
    • G. E. Coxson and J. C. Russo, Efficient exhaustive search for optimal-peak-sidelobe binary codes, IEEE Transactions on Aerospace and Electronic Systems 41, 302 (2005);
    • P. Borewein, R. Ferguson, and J. Knauer, The merit factor problem, 2000; C. J. Nunn and G. E. Coxson, Best-known autocorrelation peak sidelobe levels for binary codes oflength 71 to 105, IEEE Transactions on Aerospace and Electronic Systems 44, 392 (2008);
    • S. K. Shanmugam, C. Mongredien, J. Nielsen, and G. Lachapelle, Design of short synchronization codes for use in future GNSS system, International Journal of Navigation and Observation 2008 (2008), Article ID 246703;
    • C. Sommer: Die Größenabhängigkeit der Gleichgewichtsgeschwindigkeit von Partikeln beim Transport in Mikrokanälen, Dissertation der Technischen Universität Darmstadt, D17, Juli 2014;
    • S. Quint: Durchflusszytometrie nach dem Prinzip der räumlich modulierten Fluoreszenz. Dissertation der Universität Darmstadt, 2014, Abstract;
    • R. Foster: Manchester encoding: opposing definitions resolved, Engineering Science and Education Journal, Vol. 9, 2000, No. 6, S. 278-280. ISSN 0963-7346;
    • S. W. Kettlitz et al.: Particle detection form spatially modulated fluorescence Signals, Proceedings of SPIE, Vol. 9129, 2014, S. 1-7. - doi: 10.1117/12.2052169;
    • S. W. Kettlitz et al.: Sensitivity Improvement in Fluorescence-Based Particle Detection, Cytometry A, Vol. 85A, 2014, S. 746-755, - doi: 10.1002/cyto.a.22499;
    • J. Lindner: Binary sequences up to length 40 with best possible autocorrelation function, Electronics Letters, Vol. 11, 1975, No. 21, S. 507;
    • C. Sommer et al.: The equilibrium velocity of spherical particles in rectangular microfluidic channels for size measurement, Lab on a Chip, Vol. 14, 2014, No. 13, S. 2319-2326;
    • M. N. Cohen et al: Minimum peak sidelobe pulse compression codes, IEEE International Conference on Radar, 7-10 May 1990, S. 633-638; und
    • H. D. Lüke et al: Binary and quadriphase sequences with optimal autocorrelation properties: a survey, IEEE Transactions on Information Theory, Vol. 49, 2003, No. 12, S. 3271-3282.
  • Es besteht mithin der Bedarf, Sequenzen zu finden, mit denen sich die Erkennung der Zellen verbessern lässt. Diesem Bedarf wird durch eine Vorrichtung, ein Verfahren und eine Computerprogramm zum Bereitstellen von Information über zumindest eine Sequenz Rechnung getragen.
  • Zusammenfassung
  • Ausführungsbeispiele beruhen auf der Erkenntnis, dass Sequenzen, die in einem Zytometer genutzt werden, das auf RMF basiert, vorteilhafte Eigenschaften aufweisen, wenn sie nicht nur in einer Autokorrelationsfunktion bzw. Korrelationsfunktion mit einer amplitudenskalierten Variante reduzierte Nebenmaxima aufweisen, sondern auch Nebenmaxima in Kreuzkorrelationsfunktionen der Sequenz mit zeitlichen Skalierungen der Sequenz reduziert sind. Die Reduktion der Nebenmaxima der Kreuzkorrelationsfunktion kann dabei eine verbesserte Erkennung von verschieden schnellen Zellen ermöglichen. Ausführungsbeispiele können ferner auf der Erkenntnis basieren, dass Sequenzen, die in einem Zytometer genutzt werden, das auf RMF basiert, vorteilhafte Eigenschaften aufweisen, wenn sie unipolar sind und/oder wenn die Korrelationsfunktionen basierend auf balancierten Filtern bestimmt werden.
  • Ausführungsbeispiele schaffen eine Vorrichtung zum Bereitstellen von Information über zumindest eine Sequenz. Die zumindest eine Sequenz beschreibt zeitlich aufeinanderfolgende Signalzustände. Die Vorrichtung umfasst eine Schnittstelle, die ausgebildet ist, um Information über eine Anzahl der Signalzustände zu erhalten. Die Vorrichtung umfasst ferner ein Berechnungsmodul, welches ausgebildet ist, um eine Mehrzahl von möglichen Sequenzen zu generieren, basierend auf der Information über die Anzahl der Signalzustände. Das Berechnungsmodul ist ferner ausgebildet, um für zumindest eine Teilmenge der möglichen Sequenzen Korrelationsfunktionen zwischen einer Sequenz und zumindest einer zeitlichen Skalierung der Sequenz zu berechnen. Eine Korrelationsfunktion umfasst ein Hauptmaximum und ein oder mehrere Nebenmaxima. Das Berechnungsmodul ist ferner ausgebildet, die zumindest eine Sequenz basierend auf den Korrelationsfunktionen zu bestimmen. Die Folge der Signalzustände innerhalb der zumindest einen Sequenz ist derart gewählt, dass ein Nebenmaximum in einer Korrelationsfunktion der Sequenz mit der zumindest einen zeitlichen Skalierung der Sequenz gegenüber einem durch unterschiedliche Anordnungen der Signalzustände in der Sequenz maximal in einer Korrelationsfunktion erzielbaren Nebenmaximum reduziert ist. Das Berechnungsmodul ist ferner ausgebildet, um die Information über die zumindest eine Sequenz basierend auf der zumindest einen Sequenz zu bestimmen und über die Schnittstelle bereitzustellen. Die zumindest eine Sequenz kann in Ausführungsbeispielen genutzt werden, um die Erkennung von zeitlich skalierten Signalen mittels einer Filterbank zu verbessern, und um Nebenmaxima einer Korrelationsfunktion einer empfangenen Sequenz und zeitlich skalierten Sequenzen zu reduzieren und damit Fehlerkennungen zu reduzieren.
  • In einigen Ausführungsbeispielen kann das Berechnungsmodul ausgebildet sein, um die zumindest eine Sequenz basierend auf den Korrelationsfunktionen zu bestimmen. Die Folge der Signalzustände innerhalb der zumindest einen Sequenz könnte derart gewählt sein, dass eine Summe, die auf den Nebenmaxima einer Korrelationsfunktion der Sequenz mit der zumindest einen zeitlichen Skalierung der Sequenz basiert, gegenüber einer durch unterschiedliche Anordnungen der Signalzustände in der Sequenz maximal erzielbaren Summe, die auf aus der unterschiedlichen Anordnung resultierenden Nebenmaxima basiert, reduziert ist. Die Folge der Signalzustände innerhalb der zumindest einen Sequenz könnte derart gewählt sein, dass ein größtes Nebenmaximum in der Korrelationsfunktion der Sequenz mit der zumindest einen zeitlichen Skalierung der Sequenz gegenüber einem durch unterschiedliche Anordnungen der Signalzustände in der Sequenz maximal erzielbaren größten Nebenmaximum reduziert ist. Eine Reduzierung einer Summe, die auf den Nebenmaxima einer Korrelationsfunktion der Sequenz mit der zumindest einen zeitlichen Skalierung der Sequenz basiert, kann beispielsweise in der Erkennung mittels einer Filterbank ein Grundrauschen in der Korrelationsanalyse reduzieren und/oder bei stark verrauschten Signalen eine geringere Fehlerkennungsrate ergeben, da eine Gesamtenergie in den Nebenmaxima geringer ist als bei Sequenzen, die auf eine Reduzierung des höchsten Nebenmaximums optimiert sind. Eine Reduzierung eines höchsten Nebenmaximums kann beispielsweise in der Erkennung mittels einer Filterbank Fehlerkennungen reduzieren, und/oder könnte eine Nutzung niedriger Schwellenwerte, d.h. höhere Sensitivität, bei gleicher Fehlerrate, ermöglichen.
  • In zumindest manchen Ausführungsbeispielen könnte das Berechnungsmodul ausgebildet sein, um die zumindest eine Sequenz so zu bestimmen, dass die Summe, die auf den Nebenmaxima basiert, unter einem ersten Schwellenwert liegt. Das Berechnungsmodul könnte ausgebildet sein, um die zumindest eine Sequenz so zu bestimmen, dass das größte Nebenmaximum unter einem zweiten Schwellenwert liegt. In Ausführungsbeispielen könnte der erste und/oder der zweite Schwellenwert genutzt werden, um eine Bestimmung der Sequenz zu beschleunigen und/oder um bestimmte Sequenzen zu evaluieren.
  • In manchen Ausführungsbeispielen könnte der erste Schwellenwert einer zweitniedrigsten Summe, die auf den Nebenmaxima der Korrelationsfunktionen basiert, entsprechen. Der zweite Schwellenwert könnte einem größten Nebenmaximum einer Sequenz mit einem zweitniedrigsten größten Nebenmaximum der Korrelationsfunktionen entsprechen. Durch die Wahl der zweitniedrigsten Summe und des zweitniedrigsten größten Nebenmaximums als ersten und zweiten Schwellenwert kann das Bestimmen einer Sequenz, deren Korrelationsfunktion eine niedrigste Summe oder ein niedrigstes größtes Nebenmaximum ergibt, ermöglicht werden, beispielsweise auch rekursiv.
  • In zumindest einigen Ausführungsbeispielen könnte das Berechnungsmodul ausgebildet sein, um bereinigte Korrelationsfunktionen basierend auf den Korrelationsfunktionen zu bestimmen. Das Berechnungsmodul könnte ausgebildet sein, um einen Beitrag von Hauptmaxima in den bereinigten Korrelationsfunktionen zu reduzieren. Das Berechnungsmodul könnte ausgebildet sein, um die Information über die zumindest eine Sequenz basierend auf den bereinigten Korrelationsfunktionen zu bestimmen. Die Reduzierung des Beitrags der Hauptmaxima in den bereinigten Korrelationsfunktionen könnte eine Analyse der ein oder mehreren Nebenmaxima erleichtern und könnte einer Erkennung und Reduzierung der Hauptmaxima in einem Zytometer entsprechen.
  • In manchen Ausführungsbeispielen könnte das Berechnungsmodul ausgebildet sein, um unter den Hauptmaxima der Korrelationsfunktionen ein größtes Hauptmaximum zu bestimmen. Das Berechnungsmodul könnte ausgebildet sein, um über die zeitliche Ausdehnung des größten Hauptmaximums den Beitrag der Hauptmaxima in den bereinigten Korrelationsfunktionen zu reduzieren. Eine Reduzierung des Beitrags der Hauptmaxima über die zeitliche Ausdehnung des Hauptmaximums könnte einfacher numerisch umzusetzen sein.
  • In einigen Ausführungsbeispielen könnte das Berechnungsmodul ausgebildet sein, um basierend auf einer zeitlichen Position des größten Hauptmaximums und basierend auf der Korrelationsfunktion der zeitlich skalierten Sequenz die das größte Hauptmaximum umfasst, die Beiträge der Hauptmaxima in den bereinigten Korrelationsfunktionen der zeitlich skalierten Sequenz und benachbarten zeitlichen Skalierungen der Sequenz zu bestimmen und zu reduzieren. Durch die Reduzierung der Beiträge der Hauptmaxima in den bereinigten Korrelationsfunktionen der zeitlich skalierten Sequenz und benachbarten zeitlichen Skalierungen der Sequenz könnte beispielsweise eine genauere Isolierung der Hauptmaxima erreicht werden und Doppelereignisse könnten unterstützt werden.
  • In zumindest einigen Ausführungsbeispielen könnte eine Sequenz eine mittlere Strukturgröße haben. Eine Strukturgröße könnte auf einer Anzahl von gleichen aufeinanderfolgenden Signalzuständen basieren. Das Berechnungsmodul könnte ausgebildet sein, um die Mehrzahl von möglichen Sequenzen basierend auf einem Zielbereich für die mittlere Strukturgröße zu generieren. Das Berechnungsmodul könnte ausgebildet sein, die Teilmenge der möglichen Sequenzen basierend auf der mittleren Strukturgröße und/oder dem Zielbereich für die mittlere Strukturgröße zu bestimmen. Eine Bestimmung der Teilmenge der möglichen Sequenzen basierend auf der mittleren Strukturgröße könnte beispielsweise ein Bestimmen der zumindest einen Sequenz beschleunigen und einen Berechnungsaufwand für das Bestimmen der zumindest einen Sequenz reduzieren.
  • In zumindest manchen Ausführungsbeispielen könnte die mittlere Strukturgröße einem Durchschnitt der Strukturgrößen einer Sequenz entsprechen. Der Zielbereich für die mittlere Strukturgröße könnte beispielsweise zwischen 1,3 und 1,8 liegen. Eine Bestimmung der Teilmenge der möglichen Sequenzen basierend auf der mittleren Strukturgröße und dem Zielbereich könnte beispielsweise ein Bestimmen der zumindest einen Sequenz beschleunigen und einen Berechnungsaufwand für das Bestimmen der zumindest einen Sequenz reduzieren.
  • In zumindest manchen Ausführungsbeispielen kann das Berechnungsmodul ausgebildet sein, um die Mehrzahl von möglichen Sequenzen als unipolare Sequenzen zu bestimmen. Durch eine Nutzung von unipolaren Sequenzen könnte ein verbessertes Hauptmaximum zu Nebenmaxima-Verhältnis erreicht werden.
  • In einigen Ausführungsbeispielen könnte eine Sequenz einer Binärsequenz entsprechen. Die Binärsequenz könnte zeitlich äquidistante aufeinanderfolgende Signalzustände beschreiben. Die Nutzung einer Binärsequenz könnte beispielsweise eine einfache Notation der Sequenz ermöglichen.
  • Ausführungsbeispiele stellen ferner eine Vorrichtung für ein Zytometer bereit. Die Vorrichtung ist ausgebildet zur Bereitstellung von Information über ein oder mehrere Zellen in einem Medium in einem Kanal. Die Vorrichtung umfasst ein Sensormodul, das ausgebildet ist, um ein gefiltertes Signal zu erfassen. Das gefilterte Signal wird von den durch den Kanal strömenden Zellen beeinflusst. Die Vorrichtung umfasst ferner ein Filter, das ausgebildet ist, um eine vorgegebene Sequenz räumlich abzubilden, und ausgebildet ist, um basierend auf der Sequenz und einem von den durch den Kanal strömenden Zellen beeinflussten Signal das gefilterte Signal bereitzustellen. Die Vorrichtung umfasst ferner ein Kontrollmodul das ausgebildet ist, um die Information über die ein oder mehreren Zellen basierend auf einer Korrelationsanalyse des erfassten gefilterten Signals und zumindest einer zeitlichen Skalierung der Sequenz bereitzustellen. Die Sequenz beschreibt zeitlich aufeinanderfolgende Signalzustände. Die Folge der Signalzustände innerhalb der Sequenz ist derart gewählt, dass ein Nebenmaximum in einer Korrelationsfunktion der Sequenz mit der zumindest einen zeitlichen Skalierung der Sequenz gegenüber einem durch unterschiedliche Anordnungen der Signalzustände in der Sequenz maximal erzielbaren Nebenmaximum in der Korrelationsfunktion reduziert ist. Die Sequenz kann in Ausführungsbeispielen genutzt werden, um die Erkennung von zeitlich skalierten Signalen mittels einer Filterbank, welche beispielsweise im Kontrollmodul umfasst sein könnte, zu verbessern, und um Nebenmaxima einer Korrelationsfunktion einer empfangenen Sequenz und zeitlich skalierten Sequenzen zu reduzieren und damit Fehlerkennungen zu reduzieren.
  • In einigen Ausführungsbeispielen könnte eine Folge der Signalzustände innerhalb der zumindest einen Sequenz derart gewählt sein, dass eine Summe, die auf Nebenmaxima der Korrelationsfunktion der Sequenz mit der zumindest einen zeitlichen Skalierung der Sequenz basiert, gegenüber einer durch unterschiedliche Anordnungen der Signalzustände in der Sequenz maximal erzielbaren Summe, die auf aus der unterschiedlichen Anordnung resultierenden Nebenmaxima basiert, reduziert ist. Die Folge der Signalzustände innerhalb der zumindest einen Sequenz könnte derart gewählt sein, dass ein größtes Nebenmaximum der Korrelationsfunktion der Sequenz mit der zumindest einen zeitlichen Skalierung der Sequenz gegenüber einem durch unterschiedliche Anordnungen der Signalzustände in der Sequenz maximal erzielbaren größten Nebenmaximum reduziert ist. Eine Reduzierung einer Summe, die auf den Nebenmaxima einer Korrelationsfunktion der Sequenz mit der zumindest einen zeitlichen Skalierung der Sequenz basiert, könnte beispielsweise in der Erkennung mittels einer Filterbank ein Grundrauschen in der Korrelationsanalyse reduzieren. Eine Reduzierung eines höchsten Nebenmaximums könnte beispielsweise in der Erkennung mittels einer Filterbank Fehlerkennungen reduzieren.
  • In manchen Ausführungsbeispielen könnte die Summe, die auf den Nebenmaxima basiert, unter einem ersten Schwellenwert liegen. Das größte Nebenmaximum könnte unter einem zweiten Schwellenwert liegen. In Ausführungsbeispielen kann der erste und/oder der zweite Schwellenwert genutzt werden, um eine Bestimmung der Sequenz zu beschleunigen und/oder um bestimmte Sequenzen zu evaluieren.
  • In zumindest einigen Ausführungsbeispielen könnte der erste Schwellenwert einer zweitniedrigsten Summe, die auf Nebenmaxima von Korrelationsfunktionen möglicher Sequenzen basiert, entsprechen. Der zweite Schwellenwert könnte einem größten Nebenmaximum einer Sequenz mit einem zweitniedrigsten größten Nebenmaximum der Korrelationsfunktionen entsprechen. Durch die Wahl der zweitniedrigsten Summe und des zweitniedrigsten größten Nebenmaximums als ersten und zweiten Schwellenwert könnte das Bestimmen einer Sequenz, deren Korrelationsfunktion eine niedrigste Summe oder ein niedrigstes größtes Nebenmaximum ergibt, ermöglicht werden, beispielsweise auch rekursiv.
  • In manchen Ausführungsbeispielen könnte die Sequenz einer unipolaren Sequenz entsprechen. Durch eine Nutzung einer unipolaren Sequenz könnte ein verbessertes Hauptmaximum zu Nebenmaxima-Verhältnis erreicht werden.
  • Ausführungsbeispiele schaffen ferner ein Zytometer, welches die Vorrichtung sowie einen Kanal umfasst. Die Vorrichtung ist ausgebildet, um Information über ein oder mehrere Zellen in einem Medium in dem Kanal dem Zytometer bereitzustellen. Die Vorrichtung, welche die Sequenz nutzt, könnte in Ausführungsbeispielen genutzt werden, um die Erkennung von zeitlich skalierten Signalen, etwa optischen Signalen der Zellen im Kanal, mittels einer Filterbank, welche beispielsweise im Kontrollmodul der Vorrichtung umfasst sein könnte, zu verbessern, und um Nebenmaxima einer Korrelationsfunktion einer empfangenen Sequenz und zeitlich skalierten Sequenzen zu reduzieren und damit Fehlerkennungen zu reduzieren.
  • Ausführungsbeispiele schaffen ferner ein Verfahren zum Bereitstellen von Information über zumindest eine Sequenz. Die zumindest eine Sequenz beschreibt zeitlich aufeinanderfolgende Signalzustände. Das Verfahren umfasst Erhalten von Information über eine Anzahl der Signalzustände. Das Verfahren umfasst ferner Generieren einer Mehrzahl von möglichen Sequenzen, basierend auf der Information über die Anzahl der Signalzustände. Das Verfahren umfasst ferner Berechnen von Korrelationsfunktionen zwischen einer Sequenz und zumindest einer zeitlichen Skalierung der Sequenz für zumindest eine Teilmenge der möglichen Sequenzen. Eine Korrelationsfunktion umfasst ein Hauptmaximum und ein oder mehrere Nebenmaxima. Das Verfahren umfasst ferner Bestimmen der zumindest einen Sequenz basierend auf den Korrelationsfunktionen. Die Folge der Signalzustände innerhalb der zumindest einen Sequenz ist derart gewählt, dass ein Nebenmaximum in einer Korrelationsfunktion der Sequenz mit der zumindest einen zeitlichen Skalierung der Sequenz gegenüber einem durch unterschiedliche Anordnungen der Signalzustände in der Sequenz maximal in einer Korrelationsfunktion erzielbaren Nebenmaximum reduziert ist. Das Verfahren umfasst ferner Bestimmen der Information über die zumindest eine Sequenz basierend auf der zumindest einen Sequenz, und Bereitstellen der Information über die zumindest eine Sequenz.
  • Ausführungsbeispiele schaffen ferner ein Verfahren zur Bereitstellung von Information über ein oder mehrere Zellen in einem Medium in einem Kanal. Die Vorrichtung umfasst Erfassen eines gefilterten Signals. Das gefilterte Signal wird von den durch den Kanal strömenden Zellen beeinflusst. Das Verfahren umfasst ferner Bereitstellen des gefilterten Signals basierend auf einer vorgegebenen Sequenz, einem von den durch den Kanal strömenden Zellen beeinflussten Signal und einer räumlichen Abbildung der Sequenz. Das Verfahren umfasst ferner Bereitstellen der Information über die ein oder mehreren Zellen basierend auf einer Korrelationsanalyse des erfassten gefilterten Signals und zumindest einer zeitlichen Skalierung der Sequenz. Die Sequenz beschreibt zeitlich aufeinanderfolgende Signalzustände. Die Folge der Signalzustände innerhalb der Sequenz ist derart gewählt, dass ein Nebenmaximum in einer Korrelationsfunktion der Sequenz mit der zumindest einen zeitlichen Skalierung der Sequenz gegenüber einem durch unterschiedliche Anordnungen der Signalzustände in der Sequenz maximal erzielbaren Nebenmaximum in der Korrelationsfunktion reduziert ist.
  • Ein weiteres Ausführungsbeispiel ist ein Computerprogramm zur Durchführung zumindest eines der oben beschriebenen Verfahren, wenn das Computerprogramm auf einem Computer, einem Prozessor oder einer programmierbaren Hardwarekomponente abläuft. Ein weiteres Ausführungsbeispiele ist auch ein digitales Speichermedium, das maschinen- oder computerlesbar ist, und das elektronisch lesbare Steuersignale aufweist, die mit einer programmierbaren Hardwarekomponente so zusammenwirken können, dass eines der oben beschriebenen Verfahren ausgeführt wird.
  • Figurenliste
  • Weitere vorteilhafte Ausgestaltungen werden nachfolgend anhand der in den Zeichnungen dargestellten Ausführungsbeispiele, auf welche Ausführungsbeispiele generell jedoch nicht insgesamt beschränkt sind, näher beschrieben. Es zeigen:
    • 1 zeigt ein Blockdiagramm eines Ausführungsbeispiels einer Vorrichtung zum Bereitstellen von Information über zumindest eine Sequenz;
    • 2 zeigt ein Diagramm des Prinzips der Räumlich Modulierten Fluoreszenz;
    • 3 zeigt ein Diagramm einer zeitlichen und räumlichen Abbildung einer Sequenz;
    • 4 zeigt ein Filter und ein Korrelationssignal mit Hauptpeak, Gleichanteil und niederfrequentem Untergrund;
    • 5 zeigt ein Balanciertes Filter für eine Sequenz LABS26 und Ausgangssignal mit Hauptmaximum und Nebenmaxima;
    • 6 zeigt eine laterale Transportgeschwindigkeit von Zellen in einem fluidischen Flussprofil in einem Kanal;
    • 7 zeigt eine balancierte Filterung eines Eingangssignals mit einer Vielzahl von Filtern;
    • 8 zeigt für die Sequenz LABS26 einen bipolaren, gleichanteilfreien Signalverlauf, einen Filterverlauf und ein Ausgangssignal;
    • 9 zeigt ein unipolares, gleichanteilbehaftetes Signal, einen Filterverlauf und ein Ausgangssignal für die Sequenz Q1DP26(12,#1);
    • 10 zeigt einen Vergleich der Qualitätsmesszahl Q1DI für Sequenzen aus einer Klasse von LABS und ULABS;
    • 11 zeigt in einem beispielhaften Ausführungsbeispiel die jeweils besten untersuchten Sequenzen aus einer Klasse von LABS und ULABS bewertet nach einer Qualitätsmesszahl Q1DI;
    • 12 zeigt einen Vergleich einer Qualitätsmesszahl Q1DP für Sequenzen aus einer Klasse von LABS und ULABS;
    • 13 zeigt in einem beispielhaften Ausführungsbeispiel die jeweils besten Sequenzen aus einer Klasse von LABS und ULABS bewertet nach einer Qualitätsmesszahl Q1DP,
    • 14 zeigt einen Vergleich der jeweils besten untersuchten Sequenzen aus Ausführungsbeispielen mit L=33 bezüglich der Qualitätsmesszahlen Q1DP und Q1DI;
    • 15 zeigt einen Vergleich von 2D quadrierten Korrelationskanälen für ausgewählte Geschwindigkeitskanäle für LABS26(16,#2) und Q2DPs26(15,#2);
    • 16a und 16b zeigen ein Ausschneiden der Hauptmaxima zum Bestimmen der bereinigten Korrelationsfunktionen in Ausführungsbeispielen beispielhaft für ULABS-Q2DPs31(16,#1);
    • 17a und 17b zeigen ein Ausschneiden der Hauptmaxima zum Bestimmen der bereinigten Korrelationsfunktionen in Ausführungsbeispielen für LABS-Q2DPs31(19,#1);
    • 18 zeigt einen Vergleich der Qualitätsmesszahl Q2DPs für bekannte LABS und Beispiele aus der Klasse von ULABS aus Ausführungsbeispielen;
    • 19 zeigt eine beispielhafte Auflistung von jeweils besten untersuchten Sequenzen aus der Klasse von LABS und ULABS bewertet nach der Qualitätsmesszahl Q2DPs;
    • 20 zeigt einen Vergleich der Qualitätsmesszahl Q2DP für bekannte LABS und Beispiele aus der Klasse von Sequenzen „ULABS“ aus Ausführungsbeispielen;
    • 21 zeigt eine beispielhafte Auflistung von jeweils besten untersuchten Codes aus der Klasse von LABS und ULABS bewertet nach der Qualitätsmesszahl Q2DP;
    • 22 zeigt eine beispielhafte Eingangsverteilung;
    • 23 zeigt eine beispielhafte Ausgangsverteilung für LABS26(16,#2) mit einer Positiverkennungsrate von 99,6%;
    • 24 zeigt eine beispielhafte Ausgangsverteilung für Q2DPs26(15,#2) mit einer Positiverkennungsrate von 99,8%;
    • 25 zeigt eine beispielhafte Ausgangsverteilung für Q1DI26(14,#1) mit einer Positiverkennungsrate von 99,8%;
    • 26 zeigt eine beispielhafte Ausgangsverteilung für B13Man26(13) mit einer Positiverkennungsrate von 99,8%.
    • 27 zeigt beispielhaft eine Erkennungsrate als Funktion des mittleren Partikelabstandes in Einheiten der Maskenlänge;
    • 28 zeigt eine Fehlerkennungsrate von Ausführungsbeispielen im Rauschen als Funktion des mittleren Partikelabstandes in Einheiten der Maskenlänge
    • 29 zeigt einen beispielhaften Vergleich einer Nebenmaximahäufigkeit;
    • 30 zeigt eine mittlere Strukturgröße für die besten untersuchten 2D-Codes gemäß eines Qualitätskriteriums aus Ausführungsbeispielen und LABS;
    • 31 zeigt ein Ausführungsbeispiel einer Vorrichtung für ein Zytometer;
    • 32 zeigt ein Blockdiagramm eines Ausführungsbeispiels eines Zytometers;
    • 33 zeigt ein Flussdiagramm eines Verfahrens zum Bereitstellen von Information über zumindest eine Sequenz; und
    • 34 zeigt ein Flussdiagramm eines Verfahrens zur Bereitstellung von Information über ein oder mehrere Zellen in einem Medium in einem Kanal.
  • Beschreibung
  • Verschiedene Ausführungsbeispiele werden nun ausführlicher unter Bezugnahme auf die beiliegenden Zeichnungen beschrieben, in denen einige Ausführungsbeispiele dargestellt sind. In den Figuren können die Dickenabmessungen von Linien, Schichten und/oder Regionen um der Deutlichkeit Willen übertrieben dargestellt sein.
  • Bei der nachfolgenden Beschreibung der beigefügten Figuren, die lediglich einige exemplarische Ausführungsbeispiele zeigen, können gleiche Bezugszeichen gleiche oder vergleichbare Komponenten bezeichnen. Ferner können zusammenfassende Bezugszeichen für Komponenten und Objekte verwendet werden, die mehrfach in einem Ausführungsbeispiel oder in einer Zeichnung auftreten, jedoch hinsichtlich eines oder mehrerer Merkmale gemeinsam beschrieben werden. Komponenten oder Objekte, die mit gleichen oder zusammenfassenden Bezugszeichen beschrieben werden, können hinsichtlich einzelner, mehrerer oder aller Merkmale, beispielsweise ihrer Dimensionierungen, gleich, jedoch gegebenenfalls auch unterschiedlich ausgeführt sein, sofern sich aus der Beschreibung nicht etwas anderes explizit oder implizit ergibt.
  • Obwohl Ausführungsbeispiele auf verschiedene Weise modifiziert und abgeändert werden können, sind Ausführungsbeispiele in den Figuren als Beispiele dargestellt und werden hierin ausführlich beschrieben. Es sei jedoch klargestellt, dass nicht beabsichtigt ist, Ausführungsbeispiele auf die jeweils offenbarten Formen zu beschränken, sondern dass Ausführungsbeispiele vielmehr sämtliche funktionale und/oder strukturelle Modifikationen, Äquivalente und Alternativen, die im Bereich der Erfindung liegen, abdecken sollen. Gleiche Bezugszeichen bezeichnen in der gesamten Figurenbeschreibung gleiche oder ähnliche Elemente.
  • Man beachte, dass ein Element, das als mit einem anderen Element „verbunden“ oder „verkoppelt“ bezeichnet wird, mit dem anderen Element direkt verbunden oder verkoppelt sein kann oder dass dazwischenliegende Elemente vorhanden sein können. Wenn ein Element dagegen als „direkt verbunden“ oder „direkt verkoppelt“ mit einem anderen Element bezeichnet wird, sind keine dazwischenliegenden Elemente vorhanden. Andere Begriffe, die verwendet werden, um die Beziehung zwischen Elementen zu beschreiben, sollten auf ähnliche Weise interpretiert werden (z.B., „zwischen“ gegenüber „direkt dazwischen“, „angrenzend“ gegenüber „direkt angrenzend“ usw.).
  • Die Terminologie, die hierin verwendet wird, dient nur der Beschreibung bestimmter Ausführungsbeispiele und soll die Ausführungsbeispiele nicht beschränken. Wie hierin verwendet, sollen die Singularformen „einer,‟ „eine‟, „eines‟ und „der, die, das“ auch die Pluralformen beinhalten, solange der Kontext nicht eindeutig etwas anderes angibt. Ferner sei klargestellt, dass die Ausdrücke wie z.B. „beinhaltet“, „beinhaltend“, „aufweist“, „umfasst“, „umfassend“ und/oder „aufweisend“, wie hierin verwendet, das Vorhandensein von genannten Merkmalen, ganzen Zahlen, Schritten, Arbeitsabläufen, Elementen und/oder Komponenten angeben, aber das Vorhandensein oder die Hinzufügung von einem bzw. einer oder mehreren Merkmalen, ganzen Zahlen, Schritten, Arbeitsabläufen, Elementen, Komponenten und/oder Gruppen davon nicht ausschließen.
  • Solange nichts anderes definiert ist, haben sämtliche hierin verwendeten Begriffe (einschließlich von technischen und wissenschaftlichen Begriffen) die gleiche Bedeutung, die ihnen ein Durchschnittsfachmann auf dem Gebiet, zu dem die Ausführungsbeispiele gehören, beimisst. Ferner sei klargestellt, dass Ausdrücke, z.B. diejenigen, die in allgemein verwendeten Wörterbüchern definiert sind, so zu interpretieren sind, als hätten sie die Bedeutung, die mit ihrer Bedeutung im Kontext der einschlägigen Technik konsistent ist, und nicht in einem idealisierten oder übermäßig formalen Sinn zu interpretieren sind, solange dies hierin nicht ausdrücklich definiert ist.
  • Bei der Zellzählung basierend auf dem Prinzip der Räumlich Modulierten Fluoreszenz (RMF) werden häufig binäre Codes eingesetzt, mit Hilfe derer die räumliche Modulationsabfolge der Beleuchtung oder der Fluoreszenzemission festgelegt wird (siehe 2).
  • 2 zeigt ein Schema des Prinzips der räumlich modulierten Emission. Der Abfolge der offenen und geschlossenen Bereiche der Maske 2002 wird durch sogenannte binäre Sequenzen/Codes beschrieben, wobei üblicherweise eine „1“ im Code für eine Maskenöffnung und eine „0“ für geschlossene Bereiche verwendet wird. 2 zeigt zusätzlich einen Sensor 2004, der ein Lichtsignal von fluoreszierenden Partikeln 2006 in einem Kanal 2008 erfasst und als gefiltertes Signal 2010 bereitstellt.
  • Im nicht-prozessierten Ausgangssignal 2010 der RMF kann sich die räumliche Modulationsabfolge in einer entsprechenden zeitlichen Modulation der Signalamplitude (Nutzsignal) wiederfinden (vgl. 3). Reale Signale können aus einer additiven Überlagerung des Nutzsignales, eines Gleichanteils und von Rauschen, das aus verschiedenen Quellen stammen kann (z.B. Detektorrauschen, Schrotrauschen, Quantisierungsrauschen, etc.) bestehen. 3 zeigt ein Beispiel für ein Ausgangssignal 3002 codiert mit der Sequenz LABS26(16) (S11011011010101111110001100T). LABS steht für Low-Autocorrelation Binary Sequences, Schwach-Autokorrelierte Binäre Sequenzen, eine bekannte Gruppe von Sequenzen mit schwacher Autokorrelation. 3004 zeigt das entsprechende Filter.
  • Zur Bestimmung des Transitzeitpunktes und der Fluoreszenzintensität der Zellen, könnte bei der RMF eine kontinuierliche Korrelation 4004 des Eingangssignals mit der zugrundeliegenden Sequenz (Filter, 4002), die auf die Länge des Messsignals gestreckt wird (Kiesel et al. 2009) erfolgen. Das Korrelations-/Filterungsergebnis 4004 könnte anschließend differenziert werden, um Gleichanteile sowie den niederfrequenten Untergrund zu unterdrücken 4006. 4 zeigt Filter 4002 und Korrelationssignal 4004 mit Hauptpeak, Gleichanteil und niederfrequentem Untergrund. 4006 zeigt ein zeitlich differenziertes Korrelationssignal zur Eliminierung der Gleichanteile des niederfrequenten Untergrundes (Differentielle Filterung).
  • Im Ausgangssignal können neben einem Hauptmaximum, die das Zellereignis wiedergibt, signifikante Nebenmaxima verbleiben, die die Erkennung der Hauptsignatur im Fall von überlappenden Partikeln erschweren können. In Ausführungsbeispielen können Nebenmaxima beispielsweise lokalen Maxima in einer Korrelationsfunktion entsprechen, die nicht einem Hauptmaximum zugeordnet werden. Die lokalen Maxima können beispielsweise auf einer Betragsfunktion basieren, und ein Maximum für eine (positive oder negative) Abweichung von einem Referenzwert angeben. In manchen Ausführungsbeispielen kann ein Nebenmaximum auch einer Nebenkeule entsprechen.
  • 5 zeigt beispielhaft, wie mittels einer balancierten Filterung in einem Korrelationsschritt ein Ausgangssignal 5004 mit reduzierter Rauschverstärkung aus dem Eingangssignal der 3 3002 gewonnen werden kann. Mit dieser Methode könnte der Rechenaufwand reduziert werden und/oder das Signal-zu-Rausch-Verhältnis (auch engl. Signal-to-Noise-Ratio, SNR) verbessert werden. Allerdings könnten Nebenmaxima mit vergleichbarer Amplitude im Signal verbleiben. 5 zeigt ein Balanciertes Filter 5002 für die Sequenz LABS26 und Ausgangssignal 5004 mit Hauptmaximum und Nebenmaxima. 5006 zeigt den Absolutbetrag des Ausgangssignals. Im Vergleich zur 4 ist der Gleichanteil unmittelbar eliminiert und der niederfrequente Untergrund ist deutlich reduziert. Im Vergleich zu 4006 ist das Ausganssignal sichtbar weniger verrauscht, allerdings können Nebenmaxima vergleichbarer Amplitude im Signal verbleiben.
  • Die zeitliche Ausdehnung des Signales ist häufig nicht bekannt, da z.B. bei einer mikrofluidischen Zellzählung die laterale Transportgeschwindigkeit von Zellen 6004 von ihrer Transportposition im fluidischen Flussprofil 6002 abhängt (siehe 6). Um dennoch Zellen identifizieren zu können, könnte nun die Korrelation des Signales mit einer Vielzahl von Filtern so erfolgen, dass der gesamte erwartete Geschwindigkeits- / Laufzeitbereich der Objekte abgedeckt wird. Die Filterung könnte beispielsweise auf einer Vielzahl von sogenannten Geschwindigkeitskanälen 7002 (vgl. 7) erfolgen und wird auch zweidimensionale (2D-) Korrelation genannt.
  • 7 zeigt eine balancierte Filterung eines Eingangssignals 7000 mit einer Vielzahl von Filtern 7004 auf den Geschwindigkeitskanälen 7002. Die Filter 7004 sind zeitlich gestreckt um die tatsächliche Länge des Eingangssignals 7000 zu ermitteln. Zudem sind die Amplituden der Filter normiert, um ausgangsseitig die Amplitude des Eingangssignales zu reproduzieren. Aus dem globalen Maximum über alle Korrelations-/Ausgangsignale 7006 ergeben sich die Eingangssignalamplitude, die Transportgeschwindigkeit des Objektes und der Zeitpunkt des Transits.
  • In den Ausgangssignalen könnte beispielsweise aus dem globalen Maximum (größten Hauptmaximum) ein sogenannter resonanter Geschwindigkeitskanal ermittelt werden. In zumindest manchen Ausführungsbeispielen kann ein Geschwindigkeitskanal, oder auch Kanal, einer zeitlich skalierten Sequenz oder einem Filter basierend auf einer zeitlich skalierten Sequenz entsprechen. In diesem Kanal stimmt die Länge des Eingangssignals mit der Länge des Filters überein und es könnten beispielsweise die Signalamplitude, die Geschwindigkeit und der Zeitpunkt des Zelltransits ermittelt werden. Allerdings können die nicht resonanten Kanäle ausgeprägte Maxima zeigen, die in ihrer Amplitude z.T. deutlich höher sind als die Nebenmaxima auf dem resonanten Kanal. Bei der 2D-Korrelation können bei Verwendung der bekannten LABS auf den nicht resonanten Kanälen ausgeprägte Nebenmaxima auftreten. Dies hat drei wesentliche Folgen:
    • ■ Rauschen im Eingangssignal kann dazu führen, dass Nebenmaxima deutlich über ihren regulären Level überhöht werden, sodass Fehlerkennungen auftreten. Schwellenwerte für die Ereigniserkennung könnten hinreichend hoch gesetzt werden um Fehlerkennungen zu vermeiden. Das Detektionslimit ist somit durch die Amplitude der größten Nebenmaxima gegeben.
    • ■ Bei überlappenden Eingangssignalen, d.h. bei gleichzeitigem Transit mehrerer Zellen durch die Detektionszone könnten die Amplituden der Hauptmaxima stark durch die Nebenmaxima der benachbarten Ereignisse verfälscht werden, oder es könnten durch die konstruktive Überlagerung von Nebenmaxima Fehlerkennungen zustande kommen.
    • ■ Bei überlappenden Zellsignalen könnte die Amplitude des schwächeren Signales immer noch deutlich geringer als die größten Nebenmaxima des stärkeren Signales sein. Die Messbereichsdynamik könnte im Zeitbereich um das Hauptmaximum, in dem sich die Nebenmaxima befinden, stark eingeschränkt sein.
  • Ausführungsbeispiele basieren auf der Erkenntnis, dass es sinnvoll sein kann, die maximale Amplitude der Nebenmaxima für den Fall unterschiedlicher Zellgeschwindigkeiten - also für den 2D Fall - auch auf den nicht resonanten Geschwindigkeitskanälen zu minimieren. Zumindest manche Ausführungsbeispiele beschreiben, nach welchen Kriterien Codes auszuwählen sind, die diese Forderung erfüllen.
  • 1 zeigt ein Blockdiagramm eines Ausführungsbeispiels einer Vorrichtung 20 zum Bereitstellen von Information über zumindest eine Sequenz. Die zumindest eine Sequenz beschreibt zeitlich aufeinanderfolgende Signalzustände. Ein Signalzustand könnte einem binären, einem unipolaren, einem bipolaren oder einem multipolaren Signalzustand entsprechen. Ein Signalzustand könnte einen logischen oder diskreten Wert umfassen oder auf einer Spannung, einer Stromstärke und/oder einer Messintensität basieren. In zumindest manchen Ausführungsbeispielen könnte die Sequenz zeitlich äquidistante Signalzustände beschreiben, oder die Sequenz könnte einen variablen zeitlichen Zusammenhang zwischen Signalzuständen beschreiben. In zumindest manchen Ausführungsbeispielen kann eine Sequenz einer Binärsequenz entsprechen. Die Binärsequenz könnte beispielsweise zeitlich äquidistante aufeinanderfolgende Signalzustände beschreiben.
  • Die Vorrichtung 20 umfasst eine Schnittstelle 22, die ausgebildet ist, um Information über eine Anzahl der Signalzustände zu erhalten. Die Schnittstelle 22 kann beispielsweise einem oder mehreren Eingängen und/oder einem oder mehreren Ausgängen zum Empfangen und/oder Übertragen von Informationen entsprechen, etwa in digitalen Bitwerten, basierend auf einem Code, innerhalb eines Moduls, zwischen Modulen, oder zwischen Modulen verschiedener Entitäten. Die Information über die Anzahl der Signalzustände könnte beispielsweise Information über eine Anzahl von Zeitschritten in der Sequenz umfassen, Information über eine Anzahl von Signalzustandswechseln und/oder Information über eine Länge der Sequenz.
  • Die Vorrichtung 20 umfasst ferner ein Berechnungsmodul 24, welches ausgebildet ist, um eine Mehrzahl von möglichen Sequenzen zu generieren, basierend auf der Information über die Anzahl der Signalzustände. In einigen Ausführungsbeispielen könnte die Mehrzahl von möglichen Sequenzen alle Sequenzen umfassen, die basierend auf der Anzahl der Signalzustände und basierend auf den Signalzuständen möglich sind. Bei binären Signalzuständen und einer Anzahl von 5 Signalzuständen könnte die Mehrzahl von möglichen Sequenzen beispielsweise 32 Sequenzen umfassen. In zumindest einigen anderen Ausführungsbeispielen könnte die Mehrzahl von möglichen Sequenzen eine Untermenge von möglichen Sequenzen umfassen, die beispielsweise durch Eigenschaften der Sequenzen eingeschränkt ist, beispielsweise durch eine Periodizität oder eine mittlere Strukturlänge.
  • In Ausführungsbeispielen kann das Berechnungsmodul 24 einem beliebigen Controller oder Prozessor oder einer programmierbaren Hardwarekomponente entsprechen. Beispielsweise kann das Berechnungsmodul 24 auch als Software realisiert sein, die für eine entsprechende Hardwarekomponente programmiert ist. Insofern kann das Berechnungsmodul 24 als programmierbare Hardware mit entsprechend angepasster Software implementiert sein. Dabei können beliebige Prozessoren, wie Digitale Signalprozessoren (DSPs) zum Einsatz kommen. Ausführungsbeispiele sind dabei nicht auf einen bestimmten Typ von Prozessor eingeschränkt. Es sind beliebige Prozessoren oder auch mehrere Prozessoren zur Implementierung des Berechnungsmoduls 24 denkbar.
  • In zumindest einigen Ausführungsbeispielen könnte eine Sequenz eine mittlere Strukturgröße haben. Eine Strukturgröße könnte auf einer Anzahl von gleichen aufeinanderfolgenden Signalzuständen basieren. Das Berechnungsmodul 24 könnte ausgebildet sein, um die Mehrzahl von möglichen Sequenzen basierend auf einem Zielbereich für die mittlere Strukturgröße zu generieren.
  • Die mittlere Strukturgröße kann einem Durchschnitt der Strukturgrößen einer Sequenz entsprechen. In zumindest manchen Ausführungsbeispielen könnte der Zielbereich für die mittlere Strukturgröße zwischen 1,3 und 1,8 liegen. In einigen Ausführungsbeispielen könnte der Zielbereich beispielsweise zwischen 1,4 und 1,7, zwischen 1,5 und 1,6, zwischen 1,4 und 1,7, zwischen 1,5 und 1,7, oder zwischen 1,4 und 1,6 liegen.
  • Das Berechnungsmodul 24 ist ferner ausgebildet, für zumindest eine Teilmenge der möglichen Sequenzen Korrelationsfunktionen zwischen einer Sequenz und zumindest einer zeitlichen Skalierung der Sequenz zu berechnen. Eine Korrelationsfunktion umfasst ein Hauptmaximum und ein oder mehrere Nebenmaxima.
  • Die zumindest eine zeitliche Skalierung der Sequenz könnte beispielsweise auf unterschiedlich langen zeitlichen Ausdehnungen der Signalzustände basieren. In zumindest manchen Ausführungsbeispielen kann die zumindest eine zeitliche Skalierung der Sequenz auf einer Abbildung der Sequenz auf unterschiedliche Anzahlen von Abtastzeitpunkten oder Vergleichszeitpunkten basieren. Beispielsweise könnte die zumindest eine zeitliche Skalierung der Sequenz zeitlich skalierte Sequenzen mit 2 Abtastzeitpunkten oder Vergleichszeitpunkten (auch engl. Samples) pro Zustandsänderung, 3 Samples pro Zustandsänderung, 4 Samples pro Zustandsänderung etc. umfassen.
  • Die Korrelationsfunktionen können beispielsweise Kreuzkorrelationsfunktionen entsprechen. Die Sequenz könnte beispielsweise auf einer zeitlichen Referenzskalierung basieren, und das Berechnungsmodul 24 könnte ausgebildet sein, die Kreuzkorrelationsfunktionen zwischen der zeitlichen Referenzskalierung der Sequenz und der zumindest einen zeitlichen Skalierung der Sequenz zu berechnen. In zumindest manchen Ausführungsbeispielen könnten die Korrelationsfunktionen Autokorrelationsfunktionen umfassen. In einigen Ausführungsbeispielen könnte das Berechnungsmodul 24 ferner ausgebildet sein, die Autokorrelationsfunktion der Sequenz basierend auf der zeitlichen Referenzskalierung und/oder Autokorrelationsfunktionen der zeitlich skalierten Sequenzen zu berechnen. In manchen Ausführungsbeispielen, in denen die zumindest eine zeitliche Skalierung der Sequenz zwei oder mehr unterschiedliche zeitliche Skalierungen der Sequenz umfasst, könnte das Berechnungsmodul 24 ferner ausgebildet sein, Kreuzkorrelationsfunktionen zwischen unterschiedlichen zeitlich skalierten Sequenzen der zwei oder mehr unterschiedlichen zeitlichen Skalierungen der Sequenz zu berechnen.
  • Das Berechnungsmodul 24 könnte ausgebildet sein, die Teilmenge der möglichen Sequenzen basierend auf der mittleren Strukturgröße zu bestimmen. Alternativ oder zusätzlich könnte das Berechnungsmodul 24 ausgebildet sein, die Teilmenge der möglichen Sequenzen basierend auf dem Zielbereich für die mittlere Strukturgröße zu bestimmen.
  • Das Berechnungsmodul 24 ist ferner ausgebildet, die zumindest eine Sequenz basierend auf den Korrelationsfunktionen zu bestimmen. Die Folge der Signalzustände innerhalb der zumindest einen Sequenz ist derart gewählt, dass ein Nebenmaximum in einer Korrelationsfunktion der Sequenz mit der zumindest einen zeitlichen Skalierung der Sequenz gegenüber einem durch unterschiedliche Anordnungen der Signalzustände in der Sequenz maximal in einer Korrelationsfunktion erzielbaren Nebenmaximum reduziert ist.
  • Das Berechnungsmodul 24 ist ferner ausgebildet, die Information über die zumindest eine Sequenz basierend auf der zumindest einen Sequenz zu bestimmen und über die Schnittstelle 22 bereitzustellen. Die Schnittstelle 22 ist mit dem Berechnungsmodul 24 gekoppelt.
  • In zumindest manchen Ausführungsbeispielen könnte das Berechnungsmodul 24 ausgebildet sein, um die Mehrzahl von möglichen Sequenzen als unipolare Sequenzen zu bestimmen. In einigen Ausführungsbeispielen könnte das Berechnungsmodul 24 ferner ausgebildet sein, die Korrelationsfunktionen basierend auf den unipolaren Sequenzen zu berechnen und/oder die zumindest eine Sequenz basierend auf den unipolaren Sequenzen zu bestimmen. In zumindest einigen Ausführungsbeispielen könnte das Berechnungsmodul 24 ferner ausgebildet sein, die zeitlich skalierten Sequenzen als balancierte zeitlich-skalierte Sequenzen zu bestimmen, die zeitlich skalierten Sequenzen als normierte zeitlich-skalierte Sequenzen zu bestimmen, und/oder die Korrelationsfunktionen basierend auf den balancierten zeitlich-skalierten Sequenzen und/oder den normierten balancierten zeitlich-skalierten Sequenzen zu berechnen.
  • In einigen Ausführungsbeispielen könnte das Berechnungsmodul 24 ausgebildet sein, um die zumindest eine Sequenz basierend auf den Korrelationsfunktionen zu bestimmen. Die Folge der Signalzustände innerhalb der zumindest einen Sequenz könnte derart gewählt sein, dass eine Summe, die auf den Nebenmaxima einer Korrelationsfunktion der Sequenz mit der zumindest einen zeitlichen Skalierung der Sequenz basiert, gegenüber einer durch unterschiedliche Anordnungen der Signalzustände in der Sequenz maximal erzielbaren Summe, die auf aus der unterschiedlichen Anordnung resultierenden Nebenmaxima basiert, reduziert ist. Die Summen könnten beispielsweise auf numerischen Repräsentationen der Korrelationsfunktionen basieren. Die Summen könnten beispielsweise auf Quadrierungs- oder Beitragsfunktionen der numerischen Repräsentationen der Nebenmaxima basieren, beispielsweise geteilt durch die Anzahl von betrachteten Zeitschritten. Die Summen könnten beispielsweise normiert sein.
  • Alternativ oder zusätzlich könnte die Folge der Signalzustände innerhalb der zumindest einen Sequenz derart gewählt sein, dass ein größtes Nebenmaximum in der Korrelationsfunktion der Sequenz mit der zumindest einen zeitlichen Skalierung der Sequenz gegenüber einem durch unterschiedliche Anordnungen der Signalzustände in der Sequenz maximal erzielbaren größten Nebenmaximum reduziert ist. Das größte Nebenmaximum könnte beispielsweise auf einer Quadrierungs- oder Beitragsfunktion der numerischen Repräsentation der Nebenmaxima basieren.
  • In einigen Ausführungsbeispielen könnte das Berechnungsmodul 24 ausgebildet sein, um die zumindest eine Sequenz so zu bestimmen, dass die Summe, die auf den Nebenmaxima basiert, unter einem ersten Schwellenwert liegt.
  • Der erste Schwellenwert könnte beispielsweise einer zweitniedrigsten Summe, die auf den Nebenmaxima der Korrelationsfunktionen basiert, entsprechen. Das Berechnungsmodul 24 könnte die zweitniedrigste Summe beispielsweise rekursiv berechnen. Beispielsweise könnte das Berechnungsmodul 24 die Korrelationsfunktionen iterativ berechnen, wobei das Berechnungsmodul 24 die zu diesem Zeitpunkt niedrigste Summe, die auf den Nebenmaxima der Korrelationsfunktionen basiert, als ersten Schwellenwert nutzt.
  • In zumindest manchen Ausführungsbeispielen könnte der erste Schwellenwert auf einer niedrigsten Summe, die auf Nebenmaxima von Korrelationsfunktionen zwischen Sequenzen und zeitlich skalierten Sequenzen basierend auf zumindest einem Element der Gruppe von Low-Autocorrelation Binary Sequences (Binäre Sequenzen mit geringer Autokorrelation, LABS), Barker-Codes, CS100 codes, S100 codes, Neuman-Hofman Codes, Merten's Codes basieren.
  • Das Berechnungsmodul 24 könnte ausgebildet sein, um die zumindest eine Sequenz so zu bestimmen, dass das größte Nebenmaximum unter einem zweiten Schwellenwert liegt. Der zweite Schwellenwert könnte einem größten Nebenmaximum einer Sequenz mit einem zweitniedrigsten größten Nebenmaximum der Korrelationsfunktionen entsprechen. Das Berechnungsmodul 24 könnte das zweitniedrigste größte Nebenmaximum beispielsweise iterativ oder rekursiv berechnen. Beispielsweise könnte das Berechnungsmodul 24 die Korrelationsfunktionen iterativ berechnen, wobei das Berechnungsmodul 24 das zu diesem Zeitpunkt niedrigste größte Nebenmaximum, das auf den Nebenmaxima der Korrelationsfunktionen basiert, als zweiten Schwellenwert nutzt. In zumindest manchen Ausführungsbeispielen könnten die ein oder mehreren Nebenmaxima und/oder das Hauptmaximum normiert sein.
  • In einigen Ausführungsbeispielen könnte ein niedrigstes größtes Nebenmaximum einem Nebenmaximum entsprechen, das, in einer Mehrzahl von Korrelationsfunktionen zwischen Sequenzen und zeitlich skalierten Sequenzen einer Mehrzahl verschiedener Sequenzen, wobei die Korrelationsfunktionen jeweils ein größtes Nebenmaximum umfassen, das niedrigste größte Nebenmaximum darstellt.
  • In manchen Ausführungsbeispielen könnte der erste Schwellenwert auf einem niedrigsten größten Nebenmaximum von Korrelationsfunktionen zwischen Sequenzen und zeitlich skalierten Sequenzen basierend auf zumindest einem Element der Gruppe von Low-Autocorrelation Binary Sequences (Binäre Sequenzen mit geringer Autokorrelation, LABS), Barker-Codes, CS100 codes, S100 codes, Neuman-Hofman Codes, Merten's Codes basieren.
  • In konventionellen Systemen werden auch für unipolare Signale vorwiegend sog. „Low Autocorrelation Binary Sequences“ (LABS) zur Definition des Signalverlaufes eingesetzt. Diese Sequenzen sind bipolar, d.h. die Koeffizienten sind entweder „1“ oder „-1‟ und können beispielsweise für bipolare und gleichanteilfreie Signale optimiert sein. Die Optimierung kann beispielsweise darauf ausgelegt sein entweder die Gesamtenergie (oder Amplitude) in den Nebenmaxima oder die Energie/Amplitude des größten Nebenmaximums zu minimieren. Die oben bereits eingeführte Sequenz LABS26 ist eine solche Sequenz mit der Länge 26. Bei der Umwandlung einer bipolaren (1,-1) Sequenz in eine unipolare (0,1) Sequenz werden im Folgenden zwei Transformationsvorschriften genutzt: 1 0,     1 1
    Figure DE102015110316B4_0001
    1 1,      1 0
    Figure DE102015110316B4_0002
  • D.h. aus einer bipolaren Sequenz könnten zwei unipolare Sequenzen erzeugt werden, die sich im Sinne der im Folgenden beschriebenen Eigenschaften unterschiedlich verhalten können. Im Folgenden wurden jeweils beide Sequenzen untersucht.
  • In 8 sind für die Sequenz LABS26 der bipolare, gleichanteilfreie Signalverlauf 8002, der Filterverlauf 8004 und das Ausgangssignal 8006 dargestellt. Im Vergleich zu 4 fällt deutlich auf, dass das Hauptmaximum (AH) wesentlich größer als das höchste Nebenmaximum (max(abs(NKi)), NKi bezeichnet das i-te Nebenmaximum) ist. In diesem Fall ergibt sich für die bipolare Situation das Verhältnis V b i p = A H m a x ( a b s ( N K i ) ) = 77,4
    Figure DE102015110316B4_0003
  • Für die beiden unipolaren Fälle ergibt sich Vdifferentiell = 15,1 und Vbalanciert = 13,0.
  • Im Fall der RMF ist die Form der Signale nun unipolar und gleichanteilbehaftet, wie das Beispielsignal 3002 aus 3 zeigt. Zumindest manche Ausführungsbeispiele zeigen, dass im Fall solcher unipolarer Signale eine neue Klasse von Sequenzen, die sich von den bereits bekannten LABS unterscheidet, vorteilhaft sein kann. 2 zeigt anhand eines Beispielsequenz (Q1DP26(12,#1) [01001001010111011000111000]), dass ein wesentlich besseres Hauptmaximum zu Nebenmaximaverhältnis V erzielt werden kann als mit literaturbekannten bipolaren LABS (Levanon et al. 2004, Turyn et al. 1961, Mertens et al. 1996, Coxson et al., 2005, Borewein et al. 2000, Nunn et al. 2008, Shanmugam et al. 2008). So beträgt beispielsweise für die Codelänge 26 ein in Ausführungsbeispielen gefundenes verbessertes Verhältnis Vbal(Q1DP26(12, #1)) = 46,0 unter Verwendung der Sequenz Q1DP26(12,#1). Für die bekannten LABS ergibt sich das beste gefundene Verhältnis zu Vbal(LABS26) = 13,0. In Ausführungsbeispielen kann ein Code einer Sequenz entsprechen.
  • 9 zeigt das unipolare, gleichanteilbehaftete Signal 9002, den Filterverlauf 9004 und das Ausgangssignal 9006 für eine Sequenz Q1DP26(12,#1) aus Ausführungsbeispielen.
  • Im Folgenden wird beispielhaft ein Bewertungsmodell für die Sequenz-/Codequalität für den Fall bekannter Signallänge (1D, bzw. ein-dimensionaler Fall) definiert. Eine unipolare Sequenz (Code) C der Länge L ist definiert durch L Koeffizienten bi : C = [ b i ] ,   ( i = 1,2, ,L ) b i   { 1,   0 }
    Figure DE102015110316B4_0004
  • Die Anzahl E der Einsen im Code berechnet sich dann nach E = i = 1 L b i
    Figure DE102015110316B4_0005
    und die Anzahl Z der Nullen im Code ergibt sich zu Z = L E
    Figure DE102015110316B4_0006
  • Ein Ausschnitt Mj der Länge KM aus dem Messsignal für die Zeitschritte tj, tj+1, ... , ti+KM-1 ist definiert durch eine Folge von KM zeitlich äquidistanten Werten mj+k-1, ( k = 1, 2, ..., KM) die sich jeweils additiv zusammensetzen aus einem quasi stationären Gleichanteil gj+k-1 = g, einem Rauschanteil rj+k-t und dem eigentlichen Nutzsignal S mit den Einzelwerten aj+k-1, die die Nutzsignalamplitude A einnehmen, falls im Zeitschritt tj+k-1 gilt: bi = 1. D.h. zum Zeitpunkt tj+k-1 befindet sich das fluoreszierende Partikel 2006 an einem Ort relativ zur Maske 2002, an dem der Detektor oder Sensor 2004 das Signal des Partikels „sieht“ (vgl. 2). Ansonsten gilt aj+k-1 = 0. Also ergibt sich: M j = [ m j + k 1 ] = [ a j + k 1 + G + r j + k 1 ] , ( j = 1,2 ) ( k = 1,2, ,K M ) a j + k 1 { A ,0 }
    Figure DE102015110316B4_0007
  • Die Gesamtlänge Kv des Nutzsignals S im Messsignal Mj berechnet sich aus der mittleren Anzahl der Zeitschritten Kbv während derer sich das Partikel im Bereich eines Codekoeffizienten bi befindet und der Codelänge L nach K v = K b v L
    Figure DE102015110316B4_0008
  • In zumindest manchen Ausführungsbeispielen muss Kbv nicht ganzzahlig sein. Zur Vereinfachung werden die folgenden Betrachtungen dennoch nur für ganzzahlige Kbv durchgeführt.
  • Das Filter Fv mit den Koeffizienten fk·i,v ist mit der Länge Kv des Nutzsignals zu wählen um eine verbesserte Erkennung zu gewährleisten und seine Koeffizienten werden nach folgender Vorschrift belegt. Für jeden Codekoeffizienten bi könnten Kbv Filterkoeffizienten fk·i,v so gewählt werden, dass gilt: f k i , v = f p , v = { h = 1 E K b v , f a l l s   b i = 1 , ( p = k i l = 1 Z K b v , f a l l s   b i = 0                     = 1,2, ,K v )     
    Figure DE102015110316B4_0009
  • Durch diese Wahl kann das Filter Fv in Ausführungsbeispielen balanciert sein, d.h. es gilt p f p , v = 0,   ( p = 1,2, ,K v )
    Figure DE102015110316B4_0010
    und das Filter Fv kann normiert sein, so dass das Hauptmaximum nach der Korrelation den Wert A der Nutzsignalamplitude im Nutzsignal S annimmt (rauschfreier Fall). Falls also im Messsignalausschnitt Mj ein Nutzsignal mit der Amplitude A enthalten ist, dann kann für das (rauschfreie) Korrelationssignal C gelten m a x ( c q ) j , v = m a x ( p = 1 K v m j + q 1 + p 1 f p , v ) = A , ( q = 1,2, ,K M K v + 1 )
    Figure DE102015110316B4_0011
  • Allgemein sind dabei die Koeffizienten der Korrelation Cj,v = Mj × Fv = [cq,j,v] des Messsignals Mj der Länge KM mit dem Filter Fv folgendermaßen definiert: c q , j , v = p = 1 K v m j + q 1 + p 1 f p , v ,    ( q = 1,2, ,K M K v + 1 )
    Figure DE102015110316B4_0012
  • Damit in der Korrelation alle Koeffizienten, d.h. auch alle Nebenmaxima vollständig erfasst sind, könnte das Messsignal zumindest eine Länge KM = 3Kv - 2 haben.
  • Zur Charakterisierung verschiedener Codes könnte in Ausführungsbeispielen das Messsignal nun vereinfacht werden, indem der Gleichanteil g = 0 gesetzt wird und das Rauschen vernachlässigt wird (rj+k-1 = 0∀k), d.h. die Codecharakterisierung könnte rein anhand des Nutzsignales erfolgen. Zur weiteren Vereinfachung wird die Amplitude des Nutzsignales auf A = 1 gesetzt um ein einfaches Filtertestsignal zu erhalten. Das Filtertestsignal ST,v für ein Objekt der Geschwindigkeit v könnte dann die Form haben: S T ,v = [ s p , v ] , w o b e i   s p , v = s k i , v       = { 1, f a l l s   b i = 1 , ( p = 1,2, , K v ) 0, f a l l s   b i = 0
    Figure DE102015110316B4_0013
  • Da unter diesen Bedingungen in Ausführungsbeispielen Koeffizienten außerhalb des Filtertestsignals und des Filters zu Null angenommen werden könnten, könnte die Korrelation beider Signale folgendermaßen dargestellt werden: c q , v = p = 1 K v s p , v f p + q K v ,        ( q = 1,2, ,2K v 1 ) ,                         w o b e i   s x , v = f x = 0,   f u ¨ r   x < 1   o d e r   x > K v
    Figure DE102015110316B4_0014
  • Die Korrelation hat beispielsweise 2Kv - 1 Korrelationskoeffizienten. Davon werden per Definition die Koeffizienten mit dem Index Kv - (Kbv - 1), Kv - (Kbv - 1) + 1, Kv + (Kbv - 1), also insgesamt (2Kbv + 1) Koeffizienten, dem Hauptmaximum zugerechnet und alle anderen Koeffizienten gehören zu den Nebenmaxima. Diese Definition unterstellt, dass das Hauptmaximum doppelt so viele Korrelationskoeffizienten enthält, wie Zeitschritte in einem Codekoeffizienten vorliegen. Um positive und negative Werte gleichermaßen zu behandeln könnten die Korrelationskoeffizienten in zumindest einigen Ausführungsbeispielen quadriert werden und das Maximum 1 D P S L R = m a x ( c q , v 2 ) f u ¨ r   q K v ( K b v 1 ) ,K v ( K b v 1 ) + 1, ,K v + ( K b v 1 )
    Figure DE102015110316B4_0015
    könnten einen charakteristischen Wert für das größte Nebenmaximum des Codes angeben, wobei alle Korrelationskoeffizienten in einer Gesamtbreite von (2Kbv + 1) um das Hauptmaximum ausgelassen werden. In zumindest manchen Ausführungsbeispielen könnte das größte Nebenmaximum auf 1DPSLR basieren.
  • Ein Mittelwert 1 D I S L R = q c q , v 2 2 K v 1 2 ( K b v 1 ) f u ¨ r   q K v ( K b v 1 ) ,K v ( K b v 1 ) + 1, ,K v + ( K b v 1 )
    Figure DE102015110316B4_0016
    ist beispielshaft ein weiterer charakteristischer Wert für die Energie/Amplitude die in den Nebenmaxima enthalten ist, wobei ebenfalls alle Korrelationskoeffizienten in einer Gesamtbreite von (2Kbv + 1) um das Hauptmaximum ausgelassen werden. In zumindest manchen Ausführungsbeispielen könnte die Summe der Nebenmaxima auf 1DISLR basieren.
  • Alternativ könnten als Bereich für das Hauptmaximum alle Korrelationskoeffizienten bis zu den beiden benachbarten lokalen Minima herangezogen werden.
  • Ein weiterer Parameter für die Bewertung von Sequenzen für ihre Eignung für die RMS Technik ist in manchen Ausführungsbeispielen die Rauschverstärkung. Allgemein berechnet sich die Rauschverstärkung V zu V v = p = 1 K v f p , v 2
    Figure DE102015110316B4_0017
  • Aus dem Verhältnis von Eingangssignalamplitude (= Amplitude des Hauptmaximums bei balancierten Codes) zu Rauschverstärkung ergibt sich der SNR-Gewinn beispielsweise zu S N R G a i n , v = A p = 1 K v f p , v 2 = 2 1 A p = 1 K v f p 2 = K b E Z L
    Figure DE102015110316B4_0018
  • Der SNR-Gewinn eines balancierten Filters in Verbindung mit unipolaren Modulationssequenzen könnte für E = Z für gerades L oder für abs(E - Z) = 1 für ungerades L ein Maximum einnehmen.
  • Ein Code könnte zwei Bedingungen erfüllen:
    1. 1. Er hat ein verbessertes SNRGain und
    2. 2. er hat ein verkleinertes 1DISLR oder 1DPSLR
  • Mathematisch könnte dies z.B. zusammengefasst werden in den Qualitätsmesszahlen Q 1 D I = S N R G a i n 1 D I S L R
    Figure DE102015110316B4_0019
    und Q 1 D P = S N R G a i n 1 D P S L R
    Figure DE102015110316B4_0020
  • In zumindest manchen Ausführungsbeispielen könnte das größte Nebenmaximum auf Q1DP basieren. In zumindest manchen Ausführungsbeispielen könnte die Summe der Nebenmaxima auf Q1DI basieren.
  • In beispielhaften Ausführungsbeispielen hat eine vollständige Suche in allen Codes bis zur Länge L = 33 ergeben, dass zahlreiche unipolare Codes „ULABS“ existieren, die unter Anwendung dieser Kriterien besser für die Verwendung in der RMF (und vergleichbare Anwendungen, die unipolare Codes nutzen) geeignet sein könnten als die bekannten bipolaren LABS (siehe 10 und 12, die zugehörigen Codes und Qualitätsmesszahlen finden sich in 11 und 13).
  • In zumindest manchen Ausführungsbeispielen beziehen sich die Bezeichnungen beste Codes, optimale Codes, beste Sequenzen oder optimale Sequenzen auf Sequenzen, die bei Berücksichtigung eines oder mehrerer Kriterien oder Qualitätsmesszahlen, etwa Q1DI, Q1DP, Q2DI, Q2DP, Q2DPs, unter einer Mehrzahl von betrachteten Sequenzen vorteilhafte oder die vorteilhaftesten Eigenschaften aufweisen.
  • 10 zeigt einen Vergleich der Qualitätsmesszahl Q1DI für die bekannten LABS (Bestes LABS 1002) und Beispiele aus der Klasse von Codes „ULABS“ aus Ausführungsbeispielen (bester Code: 1004, beste 10: 1006).
  • 11 zeigt in einem beispielhaften Ausführungsbeispiel die jeweils besten betrachteten Sequenzen aus der Klasse der LABS 1116 und der ULABS 1114 bewertet nach der Qualitätsmesszahl Q1DI 1112. 1102 zeigt eine Code-ID der Codes, 1104 eine Binärrepräsentation der Codes, 1106 die Codelänge, 1108 die Anzahl der Einsen im Binärcode, 1110 zeigt SNRGain und 1112 die Qualitätsmesszahl Q1DI.
  • 12 zeigt beispielhaft einen Vergleich der Qualitätsmesszahl Q1DP für die bekannten LABS (bestes LABS: 1202) und Beispiele aus der neu identifizierten Klasse von Codes „ULABS“ (bester Code: 1204, beste 10: 1206).
  • 13 zeigt in einem beispielhaften Ausführungsbeispiel die jeweils besten betrachteten Sequenzen aus der Klasse der LABS 1316 und der ULABS 1314 bewertet nach der Qualitätsmesszahl Q1DP 1312. 1302 zeigt die Code-ID der Codes, 1304 eine Binärrepräsentation der Sequenzen, 1306 die Sequenzlänge, 1308 die Anzahl der Einsen in der Binärsequenz, 1310 zeigt SNRGain und 1312 die Qualitätsmesszahl Q1DP.
  • 14 zeigt einen Vergleich der jeweils besten Sequenzen 1402 (Q1DI) und 1404 (Q1DP) aus Ausführungsbeispielen mit L = 33 aus 11 und 13 bezüglich der Qualitätsmesszahlen Q1DP 1312 und Q1DI 1112. Für beide Kriterien findet sich das selbe LABS 1406.
  • Die bis hierhin ausgeführten Berechnungen und Auswahlkriterien beziehen sich auf den 1D-Fall, d.h. die zeitliche Ausdehnung des Eingangssignals, bzw. die Geschwindigkeit der Zelle / des Objektes ist bekannt. Im Folgenden werden die Kriterien auf den 2D-Fall erweitert. In 7 ist bereits gezeigt, dass in diesem Fall über eine Vielzahl von sogenannten Geschwindigkeitskanälen korreliert werden könnte, d.h. für jede Geschwindigkeit könnte ein speziell angepasster Filter genutzt werden, beispielsweise durch verschiedene zeitliche Skalierungen der Sequenz. Eine entsprechende Filterbank mit einem Satz von Filtern Fj für Kv Geschwindigkeitskanäle könnte in Ausführungsbeispielen unter Verwendung von Gl. (9) berechnet werden. F j = [ f ( k i ) j ] , ( i = 1,2, , L ) , ( k = 1,2, ,K b j ) , ( K b j = K u ,K u + 1, ,K o ) ,
    Figure DE102015110316B4_0021
    wobei Ku · L und Ko · L die Länge des Filters für die höchste, bzw. niedrigste Geschwindigkeit in der Filterbank beschreiben.
  • Analog zu Gl. (13) wird ein Testsignal ST = [sp] mit KT (Ku ≤ KT ≤ Ko) Abtastungen pro Codekoeffizient gewählt und die 2D Korrelation ergibt sich analog Gl. (14) zu c q j = p = 1 K j s p f p + q K j ,            ( q = 1,2, ,2 K j 1 ) , ( K j = L K b j )            w o b e i   s x = f x = 0,   f u ¨ r   x < 1   o d e r   x > K j
    Figure DE102015110316B4_0022
  • 7 zeigt die Form von Filtern 7004 und Korrelationssignalen 7006 (Ausgangssignale) für ausgewählte Geschwindigkeitskanäle 7002. Es zeigt sich, dass in zumindest manchen Ausführungsbeispielen im sogenannten resonanten Kanal das größte Hauptmaximum zu finden ist. Anhand der Parameter Ku und Ko könnte in Ausführungsbeispielen nun ein Geschwindigkeitsbereich festgelegt werden, in dem die Nebenmaxima um das Hauptmaximum bewertet werden könnten. In 7 wurde Ku = 20 (Kanal Nr. 1, vmax,opt) und Ko = 80 (Kanal Nr. 61, vmin,opt) gewählt und der resonante Kanal wurde auf KT = 40 (Kanal Nr. 21, vcenter) gesetzt. Der kürzeste Filter ist in einer beispielhaften Implementierung auf die doppelte Geschwindigkeit des Testpartikels ausgelegt und der längste Filter entsprechend auf die Hälfte und es gilt: v m a x , o p t / 2 = v c e n t e r = 2 v m i n , o p t
    Figure DE102015110316B4_0023
  • Diese Beziehung zeigt an, in welchem Geschwindigkeitsbereich um das Hauptmaximum herum beispielhaft die Nebenmaxima untersucht und minimiert werden könnten. Für die Anwendung der Technik im Fall einer unbekannten Geschwindigkeit bedeutet dies, dass der Code auf eine Geschwindigkeitsdynamik [ v m i n ; v m a x = 2 v m i n ]
    Figure DE102015110316B4_0024
    ausgelegt ist, da dies genau der Geschwindigkeitsbereich ist in dem für das schnellste Objekt die Optimierung bis zur langsamsten Geschwindigkeit ausgelegt ist und umgekehrt.
  • In einigen Ausführungsbeispielen könnte das Berechnungsmodul 24 ausgebildet sein, um bereinigte Korrelationsfunktionen basierend auf den Korrelationsfunktionen zu bestimmen. Das Berechnungsmodul 24 könnte ausgebildet sein, um einen Beitrag von Hauptmaxima in den bereinigten Korrelationsfunktionen zu reduzieren. Das Berechnungsmodul 24 könnte beispielsweise ausgebildet sein, um die Information über die zumindest eine Sequenz basierend auf den bereinigten Korrelationsfunktionen zu bestimmen.
  • In zumindest manchen Ausführungsbeispielen könnte das Berechnungsmodul 24 ausgebildet sein, um unter den Hauptmaxima der Korrelationsfunktionen ein größtes Hauptmaximum zu bestimmen. Das Berechnungsmodul 24 könnte ausgebildet sein, um über die zeitliche Ausdehnung des größten Hauptmaximums den Beitrag der Hauptmaxima in den bereinigten Korrelationsfunktionen zu reduzieren. Beispielsweise könnte das Berechnungsmodul 24 ausgebildet sein, alle anderen Haupt- und Nebenmaxima während der zeitlichen Ausdehnung des größten Hauptmaximums zu unterdrücken oder reduzieren.
  • Alternativ oder zusätzlich könnte das Berechnungsmodul 24 ausgebildet sein, um basierend auf einer zeitlichen Position des größten Hauptmaximums und basierend auf der Korrelationsfunktion der zeitlich skalierten Sequenz die das größte Hauptmaximum umfasst, die Beiträge der Hauptmaxima in den bereinigten Korrelationsfunktionen der zeitlich skalierten Sequenz und benachbarten zeitlichen Skalierungen der Sequenz zu bestimmen und zu reduzieren. Die zeitlich skalierte Sequenz, die das größte Hauptmaximum umfasst könnte beispielsweise einem resonanten Kanal entsprechen. In Ausführungsbeispielen könnte das Berechnungsmodul 24 ausgebildet sein, die Beiträge der Hauptmaxima bis zu einem nächsten (lokalen) Minimum zu reduzieren, beispielsweise basierend auf der zeitlichen Position des größten Hauptmaximums.
  • 15 zeigt einen Vergleich von 2D quadrierten Korrelationskanälen 1506 für ausgewählte Geschwindigkeitskanäle für LABS26(16,#2) 1502 und Q2DPs26(15,#2) 1504. Q2DPs26(15,#2) 1504 zeigt deutlich niedrigere Nebenmaxima über den gesamten Geschwindigkeitsbereich. Zur Ermittlung von Qualitätsmesszahlen für die Codes könnten zwei unterschiedlich scharfe Kriterien genutzt werden:
    • ■ Kriterium 1: Ausschneiden des Zeitbereiches den das größte Hauptmaximum abdeckt auf allen Kanälen (zeitlichen Skalierungen). Vorteil: Dieses Kriterium könnte numerisch einfach umzusetzen sein. Nachteil: Die Optimierung umfasst möglicherweise nicht den Fall von Doppelereignissen, bei denen ein Partikel das andere exakt in der Maskenmitte überholt.
    • ■ Kriterium 2: Ausschneiden des größten Hauptmaximums im resonanten Kanal und der Hauptmaxima in den benachbarten Kanälen (zeitlichen Skalierungen), jeweils bis zum nächsten Minimum. Dieses Kriterium könnte schärfer aber numerisch aufwendiger sein. Für Überholereignisse könnte die Optimierung ebenfalls gelten.
  • 16a und 16b zeigen ein Ausschneiden der Hauptmaxima zum Bestimmen der bereinigten Korrelationsfunktionen in Ausführungsbeispielen beispielhaft für ULABS-Q2DPs31(16,#1). 1602a und 1602b zeigen die Korrelationsfunktionen ULABS-Q2DPs31(16,#1) mit den Hauptmaxima 1604a und 1604b. 1606a und 1606b zeigen die bereinigten Korrelationsfunktionen mit den Leerstellen der entfernten Hauptmaxima 1608a und 1608b.
  • 17a und 17b zeigen ein Ausschneiden der Hauptmaxima zum Bestimmen der bereinigten Korrelationsfunktionen in Ausführungsbeispielen für LABS-Q2DPs31(19,#1). 1702a und 1702b zeigen die Korrelationsfunktionen LABS-Q2DPs31(19,#1) mit den Hauptmaxima 1704a und 1704b. 1706a und 1706b zeigen die bereinigten Korrelationsfunktionen mit den Leerstellen der entfernten Hauptmaxima 1708a und 1708b.
  • Aus der Menge an verbleibenden Korrelationswerten [cq] und [cq]s könnten analog zu Gl. (15) jeweils die charakteristische Werte 2DPSLR für Kriterium 1 und 2DPSLRs für Kriterium 2 berechnet werden, die die Amplitude des höchsten Nebenmaximums bewerten: 2 D P S L R = m a x ( [ c q 2 ] ) 2 D P S L R = m a x ( [ c q 2 ] s )
    Figure DE102015110316B4_0025
  • Für beide Kriterien könnten jeweils gezählt werden, wie viele Korrelationswerte aus der Korrelationsmatrix für das Hauptmaximum entnommen wurden um analog zu Gl. (16) die Mittelwerte als weitere charakteristische Größe zu berechnen: 2 D I S L R = q c q 2 K N K ,      c q a u s g e w a ¨ h l t   n a c h   K r i t e r i u m   1, 2 D I S L R s = q c q 2 K N K s ,      c q a u s g e w a ¨ h l t   n a c h   K r i t e r i u m  2 ,
    Figure DE102015110316B4_0026
  • Wobei KNK und KNKS jeweils die Anzahl der verbleibenden Korrelationswerte angibt.
  • Erneut könnten unter Berücksichtigung des Signal zu Rausch Gewinns SNRGain (siehe Gl. (18)) verschiedene Qualitätsmesszahlen berechnet werden Q 2 D I = S N R G a i n 2 D I S L R Q 2 D I s = S N R G a i n 2 D I S L R s
    Figure DE102015110316B4_0027
    und Q 2 D P = S N R G a i n 2 D P S L R Q 2 D P s = S N R G a i n 2 D P S L R s
    Figure DE102015110316B4_0028
  • In zumindest manchen Ausführungsbeispielen könnte das größte Nebenmaximum auf 2DPSLR, 2DPSLRs, Q2DP und/oder Q2DPs basieren. In zumindest manchen Ausführungsbeispielen könnte die Summe der Nebenmaxima auf 2DISLR, 2DISLRs, Q2DI und/oder Q2DIs basieren.
  • Eine vollständige Suche in allen Codes bis zur Länge L = 33 hat ergeben, dass zahlreiche unipolare Codes existieren, die unter Anwendung dieser Kriterien besser für die Verwendung in der RMF (und vergleichbare Anwendungen, die unipolare Codes nutzen) geeignet sind. Beispielsweise könnten diese Codes/Sequenzen in Ausführungsbeispielen in der Lichtdetektion und Abstandsmessung (auch engl. LIght Detection And Ranging, LIDAR), bei dem die Pulsfolgen zeitlich moduliert werden genutzt werden, oder für Radar-Anwendungen bei denen keine bipolare Codierung möglich ist.
  • 18 zeigt einen Vergleich der Qualitätsmesszahl Q2DPs für die bekannten LABS (bestes LABS: 1802) und die Beispiele aus der Klasse von Codes „ULABS“ aus Ausführungsbeispielen (bester Code: 1804, beste 10: 1806).
  • 19 zeigt eine beispielhafte Auflistung der jeweils besten betrachteten Sequenzen aus der Klasse der LABS und der ULABS bewertet nach der Qualitätsmesszahl Q2DPs. 1902 zeigt die Code-ID der Codes, 1904 eine Binärrepräsentation der Codes, 1906 die Codelänge, 1908 die Anzahl der Einsen im Binärcode, 1910 zeigt SNRGain und 1912 die Qualitätsmesszahl Q2DPs.
  • 20 zeigt einen Vergleich der Qualitätsmesszahl Q2DP für die bekannten LABS (bestes LABS: 2002a) und die Beispiele aus der Klasse von Codes „ULABS“ aus Ausführungsbeispielen (bester Code: 2004a, beste 10: 2006a).
  • 21 zeigt eine beispielhafte Auflistung der jeweils besten Codes aus der Klasse der LABS und der ULABS bewertet nach der Qualitätsmesszahl Q2DP. 2102 zeigt die Code-ID der Codes, 2104 eine Binärrepräsentation der Codes, 2106 die Codelänge, 2108 die Anzahl der Einsen im Binärcode, 2110 zeigt SNRGain und 2112 die Qualitätsmesszahl Q2DP.
  • Zum Vergleich der Qualität ausgewählter Codes in Ausführungsbeispielen wurden anhand von simulierten Ereignissen die Erkennungsrate sowie die Fehlerkennungsrate aufgrund von Nebenmaxima und im Rauschen ermittelt. Dieser Vergleich erfolgte an einer für alle betrachteten Codes einheitlichen Eingangsverteilung (22). Die Eingangsverteilung besteht in Ausführungsbeispielen aus 1000 Partikelsignalen 2202 mit Gaußverteilter Amplitude (0,01+-0,002) und Geschwindigkeit (450+-5mm/s. Für jeden betrachteten Code wurde so ein Rohsignal erzeugt und es wurde jeweils weißes Rauschen mit einer Standardabweichung von 0,002 addiert. Die Rohsignale unterscheiden sich also nur im Verlauf des Nutzsignales, das dem jeweils betrachteten Code folgt.
  • Die Rohsignale wurden jeweils mit balancierten Filterbänken, die einen Geschwindigkeitsbereich von 300mm/s bis 600mm/s abdecken, korreliert. Die Ereigniserkennung erfolgte durch einen einheitlichen Schwellenwertfaktor von 1,3. In einem Zeitfenster von doppelter Laufzeit eines Partikels das mit 300mm/s durch die Messzone transportiert wird die mittlere Amplitude aller Maxima in der quadrierten Korrelationsmatrix ermittelt. Ein Maximum, das um den Schwellenwertfaktor höher liegt als dieser Mittelwert wird in Ausführungsbeispielen als Ereignis angesehen. 23, 24, 25 und 26 zeigen die erkannten Ereignisse von Ausführungsbeispielen für die jeweiligen Codes. Neben den in den obigen Tabellen aufgeführten Codes wurde zusätzlich zum Vergleich der B13Man26(13) hinzugenommen. Dieser Code wird aus dem Barker 13 Code (1111100110101) über eine Manchester Transformation (1 →01, 0→ 10) gewonnen und hat folglich die Länge 26.
  • Der Schwellenwertfaktor wurde in zumindest manchen Ausführungsbeispielen einheitlich gewählt und ist so eingestellt, dass alle Codes mehr als 99% der Ereignisse richtig erkennen. Alle Sequenzen ermöglichen in Ausführungsbeispielen gleichermaßen eine Erkennung von Ereignissen an der Rauschschwelle. Dies ist unproblematisch, da diese über einen globalen Amplitudenschwellenwert aussortiert werden könnten. Es zeigt sich aber, dass der Code Q2DPs26(15,#2) (24) in Ausführungsbeispielen die weitaus wenigsten Fehlerkennungen aufweist. Für diesen Code könnte es also möglich sein weitere Ereignisse zwischen Rauschschwelle und Amplitude der Eingangsverteilung zu erkennen. Andere Sequenzen aus Ausführungsbeispielen könnten hier versagen, da die Fehlerkennungen durch Nebenmaxima weite Amplituden- und Geschwindigkeitsbereiche verdecken könnten.
  • 23 zeigt eine beispielhafte Ausgangsverteilung für das LABS26(16,#2) mit einer Positiverkennungsrate von 99,6%. 2302 (analog 2402, 2502, 2602) zeigt die Ereignisse, 2304 (analog 2404, 2504, 2604) Fehlerkennungen basierend auf einem Nebenmaximum und 2306 (analog 2406, 2506, 2606) Fehlerkennungen basierend auf Rauschen. 24 zeigt eine beispielhafte Ausgangsverteilung für das Q2DPs26(15,#2) mit einer Positiverkennungsrate von 99,8%. 25 zeigt eine beispielhafte Ausgangsverteilung für das Q1DI26(14,#1) mit einer Positiverkennungsrate von 99,8%. 26 zeigt eine beispielhafte Ausgangsverteilung für das B13Man26(13) mit einer Positiverkennungsrate von 99,8%.
  • Wird der mittlere zeitliche Abstand der Ereignisse in der Eingangsverteilung verringert, geht in Ausführungsbeispielen für alle Codes die Erkennungsrate zurück wie beispielhaft in 27 dargestellt. Am schlechtesten schneidet in manchen Ausführungsbeispielen der bereits bekannte LABS26 2706 ab. Die weiteren Codes sind in einer beispielhaften Evaluierung nahezu identisch, wobei leichte Vorteile für den Q2DPs26(15,#2) 2702 vorliegen. 27 zeigt beispielhaft eine Erkennungsrate als Funktion des mittleren Partikelabstandes in Einheiten der Maskenlänge, für Q1DI26(14, #1) 2702 (analog 2802, 2902), Q2DPs26(15, #2) 2704 (analog 2804, 2904), LABS26 (16, #2) 2706 (analog 2806, 2906) und B13Man26(13) 2708 (analog 2808, 2908).
  • Die Fehlerkennungsrate im Rauschen verhält sich weitgehende einheitlich für die betrachteten Codes in betrachteten Ausführungsbeispielen. Die Abnahme mit zunehmender Partikeldichte ist konsistent, da die zeitlichen Anteile, in denen sich kein Partikelbeitrag im Signal befindet mit zunehmender Partikeldichte weniger werden. Es könnten also anteilig weniger Fehlerkennungen im Rauschen stattfinden. 28 zeigt die Fehlerkennungsrate von Ausführungsbeispielen im Rauschen als Funktion des mittleren Partikelabstandes in Einheiten der Maskenlänge.
  • Eine wichtige Erkenntnis in Ausführungsbeispielen könnte sich nun beim beispielhaften Vergleich der Nebenmaximahäufigkeit (29) zeigen. Das konventionelle LABS26(16,#2) 2906 schneidet in beispielhaften Ausführungsbeispielen am schlechtesten ab. B13Man26(13) 2908 und Q1DI26(14,#1) 2902 sind besser. Für beide Codes ergibt sich ein besseres Q2DPs als für den LABS26(16,#2) 2906, obwohl diese Codes nicht auf Q2DPs optimiert sind. Bezüglich der Nebenmaximalhäufigkeit verhält sich Q2DPs26(15,#2) 2904 besser, der weitaus die wenigsten Nebenmaxima in beispielhaften Ausführungsbeispielen erkennt. 29 zeigt die Fehlerkennungsrate aufgrund von Nebenmaxima als Funktion des mittleren Partikelabstandes in Einheiten der Maskenlänge.
  • In einer Strukturgrößennotation für Codes wird gezählt, wie oft ein Symbol eines Codes in direkter Folge wiederholt wird. Am Beispiel des Barker Code 13 bedeutet dies: [ 1 ; 1; 1; 1; 1; 0; 0; 1; 1; 0; 1; 0; 1 ] [ 5 ; 2; 2; 1; 1; 1; 1 ]
    Figure DE102015110316B4_0029
  • Die Strukturgrößennotation ist kürzer. Der Mittelwert der Strukturgrößen wird mittlere Strukturgröße genannt. 30 vergleicht die mittlere Strukturgröße für die besten 2D-Codes 3004a und 3006a aus Ausführungsbeispielen mit den LABS 3002a. Die mittlere Strukturgröße liegt für konventionelle LABS 3002a in einem Bereich nahe 2. Für zunehmendes L (L>40, nicht dargestellt) könnte dieser Zusammenhang noch schärfer werden. Codes mit dem Besten Q2DP 3004a und Q2DPs 3006a könnten sich beispielsweise bevorzugt bei 1,5 bis 1,6 finden. Dies könnte als Kriterium genutzt werden, wonach ein Berechnungsmodul, etwa das Berechnungsmodul 24, bei hohen Codelängen nicht für alle Codes die aufwändig zu ermittelnden Werte Q2DP und Q2DPs berechnen könnte, die Mehrzahl von möglichen Sequenzen einschränken könnte, oder die Teilmenge der möglichen Sequenzen einschränken könnte. 30 zeigte einen Zusammenhang zwischen mittlerer Strukturgröße und Codelänge verglichen für die besten Codes nach den Kriterien Q2DP 3004a und Q2DPs 3006a sowie für LABS 3002a.
  • Im Fall von experimentellen Aufbauten (z.B. Array-Detektoren), in denen die mehrere Modulationsfolgen auf ein Ereignis aufgeprägt werden könnten könnte es sich anbieten, Sequenzen mit verschiedenen Eigenschaften zu nutzen um die Erkennung zu verbessern. So könnte zum Beispiel ein Partikel mit einem optimalen 2D LABS erkannt werden und anschließend mit einer streng periodischen Sequenz (etwa 010101010101010 ..) die Geschwindigkeit ermittelt werden.
  • 31 zeigt ein Ausführungsbeispiel einer Vorrichtung 10 für ein Zytometer 100. Die Vorrichtung 10 ist ausgebildet zur Bereitstellung von Information über ein oder mehrere Zellen in einem Medium in einem Kanal 110. Die Vorrichtung umfasst ein Sensormodul 12, das ausgebildet ist, um ein gefiltertes Signal zu erfassen. Das gefilterte Signal wird von den durch den Kanal 110 strömenden Zellen beeinflusst.
  • In Ausführungsbeispielen könnte das Sensormodul 12 beispielsweise einem lichtempfindlichen Sensor, einem Photodetektor, einem Photosensor, einem Strahlungsdetektor oder einem Strahlungssensor entsprechen. In zumindest manchen Ausführungsbeispielen könnte das Sensormodul 12 ausgebildet sein, Lichtsignale zu erfassen. Alternativ oder zusätzlich könnte das Sensormodul 12 ausgebildet sein, Signale aus elektromagnetische Wellen zu erfassen. In zumindest manchen Ausführungsbeispielen könnte das Sensormodul 12 ein oder mehrere Sensoreinheiten umfassen.
  • Der Kanal 110 könnte beispielsweise einem Flüssigkeitskanal oder einem Flusszytometriekanal entsprechen. Der Kanal 110 könnte beispielsweise ausgebildet sein, das Medium, welches beispielsweise einer Suspension oder einem Fluid entsprechen könnte, welches die ein oder mehreren Zellen umfasst, entlang des Sensormoduls 12 zu transportieren oder zu leiten. Das Medium könnte beispielsweise mittels eines Pumpenmoduls durch den Kanal transportiert werden. In zumindest manchen Ausführungsbeispielen könnte das Zytometer 100 einem Durchflusszytometer und/oder einem Fluoreszenz-basierten Durchflusszytometer entsprechen. In einigen Ausführungsbeispielen könnte das Zytometer 100 beispielsweise einem Zytometer 100 mit Sortierfunktionalität entsprechen, beispielsweise basierend auf Fluoreszenz-aktivierter Zellen-Sortierung (auch engl. Fluorescense-Activated Cell Sorting, FACS).
  • Die Vorrichtung 10 umfasst ferner ein Filter 14. Das Filter 14 ist ausgebildet, um eine vorgegebene Sequenz räumlich abzubilden. Das Filter 14 ist ausgebildet, um basierend auf der Sequenz und einem von den durch den Kanal strömenden Zellen beeinflussten Signal das gefilterte Signal bereitzustellen. In zumindest manchen Ausführungsbeispielen könnte das Filter 14 einer Filtermaske entsprechen. Die Filtermaske könnte beispielsweise ausgebildet sein, das von den durch den Kanal strömenden Zellen beeinflussten Signal räumlich abzuschatten um das gefilterte Signal bereitzustellen. Das Signal könnte beispielsweise auf den ein oder mehreren Zellen basieren, sie beispielsweise von einem Stimulationsmodul angeregt werden, beispielsweise mittels Fluoreszenz basierend auf einem Laser. Die Filtermaske könnte ausgebildet sein, das durch die Fluoreszenz emittierte Licht räumlich basierend auf der Sequenz abzuschatten um das gefilterte Signal bereitzustellen. Alternativ oder zusätzlich könnte das Filter 14 das Stimulationsmodul abschatten, so dass beispielsweise die ein oder mehreren Zellen basierend auf der Sequenz angeregt werden und/oder das Filter könnte ein oder mehrere Stimulationsquellen umfassen, welche ausgebildet sind, basierend auf der Sequenz die ein oder mehreren Zellen anzuregen, um das gefilterte Signal bereitzustellen.
  • Die Vorrichtung 10 umfasst ferner ein Kontrollmodul 16, das ausgebildet ist, um die Information über die ein oder mehreren Zellen basierend auf einer Korrelationsanalyse des erfassten gefilterten Signals und zumindest einer zeitlichen Skalierung der Sequenz bereitzustellen. Die Sequenz beschreibt zeitlich aufeinanderfolgende Signalzustände. Die Folge der Signalzustände innerhalb der Sequenz ist derart gewählt, dass ein Nebenmaximum in einer Korrelationsfunktion der Sequenz mit der zumindest einen zeitlichen Skalierung der Sequenz gegenüber einem durch unterschiedliche Anordnungen der Signalzustände in der Sequenz maximal erzielbaren Nebenmaximum in der Korrelationsfunktion reduziert.
  • In Ausführungsbeispielen kann das Kontrollmodul 16 einem beliebigen Controller oder Prozessor oder einer programmierbaren Hardwarekomponente entsprechen. Beispielsweise kann das Kontrollmodul 16 auch als Software realisiert sein, die für eine entsprechende Hardwarekomponente programmiert ist. Insofern kann das Kontrollmodul 16 als programmierbare Hardware mit entsprechend angepasster Software implementiert sein. Dabei können beliebige Prozessoren, wie Digitale Signalprozessoren (DSPs) zum Einsatz kommen. Ausführungsbeispiele sind dabei nicht auf einen bestimmten Typ von Prozessor eingeschränkt. Es sind beliebige Prozessoren oder auch mehrere Prozessoren zur Implementierung des Kontrollmoduls 16 denkbar.
  • In zumindest manchen Ausführungsbeispielen könnte eine Folge der Signalzustände innerhalb der zumindest einen Sequenz derart gewählt sein, dass eine Summe, die auf Nebenmaxima der Korrelationsfunktion der Sequenz mit der zumindest einen zeitlichen Skalierung der Sequenz basiert, gegenüber einer durch unterschiedliche Anordnungen der Signalzustände in der Sequenz maximal erzielbaren Summe, die auf aus der unterschiedlichen Anordnung resultierenden Nebenmaxima basiert, reduziert ist.
  • Die Folge der Signalzustände innerhalb der zumindest einen Sequenz könnte derart gewählt sein, dass ein größtes Nebenmaximum der Korrelationsfunktion der Sequenz mit der zumindest einen zeitlichen Skalierung der Sequenz gegenüber einem durch unterschiedliche Anordnungen der Signalzustände in der Sequenz maximal erzielbaren größten Nebenmaximum reduziert ist.
  • In einigen Ausführungsbeispielen könnte die Summe, die auf den Nebenmaxima basiert, unter einem ersten Schwellenwert liegen. Das größte Nebenmaximum könnte beispielsweise unter einem zweiten Schwellenwert liegen.
  • In zumindest einigen Ausführungsbeispielen könnte der erste Schwellenwert einer zweitniedrigsten Summe, die auf Nebenmaxima von Korrelationsfunktionen möglicher Sequenzen basiert, entsprechen. Der zweite Schwellenwert könnte beispielsweise einem größten Nebenmaximum einer Sequenz mit einem zweitniedrigsten größten Nebenmaximum der Korrelationsfunktionen entspricht.
  • In zumindest manchen Ausführungsbeispielen könnte die Sequenz einer unipolaren Sequenz entsprechen. In zumindest manchen Ausführungsbeispielen könnte das Filter 14 balanciert sein.
  • Mehr Details und Aspekte der Vorrichtung 10 (z.B. Sequenz, Korrelationsfunktion, zeitliche Skalierung, zeitlich aufeinanderfolgende Signalzustände, erster Schwellenwert, zweiter Schwellenwert, unipolare Sequenz) werden in Verbindung mit dem Konzept oder Beispielen genannt, die vorher (z.B. 1 bis 30) beschrieben wurden. Die Vorrichtung 10 kann ein oder mehrere zusätzliche optionale Merkmale umfassen, die ein oder mehreren Aspekten des vorgeschlagenen Konzepts oder der beschriebenen Beispiele entsprechen, wie sie vorher oder nachher beschrieben wurden.
  • 32 zeigt ein Blockdiagramm eines Ausführungsbeispiels eines Zytometers 100. Das Zytometer umfasst eine Vorrichtung 10 sowie einen Kanal 110. Die Vorrichtung 10 ist ausgebildet, Information über ein oder mehrere Zellen in einem Medium in dem Kanal 110 dem Zytometer 100 bereitzustellen.
  • Mehr Details und Aspekte des Zytometers 100 (z.B. Zytometer 100, Vorrichtung 10, Kanal 110) werden in Verbindung mit dem Konzept oder Beispielen genannt, die vorher (z.B. 1 bis 31) beschrieben wurden. Das Zytometer 100 kann ein oder mehrere zusätzliche optionale Merkmale umfassen, die ein oder mehreren Aspekten des vorgeschlagenen Konzepts oder der beschriebenen Beispiele entsprechen, wie sie vorher oder nachher beschrieben wurden.
  • 33 zeigt ein Flussdiagramm eines Verfahrens zum Bereitstellen von Information über zumindest eine Sequenz. Die zumindest eine Sequenz beschreibt zeitlich aufeinanderfolgende Signalzustände. Das Verfahren umfasst Erhalten 32 von Information über eine Anzahl der Signalzustände und Generieren 34 einer Mehrzahl von möglichen Sequenzen, basierend auf der Information über die Anzahl der Signalzustände. Das Verfahren umfasst ferner Berechnen 36 von Korrelationsfunktionen zwischen einer Sequenz und zumindest einer zeitlichen Skalierung der Sequenz für zumindest eine Teilmenge der möglichen Sequenzen. Eine Korrelationsfunktion umfasst ein Hauptmaximum und ein oder mehrere Nebenmaxima. Das Verfahren umfasst ferner Bestimmen 38 der zumindest eine Sequenz basierend auf den Korrelationsfunktionen. Die Folge der Signalzustände innerhalb der zumindest einen Sequenz ist derart gewählt, dass ein Nebenmaximum in einer Korrelationsfunktion der Sequenz mit der zumindest einen zeitlichen Skalierung der Sequenz gegenüber einem durch unterschiedliche Anordnungen der Signalzustände in der Sequenz maximal in einer Korrelationsfunktion erzielbaren Nebenmaximum reduziert ist. Das Verfahren umfasst ferner Bestimmen 40 der Information über die zumindest eine Sequenz basierend auf der zumindest einen Sequenz, und Bereitstellen 42 der Information über die zumindest eine Sequenz.
  • 34 zeigt ein Flussdiagramm eines Verfahrens zur Bereitstellung von Information über ein oder mehrere Zellen in einem Medium in einem Kanal 110. Die Vorrichtung umfasst Erfassen 52 eines gefilterten Signals. Das gefilterte Signal wird von den durch den Kanal 110 strömenden Zellen beeinflusst. Das Verfahren umfasst ferner Bereitstellen 54 des gefilterten Signals basierend auf einer vorgegebenen Sequenz, einem von den durch den Kanal strömenden Zellen beeinflussten Signal und einer räumlichen Abbildung der Sequenz. Das Verfahren umfasst ferner Bereitstellen 56 der Information über die ein oder mehreren Zellen basierend auf einer Korrelationsanalyse des erfassten gefilterten Signals und zumindest einer zeitlichen Skalierung der Sequenz. Die Sequenz beschreibt zeitlich aufeinanderfolgende Signalzustände. Die Folge der Signalzustände innerhalb der Sequenz ist derart gewählt, dass ein Nebenmaximum in einer Korrelationsfunktion der Sequenz mit der zumindest einen zeitlichen Skalierung der Sequenz gegenüber einem durch unterschiedliche Anordnungen der Signalzustände in der Sequenz maximal erzielbaren Nebenmaximum in der Korrelationsfunktion reduziert ist.
  • Ein weiteres Ausführungsbeispiel ist ein Computerprogramm zur Durchführung zumindest eines der oben beschriebenen Verfahren, wenn das Computerprogramm auf einem Computer, einem Prozessor oder einer programmierbaren Hardwarekomponente abläuft. Ein weiteres Ausführungsbeispiele ist auch ein digitales Speichermedium, das maschinen- oder computerlesbar ist, und das elektronisch lesbare Steuersignale aufweist, die mit einer programmierbaren Hardwarekomponente so zusammenwirken können, dass eines der oben beschriebenen Verfahren ausgeführt wird.
  • Die in der vorstehenden Beschreibung, den nachfolgenden Ansprüchen und den beigefügten Figuren offenbarten Merkmale können sowohl einzeln wie auch in beliebiger Kombination für die Verwirklichung eines Ausführungsbeispiels in ihren verschiedenen Ausgestaltungen von Bedeutung sein und implementiert werden.
  • Obwohl manche Aspekte im Zusammenhang mit einer Vorrichtung beschrieben wurden, versteht es sich, dass diese Aspekte auch eine Beschreibung des entsprechenden Verfahrens darstellen, sodass ein Block oder ein Bauelement einer Vorrichtung auch als ein entsprechender Verfahrensschritt oder als ein Merkmal eines Verfahrensschrittes zu verstehen ist. Analog dazu stellen Aspekte, die im Zusammenhang mit einem oder als ein Verfahrensschritt beschrieben wurden, auch eine Beschreibung eines entsprechenden Blocks oder Details oder Merkmals einer entsprechenden Vorrichtung dar.
  • Je nach bestimmten Implementierungsanforderungen können Ausführungsbeispiele der Erfindung in Hardware oder in Software implementiert sein. Die Implementierung kann unter Verwendung eines digitalen Speichermediums, beispielsweise einer Floppy-Disk, einer DVD, einer Blu-Ray Disc, einer CD, eines ROM, eines PROM, eines EPROM, eines EEPROM oder eines FLASH-Speichers, einer Festplatte oder eines anderen magnetischen oder optischen Speichers durchgeführt werden, auf dem elektronisch lesbare Steuersignale gespeichert sind, die mit einer programmierbaren Hardwarekomponente derart zusammenwirken können oder zusammenwirken, dass das jeweilige Verfahren durchgeführt wird.
  • Eine programmierbare Hardwarekomponente kann durch einen Prozessor, einen Computerprozessor (CPU = Central Processing Unit), einen Grafikprozessor (GPU = Graphics Processing Unit), einen Computer, ein Computersystem, einen anwendungsspezifischen integrierten Schaltkreis (ASIC = Application-Specific Integrated Circuit), einen integrierten Schaltkreis (IC = Integrated Circuit), ein Ein-Chip-System (SOC = System on Chip), ein programmierbares Logikelement oder ein feldprogrammierbares Gatterarray mit einem Mikroprozessor (FPGA = Field Programmable Gate Array) gebildet sein.
  • Das digitale Speichermedium kann daher maschinen- oder computerlesbar sein. Manche Ausführungsbeispiele umfassen also einen Datenträger, der elektronisch lesbare Steuersignale aufweist, die in der Lage sind, mit einem programmierbaren Computersystem oder einer programmierbare Hardwarekomponente derart zusammenzuwirken, dass eines der hierin beschriebenen Verfahren durchgeführt wird. Ein Ausführungsbeispiel ist somit ein Datenträger (oder ein digitales Speichermedium oder ein computerlesbares Medium), auf dem das Programm zum Durchführen eines der hierin beschriebenen Verfahren aufgezeichnet ist.
  • Allgemein können Ausführungsbeispiele der vorliegenden Erfindung als Programm, Firmware, Computerprogramm oder Computerprogrammprodukt mit einem Programmcode oder als Daten implementiert sein, wobei der Programmcode oder die Daten dahin gehend wirksam ist bzw. sind, eines der Verfahren durchzuführen, wenn das Programm auf einem Prozessor oder einer programmierbaren Hardwarekomponente abläuft. Der Programmcode oder die Daten kann bzw. können beispielsweise auch auf einem maschinenlesbaren Träger oder Datenträger gespeichert sein. Der Programmcode oder die Daten können unter anderem als Quellcode, Maschinencode oder Bytecode sowie als anderer Zwischencode vorliegen.
  • Ein weiteres Ausführungsbeispiel ist ferner ein Datenstrom, eine Signalfolge oder eine Sequenz von Signalen, der bzw. die das Programm zum Durchführen eines der hierin beschriebenen Verfahren darstellt bzw. darstellen. Der Datenstrom, die Signalfolge oder die Sequenz von Signalen kann bzw. können beispielsweise dahin gehend konfiguriert sein, um über eine Datenkommunikationsverbindung, beispielsweise über das Internet oder ein anderes Netzwerk, transferiert zu werden. Ausführungsbeispiele sind so auch Daten repräsentierende Signalfolgen, die für eine Übersendung über ein Netzwerk oder eine Datenkommunikationsverbindung geeignet sind, wobei die Daten das Programm darstellen.
  • Ein Programm gemäß einem Ausführungsbeispiel kann eines der Verfahren während seiner Durchführung beispielsweise dadurch umsetzen, dass dieses Speicherstellen ausliest oder in diese ein Datum oder mehrere Daten hinein schreibt, wodurch gegebenenfalls Schaltvorgänge oder andere Vorgänge in Transistorstrukturen, in Verstärkerstrukturen oder in anderen elektrischen, optischen, magnetischen oder nach einem anderen Funktionsprinzip arbeitenden Bauteile hervorgerufen werden. Entsprechend können durch ein Auslesen einer Speicherstelle Daten, Werte, Sensorwerte oder andere Informationen von einem Programm erfasst, bestimmt oder gemessen werden. Ein Programm kann daher durch ein Auslesen von einer oder mehreren Speicherstellen Größen, Werte, Messgrößen und andere Informationen erfassen, bestimmen oder messen, sowie durch ein Schreiben in eine oder mehrere Speicherstellen eine Aktion bewirken, veranlassen oder durchführen sowie andere Geräte, Maschinen und Komponenten ansteuern.
  • Ausführungsbeispiele können beispielsweise genutzt werden, um verbesserte Codes zu bestimmen, die folgende Vorteile aufweisen können:
    • ■ Signifikant verringerte Fehlererkennungsrate
    • ■ Verbesserte Erkennung im Fall überlappender Signale
    • ■ Verbesserte Sensitivität für Signale geringer Amplitude im Fall der Überlappung mit Signalen höherer Amplitude
  • Die oben beschriebenen Ausführungsbeispiele stellen lediglich eine Veranschaulichung der Prinzipien der vorliegenden Erfindung dar. Es versteht sich, dass Modifikationen und Variationen der hierin beschriebenen Anordnungen und Einzelheiten anderen Fachleuten einleuchten werden. Deshalb ist beabsichtigt, dass die Erfindung lediglich durch den Schutzumfang der nachstehenden Patentansprüche und nicht durch die spezifischen Einzelheiten, die anhand der Beschreibung und der Erläuterung der Ausführungsbeispiele hierin präsentiert wurden, beschränkt sei.
  • Die Forschungsarbeiten, die zu diesen Ergebnissen geführt haben, wurden von der Europäischen Union gefördert.

Claims (20)

  1. Vorrichtung (20) zum Bereitstellen von Information über zumindest eine Sequenz, wobei die zumindest eine Sequenz zeitlich aufeinanderfolgende Signalzustände beschreibt, die Vorrichtung (20) umfassend eine Schnittstelle (22), die ausgebildet ist, um Information über eine Anzahl der Signalzustände zu erhalten; und ein Berechnungsmodul (24), welches ausgebildet ist, um: eine Mehrzahl von möglichen Sequenzen zu generieren, basierend auf der Information über die Anzahl der Signalzustände, für zumindest eine Teilmenge der möglichen Sequenzen Korrelationsfunktionen zwischen einer Sequenz und zumindest einer zeitlichen Skalierung der Sequenz zu berechnen, wobei eine Korrelationsfunktion ein Hauptmaximum und ein oder mehrere Nebenmaxima umfasst, die zumindest eine Sequenz basierend auf den Korrelationsfunktionen zu bestimmen, wobei die Folge der Signalzustände innerhalb der zumindest einen Sequenz derart gewählt ist, dass ein Nebenmaximum in einer Korrelationsfunktion der Sequenz mit der zumindest einen zeitlichen Skalierung der Sequenz gegenüber einem durch unterschiedliche Anordnungen der Signalzustände in der Sequenz maximal in einer Korrelationsfunktion erzielbaren Nebenmaximum reduziert ist, und die Information über die zumindest eine Sequenz basierend auf der zumindest einen Sequenz zu bestimmen und über die Schnittstelle (22) bereitzustellen.
  2. Vorrichtung (20) gemäß Anspruch 1, wobei das Berechnungsmodul (24) ausgebildet ist, um die zumindest eine Sequenz basierend auf den Korrelationsfunktionen zu bestimmen, wobei die Folge der Signalzustände innerhalb der zumindest einen Sequenz derart gewählt ist, dass eine Summe, die auf den Nebenmaxima einer Korrelationsfunktion der Sequenz mit der zumindest einen zeitlichen Skalierung der Sequenz basiert, gegenüber einer durch unterschiedliche Anordnungen der Signalzustände in der Sequenz maximal erzielbaren Summe, die auf aus der unterschiedlichen Anordnung resultierenden Nebenmaxima basiert, reduziert ist, und/oder wobei die Folge der Signalzustände innerhalb der zumindest einen Sequenz derart gewählt ist, dass ein größtes Nebenmaximum in der Korrelationsfunktion der Sequenz mit der zumindest einen zeitlichen Skalierung der Sequenz gegenüber einem durch unterschiedliche Anordnungen der Signalzustände in der Sequenz maximal erzielbaren größten Nebenmaximum reduziert ist.
  3. Vorrichtung gemäß Anspruch 2, wobei das Berechnungsmodul (24) ausgebildet ist, um die zumindest eine Sequenz so zu bestimmen, dass die Summe, die auf den Nebenmaxima basiert, unter einem ersten Schwellenwert liegt, und/oder wobei das Berechnungsmodul (24) ausgebildet ist, um die zumindest eine Sequenz so zu bestimmen, dass das größte Nebenmaximum unter einem zweiten Schwellenwert liegt.
  4. Vorrichtung gemäß Anspruch 3, wobei der erste Schwellenwert einer zweitniedrigsten Summe, die auf den Nebenmaxima der Korrelationsfunktionen basiert, entspricht, und/oder wobei der zweite Schwellenwert einem größten Nebenmaximum einer Sequenz mit einem zweitniedrigsten größten Nebenmaximum der Korrelationsfunktionen entspricht.
  5. Vorrichtung (20) gemäß einem der vorangehenden Ansprüche, wobei das Berechnungsmodul (24) ausgebildet ist, um bereinigte Korrelationsfunktionen basierend auf den Korrelationsfunktionen zu bestimmen, wobei das Berechnungsmodul (24) ausgebildet ist, um einen Beitrag von Hauptmaxima in den bereinigten Korrelationsfunktionen zu reduzieren, und wobei das Berechnungsmodul (24) ausgebildet ist, um die Information über die zumindest eine Sequenz basierend auf den bereinigten Korrelationsfunktionen zu bestimmen.
  6. Vorrichtung (20) gemäß Anspruch 5, wobei das Berechnungsmodul (24) ausgebildet ist, um unter den Hauptmaxima der Korrelationsfunktionen ein größtes Hauptmaximum zu bestimmen, und wobei das Berechnungsmodul (24) ausgebildet ist, um über die zeitliche Ausdehnung des größten Hauptmaximums den Beitrag der Hauptmaxima in den bereinigten Korrelationsfunktionen zu reduzieren.
  7. Vorrichtung (20) gemäß Anspruch 6, wobei das Berechnungsmodul (24) ausgebildet ist, um basierend auf einer zeitlichen Position des größten Hauptmaximums und basierend auf der Korrelationsfunktion der zeitlich skalierten Sequenz die das größte Hauptmaximum umfasst, die Beiträge der Hauptmaxima in den bereinigten Korrelationsfunktionen der zeitlich skalierten Sequenz und benachbarten zeitlichen Skalierungen der Sequenz zu bestimmen und zu reduzieren.
  8. Vorrichtung (20) gemäß einem der vorhergehenden Ansprüche, wobei eine Sequenz eine mittlere Strukturgröße hat, wobei eine Strukturgröße auf einer Anzahl von gleichen aufeinanderfolgenden Signalzuständen basiert, und wobei das Berechnungsmodul (24) ausgebildet ist, um die Mehrzahl von möglichen Sequenzen basierend auf einem Zielbereich für die mittlere Strukturgröße zu generieren und/oder wobei das Berechnungsmodul (24) ausgebildet ist, die Teilmenge der möglichen Sequenzen basierend auf der mittleren Strukturgröße und/oder dem Zielbereich für die mittlere Strukturgröße zu bestimmen.
  9. Vorrichtung (20) gemäß Anspruch 8, wobei die mittlere Strukturgröße einem Durchschnitt der Strukturgrößen einer Sequenz entspricht, und wobei der Zielbereich für die mittlere Strukturgröße zwischen 1,3 und 1,8 liegt.
  10. Vorrichtung (20) gemäß einem der vorhergehenden Ansprüche, wobei das Berechnungsmodul (24) ausgebildet ist, um die Mehrzahl von möglichen Sequenzen als unipolare Sequenzen zu bestimmen.
  11. Vorrichtung (20) gemäß einem der vorhergehenden Ansprüche, wobei eine Sequenz einer Binärsequenz entspricht und wobei die Binärsequenz zeitlich äquidistante aufeinanderfolgende Signalzustände beschreibt.
  12. Vorrichtung (10) für ein Zytometer (100), zur Bereitstellung von Information über ein oder mehrere Zellen in einem Medium in einem Kanal (110), die Vorrichtung umfassend ein Sensormodul (12), das ausgebildet ist, um ein gefiltertes Signal zu erfassen, wobei das gefilterte Signal von den durch den Kanal (110) strömenden Zellen beeinflusst wird, ein Filter (14), das ausgebildet ist, um eine vorgegebene Sequenz räumlich abzubilden, und ausgebildet ist, um basierend auf der Sequenz und einem von den durch den Kanal strömenden Zellen beeinflussten Signal das gefilterte Signal bereitzustellen, und ein Kontrollmodul (16) das ausgebildet ist, um die Information über die ein oder mehreren Zellen basierend auf einer Korrelationsanalyse des erfassten gefilterten Signals und zumindest einer zeitlichen Skalierung der Sequenz bereitzustellen, wobei die Sequenz zeitlich aufeinanderfolgende Signalzustände beschreibt, und wobei die Folge der Signalzustände innerhalb der Sequenz derart gewählt ist, dass ein Nebenmaximum in einer Korrelationsfunktion der Sequenz mit der zumindest einen zeitlichen Skalierung der Sequenz gegenüber einem durch unterschiedliche Anordnungen der Signalzustände in der Sequenz maximal erzielbaren Nebenmaximum in der Korrelationsfunktion reduziert ist.
  13. Vorrichtung (10) gemäß Anspruch 12, wobei eine Folge der Signalzustände innerhalb der zumindest einen Sequenz derart gewählt ist, dass eine Summe, die auf Nebenmaxima der Korrelationsfunktion der Sequenz mit der zumindest einen zeitlichen Skalierung der Sequenz basiert, gegenüber einer durch unterschiedliche Anordnungen der Signalzustände in der Sequenz maximal erzielbaren Summe, die auf aus der unterschiedlichen Anordnung resultierenden Nebenmaxima basiert, reduziert ist, und/oder wobei die Folge der Signalzustände innerhalb der zumindest einen Sequenz derart gewählt ist, dass ein größtes Nebenmaximum der Korrelationsfunktion der Sequenz mit der zumindest einen zeitlichen Skalierung der Sequenz gegenüber einem durch unterschiedliche Anordnungen der Signalzustände in der Sequenz maximal erzielbaren größten Nebenmaximum reduziert ist.
  14. Vorrichtung (10) gemäß Anspruch 13, wobei die Summe, die auf den Nebenmaxima basiert, unter einem ersten Schwellenwert liegt, und/oder wobei das größte Nebenmaximum unter einem zweiten Schwellenwert liegt.
  15. Vorrichtung gemäß Anspruch 14, wobei der erste Schwellenwert einer zweitniedrigsten Summe, die auf Nebenmaxima von Korrelationsfunktionen möglicher Sequenzen basiert, entspricht, und/oder wobei der zweite Schwellenwert einem größten Nebenmaximum einer Sequenz mit einem zweitniedrigsten größten Nebenmaximum der Korrelationsfunktionen entspricht.
  16. Vorrichtung (10) gemäß einem der Ansprüche 12 bis 15, wobei die Sequenz einer unipolaren Sequenz entspricht.
  17. Zytometer (100) umfassend eine Vorrichtung (10) gemäß einem der Ansprüche 12 bis 16 sowie einen Kanal (110), wobei die Vorrichtung (10) ausgebildet ist, Information über ein oder mehrere Zellen in einem Medium in dem Kanal (110) dem Zytometer (100) bereitzustellen.
  18. Verfahren zum Bereitstellen von Information über zumindest eine Sequenz, wobei die zumindest eine Sequenz zeitlich aufeinanderfolgende Signalzustände beschreibt, das Verfahren umfassend Erhalten (32) von Information über eine Anzahl der Signalzustände, Generieren (34) einer Mehrzahl von möglichen Sequenzen, basierend auf der Information über die Anzahl der Signalzustände, Berechnen (36) von Korrelationsfunktionen zwischen einer Sequenz und zumindest einer zeitlichen Skalierung der Sequenz für zumindest eine Teilmenge der möglichen Sequenzen, wobei eine Korrelationsfunktion ein Hauptmaximum und ein oder mehrere Nebenmaxima umfasst, Bestimmen (38) der zumindest eine Sequenz basierend auf den Korrelationsfunktionen, wobei die Folge der Signalzustände innerhalb der zumindest einen Sequenz derart gewählt ist, dass ein Nebenmaximum in einer Korrelationsfunktion der Sequenz mit der zumindest einen zeitlichen Skalierung der Sequenz gegenüber einem durch unterschiedliche Anordnungen der Signalzustände in der Sequenz maximal in einer Korrelationsfunktion erzielbaren Nebenmaximum reduziert ist, und Bestimmen (40) der Information über die zumindest eine Sequenz basierend auf der zumindest einen Sequenz, und Bereitstellen (42) der Information über die zumindest eine Sequenz.
  19. Verfahren zur Bereitstellung von Information über ein oder mehrere Zellen in einem Medium in einem Kanal (110), das Verfahren umfassend Erfassen (52) eines gefilterten Signals, wobei das gefilterte Signal von den durch den Kanal (110) strömenden Zellen beeinflusst wird, Bereitstellen (54) des gefilterten Signals basierend auf einer vorgegebenen Sequenz, einem von den durch den Kanal strömenden Zellen beeinflussten Signal und einer räumlichen Abbildung der Sequenz, und Bereitstellen (56) der Information über die ein oder mehreren Zellen basierend auf einer Korrelationsanalyse des erfassten gefilterten Signals und zumindest einer zeitlichen Skalierung der Sequenz, wobei die Sequenz zeitlich aufeinanderfolgende Signalzustände beschreibt, und wobei die Folge der Signalzustände innerhalb der Sequenz derart gewählt ist, dass ein Nebenmaximum in einer Korrelationsfunktion der Sequenz mit der zumindest einen zeitlichen Skalierung der Sequenz gegenüber einem durch unterschiedliche Anordnungen der Signalzustände in der Sequenz maximal erzielbaren Nebenmaximum in der Korrelationsfunktion reduziert ist.
  20. Computerprogramm mit einem Programmcode, der ein Durchführen eines der Verfahren der Ansprüche 18 und 19 bewirkt, wenn der Programmcode auf einem Prozessor oder einer programmierbaren Hardwarekomponente abläuft.
DE102015110316.8A 2015-06-26 2015-06-26 Vorrichtungen, Zytometer, Verfahren und Computerprogramm zum Bereitstellen von Information über zumindest eine Sequenz Active DE102015110316B4 (de)

Priority Applications (4)

Application Number Priority Date Filing Date Title
DE102015110316.8A DE102015110316B4 (de) 2015-06-26 2015-06-26 Vorrichtungen, Zytometer, Verfahren und Computerprogramm zum Bereitstellen von Information über zumindest eine Sequenz
CN201680036244.2A CN107771278B (zh) 2015-06-26 2016-06-23 提供关于序列的信息的装置、细胞仪、方法和存储介质
PCT/EP2016/064492 WO2016207262A1 (de) 2015-06-26 2016-06-23 Vorrichtungen, zytometer, verfahren und computerprogramm zum bereitstellen von information über zumindest eine sequenz
US15/738,162 US10551295B2 (en) 2015-06-26 2016-06-23 Devices, cytometers, methods and computer program for providing information on at least one sequence

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
DE102015110316.8A DE102015110316B4 (de) 2015-06-26 2015-06-26 Vorrichtungen, Zytometer, Verfahren und Computerprogramm zum Bereitstellen von Information über zumindest eine Sequenz

Publications (2)

Publication Number Publication Date
DE102015110316A1 DE102015110316A1 (de) 2016-12-29
DE102015110316B4 true DE102015110316B4 (de) 2021-09-09

Family

ID=56194494

Family Applications (1)

Application Number Title Priority Date Filing Date
DE102015110316.8A Active DE102015110316B4 (de) 2015-06-26 2015-06-26 Vorrichtungen, Zytometer, Verfahren und Computerprogramm zum Bereitstellen von Information über zumindest eine Sequenz

Country Status (4)

Country Link
US (1) US10551295B2 (de)
CN (1) CN107771278B (de)
DE (1) DE102015110316B4 (de)
WO (1) WO2016207262A1 (de)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR102664396B1 (ko) * 2019-01-15 2024-05-08 삼성전자주식회사 라이다 장치 및 그 동작 방법
EP4017635A1 (de) * 2019-08-21 2022-06-29 Life Technologies Corporation Vorrichtungen mit mehrspuriger durchflusszelle

Family Cites Families (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20040206771A1 (en) * 2003-04-18 2004-10-21 Eric Junkel Water toy with two port elastic fluid bladder
KR100592547B1 (ko) * 2003-06-03 2006-06-23 서울시립대학교 산학협력단 스트리밍을 위한 패킷 스케줄링 방법
US8315307B2 (en) * 2004-04-07 2012-11-20 Qualcomm Incorporated Method and apparatus for frame prediction in hybrid video compression to enable temporal scalability
EP1817374A4 (de) * 2004-12-03 2008-03-26 Lg Chemical Ltd Rohrbehälter mit sperreigenschaft
JP4745030B2 (ja) * 2005-11-15 2011-08-10 シスメックス株式会社 血液分析装置
JP5511149B2 (ja) * 2008-03-28 2014-06-04 シスメックス株式会社 生体試料分析装置及び生体試料分析方法
EP2211164A1 (de) 2009-01-27 2010-07-28 Koninklijke Philips Electronics N.V. Fingerelektroden zur mikrofluidischen Einzelteilchenanalyse
US9423335B2 (en) * 2011-07-21 2016-08-23 Invitrox, Inc. Instrument and method for optical particle sensing
CN103084229B (zh) * 2012-01-16 2015-09-23 中国科学院深圳先进技术研究院 微流控芯片、血细胞分析系统以及血细胞分析方法
CN104641217B (zh) * 2012-05-30 2019-02-22 艾瑞斯国际有限公司 流式细胞仪
CN102867225A (zh) * 2012-09-21 2013-01-09 四川大学 一种电力小时负荷的混沌线性回归模型安全监控预测方法
CN104598998A (zh) * 2015-02-04 2015-05-06 国家电网公司 基于经济增长指标的能源需求量预测方法

Non-Patent Citations (27)

* Cited by examiner, † Cited by third party
Title
C. Sommer et al.: The equilibrium velocity of spherical particles in rectangular microfluidic channels for size measurement, Lab on a Chip, Vol. 14, 2014, No. 13, S. 2319-2326
C. Sommer: Die Größenabhängigkeit der Gleichgewichtsgeschwindigkeit von Partikeln beim Transport in Mikrokanälen, Dissertation der Technischen Universität Darmstadt, D17, Juli 2014
COHEN, Marvin N. ; FOX, Marshall R. ; BADEN, J. Michael: Minimum peak sidelobe pulse compression codes. In: IEEE International Conference on Radar, 7-10 May 1990, S. 633-638. – URL: https://ieeexplore.ieee.org/stamp/stamp.jsp?tp=&arnumber=201102 [abgerufen am 24.3.21]
FORSTER, Roger: Manchester encoding: opposing definitions resolved. In: Engineering Science and Education Journal, Vol. 9, 2000, No. 6, S. 278-280. - ISSN 0963-7346
G. E. Coxson and J. C. Russo, Efficient exhaustive search for optimal-peak-sidelobe binary codes, IEEE Transactions on Aerospace and Electronic Systems 41, 302 (2005)
H. D. Lüke et al: Binary and quadriphase sequences with optimal autocorrelation properties: a survey, IEEE Transactions on Information Theory, Vol. 49, 2003, No. 12, S. 3271-3282
J. Lindner: Binary sequences up to length 40 with best possible autocorrelation function, Electronics Letters, Vol. 11, 1975, No. 21, S. 507
KETTLITZ, Siegfried W. [et al.]: Particle detection form spatially modulated fluorescence signals. In: Proceedings of SPIE, Vol. 9129, 2014, S. 1-7. – doi: 10.1117/12.2052169
KETTLITZ, Siegfried W. [et al.]: Sensitivity Improvement in Fluorescence-Based Particle Detection. In: Cytometry A, Vol. 85A, 2014, S. 746-755. – doi: 10.1002/cyto.a.22499
KIESEL, Peter [et al.]: Spatially modulated fluorescence emission from moving particles. In. Appl. Phys. Lett., Vol. 94, 2009, S. 411071-411073.
Levanon et al. 2004, Turyn et al. 1961, Mertens et al. 1996, Coxson et al., 2005, Borewein et al. 2000, Nunn et al. 2008, Shanmugam et al. 2008
LINDNER, J.: Binary sequences up to length 40 with best possible autocorrelation function. In: Electronics Letters, Vol. 11, 1975, No. 21, S. 507. – ISSN 0013-5194
LÜKE, H. Dieter ; SCHOTTEN, Hans D. ; HADINEJAD-MAHRAM, Hafez: Binary and quadriphase sequences with optimal autocorrelation properties: a survey. In: IEEE Transactions on Information Theory, Vol. 49, 2003, No. 12, S. 3271-3282. – ISSN 0018-9448
M. N. Cohen et al: Minimum peak sidelobe pulse compression codes, IEEE International Conference on Radar, 7-10 May 1990, S. 633-638
N. Levanon and E. Mozeson, Radar Signals, Wiley, New York, 2004
P. Borewein, R. Ferguson, and J. Knauer, The merit factor problem, 2000; C. J. Nunn and G. E. Coxson, Best-known autocorrelation peak sidelobe levels for binary codes oflength 71 to 105, IEEE Transactions on Aerospace and Electronic Systems 44, 392 (2008)
P. Kiesel, M. Bassler, M. Beck and N. Johnson, Appl. Phys. Lett., 2009, 94, p. 41107
QUINT, Stephan: Durchflusszytometrie nach dem Prinzip der räumlich modulierten Fluoreszenz. Dissertation der Universität Darmstadt, 2014, Abstract
R. Foster: Manchester encoding: opposing definitions resolved, Engineering Science and Education Journal, Vol. 9, 2000, No. 6, S. 278-280. ISSN 0963-7346
R. Turyn and S. Storer, On binary sequences, Proceedings of the American Mathematical Society 12, 394 (1961)
S. K. Shanmugam, C. Mongredien, J. Nielsen, and G. Lachapelle, Design of short synchronization codes for use in future GNSS system, International Journal of Navigation and Observation 2008 (2008), Article ID 246703
S. Mertens, Exhaustive search for low-autocorrelation binary sequences, Journal of Physics A: Mathematical and General 29, L473 (1996)
S. Quint: Durchflusszytometrie nach dem Prinzip der räumlich modulierten Fluoreszenz. Dissertation der Universität Darmstadt, 2014, Abstract
S. W. Kettlitz et al.: Particle detection form spatially modulated fluorescence Signals, Proceedings of SPIE, Vol. 9129, 2014, S. 1-7. - doi: 10.1117/12.2052169
S. W. Kettlitz et al.: Sensitivity Improvement in Fluorescence-Based Particle Detection, Cytometry A, Vol. 85A, 2014, S. 746-755, - doi: 10.1002/cyto.a.22499
SOMMER, Christian [et al.]: The equilibrium velocity of spherical particles in rectangular microfluidic channels for size measurement. In: Lab on a Chip, Vol. 14, 2014, No. 13, S. 2319-2326. – ISSN 0013-5194
SOMMER, Christian: Die Größenabhängigkeit der Gleichgewichtsgeschwindigkeit von Partikeln beim Transport in Mikrokanälen. Dissertation der Technischen Universität Darmstadt, D17, Juli 2014

Also Published As

Publication number Publication date
DE102015110316A1 (de) 2016-12-29
WO2016207262A1 (de) 2016-12-29
US10551295B2 (en) 2020-02-04
US20180172575A1 (en) 2018-06-21
CN107771278B (zh) 2020-11-20
CN107771278A (zh) 2018-03-06

Similar Documents

Publication Publication Date Title
EP3410154B1 (de) Verfahren und sensorsystem zur vermessung der eigenschaften einer übertragungsstrecke eines messsystems zwischen sender und empfänger
EP3341891A1 (de) Aktives maschinelles lernen zum trainieren eines ereignisklassifikators
Pnevmatikakis et al. Sparse nonnegative deconvolution for compressive calcium imaging: algorithms and phase transitions
DE102018124117A1 (de) Verfahren zur Verfolgung eines Objekts
DE102020203796A1 (de) Verfahren und Vorrichtung zum Bestimmen einer Messinformation und LiDAR-Gerät
EP2584493A2 (de) Verfahren zur Unterscheidung zwischen einem realen Gesicht und einer zweidimensionalen Abbildung des Gesichts in einem biometrischen Erfassungsprozess
DE112016005278T5 (de) Probabilistischer single-flux-quanten-digitalumsetzer
EP2400458B1 (de) Verfahren und Vorrichtung zur Segmentierung von biologischen Zellen in einer Aufnahme
DE102011017514A1 (de) Vorrichtung und Verfahren zum Verarbeiten eines Doppler-Signals
CN105913081A (zh) 基于改进的PCAnet的SAR图像分类方法
DE102015110316B4 (de) Vorrichtungen, Zytometer, Verfahren und Computerprogramm zum Bereitstellen von Information über zumindest eine Sequenz
DE19648834A1 (de) Verfahren zum Betrieb und zur Auswertung von Signalen einer Wirbelstromsonde und Vorrichtung zur Ausübung des Verfahrens
Ghassami et al. Interaction information for causal inference: The case of directed triangle
DE102018206108A1 (de) Generieren von Validierungsdaten mit generativen kontradiktorischen Netzwerken
EP2549408A1 (de) Verfahren zur Detektion und Klassifikation von Objekten
EP3430565B1 (de) Verfahren zur untersuchung verteilter objekte unter segmentierung eines übersichtbildes
Saha et al. A classification-based approach to estimate the number of resting fMRI dynamic functional connectivity states
DE102008008571B4 (de) Hintergrundbilderzeugungsvorrichtung
Huang et al. On the access fairness of the DQDB MAN protocol
DE102015110359B4 (de) Verfahren zur Detektion von strahlungsemittierenden Partikeln
DE60221855T2 (de) Schätzung von zielsignalen durch einen dynamischen fluxgate-sensor
DE112020006661T5 (de) Merkmalsmenge-extrahierungseinrichtung, zeitsequentielle inferenzvorrichtung,zeitsequentielles lernsystem, zeitsequentielles merkmalsemengeextrahierungsverfahren,zeitsequentielles inferenzverfahren und zeitsequentielles lernverfahren
EP2550540B1 (de) Verfahren zum erkennen von magnetisch gekennzeichneten objekten sowie entsprechende vorrichtung
Bagalkote et al. Multiresolution Analysis and Implementation of Grape Species Classification Using Neural Network
Tiwari et al. Identification of topological measures of visibility graphs for analyzing transitions in complex time series

Legal Events

Date Code Title Description
R012 Request for examination validly filed
R016 Response to examination communication
R016 Response to examination communication
R018 Grant decision by examination section/examining division
R020 Patent grant now final