Skip to main content
mjmartin13
Participating Frequently
February 22, 2017
Answered

How to equalize a frequency response curve from a hydrophone with a flat frequency response (+/- 3 dB) up to 150 kHz but ranges up to 288 kHz

  • February 22, 2017
  • 1 reply
  • 3082 views

Hi,

I am working with a high-frequency recording underwater hydrophone that is not supplied with a frequency response curve. The sampling rate is 576 kHz, therefore the frequency range reaches up to 288 kHz. The hydrophone has a flat frequency response from 20 Hz - 150 kHz (+/- 3 dB) and then suffers a roll-off of 6 dB per octave above 150 kHz. The next octave would include all frequencies up to 300 kHz. I would like to know if I am applying the appropriate equalization settings in order to achieve a full normalized response over the range of the hydrophone. Ideally, I need to have a full flat response in order to calculate measurements of ultrasonic dolphin echolocation clicks in which some of the energy from each click extends well above 150 kHz.

First, I calculated the roll-off correction values from 151 kHz - 300 kHz based on the 6 dB roll-off. The correction value at 288 kHz is 5.64 dB. However, I am not sure if I need to include an additional 3 dB from the flat response portion of the range up to 150 kHz.

These are the steps I have used so far:

I selected the first filtering option under 'Filter and EQ'  called 'FFT filter effect' and created a curve from 150 kHz to 288 kHz increasing from 0 dB to 5.64 dB (from the roll-off correction value at 288 kHz). I selected the scale to be linear and selected  'spline curve' so the curve was smoother. Under the advanced settings tab, I set FFT window size to 512 and the type of window (not sure which one is best, but I scanned my .wav files using Blackman-Harris so I used this). Then clicked apply. Potentially, if the +/- 3 dB needs to be added, I could add a curve increasing from 0 to 3 dB from 20 Hz to 150 kHz and then start the second curve at 3 dB at 150 kHz and have it extend to 8.64 dB as it reaches 288 kHz.

Screenshot attached of what I did in Adobe...hope I am making sense!  In the screenshot inside the Frequency Analysis box, the blue line is from a click I scanned before apply the FFT filter and the orange line is from the same click I scanned after applying the FFT filter.

Please let me know if this is the correct method to create a full response of my hydrophone range and if it is necessary to the additional 3 dB from the flat response up to 150 kHz.

Thank you kindly!

Morgan

    This topic has been closed for replies.
    Correct answer SteveG_AudioMasters_

    Hey Steve,

    Thanks again, I have been doing some reading and realized I mistook what you meant by performing a transform. I'm so used to seeing FFT as an acronym in the literature and code that I forgot it was the main transform. Doh! I will eventually remember definitions of all this jargon...I understand what you meant by not low-pass filtering out any of the HF data from the clicks- I should not filter until it has been transformed from the time domain to the frequency domain. I was getting confused by phrasing in the literature when discussing how people corrected for their frequency response curves of hydrophones. I kept reading that authors 'applied the inverse of the system transfer function to correct the received levels according to the sensitivity of the recording system.' Please correct me if I am wrong but I think they are referring to applying the inverse of the FFT in order to revert back to the time domain in order to use the corrected waveform of their clicks.

    As I was starting this forum, I asked a scientist via email about how she corrected for her frequency response curve. (She had a similar design set-up as my study.) She was provided the frequency response curve from the hydrophone manufacturer and advised that I should try to create my own curve via calibration of the SoundTrap since it did not come with one. Her reply:

    "You should correct your signal after calculating the FFT of your waveform. Make sure that your spectrum is in dB scale so you just add the correction factors. After that you will have your spectrum levels calibrated. You should ask the SoundTrap team for the frequency response curve below 150 kHz, and if they don't have it, ideally you should measure it by calibrating your hydrophone. Once you have that done, you can calibrate your spectral parameters and if you want to measure RLp-p or RLrms you should apply the inverse FFT to get the corrected waveform. Regarding your code, you will need to modify it so that after calculating the FFT of your input signal, the code adds the frequency response curve"

    So let me get this straight- I can avoid all of this by equalizing my .wav files for the 6 dB roll-off above 150 kHz using the FFT filter setting that I sent a screenshot of previously (0-5.6 dB increase from 150kHz -300 kHz)? If I filter each .wav file in its entirety with the FFT filter setting to correct for the 6dB roll-off, then I can proceed with my analysis of clicks taken from the new modified .wav files to calculate my spectral parameters? I plan to apply a 512 point FFT, rectangular window, with a bandpass filter between 25 kHz and 275 kHz to avoid LF noise and aliasing. Please let me know your thoughts.

    Thanks as always!

    Morgan


    It's perfectly possible to equalise a file of this nature in the time domain first - if you have a correction curve. And unfortunately, SoundTrap don't seem to have produced one of those that's anything like definitive. All they've done is to provide a pistonphone calibration which will give you an accurate amplitude calibration at the repetition rate of the piston, but without a correction curve it's going to be impossible to apply any correction correctly.

    And in the absence of that, the only thing you can do is go with the generalised correction that they've indicated above 150kHz, which I think you've already correctly arrived at. The effect that will have on any transform is affect the levels of the higher frequency components of the signal - this is quite a normal EQ process. What your correspondent has suggested when she says "You should correct your signal after calculating the FFT of your waveform" is an interesting way of going about it - she seems to be suggesting that by defining what she thinks the levels of the upper components should be, you can go back and correct the original file so that this is what the FFT result is. Now, whilst that's technically possible, it's going to be rather a hit-and miss process, to say the very least; normally a correction chart would give you corrected values at different frequencies, and you'd set up a correction filter that provided a smooth curve between them. Alternatively, it would be in the form of a response curve on a graph, and you'd set up a filter that was the inverse of it.

    The thing is though that unless you're absolutely sure of the nature of the stimulus (and I don't see how you can be), you can't apply any correction based upon it. Calibrating the response of a hydrophone for yourself requires a pile of specialised equipment from B&K, and even though I know how to do this, it's not something I'd undertake lightly, and it would be expensive; you have to hire all the kit, and then create the correct environment.

    Bottom line: unless you can get an accurate correction graph or data from SoundTrap, the approach you've suggested in your last paragraph is the only one you can use - and pretty much, if it's a long-distance signal, then the variable sea conditions will be what you are analysing, not the source...

    1 reply

    SteveG_AudioMasters_
    Community Expert
    Community Expert
    February 22, 2017

    No you don't need to add 3dB at all - the response is what it says it is - a -6dB slope from 150kHz. Your real problem is that as it's not a calibrated device that you don't really know where you're starting from; you have a device that, within its normally calibrated range is operating within a 6dB window anyway, so what you think of as 'flat' is likely to be anything but.

    There is more, though. I don't know what the field response of this hydrophone is, and that in itself could cause a response issue, as the compensation for a free field response is different from a diffuse response, and that's different again from a pressure response - and all of these become significant at high frequencies. If you can point me at the spec of the hydrophone, I might be able to glean a little more from it...

    mjmartin13
    Participating Frequently
    February 22, 2017

    Hi!

    Thank you so much for the help and quick reply. I am using a high frequency

    SoundTrap 300 created by Ocean Instruments. My settings were as follows:

    Sampling rate: 576 kHz

    SoundTrap sensitivity: 171.1 dB full scale (I assume this is only up to 150

    kHz)

    Gain setting: High (no further info supplied)

    High pass filter: On = 400 Hz filter

    Flat frequency range: 400 Hz - 150 kHz

    Messages from the SoundTrap contact:

    The anti alias filter in the HF recorders is set to 150 kHz. This means the

    response of the recorder is down 3dB at 150 kHz. The response then drops 6

    dB per octave from there upwards. If we were to set the AA filter to half

    the sample rate you would end up with major distortion because filters are

    not perfect so you must allow some head room for the 6 dB 'roll-off'. The

    signals you are seeing beyond 150 kHz are very high intensity signals that

    are still visible after the above attenuation.

    (This next one corresponds to me asking about how to correct for the 6 dB

    roll-off for the purpose of measuring high frequency dolphin

    clicks-click peak energy is below 150 kHz, FYI):

    This isn't so easy to answer as it depends on exactly what you;re trying to

    achieve. in the simplest case - say you have a narrow band signal at

    or above 150kHz and you want to know the received signal strength then you

    would simply add the reduction in sensitivity to the measured signal

    strength. So, for example, if you measured a 150 kHz signal to be 85 dB,

    you would add another 3 db to get the corrected value of 88 dB. If the

    signal were at 300 Hz (+ 1 octave) then you'd need to add 3685 = 94 dB.

    What I can find on the SoundTrap spec sheet is attached in the screenshot.

    There is a very short spec sheet available online for the HF SoundTrap 300

    but I don't think it has anything not shown here.

    The manual is very short (12 double spaced pgs) and mostly filled with

    general usage and trouble shooting info. I couldn't find anything that was

    helpful for this issue.

    Please let me know if I can help you with any more spec details.

    Thanks so much,

    Morgan

    On Wednesday, 22 February 2017, SteveG(AudioMasters) <

    forums_noreply@adobe.com

    SteveG_AudioMasters_
    Community Expert
    Community Expert
    February 22, 2017

    One thing that might be useful would be the calibration data for your device. I can find that online by the looks of it, but you'd have to provide your individual unit's serial number. But as you say, any data about the actual hydraphone is rather scant. It's slightly strange that they claim that the roll-off is deliberate, and on account of anti-aliasing. Since it's such a shallow rate, I'm going to take a guess that this is based on an over-sampling approach, but it does rather indicate that they don't want you to rely particularly on any signals you get above the turnover point - which is why they are really only specifying the unit as working up to 150kHz.

    Sound transmission in water is somewhat different from what happens in air, and in fact rather more complicated - as my Masters tutor delighted in telling me on several occasions. Even though it's essentially pressure waves that a hydrophone will pick up, it's really difficult to be very specific about what you will have recorded, as the propagation depends on temperature, salinity and acidity and that's before you've considered the impact of multiple path transmission through different layers (especially if it's distant sound), and that's rather important, as it will screw with the phase responses of what you've detected - and with clicks, that will alter the impulse response significantly. The NPL here have a basic absorption calculator which will get you started - it's here: Calculation of absorption of sound in seawater and if you click on the underlying physics link, you'll begin to get some sort of idea as to why you have to be extremely cautious with the results!