ANALYSIS DATA PROCESSING METHOD AND DEVICE

Abstract

When conducting imaging mass analysis for a region to be measured on a sample, an individual reference value calculating part obtains a maximum value in P.sub.i/I.sub.i of respective measuring points, and stores the value together with measured data as an individual reference value. When performing comparison analysis for a plurality of the data obtained from different samples, a common reference value determining part reads out corresponding a plurality of the individual reference values and determines a minimum value as a common reference value Fmin. A normalization calculation processing part normalizes the respective intensity values by multiplying the intensity values read out from an external memory device by a normalization coefficient long_Max×(Fmin/P.sub.i) obtained from the common reference value Fmin, TIC values Pi at the respective measuring points, and a maximum allowable value long_Max of a variable storing the intensity values at the time of operation.

Claims

1-8. (canceled)

9. An analytical data processing method that handles a spectrum data set including spectrum data collected by performing predetermined analyses respectively on a plurality of the measuring points in a predetermined region to be measured set on the sample and indicating a relationship between a predetermined parameter and a signal intensity value, being associated with spatial position information of measuring points, the method comprising: for a plurality of the spectrum data sets different from each other, a) an individual reference value acquiring step for obtaining a maximum reference value at which intensity values after normalization do not exceed a predetermined allowable value when normalizing the intensity values for all or part of parameter values in spectrum data at respective measuring points for each of the spectrum data sets; b) a common reference value selecting step for selecting a minimum reference value from a plurality of the reference values obtained in the individual reference value acquiring step as a common reference value or setting a value smaller than the minimum value as the common reference value for all the spectrum data sets; and c) a normalization executing step for normalizing the intensity values for all or part of the parameter values in the spectrum data contained in each of the spectrum data sets based on the common reference value selected in the common reference value selecting step.

10. An analytical data processing method that handles a spectrum data set including spectrum data collected by performing predetermined analyses respectively on a plurality of the measuring points in a predetermined region to be measured set on the sample and indicating a relationship between a predetermined parameter and a signal intensity value, being associated with spatial position information of measuring points, the method comprising: for a plurality of the spectrum data sets different from each other, a) an individual reference value acquiring step for obtaining a value at which a ratio between an integrated value for an intensity value and a maximum intensity value for all or part of parameter values in the spectrum data at one measuring point becomes the largest among those at all the measuring points as a reference value for normalization for each of the spectrum data sets; b) a common reference value selecting step for selecting a minimum reference value from a plurality of the reference values obtained in the individual reference value acquiring step as a common reference value or setting a value smaller than the minimum value as the common reference value for all the spectrum data sets; and c) a normalization executing step for normalizing the intensity values for all or part of the parameter values in the spectrum data included in each of the spectrum data sets so that the intensity values after normalization do not exceed a predetermined allowable value based on the common reference value selected in the common reference value selecting step.

11. The analytical data processing method according to claim 9, wherein the predetermined allowable value is a maximum possible value that a variable storing the intensity values can take when executing a normalization operation in the normalization executing step.

12. The analytical data processing method according to claim 9, wherein acquisition of the reference value for each of the spectrum data sets in the individual reference value acquiring step is performed when the spectrum data included in the spectrum data set are collected by executing the predetermined analyses for the plurality of measuring points in the predetermined region to be measured.

13. The analytical data processing method according to claim 9, wherein in the common reference value selecting step, the plurality of reference values acquired in the individual reference value acquiring step for all the spectrum data sets are presented to a user and a value according to a user instruction based on presented information is set as the common reference value.

14. The mass analytical data processing method according to claim 9, wherein the spectrum data is mass spectrum data which indicates a relationship between a mass charge ratio and an ion intensity value.

15. An analytical data processing apparatus that handles a spectrum data set including spectrum data collected by performing predetermined analyses respectively on a plurality of the measuring points in a predetermined region to be measured set on the sample and indicating a relationship between a predetermined parameter and a signal intensity value, being associated with spatial position information of measuring points, the apparatus comprising: for a plurality of the spectrum data sets different from each other, a) an individual reference value acquiring part for obtaining a maximum reference value at which the intensity values after a normalization do not exceed a predetermined allowable value when normalizing intensity values for all or part of parameter values in the spectrum data at the respective measuring points for each of the spectrum data sets; b) a common reference value selecting part for selecting a minimum reference value from a plurality of the reference values obtained by the individual reference value acquiring part as a common reference value or setting a value smaller than the minimum value as the common reference value for all the spectrum data sets; and c) a normalization executing part for normalizing the intensity values for all or part of parameter values in the spectrum data contained in each of the spectrum data sets based on the common reference value selected by the common reference value selecting part.

16. An analytical data processing apparatus that handles a spectrum data set including spectrum data collected by performing predetermined analyses respectively on a plurality of the measuring points in a predetermined region to be measured set on the sample and indicating a relationship between a predetermined parameter and a signal intensity value, being associated with spatial position information of measuring points, the apparatus comprising: for a plurality of the spectrum data sets different from each other, a) an individual reference value acquiring part for obtaining a value at which a ratio between an integrated value for an intensity value and a maximum intensity value for all or part of parameter values in spectrum data at one measuring point becomes the largest among those at all the measuring points as a reference value for normalization for each of the spectrum data sets; b) a common reference value selecting part for selecting a minimum reference value from a plurality of the reference values obtained by the individual reference value acquiring part as a common reference value or setting a value smaller than the minimum value as the common reference value for all the spectrum data sets; and c) a normalization executing part for normalizing the intensity values for all or part of the parameter values in the spectrum data included in each of the spectrum data sets so that intensity values after normalization do not exceed a predetermined allowable value based on the common reference value selected by the common reference value selecting part.

17. The analytical data processing method according to claim 10, wherein the predetermined allowable value is the maximum possible value that a variable storing the intensity values can take when executing the normalization operation in the normalization executing step.

18. The analytical data processing method according to claim 10, wherein acquisition of the reference value for each of the spectrum data sets in the individual reference value acquiring step is performed when the spectrum data included in the spectrum data set is collected by executing the predetermined analyses for the plurality of measuring points in the predetermined region to be measured.

19. The analytical data processing method according to claim 10, wherein in the common reference value selecting step, a plurality of reference values acquired in the individual reference value acquiring step for all the spectrum data sets are presented to a user and a value according to a user instruction based on presented information is set as the common reference value.

20. The mass analysis analytical data processing method according to claim 10, wherein the spectrum data is mass spectrum data which indicates a relationship between a mass charge ratio and an ion intensity values.

21. An analytical data processing method that handles data including spectrum data collected by performing predetermined analyses respectively on a plurality of the measuring points in a predetermined region to be measured set on the sample and indicating a relationship between a predetermined parameter and a signal intensity value, being associated with spatial position information of measuring points, the method comprising: a) an individual reference value acquiring step for obtaining a reference value at which intensity values after the normalization do not exceed a predetermined allowable value when normalizing the intensity values for all or part of parameter values in spectrum data at one measuring point for each of the measuring points; b) a common reference value selecting step for selecting the minimum reference value from a plurality of the reference values for the spectrum data at the respective measuring points obtained in the individual reference value acquiring step as a common reference value or setting a value smaller than the minimum value as the common reference value; and c) a normalization executing step for normalizing the intensity values for all or part of the parameter values in the respective spectrum data based on the common reference value selected in the common reference value selecting step.

22. An analytical data processing apparatus that handles data including spectrum data collected by performing predetermined analyses respectively on a plurality of the measuring points in a predetermined region to be measured set on the sample and indicating a relationship between a predetermined parameter and a signal intensity value, being associated with spatial position information of measuring points, the apparatus comprising: a) an individual reference value acquiring part for obtaining a reference value at which intensity values after the normalization do not exceed a predetermined allowable value when normalizing the intensity values for all or part of parameter values in spectrum data at one measuring point for each of the measuring points; b) a common reference value selecting part for selecting the minimum reference value from the plurality of reference values for the spectrum data at the respective measuring points obtained in the individual reference value acquiring part as a common reference value or setting a value smaller than the minimum value as the common reference value; and c) a normalization executing part for normalizing the intensity values for all or part of the parameter values in the respective spectrum data based on the common reference value selected by the common reference value selecting part.

23. An analytical data processing method that handles data including spectrum data collected by performing predetermined analyses respectively on a plurality of the measuring points in a predetermined region to be measured set on the sample and indicating a relationship between a predetermined parameter and a signal intensity value, being associated with spatial position information of measuring points, the method comprising: a) an individual reference value acquiring step for obtaining a value at which a ratio between an integrated value for the intensity value and the maximum intensity value for all or part of the parameter values in a spectrum data at one measuring point becomes the largest among those at all the measuring points as a reference value for normalization for each of the measuring points; b) a common reference value selecting step for selecting the minimum reference value from the plurality of reference values for the spectrum data at the respective measuring points obtained in the individual reference value acquiring step as a common reference value or setting a value smaller than the minimum value as the common reference value; and c) a normalization executing step for normalizing intensity values for all or part of the parameter values in the respective spectrum data so that intensity values after normalization do not exceed a predetermined allowable value based on the common reference value selected in the common reference value selecting step.

24. An analytical data processing apparatus that handles a data including spectrum data collected by performing predetermined analyses respectively on a plurality of the measuring points in a predetermined region to be measured set on the sample and indicating a relationship between a predetermined parameter and a signal intensity value, being associated with spatial position information of measuring points, the apparatus comprising: a) an individual reference value acquiring part for obtaining a value at which a ratio between an integrated value for the intensity value and the maximum intensity value for all or part of the parameter values in a spectrum data at one measuring point becomes the largest among those at all the measuring points as a reference value for normalization for each of the measuring points; b) a common reference value selecting part for selecting the minimum reference value from a plurality of reference values for the spectrum data at the respective measuring points obtained in the individual reference value acquiring part as a common reference value or setting a value smaller than the minimum value as the common reference value; and c) a normalization executing part for normalizing the intensity values for all or part of the parameter values in the respective spectrum data so that intensity values after normalization do not exceed a predetermined allowable value based on the common reference value selected in the common reference value selec

Description

BRIEF DESCRIPTION OF DRAWINGS

[0051] FIG. 1 is a schematic configuration drawing of illustrating an imaging mass analyzing device, which is an example of an analytical data processing apparatus of the prevent invention.

[0052] FIG. 2 is a flowchart of a data processing which is characteristic at the time of executing measurement according to an imaging mass analyzing device of the example.

[0053] FIG. 3 is a flowchart of a normalization processing for a plurality of imaging mass analysis data in the imaging mass analyzing device of the example.

[0054] FIG. 4 is schematic explanatory drawings of the normalization processing for the plurality of imaging mass analysis data to be executed in the imaging mass analyzing device of the example,

[0055] FIG. 5 is a schematic explanatory drawing illustrating data acquired by imaging mass analysis and preparation and display of a two-dimensional imaging image preparation and display based on the acquired data.

DESCRIPTION OF EMBODIMENTS

[0056] Hereinafter, an analytical data processing method and apparatus of the present invention will be described with an example of mass analysis imaging data acquired by an imaging mass analyzing device.

[0057] First of all, a method of normalization processing for a plurality of the mass analysis imaging data, which is one of characteristics of the present invention, will be described with reference to FIG. 4.

[0058] Here, a case of normalizing the mass analysis imaging data acquired from three different samples A, B, C as illustrated in FIGS. 4(a), (b), and (c) will be described as an example. As described in conjunction with FIG. 5, one mass analysis imaging data is a set of mass spectrum data in a predetermined mass charge ratio range at a number of measuring points 102 in a predetermined region to be measured 101 on a sample 100. In the three mass analysis imaging data, all the size of the measuring point (interval between measuring points), the number of measuring points, the mass charge ratio range, and mass resolution do not have to be provided.

[0059] In one mass analysis imaging data, a maximum intensity value I.sub.i in a mass spectrum at i.sup.th (i=1 to N, where N is a total number of the measuring points 102 in the region to be measured 101) measuring point, and P.sub.i is a TIC value, which is a sum of all intensity values on the mass spectra. Ps.sub.j is the TIC value which is a reference of normalization in the mass analysis imaging data, that is, an individual reference value. Here, j is a sequential number allocated to the mass analysis imaging data, and a sample A is assumed to be j=1, a sample B is assumed to be j=2, and a sample C is assumed to be j=3 as illustrated in FIG. 4. In addition, long_Max is a maximum allowable value that a variable storing the intensity values may take when performing arithmetic processing for normalization on a computer.

[0060] In this case, if the following expression (1) is satisfied for all the measuring points in the one mass analysis imaging data, overflow does not occur during operation of normalization.

I.sub.i×(Ps.sub.j/P.sub.i)<long_Max   (1)

If the expression (1) is modified,

Ps.sub.j<long_Max×(P.sub.i/I.sub.i)   (2)

is satisfied.

[0061] In FIG. 4(a), mass spectra for two measuring points, i=8, i=16 are exemplified. From the former, I.sub.8 and P.sub.8 are obtained, and from the latter, I.sub.16, P.sub.16 are obtained.

[0062] When the individual reference value Ps.sub.j used in TIC normalization satisfies the expression (2), overflow does not occur at the time of normalization operation. In other words, at the time of the TIC normalization, normalization without causing overflow is achieved by performing the normalization so that a given value not larger than the maximum value from among long_Max×(P.sub.i/Ii) for all the i are determined to be an individual reference value Ps.sub.j, and the TIC value at the respective measuring points to be equalized to the individual reference value PS.sub.j. However, if the individual reference value Ps.sub.j is too small, the round-off errors at the time of normalization operation will be increased. Therefore, the individual reference value Ps.sub.j is preferably determined to a value closer to the maximum value of long_Max×(P.sub.i/I.sub.i) as much as possible. Therefore, generally, the maximum value of long_Max×(P.sub.i/I.sub.i) may be determined to the individual reference value Ps.sub.j. In this case, an element of long_Max is contained in the individual reference value Ps.sub.j, and in normalization based thereon, it is guaranteed that no overflow occurs.

[0063] Based on the above-described expression (2), individual reference values Ps.sub.j which do not have probability of causing overflow in the normalization operation are obtained respectively for the plurality of mass analysis imaging data which are desired to be compared to each other. In the example in FIG. 4, three individual reference values Ps.sub.1, Ps.sub.2, and Ps.sub.3 are obtained. Generally, the individual reference values Ps.sub.j are different from each other. Accordingly, the individual reference values Ps.sub.j obtained respectively for the plurality of mass analysis imaging data are compared and the one having the minimum value is selected. The selected individual reference value is determined to be a common reference value Pmin. The common reference value Pmin is determined by considering the TIC values and the maximum intensity values for all the measuring points on the plurality of mass analysis imaging data to be compared. Therefore, in all the mass analysis imaging data to be compared, the respective intensity values of the mass spectra are normalized so that the TIC values at the respective measuring points become the common reference value Pmin, and thus the plurality of mass analysis imaging data may be normalized by the same reference without causing overflow. The round-off errors generated at the time of the normalization operation may be minimized. Specifically, the intensity values may be corrected by multiplying the intensity values of the mass spectra at the respective measuring points by a normalization coefficient G which is obtained by the following expression (3).

G=Pmin/P.sub.i   (3)

[0064] Note that although the description given above supposes a case where the TIC normalization is performed, the principle is the same also in the case where XIC normalization is performed. The reference value used for normalization is, in the case of the TIC normalization, an integrated value for the intensity value over an entire mass charge ratio range of the mass spectra, and in the XIC normalization, the intensity value in the specific mass charge ratio or the integrated value for the intensity value in the range of the specific mass charge ratio. Therefore, a common reference value may be determined based on the description given above after defining the mass charge ratio and the range thereof, and the normalization may be performed basso on the common reference value. Furthermore, the same applies in the case of normalization based on p-norm as disclosed in PTL 2, and the TIC value P.sub.i is replaced by the o-norm.

[0065] In order to perform the above-described arithmetic processing on a computer, the plurality of mass analysis imaging data to be compared are read from, for example, an external memory device into a main memory in the computer or a resister or the like in the interior of a CPU, and the maximum value of long_Max×(P.sub.i/I.sub.i) is obtained for each of the mass analysis imaging data, so that an individual reference value Ps.sub.j is obtained, Then, Pmin, which is the minimum value, is obtained out of the individual reference value Ps.sub.j for the plurality of mass analysis imaging data, and based on the obtained Pmin, the respective intensity values are normalized according to the expression (3). However, the processing procedure descried above has the following problem in practice.

[0066] In other words, there is a case where the total data quantity of the mass analysis imaging data may become as much as several tens of GB as described above. Since a capacity of the main memory of the computer is limited, it is generally difficult to read all the plurality of mass analysis imaging data to be compared into the main memory when comparing the data. Therefore, the plurality of mass analysis imaging data are read into the main memory one by one and an individual reference value Ps.sub.j is obtained, and these procedures are repeated to obtain a plurality of the individual reference value PS.sub.j. Then, it is necessary to select a common reference value Pmin, then read the plurality of mass analysis imaging data into the main memory one by one, and perform normalization with reference to the common reference value Pmin. In this case, a procedure of reading the same data into the main memory is required at the time of calculation of the individual reference values PS.sub.j and at the time of normalization of the respective intensity values. In other words, the operation of reading the same data is required twice, and thus data analysis needs more time.

[0067] Accordingly, in order to avoid such a problem, the processing may be performed in the following procedure. In other words, an individual reference value Psj is calculated for each of the mass analysis imaging data when collecting data by executing measurement of the sample in the imaging mass analyzing device or when performing pre-processing such as noise removal for the data to be analyzed. Subsequently, the individual reference values Ps.sub.j, for example, are stored in the same data file as the mass analysis imaging data or in another data associated to the data files. When performing comparison analysis for the plurality of mass analysis imaging data, the files in which the individual reference values Ps.sub.j for the plurality of mass analysis imaging data to be compared are stored are read out, or only the individual reference values Ps.sub.j are read out from these files to determine a common reference value Pmin. Subsequently, individual mass analysis imaging data are read in the main memory and normalization with reference to the common reference value Pmin is performed. Accordingly, an operation of reading the mass analysis imaging data into the main memory for obtaining the individual reference values Ps.sub.j for normalization is no longer necessary, and time required for at least the normalization operation processing may be shortened.

[0068] Subsequently, a configuration and an operation of an example of the imaging mass analyzing device which performs a normalization processing of the data described above will be described. FIG. 1 is a configuration drawing of a principal portion of the imaging mass analyzing device of this example.

[0069] The imaging mass analyzing device includes an imaging mass analyzing part 1, a data processing part 2 configured to execute various types of data processing as described later on the obtained data, a large-capacity external memory device 3 such as a hard disk drive (HDD) or a solid state drive (SSD) for storing the mass spectrum data acquired by the imaging mass analyzing part 1, an operating unit 4 to be operated by an analyst, and a display unit 5 configured to display results of analysis and the like.

[0070] The imaging mass analyzing part 1, although not illustrated, includes a microscopic observation part for observing optical microscopic images on a sample, and a mass analysis part including an atmospheric pressure MALDI ion source, an ion trap, and a flight time type mass spectrometer (TOFMS). The imaging mass analyzing part 1 is capable of performing mass analysis respectively for a number of measuring points (minute regions) 102 in the region to be measured 101 respectively on the sample 100 specified by the analyst to obtain mass spectrum data over a predetermined mass charge ratio range for each of the measuring points. The mass spectrum data may be MS/MS spectrum (in addition, MS.sup.n spectrum having an n value not smaller than 3) data obtained by dissociating ion corresponding to, for example, given peaks appeared on the mass spectra in the ion trap and performing mass analysis on ion generated thereby in addition to those obtained by analysis without involving ion dissociation.

[0071] The data processing part 2 may include function blocks such as a data collecting part 20, an individual reference value calculating part 21, a common reference value determining part 22, a normalization calculation processing part 23, a peak matrix creating part 24, an imaging image creating part 25, a statistical analysis 26 a data file creating part 27, a main memory 28, and a display processing part 29. The materiality of the data processing part 2 is a personal computer including a CPU, a RAM, a ROM and the like (or a work station with higher performance), and functions of the respective blocks described above except for the main memory 28 are realized by executing dedicated data analysis software installed in advance in the computer.

[0072] FIG. 2 is a flowchart of characteristic data processing when performing measurement on one sample in the imaging mass analyzing device of the example, and FIG. 3 is a flowchart of normalization processing for the plurality of imaging mass analysis data. The characteristic data processing operation in the imaging mass analyzing device of the example will be described with reference to these drawings.

[0073] When a sample to be measured is set in the imaging mass analyzing part 1, the microscopic observation part, not illustrated, images a high-resolution microscopic image on the sample, and the image is displayed on a screen of the display unit 5 via the data processing part 2. The analyst specifies a region to be measured on the sample by the operating unit 4 with reference to the image. Then, a control part, not illustrated, controls the imaging mass analyzing part 1 so that the mass analysis are executed in sequence for a number of the measuring points in the specified region to be measured. In other words, a variable i which specifies a measuring point is set to 1 (Step S11), and measurement for ith measuring point is executed (Step S12). Accordingly, the mass spectrum data at the ith measuring point is obtained, and thus the data collecting part 20 stores the obtained data once in the main memory 28 (Step S13).

[0074] In the data processing part 2, the individual reference value calculating part 21 calculates the largest intensity value I.sub.i and TIC value P.sub.i that is integrated with an intensity value over the entire mass charge ratio from the mass spectrum data obtained for the ith measuring point, and these values are stored therein temporarily in correspondence with the value of i (Step S14). Therefore, when i=1, the maximum intensity value I.sub.1 and the TIC value P.sub.1 are stored.

[0075] Subsequently, the control part determines whether or not the variable i reaches N, which is a total number of the measuring points in the region to be measured (Step S15). If the variable i does not reach a total number N of the measuring points, the value of i is incremented (Step S16) and the procedure goes back to Step S12. In contrast, if the variable i reaches the total number N of the measuring points, the procedure goes from Steps S15 to S17. Therefore, the processing from steps S12 to S16 are repeated until the variable i reaches the total number N of the measuring points. Accordingly, mass analysis for all the measuring points in the region to be measured specified by the analyst is performed, and the mass spectrum data for all the measuring points, that is, the mass analysis imaging data for the sample, are stored in the main memory 28. In contrast, the maximum intensity value I.sub.i and the TIC value P.sub.i for i=1 to N are memorized in the individual reference value calculating part 21. Subsequently, the individual reference value calculating part 21 calculates P.sub.i/I.sub.i for all the measuring points (that is, i=1 to N), and selects the maximum value therefrom as a unsealed individual reference value F.sub.j (Step S17). The unsealed individual reference value F.sub.i is different from the above-described individual reference value Ps.sub.j, and does not include an element of long_Max.

[0076] The data file creating part 27 stores the mass analysis imaging data for one sample stored in the main memory 28, the unsealed individual reference value F.sub.j obtained by the individual reference value calculating part 21, and the TIC value P.sub.i obtained at all the measuring points of i=1 to N in one data file, and the data file is stored in the external memory device 3.

[0077] Although not illustrated in FIG. 2, when multiple classification analyses such as principal component analysis (PCA), principal component regression analysis (PCR), partial least square method (PLS), hierarchical cluster analysis (HCA), hypothesis test, and support vector machine (SVM) are performed for the mass analysis imaging data, it is convenient to prepare a peak matrix based on peak information extracted from the mass spectra (the mass charge ratio value and the intensity values of peaks). Accordingly, in parallel to storage of the data obtained by measurement of the sample in the main memory 28, the peak matrix creating part 24 performs peak detection for the mass spectrum data as needed to prepare a peak matrix including the respective intensity values of the detected peaks arranged in the mass charge ratio value. Then, the data file creating part 27 may store the peak matrix in the same data file as the mass analysis imaging data.

[0078] When performing peak detection for the mass spectra, by narrowing the conditions for peak detection such as peaks having intensity values not smaller than a predetermined threshold, value, or peaks appearing in a specific range of mass charge ratio, the number of detected peaks is reduced, so that the size of the peak matrix data is reduced. Therefore, if reduction of processing time required for multivariable analysis for comparison analysis described later or data quantity to be read in the main memory 28 at the time of multivariable analysis is desired, the peak detection conditions may be narrowed as described above to be more strict to reduce the peak matrix data size.

[0079] When XIC normalization is performed instead of TIC normalization, the unsealed individual reference value may be calculated in the same manner by using the intensity values in a specific mass charge ratio in the obtained mass spectra or the integrated value for the intensity value in a specific mass charge ratio range instead of the TCI value P.sub.i and stored in the data file. The same applies to a case where normalization based on the p-norm as disclosed in PTL 2. In addition, in the case where TIC normalization and XIC normalization are desired to be performed selectively at the time of comparison analysis, or in the case where switching between the mass charge ratio value and the mass charge ratio range focused also in XIC normalization is desired, the unsealed individual reference values corresponding to the respective conditions are calculated and are stored together in the data file.

[0080] The flowchart illustrated in FIG. 2 is for a case where measurement is performed for one sample, and processing according to the flowchart illustrated in FIG. 2 is performed separately for another sample. Therefore, every time when the sample is measured, or every time when measurement for another region to be measured is performed, a data file in which the mass analysis imaging data, the unsealed individual reference value, and TIC values for the respective measuring points for the sample or the region to be measured are stored is created and is stored in the external memory device 3.

[0081] Subsequently, data processing to be executed when comparison analysis of the plurality of mass analysis imaging data is performed in a state in which data files for a plurality of samples are stored in the external memory device 3 as described above will be described with reference to FIG. 3.

[0082] In this case, the analyst specifies the plurality of mass analysis imaging data to be compared by the operating unit 4 (Step S21). FIG. 4, which has already been described, illustrates an example in which the mass analysis imaging data for the three samples A, B, C are specified. Prior to comparison analysis therefor, the common reference value determining part 22 reads out unsealed individual reference values F.sub.j from the data file stored in the external memory device 3 and corresponding to the plurality of mass analysis imaging data (Step S22). Then, the minimum value is obtained from a read plurality of unsealed individual reference values F.sub.j and is stored in the interior as an unsealed common reference value Fmin (Step S23).

[0083] Subsequently, the normalization calculation processing part 23 reads out mass analysis imaging data and TIC values P.sub.i at the respective measuring points from the data file stored in the external memory device 3 (Step S24). Then, a normalization coefficient G′ is calculated from the unsealed common reference value Fmin, the TIC values P.sub.i at the respective measuring points, and the maximum allowable value long_Max of variable, which is a value specific to the apparatus according to the following expression (4),

G′=long_Max×(Fmin/P.sub.i)   (4).

Then, the intensity values of the mass spectra at the respective measuring paints included in the mass analysis imaging data is multiplied by the normalization coefficient G′, so that the intensity values are corrected and stored in the main memory 28 (Step S25). In this case, being different from the description of the normalization processing method described above, the element of long_Max is not included in the unsealed individual reference value F.sub.j, and as a matter of course, the element of long_Max is not included in the unsealed common reference value Fmin. However, since the element of long_Max is included in the normalization coefficient G′ obtained by the expression (4), scaling is performed simultaneously so as not to cause overflow at the time of normalization operation. Therefore, the same normalization as the normalization method described above is achieved.

[0084] By multiplying the respective intensity values with the normalization coefficient G′ as described above, the mass analysis imaging data having all the intensity values normalized based on the unsealed common reference value Fmin is stored in the main memory 28. Processing of Steps S24, S25 are performed for another mass analysis imaging data stored in a different data file, and the respective intensity values are normalized based on the same unsealed common reference value Fmin.

[0085] In the case where comparison of two-dimensional distributions of the ion obtained from the plurality of samples and having a specific mass charge ratio is desired, only data that indicates the intensity values at the specific mass charge ratio out of the respective mass analysis imaging data may be extracted, and the intensity values may be normalized and stored in the main memory 28. In this manner, when the plurality of mass analysis imaging data in which the intensity values are normalized under the same common reference value are prepared in the main memory 28, for example, the statistical analysis 26 performs predetermined statistical analysis based on the data in which the intensity values are normalized, and the result is output to a screen of the display unit 5 via the display processing part 29. When displaying a plurality of imaging images to be compared visually, the imaging image creating part 25 reproduces imaging images based on data stored in the main memory 28 and outputs the image onto the screen of the display unit 5 via the display processing part 29.

[0086] As described above, the peak matrix data are stored respectively in the respective data files stored in the external memory device 3, and if execution of comparison analysis of the plurality of mass analysis imaging data is desired by performing multivariable analysis such as main component analysis or the like using the peak matrix data, the normalization calculation processing part 23 may read out the peak matrix data instead of the mass spectrum data in the data files stored in the external memory device 3, and multiply the intensity values, which is respective elements in the peak matrix data, by the above-described normalization coefficient. By performing multivariable analysis using a plurality of normalized peak matrix data obtained in this manner, influence of variations in condition of measurement when measuring the plurality of samples is alleviated, so that multivariable analysis with high degree of accuracy is achieved.

[0087] As described above, by making the conditions of peak detection more strict and reducing the original size of the peak matrix data, the data quantity of the peak matrix data after normalization to be stored in the main memory 28 may be small. Therefore, an occupied capacity in the main memory 28 may be reduced, and time required for multivariable analysis is also reduced.

[0088] In the imaging mass analyzing device of this example, in parallel to execution of measurement for the sample for collecting data, an individual reference value is calculated and stored the individual reference value in the data file. However, in the case where comparison analysis between mass analysis imaging data acquired by the imaging mass analyzing device of the example and the mass analysis imaging data acquired by another apparatus is desired, the data file creased by the another apparatus may include only the mass spectrum data at the respective measuring points. In other words, the data file in which the individual reference value is not stored may be handled depending on the case.

[0089] Therefore, by using data conversion software or the like installed separately on the computer, for example, data is read in the main memory by processing by the software to detect peaks, a peak matrix is prepared and stored in a data file in which data is stored, and the individual reference value is calculated and is stored in the same data file. By creating the data file in which the peak matrix and the individual reference value are added in this manner, the processing as described in conjunction with FIG. 3 is enabled in the imaging mass analyzing device of the above-described example. As a matter of course, such a function of the data conversion software may be integrated in the data analysis software that performs the normalization processing as described above.

[0090] In the imaging mass analyzing device of the example described above, the common reference value is automatically determined based on the individual reference values determined from the respective mass analysis imaging data. However a room for determination for the analyst for setting the common reference value may also be provided. In order to do so, for example, the common reference value determining part 22 displays the plurality of unsealed individual reference value F.sub.j (or the individual reference value Ps.sub.j) read from the data files on a screen of the display unit 5 via the display processing part 29 at the time of comparison analysis of the plurality of mass analysis imaging data. The analyst confirms these values and inputs, for example, a value smaller than the minimum value therein from the operating unit 4 as the common reference value. At this time, if the input value is larger than the smallest value in the plurality of individual reference values, probability of occurrence of overflow during the normalization operation remains. Therefore, a configuration in which an alarm is issued when such a value is input or an upper limit is set so as to disable such an input is also applicable. The common reference value determining part 22 sets the input value as the common reference value. Accordingly, normalization of the intensity values with reference to the value input by the analyst instead of the value automatically determined is enabled.

[0091] As described above, in the imaging mass analyzing device of the example described above, the unsealed individual reference values F.sub.j in which the element of the long_Max is not included is stored in the data files. However, since long_Max is determined specifically for the apparatus (software), the individual reference value Ps.sub.j in which the element of long_Max is included stay be obtained based on the expression (2) and may be stored in the data file as described in the description or the normalization processing method. In this case, the common reference value Pmin is obtained from the individual reference value Ps.sub.j, the intensity values of the mass spectrum data at the respective measuring points may be normalized by using the normalization coefficient G obtained according to the expression (3) based on the common reference value Pmin, and scaling does not have to be performed at the time of normalization.

[0092] The example described above is an example in which the present invention is applied to processing of data obtained by the imaging mass analyzing device. However, the present invention may be applied not only to the imaging mass analyzing device, but also to various analyzing devices configured to acquire any spectrum information respectively for a number of measuring points set in the regions to be measured on samples.

[0093] For example, in the FTIR imaging device, Fourier transform infrared spectrophotometric measurement is performed on a number of measuring points in regions to be measured on samples, and infrared absorption spectra are obtained respectively. Alternatively, in the Raman spectroscopic imaging device, Raman scattering spectra is performed on a number of measuring points in a region to be measured on a sample, and Raman scattering spectra are acquired respectively. In the electron beam microanalyser or the scanning type electron microscope/energy dispersion type X-ray analyzing device, a number of measuring points in a region to be measured on a sample are irradiated with an electron beam as exciting beam, respectively, and spectra indicating relationships between energy and intensity of the characteristic X-ray emitted from the sample thereby are acquired, respectively. In the scanning-type probe microscope, a scanning probe microscope configured to measure a force acting on a probe while changing a distance between the probe and the sample at a number of measuring points on the sample and acquire spectra which indicates a relationship between a height from the surface of the sample and physical property data such a spectra (force curve) indicating a relationship between a distance between the probe and the sample and a force applied to the probe. Therefore, for example, even when comparison analysis of data obtained by these analyzing devices is performed for the plurality of samples, accurate comparison analysis is enabled by normalization of the intensity values as described above.

[0094] The example described above and various modifications described above are examples of the present invention only, and modifications, corrections, and additions performed within the scope of the present invention are included within the claims of this application, as a matter of course.

REFERENCE SIGNS LIST

[0095] 1 . . . imaging mass analyzing part [0096] 2 . . . data processing part [0097] 20 . . . data collecting part [0098] 21 . . . individual reference value calculating part [0099] 22 . . . common reference value determining part [0100] 23 . . . normalization calculation processing part [0101] 24 . . . peak matrix creating part [0102] 25 . . . imaging image creating part [0103] 26 . . . statistical analysis operating part [0104] 27 . . . data file creating part [0105] 28 . . . main memory [0106] 29 . . . display processing part [0107] 3 . . . external memory device [0108] 4 . . . operating unit [0109] 5 . . . display unit [0110] 100 . . . sample [0111] 101 . . . region to be measured [0112] 102 . . . measuring point