as indicated above , embodiments of the disclosure are directed to car monitoring and diagnostic systems , methods and devices . an example of one such system 10 is illustrated in fig1 . as presented here , the system includes a sensor array 12 . the sensor array 12 , an example of which is shown in detail in fig1 a and 2 is comprised of : near infrared spectroscopy ( nirs ) oximetry elements 14 for optical emission and detection of optical signals ; an applanation tonometry pressure sensing device 16 ; a sensor patch 18 of compliant material that contains the nirs sensor 14 and applanation tonometry pressure sensor 16 in a desired geometric configuration ; a mechanism for digitization , preamplification , and prefiltering of the nirs and pressure signals ; and a mechanism of data transmission of the signals ( cable , bluetooth ® transmitter , etc .) 20 to a signal processing member 22 . applanation tonometry measures relative blood pressure changes , such that pressure waveform morphology and lead time relationships can be calculated . in addition to the sensor array , the embodiment shown in fig1 includes a headset frame 24 to which the sensor array 12 is mounted . in this embodiment shown in fig1 aspects of the sensor array are mounted at selected points on the headset . for example , the nirs oximetry sensing member 14 is firmly appositioned to the forehead scalp , or other scalp area of the patient . correct positioning of the nirs sensor 14 is achieved by tensioning of an adjustable band 26 that serves the purpose of providing a firm headset fit for a range of cranium sizes . the applanation pressure sensing member 16 is positioned between two force - bearing contact pads 18 in order to mechanically isolate the apposition pressure of the applanation sensor 16 from the headset pressure due to the tensioning of the adjustable band 26 . in at least one embodiment , the nirs sensor 14 with nirs optical elements 15 and signal amplifiers is mounted in the headset and located over a region of the forehead with a secure apposition of the nirs sensor 14 to the skin . separately the applanation tonometry pressure ( atp ) sensor 16 is mounted in the same headset 24 such that the atp sensor 16 is located over a region of the superficial temporal artery for optimal sensitivity of arterial pressure variations . in some embodiments , the headset 24 includes a feature that may be located in a preferred location relative to an anatomical feature of the test subject . for example , the positioning of the atp sensor 16 is optimally located over a preferred segment of the superficial temporal artery . in at least one embodiment , the anatomical feature is the tragus , such that the feature of the headset design enables positioning of the atp sensor 16 near , or in a desired position relative to , the tragus . in at least one embodiment , the atp sensor 16 is distinct or separate from the headset 24 so that it may be placed over the radial artery in the wrist of a patient , instead of over the temporal artery . in at least one embodiment , the atp sensor 16 is placed over the brachial artery for arterial pressure applanation . alternately , the pressure sensor 16 may be placed noninvasively over any artery in which pressure variations may be transmitted to the signal processor 22 . as is shown in fig3 analysis of the nirs sensor 14 and pressure sensor 16 data is performed by the signal processing element 22 to determine , among other things , time delays between the parameters , phase differences , and related characteristics of nirs parameters and blood pressure variations . results are displayed on the graphical user interphase ( gui ) 30 . the signal processing element 22 is in communication with the sensor array 12 by direct ( wired ) or wireless communication comprised of : an input section 32 capable of receiving wireless data from the nirs sensor 14 and pressure sensor 16 ; amplifiers for bandpass filtering and amplification of the applanation blood pressure signal ; filtering and pre - conditioning of the optical signal from the nirs optical detection element ; hardware and software 34 for implementing an algorithm for analysis of the nirs and pressure signals to generate a cross correlation result representative of the lead time between the nirs and pressure signals ; and display of the analysis results in a graphic user interface ( gui ) 30 showing numerical and trending data . in at least one embodiment , the signal processor 22 includes a processing algorithm for analysis of the raw nirs signals and produces an estimate of the oxy - hemoglobin concentration , deoxy - hemoglobin concentration , total hemoglobin concentration , or other parameter derived from the raw nirs signals , prior to analysis to generate a cross correlation result . in at least one embodiment of the invention the sensor array 12 detects and measures two physiological parameters , namely blood pressure ( via the atp sensor or module 16 ) and cerebral oxygen concentration ( via the nirs sensor 14 ) over multiple periods of time or episodes . episodes of interest are identified in the time series of the first parameter , then an assessment is made of the presence and temporal relationships of corresponding episodes in the second parameter . in some embodiments , the first parameter may be arterial applanation tonometry pressure measured by applanation of an artery . the second parameter may be blood velocity measured in a cerebral artery , or regional cerebral oxygen saturation , as a surrogate for time - varying changes in cerebral blood flow . in embodiments of the disclosure , the system may be designed , with the use of signal bandpass filtering to , measure time - varying magnitude changes occurring only over relatively long intervals , such as 30 - 120 seconds or even several minutes ; or only over relatively shorter time intervals , such as 3 - 15 seconds , in order to evaluate physiological effects occurring within a certain range of time scales . cerebrovascular blood flow autoregulation may be measured with the device by an evaluation of a series of time delays of regional oxygen saturation waveforms relative to arterial applanation tonometry pressure waveforms during episodes of interest in the pressure parameter . as discussed , the system of fig1 - 3 includes a device for measuring cerebral pressure autoregulation by the use of arterial blood pressure and regional cerebral oxygen saturation ( first and second parameters ). the device detects whether slow arterial blood pressure changes are having an effect on cerebral blood flow , and by this means assesses whether cerebral pressure autoregulation is functional . a continuous measure of arterial pressure is used to track slow blood pressure changes . regional oxygen saturation is measured to track the cerebral blood flow changes . this substitution is possible because , over the time scales of interest , cerebral blood flow changes are directly reflected in regional cerebral oxygen saturation changes . example time series of monitored first 40 and second 50 physiological parameters are shown in fig4 and 5 respectively . fig4 depicts a time series of measurements of a first parameter ( blood pressure ) 40 during a 4 - minute episode . slow arterial pressure variations are not cyclic changes with a regular frequency , but are highly non - stationary . this is readily seen in the example of fig4 . the various devices and approaches described in the prior art for measurement of cerebral pressure autoregulation all lack appreciation for this characteristic . therefore it is a novel aspect of the present invention that slow arterial pressure variations are each considered as a separate event , with no requirement for stationary cyclic behavior of the parameter . in order to measure the effect of arterial blood pressure ( abp ) changes 40 on regional cerebral oxygen ( rso2 ) 50 , the system measures the correspondence of rso2 changes relative to atp that occur over time scales of 30 - 120 seconds , such as in the manner shown in fig6 and 11 . atp and rso2 time series measurements may be sampled at a higher data rate and then digitally filtered to remove energy outside the bandwidth of 0 . 008 - 0 . 033 hz , which is not of interest . a change in atp is characterized as a waveform which has a beginning at a minimum value , followed by a rising segment during which atp increases to a maximum value , followed ( in some cases ) by a maximum plateau , followed by a decreasing segment during which atp decreases and returns to a minimum value , followed ( in some cases ) by a baseline segment which precedes the beginning of the next waveform . the rising segment of a waveform is referred to herein as an “ episode of interest ”, which is an abp response to a neurogenic effect which raises blood pressure . these episodes of interest are identified as e 1 , e 2 , etc in fig4 - 9 . the time intervals e 1 - e 5 shown in fig4 identify episodes of interest in the first parameter atp . in this example , an episode of interest is the occurrence of an increasing signal magnitude beginning from a local minimum value and ending at a local maximum value , which includes a zero - crossing or crossing of a threshold value . the episodes of interest have varying temporal characteristics , each corresponding to the rising portion of a waveform . fig5 shows a time series of measurements of a second parameter of cerebral oxygen or rso2 during a 3 - minute episode . the episodes of interest e 1 - e 4 , corresponding to the rising portion of the most significant waveforms , are identified for illustration purposes . in practice the episodes used for evaluation are defined by the episodes of interest in only one parameter , and typically the first parameter . fig6 shows an example 4 - minute time series of both the first and the second parameters over a time interval which includes five episodes of interest ( e 1 - e 5 ) in the first parameter , represented by the solid line . the second parameter , represented by the dashed line , shows no similarity to the first parameter . fig7 shows an example 4 - minute time series of both the first and the second parameters over a time interval which includes four episodes of interest ( e 1 - e 4 ) in the first parameter . in this instance , the second parameter ( dashed line ) exhibits a significant similarity of pattern to the first parameter . fig8 shows an example measurement of the lead times of the waveforms of the second parameter relative to the first parameter in fig6 , as measured during each of the five episodes of interest . the lead time alternates between large negative and positive values , indicating that the second parameter episodes ( or waveforms ) are independent of the first parameter . fig9 shows a measurement of time delays of the waveforms of the second parameter relative to the first parameter in fig7 , as measured during the four episodes of interest . in this case the lead times are all within a relatively narrow range of values , indicating a consistent temporal relationship between the second parameter episodes ( or waveforms ) and those of the first parameter . by visual examination of fig7 it is evident that the second parameter waveforms do indeed show a consistent similarity of pattern in which the episodes of interest ( i . e ., the rising segments of the waveforms ) are time - shifted to the left relative to the first parameter episodes of interest . the rso2 ( second parameter ) changes depicted on the examples of fig5 - 9 are also characterized as waveforms which have the same basic segments . the rso2 waveforms may not have a similar morphology to that of the abp ( first parameter ) waveforms . a typical difference is illustrated by comparison of fig4 and 5 . slow rso2 waveforms may be the result of cerebral blood flow responses to arterial carbon dioxide or cerebral metabolic processes , as noted above , or they may be caused by a failure of the cerebral pressure autoregulation ( i . e ., car ) mechanisms to adjust cerebrovascular resistance in order to minimize the effects of arterial pressure on cerebral blood flow . system of fig1 - 3 measures the contributory effect of abp as a cause of rso2 changes by evaluating the correspondence of episodes of interest ( i . e ., the rising segments ) of rso2 waveforms relative to contemporaneous atp waveforms . absence of rso2 waveforms indicates no effect of abp variations on rso2 , and therefore intact car . in the presence of rso2 waveforms , a lack of correspondence indicates that the rso2 changes , or waveforms , are the result of other effects , and therefore that car is functional . close correspondence between the episodes of interest in the two sets of waveforms indicates that the rso2 changes are correlated with abp , and therefore that cerebral pressure autoregulation is not fully functional . correspondence is evaluated by measuring the lead time of the rso2 ( second parameter ) versus atp ( first parameter ) episodes of interest for each of a series of atp waveforms , such as is depicted in fig8 and 9 . in practice the lead times may be either positive or negative . consecutive small or zero lead times over a series of waveforms shows a high degree of correspondence , whereas large lead time variability shows lack of correspondence . referring to the figures , episodes of interest are first identified in a series of atp and rso2 waveforms , as illustrated in fig4 and 5 . for each atp waveform , a measurement is made of the average lead time of the most contemporaneous rso2 waveform during the atp episode of interest . fig6 illustrates a series of atp and rso2 waveforms which do not show correspondence with each other . fig7 illustrates a series of atp and rso2 waveforms that show correspondence , in which a rso2 ( dashed line ) waveform lead time is evident . fig8 is an example of a device output for the waveforms shown in fig6 . the rso2 lead times vary from − 20 to + 15 seconds , confirming a lack of correspondence which is visually evident in fig6 . the interpretation of the device output would be that cerebral pressure autoregulation is healthy and functional . fig9 shows the device output for the waveforms shown in fig7 , in which the last three rso2 lead times are consistently between 5 and 10 seconds . this is evidence of a corresponding effect of abp on rso2 changes , and absence of cerebrovascular resistance adjustments . the interpretation of the device output would be dysfunction of cerebral pressure autoregulation . analysis and comparison of the monitored parameters may be accomplished in a variety of ways and mechanisms by the signal processor ( see fig3 ). in some embodiments , the digitized nirs and pressure signals are bandpass filtered by the signal processor to remove signal frequencies outside the range of interest . as indicated above , the frequency range of interest may be 0 . 008 hz to 0 . 033 hz ; alternative higher frequency ranges may also be used for analysis , including 0 . 008 - 0 . 06 hz , or other bandwidth . in some embodiments , an initial analysis of the slow wave features of either the nirs or atp signals is made to determine an optimal frequency range of interest . for example , fourier analysis may be made to determine whether a significant proportion of the low frequency power is between 0 . 03 and 0 . 06 hz . alternatively , an analysis of the wavelengths of slow waves may be made to determine whether a significant proportion of slow waves occur with wavelengths from 16 - 33 seconds . from this initial analysis , a frequency range of interest for bandpass filtering and lead time calculations may be chosen . in at least one embodiment , the signal filtering is performed using a kaiser - bessel filter whose coefficients have been selected to produce a zero phase shift . in at least one embodiment , calculation of the lead time relationship between the bandpass filtered nirs parameter and atp signal is performed for each discrete pressure waveform in a time series of the data by generating a cross correlation function ( ccf ) based on a predetermined portion of the atp waveform and the corresponding nirs waveform as defined above . in such an embodiment , the time delay nτ corresponding to the maximum value of the cross correlation function f ( n ), in which τ is the time interval between successive digitized data points , is interpreted to represent a real - time estimate of the time delay between the two filtered signals . in some embodiments , a phase delay is calculated in degrees as the ratio of nτ divided by the wavelength of the atp waveform , multiplied by 360 degrees . in another embodiment , the lead time relationship for each discrete pressure waveform within a selected recent epoch of time ( the “ phase history ”), as illustrated in fig8 and 10 , may be displayed in the gui as shown in fig3 . in a related embodiment , calculation of the phase relationship is performed for each discrete nirs waveform , instead of each discrete pressure waveform , in a time series of the data . in some embodiments , the lead time relationships within the selected epoch may be further analyzed with statistical methods to determine a trend ; characterize variability ( e . g ., calculate standard deviation ); calculate a mean lead time relationship ; estimate the probability that the mean lead time is greater or less than a predetermined lead time , which may be defined as a critical or threshold value ; or calculate a confidence interval defining a range within which the mean lead time likely occurs . one or more of these outputs may be used to characterize the extent to which nirs fluctuating values are passively dependent on pressure variations , as an index of car dysfunction . in another embodiment , the analysis may be performed only on atp waveforms that meet certain criteria . alternatively , the analysis may be performed only when nirs waveforms meet certain criteria , or only when both the atp and corresponding nirs waveforms respectively meet certain criteria . for example , exclusion criteria may include low amplitude , high amplitude , presence of artifact ( e . g ., motion artifact ), large difference in wavelength relative to the wavelength of one or more previous waveforms , or large difference in wavelength between the pressure waveform and corresponding nirs waveform . in at least one embodiment , the atp sensing module and headset include structural features that prevent transfer of the headset tension to the applanation pressure sensing member . the structural features may include force bearing members on opposite sides of the atp sensing member , combined with a constant force device within the atp sensing module that applies a constant contact force for the atp sensing member against the skin . a constant force device may include a spring loading , a pressurized bladder , or other device . in at least one embodiment , a measurement procedure may be performed in the following steps : first , place the headset which includes the nirs sensing module and the atp sensing module on the subject to be tested , with the nirs sensor located over the forehead and the atp sensing member positioned over a segment of the temporal artery . fixation is facilitated by a member of the headset extending over the top of the head which has a means of adjusting the tension . second , verify that the nirs and pressure signals are free of artifacts that may result from poor placement technique , and adjust placement of the sensors as necessary . third , assure that the surroundings are suitable for stable measurements , and in particular remove noise and other cognitive distractions , while making the patient comfortable . next , begin data acquisition and analysis . in this part of the procedure , the signal processor receives the digitized data from the two sensors , and then applies bandpass filtering to successive overlapping epochs of data from each sensor . ideally the epochs will be sufficiently long to include at least two of the lowest frequency waveforms within the passband ( i . e ., at least two wavelengths ). for example , if the passband is 0 . 008 - 0 . 03 hz , the filtered epoch may be 2 / 0 . 008 = 250 seconds or more . next , the signal processor identifies the time of the beginning of the pressure rise in a first selected waveform ( t p , min ), and also the time of the peak pressure in the first selected waveform ( t p , max ). next , the signal processor calculates a non - normalized cross correlation function ( ccf ), defined as the sum of the products of the atp values and time - delayed nirs values over the interval from t p , min to t p , max for each of the predetermined time delay ( nτ ) values . therefore , for the first pressure waveform , where τ is the sampling interval ( i . e ., 1 / sample rate ), and f ( n ) is calculated over the time interval from t = t p , min to t p , max and over a predetermined range of n that is sufficient to capture a maximum value of f ( n ) ( typically the maximum nτ is about one half the wavelength of the pressure waveform , and f ( n ) is calculated for both positive and negative n values ). see fig1 . in at least some embodiments each of the two waveforms is normalized to a range of − 1 to + 1 prior to calculation of the ccf , and this normalization is performed separately for each successive waveform pair . next the signal processor determines the value of nτ for which f ( n ) is maximum , τ · n ccf max . this represents the calculated time delay of the nirs waveform relative to the pressure waveform . the value of n ccfmax may be negative or positive . next the signal processor identifies the time of the beginning of the pressure rise from its minimum level , t p , min for the next waveform . the time interval between the two minima is defined as the wavelength λ of the first pressure waveform . next the signal processor calculates a phase angle φ in degrees using the calculated τ · n ccf max and the calculated wavelength , λ , as follows : next the signal processor updates the gui with the numerical and graphical representation of the phase angle . next , after sufficient acquisition of new data , the signal processor applies bandpass filtering to a new epoch of data overlapping with the first epoch , and sufficiently long to include at least the next waveform . then using the t p , min already determined for the next waveform , the steps of determining t p , max and calculating the ccf , τ · n ccf max , λ , and φ will be repeated , and the gui will be updated . these steps will continue to repeat until a satisfactory characterization of the behavior of the lead time or phase angle φ can be made . for example , a stable , small lead time may provide a confident assessment within a few wavelengths of time , whereas a variable lead time may require a longer time for an assessment and characterization with acceptable confidence intervals . in some embodiments , the signal processing member may include several frequency passbands ; software for analysis of the relative phase delay between the nirs and pressure parameters within each passband ; an algorithm for determining the frequency dependence of the phase delay ; and an algorithm for calculating a maximum , minimum , or optimum coherence and phase delay between the two parameters . in some embodiments , the signal processing member may include an algorithm and software for determining the signal power or wavelengths in bandpass - filtered nirs and pressure parameters over a predetermined time interval . the predetermined time interval may be 5 minutes , or any other duration that is longer than 2 ÷ f m where f m is the minimum frequency within the frequency bandpass range . the algorithm and software may include a calculation of signal quality based on the signal power of either parameter . the algorithm and software may further use the signal quality calculation to determine a level of confidence in other calculated parameters . in another embodiment , a car assessment protocol includes cognitive stress to provoke a cerebral metabolic change while nirs and pressure data are acquired . cognitive stress may be induced by working memory tasks , problem solving , or other method . lead time calculations prior to and after the stress challenge may be compared for use as an index of dysfunction . in another embodiment , a car assessment protocol includes physical exertion to introduce a cardio - respiratory challenge and a cerebrovascular response while nirs and pressure data are acquired . physical exertion may be induced by isometric exercise , aerobic exercise , or other physical activity . in another embodiment , two or more nirs sensing members may be used . the sensing members may be placed bilaterally on the forehead and temporal regions . alternatively , the sensing members may be placed on two or more scalp locations that are located over different regions of the brain , with the single atp sensing member placed over the temporal artery , radial artery , or other artery . for example , one nirs sensing member may be placed on or near a location of contusion which is associated with a concussive injury , and another nirs sensing member may be placed ipsilaterally or contralaterally over a different region of the brain . the signal processing member may include software that compares the results from the sensing members . a local car dysfunction may be assessed on the basis of , for example , sufficiently large differences in analysis results between the two locations . a global car dysfunction may be assessed on the basis of similar car dysfunction indexes obtained at two or more locations . in another embodiment , a test protocol includes first having the subject breathe normally for a first time interval , preferably at a regular constant respiration rate , for a predetermined duration while nirs and pressure data are acquired and analyzed as described above , which may be 5 minutes or other duration . after some elapsed time , the subject may breathe normally for a second time interval , preferably at a regular constant respiration rate , for a predetermined duration while nirs and pressure data are acquired and analyzed as described above . the elapsed time may be 15 minutes , 30 minutes , or other time . this procedure may be repeated for an extended time or number of time intervals as may be desired for dysfunction assessment . an additional assessment protocol is then performed by the software in the signal processing member to calculate a level of dysfunction on the basis of one or more of the lead time assessments performed in the sequence of time intervals . in some embodiments , the test subject may follow a pattern of activity while nirs and pressure data are acquired , that may include normal daily routines , physical exertion , cognitive exertion , reading , eating , sleeping , or other activities . in this embodiment the test subject &# 39 ; s respiration rate may not be regular , but may vary according to what is a comfortable rate during each activity . the assessment protocol may include an algorithm and software that selects optimum epochs for car assessment . selection criteria for optimum epochs include : absence of apparent artifacts ; stable pressure and nirs waveform amplitudes ; or other criteria . in some embodiments , the input section of the signal processing member includes a means of entering event markers that may be used to identify activities of interest in order to facilitate analysis of test results . activities of interest may include physical exertion , cognitive exertion , hyperventilation , onset of neuropsychological symptoms , heightened stress response to external stimuli , or other stimulus . in another embodiment , the signal processing member includes a means of communicating the test results to an external device . the external device may be a smartphone , tablet computer , or other handheld device . communication to the external device may be either via a cable , such as a usb cable , or via a wireless protocol , such as bluetooth . the external device may include a means of communicating current test results , historical test results , test subject identification and information , and related data to a remote computer or server using an available communications protocol such as internet protocol , file transfer protocol , or other internet or telecommunication protocols . in at least one embodiment , a device such as a smartphone or tablet computer comprises at least a portion of the signal processing member and performs at least some of the functions of the signal processing member . the many features and advantages of the invention are apparent from the above description . numerous modifications and variations will readily occur to those skilled in the art . since such modifications are possible , the invention is not to be limited to the exact construction and operation illustrated and described . rather , the present invention should be limited only by the following claims .
a system device and method are presented for determining if cerebral blood flow autoregulation functionality has been compromised in a patient . the system includes sensors for detecting and measuring at least two physiological parameters such as oxygenation and blood pressure . the system analyzes these measurements and presents them as waveforms . episodes of interest in each waveform are determined and compared . any time lag between the occurrence of episodes of interest in each waveform is established and accounted for . where there is correspondence between episodes of interest in each waveform then impairment of cerebral blood flow autoregulation is likely .