# Patent

## Title
Methods and apparatus to determine locations of audience members

## Abstract

A disclosed example method to determine a location of an audience member involves generating a correlation analysis result based on correlating first audio samples from a stationary audio detector with second audio samples from a portable audio detector carried by the audience member and determining via a neural network the location of the audience member based on the correlation analysis result.

## Background
FIELD OF THE DISCLOSURE
The present disclosure relates generally to monitoring media and, more particularly, to methods and apparatus to determine locations of audience members.
BACKGROUND
Consuming media presentations generally involves listening to audio information and/or viewing video information such as, for example, radio programs, music, television programs, movies, still images, etc. Media-centric companies such as, for example, advertising companies, broadcasting networks, etc. are often interested in the viewing and listening interests of their audience to better market their products. A technique often used to measure the number and/or demographics of audience members exposed to media involves awarding media exposure credit to a media presentation each time an audience member is exposed to (e.g., is in the vicinity of) the media presentation.
The media exposure activities of audience members may be monitored using personal portable metering devices (PPMs), which are also known as portable meters, portable metering devices, and/or portable personal meters. A PPM is an electronic device that is typically worn (e.g., clipped to a belt or other apparel) or carried by an audience member. Some PPMs are configured to monitor media exposure by detecting or collecting information (e.g., program or source identifier codes, signatures, etc.) from audio and/or video signals that are emitted or presented by media delivery devices (e.g., televisions, stereos, speakers, computers, etc.).
An audience member or monitored individual performs their usual daily routine while wearing a PPM. The monitored individual's daily routine may include listening to the radio and/or other sources of audio media and/or watching television programs and/or other sources of audio/visual media. As the audience member is exposed to (e.g., views, listens to, is in the vicinity of, etc.) media, a PPM associated with (e.g., assigned to and carried by) that audience member generates monitoring data.
Unfortunately, the typical household presents unique monitoring challenges to the PPM. For example, a typical household includes multiple media delivery devices, each configured to deliver media content to specific viewing and/or listening areas located within the home. A PPM, carried by a person who is located in one of the viewing and/or listening areas, is configured to detect any media content being delivered in the viewing and/or listening area and to credit the programming associated with the media content as having been delivered to the corresponding audience member. Thus, the PPM operates on the premise that any media content detected by the PPM is associated with programming to which the person carrying the PPM was exposed. However, in some cases, a PPM may detect media content that is emitted by a media delivery device that is not located within the viewing or listening proximity of the person carrying the PPM, thereby causing the detected programming to be improperly credited.
The ability of the PPM to detect audio/video content being delivered outside of the viewing and/or listening proximity of the person carrying the PPM is an effect referred to as “spillover” because the media content being delivered outside of the viewing and/or listening proximity of the person carrying the PPM is described as “spilling over” into the area occupied by the person carrying the PPM. Spillover may occur, for example, in a case where a monitored individual in a bedroom is reading a book, but their PPM detects audio/video content delivered by a television in an adjacent living room (i.e., outside of their viewing/listening proximity), thereby causing the person carrying the PPM to be improperly credited as a member of the audience for the audio/video content.
Another effect, referred to as “hijacking” occurs when a person's PPM detects audio/video content being emitted from multiple media delivery devices at the same time. For example, an adult watching a television news program in a household kitchen may be located near a household family room in which children are watching a television cartoon program on a different television. Yet, the cartoon programming delivered by the family room television may, in some cases, have signals that overpower or “hijack” the signals associated with the news programming being emitted by the kitchen television. As a result, information collected by the adult's PPM may lead to inaccurately crediting the cartoon program as having been viewed by the adult and failing to credit the news program with any viewing. Still further, other common difficulties such as varying volume levels, varying audio/video content type (e.g., sparse, medium, rich, etc.), varying household transmission characteristics due to open/closed doors, movement and/or placement of furniture, acoustic characteristics of room layouts, wall construction, floor coverings, ceiling heights, etc. may lead to inaccurate audio/video content exposure detection by PPMs.



BRIEF DESCRIPTION OF THE DRAWINGS
FIG. 1 depicts an example monitored environment in which locations of audience members may be detected.
FIG. 2 depicts an example apparatus constructed in accordance with the teachings of this disclosure and that may be used to detect locations of audience members in the monitored environment of FIG. 1 based on time-domain analyses.
FIG. 3 depicts another example apparatus constructed in accordance with the teachings of this disclosure and that may be used to detect locations of audience members in the monitored environment of FIG. 1 based on frequency-domain analyses.
FIGS. 4A-4C depict reference and non-stationary audio sample sets at different time-shift positions relative to one another.
FIG. 5 depicts an example correlator of the example apparatus of FIG. 3 determining correlation coefficients based on frequency-domain representations of the reference and non-stationary audio sample sets of FIG. 3.
FIGS. 6A and 6B depict a flow diagram representative of example machine readable instructions that may be executed to determine locations of audience members.
FIG. 7 depicts a flow diagram representative of example machine readable instructions that may be executed to determine correlation coefficients for the process of FIG. 6A based on time-domain analyses when a DC bias is not present in any audio sample set being analyzed.
FIGS. 8A and 8B depict a flow diagram representative of example machine readable instructions that may be executed to determine correlation coefficients for the process of FIG. 6A based on time-domain analyses when a DC bias is present in an audio sample set being analyzed.
FIG. 9 depicts a flow diagram representative of example machine readable instructions that may be executed to determine correlation coefficients for the process of FIG. 6B based on frequency-domain analyses when a DC bias is not present in any audio sample set being analyzed.
FIGS. 10A-10C depict a flow diagram representative of example machine readable instructions that may be executed to determine correlation coefficients for the process of FIG. 6B based on frequency-domain analyses when a DC bias is present in an audio sample set being analyzed.
FIGS. 11A-11C depict a flow diagram representative of example machine readable instructions that may be executed to determine correlation coefficients for the process of FIG. 6B based on frequency-domain analyses when a DC bias is present in both audio sample sets being analyzed or when it is unknown whether a DC bias is present in one or both audio sample sets being analyzed.
FIG. 12 is an example processor system that can be used to execute the example instructions of FIGS. 6A, 6B, 7, 8A, 8B, 9, 10A-10C, and 11A-11C to implement the example apparatus and/or systems of FIGS. 1, 2, and 3 described herein.



DETAILED DESCRIPTION
Example methods, apparatus, systems, and articles of manufacture disclosed herein determine locations of audience members in monitored environments. Such audience member location information may be advantageously used to determine whether audience members were located sufficiently proximate to media devices (e.g., televisions, radios, computers, etc.) to have been exposed to media content presented by those media devices.
Example methods, apparatus, systems, and articles of manufacture disclosed herein involve using neural network models to determine locations of audience members. To implement location detection techniques, an example location detection system is installed at a home or other monitored environment (e.g., an office, a commercial establishment, a restaurant, a bar, a retail store, etc.) and calibrated for that environment. The location detection system includes a reference unit that operates as a stationary audio detector and a portable unit that operates as a portable audio detector. The reference unit includes a microphone and is positioned at a fixed location in the monitored environment. The portable unit also includes a microphone and is worn or carried by an audience member as the audience member moves in the monitored environment. During a calibration process, a neural network model of the location detection system is trained based on known distances between the reference unit and one or more audio emitters (e.g., speakers) of one or more media devices in the monitored environment and known distances between the portable unit and the audio emitter(s). The known distances and amplitude characteristics of audio detected at the reference and portable units are input to the neural network model.
After calibration of the location detection system, and during operation, an audience member carries or wears the portable unit. The portable unit digitizes any detected audio and communicates the digitized audio samples to the location detection system, while the reference unit also detects and digitizes the same audio. The location detection system then correlates the portable unit audio samples (i.e., non-stationary audio samples) with the reference unit audio samples (i.e., reference audio samples) to find one or more correlation coefficients (e.g., maximum and/or minimum correlation coefficients). The correlation coefficients are then provided to the neural network model, which determines a distance from the portable unit to the source of the detected audio (e.g., a speaker of a media device that emitted the audio) based on its previous training.
The resulting distance measure provided by the neural network model can then be used to determine whether the audience member carrying or wearing the portable unit was sufficiently close to the media device that emitted the audio to credit corresponding media content (e.g., a television program, a radio program, an advertisement, Internet web content, computer content, or any other media content) with exposure to the audience member.
In some examples, the location detection system is configured to process acquired audio samples from the reference unit and acquired audio samples from the portable unit using time-domain analysis techniques to determine correlation coefficients. In other examples the location detection system is configured to use frequency-domain analysis to process the acquired audio samples to determine the correlation coefficients.
FIG. 1 depicts an example monitored environment 100 having an example location detection system 102 constructed in accordance with the teachings of this disclosure located therein. In the illustrated example, the location detection system 102 includes a processing unit 104 and a stationary reference unit 106. In the illustrated example, the location detection system 102 also includes or is in communication with a portable meter 108 carried or worn by a person 110 as the person 110 moves in the monitored environment 100. In the illustrated example, the location detection system 102 is configured to determine the location(s) of the person 110 relative to a media presentation device 112 in the monitored environment 100. In the illustrated example, the media presentation device 112 is a television, but it may alternatively be a stereo, a computer, a multi-media device, or any other media presentation device capable of presenting media content sought to be monitored for exposure to the person 110.
In the illustrated example, the stationary reference unit 106 is located next to or substantially close to the media presentation device 112 such that distances between the portable meter 108 and the media presentation device 112 are the same or substantially similar to distances between the portable meter 108 and the stationary reference unit 106. In this manner, distances determined by the location detection system 102 between the portable meter 108 and the stationary reference unit 106 are regarded as the same distances between the portable meter 108 and the media presentation device 112.
In the illustrated example, the stationary reference unit 106 operates as a stationary audio detector to detect audio emitted from the media presentation device 112 of the monitored environment 100. In the illustrated example, the stationary reference unit 106 digitizes the detected audio and stores time-stamped digitized audio samples of the detected audio. In other examples, the stationary reference unit 106 sends the detected audio in analog form to the processing unit 104 and the processing unit 104 digitizes and stores time-stamped digitized audio samples of the detected audio. In some examples, the stationary reference unit 106 is provided with a microphone to detect the audio emitted by the media presentation device 112. In other examples, the stationary reference unit 106 is provided with audio input port connectors to interface with audio output ports or speaker output connections of the media presentation device 112 via wires to receive audio reproduced by the media presentation device 112.
In the illustrated example, the portable meter 108 operates as a portable audio detector that also detects audio emitted from the media presentation device 112 in the monitored environment 100. The example portable meter 108 digitizes the detected audio and stores time-stamped digitized audio samples of the detected audio.
In the illustrated example, the processing unit 104 may be configured to determine monitored distances (dM) between the person 110 and the media presentation device 112 in real-time or near real-time as the processing unit 104 receives audio samples from the stationary reference unit 106 and the portable meter 108. Additionally or alternatively, the processing unit 104 may be configured to determine the monitored distances (dM) during a post-process.
In the illustrated example, the example processing unit 104 determines locations of the person 110 by calculating monitored distances (dM) between the person 110 and the media presentation device 112 based on time-of-flight delays of audio emissions from the media presentation device 112. That is, the time-of-flight of audio emissions characterized by the speed of sound will cause an audio emission from the media presentation device 112 to be received at different times by the stationary reference unit 106 and the portable meter 108 when they are located at different distances from the media presentation device 112. In the illustrated example, the stationary reference unit 106 and the portable meter 108 timestamp their respective digitized audio samples to indicate times at which they detected corresponding audio. The processing unit 104 subsequently uses such timestamps to determine delays between the times at which the stationary reference unit 106 received audio emissions and the times at which the portable meter 108 received the same audio emissions. The processing unit 104 can then use the measured delays to determine monitored distances (dM) corresponding to the locations of the person 110.
The stationary reference unit 106 and the portable meter 108 periodically or aperiodically send their time-stamped digitized audio samples to the processing unit 104. The processing unit 104 of the illustrated example may be implemented in a number of different ways. For example, the processing unit 104 may be implemented by the example apparatus 200 of FIG. 2 to determine locations of the person 110 using time-domain analyses of the time-stamped digitized audio samples and/or by the example apparatus 300 of FIG. 3 to determine locations of the person 110 using frequency-domain analyses of the time-stamped digitized audio samples.
In the illustrated example, the processing unit 104 includes a neural network model. Example neural network models are described below in connection with FIGS. 2 and 3. During installation of the location detection system 102 in the monitored environment 100 and prior to use, the processing unit 104 is calibrated with training data associated with known training distances to train the neural network model to determine location information associated with the person 110. In particular, during a calibration process, media content is presented via the media presentation device 112 to cause audio emissions therefrom. Throughout the collection process, the stationary reference unit 106 detects the audio emissions and generates corresponding time-stamped digitized reference audio samples. In addition, the portable meter 108 is moved to different locations in the monitored environment 100 shown in FIG. 1 as located at known training distances (dT1), (dT2), and (dT3) from the media presentation device 112 (and, thus, from the stationary reference unit 106). At each of the training distances (dT1), (dT2), and (dT3), the portable meter 108 detects the audio emissions from the media presentation device 112 and generates corresponding time-stamped digitized non-stationary audio samples.
In the illustrated example, the stationary reference unit 106 and the mobile meter 108 send their respective time-stamped digitized audio samples to the processing unit 104 through wireless or wired mediums. The processing unit 104 is also provided with the training distances (dT1), (dT2), and (dT3) and indications of times at which the portable meter 108 was located at the training distances (dT1), (dT2), and (dT3). To train the neural network model to accurately determine locations of the person 110 based on subsequent audio collection processes, the processing unit 104 provides its neural network model with a distance value for each training distance (dT1), (dT2), and (dT3) and processes audio sample sets from the stationary reference unit 106 and the portable meter 108 (e.g., digitized reference and non-stationary audio sample sets collected when the portable meter 108 was located at the respective training distances (dT1), (dT2), and (dT3)) to provide the neural network with correlation coefficients for each of the training distances (dT1), (dT2), and (dT3)). In this manner, the neural network model can be trained to correlate the training distances (dT1), (dT2), and (dT3) with the correlation coefficient results associated with the digitized audio samples to subsequently determine accurate location information when the person 110 is located at different distances from the media presentation device 112.
In the illustrated example, location information refers to locations of the person 110 relative to the media presentation device 112. That is, the processing unit 104 determines distances separating the person 110 from the media presentation device 112 such that the person 110 may be located at any position along a circumference about the media presentation device 112 having a radius equal to the determined distance. In the illustrated example of FIG. 1, when the processing unit 104 subsequently determines that the person 110 was located at a monitored distance (dM) from the media presentation device 112, the person 110 may have been located anywhere along a circumference 116 having a radius equal to the monitored distance (dM) about the media presentation device 112.
FIG. 2 depicts an example apparatus 200 that may be used to implement the example processing unit 104 of FIG. 1 to detect locations of audience members in the monitored environment 100 of FIG. 1 based on time-domain signal analyses. In the illustrated example, the apparatus 200 is in communication with the stationary reference unit 106 to receive digitized reference audio samples 202 and in communication with the portable meter (PM) 108 to receive digitized non-stationary audio samples 204. To synchronize clocks in the stationary reference unit 106 and the portable meter 108 with one another, the apparatus 200 is provided with an example time synchronizer 206. In the illustrated example, the time synchronizer 206 configures the stationary reference unit 106 and the portable meter 108 to have the same time so that timestamps of the reference audio samples 202 from the stationary reference unit 106 correspond to the same timestamps (or substantially the same timestamps within an acceptable error tolerance (e.g., +/−number of milliseconds)) of the non-stationary audio samples 204 from the portable meter 108. Ensuring that reference audio samples 202 from the stationary reference unit 106 correspond in time with respective ones of the non-stationary audio samples 204 from the portable meter 108 facilitates correlation processes performed on the audio samples 202 and 204 by the apparatus 200. In the illustrated example, the time synchronizer 206 can communicate with the stationary reference unit 106 and the portable meter 108 through wired or wireless communications to set clocks therein with the same (or substantially the same) time.
In the illustrated example, to sub-sample (or down-sample) the reference audio samples 202 from the stationary reference unit 106, the apparatus 200 is provided with a sub-sampler 208. In addition, to sub-sample the non-stationary audio samples 204 from the portable meter 108, the apparatus 200 is provided with a second sub-sampler 210. In the illustrated example, the sub-samplers 208 and 210 reduce the effective sampling rate at which the received digitized audio samples are represented. For example, if the stationary reference unit 106 and the portable meter 108 sample audio at 44 kHz, the sub-samplers 208 and 210 may select every other sample of the received digitized audio samples, which results in an effective sampling rate of 22 kHz. Sub-sampling reduces the size of the audio data sets to be processed by the apparatus 200, while maintaining the duration or the period over which the samples were collected and maintaining sufficient audio quality for determining audience member location information. In addition, the sub-sampling advantageously allows the apparatus 200 to process the digitized audio samples faster by reducing the amount of data that must be processed. In the illustrated example, the sub-samplers 208 and 210 perform 20:1 sub-sampling processes to reduce digitized audio sample sets by a factor of 20.
In the illustrated example, to reduce the quantization or digitization resolution of the non-stationary audio samples 204 from the portable meter 108, the apparatus 200 is provided with an example down-scaler 212. In the illustrated example, the down-scaler 212 is configured to reduce the number of bits used to represent each audio sample of the non-stationary audio samples 204. In the illustrated example, the down-scaler 212 down-scales audio samples from 16-bit representations to 8-bit representations. The example apparatus 200 is provided with the down-scaler 212 to normalize the bit-resolutions of the non-stationary audio samples 204 from the portable meter 108 to match the bit-resolutions of the reference audio samples 202 from the stationary reference unit 106, because the stationary reference unit 106 of the illustrated example digitizes audio at the lower, 8-bit quantization resolution. In other examples in which the portable meter 108 and the stationary reference unit 106 digitize audio using the same bit resolution, the down-scaler 212 is omitted. In yet other examples in which the stationary reference unit 106 digitizes audio using a higher bit resolution than the portable meter 108, the down-scaler 212 is used to down-scale the reference audio samples 202 from the stationary reference unit 106, rather than the non-stationary audio samples 204 from the portable meter 108.
In the illustrated example, to perform correlations between the audio samples 202 and 204, the apparatus 200 is provided with an example correlator 214. In the illustrated example, the correlator 214 determines correlation coefficients (e.g., correlation analysis results) based on correlations between audio sample subsets 402 and 404 (FIGS. 4A-4C) generated by the sub-samplers 208 and 210 and the down-scaler 212 based on the reference audio samples 202 and the non-stationary audio samples 204. The correlator 214 of the illustrated example receives a subset of reference audio samples 402 from the sub-sampler 208 and a subset of non-stationary audio samples 404 from the down-scaler 212.
Turning to FIGS. 4A-4C, the correlator 214 aligns the reference audio sample subset 402 with the non-stationary audio sample subset 404 based on respective timestamps 406 and 408. As shown in FIG. 4A, the time-aligned subsets 402 and 404 exhibit corresponding similar characteristics 410 that are offset from one another by a measured distance time offset (tdist). In some examples described herein, the measured distance time offset (tdist) corresponds to the monitored distance (dM) between the person 110 and the media presentation device 112 shown in FIG. 1.
To determine the measured distance time offset (tdist), the correlator 214 determines multiple correlation coefficients by performing correlations between the reference and non-stationary audio sample subsets 402 and 404. For each correlation coefficient calculation, the correlator 214 time-shifts the audio sample subsets 402 and 404 relative to one another by incremental time shifts to different time-shift positions (tshift) as shown in FIGS. 4A, 4B, and 4C such that at one of the time-shift positions (tshift) (e.g., the time shift position (tshift) shown in FIG. 4B) the similar characteristics 410 between the audio sample subsets 402 and 404 are close in alignment or substantially aligned with one another to produce a very strong correlation coefficient. In the illustrated examples, the strongest correlation between the audio sample subsets 402 and 404 occurs at the time shift position equal to 2.5 milliseconds (ms) (tshift=2.5 ms) of FIG. 4B as represented by the notation below the timeline of the shown graph.
In some examples, the correlator 214 initially aligns the reference audio sample subset 402 to be negatively shifted in time relative to the non-stationary audio sample subset 404. The correlator 214 then performs a time-shift sweep of the audio sample subsets 402 and 404 relative to one another when calculating the correlation coefficients such that the timestamps 406 and 408 of the audio sample subsets 402 and 404 become matchingly aligned at a subsequent time-shift and, at yet a further time-shift, the reference audio sample subset 402 becomes positively shifted in time relative to the non-stationary audio sample subset 404. For each time-shift, the correlator 214 associates a corresponding correlation coefficient with the corresponding time-shift position (tshift) between the timestamps 406 and 408 of correspondingly aligned audio samples. For the time-shift at which the similar characteristics 408 are aligned with one another, as indicated by the very strongest correlation coefficient, the measured distance time offset (tdist) (e.g., the measured distance time offset (tdist) shown in FIG. 4A) is substantially equal to the time-shift position (tshift) of that time shift.
Although only three time-shifts are shown in FIGS. 4A-4C, in the illustrated example, the correlator 214 determines correlation coefficients for the audio sample subsets 402 and 404 based on more time-shifts. In addition, the time-shifts may be based on finer or coarser time steps (e.g., 0.2 ms increments, 0.5 ms increments, 1 ms increments, 2 ms increments, 5 ms increments, etc.).
In the illustrated example, to determine correlation coefficients based on time-domain analyses, the correlator 214 may elect to use one of Equation 1 or Equation 2 below based on whether the average of one or both of the audio samples 202 and 204 is equal to zero (i.e., whether the audio samples 202 and 204 have DC biases). In the illustrated example, the average of an audio sample set is equal to zero when it does not have a direct current (DC) bias. Such a DC bias may be eliminated from an audio signal using a high-pass filter or any other suitable filtering process such as, for example, a moving average process. For example, one or both of the stationary reference unit 106 and/or the portable meter 108 may be provided with a respective high-pass filter (or any other suitable filter) to remove DC bias from a detected audio signal prior to digitizing the signal to produce the audio samples 202 and/or 204.
In the illustrated example, the correlator 214 uses Equation 1 below to perform a correlation between the reference and non-stationary audio samples 202 and 204 and determine a correlation coefficient (r) when the averages of both of the audio samples 202 and 204 are equal to zero (i.e., neither of the audio samples 202 and 204 contain a DC bias). In such instances, the stationary reference unit 106 and the portable meter 108 used high-pass filters (or other suitable filters to remove DC biases) to filter audio signals corresponding to the reference and non-stationary audio samples 202 and 204.




  
    
      
        r
        =
        
          
            ∑
            xy
          
          
            
              
                (
                
                  ∑
                  
                    x
                    2
                  
                
                )
              
              ⁢
              
                (
                
                  ∑
                  
                    y
                    2
                  
                
                )
              
            
          
        
      
    
    
      
        Equation
        ⁢
        
          
        
        ⁢
        1
      
    
  




Using Equation 1 above, the correlator 214 determines a correlation coefficient (r) (e.g., a correlation analysis result) by performing three summations (Σxy, Σx2, Σy2) based on the reference and non-stationary audio samples 202 and 204. An example process that may be used by the correlator 214 to determine correlation coefficients (r) based on Equation 1 is described below in connection with the example flow diagram of FIG. 7.

In the illustrated example, the correlator 214 uses Equation 2 below to perform a correlation between the reference and non-stationary audio samples 202 and 204 and determine a correlation coefficient (r) when the average of the reference audio samples 202 is zero (i.e., the reference audio samples 202 have no DC bias) and the average of the non-stationary audio samples 204 is not equal to zero (i.e., the non-stationary audio samples 204 have a DC bias). In such instances, the stationary reference unit 106 used a high-pass filter (or other suitable filter to remove a DC bias) to filter an audio signal corresponding to the reference audio samples 202.




  
    
      
        r
        =
        
          
            n
            ⁢
            
              ∑
              xy
            
          
          
            
              
                n
                ⁡
                
                  (
                  
                    ∑
                    
                      x
                      2
                    
                  
                  )
                
              
            
            ⁢
            
              
                
                  n
                  ⁡
                  
                    (
                    
                      ∑
                      
                        y
                        2
                      
                    
                    )
                  
                
                -
                
                  
                    (
                    
                      ∑
                      y
                    
                    )
                  
                  2
                
              
            
          
        
      
    
    
      
        Equation
        ⁢
        
          
        
        ⁢
        2
      
    
  




Using Equation 2 above, the correlator 214 determines a correlation coefficient (r) (e.g., a correlation analysis result) using four summations (Σxy, Σx2, Σy2, Σy) based on the reference and non-stationary audio samples 202 and 204. An example process that may be used by the correlator 214 to determine correlation coefficients (r) based on Equation 2 is described below in connection with the example flow diagram of FIGS. 8A and 8B.

In the illustrated example, Equations 1 and 2 may be derived from known Equation 3 below when the respective conditions noted above regarding DC biases of the reference and non-stationary audio samples 202 and 204 are true.




  
    
      
        r
        =
        
          
            
              n
              ⁢
              
                ∑
                xy
              
            
            -
            
              
                (
                
                  ∑
                  x
                
                )
              
              ⁢
              
                (
                
                  ∑
                  y
                
                )
              
            
          
          
            
              
                
                  n
                  ⁡
                  
                    (
                    
                      ∑
                      
                        x
                        2
                      
                    
                    )
                  
                
                -
                
                  
                    (
                    
                      ∑
                      x
                    
                    )
                  
                  2
                
              
            
            ⁢
            
              
                
                  n
                  ⁡
                  
                    (
                    
                      ∑
                      
                        y
                        2
                      
                    
                    )
                  
                
                -
                
                  
                    (
                    
                      ∑
                      y
                    
                    )
                  
                  2
                
              
            
          
        
      
    
    
      
        Equation
        ⁢
        
          
        
        ⁢
        3
      
    
  




In some examples, Equation 3 may be used to determine correlation coefficients (r) (e.g., a correlation analysis results) when it is unknown whether a DC bias is present in either of the audio samples 202 and 204 and/or when DC biases are present in both of the audio samples 202 and 204.

As shown in FIG. 2, to analyze and select ones of the correlation coefficients (r) determined using Equation 1, Equation 2, and/or Equation 3 above, the apparatus 200 is provided with a correlation analyzer 216. In the illustrated example, the correlation analyzer 216 receives correlation coefficients (r) and corresponding time-shift positions (tshift) determined by the correlator 214 based on multiple time shifts of the reference and non-stationary audio sample subsets 402 and 404 of FIGS. 4A-4C. In the illustrated example, the correlation analyzer 216 compares the correlation coefficients (r) to one another and selects the three maximum or peak coefficients (rmax) and three minimum or valley coefficients (rmin). In the illustrated example, the three maximum coefficients (rmax) are coefficients that have the largest values relative to the other correlation coefficients (r) for corresponding reference and non-stationary audio sample subsets 402 and 404, and the three minimum coefficients (rmin) are coefficients that have the smallest values relative to the other correlation coefficients (r) for corresponding reference and non-stationary audio sample subsets 402 and 404. In other examples, fewer or more maximum coefficients (rmax) and minimum coefficients (rmin) may be selected.
In the illustrated example of FIG. 2, to select a maximum-amplitude audio sample from the reference audio samples 202, the apparatus 200 is provided with an example maximum level detector 218. The example maximum level detector 218 selects an audio sample having the largest amplitude from the reference audio samples 202 over a time period of interest and stores a value corresponding to its amplitude. In addition, to select a maximum-amplitude audio sample from the non-stationary audio samples 204, the apparatus 200 is provided with another example maximum level detector 220. The example maximum level detector 220 selects an audio sample having the largest amplitude from the non-stationary audio samples 204 over a time period of interest and stores a value corresponding to its amplitude.
In the illustrated example of FIG. 2, to determine the monitored distance (dM) of the person 110 from the media presentation device 112 of FIG. 1, the apparatus 200 is provided with a neural network 222. The neural network 222 receives the three maximum coefficients (rmax), the three minimum coefficients (rmin), and corresponding time-shift positions (tshift) from the correlation analyzer 216. In addition, the neural network 222 receives a maximum-amplitude reference audio sample value from the maximum level detector 218 and a maximum-amplitude non-stationary audio sample value from the maximum level detector 220. In the illustrated example, the neural network 222 implements a neural network model that is previously trained based on the known training distances (dT1), (dT2), and (dT3) as discussed above in connection with FIG. 1 to generate distance measures such as the monitored distance (dM). Thus, in the illustrated example, the neural network 222 applies its previously trained neural network model to the three maximum coefficients (rmax), the three minimum coefficients (rmin), corresponding time-shift positions (tshift), the maximum-amplitude reference audio sample value, and the maximum-amplitude non-stationary audio sample value to determine a distance result 224 (e.g., the monitored distance (dM) of FIG. 1).
FIG. 3 depicts another example apparatus 300 that may be used to implement the example processing unit 104 to detect locations of audience members in the monitored environment 100 of FIG. 1 based on frequency-domain analyses. In the illustrated example, the apparatus 300 is in communication with the stationary reference unit 106 to receive the digitized reference audio samples 202 and in communication with the portable meter (PM) 108 to receive the digitized non-stationary audio samples 204. To synchronize clocks in the stationary reference unit 106 and the portable meter 108 with one another as discussed above in connection with FIG. 2, the apparatus 300 is provided with an example time synchronizer 306, which is substantially similar or identical to the time synchronizer 206 of FIG. 2.
In the illustrated example, to convert the reference audio samples 202 from the time domain to the frequency domain, the apparatus 300 is provided with an example frequency-domain transformer 308. In addition, to convert the non-stationary audio samples 204 to the frequency domain, the apparatus 300 is provided with another example frequency-domain transformer 310. In the illustrated example, the frequency-domain transformers 308 and 310 determine fast Fourier transforms (FFTs) of audio samples to enable frequency domain analyses of the same.
In the illustrated example, the apparatus 300 is not provided with sub-samplers (e.g., the sub-samplers 208 and 210 of FIG. 2). The example apparatus 300 can perform frequency domain analyses based on a full audio sample sets collected during a shorter duration by the stationary reference unit 106 and the portable unit 108, whereas the apparatus 200 performs time-domain analyses on audio sample sets collected during longer durations and uses the sub-samplers 208 and 210 to reduce the size of the audio sample sets. In the illustrated examples disclosed herein, the quantity of audio samples collected by each of the stationary reference unit 106 and the portable meter 108 is based on the desired size of a correlation window and time-shift range for determining correlation coefficients (r). In examples using the example apparatus 200 of FIG. 2 to perform time-domain analyses, each of the stationary reference unit 106 and the portable meter 108 may collect 750 ms of audio samples plus an additional quantity of audio samples equal to twice the duration of the desired time-shift range. In such examples, the sub-samplers 208 and 210 are then used to reduce quantities of samples while maintaining the duration over which the audio samples were acquired. In other examples using the example apparatus 300 of FIG. 3 to perform frequency-domain analyses, each of the stationary reference unit 106 and the portable meter 108 may collect 32 ms of audio samples plus an additional quantity of audio samples equal to twice the duration of the desired time-shift range. The example apparatus 300 does not employ sub-sampling of the acquired audio samples for frequency-domain analyses, because the number of audio samples acquired during the relatively shorter duration is already of a smaller, more manageable quantity compared to the number of audio samples collected by the stationary reference unit 106 and the portable meter 108 for time-domain analyses. Frequency-domain analyses may be advantageously used to improve battery charge performance of the portable meter 108, because less audio samples are acquired for frequency-domain analyses to achieve similar results as those produced using time-domain analyses.
Referring again to FIG. 3, to perform correlations between frequency-domain representations of the audio samples 202 and 204, the apparatus 300 is provided with an example correlator 314. In the illustrated example, the frequency-domain transformers 308 and 310 and the correlator 314 work cooperatively to frequency-correlate multiple time-shifted instances of the reference and non-stationary audio samples 202 and 204. In this manner, the correlator 314 may determine multiple correlation coefficients, each corresponding to a different time-shifted instance of the audio samples 202 and 204. For example, turning briefly to FIG. 5, the frequency domain transformer 308 is shown as generating multiple reference frequency transforms 502, each corresponding to a respective one of a 0 ms time-shift position (tshift=0 ms), a 2.5 ms time-shift position (tshift=2.5 ms), and a 3.5 ms time-shift position (tshift=3.5 ms). In addition, the frequency-domain transformer 310 is shown as generating multiple non-stationary frequency transforms 504, each corresponding to a respective one of the 0 ms time-shift position (tshift=0 ms), the 2.5 ms time-shift position (tshift=2.5 ms), and the 3.5 ms time-shift position (tshift=3.5 ms). For the frequency-domain transforms performed by the frequency-domain transformers 308 and 310, the time shifting is performed in similar fashion as described above in connection with FIGS. 4A-4C. Thus, the time shift positions (tshift) of the audio samples 202 and 204 are substantially similar or identical to the time shift positions shown in FIGS. 4A-4C. Although three time shifts are noted (tshift=0 ms, tshift=2.5 ms, and tshift=3.5 ms), the correlator 314 may determine correlation coefficients (r) for any other quantity of time shifts.
In the illustrated example of FIG. 5, for each time-shift position (tshift), the correlator 314 determines frequency-based correlation coefficients (r) 506 for the multiple frequency-domain representations 502 and 504. In the illustrated example, the correlator 314 may use Equation 4 or Equation 5 below to determine the frequency-based correlation coefficients (r) 506. In particular, the correlator 314 may elect to use one of Equation 4 or Equation 5 below based on whether the average of one or both of the audio samples 202 and 204 is equal to zero (i.e., whether either or both of the audio samples 202 and 204 have DC biases). In the illustrated example, the average of an audio sample set is equal to zero when it does not have a direct current (DC) bias. Such a DC bias may be eliminated from an audio signal using a high-pass filter or any other suitable filtering process such as, for example, a moving average process. For example, one or both of the stationary reference unit 106 and/or the portable meter 108 may be provided with a respective high-pass filter (or any other suitable filter) to remove DC bias from a detected audio signal prior to digitizing the signal to produce the audio samples 202 and/or 204.
In the illustrated example, the correlator 314 uses Equation 4 below to perform frequency-based correlations between different time-shift position (tshift) instances of the reference and non-stationary audio samples 202, 204 and to determine respective correlation coefficients (r) when the averages of both of the audio samples 202 and 204 are equal to zero (i.e., neither of the audio samples 202 and 204 contains a DC bias). In such instances, the stationary reference unit 106 and the portable meter 108 used high-pass filters (or other suitable filters to remove DC biases) to filter audio signals corresponding to the reference and non-stationary audio samples 202 and 204.




  
    
      
        r
        =
        
          
            ∑
            
              XY
              *
            
          
          
            
              
                (
                
                  ∑
                  
                    XX
                    *
                  
                
                )
              
              ⁢
              
                (
                
                  ∑
                  
                    YY
                    *
                  
                
                )
              
            
          
        
      
    
    
      
        Equation
        ⁢
        
          
        
        ⁢
        4
      
    
  




Using Equation 4 above, the correlator 314 determines a correlation coefficient (r) (e.g., a correlation analysis result) by performing three summations (ΣXY*, ΣXX*, ΣYY*) based on the complex conjugate pairs (X,X*) and (Y,Y*) of frequency-domain representations 502 and 504 (FIG. 5) of the reference and non-stationary audio samples 202 and 204. An example process that may be used by the correlator 314 of FIG. 3 to determine correlation coefficients (r) based on Equation 4 is described below in connection with the example flow diagram of FIG. 9.

In the illustrated example, the correlator 314 uses Equation 5 below to perform frequency-based correlations between different time-shift position (tshift) instances of the reference and non-stationary audio samples 202, 204 and to determine respective correlation coefficients (r) when the average of the reference audio samples 202 is zero (i.e., the reference audio samples 202 have no DC bias), but the average of the non-stationary audio samples 204 is not equal to zero (i.e., the non-stationary audio samples 204 have a DC bias). In such instances, the stationary reference unit 106 used a high-pass filter (or other suitable filter to remove a DC bias) to filter an audio signal corresponding to the reference audio samples 202.




  
    
      
        r
        =
        
          
            n
            ⁢
            
              ∑
              
                XY
                *
              
            
          
          
            
              
                n
                ⁡
                
                  (
                  
                    ∑
                    
                      XX
                      *
                    
                  
                  )
                
              
            
            ⁢
            
              
                
                  n
                  ⁡
                  
                    (
                    
                      ∑
                      
                        YY
                        *
                      
                    
                    )
                  
                
                -
                
                  
                    (
                    
                      ∑
                      Y
                    
                    )
                  
                  ⁢
                  
                    
                      (
                      
                        ∑
                        Y
                      
                      )
                    
                    *
                  
                
              
            
          
        
      
    
    
      
        Equation
        ⁢
        
          
        
        ⁢
        5
      
    
  




Using Equation 5 above, the correlator 314 determines a correlation coefficient (r) (e.g., a correlation analysis result) using four summations (ΣXY*, ΣXX*, ΣYY, ΣY) based on the complex conjugate pairs (X,X*) and (Y,Y*) of frequency-domain representations 502 and 504 (FIG. 5) of the reference and non-stationary audio samples 202 and 204. An example process that may be used by the correlator 314 of FIG. 3 to determine correlation coefficients (r) based on Equation 5 is described below in connection with the example flow diagram of FIGS. 10A-10C.

In the illustrated example, Equations 4 and 5 may be derived from Equation 6 below when the respective conditions noted above regarding DC biases of the reference and non-stationary audio samples 202 and 204 are true.




  
    
      
        r
        =
        
          
            
              n
              ⁢
              
                ∑
                
                  (
                  
                    XY
                    *
                  
                  )
                
              
            
            -
            
              
                (
                
                  ∑
                  X
                
                )
              
              ⁢
              
                (
                
                  ∑
                  
                    Y
                    *
                  
                
                )
              
            
          
          
            
              
                
                  n
                  ⁡
                  
                    (
                    
                      ∑
                      
                        XX
                        *
                      
                    
                    )
                  
                
                -
                
                  
                    (
                    
                      ∑
                      X
                    
                    )
                  
                  ⁢
                  
                    
                      (
                      
                        ∑
                        X
                      
                      )
                    
                    *
                  
                
              
            
            ⁢
            
              
                
                  n
                  ⁡
                  
                    (
                    
                      ∑
                      
                        
                      
                      ⁢
                      
                        YY
                        *
                      
                    
                    )
                  
                
                -
                
                  
                    (
                    
                      ∑
                      
                        
                      
                      ⁢
                      Y
                    
                    )
                  
                  ⁢
                  
                    
                      (
                      
                        ∑
                        Y
                      
                      )
                    
                    *
                  
                
              
            
          
        
      
    
    
      
        Equation
        ⁢
        
          
        
        ⁢
        6
      
    
  




In some example implementations, Equation 6 may be used to determine correlation coefficients (r) (e.g., correlation analysis results) when it is unknown whether a DC bias is present in either of the audio samples 202 and 204 and/or when DC biases are present in both of the audio samples 202 and 204. An example process that may be used by the correlator 314 of FIG. 3 to determine correlation coefficients (r) based on Equation 6 is described below in connection with the example flow diagram of FIGS. 11A-11C.

For frequency-domain analyses using Equations 4-6 above, the correlation coefficients (r) may be represented in polar coordinate notation (i.e., r=ρ∠θ), in which the correlation coefficients (r) are expressed as a magnitude (ρ) at a phase difference (θ) between the two audio signals. In some examples in which phase distortion between audio signals is small (e.g., minimal) or can otherwise be ignored (e.g., there is no interest in taking phase distortion into account or no interest in analyzing phase delay), the correlation coefficients (r) may be represented using the magnitude (ρ) without the phase difference (θ).
As shown in FIG. 3, to analyze and select ones of the correlation coefficients (r) determined using Equation 4, Equation 5, and/or Equation 6 above, the apparatus 300 is provided with a correlation analyzer 316. In the illustrated example, the correlation analyzer 316 receives correlation coefficients (r) and corresponding time-shift positions (tshift) from the correlator 314 based on multiple time shifts of the reference and non-stationary audio samples 202 and 204 (e.g., similar to time-shift positions (tshift) of the reference and non-stationary audio sample subsets 402 and 404 of FIGS. 4A-4C). In the illustrated example, the correlation analyzer 316 compares the correlation coefficients (r) to one another and selects the six maximum or peak coefficients (rmax). In the illustrated example, the six maximum coefficients (rmax) are coefficients that have the largest values relative to the other correlation coefficients (r) for corresponding reference and non-stationary audio samples 202 and 204. In other examples, fewer or more maximum coefficients (rmax) may be selected.
In the illustrated example of FIG. 3, to select a maximum-amplitude audio sample from the reference audio samples 202, the apparatus 300 is provided with an example maximum level detector 318. The example maximum level detector 318 selects an audio sample having the largest amplitude from the reference audio samples 202 over a time period of interest and stores a value corresponding to its amplitude. In addition, to select a maximum-amplitude audio sample from the non-stationary audio samples 204, the apparatus 300 is provided with another example maximum level detector 320. The example maximum level detector 320 selects an audio sample having the largest amplitude from the non-stationary audio samples 204 over a time period of interest and stores a value corresponding to its amplitude.
In the illustrated example of FIG. 3, to determine the monitored distance (dM) of the person 110 of FIG. 1, the apparatus 300 is provided with a neural network 322. The neural network 322 receives the six maximum coefficients (rmax) and corresponding time-shift positions (tshift) from the correlation analyzer 316. In addition, the neural network 322 receives a maximum-amplitude reference audio sample value from the maximum level detector 318 and a maximum-amplitude non-stationary audio sample value from the maximum level detector 320. In the illustrated example, the neural network 322 implements a neural network model that is previously trained based on the known training distances (dT1), (dT2), and (dT3) as discussed above in connection with FIG. 1 to generate distance measures such as the monitored distance (dM). Thus, in the illustrated example, the neural network 322 applies its previously trained neural network model to the six maximum coefficients (rmax), corresponding time-shift positions (tshift), the maximum-amplitude reference audio sample value, and the maximum-amplitude non-stationary audio sample value to determine a distance result 324 (e.g., the monitored distance (dM) of FIG. 1).
In some examples, the apparatus 200 of FIG. 2 and the apparatus 300 of FIG. 3 are both implemented in the processing unit 104 of FIG. 1. In such examples, the processing unit 104 is configured to determine which of the apparatus 200 or the apparatus 300 it should use based on whether it is to perform a time-domain analysis or a frequency-domain analysis of received audio samples (e.g., the audio samples 202 and 204 of FIGS. 2 and 3). For example, user-defined configurations may be set in the processing unit 104 to select between performance of time-domain or frequency-domain analyses based on any number and/or types of factors. Such factors may include the acoustic characteristics of the monitored environment 100, the size of the monitored environment 100, the device type of the media presentation device 112 (e.g., a stereo, a computer, a television, etc.), the speaker type of the media presentation device 112 (e.g., mixed range speakers, high-fidelity speakers, low-cost speakers, etc.), the processing power of the processing unit 104, the processing power of the stationary reference unit 106, the processing power of the portable meter 108, the battery capacity or other battery characteristics of the portable unit 108, etc.
In other examples, the processing unit 104 may be provided with only one of the apparatus 200 or the apparatus 300. In such examples, the processing unit 104 is configured to perform only time-domain analyses or perform only frequency-domain analyses.
While example manners of implementing the apparatus 200 and 300 have been illustrated in FIGS. 2 and 3, one or more of the elements, processes and/or devices illustrated in FIGS. 2 and 3 may be combined, divided, re-arranged, omitted, eliminated and/or implemented in any other way. Further, the time synchronizer 206, the sub-samplers 208 and 210, the down-scaler 212, the correlation analyzer 216, the neural network 222 and/or, more generally, the example apparatus 200 of FIG. 2 may be implemented by hardware, software, firmware and/or any combination of hardware, software and/or firmware. In addition, the time synchronizer 306, the frequency-domain transformers 308 and 310, the correlator 314, the correlation analyzer 316, the neural network 322 and/or, more generally, the example apparatus 300 of FIG. 3 may be implemented by hardware, software, firmware and/or any combination of hardware, software and/or firmware. Thus, for example, any of the time synchronizer 206, the sub-samplers 208 and 210, the down-scaler 212, the correlation analyzer 216, the neural network 222, the example apparatus 200, the time synchronizer 306, the frequency-domain transformers 308 and 310, the correlator 314, the correlation analyzer 316, the neural network 322 and/or the example apparatus 300 could be implemented by one or more circuit(s), programmable processor(s), application specific integrated circuit(s) (ASIC(s)), programmable logic device(s) (PLD(s)) and/or field programmable logic device(s) (FPLD(s)), etc. When any of the appended apparatus claims are read to cover a purely software and/or firmware implementation, at least one of the time synchronizer 206, the sub-samplers 208 and 210, the down-scaler 212, the correlation analyzer 216, the neural network 222, the time synchronizer 306, the frequency-domain transformers 308 and 310, the correlator 314, the correlation analyzer 316, and/or the neural network 322 are hereby expressly defined to include a computer readable medium such as a memory, DVD, CD, etc. storing the software and/or firmware. Further still, the example apparatus 200 and 300 may include one or more elements, processes and/or devices in addition to, or instead of, those illustrated in FIGS. 2 and 3, and/or may include more than one of any or all of the illustrated elements, processes and devices.
FIGS. 6A, 6B, 7, 8A, 8B, 9, 10A-10C, and 11A-11C are flow diagrams representative of example machine readable instructions that can be executed to implement example methods, apparatus, and systems described herein. The example processes of FIGS. 6A, 6B, 7, 8A, 8B, 9, 10A-10C, and 11A-11C may be implemented using machine readable instructions that, when executed, cause a device (e.g., a programmable controller or other programmable machine or integrated circuit) to perform the operations shown in FIGS. 6A, 6B, 7, 8A, 8B, 9, 10A-10C, and 11A-11C. For instance, the example processes of FIGS. 6A, 6B, 7, 8A, 8B, 9, 10A-10C, and 11A-11C may be performed using a processor, a controller, and/or any other suitable processing device. For example, the example processes of FIGS. 6A, 6B, 7, 8A, 8B, 9, 10A-10C, and 11A-11C may be implemented using coded instructions stored on a tangible machine readable medium such as a flash memory, a read-only memory (ROM), and/or a random-access memory (RAM).
As used herein, the term tangible computer readable medium is expressly defined to include any type of computer readable storage and to exclude propagating signals. Additionally or alternatively, the example processes of FIGS. 6A, 6B, 7, 8A, 8B, 9, 10A-10C, and 11A-11C may be implemented using coded instructions (e.g., computer readable instructions) stored on a non-transitory computer readable medium such as a flash memory, a read-only memory (ROM), a random-access memory (RAM), a cache, or any other storage media in which information is stored for any duration (e.g., for extended time periods, permanently, brief instances, for temporarily buffering, and/or for caching of the information). As used herein, the term non-transitory computer readable medium is expressly defined to include any type of computer readable medium and to exclude propagating signals.
Alternatively, the example processes of FIGS. 6A, 6B, 7, 8A, 8B, 9, 10A-10C, and 11A-11C may be implemented using any combination(s) of application specific integrated circuit(s) (ASIC(s)), programmable logic device(s) (PLD(s)), field programmable logic device(s) (FPLD(s)), discrete logic, hardware, firmware, etc. Also, the example processes of FIGS. 6A, 6B, 7, 8A, 8B, 9, 10A-10C, and 11A-11C may be implemented as any combination(s) of any of the foregoing techniques, for example, any combination of firmware, software, discrete logic and/or hardware.
Although the example processes of FIGS. 6A, 6B, 7, 8A, 8B, 9, 10A-10C, and 11A-11C are described with reference to the flow diagrams of FIGS. 6A, 6B, 7, 8A, 8B, 9, 10A-10C, and 11A-11C, other methods of implementing the processes of FIGS. 6A, 6B, 7, 8A, 8B, 9, 10A-10C, and 11A-11C may be employed. For example, the order of execution of the blocks may be changed, and/or some of the blocks described may be changed, eliminated, sub-divided, or combined. Additionally, one or both of the example processes of FIGS. 6A, 6B, 7, 8A, 8B, 9, 10A-10C, and 11A-11C may be performed sequentially and/or in parallel by, for example, separate processing threads, processors, devices, discrete logic, circuits, etc.
Turning in detail to FIGS. 6A and 6B, the depicted example process may be used to determine locations of audience members (e.g., the person 110 of FIG. 1). Initially, a neural network for use in determining audience member locations is trained (block 602). For example, the apparatus 200 of FIG. 2 may train its neural network 222 for use in determining audience member locations based on time-domain analyses of the reference and non-stationary audio samples 202 and 204 and/or the apparatus 300 of FIG. 3 may train its neural network 322 for use in determining audience member locations based on frequency-domain analyses of the reference and non-stationary audio samples 202 and 204. In either instance, during a training process, the apparatus 200 and/or the apparatus 300 is/are provided with multiple sets of reference and non-stationary audio samples from the stationary reference unit 106 and the portable meter 108 corresponding to the training distances (dT1, dT2, dT3) along with the distance values of the training distances (dT1, dT2, dT3) corresponding to each audio sample set. In this manner, the apparatus 200 and/or the apparatus 300 can train its respective neural network 222 and/or 322 by inputting resulting correlation coefficients (r) and the distance values of the training distances (dT1, dT2, dT3) using any suitable neural network training technique so that the neural network(s) 222 and/or 322 can configure its/their logic to determine distance measures (e.g., the monitored distance dM of FIG. 1) based on subsequent reference and non-stationary audio samples 202 and 204 collected by the stationary reference unit 106 and the portable meter 108.
After the neural network training phase of block 602, the subsequent operations of the flow diagram of FIGS. 6A and 6B correspond to a monitoring phase in which the apparatus 200 and/or the apparatus 300 received reference and non-stationary audio samples (e.g., the reference and non-stationary audio samples 202 and 204 of FIGS. 2 and 3) to determine monitored distances (dM) of the person 110 of FIG. 1.
Initially, during a monitoring phase, the processing unit 104 determines whether it should perform a time-domain analysis (block 604) using the apparatus 200 of FIG. 2 to determine the monitored distance (dM) of the person 110. In some examples, the processing unit 104 may be pre-programmed or configured by a user (e.g., an installation technician) to always perform time-domain analyses or to always perform frequency-domain analyses. Such configuration may be implemented by setting a software register option (e.g., a flag or bit) in the processing unit 104 or by setting a hardware switch on the processing unit 104 specifying which type of analysis to perform. In instances in which the processing unit 104 is only provided with the example apparatus 200 of FIG. 2 to perform time-domain analyses or is only provided with the example apparatus 300 of FIG. 3 to perform frequency-domain analyses, the decision of block 604 may be omitted. That is, in instances in which the processing unit 104 is only provided with the apparatus 200 to perform time-domain analyses, the processing unit 104 performs the operations of blocks 604, 608, 610, 612, 614, 616, 618, 620, and 622 of FIG. 6A (without needing to perform the decision of block 604) to determine monitored distances (dM) of the person 110 after a neural network training phase in which the neural network 222 (FIG. 2) is trained. Whereas, in instances in which the processing unit 104 is only provided with the apparatus 300 to perform frequency-domain analyses, the processing unit 104 performs the operations shown in FIG. 6B (without needing to perform the decision of block 604) to determine monitored distances (dM) of the person 110 after a neural network training phase in which the neural network 322 (FIG. 3) is trained.
If the processing unit 104 determines at block 604 that it should use a time-domain analysis, the processing unit 104 of FIG. 1 receives the reference and non-stationary audio samples 202 and 204 (block 606). In the illustrated example, the quantity of audio samples collected by each of the stationary reference unit 106 and the portable meter 108 is based on the desired size of the correlation window and time-shift range for determining the correlation coefficients (r). In the illustrated example, for time-domain analyses, each of the stationary reference unit 106 and the portable meter 108 may collect 750 ms of audio samples plus an additional quantity of audio samples equal to twice the duration of the desired time-shift range.
The sub-samplers 208 and 210 of FIG. 2 perform sub-sampling operations on respective ones of the reference and non-stationary audio samples 202 and 204 (block 608) to generate the reference and audio sample subsets 402 and 404 of FIGS. 4A-4C. The down-scaler 212 of FIG. 2 down-scales the non-stationary audio sample subset 404 (block 610). The correlator 214 (FIG. 2) then determines correlation coefficients (r) (e.g., correlation analysis results) (block 612) based on the reference and audio sample subsets 402 and 404. For instances in which neither of the reference and non-stationary audio samples 202 and 204 have a DC bias, the correlator 214 may perform the operation of block 612 using Equation 1 above. An example process that may be used to implement Equation 1 is described below in connection with the flow diagram of FIG. 7. Alternatively, if the reference audio samples 202 do not have a DC bias, but the non-stationary audio samples 204 do have a DC bias, the correlator 214 may perform the operation of block 612 using Equation 2 above. An example process that may be used to implement Equation 2 is described below in connection with the flow diagram of FIGS. 8A and 8B.
The processing unit 104 then selects three maximum correlation coefficients (rmax) and three minimum correlation coefficients (rmin) and corresponding time shift positions (tshift) (block 614). For example, the processing unit 104 uses the correlation analyzer 216 of FIG. 2 to perform the operation of block 614 for the time-domain analysis. The processing unit 104 then determines the maximum-amplitude values of the reference and non-stationary audio samples 202 and 204 (block 616). For example, the processing unit 104 uses the maximum level detectors 218 and 220 of FIG. 2 to perform the operation of block 616 for the time-domain analysis.
The processing unit 104 then determines the monitored distance (dM) (block 618) based on the three maximum correlation coefficients (rmax), the three minimum correlation coefficients (rmin), the corresponding time shift positions (tshift), and the maximum-amplitude values of the reference and non-stationary audio samples 202 and 204. For example, the processing unit 104 uses the neural network 222 of FIG. 2 to perform the operation of block 618 for the time-domain analysis.
The processing unit 104 determines the location of the person 110 in the monitored environment 100 of FIG. 1 based on the monitored distance (dM) (block 620). For example, the processing unit 104 may determine that the location of the person 110 is along the circumference 116 about the media presentation device 112 having a radius equal to the monitored distance (dM).
In the illustrated example, the location of the person 110 is subsequently used to selectively credit media content (block 622) presented via the media presentation device 112. For example, the processing unit 104 may elect to award exposure credit to media content as having been exposed to the person 110 if the person 110 was within a threshold distance (e.g., 10 feet) of the media presentation device 112. Alternatively, the processing unit 104 may elect to withhold crediting the media content with an exposure credit if the determined location indicates that the person 110 was not within a threshold distance (e.g., 10 feet) of the media presentation device 112 to have been adequately exposed to the media content. The threshold distance is selected to be equal to a distance that is sufficiently proximate the media presentation device 112 to assure the person was exposed to the program. The distance selected may be dependent on the environment of use. The example process of FIG. 6A then ends.
Turning now to FIG. 6B, for instances in which the processing unit 104 determines at block 604 (FIG. 6A) that it should perform a frequency-domain analysis (i.e., it should not perform a time-domain analysis), the processing unit 104 of FIG. 1 receives the reference and non-stationary audio samples 202 and 204 (block 624). In the illustrated example, the quantity of audio samples collected by each of the stationary reference unit 106 and the portable meter 108 is based on the desired size of the correlation window and time-shift range for determining the correlation coefficients (r). In the illustrated example, to perform the frequency-domain analysis, each of the stationary reference unit 106 and the portable meter 108 may collect 32 ms of audio samples plus an additional quantity of audio samples equal to twice the duration of the desired time-shift range. Frequency-domain analyses may be advantageously used to improve battery charge performance of the portable meter 108, because less audio samples are required under frequency-domain analyses to achieve similar results as those produced using time-domain analyses.
The processing unit 104 uses the apparatus 300 of FIG. 3 to determine the correlation coefficients (e.g., correlation analysis results) at block 626. In such instances, the correlator 314 uses Equation 4 above to determine the correlation coefficients (r) at block 626 when neither of the reference and non-stationary audio samples 202 and 204 have a DC bias. An example process that may be used to implement Equation 4 is described below in connection with the flow diagram of FIG. 9. Alternatively, if the reference audio samples 202 do not have a DC bias, but the non-stationary audio samples 204 do have a DC bias, the correlator 314 performs the operation of block 626 using Equation 5 above. An example process that may be used to implement Equation 5 is described below in connection with the flow diagram of FIGS. 10A-10C. Alternatively if both of the reference and non-stationary audio samples 202 and 204 have a DC bias or if it is not known whether either of the reference and non-stationary audio samples 202 and 204 have DC biases, the correlator 314 performs the operation of block 626 using Equation 6 as described below in connection with the flow diagram of FIGS. 11A-11C.
The processing unit 104 then selects six maximum correlation coefficients (rmax) and corresponding time shift positions (tshift) (block 628). In the illustrated example, the processing unit 104 uses the correlation analyzer 316 of FIG. 3 to perform the operation of block 628 for the frequency-domain analysis. The processing unit 104 then determines the maximum-amplitude values of the reference and non-stationary audio samples 202 and 204 (block 630). In the illustrated example, the processing unit 104 uses the maximum level detectors 318 and 320 of FIG. 3 to perform the operation of block 630 for the frequency-domain analysis.
The processing unit 104 then determines the monitored distance (dM) (block 632) based on the six maximum correlation coefficients (rmax), the corresponding time shift positions (tshift), and the maximum-amplitude values of the reference and non-stationary audio samples 202 and 204. In the illustrated example, the processing unit 104 uses the neural network 322 of FIG. 3 to perform the operation of block 618 for the frequency-domain analysis.
The processing unit 104 determines the location of the person 110 in the monitored environment 100 of FIG. 1 based on the monitored distance (dM) (block 634). For example, the processing unit 104 may determine that the location of the person 110 is along the circumference 116 about the media presentation device 112 having a radius equal to the monitored distance (dM).
In the illustrated example, the location of the person 110 is subsequently used to selectively credit media content (block 636) presented via the media presentation device 112. For example, the processing unit 104 may elect to award exposure credit to media content as having been exposed to the person 110 if the person 110 was within a threshold distance (e.g., 10 feet) of the media presentation device 112. Alternatively, the processing unit 104 may elect to withhold crediting the media content with an exposure credit if the determined location indicates that the person 110 was not within a threshold distance (e.g., 10 feet) of the media presentation device 112 to have been adequately exposed to the media content. The threshold distance is selected to be equal to a distance that is sufficiently proximate the media presentation device 112 to assure the person was exposed to the program. The distance selected may be dependent on the environment of use. The example process of FIG. 6B then ends.
Turning now to FIG. 7, the depicted example flow diagram is representative of an example process that may be used by the correlator 214 of FIG. 2 to determine correlation coefficients (r) based on Equation 1 above. In the illustrated example, the example process of FIG. 7 may be used to implement the operation of block 612 of FIG. 6A.
Initially, the correlator 214 determines a sum of squared reference audio samples (Σx2) (block 702) based on the square of each audio sample in the reference audio samples 202. The correlator 214 also determines a sum of squared non-stationary audio samples (Σy2) (block 704) based on the square of each audio sample in the non-stationary audio samples 204. The correlator 214 multiplies the sums of squares (Σx2) and (Σy2) to determine a product ((Σx2)(Σy2)) (block 706) and determines the square root of the product (√{square root over ((Σx2)(Σy2))}{square root over ((Σx2)(Σy2))}) (block 708). The correlator 214 determines a sum of cross-sample products (Σxy) (block 710) based on products (xy) of each reference audio sample 202 (x) and a corresponding one of the non-stationary audio samples 204 (y). The correlator 214 then determines the correlation coefficient (r) associated with the reference and non-stationary audio samples 202 and 204 by dividing the sum of cross-sample products (Σxy) by the square root of the first product (√{square root over ((Σx2)(Σy2))}{square root over ((Σx2)(Σy2))}) (block 712) as shown in Equation 1 (block 710). The example process of FIG. 7 then ends.
Turning now to FIGS. 8A and 8B, the depicted example flow diagram is representative of an example process that may be performed by the correlator 214 of FIG. 2 to determine correlation coefficients (r) based on Equation 2 above. In the illustrated example, the example process of FIGS. 8A and 8B may be used to implement the operation of block 612 of FIG. 6A.
Initially, the correlator 214 determines a sum of squared reference audio samples (Σx2) (block 802) based on the square of each audio sample in the reference audio samples 202. The correlator 214 then multiplies the sum of squared reference audio samples (Σx2) by a total sample quantity value (n) (block 804) representing the total quantity of audio samples in the reference and non-stationary audio samples 202 and 204 being processed to determine a first product (n(Σx2)) and determines a first square root of the first product (√{square root over (n(Σx2))}) (block 806). The correlator 214 also determines a sum of squared non-stationary audio samples (Σy2) (block 808) based on the square of each audio sample in the non-stationary audio samples 204 and multiplies the sum of squared non-stationary audio samples (Σy2) by the total sample quantity value (n) to determine a second product (n(Σy2)) (block 810). The correlator 214 determines a square of a sum of the non-stationary audio samples 204 ((Σy)2) (block 812) and determines a second square root of the difference between the second product (n(Σy2)) and the square of the sum of the non-stationary audio samples 204 ((Σy)2) (i.e., √{square root over (n(Σy2)−(Σy)2)}{square root over (n(Σy2)−(Σy)2)}) (block 814). The correlator 214 multiplies the first square root result (√{square root over (n(Σx2))}) by the second square root result (√{square root over (n(Σy2)−(Σy)2)}{square root over (n(Σy2)−(Σy)2)}) to determine a third product (√{square root over (n(Σx2))}√{square root over (n(Σy2)−(Σy)2)}{square root over (n(Σy2)−(Σy)2)}) (block 816) (FIG. 8B). The correlator 214 determines a sum of cross-sample products (Σxy) (block 818) based on products (xy) of each reference audio sample 202 (x) and a corresponding one of the non-stationary audio samples 204 (y). The correlator 214 multiplies the sum of cross-sample products (Σxy) by the total sample quantity value (n) to determine a fifth product (nΣxy) (block 820). The correlator 214 then determines the correlation coefficient (r) by dividing the fifth product (nΣxy) by the third product (√{square root over (n(Σx2))}√{square root over (n(Σy2)−(Σy)2)}{square root over (n(Σy2)−(Σy)2)}) as shown in Equation 2 above (block 822). The example process of FIGS. 8A and 8B then ends.
Turning now to FIG. 9, the depicted example flow diagram is representative of an example process that may be performed by the correlator 314 of FIG. 3 to determine correlation coefficients (r) based on Equation 4 above. In the illustrated example, the example process of FIG. 9 may be used to implement the operation of block 626 of FIG. 6B.
Initially, the correlator 314 determines a sum of reference complex conjugate pairs ((ΣXX*)) of the reference frequency-domain representation 502 (block 902). The correlator 314 determines a sum of non-stationary complex conjugate pairs ((ΣYY*)) of the non-stationary frequency-domain representation 504 (block 904). The correlator 314 then multiplies the sum of reference complex conjugate pairs ((ΣXX*)) by the sum of non-stationary complex conjugate pairs ((ΣYY*)) to determine a first product ((ΣXX*)(ΣYY*)) (block 906). The correlator 314 determines the square root of the first product (√{square root over ((ΣXX*)(ΣYY*))}{square root over ((ΣXX*)(ΣYY*))}) (block 908). The correlator 314 multiplies reference complex numbers of the reference frequency-domain representation 502 by corresponding non-stationary conjugates of the non-stationary frequency-domain representation 504 to determine second products (XY*) (block 910) and determines a third sum of the second products (ΣXY*) (block 912). The correlator 314 determines a correlation coefficient (r) associated with the reference and non-stationary audio samples 202 and 204 by dividing the third sum of the second products (ΣXY*) by the square root of the first product (√{square root over ((ΣXX*)(ΣYY*))}{square root over ((ΣXX*)(ΣYY*))}) as shown in Equation 4 above (block 914). The example process of FIG. 9 then ends.
Turning now to FIGS. 10A-10C, the depicted example flow diagram is representative of an example process that may be performed by the correlator 314 of FIG. 3 to determine correlation coefficients (r) based on Equation 5 above. In the illustrated example, the example process of FIGS. 10A-10C may be used to implement the operation of block 626 of FIG. 6B.
Initially, the correlator 314 determines a sum of reference complex conjugate pairs ((ΣXX*)) of the reference frequency-domain representation 502 (block 1002). The correlator 314 multiplies the sum of reference complex conjugate pairs ((ΣXX*)) by a total sample quantity value (n) representing the total quantity of the audio samples of the reference and non-stationary audio samples 202 and 204 to determine a first product (n(ΣXX*)) (block 1004). The correlator 314 determines a square root of the first product (√{square root over (n(ΣXX*))}) (block 1006). The correlator 314 determines a sum of non-stationary complex conjugate pairs ((ΣYY*)) of the non-stationary frequency-domain representation 504 (block 1008) and multiplies the sum of non-stationary complex conjugate pairs ((ΣYY*)) by the total sample quantity value (n) to determine a second product (n(ΣYY*)) (block 1010). The correlator 314 determines a sum of non-stationary complex numbers ((ΣY)) of the non-stationary frequency-domain representation 504 (block 1012) and determines a conjugate of the sum of non-stationary complex numbers ((ΣY)) (i.e., (ΣY)*) (block 1014) (FIG. 10B). The correlator 314 multiplies the sum of non-stationary complex numbers ((ΣY)) by its conjugate ((ΣY)*) to generate a third product ((ΣY)(ΣY)*) (block 1016). The correlator 314 determines a second square root of the difference between the second product (n(ΣYY*)) and the third product ((ΣY)(ΣY)*) (i.e., √{square root over (n(ΣYY*)−(ΣY)(ΣY)*)}{square root over (n(ΣYY*)−(ΣY)(ΣY)*)}{square root over (n(ΣYY*)−(ΣY)(ΣY)*)}) (block 1018). The correlator 314 multiplies the first square root (√{square root over (n(ΣXX*))}) by the second square root (√{square root over (n(ΣYY*)−(ΣY)(ΣY)*)}{square root over (n(ΣYY*)−(ΣY)(ΣY)*)}{square root over (n(ΣYY*)−(ΣY)(ΣY)*)}) to determine a third product (√{square root over (n(ΣXX*))}√{square root over (n(ΣYY*)−(ΣY)(ΣY)*)}{square root over (n(ΣYY*)−(ΣY)(ΣY)*)}{square root over (n(ΣYY*)−(ΣY)(ΣY)*)}) (block 1020). The correlator 314 multiplies reference complex numbers (X) of the reference frequency-domain representation 502 by corresponding non-stationary conjugates (Y*) of the non-stationary frequency-domain representation 504 to determine fourth products (XY*) (block 1022) and determines a sum of the fourth products (ΣXY*) (block 1024). The correlator 314 multiplies the sum of the fourth products (ΣXY*) by the total sample quantity value (n) to determine a fifth product (nΣXY*) (block 1026) (FIG. 10C). The correlator 314 determines a correlation coefficient (r) by dividing the fifth product (nΣXY*) by the third product (√{square root over (n(ΣXX*))}√{square root over (n(ΣYY*)−(ΣY)(ΣY)*)}{square root over (n(ΣYY*)−(ΣY)(ΣY)*)}{square root over (n(ΣYY*)−(ΣY)(ΣY)*)}) as shown in Equation 5 (block 1028). The example process of FIGS. 10A-10C then ends.
Turning now to FIGS. 11A-11C, the depicted example flow diagram is representative of an example process that may be performed by the correlator 314 of FIG. 3 to determine correlation coefficients (r) based on Equation 6 above. In the illustrated example, the example process of FIGS. 11A-11C may be used to implement the operation of block 626 of FIG. 6B.
Initially, the correlator 314 determines a sum of reference complex conjugate pairs ((ΣXX*)) of the reference frequency-domain representation 502 (block 1102). The correlator 314 multiplies the sum of reference complex conjugate pairs ((ΣXX*)) by a total sample quantity value (n) representing the total quantity of the audio samples of the reference and non-stationary audio samples 202 and 204 to determine a first product (n(ΣXX*)) (block 1104). The correlator 314 determines a sum of reference complex numbers ((ΣX)) of the reference frequency-domain representation 504 (block 1106) and determines a conjugate of the sum of reference complex numbers ((ΣX)) (i.e., (ΣX)*) (block 1108). The correlator 314 multiplies the sum of reference complex numbers ((ΣX)) by its conjugate ((ΣX)*) to generate a second product ((ΣX)(ΣX)*) (block 1110). The correlator 314 determines a first square root of the difference between the first product (n(ΣXX*)) and the second product ((ΣX)(ΣX)*) (i.e., √{square root over (n(ΣXX*)−(ΣX)(ΣX)*)}{square root over (n(ΣXX*)−(ΣX)(ΣX)*)}{square root over (n(ΣXX*)−(ΣX)(ΣX)*)}) (block 1112).
The correlator 314 determines a sum of non-stationary complex conjugate pairs ((ΣYY*)) of the non-stationary frequency-domain representation 504 (block 1114) (FIG. 11B) and multiplies the sum of non-stationary complex conjugate pairs ((ΣYY*)) by the total sample quantity value (n) to determine a third product (n(ΣYY*)) (block 1116). The correlator 314 determines a sum of non-stationary complex numbers ((ΣY)) of the non-stationary frequency-domain representation 504 (block 1118) and determines a conjugate of the sum of non-stationary complex numbers ((ΣY)) (i.e., (ΣY)*) (block 1120). The correlator 314 multiplies the sum of non-stationary complex numbers ((ΣY)) by its conjugate ((ΣY)*) to generate a fourth product ((ΣY)(ΣY)*) (block 1122). The correlator 314 determines a second square root of the difference between the third product (n(ΣYY*)) and the fourth product ((ΣY)(ΣY)*) (i.e., √{square root over (n(ΣYY*)−(ΣY)(ΣY)*)}{square root over (n(ΣYY*)−(ΣY)(ΣY)*)}{square root over (n(ΣYY*)−(ΣY)(ΣY)*)}) (block 1124).
The correlator 314 multiples the first and second square roots to determine a fifth product (√{square root over (n(ΣXX*)−(ΣX)(ΣX)*)}{square root over (n(ΣXX*)−(ΣX)(ΣX)*)}{square root over (n(ΣXX*)−(ΣX)(ΣX)*)}√{square root over (n(ΣYY*)−(ΣY)(ΣY)*)}{square root over (n(ΣYY*)−(ΣY)(ΣY)*)}{square root over (n(ΣYY*)−(ΣY)(ΣY)*)}) (block 1126) (FIG. 11C). The correlator 314 multiplies reference complex numbers (X) of the reference frequency-domain representation 502 by corresponding non-stationary conjugates (Y*) of the non-stationary frequency-domain representation 504 to determine sixth products (XY*) (block 1128) and determines a sum of the sixth products (Σ(XY*)) (block 1130). The correlator 314 multiplies the sum of the sixth products (Σ(XY*)) by the total sample quantity value (n) to determine a seventh product (nΣ(XY*)) (block 1132).
The correlator 314 determines a sum of non-stationary conjugates ((ΣY*)) (block 1134). The correlator 314 multiples the sum of non-stationary conjugates by the sum of reference complex numbers to determine an eighth product ((ΣX)(ΣY*)) (block 1136). The correlator 314 determines a correlation coefficient (r) by dividing the difference between the eighth product ((ΣX)(ΣY*)) and seventh product (nΣ(XY*)) by the fifth product (√{square root over (n(ΣXX*)−(ΣX)(ΣX)*)}{square root over (n(ΣXX*)−(ΣX)(ΣX)*)}{square root over (n(ΣXX*)−(ΣX)(ΣX)*)}√{square root over (n(ΣYY*)−(ΣY)(ΣY)*)}{square root over (n(ΣYY*)−(ΣY)(ΣY)*)}{square root over (n(ΣYY*)−(ΣY)(ΣY)*)}) as shown in Equation 6 (block 1138). The example process of FIGS. 11A-11C then ends.
FIG. 12 is a block diagram of an example processor system 1210 that may be used to implement the example apparatus, methods, and systems described herein. As shown in FIG. 12, the processor system 1210 includes a processor 1212 that is coupled to an interconnection bus 1214. The processor 1212 may be any suitable processor, processing unit, or microprocessor. Although not shown in FIG. 12, the system 1210 may be a multi-processor system and, thus, may include one or more additional processors that are identical or similar to the processor 1212 and that are communicatively coupled to the interconnection bus 1214.
The processor 1212 of FIG. 12 is coupled to a chipset 1218, which includes a memory controller 1220 and an input/output (I/O) controller 1222. A chipset provides I/O and memory management functions as well as a plurality of general purpose and/or special purpose registers, timers, etc. that are accessible or used by one or more processors coupled to the chipset 1218. The memory controller 1220 performs functions that enable the processor 1212 (or processors if there are multiple processors) to access a system memory 1224 and a mass storage memory 1225.
In general, the system memory 1224 may include any desired type of volatile and/or non-volatile memory such as, for example, static random access memory (SRAM), dynamic random access memory (DRAM), flash memory, read-only memory (ROM), etc. The mass storage memory 1225 may include any desired type of mass storage device including hard disk drives, optical drives, tape storage devices, etc.
The I/O controller 1222 performs functions that enable the processor 1212 to communicate with peripheral input/output (I/O) devices 1226 and 1228 and a network interface 1230 via an I/O bus 1232. The I/O devices 1226 and 1228 may be any desired type of I/O device such as, for example, a keyboard, a video display or monitor, a mouse, etc. The network interface 1230 may be, for example, an Ethernet device, an asynchronous transfer mode (ATM) device, an 802.11 device, a digital subscriber line (DSL) modem, a cable modem, a cellular modem, etc. that enables the processor system 1210 to communicate with another processor system.
While the memory controller 1220 and the I/O controller 1222 are depicted in FIG. 12 as separate functional blocks within the chipset 1218, the functions performed by these blocks may be integrated within a single semiconductor circuit or may be implemented using two or more separate integrated circuits.
Although the above discloses example methods, apparatus, systems, and articles of manufacture including, among other components, firmware and/or software executed on hardware, it should be noted that such methods, apparatus, systems, and articles of manufacture are merely illustrative and should not be considered as limiting. Accordingly, while the above describes example methods, apparatus, systems, and articles of manufacture, the examples provided are not the only ways to implement such methods, apparatus, systems, and articles of manufacture.
Although certain example methods, apparatus, systems, and articles of manufacture have been described herein, the scope of coverage of this patent is not limited thereto. On the contrary, this patent covers all methods, apparatus and articles of manufacture fairly falling within the scope of the claims of this patent.

## Claims

1. A method to determine a location of an audience member, the method comprising:
generating a correlation analysis result based on correlating first audio samples from a stationary audio detector with second audio samples from a portable audio detector carried by the audience member;
determining via a neural network the location of the audience member based on the correlation analysis result; and
crediting media presented via a media delivery device with an exposure credit when the audience member was located within a threshold distance of the media delivery device as indicated by the location of the audience member.

2. A method as defined in claim 1, wherein generating the correlation analysis result comprises:
determining a first sum of squares of each audio sample in the first audio samples;
multiplying the first sum of squares by a quantity value representing a total audio sample quantity of the first and second audio samples to determine a first product;
determining a square root of the first product;
determining a second sum of squares of each audio sample in the second audio samples;
multiplying the second sum of squares by the quantity value to determine a second product;
determining a square of a third sum of the second audio samples;
determining a second square root of a difference between the second product and the square of the third sum;
multiplying the first square root by the second square root to determine a third product;
determining a fourth sum of fourth products between corresponding audio samples of the first and second audio samples;
multiplying the fourth sum by the quantity value to determine a fifth product; and
determining a correlation coefficient by dividing the fifth product by the third product.

3. A method as defined in claim 2, wherein the correlation analysis result is determined when an average value of the first audio samples is zero and an average value of the second audio samples is not zero.

4. A method as defined in claim 1, wherein the first audio samples are a first sub-sample set of first total audio samples collected by the stationary audio detector and the second audio samples are a second sub-sample set of second total audio samples collected by the portable audio detector.

5. A method as defined in claim 1, wherein determining the location of the audience member is based on a first maximum-amplitude audio sample from the first audio samples and a second maximum-amplitude audio sample from the second audio samples.

6. A method as defined in claim 1, further comprising:
determining a first frequency-domain representation of the first audio samples and a second frequency-domain representation of the second audio samples before generating the correlation analysis result; and
generating the correlation analysis result by correlating the first frequency-domain representation of the first audio samples with the second frequency-domain representation of the second audio samples.

7. A method as defined in claim 6, wherein generating the correlation analysis result comprises:
determining a first sum of complex conjugate pairs of the first frequency-domain representation;
determining a second sum of complex conjugate pairs of the second frequency-domain representation;
multiplying the first sum of complex conjugate pairs by the second sum of complex conjugate pairs to determine a first product;
determining a square root of the first product;
multiplying complex numbers of the first frequency-domain representation by corresponding conjugates of the second frequency-domain representation to determine second products;
determining a third sum of the second products; and
determining a correlation coefficient associated with the first and second audio samples by dividing the third sum of the second products by the square root of the first product.

8. A method as defined in claim 7, wherein the correlation analysis result is determined when an average value of the first audio samples is zero and an average value of the second audio samples is zero.

9. A method as defined in claim 6, wherein generating the correlation analysis result comprises:
determining a first sum of complex conjugate pairs of the first frequency-domain representation;
multiplying the first sum of complex conjugate pairs by a quantity value representing a total quantity of audio samples of the first and second audio samples to determine a first product;
determining a first square root of the first product;
determining a second sum of complex conjugate pairs of the second frequency-domain representation;
multiplying the second sum of complex conjugate pairs by the quantity value to determine a second product;
determining a third sum of complex numbers of the second frequency-domain representation;
determining a conjugate of the third sum of complex numbers;
multiplying the third sum of complex numbers by the conjugate to generate a third product;
determining a second square root of a difference between the second product and the third product;
multiplying the first square root by the second square root to determine a third product;
multiplying complex numbers of the first frequency-domain representation by corresponding conjugates of the second frequency-domain representation to determine fourth products;
determining a fourth sum of the fourth products;
multiplying the fourth sum by the quantity value to determine a fifth product; and
determining a correlation coefficient by dividing the fifth product by the third product.

10. A method as defined in claim 9, wherein the correlation analysis result is determined when an average value of the first audio samples is zero and an average value of the second audio samples is not zero.

11. A method as defined in claim 6, wherein generating the correlation analysis result comprises:
determining a first sum of reference complex conjugate pairs;
multiplying the first sum by a quantity value representing a total quantity of audio samples of the first and second audio samples to determine a first product;
determining a second sum of reference complex numbers;
determining a conjugate of a third sum of complex numbers of the first frequency-domain representation;
multiplying the second sum by the conjugate of the third sum to determine a second product;
determining a first square root of a difference between the first and second products;
determining a fourth sum of complex conjugate pairs of the second frequency-domain representation;
multiplying the fourth sum by the quantity value to determine a third product;
determining a fifth sum of complex numbers of the second frequency-domain representation;
determining a conjugate of a sixth sum of complex numbers of the second frequency-domain representation;
multiplying the fifth and sixth sums to determine a fourth product;
determining a second square root of a difference between the third and fourth products;
multiplying the first and second square roots to determine a fifth product;
multiplying complex numbers of the first frequency-domain representation by corresponding conjugates of the second frequency-domain representation to determine sixth products;
determining a seventh sum of the sixth products;
multiplying the seventh sum by the quantity value;
determining an eighth sum of conjugates of the second frequency-domain representation;
multiplying the second sum by the eighth sum to determine a seventh product; and
dividing a difference between the eighth product and the seventh product by the fifth product.

12. A method to determine a location of an audience member, the method comprising:
generating a correlation analysis result based on correlating first audio samples from a stationary audio detector with second audio samples from a portable audio detector carried by the audience member; and
determining via a neural network the location of the audience member based on the correlation analysis result,
wherein generating the correlation analysis result comprises:
determining a first sum of squares of each audio sample in the first audio samples;
determining a second sum of squares of each audio sample in the second audio samples;
multiplying the first sum of squares by the second sum of squares to determine a first product;
determining a square root of the first product;
determining a third sum of second products between corresponding audio samples of the first and second audio samples; and
determining a correlation coefficient associated with the first and second audio samples by dividing the third sum of the second products by the square root of the first product.

13. A method as defined in claim 12, wherein generating the correlation analysis result is performed when an average value of the first audio samples is zero and an average value of the second audio samples is zero.

14. An apparatus to determine a location of an audience member, the apparatus comprising:
a correlator structured to generate a correlation analysis result based on correlating first audio samples from a stationary audio detector with second audio samples from a portable audio detector carried by the audience member;
a neural network structured to determine the location of the audience member based on the correlation analysis result; and
a processor to credit media presented via a media delivery device with an exposure credit when the audience member was located within a threshold distance of the media delivery device as indicated by the location of the audience member.

15. An apparatus as defined in claim 14, wherein the correlator is to generate the correlation analysis result by:
determining a first sum of squares of each audio sample in the first audio samples;
multiplying the first sum of squares by a quantity value representing a total audio sample quantity of the first and second audio samples to determine a first product;
determining a square root of the first product;
determining a second sum of squares of each audio sample in the second audio samples;
multiplying the second sum of squares by the quantity value to determine a second product;
determining a square of a third sum of the second audio samples;
determining a second square root of a difference between the second product and the square of the third sum;
multiplying the first square root by the second square root to determine a third product;
determining a fourth sum of fourth products between corresponding audio samples of the first and second audio samples;
multiplying the fourth sum by the quantity value to determine a fifth product; and
determining a correlation coefficient by dividing the fifth product by the third product.

16. An apparatus to determine a location of an audience member, the apparatus comprising:
a correlator structured to generate a correlation analysis result based on correlating first audio samples from a stationary audio detector with second audio samples from a portable audio detector carried by the audience member; and
a neural network structured to determine the location of the audience member based on the correlation analysis result,
wherein the correlator is to generate the correlation analysis result by:
determining a first sum of squares of each audio sample in the first audio samples;
determining a second sum of squares of each audio sample in the second audio samples;
multiplying the first sum of squares by the second sum of squares to determine a first product;
determining a square root of the first product;
determining a third sum of second products between corresponding audio samples of the first and second audio samples; and
determining a correlation coefficient associated with the first and second audio samples by dividing the third sum of the second products by the square root of the first product.

17. A tangible machine accessible storage medium comprising instructions that, when executed, cause a machine to at least:
generate a correlation analysis result based on correlating first audio samples from a stationary audio detector with second audio samples from a portable audio detector carried by the audience member;
determine via a neural network the location of the audience member based on the correlation analysis result; and
credit media presented via a media delivery device with an exposure credit when the audience member was located within a threshold distance of the media delivery device as indicated by the location of the audience member.

18. The tangible machine accessible storage medium of claim 17, wherein the first audio samples are a first sub-sample set of first total audio samples collected by the stationary audio detector and the second audio samples are a second sub-sample set of second total audio samples collected by the portable audio detector.

19. The tangible machine accessible storage medium of claim 17, wherein the instructions further cause the machine to determine the location of the audience member based on a first maximum-amplitude audio sample from the first audio samples and a second maximum-amplitude audio sample from the second audio samples.

20. The tangible machine accessible storage medium of claim 17, wherein the instructions further cause the machine to:
determine a first frequency-domain representation of the first audio samples and a second frequency-domain representation of the second audio samples before generating the correlation analysis result, the correlation analysis result being generated based on correlating the first frequency-domain representation of the first audio samples with the second frequency-domain representation of the second audio samples.