{"metadata":{"kernelspec":{"language":"python","display_name":"Python 3","name":"python3"},"language_info":{"pygments_lexer":"ipython3","nbconvert_exporter":"python","version":"3.6.4","file_extension":".py","codemirror_mode":{"name":"ipython","version":3},"name":"python","mimetype":"text/x-python"}},"nbformat_minor":4,"nbformat":4,"cells":[{"cell_type":"code","source":"import numpy as np # linear algebra\nimport pandas as pd # data processing, CSV file I/O (e.g. pd.read_csv)\n\n# Input data files are available in the read-only \"../input/\" directory\n# For example, running this (by clicking run or pressing Shift+Enter) will list all files under the input directory\n\nimport os\n\n#for dirname, _, filenames in os.walk('/kaggle/input'):\n#    for filename in filenames:\n#        print(os.path.join(dirname, filename))\n\nimport h5py\nfrom datetime import datetime\nimport matplotlib.pyplot as plt\n\n# function to plot spectrogram given data name\ndef plot_spectrogram(filename, LIGO=\"H1\"):\n  # open file\n  f = h5py.File('/kaggle/input/g2net-detecting-continuous-gravitational-waves/train/'+filename, 'r')\n\n  # read fourier transform coefficients\n  LIGO_SFT = f[filename[:-5]][LIGO]['SFTs']\n  spectrogram = np.absolute(np.array(LIGO_SFT[:]))\n\n  # plot spectrogram\n  fig = plt.figure(num=None, figsize=(30, 4), dpi=80, facecolor='w', edgecolor='k')\n  fig.subplots_adjust(top=0.99, bottom=0.02, left=0.02, right=0.99)\n  plt.imshow(spectrogram, cmap='jet', aspect='equal', origin ='lower')\n\nfilename_list = ['001121a05.hdf5',\n                 '004f23b2d.hdf5',\n                 '02c478b09.hdf5',\n                 '03189bb3d.hdf5',\n                 '01bcf6533.hdf5']\n\nfor elt in filename_list:\n    plot_spectrogram(elt)\n\nplt.show()\n","metadata":{"_uuid":"8f2839f25d086af736a60e9eeb907d3b93b6e0e5","_cell_guid":"b1076dfc-b9ad-4769-8c92-a6c4dae69d19","execution":{"iopub.status.busy":"2022-10-05T16:32:29.814839Z","iopub.execute_input":"2022-10-05T16:32:29.815372Z","iopub.status.idle":"2022-10-05T16:32:35.399137Z","shell.execute_reply.started":"2022-10-05T16:32:29.815319Z","shell.execute_reply":"2022-10-05T16:32:35.397823Z"},"trusted":true},"execution_count":null,"outputs":[]},{"cell_type":"markdown","source":"Time along x (with some gaps but mostly one point every half hour), frequency along y.\n\nNext, need to find what the length scale of the chirp is.\n","metadata":{}},{"cell_type":"markdown","source":"","metadata":{}}]}