Hola, yo escribi un script (archivo adjunto) para calcular los pixeles en datos
de imagenes de satelite. Cuando lo corro obtengo lo siguiente:
Reading in: 200406010145.goes12ir
For file 200406010145.goes12ir count = 3126
Reading in: 200406010545.goes12ir
For file 200406010545.goes12ir count = 4455
Reading in: 200406012115.goes12ir
For file 200406012115.goes12ir count = 4690
Reading in: 200406010915.goes12ir
For file 200406010915.goes12ir count = 3373
Reading in: 200406010315.goes12ir
For file 200406010315.goes12ir count = 2390
Reading in: 200406011945.goes12ir
For file 200406011945.goes12ir count = 5902
Reading in: 200406010715.goes12ir
For file 200406010715.goes12ir count = 4644
Reading in: 200406011715.goes12ir
For file 200406011715.goes12ir count = 5613
Reading in: 200406012145.goes12ir
For file 200406012145.goes12ir count = 3976
Reading in: 200406010845.goes12ir
For file 200406010845.goes12ir count = 5269
Reading in: 200406011345.goes12ir
Como le puedo hacer para guardar estas salidas del script en un archivo de
texto.
Boris Vladimir Comi Gonzalez
Universidad Nacional Autónoma de México
Grupo de Tormentas Convecivas
# Call data from the directory(/home/mcidasv/Documentos/SSEC/20040601)
def scm(directory):
import os;
fs = os.listdir(directory);
# now fs will be a list of all the files in directory
from edu.wisc.ssec.mcidas import AreaFile;
for name in fs:
print "Reading in:",name
af = AreaFile(directory+"/"+name);
ad = af.getAreaDirectory();
count = 0;
data = af.getFloatData();
# now look through the first band y count pixels
# MCS detected when his temperature infrared (TIR) is < 219 K
for i in xrange(ad.getLines()):
for j in xrange(ad.getElements()):
if (data[0][i][j]) == 199:
count = count + 1;
print "For file",name," count = ",count
scm("/home/mcidasv/Documentos/SSEC/20040601")
_______________________________________________
Python-es mailing list
[email protected]
http://mail.python.org/mailman/listinfo/python-es
FAQ: http://python-es-faq.wikidot.com/
