JPWO2019229998A1 - 機器分析用データ処理方法及び装置 - Google Patents
機器分析用データ処理方法及び装置 Download PDFInfo
- Publication number
- JPWO2019229998A1 JPWO2019229998A1 JP2020522557A JP2020522557A JPWO2019229998A1 JP WO2019229998 A1 JPWO2019229998 A1 JP WO2019229998A1 JP 2020522557 A JP2020522557 A JP 2020522557A JP 2020522557 A JP2020522557 A JP 2020522557A JP WO2019229998 A1 JPWO2019229998 A1 JP WO2019229998A1
- Authority
- JP
- Japan
- Prior art keywords
- data
- value
- array
- intensity
- block
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000004458 analytical method Methods 0.000 title claims description 110
- 238000003672 processing method Methods 0.000 title claims description 63
- 238000007906 compression Methods 0.000 claims abstract description 151
- 230000006835 compression Effects 0.000 claims abstract description 149
- 238000001819 mass spectrum Methods 0.000 claims abstract description 91
- 238000000034 method Methods 0.000 claims description 231
- 238000012545 processing Methods 0.000 claims description 187
- 238000005259 measurement Methods 0.000 claims description 91
- 230000008569 process Effects 0.000 claims description 70
- 238000003491 array Methods 0.000 claims description 40
- 238000006243 chemical reaction Methods 0.000 claims description 34
- 230000003068 static effect Effects 0.000 claims description 15
- 238000004949 mass spectrometry Methods 0.000 claims 6
- 238000003384 imaging method Methods 0.000 description 48
- 238000013144 data compression Methods 0.000 description 42
- 239000011159 matrix material Substances 0.000 description 27
- 238000007792 addition Methods 0.000 description 22
- 238000001228 spectrum Methods 0.000 description 22
- 238000013500 data storage Methods 0.000 description 18
- 230000006837 decompression Effects 0.000 description 13
- 238000004364 calculation method Methods 0.000 description 12
- 238000001871 ion mobility spectroscopy Methods 0.000 description 12
- 238000000491 multivariate analysis Methods 0.000 description 9
- 238000010586 diagram Methods 0.000 description 8
- 238000003860 storage Methods 0.000 description 7
- 229910052770 Uranium Inorganic materials 0.000 description 6
- 229910052698 phosphorus Inorganic materials 0.000 description 6
- 239000000523 sample Substances 0.000 description 6
- 230000007423 decrease Effects 0.000 description 5
- 230000000694 effects Effects 0.000 description 5
- 238000001069 Raman spectroscopy Methods 0.000 description 4
- 229910052799 carbon Inorganic materials 0.000 description 4
- 238000013480 data collection Methods 0.000 description 4
- 238000000701 chemical imaging Methods 0.000 description 3
- 150000001875 compounds Chemical class 0.000 description 3
- 230000005484 gravity Effects 0.000 description 3
- 230000000670 limiting effect Effects 0.000 description 3
- 238000012935 Averaging Methods 0.000 description 2
- 238000000305 Fourier transform infrared microscopy Methods 0.000 description 2
- 230000008901 benefit Effects 0.000 description 2
- 238000009826 distribution Methods 0.000 description 2
- 238000010894 electron beam technology Methods 0.000 description 2
- 238000000605 extraction Methods 0.000 description 2
- 230000010354 integration Effects 0.000 description 2
- 238000002360 preparation method Methods 0.000 description 2
- 238000000513 principal component analysis Methods 0.000 description 2
- 230000002829 reductive effect Effects 0.000 description 2
- 230000002441 reversible effect Effects 0.000 description 2
- 238000002798 spectrophotometry method Methods 0.000 description 2
- 238000001269 time-of-flight mass spectrometry Methods 0.000 description 2
- 238000012546 transfer Methods 0.000 description 2
- 238000005033 Fourier transform infrared spectroscopy Methods 0.000 description 1
- 238000002083 X-ray spectrum Methods 0.000 description 1
- 239000012472 biological sample Substances 0.000 description 1
- 238000007405 data analysis Methods 0.000 description 1
- 230000002427 irreversible effect Effects 0.000 description 1
- 239000007788 liquid Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000000414 obstructive effect Effects 0.000 description 1
- 238000007781 pre-processing Methods 0.000 description 1
- 102000004169 proteins and genes Human genes 0.000 description 1
- 108090000623 proteins and genes Proteins 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 229920006395 saturated elastomer Polymers 0.000 description 1
- 238000004904 shortening Methods 0.000 description 1
- 230000003595 spectral effect Effects 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 229910052717 sulfur Inorganic materials 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
- 239000002699 waste material Substances 0.000 description 1
Images
Classifications
-
- H—ELECTRICITY
- H01—ELECTRIC ELEMENTS
- H01J—ELECTRIC DISCHARGE TUBES OR DISCHARGE LAMPS
- H01J49/00—Particle spectrometers or separator tubes
- H01J49/0004—Imaging particle spectrometry
-
- H—ELECTRICITY
- H01—ELECTRIC ELEMENTS
- H01J—ELECTRIC DISCHARGE TUBES OR DISCHARGE LAMPS
- H01J49/00—Particle spectrometers or separator tubes
- H01J49/0027—Methods for using particle spectrometers
- H01J49/0036—Step by step routines describing the handling of the data generated during a measurement
-
- H—ELECTRICITY
- H03—ELECTRONIC CIRCUITRY
- H03M—CODING; DECODING; CODE CONVERSION IN GENERAL
- H03M7/00—Conversion of a code where information is represented by a given sequence or number of digits to a code where the same, similar or subset of information is represented by a different sequence or number of digits
- H03M7/30—Compression; Expansion; Suppression of unnecessary data, e.g. redundancy reduction
- H03M7/46—Conversion to or from run-length codes, i.e. by representing the number of consecutive digits, or groups of digits, of the same kind by a code word and a digit indicative of that kind
Landscapes
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Other Investigation Or Analysis Of Materials By Electrical Means (AREA)
- Compression, Expansion, Code Conversion, And Decoders (AREA)
- Chemical Kinetics & Catalysis (AREA)
- Electrochemistry (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
Abstract
Description
例えば特許文献1には、飛行時間型質量分析装置(TOFMS)で得られる高時間分解能のマススペクトルデータを圧縮する方法が開示されている。この圧縮方法では、マススペクトルにおいて質量電荷比(m/z)値の順に強度値を並べたデータ配列を、ノイズレベルであるとみなせるデータ配列とそれ以外のデータ配列とに分ける。そして、ノイズレベルであるとみなしたデータを除いたデータ配列(つまりは有用である可能性が高いデータ配列)に対して連長符号化を行う一方、ノイズレベルであるとみなした部分に対してはエントロピー符号化の一種である静的ハフマン符号化を行う。この圧縮方法では、圧縮率を上げながら可逆圧縮を行うことができる。
マススペクトルを表すデータ配列は、マススペクトルの横軸であるm/z値の情報がそれぞれ格納されたm/z値の配列と、それと同じデータ数である強度値の配列とから構成される。例えばイメージング質量分析装置で得られる多数の測定点毎のマススペクトルデータを保存する場合、或いは、同じ質量分析装置で得られる多数の試料に対するマススペクトルデータを保存する場合、圧縮前のデータ配列は、複数のマススペクトルに共通である一つのm/z値の配列と、複数のマススペクトルにそれぞれ対応する強度値の配列とであることが多い。
a)処理対象である圧縮データの配列の先頭から順に、有意なレベルである強度値の連続数、及び無効なレベルである強度値の連続数、の情報を取得して逐次的に加算し、その加算結果が目的とする特定のパラメータ値に対応する配列番号に一致するとき又は該配列番号を超えるときの連続数を特定する連続数特定ステップと、
b)前記連続数特定ステップで特定された連続数が無効なレベルである強度値の連続数である場合には、前記特定のパラメータ値における強度値として、無効レベルを示す情報又はそれに対応付けられた特定の強度値を出力し、前記特定された連続数が有意なレベルである強度値の連続数である場合には、前記特定のパラメータ値における強度値として、圧縮データの配列上で該連続数を示すデータとその次に現れる連続数を示すデータとの間の強度値の中から該当する配列番号の強度値を見つけて出力する強度値探索ステップと、
を有することを特徴としている。
a)処理対象である圧縮データの配列の先頭から順に、有意なレベルである強度値の連続数、及び無効なレベルである強度値の連続数、の情報を取得して逐次的に加算し、その加算結果が目的とする特定のパラメータ値に対応する配列番号に一致するときの又は該配列番号を超えるときの連続数を特定する連続数特定部と、
b)前記連続数特定部により特定された連続数が無効なレベルである強度値の連続数である場合には、前記特定のパラメータ値における強度値として、無効レベルを示す情報又はそれに対応付けられた特定の強度値を出力し、前記特定された連続数が有意なレベルである強度値の連続数である場合には、前記特定のパラメータ値における強度値として、圧縮データの配列上で該連続数を示すデータとその次に現れる連続数を示すデータとの間の強度値の中から該当する配列番号の強度値を見つけて出力する強度値探索部と、
を備えることを特徴としている。
分析装置がLCやGC等のクロマトグラフ装置である場合、所定のパラメータとは時間であり、測定データとはクロマトグラムデータである。
分析装置がラマン分光分析装置やFTIRである場合、所定のパラメータとは波長(又は波数)であり、測定データとは分光スペクトルデータである。分析装置がEPMAである場合、所定のパラメータとはX線の波長又はエネルギーであり、測定データとはX線スペクトルデータである。
a)測定データの配列を所定の個数のデータ毎にブロックに区切るデータ分割ステップと、
b)ブロック毎に、当該ブロックに含まれる測定データの配列順に従って、予め決められた無効レベルである強度値が一以上連続する場合には、その強度値を一括して強度値の連続数に置き換える一方、前記無効レベル以外の有意なレベルである強度値が出現した場合にはスイッチ連長符号化を行うことにより、各ブロックにおける圧縮後のデータ配列を求める配列変換ステップと、
c)前記圧縮後の各ブロックのデータを並べた1次元配列上での各ブロックの先頭のデータの位置を示す配列番号の情報を集めてインデクスを作成し、前記配列変換ステップで得られた各ブロックの圧縮後のデータ配列と関連付けて保存するインデクス作成ステップと、
を有することを特徴としている。
a)測定データの配列を所定の個数のデータ毎にブロックに区切るデータ分割部と、
b)ブロック毎に、当該ブロックに含まれる測定データの配列順に従って、予め決められた無効レベルである強度値が一以上連続する場合には、その強度値を一括して強度値の連続数に置き換える一方、前記無効レベル以外の有意なレベルである強度値が出現した場合にはスイッチ連長符号化を行うことにより、各ブロックにおける圧縮後のデータ配列を求める配列変換部と、
c)前記圧縮後の各ブロックのデータを並べた1次元配列上での各ブロックの先頭のデータの位置を示す配列番号の情報を集めてインデクスを作成し、前記配列変換部により得られた各ブロックの圧縮後のデータ配列と関連付けて保存するインデクス作成部と、
を備えることを特徴としている。
該圧縮データから特定のパラメータ値における強度値を求めるために、
d)前記インデクスを参照して目的とする特定のパラメータ値におけるデータの配列番号に対応する圧縮前のデータを圧縮した情報が含まれるブロックを特定し、該ブロックの先頭のデータ位置の情報を求めるブロック特定ステップと、
e)前記特定されたブロックに対応する圧縮データの配列の先頭から順に、有意なレベルである強度値の連続数、及び無効なレベルである強度値の連続数、の情報を取得して逐次的に加算し、その加算結果が目的とする特定のパラメータ値に対応する配列番号に一致するとき又は該配列番号を超えるときの連続数を特定する連続数特定ステップと、
f)前記連続数特定ステップで特定された連続数が無効なレベルである強度値の連続数である場合には、前記特定のパラメータ値における強度値として、無効レベルを示す情報又はそれに対応付けられた特定の強度値を出力し、前記特定された連続数が有意なレベルである強度値の連続数である場合には、前記特定のパラメータ値における強度値として、圧縮データの配列上で該連続数を示すデータとその次に現れる連続数を示すデータとの間の強度値の中から該当する配列番号の強度値を見つけて出力する強度値探索ステップと、
を有するものとするとよい。
該圧縮データから特定のパラメータ値における強度値を求めるために、
d)前記インデクスを参照して目的とする特定のパラメータ値におけるデータの配列番号に対応する圧縮前のデータを圧縮した情報が含まれるブロックを特定し、該ブロックの先頭のデータ位置の情報を求めるブロック特定部と、
e)前記特定されたブロックに対応する圧縮データの配列の先頭から順に、有意なレベルである強度値の連続数、及び無効なレベルである強度値の連続数、の情報を取得して逐次的に加算し、その加算結果が目的とする特定のパラメータ値に対応する配列番号に一致するとき又は該配列番号を超えるときの連続数を特定する連続数特定部と、
f)前記連続数特定部により特定された連続数が無効なレベルである強度値の連続数である場合には、前記特定のパラメータ値における強度値として、無効レベルを示す情報又はそれに対応付けられた特定の強度値を出力し、前記特定された連続数が有意なレベルである強度値の連続数である場合には、前記特定のパラメータ値における強度値として、圧縮データの配列上で該連続数を示すデータとその次に現れる連続数を示すデータとの間の強度値の中から該当する配列番号の強度値を見つけて出力する強度値探索部と、
を備える構成とするとよい。
a)処理対象である複数の圧縮データ配列のうち、前記上位のビット列に対応する圧縮データ配列の先頭から順に、「0」でない値を示すビット列の連続数、及び「0」を示すビット列の連続数、の情報を取得して逐次的に加算し、その加算結果が目的とする特定のパラメータ値に対応する配列番号に一致するときの連続数又は該配列番号を超える直前の連続数を特定する連続数特定ステップと、
b)前記特定された連続数が「0」を示すビット列の連続数である場合には、前記特定のパラメータ値における上位のビット列として「0」の値を出力し、前記特定された連続数が「0」でない値を示すビット列の連続数である場合には、前記特定のパラメータ値における上位のビット列として、圧縮データの配列上で該連続数を示すデータとその次に現れる連続数を示すデータとの間のビット列の中から該当する順番のビット列を見つけて出力する上位ビット情報探索ステップと、
c)前記下位のビット列に対応する圧縮データ配列において特定のパラメータ値に対応する配列番号から下位のビット列の情報を取得し、前記上位ビット情報探索ステップにより出力されるビット列の情報と合わせて特定のパラメータ値に対する強度値情報を取得する強度値取得ステップと、
を有することを特徴としている。
a)処理対象である複数の圧縮データ配列のうち、前記上位のビット列に対応する圧縮データ配列の先頭から順に、「0」でない値を示すビット列の連続数、及び「0」を示すビット列の連続数、の情報を取得して逐次的に加算し、その加算結果が目的とする特定のパラメータ値に対応する配列番号に一致するときの連続数又は該配列番号を超える直前の連続数を特定する連続数特定ステップと、
b)前記特定された連続数が「0」を示すビット列の連続数である場合には、前記特定のパラメータ値における上位のビット列として「0」の値を出力し、前記特定された連続数が「0」でない値を示すビット列の連続数である場合には、前記特定のパラメータ値における上位のビット列として、圧縮データの配列上で該連続数を示すデータとその次に現れる連続数を示すデータとの間のビット列の中から該当する順番のビット列を見つけて出力する上位ビット情報探索ステップと、
c)処理対象である複数の圧縮データ配列のうち、前記下位のビット列に対応する圧縮データ配列を全て伸張し、その伸張後のデータ配列から特定のパラメータ値に対応する下位のビット列の情報を取得し、前記上位ビット情報探索ステップにより出力されるビット列の情報と合わせて特定のパラメータ値に対する強度値情報を取得する強度値取得ステップと、
を有することを特徴としている。
a)測定データの配列を所定の個数のデータ毎にブロックに区切るデータ分割ステップと、
b)ブロック毎に、当該ブロックに含まれる測定データのビット列をそれぞれ、所定のレベル以下の下位のビット列とそれ以外の上位のビット列とに分割することで複数のデータ配列に分け、そのブロック毎に、上位のビット列を含むデータ配列に対しては少なくともスイッチ連長符号化を行う一方、下位のビット列を含むデータ配列に対してはそのまま又は静的ハフマン符号化を行うことで、各ブロックにおける圧縮後の複数のデータ配列を求める配列変換ステップと、
c)前記圧縮後の各ブロックの複数のデータ配列について、データを並べた1次元配列上での各ブロックの先頭のデータの位置を示す配列番号の情報を集めてインデクスを作成し、前記配列変換ステップで得られた各ブロックの圧縮後の複数のデータ配列と関連付けて保存するインデクス作成ステップと、
を有することを特徴としている。
a)測定データの配列を所定の個数のデータ毎にブロックに区切るデータ分割部と、
b)ブロック毎に、当該ブロックに含まれる測定データのビット列をそれぞれ、所定のレベル以下の下位のビット列とそれ以外の上位のビット列とに分割することで複数のデータ配列に分け、そのブロック毎に、上位のビット列を含むデータ配列に対しては少なくともスイッチ連長符号化を行う一方、下位のビット列を含むデータ配列に対してはそのまま又は静的ハフマン符号化を行うことで、各ブロックにおける圧縮後の複数のデータ配列を求める配列変換部と、
c)前記圧縮後の各ブロックの複数のデータ配列について、データを並べた1次元配列上での各ブロックの先頭のデータの位置を示す配列番号の情報を集めてインデクスを作成し、前記配列変換部により得られた各ブロックの圧縮後の複数のデータ配列と関連付けて保存するインデクス作成部と、
を備えることを特徴としている。
この第6発明及び第7発明でも、第3発明及び第4発明と同じような効果を達成することができる。
該圧縮データから特定のパラメータ値における強度値を求めるために、
d)前記インデクスのうち上位のビット列に対応するインデクスを参照して目的とする特定のパラメータ値におけるデータの配列番号に対応する圧縮前の上位のビット列を圧縮した情報が含まれるブロックを特定し、該ブロックの先頭のデータ位置の情報を求める第1ブロック特定ステップと、
e)前記第1ブロック特定ステップで特定されたブロックに対応する圧縮データの配列の先頭から順に、「0」でない値を示すビット列の連続数、及び「0」を示すビット列の連続数、の情報を取得して逐次的に加算し、その加算結果が目的とする特定のパラメータ値に対応する配列番号に一致するときの連続数又は該配列番号を超える直前の連続数を特定する連続数特定ステップと、
f)前記連続数特定ステップで特定された連続数が「0」を示すビット列の連続数である場合には、前記特定のパラメータ値における上位のビット列として「0」の値を出力し、前記特定された連続数が「0」でない値を示すビット列の連続数である場合には、前記特定のパラメータ値における上位のビット列として、圧縮データの配列上で該連続数を示すデータとその次に現れる連続数を示すデータとの間のビット列の中から該当する順番のビット列を見つけて出力する上位ビット情報探索ステップと、
g)特定のパラメータ値に対応する配列番号から下位のビット列の情報を取得し、前記上位ビット情報探索ステップにより出力されるビット列の情報と合わせて特定のパラメータ値に対する強度値情報を取得する強度値取得ステップと、
を有するものとするとよい。
該圧縮データから特定のパラメータ値における強度値を求めるために、
d)前記インデクスのうち上位のビット列に対応するインデクスを参照して目的とする特定のパラメータ値におけるデータの配列番号に対応する圧縮前の上位のビット列を圧縮した情報が含まれるブロックを特定し、該ブロックの先頭のデータ位置の情報を求める第1ブロック特定ステップと、
e)前記第1ブロック特定ステップで特定されたブロックに対応する圧縮データの配列の先頭から順に、「0」でない値を示すビット列の連続数、及び「0」を示すビット列の連続数、の情報を取得して逐次的に加算し、その加算結果が目的とする特定のパラメータ値に対応する配列番号に一致するときの連続数又は該配列番号を超える直前の連続数を特定する連続数特定ステップと、
f)前記連続数特定ステップで特定された連続数が「0」を示すビット列の連続数である場合には、前記特定のパラメータ値における上位のビット列として「0」の値を出力し、前記特定された連続数が「0」でない値を示すビット列の連続数である場合には、前記特定のパラメータ値における上位のビット列として、圧縮データの配列上で該連続数を示すデータとその次に現れる連続数を示すデータとの間のビット列の中から該当する順番のビット列を見つけて出力する上位ビット情報探索ステップと、
g)前記インデクスのうち下位のビット列に対応するインデクスを参照して目的とする特定のパラメータ値におけるデータの配列番号に対応する圧縮前の下位のビット列が含まれるブロックを特定し、該ブロックの先頭のデータ位置の情報を求める第2ブロック特定ステップと、
h)前記第2ブロック特定ステップで特定されたブロックに対応する下位ビット由来の圧縮データ配列を全て伸張し、そのあと特定のパラメータ値に対応する配列番号とそのブロックの先頭の配列番号の差分をもとに、伸張後のデータ配列から特定のパラメータ値に対応する下位のビット列の情報を取得し、前記上位ビット情報探索ステップにより出力されるビット列の情報と合わせて特定のパラメータ値に対する強度値情報を取得する強度値取得ステップと、
を有するものとするとよい。
以下の例では、機器分析を実施する分析装置が質量分析装置であり、得られる測定データはマススペクトルデータである。なお、飛行時間型質量分析装置の場合には、分析を実行することでまず飛行時間スペクトルが得られ、飛行時間を質量電荷比に換算することでマススペクトルが得られる。したがって、飛行時間を質量電荷比に換算する前の飛行時間スペクトルを構成するデータもマススペクトルデータとみなすことができる。但し、後述するように、本発明を適用可能な分析装置はこれに限るものではない。
本発明に関連する第1の圧縮方法の原理について、図3を参照して説明する。
図3(C)はマススペクトルデータにより作成されるマススペクトルの一例である。横軸はm/z値、縦軸が強度値であり、このマススペクトルを構成するマススペクトルデータは、離散的なm/z値にそれぞれ対応する強度値の1次元配列で表すことができる。即ち、例えばm/z範囲が1〜2000であるとすると、(m/z 1,m/z 2,m/z 3,…,m/z 1999,m/z 2000)というm/z値の1次元配列と、これに対応する(0,123,3,…,13,0)などの強度値の1次元配列とが、マススペクトルデータとして得られる。
図4は、上述した第1のデータ圧縮方法による圧縮後のデータから、特定のm/z値に対応する強度値を求める際の概略的な手順を説明した図である。
図1は、上述した第1の圧縮方法をデータ保存に利用した、第1実施例のイメージング質量分析システムの概略構成図である。
このイメージング質量分析システムは、試料上の2次元領域内の多数の測定点についてそれぞれ質量分析を実行してマススペクトルデータを取得するイメージング質量分析部1と、データ圧縮処理等のデータ処理を実施するデータ処理部2と、該データ処理部2により圧縮されたデータを保存するデータ保存部3と、ユーザが操作する操作部4と、分析結果等を表示する表示部5と、を備える。データ処理部2の実体はパーソナルコンピュータであり、機能ブロックとして、スペクトルデータ収集部20、メインメモリ21、データ圧縮処理部22、データ伸張処理部23、データ読み出し処理部25、イメージング画像作成処理部26などを含む。また、データ保存部3は圧縮データ記憶領域30を有する。
図7は、第1実施例のイメージング質量分析システムにおいて特定のm/z位置の強度値を求める際にデータ読み出し処理部25で実施される詳細な処理を示すフローチャートである。ここで、Sは非圧縮データ配列上において強度値を取得したいデータが位置する配列番号である。これはm/z値の1次元配列の情報から既知である。Uは現在着目している非圧縮データ配列上における配列番号を示す変数である。また、Pは現在着目している圧縮データ配列上における配列番号を示す変数である。
以上のようにして、圧縮データを伸張することなく、目的とする特定のm/z値に対応する強度値を取得することができる。
イメージング画像を表示する際に、特定の一つのm/z値ではなく特定のm/z範囲に含まれる複数の強度値を積算した強度積算値のイメージング画像を作成したい場合がある。図8は、特定のm/z範囲に含まれる複数の強度値を積算した強度積算値を取得する際にデータ読み出し処理部25で実施される詳細な処理を示すフローチャートである。
以上のようにして、圧縮データを伸張することなく、目的とする特定のm/z範囲に対応する複数の強度値を積算した強度積算値を取得することができる。
次に、本発明に関連する第2の圧縮方法の原理について、図5を参照して説明する。
上記第1の圧縮方法はインデクスを使用しない圧縮方法であるため、データ量を削減するには有利であるものの、特定のm/z値に対応する強度値を取得する際に、圧縮データ配列の先頭から順に連続数を加算しなければならず、m/z範囲が広いマススペクトルデータにおいて高m/z領域の強度値を取得するときに時間が掛かる場合がある。この点を改善したのが第2の圧縮方法である。
図6は、上述した第2のデータ圧縮方法による圧縮後のデータから、特定のm/z値に対応する強度値を求める際の概略的な手順を説明した図である。
図2は、上述した第2の圧縮方法をデータ保存に利用した、第2実施例のイメージング質量分析システムの概略構成図である。
このイメージング質量分析システムにおいて上記第1実施例のイメージング質量分析装置における構成要素と同じ又は相当する構成要素については同じ符号を付して説明を省略する。このイメージング質量分析システムでは、データ処理部2にインデクス作成処理部24が付加され、データ保存部3にインデクス記憶領域31が設けられている。
図9は、第2実施例のイメージング質量分析システムにおいて特定のm/z位置の強度値を求める際にデータ読み出し処理部25で実施される詳細な処理を示すフローチャートである。ここで、S、U、Pはすでに図7で説明した定義と同じである。またNは圧縮前の1ブロックのデータ数である。図5、図6の例ではN=1000である。
図10は、第2実施例のイメージング質量分析システムにおいて特定のm/z範囲に含まれる複数の強度値を積算した強度積算値を取得する際にデータ読み出し処理部25で実施される詳細な処理を示すフローチャートである。Ss、Se、N、C、U、Pの定義はこれまでの説明と同じである。また、図10においてQで示す部分は図8中でQで示すステップと同じである。
一般に複数のマススペクトルデータに対して主成分分析などの多変量解析を行う場合、プロファイルデータの全てを解析に利用するのではなくプロファイルデータから有意なピークの情報(m/z値と強度)を抽出して作成したピーク強度情報について解析を行うことが多い。例えば、マススペクトル上のピークは複数のデータ点で構成されるが、ピークの中心(重心)から特定の許容幅の範囲内に含まれるデータ点の強度値を積算したり、横軸のm/z値を考慮した面積値(積分値)を計算したり、或いは、面積値をピークの幅で除した平均強度値などを求めたりすることで、一つのピークを一つ強度値で表現し、この強度値を行列状に配置したデータ行列に対し多変量解析を実行することがよくある。ピークを代表する一つの強度値として、単にピークトップの強度値を使用することもある。
(1)m/z値リストの作成
まず、解析対象である複数のマススペクトルの全て又は一部を平均して求めた平均スペクトルや複数のマススペクトルについてm/z値毎に最大強度の信号を抽出して作成した最大強度スペクトルを求める。そして、そのスペクトル上で観測される各ピークのm/z範囲を調べ、各ピークの始点及び終点に対応するm/z値のリストを作成する。ピークのm/z範囲は例えばマススペクトル上のピークを構成する複数のデータの中で強度値が所定の閾値を超えるm/z範囲とすればよい。或いは、ピークトップや重心位置の前後に許容幅を設け、その許容幅の範囲をピークのm/z範囲とみなしてピークの始点及び終点を決めてもよい。また、解析に用いるデータ量を削減するために、上述のようにして作成した仮のピークリストの中で、マススペクトル上で強度値が所定値以上であるピークのみを抽出して最終的なm/z値のリストを作成してもよい。また、複数の化合物の質量の理論値や過去の実測値に基づいてm/z値のリストを作成してもよい。
データ行列の各要素である強度値としては、次の積算強度値、ピークトップ強度値、セントロイド強度値のいずれかを用いることができる。
上記のように作成したm/z範囲のリストを元に、解析対象である複数のマススペクトルについての圧縮データ配列から、リストに挙げられている各m/z範囲に対応する積算強度値を求める。各マススペクトルからリストに挙げられている各m/z範囲の積算強度値を計算する際に、前述の強度値積算処理を繰り返し実施すればよい。また、強度値を積算する過程で、m/z値軸上で隣接するデータ点同士のm/z値の差を強度値に乗じて面積値としてもよい。また、積算後の値を積算したデータ点数で除した平均値を各ピークを代表する一つの強度値としてもよい。
ピークを代表する一つの強度値としてピークトップの強度値を使用する場合には、各マススペクトルにおけるリストに挙げられているm/z範囲毎に、次のような最大値計算処理を実施して各m/z範囲内に含まれるピークのピークトップの強度値を取得すればよい。図11は、第2実施例のイメージング質量分析システムにおいてm/z範囲に含まれる一つののピークのピークトップの強度値を取得する際にデータ読み出し処理部25で実施される詳細な処理を示すフローチャートである。Ss、Se、N、C、U、Pの定義はこれまでの説明と同じである。但し、ここではSs、Seは強度値を積算するm/z範囲の上限、下限ではなく、最大強度を探索するためのm/z範囲の上限、下限である。また、Iは最大強度値である。なお、図11においてQで示す部分は図8中でQでまとめて示すステップと同じである。
データ行列を作成する際に、事前に各マススペクトルの全てのピーク情報として、ピークの重心のm/z値と面積値とをセットにした、又は、ピークトップのm/z値と該ピークトップの強度値とをセットにしたセントロイドデータを求めておく。そして、各セントロイドデータに基づいて、データ行列の作成対象であるm/z値リストに挙げられているm/z範囲に含まれるセントロイドの強度値を、データ行列の強度値としてもよい。
上記データ行列は一つのピークに一つの強度値が対応しているデータ行列であるが、ピーク単位ではなく、ピークの有無とは無関係にマススペクトルのm/z軸を連続する複数の区間に区切って、その区間毎に一つの強度値を対応付けてデータ行列を作成してもよい。こうした処理をビニング(厳密にはm/z上のビニング)という。ビニングでは、各区間のm/z値幅は等しくなくてもよいから、例えば、有意なピークが存在する部分では区間を狭め、有意なピークが存在しない部分では区間を広くしてもよい。一つの区間に対応する一つの強度値は、その区間範囲内の強度値の積算値を用いることができる。
[第3の圧縮方法]
本発明に関連する第3の圧縮方法の原理について、図14を参照して説明する。
上記第3の圧縮方法で圧縮されたデータのうち、スイッチ連長符号化された上位6ビットのビット列に対して、第1の圧縮方法で圧縮されたデータに対する処理と同じ処理により、特定のm/z値に対応するデータ配列上のビット列を取得する。一方、下位2ビットがそのまま保存されている場合にはそのデータ配列において特定のm/z値に対応する配列番号からビット列を読み出し、それらのビット情報を結合することで特定のm/z値に対応する強度値を取得する。下位2ビットのビット列に対し静的ハフマン符号化を行っている場合には、スイッチ連長符号化の場合と同様に圧縮されたデータから、特定の配列番号に対応する値を取り出すことはできないため、全てのデータを一旦伸張して特定のm/z値に対応するビット列を取り出し、スイッチ連長符号化されたビット列から取り出したビット情報と結合することで特定のm/z値に対応する強度値を取得する。
次に、本発明に関連する第4の圧縮方法の原理について図15を参照して説明する。この圧縮方法は第3の圧縮方法で用いた可逆圧縮の手法と第2の圧縮方法で用いたブロック分割及びインデクス作成の手法とを組み合わせたものである。
上述したように可逆圧縮された情報のうち、スイッチ連長符号化されたビット列に対して、該圧縮データから特定のm/z値における強度値を求めるために、第2の圧縮方法で圧縮されたデータに対する処理と同様に、インデクス情報を参照して特定のm/z値に対応するビット列の値を取得する。また、下位2ビットのビット列が圧縮されていない場合には、特定のm/z値に対応する配列番号からビット列を読み出し、それらのビット情報を結合することで特定のm/zに対応する強度値を取り出す。下位2ビットのビット列に対し静的ハフマン符号化を行っている場合には、インデクス情報を参照して特定のm/z値に対応するデータが含まれるブロックを特定する。そして、そのブロックに含まれるすべてのデータを伸張し、特定のm/z値に対応する配列番号と、ブロックの先頭の配列番号の差分をもとに、伸張後の配列から特定のm/z値に対応するビット列を取得し、上記スイッチ連長符号化されたビット列から取得した特定のm/z値に対応するビット列と結合して強度値情報を取得すればよい。
(1)複数のマススペクトルを共通の一つのm/z値の配列と圧縮された複数の強度値の配列とで表すことで、圧縮前に比べてデータ容量を削減することができる。
2…データ処理部
20…スペクトルデータ収集部
21…メインメモリ
22…データ圧縮処理部
23…データ伸張処理部
24…インデクス作成処理部
25…データ読み出し処理部
26…イメージング画像作成処理部
3…データ保存部
30…圧縮データ記憶領域
31…インデクス記憶領域
4…操作部
5…表示部
該圧縮データから特定のパラメータ値における強度値を求めるために、
d)前記インデクスのうち上位のビット列に対応するインデクスを参照して目的とする特定のパラメータ値におけるデータの配列番号に対応する圧縮前の上位のビット列を圧縮した情報が含まれるブロックを特定し、該ブロックの先頭のデータ位置の情報を求める第1ブロック特定ステップと、
e)前記第1ブロック特定ステップで特定されたブロックに対応する圧縮データの配列の先頭から順に、「0」でない値を示すビット列の連続数、及び「0」を示すビット列の連続数、の情報を取得して逐次的に加算し、その加算結果が目的とする特定のパラメータ値に対応する配列番号に一致するときの連続数又は該配列番号を超える直前の連続数を特定する連続数特定ステップと、
f)前記連続数特定ステップで特定された連続数が「0」を示すビット列の連続数である場合には、前記特定のパラメータ値における上位のビット列として「0」の値を出力し、前記特定された連続数が「0」でない値を示すビット列の連続数である場合には、前記特定のパラメータ値における上位のビット列として、圧縮データの配列上で該連続数を示すデータとその次に現れる連続数を示すデータとの間のビット列の中から該当する順番のビット列を見つけて出力する上位ビット情報探索ステップと、
g)前記インデクスのうち下位のビット列に対応するインデクスを参照して目的とする特定のパラメータ値におけるデータの配列番号に対応する圧縮前の下位のビット列が含まれるブロックを特定し、該ブロックの先頭のデータ位置の情報を求める第2ブロック特定ステップと、
h)前記第2ブロック特定ステップで特定されたブロックに対応する下位ビット由来の圧縮データ配列を全て伸張し、そのあと特定のパラメータ値に対応する配列番号とそのブロックの先頭の配列番号の差分をもとに、伸張後のデータ配列から特定のパラメータ値に対応する下位のビット列の情報を取得し、前記上位ビット情報探索ステップにより出力されるビット列の情報と合わせて特定のパラメータ値に対する強度値情報を取得する強度値取得ステップと、
を有するものとするとよい。
ピークを代表する一つの強度値としてピークトップの強度値を使用する場合には、各マススペクトルにおけるリストに挙げられているm/z範囲毎に、次のような最大値計算処理を実施して各m/z範囲内に含まれるピークのピークトップの強度値を取得すればよい。図11は、第2実施例のイメージング質量分析システムにおいてm/z範囲に含まれる一つのピークのピークトップの強度値を取得する際にデータ読み出し処理部25で実施される詳細な処理を示すフローチャートである。Ss、Se、N、C、U、Pの定義はこれまでの説明と同じである。但し、ここではSs、Seは強度値を積算するm/z範囲の上限、下限ではなく、最大強度を探索するためのm/z範囲の上限、下限である。また、Iは最大強度値である。なお、図11においてQで示す部分は図8中でQでまとめて示すステップと同じである。
a)処理対象である複数の圧縮データ配列のうち、前記上位のビット列に対応する圧縮データ配列の先頭から順に、「0」でない値を示すビット列の連続数、及び「0」を示すビット列の連続数、の情報を取得して逐次的に加算し、その加算結果が目的とする特定のパラメータ値に対応する配列番号に一致するときの連続数又は該配列番号を超える直前の連続数を特定する連続数特定ステップと、
b)前記特定された連続数が「0」を示すビット列の連続数である場合には、前記特定のパラメータ値における上位のビット列として「0」の値を出力し、前記特定された連続数が「0」でない値を示すビット列の連続数である場合には、前記特定のパラメータ値における上位のビット列として、圧縮データの配列上で該連続数を示すデータとその次に現れる連続数を示すデータとの間のビット列の中から該当する配列番号のビット列を見つけて出力する上位ビット情報探索ステップと、
c)前記下位のビット列に対応する圧縮データ配列において特定のパラメータ値に対応する配列番号から下位のビット列の情報を取得し、前記上位ビット情報探索ステップにより出力されるビット列の情報と合わせて特定のパラメータ値に対する強度値情報を取得する強度値取得ステップと、
を有することを特徴としている。
a)処理対象である複数の圧縮データ配列のうち、前記上位のビット列に対応する圧縮データ配列の先頭から順に、「0」でない値を示すビット列の連続数、及び「0」を示すビット列の連続数、の情報を取得して逐次的に加算し、その加算結果が目的とする特定のパラメータ値に対応する配列番号に一致するときの連続数又は該配列番号を超える直前の連続数を特定する連続数特定ステップと、
b)前記特定された連続数が「0」を示すビット列の連続数である場合には、前記特定のパラメータ値における上位のビット列として「0」の値を出力し、前記特定された連続数が「0」でない値を示すビット列の連続数である場合には、前記特定のパラメータ値における上位のビット列として、圧縮データの配列上で該連続数を示すデータとその次に現れる連続数を示すデータとの間のビット列の中から該当する配列番号のビット列を見つけて出力する上位ビット情報探索ステップと、
c)処理対象である複数の圧縮データ配列のうち、前記下位のビット列に対応する圧縮データ配列を全て伸張し、その伸張後のデータ配列から特定のパラメータ値に対応する下位のビット列の情報を取得し、前記上位ビット情報探索ステップにより出力されるビット列の情報と合わせて特定のパラメータ値に対する強度値情報を取得する強度値取得ステップと、
を有することを特徴としている。
該圧縮データから特定のパラメータ値における強度値を求めるために、
d)前記インデクスのうち上位のビット列に対応するインデクスを参照して目的とする特定のパラメータ値におけるデータの配列番号に対応する圧縮前の上位のビット列を圧縮した情報が含まれるブロックを特定し、該ブロックの先頭のデータ位置の情報を求める第1ブロック特定ステップと、
e)前記第1ブロック特定ステップで特定されたブロックに対応する圧縮データの配列の先頭から順に、「0」でない値を示すビット列の連続数、及び「0」を示すビット列の連続数、の情報を取得して逐次的に加算し、その加算結果が目的とする特定のパラメータ値に対応する配列番号に一致するときの連続数又は該配列番号を超える直前の連続数を特定する連続数特定ステップと、
f)前記連続数特定ステップで特定された連続数が「0」を示すビット列の連続数である場合には、前記特定のパラメータ値における上位のビット列として「0」の値を出力し、前記特定された連続数が「0」でない値を示すビット列の連続数である場合には、前記特定のパラメータ値における上位のビット列として、圧縮データの配列上で該連続数を示すデータとその次に現れる連続数を示すデータとの間のビット列の中から該当する配列番号のビット列を見つけて出力する上位ビット情報探索ステップと、
g)特定のパラメータ値に対応する配列番号から下位のビット列の情報を取得し、前記上位ビット情報探索ステップにより出力されるビット列の情報と合わせて特定のパラメータ値に対する強度値情報を取得する強度値取得ステップと、
を有するものとするとよい。
該圧縮データから特定のパラメータ値における強度値を求めるために、
d)前記インデクスのうち上位のビット列に対応するインデクスを参照して目的とする特定のパラメータ値におけるデータの配列番号に対応する圧縮前の上位のビット列を圧縮した情報が含まれるブロックを特定し、該ブロックの先頭のデータ位置の情報を求める第1ブロック特定ステップと、
e)前記第1ブロック特定ステップで特定されたブロックに対応する圧縮データの配列の先頭から順に、「0」でない値を示すビット列の連続数、及び「0」を示すビット列の連続数、の情報を取得して逐次的に加算し、その加算結果が目的とする特定のパラメータ値に対応する配列番号に一致するときの連続数又は該配列番号を超える直前の連続数を特定する連続数特定ステップと、
f)前記連続数特定ステップで特定された連続数が「0」を示すビット列の連続数である場合には、前記特定のパラメータ値における上位のビット列として「0」の値を出力し、前記特定された連続数が「0」でない値を示すビット列の連続数である場合には、前記特定のパラメータ値における上位のビット列として、圧縮データの配列上で該連続数を示すデータとその次に現れる連続数を示すデータとの間のビット列の中から該当する配列番号のビット列を見つけて出力する上位ビット情報探索ステップと、
g)前記インデクスのうち下位のビット列に対応するインデクスを参照して目的とする特定のパラメータ値におけるデータの配列番号に対応する圧縮前の下位のビット列が含まれるブロックを特定し、該ブロックの先頭のデータ位置の情報を求める第2ブロック特定ステップと、
h)前記第2ブロック特定ステップで特定されたブロックに対応する下位ビット由来の圧縮データ配列を全て伸張し、そのあと特定のパラメータ値に対応する配列番号とそのブロックの先頭の配列番号の差分をもとに、伸張後のデータ配列から特定のパラメータ値に対応する下位のビット列の情報を取得し、前記上位ビット情報探索ステップにより出力されるビット列の情報と合わせて特定のパラメータ値に対する強度値情報を取得する強度値取得ステップと、
を有するものとするとよい。
Claims (18)
- 機器分析により得られた所定のパラメータにおける複数の値に対する強度値の1次元配列である測定データを圧縮するための機器分析用データ処理方法であって、
a)測定データの配列を所定の個数のデータ毎にブロックに区切るデータ分割ステップと、
b)ブロック毎に、当該ブロックに含まれる測定データの配列順に従って、予め決められた無効レベルである強度値が一以上連続する場合には、その強度値を一括して強度値の連続数に置き換える一方、前記無効レベル以外の有意なレベルである強度値が出現した場合にはスイッチ連長符号化を行うことにより、各ブロックにおける圧縮後のデータ配列を求める配列変換ステップと、
c)前記圧縮後の各ブロックのデータを並べた1次元配列上での各ブロックの先頭のデータの位置を示す配列番号の情報を集めてインデクスを作成し、前記配列変換ステップで得られた各ブロックの圧縮後のデータ配列と関連付けて保存するインデクス作成ステップと、
を有することを特徴とする機器分析用データ処理方法。 - 請求項1に記載の機器分析用データ処理方法であって、
前記機器分析は質量分析であり、前記所定のパラメータは質量電荷比であり、前記測定データはマススペクトルを構成するデータであることを特徴とする機器分析用データ処理方法。 - 請求項1に記載の機器分析用データ処理方法であり、前記データ分割ステップ、前記配列変換ステップ、及び前記インデクス作成ステップの各ステップでの処理により圧縮して得られた圧縮データを処理する機器分析用データ処理方法であって、
該圧縮データから特定のパラメータ値における強度値を求めるために、
d)前記インデクスを参照して目的とする特定のパラメータ値におけるデータの配列番号に対応する圧縮前のデータを圧縮した情報が含まれるブロックを特定し、該ブロックの先頭のデータ位置の情報を求めるブロック特定ステップと、
e)前記特定されたブロックに対応する圧縮データの配列の先頭から順に、有意なレベルである強度値の連続数、及び無効なレベルである強度値の連続数、の情報を取得して逐次的に加算し、その加算結果が目的とする特定のパラメータ値に対応する配列番号に一致するときの連続数又は該配列番号を超える直前の連続数を特定する連続数特定ステップと、
f)前記連続数特定ステップで特定された連続数が無効なレベルである強度値の連続数である場合には、前記特定のパラメータ値における強度値として、無効レベルを示す情報又はそれに対応付けられた特定の強度値を出力し、前記特定された連続数が有意なレベルである強度値の連続数である場合には、前記特定のパラメータ値における強度値として、圧縮データの配列上で該連続数を示すデータとその次に現れる連続数を示すデータとの間の強度値の中から該当する順番の強度値を見つけて出力する強度値探索ステップと、
を有することを特徴とする機器分析用データ処理方法。 - 機器分析により得られた所定のパラメータにおける複数の値に対する強度値の1次元配列である測定データを圧縮するための機器分析用データ処理装置であって、
a)測定データの配列を所定の個数のデータ毎にブロックに区切るデータ分割部と、
b)ブロック毎に、当該ブロックに含まれる測定データの配列順に従って、予め決められた無効レベルである強度値が一以上連続する場合には、その強度値を一括して強度値の連続数に置き換える一方、前記無効レベル以外の有意なレベルである強度値が出現した場合にはスイッチ連長符号化を行うことにより、各ブロックにおける圧縮後のデータ配列を求める配列変換部と、
c)前記圧縮後の各ブロックのデータを並べた1次元配列上での各ブロックの先頭のデータの位置を示す配列番号の情報を集めてインデクスを作成し、前記配列変換部により得られた各ブロックの圧縮後のデータ配列と関連付けて保存するインデクス作成部と、
を備えることを特徴とする機器分析用データ処理装置。 - 請求項4に記載の機器分析用データ処理装置であって、
前記機器分析は質量分析であり、前記所定のパラメータは質量電荷比であり、前記測定データはマススペクトルを構成するデータであることを特徴とする機器分析用データ処理装置。 - 請求項4に記載の機器分析用データ処理装置であり、前記データ分割部、前記配列変換部、及び前記インデクス作成部の各部での処理により圧縮して得られた圧縮データを処理する機器分析用データ処理装置であって、
該圧縮データから特定のパラメータ値における強度値を求めるために、
d)前記インデクスを参照して目的とする特定のパラメータ値におけるデータの配列番号に対応する圧縮前のデータを圧縮した情報が含まれるブロックを特定し、該ブロックの先頭のデータ位置の情報を求めるブロック特定部と、
e)前記特定されたブロックに対応する圧縮データの配列の先頭から順に、有意なレベルである強度値の連続数、及び無効なレベルである強度値の連続数、の情報を取得して逐次的に加算し、その加算結果が目的とする特定のパラメータ値に対応する配列番号に一致するときの連続数又は該配列番号を超えるときの連続数を特定する連続数特定部と、
f)前記連続数特定部により特定された連続数が無効なレベルである強度値の連続数である場合には、前記特定のパラメータ値における強度値として、無効レベルを示す情報又はそれに対応付けられた特定の強度値を出力し、前記特定された連続数が有意なレベルである強度値の連続数である場合には、前記特定のパラメータ値における強度値として、圧縮データの配列上で該連続数を示すデータとその次に現れる連続数を示すデータとの間の強度値の中から該当する順番の強度値を見つけて出力する強度値探索部と、
を備えることを特徴とする機器分析用データ処理装置。 - 機器分析により得られた所定のパラメータにおける複数の値に対する強度値の1次元配列である測定データを、その配列順に従って、予め決められた無効レベルの強度値が一以上連続する場合には、その強度値を一括して強度値の連続数に置き換える一方、前記無効レベル以外の有意なレベルの強度値が出現した場合にはスイッチ連長符号化を行うことにより圧縮して得られた圧縮データを処理するデータ処理方法であり、該圧縮データから特定のパラメータ値における強度値を求める機器分析用データ処理方法であって、
a)処理対象である圧縮データの配列の先頭から順に、有意なレベルである強度値の連続数、及び無効なレベルである強度値の連続数、の情報を取得して逐次的に加算し、その加算結果が目的とする特定のパラメータ値に対応する配列番号に一致するときの連続数又は該配列番号を超えるとの連続数を特定する連続数特定ステップと、
b)前記連続数特定ステップで特定された連続数が無効なレベルである強度値の連続数である場合には、前記特定のパラメータ値における強度値として、無効レベルを示す情報又はそれに対応付けられた特定の強度値を出力し、前記特定された連続数が有意なレベルである強度値の連続数である場合には、前記特定のパラメータ値における強度値として、圧縮データの配列上で該連続数を示すデータとその次に現れる連続数を示すデータとの間の強度値の中から該当する順番の強度値を見つけて出力する強度値探索ステップと、
を有することを特徴とする機器分析用データ処理方法。 - 請求項7に記載の機器分析用データ処理方法であって、
前記機器分析は質量分析であり、前記所定のパラメータは質量電荷比であり、前記測定データはマススペクトルを構成するデータであることを特徴とする機器分析用データ処理方法。 - 機器分析により得られた所定のパラメータにおける複数の値に対する強度値の1次元配列である測定データを、その配列順に従って、予め決められた無効レベルの強度値が一以上連続する場合には、その強度値を一括して強度値の連続数に置き換える一方、前記無効レベル以外の有意なレベルの強度値が出現した場合にはスイッチ連長符号化を行うことにより圧縮して得られた圧縮データを処理するデータ処理装置であり、該圧縮データから特定のパラメータ値における強度値を求める機器分析用データ処理装置であって、
a)処理対象である圧縮データの配列の先頭から順に、有意なレベルである強度値の連続数、及び無効なレベルである強度値の連続数、の情報を取得して逐次的に加算し、その加算結果が目的とする特定のパラメータ値に対応する配列番号に一致するときの連続数又は該配列番号を超えるときの連続数を特定する連続数特定部と、
b)前記連続数特定部により特定された連続数が無効なレベルである強度値の連続数である場合には、前記特定のパラメータ値における強度値として、無効レベルを示す情報又はそれに対応付けられた特定の強度値を出力し、前記特定された連続数が有意なレベルである強度値の連続数である場合には、前記特定のパラメータ値における強度値として、圧縮データの配列上で該連続数を示すデータとその次に現れる連続数を示すデータとの間の強度値の中から該当する順番の強度値を見つけて出力する強度値探索部と、
を備えることを特徴とする機器分析用データ処理装置。 - 請求項9に記載の機器分析用データ処理装置であって、
前記機器分析は質量分析であり、前記所定のパラメータは質量電荷比であり、前記測定データはマススペクトルを構成するデータであることを特徴とする機器分析用データ処理装置。 - 機器分析により得られた所定のパラメータにおける複数の値に対する強度値の1次元配列である測定データのビット列を、所定のレベル以下の下位のビット列とそれ以外の上位のビット列とに分割することで複数のデータ配列に分け、上位のビット列を含むデータ配列に対しては少なくともスイッチ連長符号化による圧縮を行う一方、下位のビット列を含むデータ配列に対しては圧縮せずに得られた、複数の圧縮データ配列から成る圧縮データを処理するデータ処理方法であり、該圧縮データから特定のパラメータ値における強度値を求める機器分析用データ処理方法であって、
a)処理対象である複数の圧縮データ配列のうち、前記上位のビット列に対応する圧縮データ配列の先頭から順に、「0」でない値を示すビット列の連続数、及び「0」を示すビット列の連続数、の情報を取得して逐次的に加算し、その加算結果が目的とする特定のパラメータ値に対応する配列番号に一致するときの連続数又は該配列番号を超える直前の連続数を特定する連続数特定ステップと、
b)前記特定された連続数が「0」を示すビット列の連続数である場合には、前記特定のパラメータ値における上位のビット列として「0」の値を出力し、前記特定された連続数が「0」でない値を示すビット列の連続数である場合には、前記特定のパラメータ値における上位のビット列として、圧縮データの配列上で該連続数を示すデータとその次に現れる連続数を示すデータとの間のビット列の中から該当する順番のビット列を見つけて出力する上位ビット情報探索ステップと、
c)前記下位のビット列に対応する圧縮データ配列において特定のパラメータ値に対応する配列番号から下位のビット列の情報を取得し、前記上位ビット情報探索ステップにより出力されるビット列の情報と合わせて特定のパラメータ値に対する強度値情報を取得する強度値取得ステップと、
を有することを特徴とする機器分析用データ処理方法。 - 機器分析により得られた所定のパラメータにおける複数の値に対する強度値の1次元配列である測定データのビット列を、所定のレベル以下の下位のビット列とそれ以外の上位のビット列とに分割することで複数のデータ配列に分け、上位のビット列を含むデータ配列に対しては少なくともスイッチ連長符号化による圧縮を行う一方、下位のビット列を含むデータ配列に対して静的ハフマン符号化による圧縮を行うことで得られた、複数の圧縮データ配列から成る圧縮データを処理するデータ処理方法であり、該圧縮データから特定のパラメータ値における強度値を求める機器分析用データ処理方法であって、
a)処理対象である複数の圧縮データ配列のうち、前記上位のビット列に対応する圧縮データ配列の先頭から順に、「0」でない値を示すビット列の連続数、及び「0」を示すビット列の連続数、の情報を取得して逐次的に加算し、その加算結果が目的とする特定のパラメータ値に対応する配列番号に一致するときの連続数又は該配列番号を超える直前の連続数を特定する連続数特定ステップと、
b)前記特定された連続数が「0」を示すビット列の連続数である場合には、前記特定のパラメータ値における上位のビット列として「0」の値を出力し、前記特定された連続数が「0」でない値を示すビット列の連続数である場合には、前記特定のパラメータ値における上位のビット列として、圧縮データの配列上で該連続数を示すデータとその次に現れる連続数を示すデータとの間のビット列の中から該当する順番のビット列を見つけて出力する上位ビット情報探索ステップと、
c)処理対象である複数の圧縮データ配列のうち、前記下位のビット列に対応する圧縮データ配列を全て伸張し、その伸張後のデータ配列から特定のパラメータ値に対応する下位のビット列の情報を取得し、前記上位ビット情報探索ステップにより出力されるビット列の情報と合わせて特定のパラメータ値に対する強度値情報を取得する強度値取得ステップと、
を有することを特徴とする機器分析用データ処理方法。 - 機器分析により得られた所定のパラメータにおける複数の値に対する強度値の1次元配列である測定データを圧縮するための機器分析用データ処理方法であって、
a)測定データの配列を所定の個数のデータ毎にブロックに区切るデータ分割ステップと、
b)ブロック毎に、当該ブロックに含まれる測定データのビット列をそれぞれ、所定のレベル以下の下位のビット列とそれ以外の上位のビット列とに分割することで複数のデータ配列に分け、そのブロック毎に、上位のビット列を含むデータ配列に対しては少なくともスイッチ連長符号化を行う一方、下位のビット列を含むデータ配列に対してはそのまま又は静的ハフマン符号化を行うことで、各ブロックにおける圧縮後の複数のデータ配列を求める配列変換ステップと、
c)前記圧縮後の各ブロックの複数のデータ配列について、データを並べた1次元配列上での各ブロックの先頭のデータの位置を示す配列番号の情報を集めてインデクスを作成し、前記配列変換ステップで得られた各ブロックの圧縮後の複数のデータ配列と関連付けて保存するインデクス作成ステップと、
を有することを特徴とする機器分析用データ処理方法。 - 請求項11に記載の機器分析用データ処理方法であり、前記データ分割ステップ、前記配列変換ステップ、及び前記インデクス作成ステップの各ステップでの処理により圧縮して得られた圧縮データで、且つ前記配列変換ステップでは下位のビット列に対する圧縮が行われていない圧縮データを処理する機器分析用データ処理方法であって、
該圧縮データから特定のパラメータ値における強度値を求めるために、
d)前記インデクスのうち上位のビット列に対応するインデクスを参照して目的とする特定のパラメータ値におけるデータの配列番号に対応する圧縮前の上位のビット列を圧縮した情報が含まれるブロックを特定し、該ブロックの先頭のデータ位置の情報を求める第1ブロック特定ステップと、
e)前記第1ブロック特定ステップで特定されたブロックに対応する圧縮データの配列の先頭から順に、「0」でない値を示すビット列の連続数、及び「0」を示すビット列の連続数、の情報を取得して逐次的に加算し、その加算結果が目的とする特定のパラメータ値に対応する配列番号に一致するときの連続数又は該配列番号を超える直前の連続数を特定する連続数特定ステップと、
f)前記連続数特定ステップで特定された連続数が「0」を示すビット列の連続数である場合には、前記特定のパラメータ値における上位のビット列として「0」の値を出力し、前記特定された連続数が「0」でない値を示すビット列の連続数である場合には、前記特定のパラメータ値における上位のビット列として、圧縮データの配列上で該連続数を示すデータとその次に現れる連続数を示すデータとの間のビット列の中から該当する順番のビット列を見つけて出力する上位ビット情報探索ステップと、
g)特定のパラメータ値に対応する配列番号から下位のビット列の情報を取得し、前記上位ビット情報探索ステップにより出力されるビット列の情報と合わせて特定のパラメータ値に対する強度値情報を取得する強度値取得ステップと、
を有することを特徴とする機器分析用データ処理方法。 - 請求項11に記載の機器分析用データ処理方法であり、前記データ分割ステップ、前記配列変換ステップ、及び前記インデクス作成ステップの各ステップでの処理により圧縮して得られた圧縮データで、且つ前記配列変換ステップでは下位のビット列に対する静的ハフマン符号化が行われている圧縮データを処理する機器分析用データ処理方法であって、
該圧縮データから特定のパラメータ値における強度値を求めるために、
d)前記インデクスのうち上位のビット列に対応するインデクスを参照して目的とする特定のパラメータ値におけるデータの配列番号に対応する圧縮前の上位のビット列を圧縮した情報が含まれるブロックを特定し、該ブロックの先頭のデータ位置の情報を求める第1ブロック特定ステップと、
e)前記第1ブロック特定ステップで特定されたブロックに対応する圧縮データの配列の先頭から順に、「0」でない値を示すビット列の連続数、及び「0」を示すビット列の連続数、の情報を取得して逐次的に加算し、その加算結果が目的とする特定のパラメータ値に対応する配列番号に一致するときの連続数又は該配列番号を超える直前の連続数を特定する連続数特定ステップと、
f)前記連続数特定ステップで特定された連続数が「0」を示すビット列の連続数である場合には、前記特定のパラメータ値における上位のビット列として「0」の値を出力し、前記特定された連続数が「0」でない値を示すビット列の連続数である場合には、前記特定のパラメータ値における上位のビット列として、圧縮データの配列上で該連続数を示すデータとその次に現れる連続数を示すデータとの間のビット列の中から該当する順番のビット列を見つけて出力する上位ビット情報探索ステップと、
g)前記インデクスのうち下位のビット列に対応するインデクスを参照して目的とする特定のパラメータ値におけるデータの配列番号に対応する圧縮前の下位のビット列が含まれるブロックを特定し、該ブロックの先頭のデータ位置の情報を求める第2ブロック特定ステップと、
h)前記第2ブロック特定ステップで特定されたブロックに対応する下位ビット由来の圧縮データ配列を全て伸張し、そのあと特定のパラメータ値に対応する配列番号とそのブロックの先頭の配列番号の差分をもとに、伸張後のデータ配列から特定のパラメータ値に対応する下位のビット列の情報を取得し、前記上位ビット情報探索ステップにより出力されるビット列の情報と合わせて特定のパラメータ値に対する強度値情報を取得する強度値取得ステップと、
を有することを特徴とする機器分析用データ処理方法。 - 請求項11〜15のいずれか1項に記載の機器分析用データ処理方法であって、
前記機器分析は質量分析であり、前記所定のパラメータは質量電荷比であり、前記測定データはマススペクトルを構成するデータであることを特徴とする機器分析用データ処理方法。 - 機器分析により得られた所定のパラメータにおける複数の値に対する強度値の1次元配列である測定データを圧縮するための機器分析用データ処理装置であって、
a)測定データの配列を所定の個数のデータ毎にブロックに区切るデータ分割部と、
b)ブロック毎に、当該ブロックに含まれる測定データのビット列をそれぞれ、所定のレベル以下の下位のビット列とそれ以外の上位のビット列とに分割することで複数のデータ配列に分け、そのブロック毎に、上位のビット列を含むデータ配列に対しては少なくともスイッチ連長符号化を行う一方、下位のビット列を含むデータ配列に対してはそのまま又は静的ハフマン符号化を行うことで、各ブロックにおける圧縮後の複数のデータ配列を求める配列変換部と、
c)前記圧縮後の各ブロックの複数のデータ配列について、データを並べた1次元配列上での各ブロックの先頭のデータの位置を示す配列番号の情報を集めてインデクスを作成し、前記配列変換部により得られた各ブロックの圧縮後の複数のデータ配列と関連付けて保存するインデクス作成部と、
を備えることを特徴とする機器分析用データ処理装置。 - 請求項17に記載の機器分析用データ処理装置であって、
前記機器分析は質量分析であり、前記所定のパラメータは質量電荷比であり、前記測定データはマススペクトルを構成するデータであることを特徴とする機器分析用データ処理装置。
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
PCT/JP2018/021282 WO2019229998A1 (ja) | 2018-06-01 | 2018-06-01 | 機器分析用データ処理方法及び装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
JPWO2019229998A1 true JPWO2019229998A1 (ja) | 2021-05-13 |
JP7024866B2 JP7024866B2 (ja) | 2022-02-24 |
Family
ID=68696867
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2020522557A Active JP7024866B2 (ja) | 2018-06-01 | 2018-06-01 | 機器分析用データ処理方法及び装置 |
Country Status (4)
Country | Link |
---|---|
US (1) | US12046463B2 (ja) |
JP (1) | JP7024866B2 (ja) |
CN (1) | CN112189136B (ja) |
WO (1) | WO2019229998A1 (ja) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113753977A (zh) * | 2021-09-06 | 2021-12-07 | 北京思源广泰科技有限公司 | 一种数据处理方法及系统 |
CN115269594A (zh) * | 2022-07-20 | 2022-11-01 | 清云智通(北京)科技有限公司 | 一种工业数据处理方法、系统及计算设备 |
CN115754416B (zh) * | 2022-11-16 | 2023-06-27 | 国能大渡河瀑布沟发电有限公司 | 基于边缘计算的水轮发电机局部放电分析系统及方法 |
CN117097906B (zh) * | 2023-10-20 | 2023-12-26 | 河北天英软件科技有限公司 | 一种区域医疗资源高效利用的方法及系统 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2004343299A (ja) * | 2003-05-14 | 2004-12-02 | Dainippon Printing Co Ltd | 時系列信号の符号化装置および記録媒体 |
WO2009069225A1 (ja) * | 2007-11-30 | 2009-06-04 | Shimadzu Corporation | 飛行時間測定装置 |
JP2012038459A (ja) * | 2010-08-04 | 2012-02-23 | Shimadzu Corp | 質量分析装置 |
JP2012169979A (ja) * | 2011-02-16 | 2012-09-06 | Shimadzu Corp | 質量分析データ処理方法及び該方法を用いた質量分析装置 |
WO2014175211A1 (ja) * | 2013-04-22 | 2014-10-30 | 株式会社島津製作所 | イメージング質量分析データ処理方法及びイメージング質量分析装置 |
JP2014215043A (ja) * | 2013-04-22 | 2014-11-17 | 株式会社島津製作所 | イメージング質量分析データ処理方法及びイメージング質量分析装置 |
JP2015513200A (ja) * | 2012-03-30 | 2015-04-30 | マイクロマス ユーケー リミテッド | 質量分析計 |
JP2017118388A (ja) * | 2015-12-25 | 2017-06-29 | 株式会社島津製作所 | 分析データ処理装置 |
Family Cites Families (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7071464B2 (en) * | 2003-03-21 | 2006-07-04 | Dana-Farber Cancer Institute, Inc. | Mass spectroscopy system |
US7228239B1 (en) * | 2004-12-22 | 2007-06-05 | The Mathworks, Inc. | Methods and systems for classifying mass spectra |
GB0608470D0 (en) * | 2006-04-28 | 2006-06-07 | Micromass Ltd | Mass spectrometer |
JP4527753B2 (ja) * | 2007-07-03 | 2010-08-18 | 株式会社エスグランツ | ビット列検索装置、検索方法及びプログラム |
US20100259631A1 (en) * | 2007-10-26 | 2010-10-14 | Fujifilm Corporation | Data compression apparatus, data compression program and image-taking apparatus |
JP4978539B2 (ja) | 2008-04-07 | 2012-07-18 | カシオ計算機株式会社 | 符号化装置、符号化方法及びプログラム。 |
JP5472610B2 (ja) * | 2009-11-26 | 2014-04-16 | 大日本印刷株式会社 | 数値データ列の符号化/復号化の方法および装置 |
US8699022B2 (en) * | 2012-02-07 | 2014-04-15 | Materialytics, LLC | Methods and systems for analyzing samples |
US8836548B1 (en) * | 2013-12-05 | 2014-09-16 | Emc Corporation | Method and system for data compression at a storage system |
KR101725224B1 (ko) * | 2014-10-06 | 2017-04-11 | 삼성디스플레이 주식회사 | 축합환 화합물 및 이를 포함한 유기 발광 소자 |
US9748972B2 (en) * | 2015-09-14 | 2017-08-29 | Leco Corporation | Lossless data compression |
WO2018037491A1 (ja) * | 2016-08-24 | 2018-03-01 | 株式会社島津製作所 | イメージング質量分析装置 |
-
2018
- 2018-06-01 WO PCT/JP2018/021282 patent/WO2019229998A1/ja active Application Filing
- 2018-06-01 US US17/059,717 patent/US12046463B2/en active Active
- 2018-06-01 CN CN201880093689.3A patent/CN112189136B/zh active Active
- 2018-06-01 JP JP2020522557A patent/JP7024866B2/ja active Active
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2004343299A (ja) * | 2003-05-14 | 2004-12-02 | Dainippon Printing Co Ltd | 時系列信号の符号化装置および記録媒体 |
WO2009069225A1 (ja) * | 2007-11-30 | 2009-06-04 | Shimadzu Corporation | 飛行時間測定装置 |
JP2012038459A (ja) * | 2010-08-04 | 2012-02-23 | Shimadzu Corp | 質量分析装置 |
JP2012169979A (ja) * | 2011-02-16 | 2012-09-06 | Shimadzu Corp | 質量分析データ処理方法及び該方法を用いた質量分析装置 |
JP2015513200A (ja) * | 2012-03-30 | 2015-04-30 | マイクロマス ユーケー リミテッド | 質量分析計 |
WO2014175211A1 (ja) * | 2013-04-22 | 2014-10-30 | 株式会社島津製作所 | イメージング質量分析データ処理方法及びイメージング質量分析装置 |
JP2014215043A (ja) * | 2013-04-22 | 2014-11-17 | 株式会社島津製作所 | イメージング質量分析データ処理方法及びイメージング質量分析装置 |
JP2017118388A (ja) * | 2015-12-25 | 2017-06-29 | 株式会社島津製作所 | 分析データ処理装置 |
Also Published As
Publication number | Publication date |
---|---|
US20210210319A1 (en) | 2021-07-08 |
JP7024866B2 (ja) | 2022-02-24 |
CN112189136A (zh) | 2021-01-05 |
WO2019229998A1 (ja) | 2019-12-05 |
US12046463B2 (en) | 2024-07-23 |
CN112189136B (zh) | 2024-01-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US9859917B2 (en) | Enhanced data compression for sparse multidimensional ordered series data | |
JP7024866B2 (ja) | 機器分析用データ処理方法及び装置 | |
US8498989B2 (en) | Mass analysis data processing method and mass spectrometer using the same | |
US9412571B2 (en) | Imaging mass spectrometric data processing method and imaging mass spectrometer | |
US10312067B2 (en) | Imaging mass analysis data processing method and imaging mass spectrometer | |
CN101878423B (zh) | 飞行时间测定装置 | |
US7667630B2 (en) | Information compression-encoding device, its decoding device, method thereof, program thereof, and recording medium storing the program | |
US10553413B2 (en) | Mass spectrometer | |
EP3444952A1 (en) | Data compression apparatus, data decompression apparatus, data compression program, data decompression program, data compression method, and data decompression method | |
JP2000252832A (ja) | データ圧縮装置、およびデータ圧縮プログラムを記録した記録媒体 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20200930 |
|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20200930 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20201105 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20211005 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20211117 |
|
TRDD | Decision of grant or rejection written | ||
A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20220111 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20220124 |
|
R151 | Written notification of patent or utility model registration |
Ref document number: 7024866 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R151 |