WO2024004281A1 - Procédé de traitement de valeur de mesure et système de traitement de valeur de mesure - Google Patents

Procédé de traitement de valeur de mesure et système de traitement de valeur de mesure Download PDF

Info

Publication number
WO2024004281A1
WO2024004281A1 PCT/JP2023/009054 JP2023009054W WO2024004281A1 WO 2024004281 A1 WO2024004281 A1 WO 2024004281A1 JP 2023009054 W JP2023009054 W JP 2023009054W WO 2024004281 A1 WO2024004281 A1 WO 2024004281A1
Authority
WO
WIPO (PCT)
Prior art keywords
value
measurement
measured
values
reference value
Prior art date
Application number
PCT/JP2023/009054
Other languages
English (en)
Japanese (ja)
Inventor
太夢 井上
Original Assignee
株式会社島津製作所
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by 株式会社島津製作所 filed Critical 株式会社島津製作所
Publication of WO2024004281A1 publication Critical patent/WO2024004281A1/fr

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01F—MEASURING VOLUME, VOLUME FLOW, MASS FLOW OR LIQUID LEVEL; METERING BY VOLUME
    • G01F23/00—Indicating or measuring liquid level or level of fluent solid material, e.g. indicating in terms of volume or indicating by means of an alarm
    • G01F23/22—Indicating or measuring liquid level or level of fluent solid material, e.g. indicating in terms of volume or indicating by means of an alarm by measuring physical variables, other than linear dimensions, pressure or weight, dependent on the level to be measured, e.g. by difference of heat transfer of steam or water
    • G01F23/28—Indicating or measuring liquid level or level of fluent solid material, e.g. indicating in terms of volume or indicating by means of an alarm by measuring physical variables, other than linear dimensions, pressure or weight, dependent on the level to be measured, e.g. by difference of heat transfer of steam or water by measuring the variations of parameters of electromagnetic or acoustic waves applied directly to the liquid or fluent solid material
    • G01F23/296—Acoustic waves
    • G01F23/2962—Measuring transit time of reflected waves
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01F—MEASURING VOLUME, VOLUME FLOW, MASS FLOW OR LIQUID LEVEL; METERING BY VOLUME
    • G01F23/00—Indicating or measuring liquid level or level of fluent solid material, e.g. indicating in terms of volume or indicating by means of an alarm
    • G01F23/80—Arrangements for signal processing

Definitions

  • the present invention relates to a measured value processing method and a measured value processing system.
  • the level of the liquid in a container containing the sample or reagent may be measured to monitor the remaining amount of the sample or reagent.
  • the liquid level height in the container can be measured with a non-contact liquid level sensor such as an ultrasonic sensor.
  • the liquid level sensor itself includes measurement errors, and it is not always possible to obtain an accurate liquid level height with a single measurement. Therefore, it is preferable to measure the liquid level height multiple times and calculate the average value of the measured values to obtain the estimated true value.
  • the liquid level in the container shakes due to vibrations caused by the operation of the motor installed in the device, which can cause some fluctuations. It is known that the measured value becomes an outlier that deviates from the true value, and furthermore, the distribution of the measured value becomes a non-normal distribution. The average value of such non-normally distributed measured values tends to have a large difference from the true value, and there is a risk that the error between the determined average value and the true value may exceed the allowable error for measuring the liquid level. be.
  • the present invention has been made in view of the above problems, and an object of the present invention is to accurately exclude outliers even when the distribution of measured values is non-normal, so that an estimated value close to the true value can be obtained. That is.
  • the measurement value processing method uses N measurement values obtained by performing a large number of N measurements on a measurement object whose measurement value distribution is a non-normal distribution.
  • a measured value processing method for estimating a true value comprising: The measured value processing method includes: a reference value searching step of searching for a reference value R from among the N measured values; Based on the reference value R searched in the reference value searching step, NM measurement values are set as outliers and estimated in descending order of the deviation from the reference value R among the N measurement values.
  • the measured value processing system includes: an information storage unit that stores N measured values obtained by performing a large number of N measurements on a measurement object whose measured value distribution is a non-normal distribution; a calculation unit that performs a process of estimating the true value of the measurement target using the N measurement values stored in the information storage unit; In the estimation process, the calculation unit includes: a reference value searching step of searching for a reference value R from among the N measured values; Based on the reference value R searched in the reference value searching step, NM measurement values are set as outliers and estimated in descending order of the deviation from the reference value R among the N measurement values.
  • the median value or average value of n measured values is determined, and the deviation from the median value or the average value among the n measured values is determined.
  • FIG. 1 is a block diagram showing an example of a measured value processing system. It is a flowchart which shows the processing procedure of the measured value in the same Example. It is a flowchart which shows the search procedure of the reference value R in the same Example. This is an example of the distribution of measurement data. This is data showing the distribution of simple 15-time average values obtained using the measurement data of FIG. 4. This is data showing the average value distribution after removing outliers, which was obtained using the measurement data of FIG. 4.
  • the measured value processing system 2 includes an information storage section 4 and a calculation section 6.
  • the measured value processing system 2 receives data on a large number of measured values regarding the height of the liquid level in the container 10 acquired by the liquid level sensor 8 .
  • the liquid level sensor 8 is, for example, an ultrasonic sensor. Data of a large number of measured values taken into the measured value processing system 2 is stored in the information storage section 4.
  • the calculation unit 6 is configured to perform estimation processing of the true value of the liquid level height in the container 10 using data of a large number of measured values stored in the information storage unit 4.
  • the measured value processing system 2 is a computer device such as a personal computer that is equipped with software for processing measurement data.
  • the information storage unit 4 is a function realized by a part of the storage area of an information storage device installed in the computer device
  • the calculation unit 6 is a function realized by a CPU (central processing unit) installed in the computer device. This is a function achieved by executing a program.
  • a measurement value to be used as the reference value R is searched among the N measurement values of interest (step 101).
  • the reference value R is a value that serves as a reference for selecting M measured values (M ⁇ N), which will be described later, and is a value different from a simple median value or average value.
  • M ⁇ N M measured values
  • the search procedure for the reference value R will be described later.
  • the reference value R is used to select M measured values from the N measured values to be used for calculating the estimated true value (step 102).
  • the selection of the M measured values is performed by excluding NM measured values having large deviations from the reference value R from among the N measured values as outliers from the targets for calculating the estimated value.
  • the ultrasonic sensor 8 measures the distance to the liquid surface based on the time it takes the ultrasonic waves to bounce back from the liquid surface. If the liquid surface to be measured is shaking, the ultrasonic sensor 8 is the most effective ultrasonic sensor. There is a tendency to measure the position close to 8, that is, the highest position of the liquid level. Therefore, the measured value of the liquid level height in the container 10 tends to be higher than the true value.
  • the deviations are the same includes not only the deviations completely matching, but also the deviations being approximate enough to be considered to be matching.
  • FIG. 4 shows measurement data obtained by measuring the liquid level height in the container 10 at 300 ⁇ sec intervals using the ultrasonic sensor 8.
  • the horizontal axis of this data is the number of measurements, and the vertical axis is the measured value (unit: 0.1 mm).
  • the true value of the liquid level height in the container 10 is 65.4 mm. As is clear from this data, almost all of the measured values that deviate from the true value (outliers) greatly exceed the true value.
  • the estimated value is calculated by taking the average of the remaining M measured values (step 103). Note that it is not excluded that M is 1. When M is 1, the reference value R becomes the estimated value of the true value.
  • the n measurement values that are candidates are sorted by the magnitude of the values (step 202), and the median value of the n measurement values is determined (step 203). If n is an even number, two measurement values will be applicable as candidates for the median value. However, as mentioned above, considering that outliers are almost always higher than the true value, two measurement values are considered as candidates for the median value. Among the values, the measured value with the smaller value is taken as the median value. Thereby, a measured value closer to the true value can be set as the median value.
  • the measurement value with the largest deviation from the median value is excluded from the candidates (step 204). At this time, if there are a plurality of measured values with the largest deviation from the median value, the measured value with the larger value is preferentially excluded from the candidates for the same reason as the selection of outliers.
  • the above steps 202 to 204 are repeatedly executed until the number n of measurement values remaining as candidates becomes 1 (step 205), and the last measurement value remaining as a candidate is set as the reference value R (step 206). ).
  • the median value of n measured values is calculated and the measured value with the largest deviation from the median value is excluded from the candidates for the reference value R, but the present invention is limited to this. It's not a thing.
  • the average value of n measured values may be calculated, and the measured value with the largest deviation from the average value may be excluded from the candidates for the reference value R. In that case, step 201 of sorting the n measurements is not necessary.
  • the distribution of estimated values for each time was calculated by taking the average of the measured value of each time and the seven measurements before and after it (15 in total).
  • Figure 6 shows that based on the measurement data in Figure 4, 10 outliers are removed from the measurement value of each time and 7 measurement values before and after it (15 measurements in total), and the remaining 5 measurement values are The distribution of estimated values at each time was determined by taking the average. For example, the estimated true value at the 50th measurement is calculated using the 43rd to 57th measurement values.
  • Embodiments of the measured value processing method and measured value processing system according to the present invention are as follows.
  • An embodiment of the measured value processing method uses N measured values obtained by performing a large number of N measurements on a measurement object whose measured value distribution is a non-normal distribution.
  • a measurement value processing method for estimating the true value of the measurement target comprising: The measured value processing method includes: a reference value searching step of searching for a reference value R from among the N measured values; Based on the reference value R searched in the reference value searching step, NM measurement values are set as outliers and estimated in descending order of the deviation from the reference value R among the N measurement values.
  • the other measured value is the reference value R.
  • the median value of the n measurement values is determined, and the one having the largest deviation from the median value among the n measurement values is determined.
  • the two measured values are excluded and the remaining measured values are used as candidates for the reference value R.
  • the measurement target is the liquid level height in the container
  • the measurement value is obtained by an ultrasonic sensor
  • the reference value searching step is performed by determining the median value or the average value.
  • the measurement target is the liquid level height in the container
  • the measurement value is obtained by an ultrasonic sensor
  • the selection step includes NM measurement values.
  • information that stores N measured values obtained by performing a large number of N measurements on a measurement object whose measured value distribution is a non-normal distribution.
  • storage section a calculation unit that performs a process of estimating the true value of the measurement target using the N measurement values stored in the information storage unit;
  • the calculation unit includes: a reference value searching step of searching for a reference value R from among the N measured values; Based on the reference value R searched in the reference value searching step, NM measurement values are set as outliers and estimated in descending order of the deviation from the reference value R among the N measurement values.
  • the reference value R is configured to be one of the measured values.
  • the calculation unit calculates a median value of the n measurement values in the loop of the reference value searching step, and calculates a median value of the n measurement values from the median value. It is configured to exclude one measurement value with the largest deviation and use the remaining measurement values as candidates for the reference value R.
  • the measurement target is a liquid level height in a container, and the measurement value is obtained by an ultrasonic sensor,
  • the calculation unit excludes one measurement value having the largest deviation from the median value or the average value in the reference value searching step, if there are a plurality of measurement values having the same deviation from the reference value R, is configured to preferentially exclude the measured value having a large value.
  • the measurement target is the liquid level height in the container, and the measurement value is obtained by an ultrasonic sensor, When excluding the NM measured values from the estimation target as the outliers in the selection step, if there are a plurality of measured values having the same deviation from the reference value R, The measurement value having a large value is preferentially excluded from the estimation target.
  • Measured value processing system 4 Information storage section 6
  • Arithmetic section 8 Ultrasonic sensor 10

Landscapes

  • Physics & Mathematics (AREA)
  • Fluid Mechanics (AREA)
  • General Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Signal Processing (AREA)
  • Acoustics & Sound (AREA)
  • Electromagnetism (AREA)
  • Thermal Sciences (AREA)
  • Indication And Recording Devices For Special Purposes And Tariff Metering Devices (AREA)
  • Testing Or Calibration Of Command Recording Devices (AREA)

Abstract

La présente invention comprend : une étape de recherche de valeur de référence pour rechercher une valeur de référence R parmi N valeurs de mesure ; une étape de sélection pour sélectionner M valeurs de mesure en excluant, à partir de la cible d'estimation, N-M valeurs de mesure, en tant que valeurs aberrantes, dans un ordre décroissant d'écart par rapport à la valeur de référence R, parmi les N valeurs de mesure, sur la base de la valeur de référence R recherchée dans l'étape de recherche de valeur de référence ; et une étape de calcul de valeur d'estimation pour définir une valeur moyenne des M valeurs de mesure sélectionnées dans l'étape de sélection en tant que valeur d'estimation d'une valeur réelle. L'étape de recherche de valeur de référence comprend une boucle se répétant d'un état n = N jusqu'à n = 1, un processus de détermination de la valeur médiane ou moyenne des n valeurs de mesure, à l'exclusion d'une valeur de mesure ayant le plus grand écart par rapport à la valeur médiane ou à la valeur moyenne parmi les n valeurs de mesure, et la définition de n-1 valeurs de mesure en tant que candidats pour la valeur de référence R, et après les extrémités de boucle, une dernière valeur de mesure restante lorsque le candidat est défini comme valeur de référence R.
PCT/JP2023/009054 2022-06-28 2023-03-09 Procédé de traitement de valeur de mesure et système de traitement de valeur de mesure WO2024004281A1 (fr)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
JP2022103250 2022-06-28
JP2022-103250 2022-06-28

Publications (1)

Publication Number Publication Date
WO2024004281A1 true WO2024004281A1 (fr) 2024-01-04

Family

ID=89381935

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/JP2023/009054 WO2024004281A1 (fr) 2022-06-28 2023-03-09 Procédé de traitement de valeur de mesure et système de traitement de valeur de mesure

Country Status (1)

Country Link
WO (1) WO2024004281A1 (fr)

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101330030A (zh) * 2007-06-21 2008-12-24 中芯国际集成电路制造(上海)有限公司 检测数据中异常点的去除方法
US20150169706A1 (en) * 2013-12-16 2015-06-18 Unmesh Sreedharan Iterative approach to detect outliers
JP2021523370A (ja) * 2018-05-09 2021-09-02 トリナミクス ゲゼルシャフト ミット ベシュレンクテル ハフツング 少なくとも１つの保管ユニットの充填レベルを決定する方法及び装置

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101330030A (zh) * 2007-06-21 2008-12-24 中芯国际集成电路制造(上海)有限公司 检测数据中异常点的去除方法
US20150169706A1 (en) * 2013-12-16 2015-06-18 Unmesh Sreedharan Iterative approach to detect outliers
JP2021523370A (ja) * 2018-05-09 2021-09-02 トリナミクス ゲゼルシャフト ミット ベシュレンクテル ハフツング 少なくとも１つの保管ユニットの充填レベルを決定する方法及び装置

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
GRUBBS FRANK E.: "Sample Criteria for Testing Outlying Observations", THE ANNALS OF MATHEMATICAL STATISTICS, vol. 21, no. 1, 1 March 1950 (1950-03-01), pages 27 - 58, XP093121966, DOI: 10.1214/aoms/1177729885 *

Similar Documents

Publication Publication Date Title
JP5155533B2 (ja) 補正プログラム、及び測定装置
JP6776319B2 (ja) 体液試料内の一時誤りを検出すること
US10043276B2 (en) Image processing method, image processing apparatus, robot apparatus, program, and recording medium
JP6184344B2 (ja) 観測値処理装置
JPWO2019235614A1 (ja) 関係性分析装置、関係性分析方法およびプログラム
EP1855085B1 (fr) Procédé de correction pour un instrument mesurant un contour et instrument mesurant un contour
WO2024004281A1 (fr) Procédé de traitement de valeur de mesure et système de traitement de valeur de mesure
JP4531685B2 (ja) 形状測定装置、形状測定方法
WO2019235603A1 (fr) Dispositif d'analyse de relation, procédé d'analyse de relation et support d'enregistrement
EP3350604B1 (fr) Mesures d'inertie à haute performance à l'aide d'un réseau redondant de capteurs inertiels peu coûteux
CN112797967A (zh) 一种基于图优化的mems陀螺仪随机漂移误差补偿方法
US20040039473A1 (en) Integrated circuit profile value determination
JP7059050B2 (ja) 応力推定装置、応力推定方法およびプログラム
US9664604B2 (en) Measurement apparatus, measurement method, and method of manufacturing article
CN115878986A (zh) 一种变工况轴承的退化趋势评估方法
CN114396881A (zh) 一种光谱测量分析中快速傅里叶变换拟合的方法及装置
Rossi A form of deviation-based method for coordinate measuring machine sampling optimization in an assessment of roundness
CN113255851A (zh) 一种机器人关节光磁混合编码器的数据融合方法及装置
CN108375337B (zh) 一种机器人及其工艺设备相对位姿的测定方法及测定装置
KR20190128565A (ko) 선형 단면 형상을 갖는 대상물의 표면 측정장치 및 측정방법
CN111090281A (zh) 一种基于改进粒子滤波算法估算移动机器人精确方位的方法和装置
JP2023504001A (ja) 計測ランドスケープに基づく計測最適化のためのシステムおよび方法
US20220146452A1 (en) Detector, detection method, and program
CN114624725B (zh) 目标物体监测方法及装置
CN109425634A (zh) 热膨胀系数测定装置和热膨胀系数测定方法

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 23830736

Country of ref document: EP

Kind code of ref document: A1