转换浮动与NumPy阵列大端

我写了一个脚本，并在最后，我需要这个数组是在类型float64于大端INT（> 2I）转换：转换浮动与NumPy阵列大端

[[ 0.92702157 1.03092008 0.9072934 ..., 0.71617331 1.02524888 
    1.07284994] 
[ 0.99573712 0.96416766 0.9230931 ..., 0.66935196 0.64930711 
0.5357821 ] 
[ 0.98846306 1.03608056 0.79976885 ..., 0.69383804 0.62434976 
0.88219911] 
..., 
[ 0.91196013 0.87880101 0.97145563 ..., 0.79110817 1.19651477 
0.98244941] 
[ 1.0129829 0.81045263 0.95434107 ..., 0.99752385 1.08271169 
1.12872492] 
[ 0.94037117 0.81365084 0.94384051 ..., 0.82754351 1.03742172 1.]]

我怎么能这样做？

这里是整个脚本

import numpy as np 
import pyfits 
from matplotlib import pyplot as plt 
import glob 
import os 
import re 
from struct import unpack,pack 

global numbers 
numbers=re.compile(r'(\d+)') 


def numericalSort(value): 
    parts = numbers.split(value) 
    parts[1::2] = map(int, parts[1::2]) 
    return parts 

dark=sorted(glob.glob(' *.fits'),key=numericalSort) 
flat=sorted(glob.glob('/ *.fits'),key=numericalSort) 
img=sorted(glob.glob('/ *.fits'),key=numericalSort) 

#dark 
sumd0 = pyfits.open(dark[0]) 
sumdd=sumd0[0].data 
sumdd = sumdd.astype(float,copy=False) 

for i in range(1,len(dark)): 
    sumdi=pyfits.open(dark[i]) 
    sumdi=sumdi[0].data 
    sumdd=sumdd+sumdi.astype(float, copy=False) 


dd=sumdd/len(dark) 


#flat 
sumf0 = pyfits.open(flat[0]) 
sumff=sumf0[0].data 
sumff = sumff.astype(float, copy=False) 

for i in range(1,len(flat)): 
    sumfi=pyfits.open(flat[i]) 
    sumfi=sumfi[0].data 
    sumff=sumff+sumfi.astype(float,copy=False) 

ff=sumff/len(flat) 


df=(ff-dd) 
maxx=np.max(df) 
df=np.clip(df,1,maxx) 


for n in range(len(img)): 
    im=pyfits.open(img[n]) 
    imgg=im[0].data 
    header=im[0].header 
    imgg=imgg.astype(float,copy=False) 
    x,y=im[0].shape 
    m=np.max(imgg) 
    imgg=np.clip((imgg-dd),1,m) 
    imgg=imgg/df 
    imgg=np.clip(imgg,0.5,1.5) 


    #print imgg.dtype 
    #imgg=imgg[200:950,150:1250] 
    #imgg=imgg[::-1,:y] 
    hdu = pyfits.PrimaryHDU(imgg,header) 
    hdulist = pyfits.HDUList([hdu]) 
    hdulist.writeto('/c'+img[n][48:]) 
    plt.imshow(imgg,cmap=plt.cm.Greys_r) 
#plt.savefig('/+'.png') 
plt.show()

在过去的9号线，我需要从float64转换阵列imgg到> I2

来源

2015-07-21 Franta Konopnik

你可以试试'a.view （'> 2i'）'（其中'a'是你数组的名称） - 你将得到一个数组视图，其中元素的数量是两倍，并且是一个额外的维度。这是你想要的吗？ –

不......我需要将这个数组保存到数据文件中......它是类型> 2i中的图像...在开始时我在适合文件（2.8MiB）中有图像数据...我的脚本进行一些操作（+ - /）与此数据数组，但此操作需要类型浮法...但是，当我得到11 MiB的大小后保存此数组...我需要将此数组转换回> 2i –

>>> import numpy as np 
>>> big_end = bytes(chr(0) + chr(1) + chr(3) + chr(2), 'utf-8') 
>>> np.ndarray(shape=(2,),dtype='>i2', buffer=big_end) 
array([ 1, 770], dtype=int16)

还看到：Byte-swapping

来源

2015-07-21 18:18:59 Robert

我得到此错误：big_end = bytes （0）+ chr（1）+ chr（3）+ chr（2），'utf-8'） TypeError：str（）最多支持1个参数（给出2） –

您是否正在使用Python 2.X ？如果你是，你可以完全省略bytes（）。 – Robert

我使用Python 2.7和debian 7 64x –

转换浮动与NumPy阵列大端

回答

相关问题