import asyncio import os import sys import logging import numpy as np from scipy.signal import blackmanharris from quamash import QEventLoop, QtWidgets, QtCore import pyqtgraph as pg from sipyco.pipe_ipc import AsyncioChildComm from sipyco import pyon class SpectrogramItem(pg.ImageItem): def __init__(self, freq_sample, freq_base, block_size): pg.ImageItem.__init__(self) depth = 100 self.setOpts(axisOrder="row-major") self.img_array = np.zeros((depth, block_size)) self.setImage(self.img_array, autoLevels=True, autoDownsample=True) self.setRect(QtCore.QRectF((freq_base-freq_sample/2)/1e6, -float(depth), freq_sample/1e6, float(depth))) def add_block(self, block): self.img_array = np.roll(self.img_array, -1, 0) self.img_array[-1:] = np.fft.fftshift(block) self.setImage(self.img_array, autoLevels=True, autoDownsample=True) class MainWindow(pg.GraphicsLayoutWidget): def __init__(self, freq_sample, freq_base, block_size): pg.GraphicsLayoutWidget.__init__(self) self.setWindowTitle("NOPTICA Wavemeter") self.freq_sample = freq_sample self.freq_base = freq_base self.block_size = block_size self.text_ref = pg.LabelItem(size="24pt", justify="left") self.addItem(self.text_ref, row=0, col=0) self.text_locked = pg.LabelItem(size="24pt", justify="right") self.addItem(self.text_locked, row=0, col=1) self.update_ref(None, None, False) p1 = self.addPlot(row=1, col=0, colspan=2, title="REF spectrum") self.ref_spectrum = SpectrogramItem(freq_sample, freq_base, block_size) p1.addItem(self.ref_spectrum) p2 = self.addPlot(row=2, col=0, colspan=2, title="MEAS spectrum") self.meas_spectrum = SpectrogramItem(freq_sample, freq_base, block_size) p2.addItem(self.meas_spectrum) self.position = self.addPlot(row=3, col=0, colspan=2, title="Position (nm)") self.position_history = np.zeros(300) def update_ref(self, block, peak_freq, locked): if block is not None: self.ref_spectrum.add_block(block) if peak_freq is None: self.text_ref.setText("REF: NO SIGNAL") else: self.text_ref.setText("REF: {:.3f}MHz".format((self.freq_base + peak_freq)/1e6)) if locked: self.text_locked.setText("REF LASER LOCKED", color="00FF00") else: self.text_locked.setText("REF LASER UNLOCKED", color="FF0000") def update_meas(self, block): assert len(block) == self.block_size spectrum = np.abs(np.fft.fft(block*blackmanharris(self.block_size))) self.meas_spectrum.add_block(spectrum) def update_position(self, position): self.position_history = np.roll(self.position_history, -1) self.position_history[-1] = position*632.816/(4.0*np.pi) self.position.clear() self.position.plot(self.position_history) class IPCClient(AsyncioChildComm): def set_close_cb(self, close_cb): self.close_cb = close_cb async def read_pyon(self): line = await self.readline() return pyon.decode(line.decode()) async def listen(self, main_window): while True: obj = await self.read_pyon() try: action = obj["action"] del obj["action"] if action == "terminate": self.close_cb() return else: getattr(main_window, action)(**obj) except: logging.error("error processing parent message", exc_info=True) self.close_cb() def main(): freq_sample = float(sys.argv[1]) freq_base = float(sys.argv[2]) block_size = int(sys.argv[3]) app = QtWidgets.QApplication([]) loop = QEventLoop(app) asyncio.set_event_loop(loop) try: ipc = IPCClient(os.getenv("NOPTICA2_IPC")) loop.run_until_complete(ipc.connect()) try: main_window = MainWindow(freq_sample, freq_base, block_size) main_window.showFullScreen() ipc.set_close_cb(main_window.close) asyncio.ensure_future(ipc.listen(main_window)) loop.run_forever() finally: ipc.close() finally: loop.close() if __name__ == "__main__": main()