wavspa 0.1.0

checksums.yaml

@@ -0,0 +1,7 @@
+---
+SHA256:
+  metadata.gz: 3619e3dfdcaa48cc41eb51037b91011391cafa4c7e2fe9aa066790bc271d73d7
+  data.tar.gz: 320b9f33e00b3f3b367cd265f04835702585c6a5aa8d2e4a610022930adbf1ac
+SHA512:
+  metadata.gz: 5b3d0e69c4e5db5dc24d233769728e1bf30da2581ab2affbbf80eda56d495a554212aa84d70b76cffed0351572bf05ac2cb75b333eb3421ee624a01f7378adea
+  data.tar.gz: 747b1bd4ba30bab6543d5add4f4c2a6b3dbb9e6f3a0451d8e23ac1cf6db7256056b6d0a30e733295b95aa5a9b6b20dce41a10782ddbcd22ace64ca68b0fed77d

data/.gitignore

@@ -0,0 +1,8 @@
+/.bundle/
+/.yardoc
+/_yardoc/
+/coverage/
+/doc/
+/pkg/
+/spec/reports/
+/tmp/

data/Gemfile

@@ -0,0 +1,4 @@
+source "https://rubygems.org"
+
+# Specify your gem's dependencies in wavspa.gemspec
+gemspec

data/LICENSE.txt

@@ -0,0 +1,21 @@
+The MIT License (MIT)
+
+Copyright (c) 2019 Hirosho Kuwagata
+
+Permission is hereby granted, free of charge, to any person obtaining a copy
+of this software and associated documentation files (the "Software"), to deal
+in the Software without restriction, including without limitation the rights
+to use, copy, modify, merge, publish, distribute, sublicense, and/or sell
+copies of the Software, and to permit persons to whom the Software is
+furnished to do so, subject to the following conditions:
+
+The above copyright notice and this permission notice shall be included in
+all copies or substantial portions of the Software.
+
+THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR
+IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
+FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE
+AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER
+LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM,
+OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN
+THE SOFTWARE.

data/README.md

@@ -0,0 +1,196 @@
+# wavspa
+
+"wavspa" is a WAV file spectrum analyzer. Includes FFT and wavelet analyzers. 
+
+## Installation
+
+Add this line to your application's Gemfile:
+
+```ruby
+gem 'wavspa'
+```
+
+And then execute:
+
+    $ bundle
+
+Or install it yourself as:
+
+    $ gem install wavspa
+
+This application uses OpenMP, but by default tries to link the libgomp. When using another implementation (such as intel OpenMP) with macOS etc, specify the libary path and library name by command line options as follows:
+
+    $ gem insatll wavspa -- --with-openomp ~/env/opts/lib --omp-name omp
+
+## Usage
+
+### FFT analyzer
+
+```
+wavfft [options] WAV-FILE
+
+options:
+    -o, --output=FILE
+    -p, --preset=NAME
+    -f, --fft-size=SIZE
+    -u, --unit-time=CS
+    -W, --output-width=SIZE
+    -w, --window-function=FUNCTION
+    -r, --frequency-range=LO,HI,
+    -m, --scale-mode=MODE
+        --show-params
+    -c, --col-steps=SIZE
+    -F, --draw-freq-line
+    -T, --draw-time-line
+    -v, --verbose
+    -h, --help
+```
+
+#### options
+
+<dl>
+  <dt>-o, --output</dt>
+  <dd>specify the output file name. by default, output to the file that extension of input file changed to ".png".</dd>
+
+  <dt>-p, --preset=NAME</dt>
+  <dd>specify preset settings. you can specify one of "default", "32k" or "cd".</dd>
+
+  <dt>-f, --fft-size=FFT</dt>
+  <dd>specify the FFT size (by number of samples, this value must spcify a power of 2).</dd>
+
+  <dt>-u, --unit-time=CS</dt>
+  <dd>specify the unit time of short time FFT in centiseconds (1/100 sec). This value is the time per pixel in horizontal direction of the output PNG.</dt>
+
+  <dt>-W, --output-width=SIZE</dt>
+  <dd>specify the hight of output PNG (by number of pixels).</dd>
+
+  <dt>-r, --frequency-range=LO,HI</dt>
+  <dd>specify the frequency band show on the output PNG (upper limit to "HI", and lower limit to "LO").<dd>
+
+  <dt>-m, --scale=mode</dt>
+  <dd>specify the mode of grid scale in vertical direction of th output PNG. you can specify one of "LOGSCALE" or "LINEARSCALE" (default is "LOGSCALE").</dd>
+
+  <dt>-c, --col-steps=SIZE</dt>
+  <dd>specify the horizontal magnify ratio of th output file.</dd>
+
+  <dt>--show-params</dt>
+  <dd>show sumarry of settings.</dd>
+
+  <dt>-F, --draw-freq-line</dt>
+  <dd>enable frequency line (vertical grid).</dd>
+
+  <dt>-F, --draw-freq-line</dt>
+  <dd>enable time line (horizontal grid).</dd>
+
+   <dt>-v, --verbose</dt>
+   <dd>enable verbose mode</dd>
+
+   <dt>-h, --help</dt>
+   <dd>show help message</dd>
+</dl>
+
+### Wavelet analyzer
+
+```
+wavlet [options] WAV-FILE
+
+options:
+    -o, --output=FILE
+    -p, --preset=NAME
+    -s, --sigma=NUMBER
+    -a, --amplitude-mode
+    -u, --unit-time=CENTISECOND
+    -W, --output-width=SIZE
+    -r, --frequency-range=LO,HI
+        --floor-gain=DB
+        --ceil-gain=DB
+    -m, --scale-mode=STRING
+    -c, --col-steps=SIZE
+        --show-params
+    -F, --draw-freq-line
+    -T, --draw-time-line
+    -v, --verbose
+```
+
+<dl>
+  <dt>-o, --output</dt>
+  <dd>specify the output file name. by default, output to the file that extension of input file changed to ".png".</dd>
+
+  <dt>-p, --preset=NAME</dt>
+  <dd>specify preset settings. you can specify one of "default" or "cd".</dd>
+
+  <dt>-s, --sigma=FLOAT</dt>
+  <dd>specify the sigma value. Increasing this value improves the accuracy but increases the processing load.</dd>
+
+  <dt>-a, --amplitude-mode</dt>
+  <dd>When this option is specified, the amplitude spectrum is output (otherwise, the power spectrum is output).</dd>
+
+  <dt>-u, --unit-time=CS</dt>
+  <dd>specify the unit time of wavelet stransform in centiseconds (1/100 sec). This value is the time per pixel in horizontal direction of the output PNG.</dt>
+
+  <dt>-W, --output-width=SIZE</dt>
+  <dd>specify the hight of output PNG (by number of pixels).</dd>
+
+  <dt>-r, --frequency-range=LO,HI</dt>
+  <dd>specify the frequency band show on the output PNG (upper limit to "HI", and lower limit to "LO").<dd>
+
+  <dt>--floor-gain=DB</dt>
+  <dd>specify the upper limit value of the gain to be displayed. values exceeding this number are displayed as saturated.</dd>
+
+  <dt>--ceil-gain=DB</dt>
+  <dd>specify the lower limit value of the gain to be displayed. values less than this number will be masked.</dd>
+
+  <dt>-m, --scale=mode</dt>
+  <dd>specify the mode of grid scale in vertical direction of th output PNG. you can specify one of "LOGSCALE" or "LINEARSCALE" (default is "LOGSCALE").</dd>
+
+  <dt>-c, --col-steps=SIZE</dt>
+  <dd>specify the horizontal magnify ratio of th output file.</dd>
+
+  <dt>--show-params</dt>
+  <dd>show sumarry of settings.</dd>
+
+  <dt>-F, --draw-freq-line</dt>
+  <dd>enable frequency line (vertical grid).</dd>
+
+  <dt>-F, --draw-freq-line</dt>
+  <dd>enable time line (horizontal grid).</dd>
+
+   <dt>-v, --verbose</dt>
+   <dd>enable verbose mode</dd>
+
+   <dt>-h, --help</dt>
+   <dd>show help message</dd>
+</dl>
+
+## Output example
+As a sample data, transformed from "Call to Quarters" (https://archive.org/details/CallToQuarters).
+
+### FFT
+```
+ wavfft -p cd -F -T -v -o Call_To_Quarters-fft.png Call_To_Quarters.wav
+```
+
+![FFT result](example/Call_To_Quarters-fft.png)
+
+
+### Wavelet
+```
+ wavlet -p cd -F -T -c 2 -v -o Call_To_Quarters-wavelet.png Call_To_Quarters.wav
+```
+
+![wavelet result](example/Call_To_Quarters-wavelet.png)
+
+
+## Contributing
+
+Bug reports and pull requests are welcome on GitHub at https://github.com/kwgt/wavspa.
+
+## その他
+以下のソースコード及びデータを流用しています。各々の開発者の方に感謝いたします。
+
+* wavfftで使用しているFFT処理は大浦氏が配布されているコード(通称「大浦FFT」 http://www.kurims.kyoto-u.ac.jp/~ooura/fft.html) を流用しています(ext/wavspa/fft/fftsg.c)。
+* プロット時に使用しているフォントはM+ BITMAP (http://mplus-fonts.osdn.jp/mplus-bitmap-fonts/) で配布されているm+ gothic 10rを変換し使用しています(ext/wavspa/fb/rb_fb.c)。
+
+## License
+
+The gem is available as open source under the terms of the [MIT License](https://opensource.org/licenses/MIT).

data/Rakefile

@@ -0,0 +1,2 @@
+require "bundler/gem_tasks"
+task :default => :spec

data/bin/wavfft

@@ -0,0 +1,167 @@
+#! /usr/bin/env ruby
+# coding: utf-8
+
+#
+# Spectrum analyzer for WAV file
+#
+#   Copyright (C) 2019 Hiroshi Kuwagata <kgt9221@gmail.com>
+#
+
+require 'optparse'
+require 'pathname'
+require 'pp'
+
+#
+# パッケージ内で参照する定数の定義
+#
+
+TRADITIONAL_NAME = "Spectrum analyzer for WAV file"
+APP_NAME         = "wavfft"
+
+BASE_DIR         = Pathname.new(File.expand_path(__FILE__)).dirname.parent
+LIB_DIR          = BASE_DIR + "lib"
+PKG_LIB_DIR      = LIB_DIR + "wavspa"
+APP_LIB_DIR      = PKG_LIB_DIR + APP_NAME
+
+#
+# ライブラリロードパスの差し込み
+#
+
+$LOAD_PATH.unshift(LIB_DIR)
+
+#
+# アプリケーション本体の読み込み
+#
+
+require "#{PKG_LIB_DIR + "version"}"
+require "#{PKG_LIB_DIR + "common"}"
+require "#{APP_LIB_DIR + "preset"}"
+require "#{APP_LIB_DIR + "main"}"
+
+#
+# デフォルトパラメータのセット
+#
+
+include WavSpectrumAnalyzer
+
+params = FFTApp::PRESET_TABLE["default"]
+output = nil
+
+#
+# コマンドラインオプションのパース
+#
+
+OptionParser.new {|opt|
+  window_funcs = %w{
+    RECTANGULAR HAMMING HANN BLACKMAN BLACKMAN_NUTTALL FLAT_TOP
+  }
+
+  opt.banner += " WAV-FILE"
+  opt.version = VERSION
+
+  opt.on("-o", "--output=FILE", String) { |name|
+    output = name
+  }
+
+  opt.on("-p", "--preset=NAME", String) { |name|
+    name.downcase!
+
+    if name == "help"
+      FFTApp::PRESET_TABLE.keys.each { |str|
+        printf("    %s\n", str)
+      }
+      exit
+    end
+
+    if FFTApp::PRESET_TABLE.include?(name)
+      params = FFTApp::PRESET_TABLE[name]
+    else
+      STDERR.print("error: unknown preset name.\n")
+      exit(1)
+    end
+  }
+
+  opt.on("-f", "--fft-size=SIZE", Integer) {|val|
+    params[:fft_size] = val
+  }
+
+  opt.on("-u", "--unit-time=MSEC", Integer) {|val|
+    params[:unit_time] = val
+  }
+
+  opt.on("-W", "--output-width=SIZE", Integer) {|val|
+    params[:output_width] = val
+  }
+
+  opt.on("-w", "--window-function=FUNCTION", String) {|name|
+    name.upcase!
+
+    if name == "HELP"
+      window_funcs.each {|func| printf("    %s\n", func)}
+      exit
+    end
+
+    if window_funcs.include?(name)
+      params[:window_function] = name
+    else
+      STDERR.print("unknown windows function.\n")
+      exit(1)
+    end
+  }
+
+  opt.on("-r", "--frequency-range=LO,HI", Array) { |val|
+    params[:range] = [val[0].to_f, val[1].to_f]
+  }
+
+  opt.on("-m", "--scale-mode=MODE", String) { |val|
+    name = name.upcase.to_sym
+    if not ["LINEAR", "LOGSCALE"].include?(name)
+      STDERR.print("error: unknown scale mode.\n")
+      exit(1)
+    end
+
+    params[:scale_mode] = name
+  }
+
+  opt.on("-c", "--col-steps=SIZE", Integer) {|val|
+    params[:col_step] = val
+  }
+
+  opt.on("--show-params") {
+    printf("FFT size        %20d entries\n", params[:fft_size])
+    printf("unit time       %20d msec\n", params[:unit_time])
+    printf("window function %20s\n", params[:window_function].to_s)
+    printf("scale mode      %20s\n", params[:scale_mode].to_s)
+    printf("frequency range %20s Hz\n",
+                  ("%.0f - %.0f" % [params[:range][0], params[:range][1]]))
+    printf("column steps    %20d pixels\n", params[:col_step])
+    exit
+  }
+
+  opt.on("-F", "--draw-freq-line") {
+    $draw_freq_line = true
+  }
+
+  opt.on("-T", "--draw-time-line") {
+    $draw_time_line = true
+  }
+
+  opt.on('-v', "--verbose") {
+    $verbose = true
+  }
+
+  opt.parse!(ARGV)
+
+  if ARGV.empty?
+    STDERR.print("error: target file not specified.\n")
+    exit(1)
+  end
+
+  output ||= File.basename(ARGV[0], ".wav") + ".png"
+}
+
+#
+# アプリケーションの起動
+#
+
+FFTApp.main(ARGV[0], params, output)

data/bin/wavlet

@@ -0,0 +1,160 @@
+#! /usr/bin/env ruby
+# coding: utf-8
+
+#
+# Spectrum analyzer for WAV file
+#
+#   Copyright (C) 2019 Hiroshi Kuwagata <kgt9221@gmail.com>
+#
+
+require 'optparse'
+require 'pathname'
+require 'pp'
+
+#
+# パッケージ内で参照する定数の定義
+#
+
+TRADITIONAL_NAME = "Spectrum analyzer for WAV file"
+APP_NAME         = "wavlet"
+
+BASE_DIR         = Pathname.new(File.expand_path(__FILE__)).dirname.parent
+LIB_DIR          = BASE_DIR + "lib"
+PKG_LIB_DIR      = LIB_DIR + "wavspa"
+APP_LIB_DIR      = PKG_LIB_DIR + APP_NAME
+
+#
+# ライブラリロードパスの差し込み
+#
+
+$LOAD_PATH.unshift(LIB_DIR)
+
+#
+# アプリケーション本体の読み込み
+#
+
+require "#{PKG_LIB_DIR + "version"}"
+require "#{PKG_LIB_DIR + "common"}"
+require "#{APP_LIB_DIR + "preset"}"
+require "#{APP_LIB_DIR + "main"}"
+
+#
+# デフォルトパラメータのセット
+#
+
+include WavSpectrumAnalyzer
+
+params = WaveLetApp::PRESET_TABLE["default"]
+output = nil
+
+#
+# コマンドラインオプションのパース
+#
+
+OptionParser.new {|opt|
+  opt.banner += " WAV-FILE"
+  opt.version = VERSION
+
+  opt.on("-o", "--output=FILE", String) { |name|
+    output = name
+  }
+
+  opt.on("-p", "--preset=NAME", String) { |name|
+    name.downcase!
+
+    if name == "help"
+      FFTApp::PRESET_TABLE.keys.each { |str|
+        printf("    %s\n", str)
+      }
+      exit
+    end
+
+    if WaveLetApp::PRESET_TABLE.include?(name)
+      params = WaveLetApp::PRESET_TABLE[name]
+    else
+      STDERR.print("error: unknown preset name.\n")
+      exit(1)
+    end
+  }
+
+  opt.on("-s", "--sigma=NUMBER", Float) {|val|
+    params[:sigma] = val
+  }
+
+  opt.on("-a", "--amplitude-mode", Integer) {|val|
+    params[:transform_mode] = :AMPLITUDE
+  }
+
+  opt.on("-u", "--unit-time=CENTISECOND", Integer) {|val|
+    params[:unit_time] = val * 10
+  }
+
+  opt.on("-W", "--output-width=SIZE", Integer) {|val|
+    params[:output_width] = val
+  }
+
+  opt.on("-r", "--frequency-range=LO,HI", Array) { |val|
+    params[:range] = [val[0].to_f, val[1].to_f]
+  }
+
+  opt.on("--floor-gain=DB", Float) {|val|
+    params[:floor] = val
+  }
+
+  opt.on("--ceil-gain=DB", Float) {|val|
+    params[:ceil] = val
+  }
+
+  opt.on("-m", "--scale-mode=STRING", String) {|name|
+    name = name.upcase.to_sym
+    if not ["LINEAR", "LOGSCALE"].include?(name)
+      STDERR.print("error: unknown scale mode.\n")
+      exit(1)
+    end
+
+    params[:scale_mode] = name
+  }
+
+  opt.on("-c", "--col-steps=SIZE", Integer) {|val|
+    params[:col_step] = val
+  }
+
+  opt.on("--show-params") {
+    printf("sigma           %20f\n", params[:sigma])
+    printf("unit time       %20d msec\n", params[:unit_time])
+    printf("scale mode      %20s\n", params[:scale_mode].to_s)
+    printf("ceiling         %20d\n", params[:ceiling].to_s)
+    printf("floor           %20d\n", params[:floor].to_s)
+    printf("frequency range %20s Hz\n",
+                  ("%.0f - %.0f" % [params[:range][0], params[:range][1]]))
+    printf("column steps    %20d pixels\n", params[:col_step])
+    exit
+  }
+
+  opt.on("-F", "--draw-freq-line") {
+    $draw_freq_line = true
+  }
+
+  opt.on("-T", "--draw-time-line") {
+    $draw_time_line = true
+  }
+
+  opt.on('-v', "--verbose") {
+    $verbose = true
+  }
+
+  opt.parse!(ARGV)
+
+  if ARGV.empty?
+    STDERR.print("error: target file not specified.\n")
+    exit(1)
+  end
+
+  output ||= File.basename(ARGV[0], ".wav") + ".png"
+}
+
+#
+# アプリケーションの起動
+#
+
+WaveLetApp.main(ARGV[0], params, output)

data/ext/wavspa/fb/extconf.rb

@@ -0,0 +1,3 @@
+require 'mkmf'
+
+create_makefile( "wavspa/fb")