AskOverflow.Dev

AskOverflow.Dev Logo AskOverflow.Dev Logo

AskOverflow.Dev Navigation

  • Início
  • system&network
  • Ubuntu
  • Unix
  • DBA
  • Computer
  • Coding
  • LangChain

Mobile menu

Close
  • Início
  • system&network
    • Recentes
    • Highest score
    • tags
  • Ubuntu
    • Recentes
    • Highest score
    • tags
  • Unix
    • Recentes
    • tags
  • DBA
    • Recentes
    • tags
  • Computer
    • Recentes
    • tags
  • Coding
    • Recentes
    • tags
Início / coding / Perguntas / 79552366
Accepted
ed190
ed190
Asked: 2025-04-03 16:26:26 +0800 CST2025-04-03 16:26:26 +0800 CST 2025-04-03 16:26:26 +0800 CST

Bloco Python GNU Radio: como armazenar um arquivo criado após executar o fluxograma?

  • 772

Criei um bloco Python que recebe cinco entradas do usuário e gera um arquivo .txt contendo registros de data e hora na coluna 1 e dados de desvio Doppler na coluna 2 com base nas entradas.

As duas primeiras entradas são elementos de duas linhas, a terceira é a frequência de um CubeSat, a quarta é um registro de data e hora e a última especifica o destino do arquivo onde o arquivo .txt será armazenado, semelhante ao bloco "File Sink".

Por exemplo, no bloco padrão "File Sink" do GNU Radio, clicar nos três pontos abre uma janela para navegar e selecionar o destino final do arquivo. No entanto, não consegui replicar esse comportamento no meu código.

Além disso, é possível executar este bloco antes que todo o fluxograma comece ao clicar em "Executar o fluxograma"?

Exemplo: Etapa 1. Entrada 1: "1 49263U 21088D 24308.92471420 .00079255 00000+0 26750-2 0 9998"

Entrada 2: "2 49263 97.4796 14.7363 0010511 313.1417 46.8944 15.30328253170926"

Entrada 3: 437.250e6

Entrada 4: 1730713098

Entrada 5: ---clique nos três pontos para abrir o navegador de arquivos e salvar o arquivo---

Etapa 2. Clique em "Executar o fluxograma" ---executa o bloco python antes do fluxograma principal---

import scipy.constants
import skyfield.api
from skyfield.api import wgs84, EarthSatellite
import numpy as np
import datetime
from tkinter import filedialog
import tkinter as tk
from gnuradio import gr

class DopplerBlock(gr.sync_block):
  
    
    def __init__(self, tle_line1='', tle_line2='', frequency=0, timestamp=0):
        gr.sync_block.__init__(
            self,
            name='Doppler Shift Generator',
            in_sig=None,
            out_sig=None
        )
        
        self.tle_line1 = tle_line1
        self.tle_line2 = tle_line2
        self.frequency = frequency
        self.timestamp = timestamp  # User-provided timestamp
        self.ts = skyfield.api.load.timescale()
        self.groundstation = wgs84.latlon(53.1111, 8.8583, 0)  # ESTEC ground station
        self.output_file_path = self.get_file_path()

    def get_file_path(self):
        """Opens a file dialog to let the user choose a file."""
        root = tk.Tk()
        root.withdraw()  # Hide the root window
        file_path = filedialog.asksaveasfilename(title="Select file to save data")
        return file_path if file_path else None

    def compute_doppler(self):
        if not self.tle_line1 or not self.tle_line2:
            raise RuntimeError("Both TLE lines must be provided.")
        
        satellite = EarthSatellite(self.tle_line1, self.tle_line2, 'satellite', self.ts)
        
        unix_epoch = datetime.datetime(1970, 1, 1, tzinfo=datetime.timezone.utc)
        t0 = unix_epoch + datetime.timedelta(seconds=self.timestamp)
        t0 = self.ts.from_datetime(t0)

        # Generate time steps
        duration_s = 25 * 60  # 25 minutes in seconds
        time_steps = np.arange(0, duration_s, 0.1)  # Every 0.1 sec
        t = t0 + time_steps / (24 * 3600)  # Convert seconds to fractional days
        
        # Compute Doppler shift
        topocentric = satellite.at(t).observe(self.groundstation)
        _, _, range_rate = topocentric.apparent().velocity.km_per_s
        range_rate *= 1e3  # Convert to m/s

        doppler_shift = -range_rate / scipy.constants.c * self.frequency

        # Save to file if path is provided
        if self.output_file_path:
            with open(self.output_file_path, 'w') as output_file:
                for time, shift in zip(time_steps, doppler_shift):
                    output_file.write(f'{self.timestamp + time}\t{shift}\n')

    def work(self, input_items, output_items):
        """This block does not process streaming signals, so work does nothing."""
        return 0
gnuradio
  • 2 2 respostas
  • 25 Views

2 respostas

  • Voted
  1. Best Answer
    Marcus Müller
    2025-04-03T16:51:48+08:002025-04-03T16:51:48+08:00

    Esse não é o comportamento do seu código python, mas é especificado no arquivo YAML que define o bloco; tutorial explicando isso . Se você ainda não fez isso, vale a pena ler os tutoriais em https://tutorials.gnuradio.org de "O que é GNU Radio" a "Criando um OOT (exemplo de bloco Python)" em ordem linear; escrevemos tudo isso pensando em usuários como você :)

    Você pode encontrar o arquivo YAML "Definição de bloco" do File Sink no link abaixo da página de documentação do bloco ; procure dtype:no primeiro elemento de parameters:.

    Acho que já mencionei isso antes, mas você definitivamente não quer uma subclasse real de gr.blockse não estiver fazendo processamento de sinal. E, novamente, nenhuma criação de interface gráfica de usuário em seus próprios blocos :)

    • 1
  2. Praveen Wijesinghe
    2025-04-03T17:42:22+08:002025-04-03T17:42:22+08:00

    Para executar um código python antes de começar, uma maneira fácil é usar o bloco Python Snippet no GNU radio. Selecione section of Flow graphpara After Initou After Start conforme necessário. Coloque o código paraCode Snippethttps://wiki.gnuradio.org/index.php/Python_Snippet

    • 0

relate perguntas

Sidebar

Stats

  • Perguntas 205573
  • respostas 270741
  • best respostas 135370
  • utilizador 68524
  • Highest score
  • respostas
  • Marko Smith

    Reformatar números, inserindo separadores em posições fixas

    • 6 respostas
  • Marko Smith

    Por que os conceitos do C++20 causam erros de restrição cíclica, enquanto o SFINAE antigo não?

    • 2 respostas
  • Marko Smith

    Problema com extensão desinstalada automaticamente do VScode (tema Material)

    • 2 respostas
  • Marko Smith

    Vue 3: Erro na criação "Identificador esperado, mas encontrado 'import'" [duplicado]

    • 1 respostas
  • Marko Smith

    Qual é o propósito de `enum class` com um tipo subjacente especificado, mas sem enumeradores?

    • 1 respostas
  • Marko Smith

    Como faço para corrigir um erro MODULE_NOT_FOUND para um módulo que não importei manualmente?

    • 6 respostas
  • Marko Smith

    `(expression, lvalue) = rvalue` é uma atribuição válida em C ou C++? Por que alguns compiladores aceitam/rejeitam isso?

    • 3 respostas
  • Marko Smith

    Um programa vazio que não faz nada em C++ precisa de um heap de 204 KB, mas não em C

    • 1 respostas
  • Marko Smith

    PowerBI atualmente quebrado com BigQuery: problema de driver Simba com atualização do Windows

    • 2 respostas
  • Marko Smith

    AdMob: MobileAds.initialize() - "java.lang.Integer não pode ser convertido em java.lang.String" para alguns dispositivos

    • 1 respostas
  • Martin Hope
    Fantastic Mr Fox Somente o tipo copiável não é aceito na implementação std::vector do MSVC 2025-04-23 06:40:49 +0800 CST
  • Martin Hope
    Howard Hinnant Encontre o próximo dia da semana usando o cronógrafo 2025-04-21 08:30:25 +0800 CST
  • Martin Hope
    Fedor O inicializador de membro do construtor pode incluir a inicialização de outro membro? 2025-04-15 01:01:44 +0800 CST
  • Martin Hope
    Petr Filipský Por que os conceitos do C++20 causam erros de restrição cíclica, enquanto o SFINAE antigo não? 2025-03-23 21:39:40 +0800 CST
  • Martin Hope
    Catskul O C++20 mudou para permitir a conversão de `type(&)[N]` de matriz de limites conhecidos para `type(&)[]` de matriz de limites desconhecidos? 2025-03-04 06:57:53 +0800 CST
  • Martin Hope
    Stefan Pochmann Como/por que {2,3,10} e {x,3,10} com x=2 são ordenados de forma diferente? 2025-01-13 23:24:07 +0800 CST
  • Martin Hope
    Chad Feller O ponto e vírgula agora é opcional em condicionais bash com [[ .. ]] na versão 5.2? 2024-10-21 05:50:33 +0800 CST
  • Martin Hope
    Wrench Por que um traço duplo (--) faz com que esta cláusula MariaDB seja avaliada como verdadeira? 2024-05-05 13:37:20 +0800 CST
  • Martin Hope
    Waket Zheng Por que `dict(id=1, **{'id': 2})` às vezes gera `KeyError: 'id'` em vez de um TypeError? 2024-05-04 14:19:19 +0800 CST
  • Martin Hope
    user924 AdMob: MobileAds.initialize() - "java.lang.Integer não pode ser convertido em java.lang.String" para alguns dispositivos 2024-03-20 03:12:31 +0800 CST

Hot tag

python javascript c++ c# java typescript sql reactjs html

Explore

  • Início
  • Perguntas
    • Recentes
    • Highest score
  • tag
  • help

Footer

AskOverflow.Dev

About Us

  • About Us
  • Contact Us

Legal Stuff

  • Privacy Policy

Language

  • Pt
  • Server
  • Unix

© 2023 AskOverflow.DEV All Rights Reserve