SOFTWARE

 

SOL

Descargar

Instalación.-

Se descarga el archivo sol.rar y se descomprime en una carpeta.  El programa no instala ningún archivo adicional ni modifica el registro.

El archivo comprimido incluye el ejecutable, una librería gráfica y la documentación del programa.

Para desinstalarlo, simplemente hay que borrar la carpeta donde se encuentra el ejecutable.

 

Requisitos del sistema.-

- Microsoft .NET Framework 3.5  SP1

Si no se encuentra instalado, puede hacerse desde: http://www.microsoft.com/downloads/details.aspx?familyid=333325FD-AE52-4E35-B531-508D977D32A6&displaylang=es

- Windows Installer 3.1

Puede descargarse desde: http://www.microsoft.com/downloads/details.aspx?FamilyID=889482FC-5F56-4A38-B838-DE776FD4138C&displaylang=es

El programa se ha probado con buenos resultados en Windows XP y Windows 7.

 

Licencia.-

Este programa se distribuye sin licencia ni garantía de ningún tipo. Su utilización se hace bajo la única responsabilidad del usuario.


 

Introducción.-

La medida de posiciones a partir de fotografías se realiza con el programa Iris, según el método descrito en la siguiente página web:

http://www.parhelio.com/docposicfotos.html

El procedimiento es largo y por eso, el programa Sol se diseñó inicialmente para automatizar determinadas tareas (especialmente la orientación de la imagen).

Con posterioridad, el programa se amplió para permitir la medida  de posiciones con retículo, de áreas, de distintos parámetros geométricos de las protuberancias, así como el cálculo de efemérides y movimientos propios de manchas.

En su estado actual, el proceso de medir posiciones con Sol es casi totalmente automático, por lo que apenas es necesario conocer el funcionamiento de Iris.

No obstante, dado que en la mayoría de sus funciones, el programa se pensó para trabajar conjuntamente con Iris, puede ser conveniente que haya una cierta familiaridad con los métodos de trabajo de este programa:

http://www.astrosurf.com/buil/us/iris/iris.htm
http://www.astrosurf.com/buil/iris/iris.htm

En Parhelio pueden encontrarse informaciones más detalladas sobre los procesos realizados con Iris.

 

Configuración del programa.-

Ante todo es necesario configurar algunos parámetros, dependiendo de la parte del programa que vayamos a utilizar. Si solo vamos a usar el programa para generar efemérides o movimientos propios, no es necesario configurar nada.

 

 

La primera vez que lo abramos, el programa creará el archivo “sol.ini” con algunos valores por defecto. Para modificarlos vamos a la pestaña “Configuración” y allí introducimos los siguientes datos:

Direcciones:

- Carpeta de trabajo.- Aquí debemos poner la dirección de la carpeta de trabajo de Iris. Dada la interacción entre Sol e Iris, conviene que ambos compartan la misma carpeta.

- Scripts.- Para realizar algunas funciones, Sol utiliza la capacidad de Iris para trabajar con scripts.  Aquí debe figurar la misma dirección que usa Iris para almacenar los scripts.

Plantilla:

- Tamaño de la imagen (píxeles).- Debemos introducir las dimensiones horizontal y vertical de la imagen, en píxeles. Estas dimensiones son constantes para cada modelo de cámara. Si cargamos una imagen en Iris, el comando >info, nos proporcionará dichas dimensiones. Hay que utilizar los valores de Iris, no los que vienen en las especificaciones de la cámara, ya que en ocasiones no coinciden.

-  Intervalos de longitud y latitud (º).- Cuando Iris dibuja una plantilla sobre el disco del Sol, estos valores proporcionan el espacio en grados, entre meridianos y paralelos adyacentes. Un valor adecuado en la mayoría de casos es de 10º.

Posiciones:

- Código observador.- El código de 3 letras asignado en Parhelio a cada observador.

Areas:

- Tamaño de píxel (micrómetros).- El tamaño de los píxeles de la cámara en micrómetros.  Este dato se suele encontrar en la documentación de la cámara. El dato se utiliza también en los cálculos sobre protuberancias.

Es necesario si vamos a medir áreas de grupos individuales, pero no si solo queremos obtener el área total medida sobre una imagen de todo el disco.

Efemérides:

La fecha y hora se puede introducir en cuatro lugares diferentes: Orientación, Otros métodos, Areas y Efemérides. Como el módulo de efemérides se puede usar en combinación con los otros tres, podemos elegir si queremos que los datos introducidos en Orientación, Otros métodos, o Areas, aparezcan automáticamente en Efemérides, o si queremos que las pestañas permanezcan independientes entre sí.

La opción que marquemos solo será efectiva después de reiniciar el programa.

Archivo manual:

Marcar el método de observación utilizado. Esta opción se usa para definir la precisión de las medidas cuando creamos manualmente el archivo de posiciones.

Iris control :

- Versión Iris.- Marcar si estamos utilizando la versión en francés o en inglés.

- Proyección rectángular.- Durante el proceso para eliminar el oscurecimiento del limbo, el programa tiene la opción de pasar la imagen a proyección rectángular. Los valores introducidos aquí se utilizan para recortar la imagen resultante.  Con Longitud=180 y Latitud=90 (valores máximos), Iris muestra toda la superficie solar, por lo que la mitad aparecerá negra. Con valores Longitud=70 y Latitud=60, el programa nos muestra la superficie a +/- 70º del meridiano central y a +/- 60º de latitud.

- Retardo (ms).- Iris necesita un cierto tiempo para ejecutar los comandos y puede ocurrir que ese tiempo sea mayor de lo que tarda Sol en enviar nuevos comandos. En esos casos, se produce una sobrecarga, originando errores de distinto tipo. El tiempo introducido aquí permite ajustar la velocidad a la que Sol se comunica con Iris. Si el proceso se interrumpe o vemos aparecer errores, aumentemos el retardo.

Las tareas que más tiempo requieren, con diferencia, son la conversión a formato fit y la alineación de las imágenes, así que subir el retardo no va a suponer un aumento importante en el tiempo total, pero así conseguiremos que el proceso sea más fluido.

El valor por defecto es de 800 ms. Con un ordenador rápido se puede disminuir, aunque seguramente necesitaremos procesar varias secuencias de imágenes en diferentes días hasta encontrar un valor óptimo para nuestro equipo.

El tiempo de retardo se utiliza también en otras funciones del programa, cada vez que Sol accede a la ventana de salida o la de estadísticas.

 

Procesado

 

 

La ventana se divide en dos partes: Alineación manual e Iris control. Esta última es la sección  que utilizaremos en condiciones normales.

El material de partida es una secuencia de imágenes de todo el disco, obtenida sin seguimiento. Previamente, con Iris, habremos pasado el formato raw de la cámara a formato fit.

Supondremos primero que la secuencia es suficientemente homogénea, es decir, que no hay fuertes variaciones de intensidad en el fondo de cielo ni en el disco entre las diferentes imágenes (estas variaciones suelen estar producidas por el paso de nubes).

Después veremos como tratar una secuencia con nubes, siempre que éstas no sean demasiado densas.

 

Iris control.-

Este módulo permite automatizar casi totalmente el trabajo con Iris a partir de unos pocos parámetros de partida. No obstante, también necesita de unas condiciones muy precisas, que si no se cumplen pueden hacer surgir errores o interrumpir el proceso:

1.- Ajustar correctamente el tiempo de retardo en la pestaña Configuración.

2.- Al principio, Iris no debe estar minimizado en la barra de tareas, aunque tampoco es necesario que esté en primer plano. Ahora bien, es imprescindible que Iris se encuentre en primer plano mientras esté realizando operaciones. De lo contrario, Sol no podría enviar comandos o leer la ventana de salida.

3.- En Iris, deben estar abiertas la ventana de comandos y la ventana de salida

4.- En la ventana de comandos, la última línea debe tener solo un signo “>”.

5.- Antes de empezar el procesado, tanto en luz blanca como en H-alfa, conviene que en la carpeta de trabajo no haya imágenes creadas en tratamientos anteriores. Se puede utilizar el botón “Borrar archivos temporales”, para eliminarlas.

 

A continuación se describe el proceso para imágenes en luz blanca.

En “Iris control”, ponemos el nombre genérico de las imágenes y podemos indicar si queremos aplicar wavelets. El cursor 1 realzará la capa con los detalles más finos, y el 2, la siguiente. Si los dejamos al mínimo, el programa saltará ese paso.

Aunque habitualmente no modifican el centro fotométrico, es aconsejable no utilizar wavelets salvo que sea necesario, en caso de que haya manchas muy tenues, por ejemplo. Siempre es posible aplicar wavelet posteriormente con Iris, a partir de la imagen “k.fit” que se creará en la carpeta de trabajo.

Una vez comprobado que se cumplen todas las condiciones anteriores, pinchamos en el botón “Alinear”. Las operaciones que realiza este botón son las siguientes:

- Calcula la intensidad del limbo en la imagen central de la secuencia.

- Alinea,

- Suma.

- Ajusta el histograma,

- Obtiene el centro y radio del disco de la imagen resultante.

Al terminar, vamos a la pestaña “Orientación”, y determinamos el ángulo alfa (ver más adelante).

Volvemos a la pestaña “Procesado”, y pinchamos en el botón Orientar. Las operaciones que realiza este botón son las siguientes:

- Rota la imagen y si es necesario, la invierte para colocar el norte solar arriba.

- Calcula el centro y radio del disco de esta nueva imagen.

- Crea una plantilla.

- Aplica wavelet si así lo hemos indicado.

- Recorta y guarda en formato bmp, dos imágenes del disco orientado, una con la plantilla superpuesta y otra sin ella.

Es conveniente examinar la imagen que tiene la plantilla para comprobar que ésta efectivamente se ajusta bien al limbo. Entonces, ya podemos pasar a la pestaña “Posiciones” para realizar las medidas.

 

Con imágenes H-alfa utilizamos los dos botones de la parte inferior, y el proceso es similar, con las siguientes diferencias:

1.- La conversión Raw -> Fit se hace en color, así que, antes de alinear las imágenes el programa separará los 3 canales y eliminará los canales G y B. Este paso puede tardar unos minutos.

2.- El disco y el limbo se procesan por separado; el primero igual que las imágenes en luz blanca, y el segundo, enmascarando el disco.

Al final del procesado H-alfa, se muestra una imagen del limbo que podemos usar para medir latitudes de protuberancias (ver más adelante).

 

En luz blanca, una vez medidas las posiciones podemos eliminar el oscurecimiento del limbo siguiendo el proceso que se describe en: http://www.parhelio.com/doctecnicas.html  Para ello, simplemente hay que introducir el valor de la longitud de onda y pinchar en “Eliminar gradiente”.  En una cámara digital un valor adecuado puede ser 550 nm, aunque si el ajuste no es bueno, probablemente tengamos que probar otros valores.

Si está marcada la casilla “Proyección rectangular”, al final el programa también pasará la imagen a coordenadas rectangulares y la guardará en formato fit y bmp. La imagen tendrá una escala de 0.1º/px y se puede utilizar para hacer planisferios. Las casillas Longitud y Latitud de la pestaña Configuración establecen los límites de esta imagen.

 

El botón “Borrar archivos temporales” elimina todos los archivos creados durante el proceso, excepto los siguientes:

- Las imágenes raw originales.
- El archivo Shift.lst.
- El archivo xy.txt.
- Las dos plantillas sol.lst y solp.lst
- k.fit.- La imagen orientada.
- k2.fit.- La imagen orientada después de aplicar wavelet.
- k3.fit.- La imagen sin oscurecimiento del limbo.
- aammdd-hhmm.bmp.- Imagen orientada después de aplicar wavelet y recortada.
- aammddcoord.bmp.- Igual que la anterior, pero con plantilla superpuesta.
- aammddsynthe.bmp.- Imagen sin oscurecimiento del limbo y recortada.
- aammdd-hhmmdisco.bmp.- Imagen Ha del disco, orientada, después de aplicar wavelet y recortada.
- aammddcoorddisco.bmp.- Igual que la anterior, pero con plantilla superpuesta.
- aammdd-hhmmlimbo.bmp.- Imagen Ha del limbo, orientada, después de aplicar wavelet y recortada.
- aammddcoordlimbo.bmp.- Igual que la anterior, pero con plantilla superpuesta.
- aammddpl.bmp, aammddpl.fit.- Imágenes en coordenadas rectangulares.
- aammdd medidas OBS.txt.- El archivo de medidas.

 

 “Iris control” puede fallar ocasionalmente con otros procesos activos. Si esto no ocurre de forma sistemática, a veces es suficiente con volver a repetirlo desde el último paso antes de producirse el fallo. Otra opción consiste en aumentar el tiempo de retardo, en la pestaña “Configuración”.

Si volvemos a ejecutar el botón Orientar, recordemos que es necesario obtener previamente el ángulo alfa.

 

Alineación manual .-

En ocasiones, debido sobre todo a la presencia de nubes, Iris puede encontrar dificultades para alinear las imágenes. En estos casos, una opción es calcular el centro fotométrico de una mancha que aparezca en todas las imágenes, y con dichos valores, construir manualmente el archivo Shift.lst. (Ver “Medidas de posición.pdf”). Aunque, en teoría, dos medidas son suficientes, lo aconsejable es utilizar más de tres, y en cualquier caso, hay que usar el máximo posible. El programa solo realizará los cálculos a partir de tres medidas.

El método consiste en abrir la primera imagen, elegir una mancha que se vea en toda la secuencia, recuadrarla, y con el botón derecho del ratón abrir la ventana de estadísticas. Sin cerrarla, pinchamos en el botón “Adquirir coordenadas”, y automáticamente, obtiene las coordenadas del centro fotométrico y las incorpora al listado. Este método requiere un ajuste del tiempo de retardo, en la pestaña Configuración (ver también la explicación sobre “Iris control”).

También podemos introducirlas manualmente en las casillas X, Y, y con el botón “Añadir” se incorporan a la lista. El proceso se repite con todas las imágenes.

Las coordenadas que introducimos son “absolutas” dado que el origen se encuentra en la esquina inferior izquierda de las imágenes, y aunque podemos usarlas para calcular a, Iris no puede usarlas para alinear.

El botón “Crear Shift.lst” sirve para crear un archivo Shift válido, con coordenadas “relativas”, utilizando como referencia la primera imagen. Después, el botón “Alinear” realiza las mismas operaciones que el de “Iris control”, salvo que utiliza un método diferente para alinear las imágenes

A partir de este punto, ya podemos continuar el proceso calculando primero el ángulo alfa en la pestaña “Orientación”, y después con el botón orientar de “Iris control”.

 

Orientación

Cuando Iris se dispone a alinear una secuencia de imágenes, en su carpeta de trabajo se crea el archivo Shift.lst, que contiene los desplazamientos sufridos por cada imagen para hacerla coincidir con la de referencia. Iris utiliza dicho archivo para hacer la alineación. Si la alineación se ha realizado con "Iris control" (comando >cregister), esos desplazamientos cambiados de signo, coinciden con las coordenadas del centro del disco solar en cada imagen (considerando como origen la imagen de referencia). Por tanto, la línea definida por dichos puntos nos marca la dirección E-W, y podemos usarla para orientar la imagen.

 

 

Una vez que Iris ha alineado la secuencia de imágenes, en el programa Sol vamos a la pestaña Orientación, e introducimos la fecha y hora de la observación. A medida que escribimos los datos, el programa calcula las efemérides para ese momento:

D.- Diámetro solar en segundos de arco.

P.- Inclinación del eje de rotación solar, respecto a la dirección N-S.

B.- Latitud heliográfica del centro del disco solar.

L.- Longitud heliográfica del meridiano central.

Si a continuación pinchamos en el botón “Shift.lst”, ocurren varias cosas. A la izquierda aparece el contenido del archivo y, simultáneamente, se hace un ajuste lineal y se calcula el ángulo que hemos de girar la imagen para que el norte solar quede arriba.

Asimismo, en la gráfica se representa el centro del disco solar en toda la secuencia de imágenes. Los puntos azules indican aquellos cuya diferencia con el ajuste es menor que la desviación típica (e < s); los de color naranja son puntos cuya diferencia con el ajuste se encuentra entre la desviación típica y el doble de dicho valor (s < e < 2s); y los de color rojo, aquellos cuya diferencia es más del doble de la desviación típica (e > 2s). La recta es el resultado del ajuste por mínimos cuadrados.

La escala de los ejes se ajusta automáticamente a cada caso. Esto puede suponer que si el ángulo entre la dirección E-W y el eje X es alto, la dispersión aparente en la gráfica puede parecer menor que con un ángulo bajo. Por ello, para constatar la fiabilidad de los datos, no debemos guiarnos por la gráfica, sino por los valores del coeficiente de correlación.

Inicialmente se incluyen todos los valores en el cálculo. No obstante, en previsión de que alguna de las imágenes haya quedado mal alineada, podemos descartar los valores más discrepantes en función de la desviación típica. Debajo del botón “Shift.lst” tenemos:

<1s.- Solo se consideran aquellos valores cuya diferencia con el ajuste sea menor que la desviación típica.

<2s.- Solo se consideran aquellos valores cuya diferencia con el ajuste sea menor que el doble de la desviación típica.

<3s.- Solo se consideran aquellos valores cuya diferencia con el ajuste sea menor que el triple de la desviación típica.

Generalmente, un valor de 2s suele ir bien en la mayoría de los casos.

Además aparecen los siguientes datos:

N.- Número de datos considerado.

s.- Desviación típica.

r.- Coeficiente de correlación.

a.- Angulo que hemos de girar la imagen para que el norte solar quede arriba. Este ángulo incluye 2 correcciones: la desviación de la cámara respecto a la dirección E-W, y el ángulo P.  a  es el valor que hemos de introducir como parámetro en el comando >rot de Iris para orientar la imagen adecuadamente.

Se considera por defecto, que en la secuencia original de imágenes, el Norte está arriba y el Este a la izquierda, es decir, que el Sol recorre el campo de izquierda a derecha. Ahora bien, puede ocurrir que, debido a la configuración óptica o a la orientación de la cámara, esto no sea así, de manera que es posible que tengamos que invertir la imagen vertical u horizontalmente.

En la parte inferior podemos marcar en qué caso estamos, y según la opción, veremos aparecer, debajo de a, uno o dos comandos que tendremos que aplicar en Iris después de rotar la imagen, para conseguir la orientación correcta.

 

Como hemos visto anteriormente, una vez orientada la imagen, se calculan las coordenadas del centro y el radio del disco. Esos valores aparecen en las casillas Xc, Yc, R.

Después, el programa obtiene automáticamente un archivo de texto con los parámetros que definen la red de coordenadas. No obstante, con el botón “Crear plantilla” tenemos la opción de crear manualmente ese archivo.

Finalmente, con el botón “Borrar” podemos eliminar el contenido de las ventanas, lo cual permite tratar varias secuencias de imágenes sin necesidad de reiniciar el programa en cada ocasión.

 

Posiciones

Esta pestaña se utiliza para medir las posiciones una vez obtenida la imagen orientada.

 

 

En dicha imagen, se elige una mancha, se recuadra, y con el botón derecho del ratón podemos abrir la ventana de estadísticas. Si pinchamos en el botón “Obtener posiciones”, Sol leerá en esa ventana las coordenadas del centro fotométrico y calculará automáticamente las coordenadas heliográficas.

También se pueden introducir manualmente las coordenadas (x, y) de las manchas, en las casillas correspondientes, y con el botón “Calcular” se obtienen las coordenadas heliográficas.

A continuación, hay que asignar cada posición a una mancha. Para ello, en la segunda ventana se introduce manualmente la numeración NOAA del grupo, seguida de una o dos letras que identifican cada mancha.

El número de caracteres que admite el programa es de 5 dígitos para el grupo y entre cero y tres letras para la mancha. Si el grupo aún no tiene identificación, cabe la posibilidad de introducir solo las letras, o dejar la entrada vacía.

Una vez que hayamos acabado, con el botón “Crear archivo” las medidas aparecen en la ventana de la derecha con el formato utilizado en Parhelio. Se añade en la parte inferior, el número de rotación y las efemérides correspondientes. El botón “Guardar” nos guarda un archivo de texto con el contenido de esa ventana.

Con el botón “Borrar” eliminamos el contenido de las ventanas, lo cual permite tratar varias secuencias de imágenes sin necesidad de reiniciar el programa en cada ocasión.

Si hemos usado un retículo, con el botón “Medidas” cargamos en la ventana izquierda, el archivo creado en la pestaña “Otros métodos”.

Cada vez que calculamos una posición, las coordenadas (X, Y) de la mancha se escriben en un archivo denominado “xy.txt”, en la carpeta de trabajo. El programa no utiliza dicho archivo, pero contiene todos los datos necesarios para calcular las posiciones y resulta de utilidad si más adelante queremos reprocesar la imagen o identificar algún posible error en las medidas. El botón “Importar xy” lee dicho archivo y muestra las posiciones en la ventana izquierda.

 

En la parte inferior, podemos hacer unos cálculos adicionales en grupos de manchas bipolares.

Si en la ventana derecha seleccionamos dos manchas y pinchamos en el botón “Calcular”, el programa obtiene la posición media, la distancia entre ambas, en grados y km, así como la inclinación respecto al ecuador de la línea que las une.

Los cálculos se hacen a lo largo del círculo máximo que pasa por ambas manchas. La inclinación i, se considera positiva si la latitud de la mancha f es mayor que la de la mancha p, y negativa en el caso contrario.

Con el botón “Abrir archivo” podemos cargar en la ventana un archivo de medidas creado anteriormente por Sol, o con su mismo formato, para realizar sobre él esas operaciones.

El botón “Incluir” añade estos cálculos al archivo final.

 

Areas

El método para obtener áreas se explica en http://www.parhelio.com/docareas.html

 

 

A la izquierda se pueden introducir 8 focales diferentes que se guardan automáticamente entre las preferencias. De esa manera, si un día usamos una configuración óptica diferente, solo tendremos que elegir la que corresponda.

Una vez introducidos el día y la hora, seleccionamos la focal,  identificamos el grupo (numeración NOAA),  su longitud y latitud media, y el área en píxeles. Este último dato es el que nos proporciona Iris. Con el botón “Calcular “ obtenemos el área en millonésimas de hemisferio (mh); y con “Añadir” incorporamos la medida al listado, con el formato utilizado en Parhelio. Finalmente, guardamos un archivo de texto con el contenido de la ventana.

 

No es recomendable medir áreas de grupos individuales sobre imágenes de todo el disco, porque los errores pueden ser importantes. La medida del área total es mucho más fiable, especialmente si se combinan los datos de varios observadores y se obtienen promedios mensuales.

En la parte inferior de la ventana se puede calcular el área total cubierta por las manchas, a partir de una imagen de todo el disco. Para ello, es necesario que en la carpeta de trabajo se encuentren la imagen sin el oscurecimiento del limbo (k3.fit), y el archivo con los parámetros de la plantilla asociada a esa imagen (sol.lst).

Con el botón “Procesar”, el programa realiza una serie de operaciones cuyo resultado es una imagen con el fondo blanco y unas manchas negras. Esas manchas coinciden con las manchas solares, aunque también pueden aparecer otras manchas (generalmente simples puntos negros) debidas a la granulación y a una cierta componente de ruido. Si el valor que hay en la casilla es muy alto, las zonas negras pueden cubrir gran parte del disco.

La forma de operar depende de cómo sea la imagen:

- Si aparecen los dos tipos de manchas, disminuimos un paso el valor de la casilla, y el botón “Ajustar” vuelve a calcular la imagen. Notaremos que los puntos debidos a la granulación disminuyen. Hay que ir reduciendo paso a paso ese valor hasta que solo queden las manchas solares (no importa si persiste algún punto aislado). En ese momento, con el botón “Calcular” obtenemos el área.

- Si solo aparecen las manchas solares, vamos subiendo el valor de la casilla hasta que aparezca la granulación. Entonces lo reducimos un paso y calculamos el área.

Por ejemplo, vamos subiendo (10000-10500-11000-11500-12000), y al llegar a 12000 surge la granulación. Entonces, el valor que nos sirve es 11500.

En la casilla, el valor se modifica en pasos de 500, pero si se escribe manualmente, se puede poner cualquier otro valor y ajustarlo de forma más precisa; aunque normalmente, esto no es necesario.

En ocasiones, para facilitar el ajuste es mejor haber aplicado wavelets. Esto puede hacerse con los dos cursores de la pestaña “Procesado”; o bien, antes de eliminar el oscurecimiento del limbo. En este caso, se abre la imagen k.fit, se aplican con el menú Processing, Wavelet (Traitement, Ondelettes), y se graba el resultado con el nombre k2.fit. Después ya se puede eliminar el gradiente y calcular el área de la manera habitual.

 

El área total se calcula en millonésimas de disco (md), al contrario que el área de grupos individuales, que se obtiene en millonésimas de hemisferio (mh); y por tanto, no son comparables.

 

Protuberancias

Esta pestaña permite realizar  diversos cálculos geométricos sobre las protuberancias.

Para obtener latitudes y alturas, utilizaremos imágenes H-alfa de todo el disco, obtenidas y procesadas con el mismo método que las imágenes en luz blanca.

El apartado “Velocidades” está pensado para ser utilizado con animaciones, o al menos, con imágenes que hayan sido obtenidas con la misma configuración óptica, con las mismas dimensiones y orientación, y que hayan sido previamente alineadas.

Se supone por defecto, que para el cálculo de velocidades se utiliza una cámara de vídeo cuyos píxeles tienen la dimensión especificada en la pestaña Configuración, apartado Areas.

 

 

Latitudes.-

Una protuberancia es una estructura extensa que generalmente abarca un sector más o menos grande del limbo. En Iris, utilizaremos el cursor para determinar las coordenadas (X, Y) de los extremos de dicho sector, y el programa calculará la latitud media, el ángulo de posición (medido desde el N, en sentido N-E-S-W), y la extensión de la protuberancia en grados heliográficos. Estos valores se añaden a la ventana, pudiendo posteriormente guardarlos en un archivo de texto.

El método consiste en aprovechar los datos que aparecen en la ventana de salida (“Output”) de Iris. Para operar de esta manera, primero abrimos la ventana (si no lo está) en el menú Analysis, Display data  (Analyse, Afficher les données), y borramos su contenido en caso de que haya algún dato presente (menú Edit, Erase, o en francés, Edition, Effacer, en la propia ventana).

Cada vez que pinchamos en la imagen, en la ventana aparecen tres datos: las coordenadas del píxel donde hemos pinchado, y su intensidad. El método consiste en ir marcando uno tras otro, con el ratón, los extremos de cada protuberancia. Cuando hayamos acabado, grabaremos el contenido de la ventana, con el nombre “medidas.dat”.

Cuando en Sol, pinchemos en el botón “Medidas.dat”, el programa leerá el contenido del archivo, y calculará el ángulo de posición, la latitud y la extensión de cada protuberancia.

Otro método, más rápido, es similar al anterior, pero en vez de grabar el archivo “Medidas.dat”, hay que dar al botón “Obtener medidas”. De esta manera, el programa lee el contenido de la ventana de salida, y automáticamente calculará las latitudes y extensiones de las protuberancias. Para utilizar esta opción, es necesario que delante de Iris solo se encuentre la ventana de Sol. Asimismo, requiere un ajuste del tiempo de retardo, en la pestaña Configuración (ver también la explicación sobre “Iris control”)

on ambos métodos, el programa creará automáticamente un archivo denominado “xyp.txt”, con los datos necesarios para volver a realizar los cálculos posteriormente si fuese necesario. El botón “Importar xyp” permite realizar dicha operación.

 

Alturas.-

De la misma manera que en el caso anterior, podemos obtener en Iris las coordenadas (X, Y) del punto más alto de la protuberancia. Dichas coordenadas se pueden introducir en Sol y obtener la altura sobre la superficie, en km y diámetros terrestres.

Tanto en el cálculo de latitudes como de alturas, es preciso que en la pestaña “Orientación” figuren las coordenadas del centro y el radio del disco, correspondientes a la imagen que estamos analizando. Esto no es necesario si importamos el archivo xyp.txt, porque éste ya contiene dichos valores.

 

Velocidades.-

Lo ideal para determinar velocidades, es aprovechar una secuencia de imágenes destinada a elaborar una animación de una protuberancia. Las imágenes, por tanto, tienen que cumplir las condiciones mencionadas anteriormente: la configuración óptica, dimensiones y orientación, deben ser las mismas; y es necesario que hayan sido alineadas tomando como referencia el limbo.

Al igual que en el cálculo de áreas, el programa guarda varias focales para que solo tengamos que marcar la que hayamos utilizado en la observación. La fecha, la focal y el tamaño de los píxeles (en la pestaña Configuración), son los datos que nos permiten determinar la escala de la imagen y transformar píxeles en km.

Con Iris se mide la posición de un detalle en dos de las imágenes, y las coordenadas (X1, Y1), (X2, Y2), las introduciremos en el programa junto con el intervalo de tiempo transcurrido entre ambas imágenes (en horas).

Dado que desconocemos el ángulo entre la dirección de movimiento y la visual, no podemos saber cual es la velocidad real, pero al menos podemos acotarla. Si suponemos que el desplazamiento es perpendicular a la visual, el valor obtenido es una velocidad mínima.

Cuando pinchamos en el botón Calcular, el programa determina la velocidad suponiendo que el desplazamiento y la visual son perpendiculares. No obstante, tenemos la opción de ir variando el ángulo entre 5º y 90º y obtener la velocidad en cada caso. Lógicamente, para un ángulo de 0º la velocidad sería infinita.

 

 

Otros métodos

En esta pestaña están agrupados tres métodos de trabajo diferentes: la obtención de medidas con retículo, la creación de un archivo de posiciones de forma manual, y la utilización de Sol para controlar las operaciones con Iris.

 

 

Retículo:

El método para determinar posiciones con retículo se explica en http://www.parhelio.com/docposic.html

Hay dos formas de realizar los cronometrajes: con tiempos “absolutos” y “relativos”. Obtenemos tiempos absolutos cuando cada mancha por separado, la medimos respecto al limbo. Por el contrario, obtenemos tiempos relativos, cuando la primera mancha la medimos respecto al limbo, la segunda respecto a la primera, la tercera respecto a la segunda, y así sucesivamente. Aunque hay muchas opciones posibles, el programa solo funcionará correctamente si la sucesión de cronometrajes es la indicada. Las posiciones serían erróneas si, por ejemplo, la tercera mancha la medimos respecto a la primera.

Una vez introducidos la fecha y hora, los ángulos alfa y -alfa, y el limbo de referencia, hay que indicar el tipo de cronometraje que hemos empleado. En función de ello, la forma de introducir los tiempos es diferente:

Tiempos absolutos.-

Cada mancha se calcula por separado. Se introducen los tiempos con cada ángulo, y dando al botón “Calcular”, debajo aparecen los tiempos medios, las posiciones con ambos ángulos, y la posición media.

“Añadir medida” va incorporando cada medida al listado y, cuando acabemos, con el botón “Guardar”, guardamos un archivo de texto llamado “medidas.dat”, con todas ellas.

Tiempos relativos.-

Todas las manchas se calculan a la vez. En este caso, solo se introduce un tiempo por cada mancha. El primer tiempo es el de la mancha que se mide respecto al limbo, y debajo se escriben los tiempos relativos. Dando al botón “Calcular”, el programa primero suma la secuencia de tiempos obteniendo los valores absolutos para cada mancha, y con ellos calcula las posiciones.

Tanto en un caso como en otro, a partir de aquí el proceso es el mismo que con las medidas hechas a partir de imágenes.  En la pestaña “Posiciones” cargamos “medidas.dat”, asignamos cada posición a cada mancha, y creamos el archivo con el formato usado en Parhelio.

 

Archivo manual:

El archivo de posiciones creado por Sol, tiene el formato adecuado para incorporar las medidas a Parhelio. Si no hemos utilizado Sol para realizar dichas medidas, podemos crear manualmente un archivo con el formato correcto. Para ello, ponemos a la izquierda, la fecha y hora, y a la derecha vamos introduciendo los números de las manchas y sus coordenadas.

 

Movimientos

La pestaña “Movimientos” permite calcular el movimiento propio de una mancha a partir de sus medidas de posición. Este apartado analiza las medidas de un único observador, y no está pensado para combinar medidas de varios observadores. No obstante, si previamente se promedian las posiciones de varios observadores, puede utilizarse el listado resultante como punto de partida para obtener el movimiento propio con Sol.

En principio , es necesario disponer de un archivo de texto que incluya todas las medidas de la mancha, como las que se encuentran en http://www.parhelio.com/posic.html  También podemos crearlo copiando y pegando en un archivo de texto todas las medidas que hayamos realizado de una mancha, o las que hayamos hecho durante el tránsito de la mancha o durante una rotación, por ejemplo. La única condición es conservar el formato empleado por Sol.

 

 

El botón Archivo abre el archivo de texto con las medidas, y podemos filtrar las entradas por nombre de mancha y por código de observador. En caso de que alguna medida sea muy discrepante o contenga algún error, puede eliminarse mediante el botón Eliminar entrada.

Una vez que tengamos las posiciones correspondientes a una única mancha y un solo observador, realizamos los cálculos mediante el botón Movimiento propio. El resultado aparece en la ventana inferior.

Las columnas corresponden a la fecha, hora, longitud y latitud. DLong y DLat son el movimiento propio en longitud y latitud, tomando como origen la primera medida. Estas cantidades se calculan como se explica en http://www.parhelio.com/mp.html

DLongS y DLatS corresponden a la trayectoria suavizada, calculada mediante una media móvil, promediando 3 medidas consecutivas.

“Movimiento propio” también hace aparecer las siguientes cantidades en la pequeña ventana, situada abajo a la derecha:

N.- Número de días entre la primera y última posición.
D.- Distancia en grados recorrida a lo largo de la trayectoria.
V.- Velocidad media calculada a lo largo de la trayectoria.
W.- Velocidad angular.
T.- Periodo de rotación.

Si seleccionamos en la ventana de la izquierda otras dos fechas diferentes, el botón Cálculos determina esas cantidades entre las nuevas fechas. Cálculos suavizados hace lo mismo, pero con las posiciones de la trayectoria suavizada.

El botón Gráficas abre una nueva ventana donde podemos representar las posiciones medidas, el movimiento propio y la trayectoria suavizada. Este botón solo funcionará si previamente hemos calculado el movimiento propio y disponemos de un mínimo de 3 medidas.

 

 

 

Efemérides

Esta pestaña permite el cálculo de efemérides con independencia del resto de funciones.

 

 

Una vez introducidos el día y la hora, el botón “Calcular” nos muestra las efemérides para ese instante:

DJ.- Día juliano

D.- Diámetro solar en segundos de arco.

P.- Inclinación del eje de rotación solar, respecto a la dirección N-S.

B.- Latitud heliográfica del centro del disco solar.

L.- Longitud heliográfica del meridiano central.

Rot.- Número de rotación.

En la parte inferior, si introducimos el número de una rotación, obtenemos la fecha de inicio de la misma.

Los controles de la derecha permiten obtener listados con un número de entradas determinado. “Intervalo” es la diferencia en días, horas, o minutos, entre dos entradas consecutivas. Con el botón “Crear” obtenemos los datos.

El listado está formado por 6 columnas: Fecha, hora, P B, L y rotación. Si el número de rotación es negativo, no aparecerá en la lista.

Con el botón “Rotaciones” calculamos la fecha y hora de inicio de tantas rotaciones como indiquemos en el “Nº de entradas”, obtenidas a partir de la rotación introducida abajo a la izquierda.

 

 

Medida de posiciones paso a paso

- Iris.- Se convierten las imágenes a formato fit b&n, con el menú "Digital photo" ("Photo numérique"), "Decode RAW files..." ("Décodage des fichiers RAW...").

- Sol.- Botón Alinear.

- Sol.- Botón Sumar.

- Sol.- En la pestaña Orientación introducimos Día/Hora y obtenemos el ángulo alfa.

- Sol.- Botón Orientar para terminar de preparar la imagen.

- Iris.- Se calcula el centro fotométrico de una mancha seleccionándola en un recuadro, pinchando con el botón derecho y después en "Statistics" ("Statistiques").

- Sol.- Sin cerrar la ventana "Statistics" ("Statistiques"), en la pestaña Posiciones, con el botón Obtener posiciones calculamos las coordenadas heliográficas.

- Iris / Sol.- Se repiten los dos últimos pasos con todas las manchas.

- Sol.- En la segunda ventana se escribe la identificación de cada mancha.

- Sol.- Con el botón Crear archivo, ponemos las medidas con el formato adecuado.

- Sol.- Guardamos el resultado con el botón Guardar.