Hello !<br><br>I am trying to iterate through the pixels  of  an FITS image read previusly with PyFits,<br>but it is too slow (more that 1CPU minute), so my question is if is there any other way to do that in a <br>more efficient/fast  way. <br>
<br>That&#39;s the code I&#39;m using:<br>------------------------------------------------------------------<br>f=pyfits.open(&#39;/tmp/c1.fits&#39;,memmap=1)<br>b=numpy.zeros([2048,2048],dtype=&#39;float32&#39;)<br>cdata=f[0].data.copy()<br>
f.close()<br>print &#39;start loop...&#39;<br>for i in range(0,2048):<br>    for j in range(0,2048):<br>        if cdata[i,j]&lt;1 or cdata[i,j]&gt;100000:<br>            b[i,j]=1<br><br><br>----------------------------------------------------------------<br>
<br>Any idea ?<br><br>Thanks very much,<br>JM<br>