0

Matlabに次のコードスニペットがあり、Pythonに移植したいと思います。

for i=1:runs;
   tic;
   datmat=importdata('data.txt',' ');
    Limits_Area=[1,261,522,784,1045,1305,1565,1827,2088,2349,2610,2871,3131,3158];
   for k=1:2000;
      for j=1:13;
         Year_Data=datmat.data(Limits_Area(j):Limits_Area(j+1)-1,4:37);
         Min_Year=min(Year_Data);
         Max_Year=max(Year_Data);
         Mean_Year=mean(Year_Data);
         GainLoss_Year=100-(Year_Data(1,:)+0.00000001)./(Year_Data(end,:)+0.00000001)*100;
     end;
   end;

私は本当に苦労しています

Year_Data=datmat.data(Limits_Area(j):Limits_Area(j+1)-1,4:37);

一部....何か方向性はありますか？

ありがとうございました

4

1 に答える 1

2

numpyを試しましたか？これは、Matlabのように見える科学計算用のライブラリです。例えば ​​：

Matlabで： a(1:5,:)

Numpyで： a[0:5] or a[:5] or a[0:5,:]

チェックアウト：Matlabユーザー向けのNumpy

Numpyを使用したくない場合は、理解リストを試してください。

Year_Data = [ [datmat.data(i,j) for j in range (4,38) ] for i in range(j,j+2) ]

編集： 

for i in range(runs) :
    datamat = numpy.genfromtxt('data.txt',delimiter=' ', newline ='\n' ) 
    // Adapt the previous line to the format of your txt file
    // at this point you should have a numpy.array object with the right shape
    Limits_Area= numpy.array( [1,261,522,784,1045,1305,1565,1827,2088,2349,2610,2871,3131,3158] )
    for k in range(2000):
        for j in range(13):
            Year_Data = datmat[ Limits_Area(j):Limits_Area(j+1)-1 , 4:37 ]
            etc etc ...

注意：Matlab配列インデックスは1からnになりますが、numpy配列インデックスは0からn-1になります

于 2012-10-22T14:15:42.977 に答える