WO2023026365A1 - Dispositif d'estimation de valeur de maillage, procédé d'estimation de valeur de maillage et programme - Google Patents

Dispositif d'estimation de valeur de maillage, procédé d'estimation de valeur de maillage et programme Download PDF

Info

Publication number
WO2023026365A1
WO2023026365A1 PCT/JP2021/031012 JP2021031012W WO2023026365A1 WO 2023026365 A1 WO2023026365 A1 WO 2023026365A1 JP 2021031012 W JP2021031012 W JP 2021031012W WO 2023026365 A1 WO2023026365 A1 WO 2023026365A1
Authority
WO
WIPO (PCT)
Prior art keywords
mesh
value
block
estimated
theoretical variogram
Prior art date
Application number
PCT/JP2021/031012
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 日本電信電話株式会社
Priority to JP2023543524A priority Critical patent/JPWO2023026365A1/ja
Priority to PCT/JP2021/031012 priority patent/WO2023026365A1/fr
Publication of WO2023026365A1 publication Critical patent/WO2023026365A1/fr

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N5/00—Computing arrangements using knowledge-based models
    • G06N5/04—Inference or reasoning models
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04W—WIRELESS COMMUNICATION NETWORKS
    • H04W24/00—Supervisory, monitoring or testing arrangements
    • H04W24/08—Testing, supervising or monitoring using real traffic

Definitions

  • the present invention relates to techniques for estimating mesh values based on block observations.
  • observation point and the estimation point are points at which the observation value and the estimation value are obtained, and are, for example, positional coordinates or a small area centered on the positional coordinates (hereinafter referred to as a mesh).
  • Observed values and estimated values are, for example, the number of users, the amount of underground resources, and the amount of precipitation.
  • Kriging for example, is known as a typical technique related to the spatial interpolation method for estimating a value at an estimation point using observed values obtained at a plurality of observation points existing within a given target area. ing.
  • kriging involves creating an empirical variogram from observations obtained at multiple observation points, fitting a theoretical variogram to the empirical variogram using a constrained nonlinear least-squares method, and using the fitted theoretical variogram. to estimate the value at the estimation point.
  • Non-Patent Document 1 proposes a method using block kriging to create a distribution of indices related to wheat yield.
  • observation values obtained at observation points are used to estimate values at estimation points. or estimating the value of an unobserved block using the observed value of the block, or estimating the value of an unobserved block using the observed value obtained at the observation point.
  • the observed value of the block is the average value of the values distributed within the block. is N/S.
  • the observed value of the block is the average value of the 25 meshes existing in the block. good too.
  • Keshavarzi Ali, Fereydoon Sarmadian, and Abbas Ahmadi. ⁇ Spatially-based model of land suitability analysis using Block Kriging,'' Australian Journal of Crop Science 2011.
  • block kriging can estimate the value at the estimation point using observations obtained from multiple blocks within a certain region of interest, some or all of the blocks overlap with other blocks, resulting in overlapping If the values of the area are observed dispersedly according to the accommodation rate for each of the overlapping blocks, the observed value of the block will be smaller than the value that should be observed, so it cannot always be estimated correctly.
  • a case where a part or all of a block overlaps with another block, and the value of the overlapping area is observed dispersedly according to the accommodation rate for each of the overlapping blocks is, for example, a plurality of eci ( E-UTRAN Cell ID) overlaps partially or entirely, and multiple users existing in the overlapped range are each accommodated in separate eci.
  • the eci coverage range corresponds to a block
  • the number of users existing in the overlapping range corresponds to the value of the overlapping region.
  • the accommodation rate of a certain block is the ratio of the number of users observed in that block to the number of users existing in the area of that block.
  • the eci cover range differs for each eci, and the cover range overlaps with the cover range of other eci existing in the same base station and the cover range of eci existing in surrounding base stations.
  • block kriging can estimate the value of the estimation point using the theoretical variogram fitted to the empirical variogram, it can be estimated correctly if a sufficient number of observations cannot be secured, because a highly reliable empirical variogram cannot be obtained. Not necessarily.
  • a case where a sufficient number of observations cannot be secured is, for example, a case where the eci existing in the target region are set as blocks and the number of users obtained in each eci is set as an observation value.
  • the present invention has been made in view of the above points, and considering the effect of the accommodation rate on overlapping blocks, even under conditions where a sufficient number of observations cannot be obtained, from the observation value of the block, We aim to estimate the mesh values in
  • a theoretical variogram estimating unit that estimates a theoretical variogram model and parameters necessary for estimating mesh values;
  • a mesh value estimating unit that estimates a mesh value using the model and parameters of the theoretical variogram estimated by the theoretical variogram estimating unit, and the accommodation rate of the block with respect to the mesh;
  • an output unit that outputs a final estimated value after repeatedly performing the processing of the theoretical variogram estimating unit and the processing of the mesh value estimating unit until a preset objective function is optimized;
  • a mesh value estimator is provided having:
  • the disclosed technique it is possible to estimate the mesh value in the target area from the block observation value even under the condition that a sufficient number of observations cannot be obtained, considering the influence of the accommodation rate for overlapping blocks. It becomes possible.
  • FIG. 1 is a functional configuration diagram of a mesh value estimating device according to an embodiment of the present invention
  • FIG. 4 is a flowchart for explaining the operation of the mesh value estimating device
  • 4 is a flowchart for explaining the operation of the mesh value estimating device
  • It is a figure which shows the hardware structural example of a mesh value estimation apparatus.
  • the mesh value estimating device 100 it is possible to estimate the mesh value in the target area from the observed value of the block, considering the effect of the accommodation rate for overlapping blocks, even under the condition that sufficient number of observations cannot be obtained.
  • the mesh value estimating device 100 will be described.
  • each base station serves one or more cells, and each user residing in a cell communicates with the base station serving that cell. This communication allows the user to be observed.
  • the ID of the cell is eci.
  • the eci coverage is called a block.
  • the number of users represents the number of users observed as users who communicate via eci.
  • the mesh value estimation device 100 estimates the number of users z(u i ) of the mesh u i from the number of users z(B k ) observed in eciB k which is a block of eci. The configuration and operation of the mesh value estimation device 100 will be described in detail below.
  • FIG. 1 is a diagram showing an example of the functional configuration of a mesh value estimation device 100 according to this embodiment.
  • the mesh value estimation device 100 has an input unit 110, a theoretical variogram estimation unit 120, a mesh value estimation unit 130, and an output unit 140.
  • FIG. 2 is a flow chart showing the operation of the mesh value estimation device 100. As shown in FIG. The operation of each unit constituting the mesh value estimation device 100 will be described along the procedure shown in FIG.
  • ⁇ S101 Input>
  • i 1 , 2 , .
  • B denote the number of meshes and the number of blocks that exist in the region of interest, respectively.
  • the input unit 110 may input the above data from any input source. For example, with respect to the position (x i , y i ) of the mesh u i , the input unit 110 may divide the target region into N [m] square mesh units and input the respective center positions.
  • the observation value z(B k ) of B k may be input by obtaining a log of eci, and the mesh set U k belonging to block B k may be input by obtaining the distribution of meshes in the target region
  • the coverage of block B k may be mapped and the meshes overlapping the coverage may be input as mesh set U k .
  • the data may be input by reading data stored in an auxiliary storage device or the like included in the mesh value estimation device 100 .
  • the theoretical variogram estimating unit 120 estimates the model parameters of the theoretical variogram necessary for estimating the mesh value and the accommodation rate ⁇ i,k of the block B k with respect to the mesh u i .
  • the model and parameters of the theoretical variogram may be estimated from a list of models of the theoretical variogram given in advance and a parameter search space.
  • a model of a theoretical variogram is, for example, a spherical function model, an exponential model, a Gaussian model, a nugget model, a power model, or a Hall effect model, and a parameter is a parameter contained in these models.
  • the estimation of the model parameters and capacity of the theoretical variogram may be manually estimated (set) or estimated using any optimization solver.
  • the objective function to be optimized when estimating using the optimization solver is the observed value z(B k ) of the block B k and the estimated value z(u i ) of the mesh u i ⁇ U k belonging to the block B k .
  • Mean absolute error, mean squared error, square root of mean squared error, mean squared logarithmic error, square root of mean squared logarithmic error, mean absolute percent error, square root of mean squared percent error may be used and penalties on these errors
  • a function with additional terms may also be used.
  • the penalty term may be a term that includes the difference between the empirical variogram calculated from the estimated values and the estimated theoretical variogram.
  • the theoretical variogram model is a spherical function model
  • i,k 1/(
  • the parameters a, b, and c are parameters included in the spherical function model
  • represents the number of blocks covering the mesh u i .
  • the search space of the model may be [spherical function model, exponential model, Gaussian model, nugget model, power model, Hall effect model], and the parameters of each model
  • the search space for a, b, and c may be set to [0, R max ], [0, 10], and [0, z(B) max ], respectively, and the search space for the accommodation rate may be set to [0, 1]. good.
  • R max represents the maximum distance between any two meshes in the target region
  • z(B) max represents the maximum observed value z(B k ).
  • the model and parameters of the theoretical variogram may be estimated using an optimization solver, and the accommodation rate may be manually set (for example, input from the input unit 110).
  • the theoretical variogram estimation unit 120 in the case of manually setting the model parameters and the accommodation rate of the theoretical variogram, for example, inputs the model parameters and the accommodation rate of the theoretical variogram, and passes the input data to the mesh value estimation unit 130. have a function.
  • the theoretical variogram estimating unit 120 in the case of estimating the model parameters and accommodation rate of the theoretical variogram with the optimization solver is equipped with the optimization solver.
  • the mesh value estimator 130 estimates the value of each mesh using the theoretical variogram estimated in S102 and the formula formulated in consideration of the influence of the accommodation rate on overlapping blocks. Details of the estimation method will be described later.
  • the output unit 140 outputs the mesh estimated values finally estimated by the theoretical variogram estimating unit 120 and the mesh value estimating unit 130 .
  • the output unit 140 may output the estimated value of the mesh to an arbitrary output destination.
  • the output unit 140 may output the mesh estimated value to a server device or the like via a communication network, may output the mesh estimated value to a display or the like, or may output the mesh estimated value to an auxiliary storage device or the like. Estimates may be output.
  • the mesh value estimating unit 130 estimates the estimated value ⁇ z(u i ) of the mesh u i using the formula formulated in consideration of the accommodation rate and the theoretical variogram estimated in S102 described above.
  • the estimated value ⁇ z(u i ) is formulated as follows.
  • the second numerator term on the right side is a correction term considering the influence of the accommodation rate, and plays a role of correcting the value of the observed value z(B k ) affected by the accommodation rate.
  • the above correction term includes the mesh true value z(u i ).
  • FIG. 3 is a flow chart showing an example of processing of the mesh value estimation unit 130 according to this embodiment.
  • ⁇ S201 Kriging coefficient calculation>
  • the mesh value estimation unit 130 first calculates a Kriging coefficient.
  • the Kriging coefficients are obtained by solving the following Kriging equations.
  • ⁇ (u i , u j ) represents the inter-mesh variogram. Since the values of ⁇ (u i , u j ) can be calculated from the theoretical variogram estimated in S102, the Kriging coefficients can be calculated using these values.
  • the initial value for example, an arbitrary value may be uniformly given, or a different value may be given for each ⁇ z(u i ).
  • the correction term is calculated for each block B k , a value obtained by dividing the observed value z(B k ) of the block B k by the number of meshes
  • step S205 the mesh value estimation unit 130 determines whether the estimated value ⁇ z calculated in S204 described above satisfies the condition.
  • the condition may be, for example, that the absolute error between ⁇ z and ⁇ z prev is below a threshold or minimum, that the relative error between ⁇ z and ⁇ z prev is below a threshold or minimum, or that ⁇ z and
  • the error of z(B k ) (B k ⁇ B) and the difference between the error of z prev and z(B k ) (B k ⁇ B) may be below a threshold or minimal. If the estimated value ⁇ z does not satisfy the condition, proceed to S206; if the condition is satisfied, proceed to S207.
  • ⁇ S207 Output>
  • the mesh value estimator 130 outputs ⁇ z, which satisfies the condition in S205, as the mesh estimated value for the model parameters and accommodation rate of the estimated theoretical variogram.
  • the mesh value estimation device 100 can be realized by, for example, causing a computer to execute a program.
  • This computer may be a physical computer or a virtual machine on the cloud.
  • the mesh value estimation device 100 can be realized by executing a program corresponding to the processing performed by the mesh value estimation device 100 using hardware resources such as a CPU and memory built into a computer. is.
  • the above program can be recorded in a computer-readable recording medium (portable memory, etc.), saved, or distributed. It is also possible to provide the above program through a network such as the Internet or e-mail.
  • FIG. 4 is a diagram showing a hardware configuration example of the computer.
  • the computer of FIG. 4 has a drive device 1000, an auxiliary storage device 1002, a memory device 1003, a CPU 1004, an interface device 1005, a display device 1006, an input device 1007, an output device 1008, etc., which are interconnected by a bus BS.
  • a program that implements the processing in the computer is provided by a recording medium 1001 such as a CD-ROM or memory card, for example.
  • a recording medium 1001 such as a CD-ROM or memory card
  • the program is installed from the recording medium 1001 to the auxiliary storage device 1002 via the drive device 1000 .
  • the program does not necessarily need to be installed from the recording medium 1001, and may be downloaded from another computer via the network.
  • the auxiliary storage device 1002 stores installed programs, as well as necessary files and data.
  • the memory device 1003 reads and stores the program from the auxiliary storage device 1002 when a program activation instruction is received. Input data is stored in the memory device 1003 or the auxiliary storage device 1002 .
  • the CPU 1004 implements functions related to the mesh value estimation device 100 according to programs stored in the memory device 1003 . For example, the CPU 1004 reads input data stored in the memory device 1003 or the auxiliary storage device 1002, and executes the arithmetic processing described with reference to FIGS.
  • the interface device 1005 is used as an interface for connecting to a network or the like.
  • a display device 1006 displays a GUI (Graphical User Interface) or the like by a program.
  • An input device 1007 is composed of a keyboard, a mouse, buttons, a touch panel, or the like, and is used to input various operational instructions.
  • the output device 1008 outputs the calculation result.
  • the mesh value in the target area can be estimated from the observed values of the blocks, even under the condition that a sufficient number of observations cannot be obtained, considering the influence of the accommodation rate for overlapping blocks. be able to.
  • a correction term that takes into account the storage rate is used in the mesh value estimation formula, and the initial value is given to the correction term and estimation is performed by repeatedly estimating, so the influence of the storage rate on overlapping blocks is considered. Estimates can be made.
  • This specification discloses at least a mesh value estimation device, a mesh value estimation method, and a program for each of the following items.
  • (Section 1) an input unit for inputting the position of the mesh to be estimated, the observation value of the block, and the mesh set belonging to the block;
  • a theoretical variogram estimating unit that estimates a theoretical variogram model and parameters necessary for estimating mesh values;
  • a mesh value estimating unit that estimates a mesh value using the model and parameters of the theoretical variogram estimated by the theoretical variogram estimating unit, and the accommodation rate of the block with respect to the mesh; an output unit that outputs a final estimated value after repeatedly performing the processing of the theoretical variogram estimating unit and the processing of the mesh value estimating unit until a preset objective function is optimized;
  • the objective function preset for repeatedly performing the theoretical variogram estimator and the mesh value estimator is an objective function based on the difference between the block estimate calculated from the mesh estimate and the block observed value, or An objective function with a term
  • (Section 3) an input unit for inputting the position of the mesh to be estimated, the observation value of the block, and the mesh set belonging to the block; While updating the mesh value with the estimated value, the process of calculating the estimated value of the mesh based on the observed value corrected using the block accommodation rate and mesh value for the mesh and the Kriging coefficient calculated from the theoretical variogram is performed.
  • a mesh value estimator having (Section 4) A mesh value estimation method executed by a mesh value estimation device, an input step of inputting the position of the mesh to be estimated, the observed value of the block, and the mesh set belonging to the block; A theoretical variogram estimation step of estimating the model and parameters of the theoretical variogram necessary for estimating the mesh value; a mesh value estimation step of estimating a mesh value using the model and parameters of the theoretical variogram estimated in the theoretical variogram estimation step, and the housing ratio of blocks to the mesh; an output step of outputting a final estimated value after repeatedly performing the theoretical variogram estimation step and the mesh value estimation step until a preset objective function is optimized; A mesh value estimation method with (Section 5) A program for causing a computer to function as the mesh value estimation device according to any one of items 1 to 3.
  • mesh value estimation device 110 input unit 120 theoretical variogram estimation unit 130 mesh value estimation unit 140 output unit 1000 drive device 1001 recording medium 1002 auxiliary storage device 1003 memory device 1004 CPU 1005 interface device 1006 display device 1007 input device 1008 output device

Landscapes

  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Evolutionary Computation (AREA)
  • Artificial Intelligence (AREA)
  • Computational Linguistics (AREA)
  • Data Mining & Analysis (AREA)
  • Signal Processing (AREA)
  • Physics & Mathematics (AREA)
  • Computing Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Software Systems (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

Selon la présente invention, un dispositif d'estimation de valeur de maillage comprend : une unité d'entrée qui entre la position d'un maillage à estimer, une valeur observée d'un bloc, et un ensemble de maillages appartenant au bloc ; une unité d'estimation de variogramme théorique qui estime un modèle et des paramètres pour un variogramme théorique qui est nécessaire pour estimer une valeur de maillage ; une unité d'estimation de valeur de maillage qui estime la valeur de maillage à l'aide du modèle et des paramètres pour le variogramme théorique estimé par l'unité d'estimation de variogrammes théorique ainsi qu'un taux d'adaptation de blocs par rapport au maillage ; et une unité de sortie qui émet une valeur estimée finale après que le processus de l'unité d'estimation de variogramme théorique et le processus de l'unité d'estimation de valeur de maillage aient été mis en œuvre de manière répétée jusqu'à ce qu'une fonction objective prédéfinie soit optimisée.
PCT/JP2021/031012 2021-08-24 2021-08-24 Dispositif d'estimation de valeur de maillage, procédé d'estimation de valeur de maillage et programme WO2023026365A1 (fr)

Priority Applications (2)

Application Number Priority Date Filing Date Title
JP2023543524A JPWO2023026365A1 (fr) 2021-08-24 2021-08-24
PCT/JP2021/031012 WO2023026365A1 (fr) 2021-08-24 2021-08-24 Dispositif d'estimation de valeur de maillage, procédé d'estimation de valeur de maillage et programme

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
PCT/JP2021/031012 WO2023026365A1 (fr) 2021-08-24 2021-08-24 Dispositif d'estimation de valeur de maillage, procédé d'estimation de valeur de maillage et programme

Publications (1)

Publication Number Publication Date
WO2023026365A1 true WO2023026365A1 (fr) 2023-03-02

Family

ID=85321842

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/JP2021/031012 WO2023026365A1 (fr) 2021-08-24 2021-08-24 Dispositif d'estimation de valeur de maillage, procédé d'estimation de valeur de maillage et programme

Country Status (2)

Country Link
JP (1) JPWO2023026365A1 (fr)
WO (1) WO2023026365A1 (fr)

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20080243447A1 (en) * 2007-03-30 2008-10-02 Frederic Roggero Method for Gradually Modifying Lithologic Facies Proportions of a Geological Model
US20110282635A1 (en) * 2010-05-14 2011-11-17 Conocophillips Company Stochastic downscaling algorithm and applications to geological model downscaling
JP2014155105A (ja) * 2013-02-12 2014-08-25 Ntt Docomo Inc 情報処理装置及び情報処理方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20080243447A1 (en) * 2007-03-30 2008-10-02 Frederic Roggero Method for Gradually Modifying Lithologic Facies Proportions of a Geological Model
US20110282635A1 (en) * 2010-05-14 2011-11-17 Conocophillips Company Stochastic downscaling algorithm and applications to geological model downscaling
JP2014155105A (ja) * 2013-02-12 2014-08-25 Ntt Docomo Inc 情報処理装置及び情報処理方法

Also Published As

Publication number Publication date
JPWO2023026365A1 (fr) 2023-03-02

Similar Documents

Publication Publication Date Title
Delignette-Muller et al. fitdistrplus: An R package for fitting distributions
Augustin et al. Modeling spatiotemporal forest health monitoring data
JP6093792B2 (ja) 測位装置、測位方法、測位プログラム、および、測位システム
KR20170042260A (ko) 목표의 실시간 위치 측정과 지도 구축 방법 및 장치
Shiri et al. Data splitting strategies for improving data driven models for reference evapotranspiration estimation among similar stations
US9749873B1 (en) Estimation devices and methods for estimating communication quality of wireless network and method for installing meters thereof
CN109283562A (zh) 一种车联网中车辆三维定位方法及装置
CN104462015B (zh) 处理非高斯Lévy噪声的分数阶线性离散系统状态更新方法
CN103983278B (zh) 一种测量影响卫星姿态确定系统精度的方法
Cao et al. Penalized spline estimation for functional coefficient regression models
JP6159594B2 (ja) 観測値処理装置、観測値処理方法及び観測値処理プログラム
WO2014093234A2 (fr) Calculateur probabiliste de crédits carbone
CN116578660A (zh) 一种基于稀疏数据的电磁目标态势地图构建方法
WO2023026365A1 (fr) Dispositif d'estimation de valeur de maillage, procédé d'estimation de valeur de maillage et programme
Vrugt Multi-criteria optimization using the AMALGAM software package: theory, concepts, and MATLAB implementation
Matarazzo et al. Sensitivity metrics for maximum likelihood system identification
JP5592927B2 (ja) 伝搬特性推定装置、伝搬特性推定方法及び伝搬特性推定プログラム
CN115792796B (zh) 基于相对观测等效模型的协同定位方法、装置及终端
KR100961875B1 (ko) 수직 비행경로 자동생성 방법
Tan et al. Estimation of wildfire wind conditions via perimeter and surface area optimization
JP5792502B2 (ja) 解析装置、その方法及びそのプログラム
CN110728289B (zh) 一种家庭宽带用户的挖掘方法及设备
Silva et al. Integration of a low-cost global navigation satellite system to a single-board computer using Kalman filtering
WO2023026364A1 (fr) Dispositif d'estimation de valeur de maillage, procédé d'estimation de valeur de maillage et programme
CN115307687A (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: 21954980

Country of ref document: EP

Kind code of ref document: A1

WWE Wipo information: entry into national phase

Ref document number: 2023543524

Country of ref document: JP

NENP Non-entry into the national phase

Ref country code: DE

122 Ep: pct application non-entry in european phase

Ref document number: 21954980

Country of ref document: EP

Kind code of ref document: A1