[Python-es] Hola a todos // Software para Radiointerferometro

José Sabater Montes jsm en iaa.es
Vie Abr 4 02:34:06 CEST 2014 

Hola Rolando,

Un proyecto muy bonito y también muy complejo.

Antes de mirar lo del software tengo una serie de preguntas que me parecen relevantes:
¿Cuántas antenas? ¿de qué tipo? ¿en qué configuración? ¿a qué frecuencia? ¿qué ancho de banda? ¿una o dos polarizaciones? ¿está el correlador ya listo? ¿qué correlador es?


El 03/04/2014, a las 23:14, Rolando Paz <flxpaz en gmail.com> escribió:

> Aun no entiendo como se organiza esto como proyecto. Alguien tiene idea?
> 
> En el siguiente link http://www.dmradas.co.uk/Downloads.html pueden descargar lo que considero son las mejores presentaciones que he visto respecto al tema de la radioastronomia e radiointerferometria.
> 
> Específicamente les sugiero que vean estas presentaciones:
> 
> http://www.dmradas.co.uk/Presentation%20Part%202A.ppt
> http://www.dmradas.co.uk/Presentation%20Part%202B.ppt
> 
> El script de python que se observa aca  https://casper.berkeley.edu/svn/trunk/projects/pocketcorrelator/poco-0.1.1/scripts/poco_rx_i8_c256.py, obtiene los datos del correlador en un formato de archivos llamado MIRIAD.

El formato MIRIAD se usa principalmente en CARMA (https://en.wikipedia.org/wiki/Combined_Array_for_Research_in_Millimeter-wave_Astronomy) y ATCA (https://en.wikipedia.org/wiki/Australia_Telescope_Compact_Array). En el archivo se guardan las “visibilidades” que obtiene el correlador (el resultado de la correlación) e información sobre las antenas, los receptores, etc.
http://www.atnf.csiro.au/computing/software/miriad/userguide/node224.html#app:uvvars

> 
> Estos datos los almaceno, para luego ser procesados con un software como , http://www.atnf.csiro.au/computing/software/miriad/, o como https://github.com/AaronParsons/aipy, o https://casper.berkeley.edu/astrobaki/index.php/AIPY

Tengo entendido que Miriad se usa para calibrar datos de interferómetros de disco y AIPY está preparado para calibrar datos de “arrays” de dipolos. Son programas bastante complejos y específicos y es posible que no funcionen directamente para este proyecto (o quizás sí).

> 
> También me gustaría poder desarrollar una interfaz gráfica para poder ver los espectros de las señales analógicas por medio de matplotlib o alguna otra herramienta en tiempo real. Aca se puede ver un ejemplo de script para esto, sin embargo esta programado para obtener los datos por el puerto RS232 y no ethernet: https://casper.berkeley.edu/svn/trunk/projects/pocketcorrelator/wb_poco_ibob/Python_Script/

Por el mensaje de error de antes veo que se requieren 73728 bytes de almacenamiento por cada secuencia de integración de poco menos de un segundo. Parece que eso está por encima de los 65536 bytes máximos de algún buffer y por eso puede que aparezca un fallo como indica José Sabastizagal.

Por otro lado, 73728 bytes se corresponden con 18432 parámetros de 4 bytes. Y esto encaja con 256 canales, 2 polarizaciones y 8 antenas con sus cross-correlaciones (28 lineas de base + 8 cross-correlaciones = 36). ¿Es esto correcto? Quizás podrías reducir el número de antenas para la prueba para ver si sigue dando el mismo fallo (?). Con 7 antenas deberías tener integraciones de 57344 bytes.

Un saludo,
Jose Sabater

> 
> Saludos y gracias por todo el apoyo mostrado.
> 
> Rolando Paz
> 
> 
> 
> 
> 
> 
> 
> 
> 
> 
> 
> 
> El 3 de abril de 2014, 15:24, Flavio Danesse <fdanesse en gmail.com> escribió:
> Yo me anoto en el proyecto.
> 
> No tengo la mas mínima idea de que es un interferómetro ni de radioimágenes ni nada de eso, pero programo sobre todo en python, gtk y gstreamer y a veces en otros lenguajes, capaz que puedo ayudar con el código.
> 
> De momento lo voy a mirar nada más.
> Creo que vas a necesitar organizar como proyecto esto porque ya son varios anotados para participar, alguien tiene que organizar la cosa.
> 
> 
> El 3 de abril de 2014, 13:44, Rolando Paz <flxpaz en gmail.com> escribió:
> 
> El vídeo en donde se explica el funcionamiento del correlador es el siguiente:
> 
> https://casper.berkeley.edu/presentations/workshop_october09/m4v/gfoster.m4v
> 
> Griffin Foster fue quien diseño el software del correlador.
> 
> Saludos
> 
> Rolando Paz
> 
> 
> El 3 de abril de 2014, 10:35, Rolando Paz <flxpaz en gmail.com> escribió:
> Muchas gracias a todos!
> 
> El problema que tengo es el siguiente:
> 
> Yo instale wiresahark y estoy usando ubuntu 12.04 LST
> 
> Ejecuto el comando:  init_poco.py i8_c256.py
> 
> rolando en rolando-MS-7815:~/ibob/poco011$ init_poco.py i8_c256.config 
> 
> Parsing config file i8_c256.config...OK
> 
> Writing override variables... OK
> 
> Connecting to iBOB... OK
> 
> Adjunto el archivo "init_poco.py - i8_c256.config", que es la captura de los datos de wireshark (una imagen wireshark.png), en donde se observa la información que se transmite desde mi PC hacia el correlador llamado IBOB. Sin embargo, luego de ejecutar este comando, creo que la IBOB debería empezar a enviar automáticamente los paquetes, pero no lo hace pues no se observa nada en el wireshark después.
> 
> Aun así, yo ejecuto el siguiente comando para obtener los datos en formato MIRIAD, sudo poco_rx_i8_c256.py i8_c256.config
> 
> rolando en rolando-MS-7815:~/ibob/poco011$ sudo poco_rx_i8_c256.py i8_c256.config 
> 
> Parsing config file i8_c256.config...OK
> Listening on port 7
> Expecting integration sizes of 18432 vectors, or 73728 bytes
> C2M Parameters:
> 	N Antennas: 8
> 	Bandwidth: 0.200000 GHz
> 	SDF: 0.000781 GHz
> 	Int Time: 0.747520 s
> 	Array location:  ['38:25:59.24', '-79:51:02.1']
> Recording Bandpass to file... OK
> Starting file: zen.uv.tmp
> Beginning RX thread...
> Expecting total dump size from each x engine: 65536 bytes
> ERR: Buffer is too small for header unpack
> ERR: could not unpack header
> Exception in thread Thread-1:
> Traceback (most recent call last):
>   File "/usr/lib/python2.7/threading.py", line 551, in __bootstrap_inner
>     self.run()
>   File "/usr/lib/python2.7/threading.py", line 504, in run
>     self.__target(*self.__args, **self.__kwargs)
>   File "/usr/local/bin/poco_rx_i8_c256.py", line 243, in _process_packets
>     last_offset = p['offset']
> TypeError: 'NoneType' object has no attribute '__getitem__'
> 
> Mi duda es la siguiente:
> 
> a) Al ejecutar el comando init_poco.py i8_c256.config, yo le digo a la IBOB que empiece a enviar los datos de forma automática por medio del puerto ethernet? (los archivos están en el link: https://casper.berkeley.edu/svn/trunk/projects/pocketcorrelator/poco-0.1.1/)
> 
> b) Qué significa el error que se produce cuando ejecuto el comando sudo poco_rx_i8_c256.py i8_c256.config?
> 
> Saludos y gracias.
> 
> Rolando Paz 
> 
> 
> 
> 
> 
> 
> 
> 
> 
> El 3 de abril de 2014, 9:50, Alfonso Ángel Rodríguez Escudero <arodesc en gmail.com> escribió:
> 
> Soy Ingeniero en Telecomunicación y entre otras cosas me dedico a trastear con estas cosillas de antenas. También le echaré un ojo y vamos viendo a ver que sale. ;D
> Saludos!!!
> 
> El 03/04/2014, a las 17:35, Rolando Paz <flxpaz en gmail.com> escribió:
> 
>> Ahora mismos he dedicado todo el tiempo para construir el back-end del interferómetro, y esto incluye el correlador y el software. Luego de esto debo diseñar el front-end, que en esencia es toda la parte analógica del interferómetro: amplificadores de bajo ruido, convertidores de frecuencia, filtros, alimentadores de las antena, cables y otros...
>> 
>> Me falta mucho, pero sigo adelante...
>> 
>> 
>> 
>> 
>> 
>> El 3 de abril de 2014, 8:41, Daπid <davidmenhur en gmail.com> escribió:
>> A mí también me interesa. Soy físico, con algunos "antecedentes" en astrofísica. Y si tú tienes antenas, o experiencia construyendo, mejor. :)
>> 
>> 
>> 2014-04-02 23:08 GMT+02:00 Rolando Paz <flxpaz en gmail.com>:
>> Es un gusto saludarlos.
>> 
>> Desde hace ya algún tiempo estoy diseñando un correlador para obtener las señales de varias antenas y lograr realizar radioimagenes del universo. 
>> 
>> Ahora mismo estoy tratando de entender el software desarrollado en el 2008 para dicho fin. Desde entonces ya no hubo actualizaciones y ya no existe soporte. El software son unos script escritos en python. 
>> 
>> Alguno de ustedes estaría interesado en ayudarme a entender y a actualizar los script?
>> 
>> Yo sé qué es lo que deben hacer los script, pero no logro hacer que funcionen.
>> No está demás decir que yo no sé programar :-)
>> 
>> Estaré muy agradecido por su ayuda. 
>> 
>> Saludos
>> 
>> Rolando Paz
>> 
>> _______________________________________________
>> Python-es mailing list
>> Python-es en python.org
>> https://mail.python.org/mailman/listinfo/python-es
>> FAQ: http://python-es-faq.wikidot.com/
>> 
>> 
>> 
>> _______________________________________________
>> Python-es mailing list
>> Python-es en python.org
>> https://mail.python.org/mailman/listinfo/python-es
>> FAQ: http://python-es-faq.wikidot.com/
>> 
>> 
>> _______________________________________________
>> Python-es mailing list
>> Python-es en python.org
>> https://mail.python.org/mailman/listinfo/python-es
>> FAQ: http://python-es-faq.wikidot.com/
> 
> _______________________________________________
> Python-es mailing list
> Python-es en python.org
> https://mail.python.org/mailman/listinfo/python-es
> FAQ: http://python-es-faq.wikidot.com/
> 
> 
> 
> 
> 
> _______________________________________________
> Python-es mailing list
> Python-es en python.org
> https://mail.python.org/mailman/listinfo/python-es
> FAQ: http://python-es-faq.wikidot.com/
> 
> 
> 
> _______________________________________________
> Python-es mailing list
> Python-es en python.org
> https://mail.python.org/mailman/listinfo/python-es
> FAQ: http://python-es-faq.wikidot.com/
> 
> 
> _______________________________________________
> Python-es mailing list
> Python-es en python.org
> https://mail.python.org/mailman/listinfo/python-es
> FAQ: http://python-es-faq.wikidot.com/

Más información sobre la lista de distribución Python-es