MEDICIONES SUIZAS DE RADIOACTIVIDAD

[maria urizar]

Mediciones de radioactividad de las estaciones suizas: en tiempo real. Se actualizan automaticamente

https://www.naz.ch/en/index.html

[py]

Lo primero gracias por mantener éstos hilos.

Lo segundo, cuando tenga un ratillo me comprometo a hacer un programilla para sacar la información de la página web y que genere un fichero .csv, así se puede almacenar y tratar luego como se quiera.

[maria urizar]

PY te agradezco mucho porque el incoveniente de estos graficos que se actualizan automaticamente es  que la información desaparece.

por ejemplo hace unas semanas, Tecnico preocupado pudo congelar estas graficas que indicaban la presencia de un alto nivel de radionucleidos:



[/quote]


Convendría que pudieramos guardar la información para comparar, porque lo que esta claro es que los organismos oficiales no nos van a informar, protegen al lobby nuclear

[py]

Allá va la primera versión, con los siguientes detalles:
- Hecha en linux y con python.
- Me descargo las páginas web (con wget).
- Limpio de polvo y paja los datos (a pelo).
- guardo todo en un fichero de texto (que ocupa 122kb, con lo que son +- 45mb/año, si no me equivoco), identificando de que estación son los datos con los tres caracteres que las identifican.
- Creo que he puesto todas las estaciones, que son 66.
- Tengo en mente a la gente que usa windows (con wget para windows: http://users.ugent.be/~bpuype/cgi-bin/fetch.pl?dl=wget/wget.exe) y python 2.7 (http://python.org/ftp/python/2.7/python-2.7.msi), aunque todavía está por adaptar para que funcione bien.

El resultado es un fichero de texto con la siguiente forma (importable y manejable facilmente):

Código: [Seleccionar]

TAE 22.06.2013 02:00 106 0 -
TAE 22.06.2013 03:00 115 2 -
TAE 22.06.2013 04:00 113 1 -
TAE 22.06.2013 05:00 111 1 -
TAE 22.06.2013 06:00 108 0 -
TAE 22.06.2013 07:00 102 0 -
TAE 22.06.2013 08:00 100 0 -
TAE 22.06.2013 09:00 98 0 -
TAE 22.06.2013 10:00 98 0 -
TAE 22.06.2013 11:00 98 0 -
TAE 22.06.2013 12:00 97 0 -
TAE 22.06.2013 13:00 101 0 -
ABO 20.06.2013 00:00 116 0 -
ABO 20.06.2013 01:00 116 0 -
ABO 20.06.2013 02:00 118 0 -
ABO 20.06.2013 03:00 119 0 -
ABO 20.06.2013 04:00 117 0 -
ABO 20.06.2013 05:00 119 0 -
ABO 20.06.2013 06:00 129 1 -
ABO 20.06.2013 07:00 125 0 -


[py]

El código es el siguiente (presupongo que donde se copie el fichero Estaciones_Suizas.py tiene dos subdirectorios: paginas y datos)

Si se fastidian muchos las tabulaciones, subo el fichero donde queráis.
Los nómbres de las estaciones, como los he sacado del html de la página, ni los he "limpiado".

Código: [Seleccionar]

#!/usr/bin/env python
# -*- coding: latin -*-

import os

class Estaciones:
def __init__(self):

self.URL="https://www.naz.ch/en/aktuell/"
self.Estaciones_Suizas=[("zeitverlauf_TAE.shtml","Aadorf / Tänikon TG"),
("zeitverlauf_ABO.shtml","Adelboden BE"),
("zeitverlauf_AIG.shtml","Aigle VD"),
("zeitverlauf_ALT.shtml","Altdorf UR"),
("zeitverlauf_BAS.shtml","Basel / Binningen BL"),
("zeitverlauf_BER.shtml","Bern / Zollikofen BE"),
("zeitverlauf_BEZ.shtml","Beznau AG"),
("zeitverlauf_BUS.shtml","Buchs / Aarau AG"),
("zeitverlauf_FRE.shtml","Bullet / La Frêtaz VD"),
("zeitverlauf_CHA.shtml","Chasseral BE"),
("zeitverlauf_CHU.shtml","Chur GR"),
("zeitverlauf_GSB.shtml","Col du Grand St-Bernhard VS"),
("zeitverlauf_DAV.shtml","Davos GR"),
("zeitverlauf_DIS.shtml","Disentis / Sedrun GR"),
("zeitverlauf_EGO.shtml","Egolzwil LU"),
("zeitverlauf_ENG.shtml","Engelberg OW"),
("zeitverlauf_FAH.shtml","Fahy / Boncourt JU"),
("zeitverlauf_GVE.shtml","Genève-Cointrin GE"),
("zeitverlauf_GLA.shtml","Glarus GL"),
("zeitverlauf_GOE.shtml","Gösgen SO"),
("zeitverlauf_GRE.shtml","Grenchen SO"),
("zeitverlauf_GUE.shtml","Gütsch ob Andermatt UR"),
("zeitverlauf_GUT.shtml","Güttingen TG"),
("zeitverlauf_HOE.shtml","Hörnli ZH"),
("zeitverlauf_JUN.shtml","Jungfraujoch VS"),
("zeitverlauf_CDF.shtml","La Chaux-de-Fonds NE"),
("zeitverlauf_DOL.shtml","La Dôle VD"),
("zeitverlauf_MLS.shtml","Le Moléson FR"),
("zeitverlauf_LEI.shtml","Leibstadt AG"),
("zeitverlauf_OTL.shtml" ,"Locarno / Monti TI"),
        ("zeitverlauf_LUG.shtml" ,"Lugano TI"),
        ("zeitverlauf_LUZ.shtml" ,"Luzern LU"),
        ("zeitverlauf_MAG.shtml" ,"Magadino / Cadenazzo TI"),
        ("zeitverlauf_MER.shtml" ,"Meiringen BE"),
        ("zeitverlauf_MOE.shtml" ,"Möhlin AG"),
        ("zeitverlauf_MUB.shtml" ,"Mühleberg BE"),
        ("zeitverlauf_NAP.shtml" ,"Napf BE"),
        ("zeitverlauf_NEU.shtml" ,"Neuchâtel NE"),
        ("zeitverlauf_CGI.shtml" ,"Nyon / Changins VD"),
        ("zeitverlauf_ORO.shtml" ,"Oron VD"),
        ("zeitverlauf_PAY.shtml" ,"Payerne VD"),
        ("zeitverlauf_PIO.shtml" ,"Piotta TI"),
        ("zeitverlauf_PLF.shtml" ,"Plaffeien FR"),
        ("zeitverlauf_ROB.shtml" ,"Poschiavo / Robbia GR"),
        ("zeitverlauf_PUY.shtml" ,"Pully VD"),
        ("zeitverlauf_RUE.shtml" ,"Rünenberg BL"),
        ("zeitverlauf_SAM.shtml" ,"Samedan GR"),
        ("zeitverlauf_SBE.shtml" ,"San Bernardino GR"),
        ("zeitverlauf_SMM.shtml" ,"Santa Maria GR"),
        ("zeitverlauf_SAE.shtml" ,"Säntis AI"),
        ("zeitverlauf_SHA.shtml" ,"Schaffhausen SH"),
        ("zeitverlauf_SCU.shtml" ,"Scuol GR"),
        ("zeitverlauf_SIO.shtml" ,"Sion VS"),
        ("zeitverlauf_SBO.shtml" ,"Stabio TI"),
        ("zeitverlauf_STG.shtml" ,"St. Gallen SG"),
        ("zeitverlauf_ULR.shtml" ,"Ulrichen VS"),
        ("zeitverlauf_VAD.shtml" ,"Vaduz"),
        ("zeitverlauf_VIO.shtml" ,"Vicosoprano GR"),
        ("zeitverlauf_PSI.shtml" ,"Villigen, PSI AG "),
        ("zeitverlauf_VIS.shtml" ,"Visp VS"),
        ("zeitverlauf_WAE.shtml" ,"Wädenswil ZH"),
        ("zeitverlauf_WYN.shtml" ,"Wynau BE"),
        ("zeitverlauf_ZER.shtml" ,"Zermatt VS"),
        ("zeitverlauf_REH.shtml" ,"Zürich / Affoltern ZH"),
        ("zeitverlauf_SMA.shtml" ,"Zürich / Fluntern ZH"),
        ("zeitverlauf_KLO.shtml" ,"Zürich-Kloten ZH")
]

self.Destino_Paginas="paginas"
self.Destino_Datos="datos"

def Procesa(self):
#!!!!!!!!! OJO QUE ME CARGO EL .SHTML !!!!!!!!!!!
comando="rm "+self.Destino_Paginas+"/*"
os.system(comando)

total=len(self.Estaciones_Suizas)
for i in range(0,len(self.Estaciones_Suizas)):
#Obtención de la página web sin complicaciones con wget.
URL=self.URL+self.Estaciones_Suizas[i][0]
print "(%d de %d) Descargando %s" % (i,total,self.Estaciones_Suizas[i][0])
os.system("wget -q "+URL+" -O "+self.Destino_Paginas+"/"+self.Estaciones_Suizas[i][0])
print "Descarga finalizada."

def Extraer_Datos(self):
#A lo bruto...

print "Procesando descargas..."
for j in range(0,len(self.Estaciones_Suizas)):
estacion=self.Estaciones_Suizas[j][0][-9:-6]
#print estacion

#Extraer la tabla...
A_ello=0
Cadena=""
fichero=open(self.Destino_Paginas+"/"+self.Estaciones_Suizas[j][0])
for linea in fichero.readlines():
if A_ello:
Cadena+=linea[:-1]

#if linea.find("<table cellspacing='0'>")>0:
if linea.find("Bemerkungen")>0:
A_ello=1


if linea.find("</table>")>0:
A_ello=0

#print Cadena
fichero.close()

#Extraer datos de xml...
Cadena2=""
Borrando=0
for i in range(0,len(Cadena)):
if Cadena[i]=="<":
Borrando=1
if not Borrando:
if ord(Cadena[i])>31 and ord(Cadena[i])<123:
Cadena2+=Cadena[i]
if Cadena[i]==">":
Borrando=0

#print Cadena2

#Limpiar un poco...
Cadena3=""
for i in range(0,len(Cadena2)-1):
if Cadena2[i]==" " and Cadena2[i+1]<>" ":
Cadena3+=Cadena2[i]
if Cadena2[i]<>" ":
Cadena3+=Cadena2[i]

#Palabras...
Tabla=[]
palabra=""
en_palabra=0
for i in range(0,len(Cadena3)):
if Cadena3[i]<>" ":
en_palabra=1
else:
if en_palabra:
en_palabra=0
Tabla.append(palabra)
palabra=""

if en_palabra:
palabra+=Cadena3[i]
if en_palabra:
Tabla.append(palabra)

#Probemos primero con un fichero.
fichero=open(self.Destino_Datos+"/prueba.txt","a")
#print Tabla,len(Tabla)
linea=""
for i in range(0,len(Tabla)/5):
linea+=estacion+" "
for ii in range(0,5):
linea+=Tabla[i*5+ii]+" "
linea+="\n"
fichero.write(linea)
linea=""
fichero.close()

print "Descargas procesadas."

if __name__=="__main__":
E=Estaciones()
E.Procesa()
E.Extraer_Datos()



[maria urizar]

Muchas gracias PY.

Yo lo intente de un modo casero en burbuja pero no se podían manejar lso datos.

Estos datos que vayamos recogiendo van a tener un valor incomensurable; porque los niveles de radioactividad en Japón están subiendo.

Visto lo que a ocurrido con Chernobil y Fukushima, queda claro que ningún gobierno ni organismo nos va a informar si las cosas van mal.

Tenemos que informarnos y sacar nuestras propias conclusiones nosostros, y a la contra de toda la propaganda pronuclear

Excelente trabajo PY

[py]

Modificado un poco, porque el tener el fichero es poco práctico, así que he añadido que lo guarde automáticamente en una base de datos sqlite3, que es más manejable.

Para crear la base de datos (dentro del directorio datos):
sqlite3 radioactividad.db
Y creamos la tabla (índice puesto para que no se repitan los registros si descargamos más de una vez lo mismo):
create table mediciones (estacion varchar(10), fecha varchar(10), hora varchar(10), valor1 varchar(10), valor2 varchar(10), valor3 varchar(10), primary key (estacion, fecha, hora));

Código:

Código: [Seleccionar]

#!/usr/bin/env python
# -*- coding: latin -*-

import os
import sqlite3

class Estaciones:
def __init__(self):

self.URL="https://www.naz.ch/en/aktuell/"
self.Estaciones_Suizas=[("zeitverlauf_TAE.shtml","Aadorf / Tänikon TG"),
("zeitverlauf_ABO.shtml","Adelboden BE"),
("zeitverlauf_AIG.shtml","Aigle VD"),
("zeitverlauf_ALT.shtml","Altdorf UR"),
("zeitverlauf_BAS.shtml","Basel / Binningen BL"),
("zeitverlauf_BER.shtml","Bern / Zollikofen BE"),
("zeitverlauf_BEZ.shtml","Beznau AG"),
("zeitverlauf_BUS.shtml","Buchs / Aarau AG"),
("zeitverlauf_FRE.shtml","Bullet / La Frêtaz VD"),
("zeitverlauf_CHA.shtml","Chasseral BE"),
("zeitverlauf_CHU.shtml","Chur GR"),
("zeitverlauf_GSB.shtml","Col du Grand St-Bernhard VS"),
("zeitverlauf_DAV.shtml","Davos GR"),
("zeitverlauf_DIS.shtml","Disentis / Sedrun GR"),
("zeitverlauf_EGO.shtml","Egolzwil LU"),
("zeitverlauf_ENG.shtml","Engelberg OW"),
("zeitverlauf_FAH.shtml","Fahy / Boncourt JU"),
("zeitverlauf_GVE.shtml","Genève-Cointrin GE"),
("zeitverlauf_GLA.shtml","Glarus GL"),
("zeitverlauf_GOE.shtml","Gösgen SO"),
("zeitverlauf_GRE.shtml","Grenchen SO"),
("zeitverlauf_GUE.shtml","Gütsch ob Andermatt UR"),
("zeitverlauf_GUT.shtml","Güttingen TG"),
("zeitverlauf_HOE.shtml","Hörnli ZH"),
("zeitverlauf_JUN.shtml","Jungfraujoch VS"),
("zeitverlauf_CDF.shtml","La Chaux-de-Fonds NE"),
("zeitverlauf_DOL.shtml","La Dôle VD"),
("zeitverlauf_MLS.shtml","Le Moléson FR"),
("zeitverlauf_LEI.shtml","Leibstadt AG"),
("zeitverlauf_OTL.shtml" ,"Locarno / Monti TI"),
        ("zeitverlauf_LUG.shtml" ,"Lugano TI"),
        ("zeitverlauf_LUZ.shtml" ,"Luzern LU"),
        ("zeitverlauf_MAG.shtml" ,"Magadino / Cadenazzo TI"),
        ("zeitverlauf_MER.shtml" ,"Meiringen BE"),
        ("zeitverlauf_MOE.shtml" ,"Möhlin AG"),
        ("zeitverlauf_MUB.shtml" ,"Mühleberg BE"),
        ("zeitverlauf_NAP.shtml" ,"Napf BE"),
        ("zeitverlauf_NEU.shtml" ,"Neuchâtel NE"),
        ("zeitverlauf_CGI.shtml" ,"Nyon / Changins VD"),
        ("zeitverlauf_ORO.shtml" ,"Oron VD"),
        ("zeitverlauf_PAY.shtml" ,"Payerne VD"),
        ("zeitverlauf_PIO.shtml" ,"Piotta TI"),
        ("zeitverlauf_PLF.shtml" ,"Plaffeien FR"),
        ("zeitverlauf_ROB.shtml" ,"Poschiavo / Robbia GR"),
        ("zeitverlauf_PUY.shtml" ,"Pully VD"),
        ("zeitverlauf_RUE.shtml" ,"Rünenberg BL"),
        ("zeitverlauf_SAM.shtml" ,"Samedan GR"),
        ("zeitverlauf_SBE.shtml" ,"San Bernardino GR"),
        ("zeitverlauf_SMM.shtml" ,"Santa Maria GR"),
        ("zeitverlauf_SAE.shtml" ,"Säntis AI"),
        ("zeitverlauf_SHA.shtml" ,"Schaffhausen SH"),
        ("zeitverlauf_SCU.shtml" ,"Scuol GR"),
        ("zeitverlauf_SIO.shtml" ,"Sion VS"),
        ("zeitverlauf_SBO.shtml" ,"Stabio TI"),
        ("zeitverlauf_STG.shtml" ,"St. Gallen SG"),
        ("zeitverlauf_ULR.shtml" ,"Ulrichen VS"),
        ("zeitverlauf_VAD.shtml" ,"Vaduz"),
        ("zeitverlauf_VIO.shtml" ,"Vicosoprano GR"),
        ("zeitverlauf_PSI.shtml" ,"Villigen, PSI AG "),
        ("zeitverlauf_VIS.shtml" ,"Visp VS"),
        ("zeitverlauf_WAE.shtml" ,"Wädenswil ZH"),
        ("zeitverlauf_WYN.shtml" ,"Wynau BE"),
        ("zeitverlauf_ZER.shtml" ,"Zermatt VS"),
        ("zeitverlauf_REH.shtml" ,"Zürich / Affoltern ZH"),
        ("zeitverlauf_SMA.shtml" ,"Zürich / Fluntern ZH"),
        ("zeitverlauf_KLO.shtml" ,"Zürich-Kloten ZH")
]

self.Destino_Paginas="paginas"
self.Destino_Datos="datos"

def Procesa(self):
#!!!!!!!!! OJO QUE ME CARGO EL .SHTML !!!!!!!!!!!
comando="rm "+self.Destino_Paginas+"/*"
os.system(comando)

total=len(self.Estaciones_Suizas)
for i in range(0,len(self.Estaciones_Suizas)):
#Obtención de la página web sin complicaciones con wget.
URL=self.URL+self.Estaciones_Suizas[i][0]
print "(%d de %d) Descargando %s" % (i,total,self.Estaciones_Suizas[i][0])
os.system("wget -q "+URL+" -O "+self.Destino_Paginas+"/"+self.Estaciones_Suizas[i][0])
print "Descarga finalizada."

def Extraer_Datos(self):
#A lo bruto...

print "Procesando descargas..."
for j in range(0,len(self.Estaciones_Suizas)):
estacion=self.Estaciones_Suizas[j][0][-9:-6]
#print estacion

#Extraer la tabla...
A_ello=0
Cadena=""
fichero=open(self.Destino_Paginas+"/"+self.Estaciones_Suizas[j][0])
for linea in fichero.readlines():
if A_ello:
Cadena+=linea[:-1]

#if linea.find("<table cellspacing='0'>")>0:
if linea.find("Bemerkungen")>0:
A_ello=1


if linea.find("</table>")>0:
A_ello=0

#print Cadena
fichero.close()

#Extraer datos de xml...
Cadena2=""
Borrando=0
for i in range(0,len(Cadena)):
if Cadena[i]=="<":
Borrando=1
if not Borrando:
if ord(Cadena[i])>31 and ord(Cadena[i])<123:
Cadena2+=Cadena[i]
if Cadena[i]==">":
Borrando=0

#print Cadena2

#Limpiar un poco...
Cadena3=""
for i in range(0,len(Cadena2)-1):
if Cadena2[i]==" " and Cadena2[i+1]<>" ":
Cadena3+=Cadena2[i]
if Cadena2[i]<>" ":
Cadena3+=Cadena2[i]

#Palabras...
Tabla=[]
palabra=""
en_palabra=0
for i in range(0,len(Cadena3)):
if Cadena3[i]<>" ":
en_palabra=1
else:
if en_palabra:
en_palabra=0
Tabla.append(palabra)
palabra=""

if en_palabra:
palabra+=Cadena3[i]
if en_palabra:
Tabla.append(palabra)

#Probemos primero con un fichero.
fichero=open(self.Destino_Datos+"/prueba.txt","a")
#print Tabla,len(Tabla)
linea=""
for i in range(0,len(Tabla)/5):
linea+=estacion+" "
for ii in range(0,5):
linea+=Tabla[i*5+ii]+" "
linea+="\n"
fichero.write(linea)
linea=""
fichero.close()

#Ahora, guardamos en SQlite3.
#Nombre bbdd: radioactividad.db
#Definicion de la tabla:
#create table mediciones (estacion varchar(10), fecha varchar(10), hora varchar(10), valor1 varchar(10), valor2 varchar(10), valor3 varchar(10), primary key (estacion, fecha, hora)); 
conexion=sqlite3.connect(self.Destino_Datos+"/"+"radioactividad.db")
c=conexion.cursor()
for i in range(0,len(Tabla)/5):
sql="insert into mediciones values ('"+estacion+"','"+Tabla[i*5]+"','"+Tabla[i*5+1]+"','"+Tabla[i*5+2]+"','"+Tabla[i*5+3]+"','"+Tabla[i*5+4]+"')"
try:
c.execute(sql)
except:
#print "Datos ya existentes: ", sql
pass
conexion.commit()
c.close()

print "Descargas procesadas."




if __name__=="__main__":
E=Estaciones()
E.Procesa()
E.Extraer_Datos()



[maria urizar]

Py, cogí mediciones en estos 2 hilos, desde el 2 de noviembre en unas semanas para tenerlas como rteferencias, porque la contaminación radioactiva ira subiendo y la propaganda intentara hacer creer qeu es "normal", sólo si tenemos valores anteriores podremos comparar:

http://www.burbuja.info/inmobiliaria/temas-calientes/359887-mediciones-suizas-de-radioactividad.html

http://www.burbuja.info/inmobiliaria/temas-calientes/359890-seguimiento-valores-radioactividad-ciudades-suizas.html

[maria urizar]

Py hay picos significativos en als mediciones puedes registralarlas??

[py]

Py hay picos significativos en als mediciones puedes registralarlas??

Las estoy registrando.
El gráfico de ejemplo muestra todas las mediciones de todas las estaciones por cada día, da una idea de los rangos por los que suele oscilar, vean por ejemplo el día 08/07/13.

[maria urizar]

El día 3 del 7 se registraron estos picos.
Estos picos tan altos que estamos registrando no los había visto hasta ahora

[py]

Actualización:
-Puntos amarillos corresponden a mediciones individuales de todas las estaciones.
-Cuadrados azules es la media diaria de todas las estaciones.
-Línea azul es la media del mes de todas las estaciones.

Los días que faltan no los tengo registrados.

[maria urizar]

ATENCION:



Los valores medidos de la estación Nadam en Leibstadt AG

Valor promedio por hora desde 07/10/2013 hasta 09/10/2013 El área de distribución natural de los valores de la estación de Leibstadt AG es 75 nSv / h y 195 nSv / h.

08.10.2013 14:00   93   0    -
08.10.2013 15:00   126   2    -
08.10.2013 16:00   194   5    -
08.10.2013 17:00   244   2   
08.10.2013 18:00   161   0    -
08.10.2013 19:00   111   0    -
08.10.2013 20:00   97   0    -
08.10.2013 21:00   98   0    -
08.10.2013 22:00   111   1    -
08.10.2013 23:00   117   0    -
09/10/2013 00:00   101   0    -
10.09.2013 01:00   116   1    -
09/10/2013 02:00   104   0    -
09/10/2013 03:00   119   3    -
09/10/2013 04:00   158   2    -
09/10/2013 05:00   189   2    -
09/10/2013 06:00   356   0   
09/10/2013 07:00   442   -   
09/10/2013 08:00   274   -   
09/10/2013 09:00   140   -    -





Como en  la mayor parte de estaciones no llueve no hay picos en las mediciones de radiactividad, en las pocos estaciones que llueve, con muy poca cantidad de precipitación los picos se elevan mucho.

Esto puede interpretarse como que hay una gran concentración de radionucleidos

Poca precipitación y grandes picos = gran concentracion de radionucleidos[COLOR="Silver"]

[Datum]

Py ; Si  hiciera falta algo en ansi C comentamelo.

[maria urizar]

¿Que es eso?
En este hilo haría falta ayuda, alguien que sepa congelar las graficas y colgarlas cuando estan altas;