A method and device for identifying optical artifacts in electrophysiological-optogenetic synchronous recording signals

Abstract

The present application relates to a kind of electro-physiological-optical genetic synchronous recording signal light artifact identification method and device, wherein, method includes: obtaining electro-physiological-optical genetic synchronous recording signal, and extracting the sharp peak event in the electro-physiological-optical genetic synchronous recording signal;According to the sharp peak event extracted obtains multiple sharp peak waveforms, extracts waveform morphological features to each sharp peak waveform, and according to waveform morphological features, filters out the first suspected light artifact sharp peak event;Utilize the physical time association of light artifact and laser pulse to the sharp peak event verification, obtain the second suspected light artifact sharp peak event;The intersection of first suspected light artifact sharp peak event and second suspected light artifact sharp peak event is sought, and light artifact is obtained.The present application can accurately identify the light artifact in electro-physiological-optical genetic synchronous recording signal.

Description

Technical Field

[0001] This invention relates to the field of optical artifact recognition technology, and in particular to a method and apparatus for optical artifact recognition of electrophysiological-optogenetic synchronous recording signals. Background Technology

[0002] The synergy between optogenetic stimulation and multichannel electrophysiological recording has become an important direction for circuit causality testing in basic neuroscience research. Optogenetics, by delivering genes encoding photosensitive ion channels or photosensitive ion pumps to target neuronal populations for expression, allows researchers to trigger excitation and inhibition within milliseconds, establishing a clearer chain of causal evidence between stimulus input, neural response, and behavioral output. However, during simultaneous recording, optical artifacts are prevalent, severely affecting the identification of genuine neural signals. In high-impedance channels, the amplitude of artifacts can even drown out genuine neuronal signals. Because the Becquerel effect is inherently unavoidable in optogenetic-electrophysiological recording, effectively identifying and removing artifacts from the recorded signals is crucial for subsequent signal analysis.

[0003] Currently, the suppression of optical artifacts mainly involves the following approaches:

[0004] Hardware shielding: The aim is to block interference at the physical level by adding a metal shielding layer to the probe or changing the distance and angle between the electrode and the light emission point to control artifacts. These methods can reduce the artifact amplitude to some extent, but they cannot completely eliminate artifacts caused by photoinduced transient currents (Becquerel effect), and they also increase the complexity of the device structure and the difficulty of the manufacturing process.

[0005] Post-processing of the algorithm: When hardware suppression still has residual effects, efficient algorithms can effectively separate and identify the signal. For example, template subtraction algorithms require acquiring a pure artifact template without neural signals and subtracting it from the original signal. This method requires highly repetitive artifact waveforms and a stable recording environment, but in long-term in vivo experiments, the electrode-tissue interface impedance will drift over time, causing template failure.

[0006] Materials and Surface Modification: These approaches suppress artifacts at the electrochemical level by altering electrode materials or adding coatings to the electrode surface. One of the core sources of optical artifacts is the Becquerel effect caused by light irradiating the metal electrode surface, so modifying the electrode material is an important way to solve this problem. However, this type of method requires consideration of the potential impact of different materials on the signal-to-noise ratio, long-term in vivo stability, and safety. Summary of the Invention

[0007] The technical problem to be solved by the present invention is to provide a method and device for identifying optical artifacts in electrophysiological-optogenetic synchronous recording signals, which can accurately identify optical artifacts.

[0008] The technical solution adopted by this invention to solve its technical problem is: to provide a method for identifying optical artifacts in electrophysiological-optogenetic synchronous recording signals, comprising the following steps:

[0009] Acquire electrophysiological-optical genetic synchronous recording signals and extract spike-like events from the electrophysiological-optical genetic synchronous recording signals;

[0010] Multiple peak waveforms are obtained from the extracted peak-like events. Waveform morphology features are extracted from each peak waveform, and the peak-like events suspected to be optical artifacts are selected based on the waveform morphology features.

[0011] The spike-like event was verified by using the physical time correlation between optical artifacts and laser pulses to obtain the spike-like event of the second suspected optical artifact.

[0012] Find the intersection of the spike-like events of the first suspected optical artifact and the spike-like events of the second suspected optical artifact to obtain the optical artifact.

[0013] The extraction of spike-like events from the electrophysiological-optical genetic synchronous recording signal specifically includes:

[0014] The electrophysiological-optogenetic synchronous recording signal was bandpass filtered to remove low-frequency field potentials and high-frequency noise;

[0015] Threshold detection was used to extract all spike-like events exceeding the threshold in the electrophysiological-optogenetic synchronous recording signal after bandpass filtering.

[0016] The process involves obtaining multiple peak waveforms based on the extracted peak-like events, extracting waveform morphology features from each peak waveform, and filtering out the peak-like events that are suspected to be optical artifacts based on these features. Specifically, this includes:

[0017] Using the peak point of each spike-like event as the center, waveform segments of preset time periods before and after are extracted, and the time base of all spike-like events is uniformly aligned to obtain multiple spike waveforms;

[0018] For each peak waveform, calculate the post-hyperpolarization depth, peak-to-valley amplitude, and half-peak width; where the post-hyperpolarization depth is the difference between the peak point and the subsequent lowest valley point; the peak-to-valley amplitude is the difference between the peak point and the previous depolarization valley point; and the half-peak width is the time interval between the rising and falling edges of the waveform reaching half of the peak value.

[0019] Spike waveforms that simultaneously satisfy the following conditions are selected as the first suspected optical artifact spike-like events: post-hyperpolarization depth less than post-hyperpolarization depth threshold, peak-to-valley amplitude less than peak-to-valley amplitude threshold, and half-width at half-maximum (WWHM) exceeding the WWHM threshold.

[0020] The post-hyperpolarization depth threshold is 50 μV, the peak-to-valley amplitude threshold is 150 μV, and the half-peak width threshold is 0.2 ms.

[0021] The method of verifying the spike-like event by utilizing the physical time correlation between optical artifacts and laser pulses to obtain the spike-like event of the second suspected optical artifact includes:

[0022] Align the timestamp of the spike-like event with the timestamp of the rising edge of the laser pulse, and calculate the time delay of each spike-like event relative to the rising edge of the most recent laser pulse; if the time delay is less than the time delay threshold, the spike-like event is a spike-like event of the second suspected optical artifact.

[0023] Calculate the interpeak interval sequence for each spike-like event, and calculate the phase-locking index of the interpeak interval sequence with the photostimulation period; when the phase-locking index of the photostimulation period exceeds the phase-locking index threshold, and the distribution of the interpeak interval sequence conforms to an integer multiple of the photostimulation period, the spike-like event is a spike-like event of the second suspected photoartifact.

[0024] The time delay threshold is 0.5ms; the phase-locked index threshold is 0.8.

[0025] The technical solution adopted by this invention to solve its technical problem is: to provide a device for identifying optical artifacts in electrophysiological-optogenetic synchronous recording signals, comprising:

[0026] The acquisition and extraction module is used to acquire electrophysiological-optical genetic synchronous recording signals and extract spike-like events from the electrophysiological-optical genetic synchronous recording signals;

[0027] The filtering module is used to obtain multiple peak waveforms based on the extracted peak-like events, extract waveform morphology features for each peak waveform, and filter out the peak-like events that are suspected to be optical artifacts based on the waveform morphology features.

[0028] The verification module is used to verify the spike-like event by utilizing the physical time correlation between optical artifacts and laser pulses, and to obtain the spike-like event of the second suspected optical artifact.

[0029] The intersection module is used to find the intersection of the spike-like events of the first suspected optical artifact and the spike-like events of the second suspected optical artifact, thus obtaining the optical artifact.

[0030] The technical solution adopted by the present invention to solve its technical problem is: to provide an electronic device, including a memory, a processor, and a computer program stored in the memory and capable of running on the processor, wherein when the processor executes the computer program, it implements the steps of the above-mentioned method for identifying optical artifacts in electrophysiological-optogenetic synchronous recording signals.

[0031] The technical solution adopted by the present invention to solve its technical problem is: to provide a computer-readable storage medium on which a computer program is stored, wherein when the computer program is executed by a processor, the steps of the above-mentioned method for identifying optical artifacts in electrophysiological-optogenetic synchronous recording signals are implemented.

[0032] By employing the above-mentioned technical solutions, this invention has the following advantages and positive effects compared with existing technologies: This invention uses dual criteria for identification, where waveform features are based on the physiological basis of neuronal ion channel dynamics, and temporal features are based on the physical characteristics of the Becquerel effect, significantly improving identification accuracy and effectively solving the industry problem of artifacts obscuring the true signal in high-impedance channels. This invention uses a pure algorithm post-processing scheme, requiring no modification to any hardware or electrode materials, and can be directly integrated into existing electrophysiological acquisition and analysis systems, greatly reducing application costs and technical barriers. Furthermore, this invention is unaffected by electrode-tissue interface impedance drift and changes in the recording environment, exhibiting excellent long-term stability. It supports offline batch processing and online real-time processing, significantly improving the objectivity and reproducibility of optogenetic electrophysiological data and reducing the workload of manual verification. Attached Figure Description

[0033] Figure 1 This is a flowchart of the optical artifact identification method for electrophysiological-optogenetic synchronous recording signals according to the first embodiment of the present invention;

[0034] Figure 2 This is a schematic diagram comparing the key waveform features of photoelectric artifacts and real Spikes in the first embodiment of the present invention;

[0035] Figure 3 This is a comparison diagram of the temporal characteristics of spontaneous and light-induced real neuronal firing in the first embodiment of the present invention;

[0036] Figure 4 This is an extreme timing characteristic diagram of pure optical artifacts in a high-impedance channel according to the first embodiment of the present invention. Detailed Implementation

[0037] The present invention will be further illustrated below with reference to specific embodiments. It should be understood that these embodiments are for illustrative purposes only and are not intended to limit the scope of the invention. Furthermore, it should be understood that after reading the teachings of this invention, those skilled in the art can make various alterations or modifications to the invention, and these equivalent forms also fall within the scope defined by the appended claims.

[0038] The first embodiment of the present invention relates to a method for identifying optical artifacts in electrophysiological-optogenetic synchronous recording signals. This method constructs a dual-criteria rule based on waveform morphological features and temporal statistics, such as... Figure 1 As shown, the specific steps include:

[0039] Step 1: Acquire the electrophysiological-optical genetic synchronous recording signal and extract the spike-like events from the electrophysiological-optical genetic synchronous recording signal.

[0040] This step involves bandpass filtering (typically 300Hz-3000Hz) of the acquired electrophysiological-optical genetic synchronous recording signal to remove low-frequency field potentials and high-frequency noise. Then, a threshold detection method is used to extract all spike-like events exceeding a threshold value from the bandpass-filtered electrophysiological-optical genetic synchronous recording signal. The threshold value is typically 3-5 times the standard deviation of the electrophysiological-optical genetic synchronous recording signal. This process allows for the extraction of all spike-like events from the electrophysiological-optical genetic synchronous recording signal.

[0041] Step 2: Based on the extracted spike-like events, multiple spike waveforms are obtained. Waveform morphology features are extracted from each spike waveform, and the spike-like events suspected to be optical artifacts are selected based on these features. This step is detailed below:

[0042] First, taking the peak point of each spike-like event as the center, extract waveform segments with a preset time period (e.g., 1ms) before and after (total duration 2ms), and align the time base of all spike-like events to obtain multiple spike waveforms.

[0043] Then, for each peak waveform, the post-hyperpolarization depth, peak-to-valley amplitude, and half-peak width are calculated; where the post-hyperpolarization depth is the difference between the peak point and the subsequent lowest valley point; the peak-to-valley amplitude is the difference between the peak point and the previous depolarization valley point; and the half-peak width is the time interval between the rising edge and the falling edge of the waveform reaching half of the peak value.

[0044] like Figure 2 As shown, the inventors of this invention discovered that normal action potentials produce significant post-hyperpolarization with an amplitude of approximately 297 μV, while photoelectric artifacts exhibit almost no identifiable post-hyperpolarization component. This difference indicates that photoelectric artifacts lack ion channel-mediated repolarization and post-hyperpolarization processes, consistent with their non-biological origin. Furthermore, the peak-to-trough amplitude of a normal action potential is approximately 354 μV, while the amplitude of photoelectric artifacts is significantly lower, at approximately 142 μV. This may be because photoelectric artifacts are primarily caused by transient currents generated at the electrode interface due to light irradiation, rather than transmembrane potential changes driven by ion channel opening and closing. The half-width at half-maximum (HWHM) of photoelectric artifacts is approximately 0.35 ms, significantly wider than the approximately 0.13 ms of a normal action potential, indicating that its waveform is more diffuse on a time scale and lacks the sharp characteristics formed by the rapid sequential activation of sodium and potassium channels in neuronal action potentials. Based on the above findings, this embodiment sets the post-hyperpolarization depth threshold to 50 μV, the peak-to-trough amplitude threshold to 150 μV, and the HWHM threshold to 0.2 ms.

[0045] Finally, peak waveforms that simultaneously satisfy the following conditions—post-hyperpolarization depth less than a post-hyperpolarization depth threshold, peak-to-valley amplitude less than a peak-to-valley amplitude threshold, and full width at half maximum (FWHM) exceeding a FWHM threshold—are selected as peak-like events of the first suspected optical artifact. Peak waveforms that do not simultaneously satisfy the above conditions are selected as first candidate real peaks.

[0046] Step 3: Verify the spike-like event by using the physical time correlation between the optical artifact and the laser pulse to obtain the spike-like event of the second suspected optical artifact.

[0047] like Figure 3 and Figure 4 As shown, the inventors of this invention also discovered that the action potentials of real neurons exhibit a relatively consistent time lock near the rising edge of each laser pulse, showing a significant instantaneous peak within a short latency window (approximately 2 ms) after the laser pulse begins. In contrast, photoartifacts are directly generated by the Becquerel effect, forming a photocurrent at the instant the laser is turned on or off, a process with almost no delay. Furthermore, photoartifacts only appear at the moment of light stimulation, thus exhibiting strict periodicity, while the ISI of real action potentials is more dispersed, consistent with the characteristics of physiological refractory period and random firing. Based on the above findings, this step specifically includes:

[0048] Align the timestamps of the spike-like events with the timestamps of the rising edges of the laser pulses, and calculate the time delay of each spike-like event relative to the rising edge of the most recent laser pulse. If the time delay is less than the time delay threshold (0.5ms), the spike-like event is a spike-like event of the second suspected optical artifact. If the spike-like events are concentrated within a window of 1.5-2.5ms after the rising edge of the laser (there is a neuronal activation latency), then the spike-like event is marked as the second candidate real spike.

[0049] This step also calculates the interpeak interval (ISI) sequence for each spike-like event and the phase-locking index of the ISI sequence with respect to the photostimulation period. The phase-locking index is used to assess whether there is a significant multiple relationship between the ISI and the photostimulation period. Its value ranges from 0 to 1, where 1 indicates complete phase-locking. When the phase-locking index of the photostimulation period exceeds the phase-locking index threshold (e.g., 0.8) and the distribution of the interpeak interval sequence conforms to an integer multiple of the photostimulation period, the spike-like event is a spike-like event of the second suspected photoartifact. When the phase-locking index of the photostimulation period is less than the phase-locking index threshold and the ISI sequence is dispersed, the spike-like event is a second candidate true spike.

[0050] It is worth mentioning that the above-mentioned time delay judgment and phase-locked index judgment can be performed together or one of them can be performed separately.

[0051] Step 4: Find the intersection of the spike-like events of the first suspected optical artifact and the second suspected optical artifact to obtain the optical artifact. In this step, spike-like events simultaneously marked as the first suspected optical artifact and the second suspected optical artifact are determined to be optical artifacts. Similarly, this step can also determine spike-like events simultaneously marked as the first candidate real spike and the second candidate real spike as real spikes. Spike-like events that satisfy only a single criterion (accounting for less than 5%) can be manually reviewed.

[0052] After obtaining the optical artifacts, they can be removed, thereby improving signal quality.

[0053] It is easy to see that this invention employs dual criteria for identification: waveform features are based on the physiological basis of neuronal ion channel dynamics, while temporal features are based on the physical properties of the Becquerel effect. This significantly improves identification accuracy and effectively solves the industry problem of artifacts obscuring the true signal in high-impedance channels. This invention uses a pure algorithm-based post-processing scheme, requiring no modification to any hardware or electrode materials. It can be directly integrated into existing electrophysiological acquisition and analysis systems, greatly reducing application costs and technical barriers. Furthermore, this invention is unaffected by electrode-tissue interface impedance drift or changes in the recording environment, exhibiting excellent long-term stability. It supports offline batch processing and online real-time processing, significantly improving the objectivity and reproducibility of optogenetic electrophysiological data and reducing the workload of manual verification.

[0054] The second embodiment of the present invention relates to an optical artifact recognition device for electrophysiological-optogenetic synchronous recording signals, comprising:

[0055] The acquisition and extraction module is used to acquire electrophysiological-optical genetic synchronous recording signals and extract spike-like events from the electrophysiological-optical genetic synchronous recording signals;

[0056] The filtering module is used to obtain multiple peak waveforms based on the extracted peak-like events, extract waveform morphology features for each peak waveform, and filter out the peak-like events that are suspected to be optical artifacts based on the waveform morphology features.

[0057] The verification module is used to verify the spike-like event by utilizing the physical time correlation between optical artifacts and laser pulses, and to obtain the spike-like event of the second suspected optical artifact.

[0058] The intersection module is used to find the intersection of the spike-like events of the first suspected optical artifact and the spike-like events of the second suspected optical artifact, thus obtaining the optical artifact.

[0059] The acquisition / extraction module includes:

[0060] The filtering unit is used to perform bandpass filtering on the electrophysiological-optical genetic synchronous recording signal to remove low-frequency field potentials and high-frequency noise;

[0061] The threshold detection unit is used to extract all spike-like events exceeding the threshold in the electrophysiological-optogenetic synchronous recording signal after bandpass filtering using the threshold detection method.

[0062] The filtering module includes:

[0063] The segmentation unit is used to extract waveform segments of a preset time period before and after each peak event, centering on the peak point of each peak event, and uniformly aligning the time base of all peak events to obtain multiple peak waveforms.

[0064] The calculation unit is used to calculate the hyperpolarization depth, peak-to-valley amplitude, and half-peak width for each peak waveform; where the hyperpolarization depth is the difference between the peak point and the subsequent lowest valley point; the peak-to-valley amplitude is the difference between the peak point and the previous depolarization valley point; and the half-peak width is the time interval between the rising edge and the falling edge of the waveform reaching half of the peak value.

[0065] The filtering unit is used to filter out peak waveforms that simultaneously satisfy the following conditions: post-hyperpolarization depth is less than the post-hyperpolarization depth threshold, peak-to-valley amplitude is less than the peak-to-valley amplitude threshold, and half-width exceeds the half-width threshold, as peak-like events of the first suspected optical artifact.

[0066] The post-hyperpolarization depth threshold is 50 μV, the peak-to-valley amplitude threshold is 150 μV, and the half-peak width threshold is 0.2 ms.

[0067] The verification module includes:

[0068] A time-locking unit is used to align the timestamp of the spike-like event with the rising edge timestamp of the laser pulse, calculate the time delay of each spike-like event relative to the rising edge of the most recent laser pulse; when the time delay is less than a time delay threshold, the spike-like event is a spike-like event of the second suspected optical artifact; and / or,

[0069] The periodic phase-locked index unit is used to calculate the inter-peak interval sequence of each spike-like event and to calculate the phase-locked index of the inter-peak interval sequence with the photostimulation period. When the phase-locked index of the photostimulation period exceeds the phase-locked index threshold and the distribution of the inter-peak interval sequence conforms to an integer multiple of the photostimulation period, the spike-like event is a spike-like event of the second suspected photoartifact.

[0070] The time delay threshold is 0.5ms; the phase-locked index threshold is 0.8.

[0071] The third embodiment of the present invention relates to an electronic device, including a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein the processor executes the computer program to implement the steps of the optical artifact identification method for electrophysiological-optogenetic synchronous recording signals of the first embodiment.

[0072] The fourth embodiment of the present invention relates to a computer-readable storage medium having a computer program stored thereon, wherein the computer program, when executed by a processor, implements the steps of the optical artifact identification method for electrophysiological-optogenetic synchronous recording signals of the first embodiment.

[0073] Those skilled in the art will understand that embodiments of this application can be provided as methods, systems, or computer program products. Therefore, this application can take the form of a completely hardware embodiment, a completely software embodiment, or an embodiment combining software and hardware aspects. Furthermore, this application can take the form of a computer program product implemented on one or more computer-usable storage media (including, but not limited to, disk storage and optical storage) containing computer-usable program code.

[0074] This application is described with reference to flowchart illustrations and / or block diagrams of methods, apparatus (systems), and computer program products according to embodiments of this application. It will be understood that each block of the flowchart illustrations and / or block diagrams, and combinations of blocks in the flowchart illustrations and / or block diagrams, can be implemented by computer program instructions. These computer program instructions can be provided to a processor of a general-purpose computer, special-purpose computer, embedded processor, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, generate instructions for implementing the flowchart... Figure 1 One or more processes and / or boxes Figure 1 A device that provides the functions specified in one or more boxes.

[0075] These computer program instructions may also be stored in a computer-readable storage medium that can direct a computer or other programmable data processing device to operate in a particular manner, such that the instructions stored in the computer-readable storage medium produce an article of manufacture including instruction methods implemented in a process. Figure 1 One or more processes and / or boxes Figure 1 The function specified in one or more boxes.

[0076] These computer program instructions may also be loaded onto a computer or other programmable data processing equipment to cause a series of operational steps to be performed on the computer or other programmable equipment to produce a computer-implemented process, thereby providing instructions that execute on the computer or other programmable equipment for implementing the process. Figure 1 One or more processes and / or boxes Figure 1 The steps of the function specified in one or more boxes.

[0077] The above description is merely a specific embodiment of the present invention, but the scope of protection of the present invention is not limited thereto. Any variations or substitutions that can be easily conceived by those skilled in the art within the technical scope disclosed in the present invention should be included within the scope of protection of the present invention. Therefore, the scope of protection of the present invention should be determined by the scope of the claims.

Claims

1. A method for identifying optical artifacts in electrophysiological-optogenetic synchronous recording signals, characterized in that, Includes the following steps: Acquire electrophysiological-optical genetic synchronous recording signals and extract spike-like events from the electrophysiological-optical genetic synchronous recording signals; Multiple peak waveforms are obtained from the extracted peak-like events. Waveform morphology features are extracted from each peak waveform, and the peak-like events suspected to be optical artifacts are selected based on the waveform morphology features. The spike-like event was verified by using the physical time correlation between optical artifacts and laser pulses to obtain the spike-like event of the second suspected optical artifact. Find the intersection of the spike-like events of the first suspected optical artifact and the spike-like events of the second suspected optical artifact to obtain the optical artifact.

2. The method for identifying optical artifacts in electrophysiological-optogenetic synchronous recording signals according to claim 1, characterized in that, The extraction of spike-like events from the electrophysiological-optical genetic synchronous recording signal specifically includes: The electrophysiological-optogenetic synchronous recording signal was bandpass filtered to remove low-frequency field potentials and high-frequency noise; Threshold detection was used to extract all spike-like events exceeding the threshold in the electrophysiological-optogenetic synchronous recording signal after bandpass filtering.

3. The method for identifying optical artifacts in electrophysiological-optogenetic synchronous recording signals according to claim 1, characterized in that, The process involves obtaining multiple peak waveforms based on the extracted peak-like events, extracting waveform morphology features from each peak waveform, and filtering out the peak-like events that are suspected to be optical artifacts based on these features. Specifically, this includes: Using the peak point of each spike-like event as the center, waveform segments of preset time periods before and after are extracted, and the time base of all spike-like events is uniformly aligned to obtain multiple spike waveforms; For each peak waveform, calculate the post-hyperpolarization depth, peak-to-valley amplitude, and half-peak width; where the post-hyperpolarization depth is the difference between the peak point and the subsequent lowest valley point; the peak-to-valley amplitude is the difference between the peak point and the previous depolarization valley point; and the half-peak width is the time interval between the rising and falling edges of the waveform reaching half of the peak value. Spike waveforms that simultaneously satisfy the following conditions are selected as the first suspected optical artifact spike-like events: post-hyperpolarization depth less than post-hyperpolarization depth threshold, peak-to-valley amplitude less than peak-to-valley amplitude threshold, and half-width at half-maximum (WWHM) exceeding the WWHM threshold.

4. The method for identifying optical artifacts in electrophysiological-optogenetic synchronous recording signals according to claim 3, characterized in that, The post-hyperpolarization depth threshold is 50 μV, the peak-to-valley amplitude threshold is 150 μV, and the half-peak width threshold is 0.2 ms.

5. The method for identifying optical artifacts in electrophysiological-optogenetic synchronous recording signals according to claim 1, characterized in that, The method of verifying the spike-like event by utilizing the physical time correlation between optical artifacts and laser pulses to obtain the spike-like event of the second suspected optical artifact includes: Align the timestamps of the spike-like events with the rising edge timestamps of the laser pulses, and calculate the time delay of each spike-like event relative to the rising edge of the most recent laser pulse; if the time delay is less than a time delay threshold, the spike-like event is a spike-like event of the second suspected optical artifact; and / or, Calculate the interpeak interval sequence for each spike-like event, and calculate the phase-locking index of the interpeak interval sequence with the photostimulation period; when the phase-locking index of the photostimulation period exceeds the phase-locking index threshold, and the distribution of the interpeak interval sequence conforms to an integer multiple of the photostimulation period, the spike-like event is a spike-like event of the second suspected photoartifact.

6. The method for identifying optical artifacts in electrophysiological-optogenetic synchronous recording signals according to claim 5, characterized in that, The time delay threshold is 0.5ms; the phase-locked index threshold is 0.

8.

7. A device for identifying optical artifacts in electrophysiological-optogenetic synchronous recording signals, characterized in that, include: The acquisition and extraction module is used to acquire electrophysiological-optical genetic synchronous recording signals and extract spike-like events from the electrophysiological-optical genetic synchronous recording signals; The filtering module is used to obtain multiple peak waveforms based on the extracted peak-like events, extract waveform morphology features for each peak waveform, and filter out the peak-like events that are suspected to be optical artifacts based on the waveform morphology features. The verification module is used to verify the spike-like event by utilizing the physical time correlation between optical artifacts and laser pulses, and to obtain the spike-like event of the second suspected optical artifact. The intersection module is used to find the intersection of the spike-like events of the first suspected optical artifact and the spike-like events of the second suspected optical artifact, thus obtaining the optical artifact.

8. An electronic device comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, characterized in that, When the processor executes the computer program, it implements the steps of the optical artifact recognition method for electrophysiological-optogenetic synchronous recording signals as described in any one of claims 1-6.

9. A computer-readable storage medium having a computer program stored thereon, characterized in that, When the computer program is executed by the processor, it implements the steps of the optical artifact recognition method for electrophysiological-optogenetic synchronous recording signals as described in any one of claims 1-6.

Images