Suddenly I cannot read Shake data after day 100

Technical Support
1

using ReadMSEEDFast. I’ve been using my code for three years now with no read problems. Now, starting 2020 day 100 I get this read error. No problem reading 2019 day 100 or any day in 2020 before 100. Was there a software change resulting in a change in write format to seed files?

Index in position 1 exceeds array bounds.

Error in ReadMSEEDFast>DecodeSignal (line 314)
            signalMatrix=reshape(signalMatrix,numel(encodedSignalMatrix(1,:))/4,[]);

Error in ReadMSEEDFast (line 179)
                [signalMatrix,timeVector]=DecodeSignal(encodedSignalMatrix,ENCODING,isLittleEndian,BLOCK_RECORD_SIZE,headerInfoAllVolumes.startTime(sameChannelRows,:),headerInfoAllVolumes.sampleRate(sameChannelRows,:));
2

hi jay,

i don’t know what ReadmSEEDFast is so i can’t really comment on that code itself. (on a side note, i do find that function name interesting: is there a slow version?)

and i can confirm that there’s been no change to how mini-SEED format files are written, this is code that is probably more than 20 years old.

i would suggest contacting the person responsible for the function you refer to and asking them if they can provide any further insight into why it’s failing.

cheers,

richard

3

Standard Matlab code for reading seed, been using it for many years.

4

I use ReadMSEEDFast on a daily basis and I’m having no problems.

1 Like
5

Maybe you can try to read a sample file from my Shake.AM.RAC22.00.SHZ.D.2020.101 (6.0 MB)

6

Yes, I downloaded the mseed file from FDSN, then:

w=waveform(outfile,‘seed’);
ReadMSEEDFast: d:\RPi\RaspberryShake\FDSN\temp.mseed

Worked instantly.

1 Like