/* FunctionEditor_SoundAnalysis.c * * Copyright (C) 1992-2005 Paul Boersma * * This program is free software; you can redistribute it and/or modify * it under the terms of the GNU General Public License as published by * the Free Software Foundation; either version 2 of the License, or (at * your option) any later version. * * This program is distributed in the hope that it will be useful, but * WITHOUT ANY WARRANTY; without even the implied warranty of * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. * See the GNU General Public License for more details. * * You should have received a copy of the GNU General Public License * along with this program; if not, write to the Free Software * Foundation, Inc., 675 Mass Ave, Cambridge, MA 02139, USA. */ /* * pb 2002/07/16 GPL * pb 2002/10/11 added a screen text for unavailable pitch * pb 2002/11/19 default log file names include "~" on Unix and Mach * pb 2002/11/19 added pulses and five separate analysis menus * pb 2003/02/19 clearer wording * pb 2003/02/24 spectral slices * pb 2003/03/03 Extract visible spectrogram: higher degree of oversampling than in editor * pb 2003/03/10 undid previous change because our PostScript code now does image interpolation * pb 2003/03/12 queriable * pb 2003/04/03 Get power at cursor cross * pb 2003/04/11 smaller log settings dialog * pb 2003/05/18 more shimmer measurements * pb 2003/05/20 longestAnalysis replaces the pitch/formant/intensity time steps, * pitch.speckle, and formant.maximumDuration * pb 2003/05/21 pitch floor and ceiling replace the separate ranges for viewing and analysis * pb 2003/05/27 spectrogram maximum and autoscaling * pb 2003/07/20 moved voice report to VoiceAnalysis.c * pb 2003/08/23 reintroduced formant.numberOfTimeSteps * pb 2003/09/11 make sure that analyses objects can be extracted even before first drawing * pb 2003/09/16 advanced pitch settings: pitch.timeStep, pitch.viewFrom, pitch.viewTo, pitch.timeStepsPerView * pb 2003/10/01 time step settings: timeStepStrategy, fixedTimeStep, numberOfTimeStepsPerView * pb 2003/11/30 Sound_to_Spectrogram_windowShapeText * pb 2003/12/09 moved Spectrogram settings back in sync with preferences * pb 2004/04/16 introduced pulses.maximumPeriodFactor * pb 2004/05/12 extended voice report with mean F0 and harmonic-to-noise ratios * pb 2004/07/14 introduced pulses.maximumAmplitudeFactor * pb 2004/10/16 C++ compatible struct tags * pb 2004/10/24 intensity.averagingMethod * pb 2004/10/27 intensity.subtractMeanPressure * pb 2004/11/22 simplified Sound_to_Spectrum () * pb 2004/11/28 improved screen text for unavailable pitch * pb 2004/11/28 warning in settings dialogs for non-standard tiems tep strategies * pb 2005/03/02 all pref string buffers are 260 bytes long * pb 2005/03/07 'intensity' logging sensitive to averaging method * pb 2005/06/16 units * pb 2005/08/18 editor name in log files */ #include #include "FunctionEditor_SoundAnalysis.h" #include "Sound_and_Spectrogram.h" #include "Sound_and_Spectrum.h" #include "Sound_to_Pitch.h" #include "Sound_to_Intensity.h" #include "Sound_to_Formant.h" #include "Pitch_to_PointProcess.h" #include "VoiceAnalysis.h" #include "Resources.h" #include "EditorM.h" #include "praat_script.h" struct logInfo { int toInfoWindow, toLogFile; char fileName [Resources_STRING_BUFFER_SIZE], format [Resources_STRING_BUFFER_SIZE]; }; #if defined (macintosh) #if defined __MACH__ #define LOG_1_FILE_NAME "~/Desktop/Pitch Log" #define LOG_2_FILE_NAME "~/Desktop/Formant Log" #define LOG_3_FILE_NAME "~/Desktop/Log script 3" #define LOG_4_FILE_NAME "~/Desktop/Log script 4" #else #define LOG_1_FILE_NAME "Macintosh HD:Desktop Folder:Pitch Log" #define LOG_2_FILE_NAME "Macintosh HD:Desktop Folder:Formant Log" #define LOG_3_FILE_NAME "Macintosh HD:Desktop Folder:Log script 3" #define LOG_4_FILE_NAME "Macintosh HD:Desktop Folder:Log script 4" #endif #elif defined (WIN32) #define LOG_1_FILE_NAME "C:\\WINDOWS\\DESKTOP\\Pitch Log.txt" #define LOG_2_FILE_NAME "C:\\WINDOWS\\DESKTOP\\Formant Log.txt" #define LOG_3_FILE_NAME "C:\\WINDOWS\\DESKTOP\\Log script 3.praat" #define LOG_4_FILE_NAME "C:\\WINDOWS\\DESKTOP\\Log script 4.praat" #else #define LOG_1_FILE_NAME "~/pitch_log" #define LOG_2_FILE_NAME "~/formant_log" #define LOG_3_FILE_NAME "~/log_script3" #define LOG_4_FILE_NAME "~/log_script4" #endif #define LOG_1_FORMAT "Time 'time:6' seconds, pitch 'f0:2' Hertz" #define LOG_2_FORMAT "'t1:4''tab$''t2:4''tab$''f1:0''tab$''f2:0''tab$''f3:0'" static struct { double longestAnalysis; int timeStepStrategy; double fixedTimeStep; long numberOfTimeStepsPerView; struct FunctionEditor_spectrogram spectrogram; struct FunctionEditor_pitch pitch; struct FunctionEditor_intensity intensity; struct FunctionEditor_formant formant; struct FunctionEditor_pulses pulses; struct logInfo log [2]; char logScript3 [Resources_STRING_BUFFER_SIZE], logScript4 [Resources_STRING_BUFFER_SIZE]; } preferences = { 10.0, 1 /* automatic time step */, 0.01 /* in case of a fixed step */, 100 /* in case of a fixed number per view */, { /* Spectrogram settings: */ NULL /* data */, TRUE /* show */, 0.0, 5000.0, /* Hertz */ 0.005 /* seconds */, 70.0 /* dB */, /* Advanced spectrogram settings: */ 1000, 250, 1 /* Fourier */, 5 /* Gaussian */, TRUE /* autoscaling on */, 100.0 /* dB/Hz */, 6.0 /* dB/octave */, 0.0 /* no dynamic compression */, /* Dynamic information: */ 0.0 }, { NULL, TRUE, /* Pitch settings: */ 75.0, 500.0, Pitch_UNIT_HERTZ, /* Advanced pitch settings: */ 0.0 /* auto view from */, 0.0 /* auto view to */, 1, FALSE, 15, 0.03, 0.45, 0.01, 0.35, 0.14 }, { NULL, FALSE, /* Intensity settings: */ 50.0, 100.0, Intensity_averaging_ENERGY, TRUE }, { NULL, TRUE, 5500.0, 10, 0.025, 30.0, 1.0, 1, 50.0, }, { /* Pulses settings: */ NULL, TRUE, 1.3, /* maximum period factor */ 1.6 /* maximum amplitude factor */ }, { TRUE, TRUE, LOG_1_FILE_NAME, LOG_1_FORMAT, TRUE, TRUE, LOG_2_FILE_NAME, LOG_2_FORMAT, }, LOG_3_FILE_NAME, LOG_4_FILE_NAME }; void FunctionEditor_SoundAnalysis_prefs (void) { Resources_addDouble ("FunctionEditor.longestAnalysis", & preferences.longestAnalysis); Resources_addInt ("FunctionEditor.timeStepStrategy", & preferences.timeStepStrategy); Resources_addDouble ("FunctionEditor.fixedTimeStep", & preferences.fixedTimeStep); Resources_addLong ("FunctionEditor.numberOfTimeStepsPerView", & preferences.numberOfTimeStepsPerView); Resources_addInt ("FunctionEditor.spectrogram.show", & preferences.spectrogram.show); Resources_addLong ("FunctionEditor.spectrogram.timeSteps2", & preferences.spectrogram.timeSteps); Resources_addLong ("FunctionEditor.spectrogram.frequencySteps2", & preferences.spectrogram.frequencySteps); Resources_addDouble ("FunctionEditor.spectrogram.viewFrom2", & preferences.spectrogram.viewFrom); Resources_addDouble ("FunctionEditor.spectrogram.viewTo2", & preferences.spectrogram.viewTo); Resources_addInt ("FunctionEditor.spectrogram.method2", & preferences.spectrogram.method); Resources_addDouble ("FunctionEditor.spectrogram.windowLength2", & preferences.spectrogram.windowLength); Resources_addInt ("FunctionEditor.spectrogram.windowShape2", & preferences.spectrogram.windowShape); Resources_addDouble ("FunctionEditor.spectrogram.maximum2", & preferences.spectrogram.maximum); Resources_addInt ("FunctionEditor.spectrogram.autoscaling2", & preferences.spectrogram.autoscaling); Resources_addDouble ("FunctionEditor.spectrogram.dynamicRange2", & preferences.spectrogram.dynamicRange); Resources_addDouble ("FunctionEditor.spectrogram.preemphasis2", & preferences.spectrogram.preemphasis); Resources_addDouble ("FunctionEditor.spectrogram.dynamicCompression2", & preferences.spectrogram.dynamicCompression); Resources_addInt ("FunctionEditor.pitch.show", & preferences.pitch.show); Resources_addDouble ("FunctionEditor.pitch.floor", & preferences.pitch.floor); Resources_addDouble ("FunctionEditor.pitch.ceiling", & preferences.pitch.ceiling); Resources_addInt ("FunctionEditor.pitch.unit", & preferences.pitch.unit); Resources_addDouble ("FunctionEditor.pitch.viewFrom", & preferences.pitch.viewFrom); Resources_addDouble ("FunctionEditor.pitch.viewTo", & preferences.pitch.viewTo); Resources_addInt ("FunctionEditor.pitch.method", & preferences.pitch.method); Resources_addInt ("FunctionEditor.pitch.veryAccurate", & preferences.pitch.veryAccurate); Resources_addLong ("FunctionEditor.pitch.maximumNumberOfCandidates", & preferences.pitch.maximumNumberOfCandidates); Resources_addDouble ("FunctionEditor.pitch.silenceThreshold", & preferences.pitch.silenceThreshold); Resources_addDouble ("FunctionEditor.pitch.voicingThreshold", & preferences.pitch.voicingThreshold); Resources_addDouble ("FunctionEditor.pitch.octaveCost", & preferences.pitch.octaveCost); Resources_addDouble ("FunctionEditor.pitch.octaveJumpCost", & preferences.pitch.octaveJumpCost); Resources_addDouble ("FunctionEditor.pitch.voicedUnvoicedCost", & preferences.pitch.voicedUnvoicedCost); Resources_addInt ("FunctionEditor.intensity.show", & preferences.intensity.show); Resources_addDouble ("FunctionEditor.intensity.viewFrom", & preferences.intensity.viewFrom); Resources_addDouble ("FunctionEditor.intensity.viewTo", & preferences.intensity.viewTo); Resources_addInt ("FunctionEditor.intensity.averagingMethod", & preferences.intensity.averagingMethod); Resources_addInt ("FunctionEditor.intensity.subtractMeanPressure", & preferences.intensity.subtractMeanPressure); Resources_addInt ("FunctionEditor.formant.show", & preferences.formant.show); Resources_addDouble ("FunctionEditor.formant.maximumFormant", & preferences.formant.maximumFormant); Resources_addLong ("FunctionEditor.formant.numberOfPoles", & preferences.formant.numberOfPoles); Resources_addDouble ("FunctionEditor.formant.windowLength", & preferences.formant.windowLength); Resources_addDouble ("FunctionEditor.formant.dynamicRange", & preferences.formant.dynamicRange); Resources_addDouble ("FunctionEditor.formant.dotSize", & preferences.formant.dotSize); Resources_addInt ("FunctionEditor.formant.method", & preferences.formant.method); Resources_addDouble ("FunctionEditor.formant.preemphasisFrom", & preferences.formant.preemphasisFrom); Resources_addInt ("FunctionEditor.pulses.show", & preferences.pulses.show); Resources_addDouble ("FunctionEditor.pulses.maximumPeriodFactor", & preferences.pulses.maximumPeriodFactor); Resources_addDouble ("FunctionEditor.pulses.maximumAmplitudeFactor", & preferences.pulses.maximumAmplitudeFactor); Resources_addInt ("FunctionEditor.log1.toInfoWindow", & preferences.log[0].toInfoWindow); Resources_addInt ("FunctionEditor.log1.toLogFile", & preferences.log[0].toLogFile); Resources_addString ("FunctionEditor.log1.fileName", & preferences.log[0].fileName [0]); Resources_addString ("FunctionEditor.log1.format", & preferences.log[0].format [0]); Resources_addInt ("FunctionEditor.log2.toInfoWindow", & preferences.log[1].toInfoWindow); Resources_addInt ("FunctionEditor.log2.toLogFile", & preferences.log[1].toLogFile); Resources_addString ("FunctionEditor.log2.fileName", & preferences.log[1].fileName [0]); Resources_addString ("FunctionEditor.log2.format", & preferences.log[1].format [0]); Resources_addString ("FunctionEditor.logScript3", & preferences.logScript3 [0]); Resources_addString ("FunctionEditor.logScript4", & preferences.logScript4 [0]); } void FunctionEditor_SoundAnalysis_forget (I) { iam (FunctionEditor); forget (my spectrogram.data); forget (my pitch.data); forget (my intensity.data); forget (my formant.data); forget (my pulses.data); } static Sound extractSound (FunctionEditor me, double tmin, double tmax) { Sound sound = NULL; if (my longSound.data) { if (tmin < my longSound.data -> xmin) tmin = my longSound.data -> xmin; if (tmax > my longSound.data -> xmax) tmax = my longSound.data -> xmax; sound = LongSound_extractPart (my longSound.data, tmin, tmax, TRUE); } else if (my sound.data) { if (tmin < my sound.data -> xmin) tmin = my sound.data -> xmin; if (tmax > my sound.data -> xmax) tmax = my sound.data -> xmax; sound = Sound_extractPart (my sound.data, tmin, tmax, enumi (Sound_WINDOW, Rectangular), 1.0, TRUE); } return sound; } static void computeSpectrogram (FunctionEditor me) { Melder_progressOff (); if (my spectrogram.show && my endWindow - my startWindow <= my longestAnalysis && (my spectrogram.data == NULL || my spectrogram.data -> xmin != my startWindow || my spectrogram.data -> xmax != my endWindow)) { Sound sound = NULL; double margin = my spectrogram.windowShape == 5 ? my spectrogram.windowLength : 0.5 * my spectrogram.windowLength; forget (my spectrogram.data); sound = extractSound (me, my startWindow - margin, my endWindow + margin); if (sound != NULL) { my spectrogram.data = Sound_to_Spectrogram (sound, my spectrogram.windowLength, my spectrogram.viewTo, (my endWindow - my startWindow) / my spectrogram.timeSteps, my spectrogram.viewTo / my spectrogram.frequencySteps, my spectrogram.windowShape, 8.0, 8.0); if (my spectrogram.data != NULL) my spectrogram.data -> xmin = my startWindow, my spectrogram.data -> xmax = my endWindow; else Melder_clearError (); forget (sound); } else Melder_clearError (); } Melder_progressOn (); } static void computePitch_inside (FunctionEditor me) { Sound sound = NULL; double margin = my pitch.veryAccurate ? 3.0 / my pitch.floor : 1.5 / my pitch.floor; forget (my pitch.data); sound = extractSound (me, my startWindow - margin, my endWindow + margin); if (sound != NULL) { double pitchTimeStep = my timeStepStrategy == 2 ? my fixedTimeStep : my timeStepStrategy == 3 ? (my endWindow - my startWindow) / my numberOfTimeStepsPerView : 0.0; /* The default: determined by pitch floor. */ my pitch.data = Sound_to_Pitch_any (sound, pitchTimeStep, my pitch.floor, my pitch.method == 1 ? 3.0 : 1.0, my pitch.maximumNumberOfCandidates, (my pitch.method - 1) * 2 + my pitch.veryAccurate, my pitch.silenceThreshold, my pitch.voicingThreshold, my pitch.octaveCost, my pitch.octaveJumpCost, my pitch.voicedUnvoicedCost, my pitch.ceiling); if (my pitch.data != NULL) my pitch.data -> xmin = my startWindow, my pitch.data -> xmax = my endWindow; else Melder_clearError (); forget (sound); } else Melder_clearError (); } static void computePitch (FunctionEditor me) { double margin = my pitch.veryAccurate ? 3.0 / my pitch.floor : 1.5 / my pitch.floor; Melder_progressOff (); if (my pitch.show && my endWindow - my startWindow <= my longestAnalysis && (my pitch.data == NULL || my pitch.data -> xmin != my startWindow || my pitch.data -> xmax != my endWindow)) { computePitch_inside (me); } Melder_progressOn (); } static void computeIntensity (FunctionEditor me) { Melder_progressOff (); if (my intensity.show && my endWindow - my startWindow <= my longestAnalysis && (my intensity.data == NULL || my intensity.data -> xmin != my startWindow || my intensity.data -> xmax != my endWindow)) { Sound sound = NULL; double margin = 3.2 / my pitch.floor; forget (my intensity.data); sound = extractSound (me, my startWindow - margin, my endWindow + margin); if (sound != NULL) { my intensity.data = Sound_to_Intensity (sound, my pitch.floor, my endWindow - my startWindow > my longestAnalysis ? (my endWindow - my startWindow) / 100 : 0.0, my intensity.subtractMeanPressure); if (my intensity.data != NULL) my intensity.data -> xmin = my startWindow, my intensity.data -> xmax = my endWindow; else Melder_clearError (); forget (sound); } else Melder_clearError (); } Melder_progressOn (); } static void computeFormants (FunctionEditor me) { Melder_progressOff (); if (my formant.show && my endWindow - my startWindow <= my longestAnalysis && (my formant.data == NULL || my formant.data -> xmin != my startWindow || my formant.data -> xmax != my endWindow)) { Sound sound = NULL; double margin = my formant.windowLength; forget (my formant.data); if (my endWindow - my startWindow > my longestAnalysis) sound = extractSound (me, 0.5 * (my startWindow + my endWindow - my longestAnalysis) - margin, 0.5 * (my startWindow + my endWindow + my longestAnalysis) + margin); else sound = extractSound (me, my startWindow - margin, my endWindow + margin); if (sound != NULL) { double formantTimeStep = my timeStepStrategy == 2 ? my fixedTimeStep : my timeStepStrategy == 3 ? (my endWindow - my startWindow) / my numberOfTimeStepsPerView : 0.0; /* The default: determined by analysis window length. */ my formant.data = Sound_to_Formant_any (sound, formantTimeStep, my formant.numberOfPoles, my formant.maximumFormant, my formant.windowLength, my formant.method, my formant.preemphasisFrom, 50.0); if (my formant.data != NULL) my formant.data -> xmin = my startWindow, my formant.data -> xmax = my endWindow; else Melder_clearError (); forget (sound); } else Melder_clearError (); } Melder_progressOn (); } static void computePulses (FunctionEditor me) { double margin = my pitch.veryAccurate ? 3.0 / my pitch.floor : 1.5 / my pitch.floor; Melder_progressOff (); if (my pulses.show && my endWindow - my startWindow <= my longestAnalysis && (my pulses.data == NULL || my pulses.data -> xmin != my startWindow || my pulses.data -> xmax != my endWindow)) { if (my pitch.data == NULL || my pitch.data -> xmin != my startWindow || my pitch.data -> xmax != my endWindow) { computePitch_inside (me); } if (my pitch.data != NULL) { Sound sound = NULL; forget (my pulses.data); sound = extractSound (me, my startWindow, my endWindow); if (sound != NULL) { my pulses.data = Sound_Pitch_to_PointProcess_cc (sound, my pitch.data); if (my pulses.data == NULL) Melder_clearError (); forget (sound); } else Melder_clearError (); } else Melder_clearError (); } Melder_progressOn (); } #define FunctionEditor_PART_CURSOR 1 #define FunctionEditor_PART_SELECTION 2 static const char *FunctionEditor_partString (int part) { static const char *strings [] = { "", "CURSOR", "SELECTION" }; return strings [part]; } static const char *FunctionEditor_partString_locative (int part) { static const char *strings [] = { "", "at CURSOR", "in SELECTION" }; return strings [part]; } static int makeQueriable (FunctionEditor me, int allowCursor, double *tmin, double *tmax) { if (my endWindow - my startWindow > my longestAnalysis) { return Melder_error ("Window too long to show analyses. Zoom in to at most %s seconds or set the \"longest analysis\" " "to at least %s seconds.", Melder_half (my longestAnalysis), Melder_half (my endWindow - my startWindow)); } if (my startSelection == my endSelection) { if (allowCursor) { *tmin = *tmax = my startSelection; return FunctionEditor_PART_CURSOR; } else { return Melder_error ("Make a selection first."); } } else if (my startSelection < my startWindow || my endSelection > my endWindow) { return Melder_error ("Command ambiguous: a part of the selection is out of view. Either zoom or re-select."); } *tmin = my startSelection; *tmax = my endSelection; return FunctionEditor_PART_SELECTION; } /***** VIEW MENU *****/ FORM (FunctionEditor, cb_showAnalyses, "Show analyses", 0) BOOLEAN ("Show spectrogram", 1) BOOLEAN ("Show pitch", 1) BOOLEAN ("Show intensity", 0) BOOLEAN ("Show formants", 1) BOOLEAN ("Show pulses", 0) POSITIVE ("Longest analysis (s)", "5.0") OK SET_INTEGER ("Show spectrogram", my spectrogram.show) SET_INTEGER ("Show pitch", my pitch.show) SET_INTEGER ("Show intensity", my intensity.show) SET_INTEGER ("Show formants", my formant.show) SET_INTEGER ("Show pulses", my pulses.show) SET_REAL ("Longest analysis", my longestAnalysis) DO XmToggleButtonGadgetSetState (my spectrogramToggle, preferences.spectrogram.show = my spectrogram.show = GET_INTEGER ("Show spectrogram"), False); XmToggleButtonGadgetSetState (my pitchToggle, preferences.pitch.show = my pitch.show = GET_INTEGER ("Show pitch"), False); XmToggleButtonGadgetSetState (my intensityToggle, preferences.intensity.show = my intensity.show = GET_INTEGER ("Show intensity"), False); XmToggleButtonGadgetSetState (my formantToggle, preferences.formant.show = my formant.show = GET_INTEGER ("Show formants"), False); XmToggleButtonGadgetSetState (my pulsesToggle, preferences.pulses.show = my pulses.show = GET_INTEGER ("Show pulses"), False); preferences.longestAnalysis = my longestAnalysis = GET_REAL ("Longest analysis"); FunctionEditor_redraw (me); END FORM (FunctionEditor, cb_timeStepSettings, "Time step settings", "Time step settings...") OPTIONMENU ("Time step strategy", 1) OPTION ("Automatic") OPTION ("Fixed") OPTION ("View-dependent") LABEL ("", "") LABEL ("", "If the time step strategy is \"fixed\":") POSITIVE ("Fixed time step (s)", "0.01") LABEL ("", "") LABEL ("", "If the time step strategy is \"view-dependent\":") NATURAL ("Number of time steps per view", "100") OK SET_INTEGER ("Time step strategy", my timeStepStrategy) SET_REAL ("Fixed time step", my fixedTimeStep) SET_INTEGER ("Number of time steps per view", my numberOfTimeStepsPerView) DO preferences.timeStepStrategy = my timeStepStrategy = GET_INTEGER ("Time step strategy"); preferences.fixedTimeStep = my fixedTimeStep = GET_REAL ("Fixed time step"); preferences.numberOfTimeStepsPerView = my numberOfTimeStepsPerView = GET_INTEGER ("Number of time steps per view"); forget (my pitch.data); forget (my formant.data); forget (my intensity.data); forget (my pulses.data); FunctionEditor_redraw (me); END /***** SPECTROGRAM MENU *****/ DIRECT (FunctionEditor, cb_showSpectrogram) preferences.spectrogram.show = my spectrogram.show = ! my spectrogram.show; FunctionEditor_redraw (me); END FORM (FunctionEditor, cb_spectrogramSettings, "Spectrogram settings", "Intro 3.2. Configuring the spectrogram") REAL ("left View range (Hz)", "0.0") POSITIVE ("right View range (Hz)", "5000.0") POSITIVE ("Window length (s)", "0.005") POSITIVE ("Dynamic range (dB)", "50.0") LABEL ("note1", "") LABEL ("note2", "") OK SET_REAL ("left View range", my spectrogram.viewFrom) SET_REAL ("right View range", my spectrogram.viewTo) SET_REAL ("Window length", my spectrogram.windowLength) SET_REAL ("Dynamic range", my spectrogram.dynamicRange) if (my spectrogram.timeSteps != 1000 || my spectrogram.frequencySteps != 250 || my spectrogram.method != 1 || my spectrogram.windowShape != 5 || my spectrogram.maximum != 100.0 || my spectrogram.autoscaling != TRUE || my spectrogram.preemphasis != 6.0 || my spectrogram.dynamicCompression != 0.0) { SET_STRING ("note1", "Warning: you have non-standard \"advanced settings\".") } else { SET_STRING ("note1", "(all of your \"advanced settings\" have their standard values)") } if (my timeStepStrategy != 1) { SET_STRING ("note2", "Warning: you have a non-standard \"time step strategy\".") } else { SET_STRING ("note2", "(your \"time step strategy\" has its standard value: automatic)") } DO preferences.spectrogram.viewFrom = my spectrogram.viewFrom = GET_REAL ("left View range"); preferences.spectrogram.viewTo = my spectrogram.viewTo = GET_REAL ("right View range"); preferences.spectrogram.windowLength = my spectrogram.windowLength = GET_REAL ("Window length"); preferences.spectrogram.dynamicRange = my spectrogram.dynamicRange = GET_REAL ("Dynamic range"); forget (my spectrogram.data); FunctionEditor_redraw (me); END FORM (FunctionEditor, cb_advancedSpectrogramSettings, "Advanced spectrogram settings", "Advanced spectrogram settings...") LABEL ("", "Time and frequency resolutions:") NATURAL ("Number of time steps", "1000") NATURAL ("Number of frequency steps", "250") LABEL ("", "Spectrogram analysis settings:") OPTIONMENU ("Method", 1) OPTION ("Fourier") OPTIONMENU ("Window shape", 6) { int i; for (i = 0; i < 6; i ++) { OPTION (Sound_to_Spectrogram_windowShapeText (i)) } } LABEL ("", "Spectrogram view settings:") BOOLEAN ("Autoscaling", 1) REAL ("Maximum (dB/Hz)", "100.0") REAL ("Pre-emphasis (dB/oct)", "6.0") REAL ("Dynamic compression (0-1)", "0.0") OK SET_INTEGER ("Number of time steps", my spectrogram.timeSteps) SET_INTEGER ("Number of frequency steps", my spectrogram.frequencySteps) SET_INTEGER ("Method", my spectrogram.method) SET_INTEGER ("Window shape", my spectrogram.windowShape + 1) SET_REAL ("Maximum", my spectrogram.maximum) SET_INTEGER ("Autoscaling", my spectrogram.autoscaling) SET_REAL ("Pre-emphasis", my spectrogram.preemphasis) SET_REAL ("Dynamic compression", my spectrogram.dynamicCompression) DO preferences.spectrogram.timeSteps = my spectrogram.timeSteps = GET_INTEGER ("Number of time steps"); preferences.spectrogram.frequencySteps = my spectrogram.frequencySteps = GET_INTEGER ("Number of frequency steps"); preferences.spectrogram.method = my spectrogram.method = GET_INTEGER ("Method"); preferences.spectrogram.windowShape = my spectrogram.windowShape = GET_INTEGER ("Window shape") - 1; preferences.spectrogram.maximum = my spectrogram.maximum = GET_REAL ("Maximum"); preferences.spectrogram.autoscaling = my spectrogram.autoscaling = GET_INTEGER ("Autoscaling"); preferences.spectrogram.preemphasis = my spectrogram.preemphasis = GET_REAL ("Pre-emphasis"); preferences.spectrogram.dynamicCompression = my spectrogram.dynamicCompression = GET_REAL ("Dynamic compression"); forget (my spectrogram.data); FunctionEditor_redraw (me); END DIRECT (FunctionEditor, cb_extractVisibleSpectrogram) Spectrogram publish; if (! my spectrogram.show) return Melder_error ("No spectrogram is visible.\nFirst choose \"Show spectrogram\" from the Spectrum menu."); if (! my spectrogram.data) { computeSpectrogram (me); if (! my spectrogram.data) return Melder_error ("No spectrogram available (out of memory?)."); } publish = Data_copy (my spectrogram.data); if (publish == NULL) return 0; if (my publishCallback) my publishCallback (me, my publishClosure, publish); END DIRECT (FunctionEditor, cb_viewSpectralSlice) double start = my startSelection == my endSelection ? my spectrogram.windowShape == 5 ? my startSelection - my spectrogram.windowLength : my startSelection - my spectrogram.windowLength / 2 : my startSelection; double finish = my startSelection == my endSelection ? my spectrogram.windowShape == 5 ? my endSelection + my spectrogram.windowLength : my endSelection + my spectrogram.windowLength / 2 : my endSelection; Sound sound = extractSound (me, start, finish); Spectrum publish; if (sound == NULL) return 0; Sound_multiplyByWindow (sound, my spectrogram.windowShape == 0 ? enumi (Sound_WINDOW, Rectangular) : my spectrogram.windowShape == 1 ? enumi (Sound_WINDOW, Hamming) : my spectrogram.windowShape == 2 ? enumi (Sound_WINDOW, Triangular) : my spectrogram.windowShape == 3 ? enumi (Sound_WINDOW, Parabolic) : my spectrogram.windowShape == 4 ? enumi (Sound_WINDOW, Hanning) : my spectrogram.windowShape == 5 ? enumi (Sound_WINDOW, Gaussian2) : 0); publish = Sound_to_Spectrum (sound, TRUE); forget (sound); if (! publish) return 0; Thing_setName (publish, "slice"); if (my publishCallback) my publishCallback (me, my publishClosure, publish); END DIRECT (FunctionEditor, cb_getFrequency) Melder_informationReal (my spectrogram.cursor, "Hertz"); END DIRECT (FunctionEditor, cb_getSpectralPowerAtCursorCross) double tmin, tmax; int part = makeQueriable (me, TRUE, & tmin, & tmax); iferror return 0; if (! my spectrogram.show) return Melder_error ("No spectrogram is visible.\nFirst choose \"Show spectrogram\" from the Spectrum menu."); if (! my spectrogram.data) { computeSpectrogram (me); if (! my spectrogram.data) return Melder_error ("No spectrogram available (out of memory?)."); } if (part != FunctionEditor_PART_CURSOR) return Melder_error ("Click inside the spectrogram first."); MelderInfo_open (); MelderInfo_write1 (Melder_double (Matrix_getValueAtXY (my spectrogram.data, tmin, my spectrogram.cursor))); MelderInfo_write5 (" Pa2/Hz (at time = ", Melder_double (tmin), " seconds and frequency = ", Melder_double (my spectrogram.cursor), " Hz)"); MelderInfo_close (); END /***** PITCH MENU *****/ DIRECT (FunctionEditor, cb_showPitch) preferences.pitch.show = my pitch.show = ! my pitch.show; FunctionEditor_redraw (me); END FORM (FunctionEditor, cb_pitchSettings, "Pitch settings", "Intro 4.2. Configuring the pitch contour") POSITIVE ("left Pitch range (Hz)", "75.0") POSITIVE ("right Pitch range (Hz)", "500.0") OPTIONMENU ("Unit", 1) OPTIONS_ENUM (ClassFunction_getUnitText (classPitch, Pitch_LEVEL_FREQUENCY, itext, Function_UNIT_TEXT_MENU), Pitch_UNIT_min, Pitch_UNIT_max) RADIO ("Optimize for", 1) RADIOBUTTON ("Intonation (AC method)") RADIOBUTTON ("Voice analysis (CC method)") LABEL ("note1", "") LABEL ("note2", "") OK SET_REAL ("left Pitch range", my pitch.floor) SET_REAL ("right Pitch range", my pitch.ceiling) SET_INTEGER ("Unit", my pitch.unit + 1 - Pitch_UNIT_min) SET_INTEGER ("Optimize for", my pitch.method) if (my pitch.viewFrom != 0.0 || my pitch.viewTo != 0.0 || my pitch.veryAccurate != FALSE || my pitch.maximumNumberOfCandidates != 15 || my pitch.silenceThreshold != 0.03 || my pitch.voicingThreshold != 0.45 || my pitch.octaveCost != 0.01 || my pitch.octaveJumpCost != 0.35 || my pitch.voicedUnvoicedCost != 0.14) { SET_STRING ("note1", "Warning: you have some non-standard \"advanced settings\".") } else { SET_STRING ("note1", "(all of your \"advanced settings\" have their standard values)") } if (my timeStepStrategy != 1) { SET_STRING ("note2", "Warning: you have a non-standard \"time step strategy\".") } else { SET_STRING ("note2", "(your \"time step strategy\" has its standard value: automatic)") } DO preferences.pitch.floor = my pitch.floor = GET_REAL ("left Pitch range"); preferences.pitch.ceiling = my pitch.ceiling = GET_REAL ("right Pitch range"); preferences.pitch.unit = my pitch.unit = GET_INTEGER ("Unit") - 1 + Pitch_UNIT_min; preferences.pitch.method = my pitch.method = GET_INTEGER ("Optimize for"); forget (my pitch.data); forget (my intensity.data); forget (my pulses.data); FunctionEditor_redraw (me); END FORM (FunctionEditor, cb_advancedPitchSettings, "Advanced pitch settings", "Advanced pitch settings...") LABEL ("", "Make view range different from analysis range:") REAL ("left View range (units)", "0.0 (= auto)") REAL ("right View range (units)", "0.0 (= auto)") LABEL ("", "Analysis settings:") BOOLEAN ("Very accurate", 0) NATURAL ("Max. number of candidates", "15") REAL ("Silence threshold", "0.03") REAL ("Voicing threshold", "0.45") REAL ("Octave cost", "0.01") REAL ("Octave-jump cost", "0.35") REAL ("Voiced / unvoiced cost", "0.14") OK SET_REAL ("left View range", my pitch.viewFrom) SET_REAL ("right View range", my pitch.viewTo) SET_INTEGER ("Very accurate", my pitch.veryAccurate) SET_INTEGER ("Max. number of candidates", my pitch.maximumNumberOfCandidates) SET_REAL ("Silence threshold", my pitch.silenceThreshold) SET_REAL ("Voicing threshold", my pitch.voicingThreshold) SET_REAL ("Octave cost", my pitch.octaveCost) SET_REAL ("Octave-jump cost", my pitch.octaveJumpCost) SET_REAL ("Voiced / unvoiced cost", my pitch.voicedUnvoicedCost) DO long maxnCandidates = GET_INTEGER ("Max. number of candidates"); if (maxnCandidates < 2) return Melder_error ("Maximum number of candidates must be greater than 1."); preferences.pitch.viewFrom = my pitch.viewFrom = GET_REAL ("left View range"); preferences.pitch.viewTo = my pitch.viewTo = GET_REAL ("right View range"); preferences.pitch.veryAccurate = my pitch.veryAccurate = GET_INTEGER ("Very accurate"); preferences.pitch.maximumNumberOfCandidates = my pitch.maximumNumberOfCandidates = GET_INTEGER ("Max. number of candidates"); preferences.pitch.silenceThreshold = my pitch.silenceThreshold = GET_REAL ("Silence threshold"); preferences.pitch.voicingThreshold = my pitch.voicingThreshold = GET_REAL ("Voicing threshold"); preferences.pitch.octaveCost = my pitch.octaveCost = GET_REAL ("Octave cost"); preferences.pitch.octaveJumpCost = my pitch.octaveJumpCost = GET_REAL ("Octave-jump cost"); preferences.pitch.voicedUnvoicedCost = my pitch.voicedUnvoicedCost = GET_REAL ("Voiced / unvoiced cost"); forget (my pitch.data); forget (my intensity.data); forget (my pulses.data); FunctionEditor_redraw (me); END DIRECT (FunctionEditor, cb_pitchListing) double tmin, tmax; int part = makeQueriable (me, TRUE, & tmin, & tmax); iferror return 0; if (! my pitch.show) return Melder_error ("No pitch contour is visible.\nFirst choose \"Show pitch\" from the Pitch menu."); if (! my pitch.data) { computePitch (me); if (! my pitch.data) return Melder_error ("No pitch contour available (out of memory?)."); } MelderInfo_open (); MelderInfo_writeLine2 ("Time_s F0_", ClassFunction_getUnitText (classPitch, Pitch_LEVEL_FREQUENCY, my pitch.unit, Function_UNIT_TEXT_SHORT)); if (part == FunctionEditor_PART_CURSOR) { double f0 = Pitch_getValueAtTime (my pitch.data, tmin, my pitch.unit, TRUE); f0 = ClassFunction_convertToNonlogarithmic (classPitch, f0, Pitch_LEVEL_FREQUENCY, my pitch.unit); MelderInfo_writeLine3 (Melder_fixed (tmin, 6), " ", Melder_fixed (f0, 6)); } else { long i, i1, i2; Sampled_getWindowSamples (my pitch.data, tmin, tmax, & i1, & i2); for (i = i1; i <= i2; i ++) { double t = Sampled_indexToX (my pitch.data, i); double f0 = Sampled_getValueAtSample (my pitch.data, i, Pitch_LEVEL_FREQUENCY, my pitch.unit); f0 = ClassFunction_convertToNonlogarithmic (classPitch, f0, Pitch_LEVEL_FREQUENCY, my pitch.unit); MelderInfo_writeLine3 (Melder_fixed (t, 6), " ", Melder_fixed (f0, 6)); } } MelderInfo_close (); return 1; END DIRECT (FunctionEditor, cb_getPitch) double tmin, tmax; int part = makeQueriable (me, TRUE, & tmin, & tmax); iferror return 0; if (! my pitch.show) return Melder_error ("No pitch contour is visible.\nFirst choose \"Show pitch\" from the Pitch menu."); if (! my pitch.data) { computePitch (me); if (! my pitch.data) return Melder_error ("No pitch contour available (out of memory?)."); } if (part == FunctionEditor_PART_CURSOR) { double f0 = Pitch_getValueAtTime (my pitch.data, tmin, my pitch.unit, TRUE); f0 = ClassFunction_convertToNonlogarithmic (classPitch, f0, Pitch_LEVEL_FREQUENCY, my pitch.unit); Melder_information ("%s %s (interpolated pitch at CURSOR)", Melder_double (f0), ClassFunction_getUnitText (classPitch, Pitch_LEVEL_FREQUENCY, my pitch.unit, 0)); } else { double f0 = Pitch_getMean (my pitch.data, tmin, tmax, my pitch.unit); f0 = ClassFunction_convertToNonlogarithmic (classPitch, f0, Pitch_LEVEL_FREQUENCY, my pitch.unit); Melder_information ("%s %s (mean pitch %s)", Melder_double (f0), ClassFunction_getUnitText (classPitch, Pitch_LEVEL_FREQUENCY, my pitch.unit, 0), FunctionEditor_partString_locative (part)); } END DIRECT (FunctionEditor, cb_getMinimumPitch) double tmin, tmax, f0; int part = makeQueriable (me, FALSE, & tmin, & tmax); iferror return 0; if (! my pitch.show) return Melder_error ("No pitch contour is visible.\nFirst choose \"Show pitch\" from the Pitch menu."); if (! my pitch.data) { computePitch (me); if (! my pitch.data) return Melder_error ("No pitch contour available (out of memory?)."); } f0 = Pitch_getMinimum (my pitch.data, tmin, tmax, my pitch.unit, TRUE); f0 = ClassFunction_convertToNonlogarithmic (classPitch, f0, Pitch_LEVEL_FREQUENCY, my pitch.unit); Melder_information ("%s %s (minimum pitch %s)", Melder_double (f0), ClassFunction_getUnitText (classPitch, Pitch_LEVEL_FREQUENCY, my pitch.unit, 0), FunctionEditor_partString_locative (part)); END DIRECT (FunctionEditor, cb_getMaximumPitch) double tmin, tmax, f0; int part = makeQueriable (me, FALSE, & tmin, & tmax); iferror return 0; if (! my pitch.show) return Melder_error ("No pitch contour is visible.\nFirst choose \"Show pitch\" from the Pitch menu."); if (! my pitch.data) { computePitch (me); if (! my pitch.data) return Melder_error ("No pitch contour available (out of memory?)."); } f0 = Pitch_getMaximum (my pitch.data, tmin, tmax, my pitch.unit, TRUE); f0 = ClassFunction_convertToNonlogarithmic (classPitch, f0, Pitch_LEVEL_FREQUENCY, my pitch.unit); Melder_information ("%s %s (maximum pitch %s)", Melder_double (f0), ClassFunction_getUnitText (classPitch, Pitch_LEVEL_FREQUENCY, my pitch.unit, 0), FunctionEditor_partString_locative (part)); END DIRECT (FunctionEditor, cb_extractVisiblePitchContour) Pitch publish; if (! my pitch.show) return Melder_error ("No pitch contour is visible.\nFirst choose \"Show pitch\" from the Pitch menu."); if (! my pitch.data) { computePitch (me); if (! my pitch.data) return Melder_error ("No pitch contour available (out of memory?)."); } publish = Data_copy (my pitch.data); if (! publish) return 0; if (my publishCallback) my publishCallback (me, my publishClosure, publish); END /***** INTENSITY MENU *****/ DIRECT (FunctionEditor, cb_showIntensity) preferences.intensity.show = my intensity.show = ! my intensity.show; FunctionEditor_redraw (me); END FORM (FunctionEditor, cb_intensitySettings, "Intensity settings", "Intro 6.2. Configuring the intensity contour") REAL ("left View range (dB)", "50.0") REAL ("right View range (dB)", "100.0") RADIO ("Averaging method", Intensity_averaging_ENERGY + 1) RADIOBUTTON ("median") RADIOBUTTON ("mean energy") RADIOBUTTON ("mean sones") RADIOBUTTON ("mean dB") BOOLEAN ("Subtract mean pressure", 1) LABEL ("", "Note: the pitch floor is taken from the pitch settings.") LABEL ("note2", "") OK SET_REAL ("left View range", my intensity.viewFrom) SET_REAL ("right View range", my intensity.viewTo) SET_INTEGER ("Averaging method", my intensity.averagingMethod + 1) SET_INTEGER ("Subtract mean pressure", my intensity.subtractMeanPressure) if (my timeStepStrategy != 1) { SET_STRING ("note2", "Warning: you have a non-standard \"time step strategy\".") } else { SET_STRING ("note2", "(your \"time step strategy\" has its standard value: automatic)") } DO preferences.intensity.viewFrom = my intensity.viewFrom = GET_REAL ("left View range"); preferences.intensity.viewTo = my intensity.viewTo = GET_REAL ("right View range"); preferences.intensity.averagingMethod = my intensity.averagingMethod = GET_INTEGER ("Averaging method") - 1; preferences.intensity.subtractMeanPressure = my intensity.subtractMeanPressure = GET_INTEGER ("Subtract mean pressure"); forget (my intensity.data); FunctionEditor_redraw (me); END DIRECT (FunctionEditor, cb_extractVisibleIntensityContour) Intensity publish; if (! my intensity.show) return Melder_error ("No intensity contour is visible.\nFirst choose \"Show intensity\" from the Intensity menu."); if (! my intensity.data) { computeIntensity (me); if (! my intensity.data) return Melder_error ("No intensity contour available (out of memory?)."); } publish = Data_copy (my intensity.data); if (! publish) return 0; if (my publishCallback) my publishCallback (me, my publishClosure, publish); END DIRECT (FunctionEditor, cb_intensityListing) double tmin, tmax; int part = makeQueriable (me, TRUE, & tmin, & tmax); iferror return 0; if (! my intensity.show) return Melder_error ("No intensity contour is visible.\nFirst choose \"Show intensity\" from the Intensity menu."); if (! my intensity.data) { computeIntensity (me); if (! my intensity.data) return Melder_error ("No intensity contour available (out of memory?)."); } MelderInfo_open (); MelderInfo_writeLine1 ("Time_s Intensity_dB"); if (part == FunctionEditor_PART_CURSOR) { double intensity = Vector_getValueAtX (my intensity.data, tmin, 1); MelderInfo_writeLine3 (Melder_fixed (tmin, 6), " ", Melder_fixed (intensity, 6)); } else { long i, i1, i2; Sampled_getWindowSamples (my intensity.data, tmin, tmax, & i1, & i2); for (i = i1; i <= i2; i ++) { double t = Sampled_indexToX (my intensity.data, i); double intensity = Vector_getValueAtX (my intensity.data, t, 0); MelderInfo_writeLine3 (Melder_fixed (t, 6), " ", Melder_fixed (intensity, 6)); } } MelderInfo_close (); return 1; END DIRECT (FunctionEditor, cb_getIntensity) double tmin, tmax; int part = makeQueriable (me, TRUE, & tmin, & tmax); iferror return 0; if (! my intensity.show) return Melder_error ("No intensity contour is visible.\nFirst choose \"Show intensity\" from the Intensity menu."); if (! my intensity.data) { computeIntensity (me); if (! my intensity.data) return Melder_error ("No intensity contour available (out of memory?)."); } if (part == FunctionEditor_PART_CURSOR) { Melder_information ("%s dB (intensity at CURSOR)", Melder_double (Vector_getValueAtX (my intensity.data, tmin, 1))); } else { static const char *methodString [] = { "median", "mean-energy", "mean-sones", "mean-dB" }; Melder_information ("%s dB (%s intensity %s)", Melder_double (Intensity_getAverage (my intensity.data, tmin, tmax, my intensity.averagingMethod)), methodString [my intensity.averagingMethod], FunctionEditor_partString_locative (part)); } END /***** FORMANT MENU *****/ DIRECT (FunctionEditor, cb_showFormants) preferences.formant.show = my formant.show = ! my formant.show; FunctionEditor_redraw (me); END FORM (FunctionEditor, cb_formantSettings, "Formant settings", "Intro 5.2. Configuring the formant contours") POSITIVE ("Maximum formant (Hz)", "5500.0") POSITIVE ("Number of formants", "5.0") POSITIVE ("Window length (s)", "0.025") REAL ("Dynamic range (dB)", "30.0") POSITIVE ("Dot size (mm)", "1.0") LABEL ("note1", "") LABEL ("note2", "") OK SET_REAL ("Maximum formant", my formant.maximumFormant) SET_REAL ("Number of formants", 0.5 * my formant.numberOfPoles) SET_REAL ("Window length", my formant.windowLength) SET_REAL ("Dynamic range", my formant.dynamicRange) SET_REAL ("Dot size", my formant.dotSize) if (my formant.method != 1 || my formant.preemphasisFrom != 50.0) { SET_STRING ("note1", "Warning: you have non-standard \"advanced settings\".") } else { SET_STRING ("note1", "(all of your \"advanced settings\" have their standard values)") } if (my timeStepStrategy != 1) { SET_STRING ("note2", "Warning: you have a non-standard \"time step strategy\".") } else { SET_STRING ("note2", "(your \"time step strategy\" has its standard value: automatic)") } DO preferences.formant.maximumFormant = my formant.maximumFormant = GET_REAL ("Maximum formant"); preferences.formant.numberOfPoles = my formant.numberOfPoles = 2.0 * GET_REAL ("Number of formants"); preferences.formant.windowLength = my formant.windowLength = GET_REAL ("Window length"); preferences.formant.dynamicRange = my formant.dynamicRange = GET_REAL ("Dynamic range"); preferences.formant.dotSize = my formant.dotSize = GET_REAL ("Dot size"); forget (my formant.data); FunctionEditor_redraw (me); END FORM (FunctionEditor, cb_advancedFormantSettings, "Advanced formant settings", "Advanced formant settings...") RADIO ("Method", 1) RADIOBUTTON ("Burg") POSITIVE ("Pre-emphasis from (Hz)", "50.0") OK SET_INTEGER ("Method", my formant.method) SET_REAL ("Pre-emphasis from", my formant.preemphasisFrom) DO preferences.formant.method = my formant.method = GET_INTEGER ("Method"); preferences.formant.preemphasisFrom = my formant.preemphasisFrom = GET_REAL ("Pre-emphasis from"); forget (my formant.data); FunctionEditor_redraw (me); END DIRECT (FunctionEditor, cb_extractVisibleFormantContour) Formant publish; if (! my formant.show) return Melder_error ("No formant contour is visible.\nFirst choose \"Show formants\" from the Formant menu."); if (! my formant.data) { computeFormants (me); if (! my formant.data) return Melder_error ("No formant contour available (out of memory?)."); } publish = Data_copy (my formant.data); if (! publish) return 0; if (my publishCallback) my publishCallback (me, my publishClosure, publish); END DIRECT (FunctionEditor, cb_formantListing) double tmin, tmax; int part = makeQueriable (me, TRUE, & tmin, & tmax); iferror return 0; if (! my formant.show) return Melder_error ("No formant contour is visible.\nFirst choose \"Show formants\" from the Formant menu."); if (! my formant.data) { computeFormants (me); if (! my formant.data) return Melder_error ("No formant contour available (out of memory?)."); } MelderInfo_open (); MelderInfo_writeLine1 ("Time_s F1_Hz F2_Hz F3_Hz F4_Hz"); if (part == FunctionEditor_PART_CURSOR) { double f1 = Formant_getValueAtTime (my formant.data, 1, tmin, 0); double f2 = Formant_getValueAtTime (my formant.data, 2, tmin, 0); double f3 = Formant_getValueAtTime (my formant.data, 3, tmin, 0); double f4 = Formant_getValueAtTime (my formant.data, 4, tmin, 0); MelderInfo_write5 (Melder_fixed (tmin, 6), " ", Melder_fixed (f1, 6), " ", Melder_fixed (f2, 6)); MelderInfo_writeLine4 (" ", Melder_fixed (f3, 6), " ", Melder_fixed (f4, 6)); } else { long i, i1, i2; Sampled_getWindowSamples (my formant.data, tmin, tmax, & i1, & i2); for (i = i1; i <= i2; i ++) { double t = Sampled_indexToX (my formant.data, i); double f1 = Formant_getValueAtTime (my formant.data, 1, t, 0); double f2 = Formant_getValueAtTime (my formant.data, 2, t, 0); double f3 = Formant_getValueAtTime (my formant.data, 3, t, 0); double f4 = Formant_getValueAtTime (my formant.data, 4, t, 0); MelderInfo_write5 (Melder_fixed (t, 6), " ", Melder_fixed (f1, 6), " ", Melder_fixed (f2, 6)); MelderInfo_writeLine4 (" ", Melder_fixed (f3, 6), " ", Melder_fixed (f4, 6)); } } MelderInfo_close (); return 1; END static int getFormant (FunctionEditor me, int iformant) { double tmin, tmax; int part = makeQueriable (me, TRUE, & tmin, & tmax); iferror return 0; if (! my formant.show) return Melder_error ("No formant contour is visible.\nFirst choose \"Show formants\" from the Formant menu."); if (! my formant.data) { computeFormants (me); if (! my formant.data) return Melder_error ("No formant contour available (out of memory?)."); } if (part == FunctionEditor_PART_CURSOR) { Melder_information ("%s Hertz (nearest F%d to CURSOR)", Melder_double (Formant_getValueAtTime (my formant.data, iformant, tmin, 0)), iformant); } else { Melder_information ("%s Hertz (mean F%d %s)", Melder_double (Formant_getMean (my formant.data, iformant, tmin, tmax, 0)), iformant, FunctionEditor_partString_locative (part)); } return 1; } static int getBandwidth (FunctionEditor me, int iformant) { double tmin, tmax; int part = makeQueriable (me, TRUE, & tmin, & tmax); iferror return 0; if (! my formant.show) return Melder_error ("No formant contour is visible.\nFirst choose \"Show formants\" from the Formant menu."); if (! my formant.data) { computeFormants (me); if (! my formant.data) return Melder_error ("No formant contour available (out of memory?)."); } if (part == FunctionEditor_PART_CURSOR) { Melder_information ("%s Hertz (nearest B%d to CURSOR)", Melder_double (Formant_getBandwidthAtTime (my formant.data, iformant, tmin, 0)), iformant); } else { Melder_information ("%s Hertz (B%d in centre of %s)", Melder_double (Formant_getBandwidthAtTime (my formant.data, iformant, 0.5 * (tmin + tmax), 0)), iformant, FunctionEditor_partString (part)); } return 1; } DIRECT (FunctionEditor, cb_getFirstFormant) if (! getFormant (me, 1)) return 0; END DIRECT (FunctionEditor, cb_getFirstBandwidth) if (! getBandwidth (me, 1)) return 0; END DIRECT (FunctionEditor, cb_getSecondFormant) if (! getFormant (me, 2)) return 0; END DIRECT (FunctionEditor, cb_getSecondBandwidth) if (! getBandwidth (me, 2)) return 0; END DIRECT (FunctionEditor, cb_getThirdFormant) if (! getFormant (me, 3)) return 0; END DIRECT (FunctionEditor, cb_getThirdBandwidth) if (! getBandwidth (me, 3)) return 0; END DIRECT (FunctionEditor, cb_getFourthFormant) if (! getFormant (me, 4)) return 0; END DIRECT (FunctionEditor, cb_getFourthBandwidth) if (! getBandwidth (me, 4)) return 0; END FORM (FunctionEditor, cb_getFormant, "Get formant", 0) NATURAL ("Formant number", "5") OK DO if (! getFormant (me, GET_INTEGER ("Formant number"))) return 0; END FORM (FunctionEditor, cb_getBandwidth, "Get bandwidth", 0) NATURAL ("Formant number", "5") OK DO if (! getBandwidth (me, GET_INTEGER ("Formant number"))) return 0; END /***** PULSE MENU *****/ DIRECT (FunctionEditor, cb_showPulses) preferences.pulses.show = my pulses.show = ! my pulses.show; FunctionEditor_redraw (me); END FORM (FunctionEditor, cb_advancedPulsesSettings, "Advanced pulses settings", "Advanced pulses settings...") POSITIVE ("Maximum period factor", "1.3") POSITIVE ("Maximum amplitude factor", "1.6") OK SET_REAL ("Maximum period factor", my pulses.maximumPeriodFactor) SET_REAL ("Maximum amplitude factor", my pulses.maximumAmplitudeFactor) DO preferences.pulses.maximumPeriodFactor = my pulses.maximumPeriodFactor = GET_REAL ("Maximum period factor"); preferences.pulses.maximumAmplitudeFactor = my pulses.maximumAmplitudeFactor = GET_REAL ("Maximum amplitude factor"); forget (my pulses.data); FunctionEditor_redraw (me); END DIRECT (FunctionEditor, cb_extractVisiblePulses) Pitch publish; if (! my pulses.show) return Melder_error ("No pulses are visible.\nFirst choose \"Show pulses\" from the Pulses menu."); if (! my pulses.data) { computePulses (me); if (! my pulses.data) return Melder_error ("No pulses available (out of memory?)."); } publish = Data_copy (my pulses.data); if (! publish) return 0; if (my publishCallback) my publishCallback (me, my publishClosure, publish); END DIRECT (FunctionEditor, cb_voiceReport) time_t today = time (NULL); Sound sound = NULL; double tmin, tmax; int part = makeQueriable (me, FALSE, & tmin, & tmax); iferror return 0; if (! my pulses.show) return Melder_error ("No pulses are visible.\nFirst choose \"Show pulses\" from the Pulses menu."); if (! my pulses.data) { computePulses (me); if (! my pulses.data) return Melder_error ("No pulses available (out of memory?)."); } sound = extractSound (me, tmin, tmax); if (! sound) return Melder_error ("Selection too small (or out of memory)."); MelderInfo_open (); MelderInfo_writeLine4 ("-- Voice report for ", my name, " --\nDate: ", ctime (& today)); if ((my pitch.method & 2) == 0) MelderInfo_writeLine1 ("WARNING: some of the following measurements may be imprecise.\n" "For more precision, go to \"Pitch settings\" and choose \"Optimize for voice analysis\".\n"); MelderInfo_writeLine2 ("Time range of ", FunctionEditor_partString (part)); Sound_Pitch_PointProcess_voiceReport (sound, my pitch.data, my pulses.data, tmin, tmax, my pitch.floor, my pitch.ceiling, my pulses.maximumPeriodFactor, my pulses.maximumAmplitudeFactor, my pitch.silenceThreshold, my pitch.voicingThreshold); MelderInfo_close (); forget (sound); END DIRECT (FunctionEditor, cb_pulseListing) long i, i1, i2; double tmin, tmax; int part = makeQueriable (me, FALSE, & tmin, & tmax); iferror return 0; if (! my pulses.show) return Melder_error ("No pulses are visible.\nFirst choose \"Show pulses\" from the Pulses menu."); if (! my pulses.data) { computePulses (me); if (! my pulses.data) return Melder_error ("No pulses available (out of memory?)."); } MelderInfo_open (); MelderInfo_writeLine1 ("Time_s"); i1 = PointProcess_getHighIndex (my pulses.data, tmin); i2 = PointProcess_getLowIndex (my pulses.data, tmax); for (i = i1; i <= i2; i ++) { double t = my pulses.data -> t [i]; MelderInfo_writeLine1 (Melder_fixed (t, 12)); } MelderInfo_close (); return 1; END /* static int cb_getJitter_xx (FunctionEditor me, double (*PointProcess_getJitter_xx) (PointProcess, double, double, double, double, double)) { double minimumPeriod = 0.8 / my pitch.ceiling, maximumPeriod = 1.25 / my pitch.floor; if (! my pulses.show) return Melder_error ("No pulses are visible.\nFirst choose \"Show pulses\" from the Pulses menu."); if (! my pulses.data) { computePulses (me); if (! my pulses.data) return Melder_error ("No pulses available (out of memory?)."); } if (my startSelection == my endSelection) return Melder_error ("Make a selection first."); if (! queriable (me)) return 0; Melder_informationReal (PointProcess_getJitter_xx (my pulses.data, my startSelection, my endSelection, minimumPeriod, maximumPeriod, my pulses.maximumPeriodFactor), NULL); return 1; } DIRECT (FunctionEditor, cb_getJitter_local) if (! cb_getJitter_xx (me, PointProcess_getJitter_local)) return 0; END DIRECT (FunctionEditor, cb_getJitter_local_absolute) if (! cb_getJitter_xx (me, PointProcess_getJitter_local_absolute)) return 0; END DIRECT (FunctionEditor, cb_getJitter_rap) if (! cb_getJitter_xx (me, PointProcess_getJitter_rap)) return 0; END DIRECT (FunctionEditor, cb_getJitter_ppq5) if (! cb_getJitter_xx (me, PointProcess_getJitter_ppq5)) return 0; END DIRECT (FunctionEditor, cb_getJitter_ddp) if (! cb_getJitter_xx (me, PointProcess_getJitter_ddp)) return 0; END static int cb_getShimmer_xx (FunctionEditor me, double (*PointProcess_Sound_getShimmer_xx) (PointProcess, Sound, double, double, double, double, double)) { Sound sound = NULL; double minimumPeriod = 0.8 / my pitch.ceiling, maximumPeriod = 1.25 / my pitch.floor; if (! my pulses.show) return Melder_error ("No pulses are visible.\nFirst choose \"Show pulses\" from the Pulses menu."); if (! my pulses.data) { computePulses (me); if (! my pulses.data) return Melder_error ("No pulses available (out of memory?)."); } if (my startSelection == my endSelection) return Melder_error ("Make a selection first."); if (! queriable (me)) return 0; sound = extractSound (me, my startSelection, my endSelection); if (! sound) return Melder_error ("Selection too small (or out of memory)."); Melder_informationReal (PointProcess_Sound_getShimmer_xx (my pulses.data, sound, my startSelection, my endSelection, minimumPeriod, maximumPeriod, my pulses.maximumAmplitudeFactor), NULL); forget (sound); return 1; } DIRECT (FunctionEditor, cb_getShimmer_local) if (! cb_getShimmer_xx (me, PointProcess_Sound_getShimmer_local)) return 0; END DIRECT (FunctionEditor, cb_getShimmer_local_dB) if (! cb_getShimmer_xx (me, PointProcess_Sound_getShimmer_local_dB)) return 0; END DIRECT (FunctionEditor, cb_getShimmer_apq3) if (! cb_getShimmer_xx (me, PointProcess_Sound_getShimmer_apq3)) return 0; END DIRECT (FunctionEditor, cb_getShimmer_apq5) if (! cb_getShimmer_xx (me, PointProcess_Sound_getShimmer_apq5)) return 0; END DIRECT (FunctionEditor, cb_getShimmer_apq11) if (! cb_getShimmer_xx (me, PointProcess_Sound_getShimmer_apq11)) return 0; END DIRECT (FunctionEditor, cb_getShimmer_dda) if (! cb_getShimmer_xx (me, PointProcess_Sound_getShimmer_dda)) return 0; END */ /***** SELECT MENU *****/ DIRECT (FunctionEditor, cb_moveCursorToMinimumPitch) if (! my pitch.show) return Melder_error ("No pitch contour is visible.\nFirst choose \"Show pitch\" from the View menu."); if (! my pitch.data) { computePitch (me); if (! my pitch.data) return Melder_error ("No pitch contour available (out of memory?)."); } if (my startSelection == my endSelection) { return Melder_error ("Empty selection."); } else { double time; Pitch_getMinimumAndTime (my pitch.data, my startSelection, my endSelection, my pitch.unit, 1, NULL, & time); if (! NUMdefined (time)) return Melder_error ("Selection is voiceless."); my startSelection = my endSelection = time; FunctionEditor_marksChanged (me); } END DIRECT (FunctionEditor, cb_moveCursorToMaximumPitch) if (! my pitch.show) return Melder_error ("No pitch contour is visible.\nFirst choose \"Show pitch\" from the View menu."); if (! my pitch.data) { computePitch (me); if (! my pitch.data) return Melder_error ("No pitch contour available (out of memory?)."); } if (my startSelection == my endSelection) { return Melder_error ("Empty selection."); } else { double time; Pitch_getMaximumAndTime (my pitch.data, my startSelection, my endSelection, my pitch.unit, 1, NULL, & time); if (! NUMdefined (time)) return Melder_error ("Selection is voiceless."); my startSelection = my endSelection = time; FunctionEditor_marksChanged (me); } END void FunctionEditor_SoundAnalysis_draw (I) { iam (FunctionEditor); double pitchFloor_hidden = ClassFunction_convertStandardToSpecialUnit (classPitch, my pitch.floor, Pitch_LEVEL_FREQUENCY, my pitch.unit); double pitchCeiling_hidden = ClassFunction_convertStandardToSpecialUnit (classPitch, my pitch.ceiling, Pitch_LEVEL_FREQUENCY, my pitch.unit); double pitchFloor_overt = ClassFunction_convertToNonlogarithmic (classPitch, pitchFloor_hidden, Pitch_LEVEL_FREQUENCY, my pitch.unit); double pitchCeiling_overt = ClassFunction_convertToNonlogarithmic (classPitch, pitchCeiling_hidden, Pitch_LEVEL_FREQUENCY, my pitch.unit); double pitchViewFrom_overt = my pitch.viewFrom < my pitch.viewTo ? my pitch.viewFrom : pitchFloor_overt; double pitchViewTo_overt = my pitch.viewFrom < my pitch.viewTo ? my pitch.viewTo : pitchCeiling_overt; double pitchViewFrom_hidden = ClassFunction_isUnitLogarithmic (classPitch, Pitch_LEVEL_FREQUENCY, my pitch.unit) ? log10 (pitchViewFrom_overt) : pitchViewFrom_overt; double pitchViewTo_hidden = ClassFunction_isUnitLogarithmic (classPitch, Pitch_LEVEL_FREQUENCY, my pitch.unit) ? log10 (pitchViewTo_overt) : pitchViewTo_overt; Graphics_setWindow (my graphics, 0.0, 1.0, 0.0, 1.0); Graphics_setColour (my graphics, Graphics_WHITE); Graphics_fillRectangle (my graphics, 0.0, 1.0, 0.0, 1.0); Graphics_setColour (my graphics, Graphics_BLACK); Graphics_rectangle (my graphics, 0.0, 1.0, 0.0, 1.0); if (my endWindow - my startWindow > my longestAnalysis) { Graphics_setFont (my graphics, Graphics_HELVETICA); Graphics_setFontSize (my graphics, 9); Graphics_setTextAlignment (my graphics, Graphics_CENTRE, Graphics_HALF); Graphics_printf (my graphics, 0.5, 0.67, "To see the analyses, zoom in to at most %s seconds,", Melder_half (my longestAnalysis)); Graphics_text (my graphics, 0.5, 0.33, "or raise the \"longest analysis\" setting in the View menu."); Graphics_setFontSize (my graphics, 10); return; } computeSpectrogram (me); if (my spectrogram.show && my spectrogram.data != NULL) { Spectrogram_paintInside (my spectrogram.data, my graphics, my startWindow, my endWindow, my spectrogram.viewFrom, my spectrogram.viewTo, my spectrogram.maximum, my spectrogram.autoscaling, my spectrogram.dynamicRange, my spectrogram.preemphasis, my spectrogram.dynamicCompression); } computePitch (me); if (my pitch.show && my pitch.data != NULL) { double periodsPerAnalysisWindow = my pitch.method == 1 ? 3.0 : 1.0; double greatestNonUndersamplingTimeStep = 0.5 * periodsPerAnalysisWindow / my pitch.floor; double defaultTimeStep = 0.5 * greatestNonUndersamplingTimeStep; double timeStep = my timeStepStrategy == 2 ? my fixedTimeStep : my timeStepStrategy == 3 ? (my endWindow - my startWindow) / my numberOfTimeStepsPerView : defaultTimeStep; int undersampled = timeStep > greatestNonUndersamplingTimeStep; long numberOfVisiblePitchPoints = (long) ((my endWindow - my startWindow) / timeStep); Graphics_setColour (my graphics, Graphics_CYAN); Graphics_setLineWidth (my graphics, 3.0); if (undersampled || numberOfVisiblePitchPoints < 101) { Pitch_drawInside (my pitch.data, my graphics, my startWindow, my endWindow, pitchViewFrom_overt, pitchViewTo_overt, 2, my pitch.unit); } if (! undersampled) { Pitch_drawInside (my pitch.data, my graphics, my startWindow, my endWindow, pitchViewFrom_overt, pitchViewTo_overt, FALSE, my pitch.unit); } Graphics_setColour (my graphics, Graphics_BLUE); Graphics_setLineWidth (my graphics, 1.0); if (undersampled || numberOfVisiblePitchPoints < 101) { Pitch_drawInside (my pitch.data, my graphics, my startWindow, my endWindow, pitchViewFrom_overt, pitchViewTo_overt, 1, my pitch.unit); } if (! undersampled) { Pitch_drawInside (my pitch.data, my graphics, my startWindow, my endWindow, pitchViewFrom_overt, pitchViewTo_overt, FALSE, my pitch.unit); } Graphics_setColour (my graphics, Graphics_BLACK); } computeIntensity (me); if (my intensity.show && my intensity.data != NULL) { Graphics_setColour (my graphics, my spectrogram.show ? Graphics_YELLOW : Graphics_LIME); Graphics_setLineWidth (my graphics, my spectrogram.show ? 1.0 : 3.0); Intensity_drawInside (my intensity.data, my graphics, my startWindow, my endWindow, my intensity.viewFrom, my intensity.viewTo); Graphics_setLineWidth (my graphics, 1.0); Graphics_setColour (my graphics, Graphics_BLACK); } computeFormants (me); if (my formant.show && my formant.data != NULL) { Graphics_setColour (my graphics, Graphics_RED); Formant_drawSpeckles_inside (my formant.data, my graphics, my startWindow, my endWindow, my spectrogram.viewFrom, my spectrogram.viewTo, my formant.dynamicRange, my formant.dotSize); Graphics_setColour (my graphics, Graphics_BLACK); } /* * Draw vertical scales. */ if (my pitch.show) { double pitchCursor_overt = NUMundefined, pitchCursor_hidden = NUMundefined; Graphics_setWindow (my graphics, my startWindow, my endWindow, pitchViewFrom_hidden, pitchViewTo_hidden); Graphics_setColour (my graphics, Graphics_NAVY); if (my pitch.data) { if (my startSelection == my endSelection) pitchCursor_hidden = Pitch_getValueAtTime (my pitch.data, my startSelection, my pitch.unit, 1); else pitchCursor_hidden = Pitch_getMean (my pitch.data, my startSelection, my endSelection, my pitch.unit); pitchCursor_overt = ClassFunction_convertToNonlogarithmic (classPitch, pitchCursor_hidden, Pitch_LEVEL_FREQUENCY, my pitch.unit); if (NUMdefined (pitchCursor_hidden)) { Graphics_setTextAlignment (my graphics, Graphics_LEFT, Graphics_HALF); Graphics_printf (my graphics, my endWindow, pitchCursor_hidden, "%.5g %s", pitchCursor_overt, ClassFunction_getUnitText (classPitch, Pitch_LEVEL_FREQUENCY, my pitch.unit, Function_UNIT_TEXT_SHORT | Function_UNIT_TEXT_GRAPHICAL)); } if (! NUMdefined (pitchCursor_hidden) || Graphics_dyWCtoMM (my graphics, pitchCursor_hidden - pitchViewFrom_hidden) > 5.0) { Graphics_setTextAlignment (my graphics, Graphics_LEFT, Graphics_BOTTOM); Graphics_printf (my graphics, my endWindow, pitchViewFrom_hidden - Graphics_dyMMtoWC (my graphics, 0.5), "%.4g %s", pitchViewFrom_overt, ClassFunction_getUnitText (classPitch, Pitch_LEVEL_FREQUENCY, my pitch.unit, Function_UNIT_TEXT_SHORT | Function_UNIT_TEXT_GRAPHICAL)); } if (! NUMdefined (pitchCursor_hidden) || Graphics_dyWCtoMM (my graphics, pitchViewTo_hidden - pitchCursor_hidden) > 5.0) { Graphics_setTextAlignment (my graphics, Graphics_LEFT, Graphics_TOP); Graphics_printf (my graphics, my endWindow, pitchViewTo_hidden, "%.4g %s", pitchViewTo_overt, ClassFunction_getUnitText (classPitch, Pitch_LEVEL_FREQUENCY, my pitch.unit, Function_UNIT_TEXT_SHORT | Function_UNIT_TEXT_GRAPHICAL)); } } else { Graphics_setTextAlignment (my graphics, Graphics_CENTRE, Graphics_HALF); Graphics_text (my graphics, 0.5 * (my startWindow + my endWindow), 0.5 * (pitchViewFrom_hidden + pitchViewTo_hidden), "(Cannot show pitch contour. Zoom out or change bottom of pitch range in pitch settings.)"); } Graphics_setColour (my graphics, Graphics_BLACK); } if (my intensity.show) { double intensityCursor = NUMundefined; int intensityCursorVisible; int textColour, alignment; double y; if (! my pitch.show) textColour = Graphics_GREEN, alignment = Graphics_LEFT, y = my endWindow; else if (! my spectrogram.show && ! my formant.show) textColour = Graphics_GREEN, alignment = Graphics_RIGHT, y = my startWindow; else textColour = my spectrogram.show ? Graphics_LIME : Graphics_GREEN, alignment = Graphics_RIGHT, y = my endWindow; Graphics_setWindow (my graphics, my startWindow, my endWindow, my intensity.viewFrom, my intensity.viewTo); if (my intensity.data) { if (my startSelection == my endSelection) { intensityCursor = Vector_getValueAtX (my intensity.data, my startSelection, 1); } else { intensityCursor = Intensity_getAverage (my intensity.data, my startSelection, my endSelection, my intensity.averagingMethod); } } Graphics_setColour (my graphics, textColour); intensityCursorVisible = NUMdefined (intensityCursor) && intensityCursor > my intensity.viewFrom && intensityCursor < my intensity.viewTo; if (intensityCursorVisible) { static const char *methodString [] = { " (.5)", " (\\muE)", " (\\muS)", " (\\mu)" }; Graphics_setTextAlignment (my graphics, alignment, Graphics_HALF); Graphics_printf (my graphics, y, intensityCursor, "%.4g dB%s", intensityCursor, my startSelection == my endSelection ? "" : methodString [my intensity.averagingMethod]); } if (! intensityCursorVisible || Graphics_dyWCtoMM (my graphics, intensityCursor - my intensity.viewFrom) > 5.0) { Graphics_setTextAlignment (my graphics, alignment, Graphics_BOTTOM); Graphics_printf (my graphics, y, my intensity.viewFrom - Graphics_dyMMtoWC (my graphics, 0.5), "%.4g dB", my intensity.viewFrom); } if (! intensityCursorVisible || Graphics_dyWCtoMM (my graphics, my intensity.viewTo - intensityCursor) > 5.0) { Graphics_setTextAlignment (my graphics, alignment, Graphics_TOP); Graphics_printf (my graphics, y, my intensity.viewTo, "%.4g dB", my intensity.viewTo); } Graphics_setColour (my graphics, Graphics_BLACK); } if (my spectrogram.show || my formant.show) { char text [100]; double textWidth; int frequencyCursorVisible = my spectrogram.cursor > my spectrogram.viewFrom && my spectrogram.cursor < my spectrogram.viewTo; Graphics_setWindow (my graphics, my startWindow, my endWindow, my spectrogram.viewFrom, my spectrogram.viewTo); /* * Range marks. */ Graphics_setLineType (my graphics, Graphics_DRAWN); Graphics_setColour (my graphics, Graphics_BLACK); if (! frequencyCursorVisible || Graphics_dyWCtoMM (my graphics, my spectrogram.cursor - my spectrogram.viewFrom) > 5.0) { sprintf (text, "%.4g Hz", my spectrogram.viewFrom); textWidth = Graphics_textWidth (my graphics, text) + Graphics_dxMMtoWC (my graphics, 0.5); /*Graphics_line (my graphics, my startWindow - textWidth, my spectrogram.viewFrom, my startWindow, my spectrogram.viewFrom);*/ Graphics_setTextAlignment (my graphics, Graphics_RIGHT, Graphics_BOTTOM); Graphics_text (my graphics, my startWindow, my spectrogram.viewFrom - Graphics_dyMMtoWC (my graphics, 0.5), text); } if (! frequencyCursorVisible || Graphics_dyWCtoMM (my graphics, my spectrogram.viewTo - my spectrogram.cursor) > 5.0) { sprintf (text, "%.4g Hz", my spectrogram.viewTo); textWidth = Graphics_textWidth (my graphics, text) + Graphics_dxMMtoWC (my graphics, 0.5); /*Graphics_line (my graphics, my startWindow - textWidth, my spectrogram.viewTo, my startWindow, my spectrogram.viewTo);*/ Graphics_setTextAlignment (my graphics, Graphics_RIGHT, Graphics_TOP); Graphics_text (my graphics, my startWindow, my spectrogram.viewTo, text); } /* * Cursor lines. */ Graphics_setLineType (my graphics, Graphics_DOTTED); Graphics_setColour (my graphics, Graphics_RED); if (frequencyCursorVisible) { double x = my startWindow, y = my spectrogram.cursor; Graphics_setTextAlignment (my graphics, Graphics_RIGHT, Graphics_HALF); Graphics_printf (my graphics, x, y, "%.5g Hz", y); Graphics_line (my graphics, x, y, my endWindow, y); } /* if (my startSelection >= my startWindow && my startSelection <= my endWindow) Graphics_line (my graphics, my startSelection, my spectrogram.viewFrom, my startSelection, my spectrogram.viewTo); if (my endSelection > my startWindow && my endSelection < my endWindow && my endSelection != my startSelection) Graphics_line (my graphics, my endSelection, my spectrogram.viewFrom, my endSelection, my spectrogram.viewTo);*/ /* * Cadre. */ Graphics_setLineType (my graphics, Graphics_DRAWN); Graphics_setColour (my graphics, Graphics_BLACK); Graphics_rectangle (my graphics, my startWindow, my endWindow, my spectrogram.viewFrom, my spectrogram.viewTo); } } void FunctionEditor_SoundAnalysis_drawPulses (I) { iam (FunctionEditor); computePulses (me); if (my pulses.show && my endWindow - my startWindow <= my longestAnalysis && my pulses.data != NULL) { long i; PointProcess point = my pulses.data; Graphics_setWindow (my graphics, my startWindow, my endWindow, -1.0, 1.0); Graphics_setColour (my graphics, Graphics_BLUE); if (point -> nt < 2000) for (i = 1; i <= point -> nt; i ++) { double t = point -> t [i]; if (t >= my startWindow && t <= my endWindow) Graphics_line (my graphics, t, -0.9, t, 0.9); } Graphics_setColour (my graphics, Graphics_BLACK); } } FORM (FunctionEditor, cb_logSettings, "Log settings", "Log files") OPTIONMENU ("Write log 1 to", 3) OPTION ("Log file only") OPTION ("Info window only") OPTION ("Log file and Info window") LABEL ("", "Log file 1:") TEXTFIELD ("Log file 1", LOG_1_FILE_NAME) LABEL ("", "Log 1 format:") TEXTFIELD ("Log 1 format", LOG_1_FORMAT) OPTIONMENU ("Write log 2 to", 3) OPTION ("Log file only") OPTION ("Info window only") OPTION ("Log file and Info window") LABEL ("", "Log file 2:") TEXTFIELD ("Log file 2", LOG_2_FILE_NAME) LABEL ("", "Log 2 format:") TEXTFIELD ("Log 2 format", LOG_2_FORMAT) LABEL ("", "Log script 3:") TEXTFIELD ("Log script 3", LOG_3_FILE_NAME) LABEL ("", "Log script 4:") TEXTFIELD ("Log script 4", LOG_4_FILE_NAME) OK SET_INTEGER ("Write log 1 to", preferences.log[0].toLogFile + 2 * preferences.log[0].toInfoWindow) SET_STRING ("Log file 1", preferences.log[0].fileName) SET_STRING ("Log 1 format", preferences.log[0].format) SET_INTEGER ("Write log 2 to", preferences.log[1].toLogFile + 2 * preferences.log[1].toInfoWindow) SET_STRING ("Log file 2", preferences.log[1].fileName) SET_STRING ("Log 2 format", preferences.log[1].format) SET_STRING ("Log script 3", preferences.logScript3) SET_STRING ("Log script 4", preferences.logScript4) DO preferences.log[0].toLogFile = (GET_INTEGER ("Write log 1 to") & 1) != 0; preferences.log[0].toInfoWindow = (GET_INTEGER ("Write log 1 to") & 2) != 0; strcpy (preferences.log[0].fileName, GET_STRING ("Log file 1")); strcpy (preferences.log[0].format, GET_STRING ("Log 1 format")); preferences.log[1].toLogFile = (GET_INTEGER ("Write log 2 to") & 1) != 0; preferences.log[1].toInfoWindow = (GET_INTEGER ("Write log 2 to") & 2) != 0; strcpy (preferences.log[1].fileName, GET_STRING ("Log file 2")); strcpy (preferences.log[1].format, GET_STRING ("Log 2 format")); strcpy (preferences.logScript3, GET_STRING ("Log script 3")); strcpy (preferences.logScript4, GET_STRING ("Log script 4")); END static int cb_deleteLogFile (FunctionEditor me, int which) { structMelderFile file; (void) me; if (! Melder_pathToFile (preferences.log[which].fileName, & file)) return 0; MelderFile_delete (& file); return 1; } DIRECT (FunctionEditor, cb_deleteLogFile1) if (! cb_deleteLogFile (me, 0)) return 0; END DIRECT (FunctionEditor, cb_deleteLogFile2) if (! cb_deleteLogFile (me, 0)) return 0; END static int cb_log (FunctionEditor me, int which) { char format [1000], *p; double tmin, tmax; int part = makeQueriable (me, TRUE, & tmin, & tmax); iferror return 0; strcpy (format, preferences.log[which].format); for (p = format; *p !='\0'; p ++) if (*p == '\'') { /* * Found a left quote. Search for a matching right quote. */ char *q = p + 1, varName [300], *r, *s, *colon; int precision = -1; double value = NUMundefined; char *stringValue = NULL; while (*q != '\0' && *q != '\'') q ++; if (*q == '\0') break; /* No matching right quote: done with this line. */ if (q - p == 1) continue; /* Ignore empty variable names. */ /* * Found a right quote. Get potential variable name. */ for (r = p + 1, s = varName; q - r > 0; r ++, s ++) *s = *r; *s = '\0'; /* Trailing null byte. */ colon = strchr (varName, ':'); if (colon) { precision = atoi (colon + 1); *colon = '\0'; } if (strequ (varName, "time")) { value = 0.5 * (tmin + tmax); } else if (strequ (varName, "t1")) { value = tmin; } else if (strequ (varName, "t2")) { value = tmax; } else if (strequ (varName, "dur")) { value = tmax - tmin; } else if (strequ (varName, "freq")) { value = my spectrogram.cursor; } else if (strequ (varName, "tab$")) { stringValue = "\t"; } else if (strequ (varName, "editor$")) { stringValue = my name; } else if (strequ (varName, "f0")) { if (! my pitch.show) return Melder_error ("No pitch contour is visible.\nFirst choose \"Show pitch\" from the View menu."); if (! my pitch.data) { return Melder_error ("No pitch contour available (out of memory?)."); } if (part == FunctionEditor_PART_CURSOR) { value = Pitch_getValueAtTime (my pitch.data, tmin, my pitch.unit, 1); } else { value = Pitch_getMean (my pitch.data, tmin, tmax, my pitch.unit); } } else if (varName [0] == 'f' && varName [1] >= '1' && varName [1] <= '5' && varName [2] == '\0') { if (! my formant.show) return Melder_error ("No formant contour is visible.\nFirst choose \"Show formant\" from the View menu."); if (! my formant.data) { return Melder_error ("No formant contour available (out of memory?)."); } if (part == FunctionEditor_PART_CURSOR) { value = Formant_getValueAtTime (my formant.data, varName [1] - '0', tmin, 0); } else { value = Formant_getMean (my formant.data, varName [1] - '0', tmin, tmax, 0); } } else if (varName [0] == 'b' && varName [1] >= '1' && varName [1] <= '5' && varName [2] == '\0') { if (! my formant.show) return Melder_error ("No formant contour is visible.\nFirst choose \"Show formant\" from the View menu."); if (! my formant.data) { return Melder_error ("No formant contour available (out of memory?)."); } value = Formant_getBandwidthAtTime (my formant.data, varName [1] - '0', 0.5 * (tmin + tmax), 0); } else if (strequ (varName, "intensity")) { if (! my intensity.show) return Melder_error ("No intensity contour is visible.\nFirst choose \"Show intensity\" from the View menu."); if (! my intensity.data) { return Melder_error ("No intensity contour available (out of memory?)."); } if (part == FunctionEditor_PART_CURSOR) { value = Vector_getValueAtX (my intensity.data, tmin, 1); } else { value = Intensity_getAverage (my intensity.data, tmin, tmax, my intensity.averagingMethod); } } if (NUMdefined (value)) { int varlen = (q - p) - 1, headlen = p - format; char formattedNumber [400]; int arglen; if (precision >= 0) { sprintf (formattedNumber, "%.*f", precision, value); } else { sprintf (formattedNumber, "%.17g", value); } arglen = strlen (formattedNumber); strncpy (Melder_buffer1, format, headlen); strcpy (Melder_buffer1 + headlen, formattedNumber); strcpy (Melder_buffer1 + headlen + arglen, p + varlen + 2); strcpy (format, Melder_buffer1); p += arglen - 1; } else if (stringValue != NULL) { int varlen = (q - p) - 1, headlen = p - format, arglen = strlen (stringValue); strncpy (Melder_buffer1, format, headlen); strcpy (Melder_buffer1 + headlen, stringValue); strcpy (Melder_buffer1 + headlen + arglen, p + varlen + 2); strcpy (format, Melder_buffer1); p += arglen - 1; } else { p = q - 1; /* Go to before next quote. */ } } if (preferences.log[which].toInfoWindow) { Melder_info ("%s", format); } if (preferences.log[which].toLogFile) { structMelderFile file; strcat (format, "\n"); if (! Melder_relativePathToFile (preferences.log[which].fileName, & file)) return 0; if (! MelderFile_appendText (& file, format)) return 0; } return 1; } DIRECT (FunctionEditor, cb_log1) if (! cb_log (me, 0)) return 0; END DIRECT (FunctionEditor, cb_log2) if (! cb_log (me, 1)) return 0; END DIRECT (FunctionEditor, cb_logScript3) if (! DO_RunTheScriptFromAnyAddedEditorCommand (me, preferences.logScript3)) return 0; END DIRECT (FunctionEditor, cb_logScript4) if (! DO_RunTheScriptFromAnyAddedEditorCommand (me, preferences.logScript4)) return 0; END void FunctionEditor_SoundAnalysis_viewMenus (I) { iam (FunctionEditor); Editor_addCommand (me, "View", "Analysis window:", motif_INSENSITIVE, cb_showAnalyses); Editor_addCommand (me, "View", "Show analyses...", 0, cb_showAnalyses); Editor_addCommand (me, "View", "Time step settings...", 0, cb_timeStepSettings); Editor_addCommand (me, "View", "-- sound analysis --", 0, 0); } void FunctionEditor_SoundAnalysis_selectionQueries (I) { iam (FunctionEditor); Editor_addCommand (me, "Query", "-- query log --", 0, NULL); Editor_addCommand (me, "Query", "Log settings...", 0, cb_logSettings); Editor_addCommand (me, "Query", "Delete log file 1", 0, cb_deleteLogFile1); Editor_addCommand (me, "Query", "Delete log file 2", 0, cb_deleteLogFile2); Editor_addCommand (me, "Query", "Log 1", motif_F12, cb_log1); Editor_addCommand (me, "Query", "Log 2", motif_F12 + motif_SHIFT, cb_log2); Editor_addCommand (me, "Query", "Log script 3 (...)", motif_F12 + motif_OPTION, cb_logScript3); Editor_addCommand (me, "Query", "Log script 4 (...)", motif_F12 + motif_COMMAND, cb_logScript4); } void FunctionEditor_SoundAnalysis_addMenus (I) { iam (FunctionEditor); EditorMenu menu; menu = Editor_addMenu (me, "Spectrum", 0); my spectrogramToggle = EditorMenu_addCommand (menu, "Show spectrogram", motif_CHECKABLE | (preferences.spectrogram.show ? motif_CHECKED : 0), cb_showSpectrogram); EditorMenu_addCommand (menu, "Spectrogram settings...", 0, cb_spectrogramSettings); EditorMenu_addCommand (menu, "Advanced spectrogram settings...", 0, cb_advancedSpectrogramSettings); EditorMenu_addCommand (menu, "Extract visible spectrogram", 0, cb_extractVisibleSpectrogram); EditorMenu_addCommand (menu, "View spectral slice", 'L', cb_viewSpectralSlice); EditorMenu_addCommand (menu, "-- query spectrogram --", 0, NULL); EditorMenu_addCommand (menu, "Query:", motif_INSENSITIVE, cb_getFrequency /* dummy */); EditorMenu_addCommand (menu, "Get frequency at frequency cursor", 0, cb_getFrequency); EditorMenu_addCommand (menu, "Get spectral power at cursor cross", motif_F7, cb_getSpectralPowerAtCursorCross); menu = Editor_addMenu (me, "Pitch", 0); my pitchToggle = EditorMenu_addCommand (menu, "Show pitch", motif_CHECKABLE | (preferences.pitch.show ? motif_CHECKED : 0), cb_showPitch); EditorMenu_addCommand (menu, "Pitch settings...", 0, cb_pitchSettings); EditorMenu_addCommand (menu, "Advanced pitch settings...", 0, cb_advancedPitchSettings); EditorMenu_addCommand (menu, "Extract visible pitch contour", 0, cb_extractVisiblePitchContour); EditorMenu_addCommand (menu, "-- query pitch --", 0, NULL); EditorMenu_addCommand (menu, "Query:", motif_INSENSITIVE, cb_getFrequency /* dummy */); EditorMenu_addCommand (menu, "Pitch listing", 0, cb_pitchListing); EditorMenu_addCommand (menu, "Get pitch", motif_F5, cb_getPitch); EditorMenu_addCommand (menu, "Get minimum pitch", motif_F5 + motif_COMMAND, cb_getMinimumPitch); EditorMenu_addCommand (menu, "Get maximum pitch", motif_F5 + motif_SHIFT, cb_getMaximumPitch); EditorMenu_addCommand (menu, "-- select pitch --", 0, NULL); EditorMenu_addCommand (menu, "Select:", motif_INSENSITIVE, cb_moveCursorToMinimumPitch /* dummy */); EditorMenu_addCommand (menu, "Move cursor to minimum pitch", motif_COMMAND + motif_SHIFT + 'L', cb_moveCursorToMinimumPitch); EditorMenu_addCommand (menu, "Move cursor to maximum pitch", motif_COMMAND + motif_SHIFT + 'H', cb_moveCursorToMaximumPitch); menu = Editor_addMenu (me, "Intensity", 0); my intensityToggle = EditorMenu_addCommand (menu, "Show intensity", motif_CHECKABLE | (preferences.intensity.show ? motif_CHECKED : 0), cb_showIntensity); EditorMenu_addCommand (menu, "Intensity settings...", 0, cb_intensitySettings); EditorMenu_addCommand (menu, "Extract visible intensity contour", 0, cb_extractVisibleIntensityContour); EditorMenu_addCommand (menu, "-- query intensity --", 0, NULL); EditorMenu_addCommand (menu, "Query:", motif_INSENSITIVE, cb_getFrequency /* dummy */); EditorMenu_addCommand (menu, "Intensity listing", 0, cb_intensityListing); EditorMenu_addCommand (menu, "Get intensity", motif_F8, cb_getIntensity); menu = Editor_addMenu (me, "Formant", 0); my formantToggle = EditorMenu_addCommand (menu, "Show formants", motif_CHECKABLE | (preferences.formant.show ? motif_CHECKED : 0), cb_showFormants); EditorMenu_addCommand (menu, "Formant settings...", 0, cb_formantSettings); EditorMenu_addCommand (menu, "Advanced formant settings...", 0, cb_advancedFormantSettings); EditorMenu_addCommand (menu, "Extract visible formant contour", 0, cb_extractVisibleFormantContour); EditorMenu_addCommand (menu, "-- query formants --", 0, NULL); EditorMenu_addCommand (menu, "Query:", motif_INSENSITIVE, cb_getFrequency /* dummy */); EditorMenu_addCommand (menu, "Formant listing", 0, cb_formantListing); EditorMenu_addCommand (menu, "Get first formant", motif_F1, cb_getFirstFormant); EditorMenu_addCommand (menu, "Get first bandwidth", 0, cb_getFirstBandwidth); EditorMenu_addCommand (menu, "Get second formant", motif_F2, cb_getSecondFormant); EditorMenu_addCommand (menu, "Get second bandwidth", 0, cb_getSecondBandwidth); EditorMenu_addCommand (menu, "Get third formant", motif_F3, cb_getThirdFormant); EditorMenu_addCommand (menu, "Get third bandwidth", 0, cb_getThirdBandwidth); EditorMenu_addCommand (menu, "Get fourth formant", motif_F4, cb_getFourthFormant); EditorMenu_addCommand (menu, "Get fourth bandwidth", 0, cb_getFourthBandwidth); EditorMenu_addCommand (menu, "Get formant...", 0, cb_getFormant); EditorMenu_addCommand (menu, "Get bandwidth...", 0, cb_getBandwidth); menu = Editor_addMenu (me, "Pulses", 0); my pulsesToggle = EditorMenu_addCommand (menu, "Show pulses", motif_CHECKABLE | (preferences.pulses.show ? motif_CHECKED : 0), cb_showPulses); EditorMenu_addCommand (menu, "Extract visible pulses", 0, cb_extractVisiblePulses); EditorMenu_addCommand (menu, "Advanced pulses settings...", 0, cb_advancedPulsesSettings); EditorMenu_addCommand (menu, "-- query pulses --", 0, NULL); EditorMenu_addCommand (menu, "Query:", motif_INSENSITIVE, cb_getFrequency /* dummy */); EditorMenu_addCommand (menu, "Voice report", 0, cb_voiceReport); EditorMenu_addCommand (menu, "Pulse listing", 0, cb_pulseListing); /* EditorMenu_addCommand (menu, "Get jitter (local)", 0, cb_getJitter_local); EditorMenu_addCommand (menu, "Get jitter (local, absolute)", 0, cb_getJitter_local_absolute); EditorMenu_addCommand (menu, "Get jitter (rap)", 0, cb_getJitter_rap); EditorMenu_addCommand (menu, "Get jitter (ppq5)", 0, cb_getJitter_ppq5); EditorMenu_addCommand (menu, "Get jitter (ddp)", 0, cb_getJitter_ddp); EditorMenu_addCommand (menu, "Get shimmer (local)", 0, cb_getShimmer_local); EditorMenu_addCommand (menu, "Get shimmer (local_dB)", 0, cb_getShimmer_local_dB); EditorMenu_addCommand (menu, "Get shimmer (apq3)", 0, cb_getShimmer_apq3); EditorMenu_addCommand (menu, "Get shimmer (apq5)", 0, cb_getShimmer_apq5); EditorMenu_addCommand (menu, "Get shimmer (apq11)", 0, cb_getShimmer_apq11); EditorMenu_addCommand (menu, "Get shimmer (dda)", 0, cb_getShimmer_dda); */ } void FunctionEditor_SoundAnalysis_init (I) { iam (FunctionEditor); my longestAnalysis = preferences.longestAnalysis; if (preferences.timeStepStrategy < 1 || preferences.timeStepStrategy > 3) preferences.timeStepStrategy = 1; if (preferences.log[0].toLogFile == FALSE && preferences.log[0].toInfoWindow == FALSE) preferences.log[0].toLogFile = TRUE, preferences.log[0].toInfoWindow = TRUE; if (preferences.log[1].toLogFile == FALSE && preferences.log[1].toInfoWindow == FALSE) preferences.log[1].toLogFile = TRUE, preferences.log[1].toInfoWindow = TRUE; my timeStepStrategy = preferences.timeStepStrategy; my fixedTimeStep = preferences.fixedTimeStep; my numberOfTimeStepsPerView = preferences.numberOfTimeStepsPerView; my spectrogram = preferences.spectrogram; my pitch = preferences.pitch; my intensity = preferences.intensity; my formant = preferences.formant; my pulses = preferences.pulses; } /* End of file FunctionEditor_SoundAnalysis.c */