import asyncio import os import sys import logging import numpy as np 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.img_array = np.zeros((depth, block_size)) self.setImage(self.img_array, autoLevels=True) pos = np.array([0., 1., 0.5, 0.25, 0.75]) color = np.array([[0,255,255,255], [255,255,0,255], [0,0,0,255], (0, 0, 255, 255), (255, 0, 0, 255)], dtype=np.ubyte) cmap = pg.ColorMap(pos, color) lut = cmap.getLookupTable(0.0, 1.0, 256) self.setLookupTable(lut) self.setRect(QtCore.QRectF(0.0, (freq_base-freq_sample/2)/1e6, float(depth), freq_sample/1e6)) 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) class MainWindow(pg.GraphicsWindow): def __init__(self, freq_sample, freq_base, block_size): pg.GraphicsWindow.__init__(self) self.setWindowTitle("NOPTICA Wavemeter") self.freq_sample = freq_sample self.freq_base = freq_base self.text = pg.LabelItem() self.addItem(self.text) p1 = self.addPlot(row=1, col=0) self.beat_spectrum = SpectrogramItem(freq_sample, freq_base, block_size) p1.addItem(self.beat_spectrum) def update_params(self, peak_freq): if peak_freq is None: self.text.setText("REF: NO BEAT") else: self.text.setText("REF: {:.3f}MHz".format((self.freq_base + self.freq_sample*peak_freq)/1e6)) 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"] if action == "update_beat_spectrum": main_window.beat_spectrum.add_block(obj["block"]) main_window.update_params(obj["peak_freq"]) if action == "terminate": self.close_cb() return 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.show() 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()