Ocean Observation Plan Formulation Support Program, Ocean Observation Plan Formulation Support Device, Ocean Observation Plan Formulation Support Method

Abstract

To perform ocean observation planning support by preliminarily presenting an area from which observation data suitable to improve the accuracy of predicting an ocean environment can be acquired to a user.SOLUTION: An ocean observation planning support program 121 to be executed by an ocean observation planning support device causes an arithmetic processing part to execute segmentation processing for dividing an observation object sea area where ocean observation is performed into a plurality of areas on the basis of a mutual relation between respective areas within the sea area, and display processing of a GUI screen for displaying the plurality of areas obtained by dividing the observation object sea area by the segmentation processing on the screen by an input-output part 110.SELECTED DRAWING: Figure 4

Description

Technical Field

[0001] The present invention relates to a program, an apparatus, and a method for assisting in the planning of ocean observations.

Background Art

[0002] In order to accurately predict the ocean environment, it is necessary to conduct ocean observations to observe physical quantities related to the ocean environment, such as the water temperature and salinity of the ocean, and obtain the observation data. On the other hand, ocean observations are costly and time-consuming. Therefore, in order to efficiently conduct ocean observations, it is required to set an area where it is possible to obtain observation data suitable for improving the prediction accuracy of the ocean environment within the target observation sea area, and determine observation points within that area.

[0003] Regarding the setting of observation points in ocean observations, the technique of Patent Document 1 is known. Patent Document 1 discloses the following technique. A physical correlation coefficient regarding the ocean environment to be forecast is calculated based on a numerical ocean model, and a physical correlation distance is calculated by assuming a Gaussian-type correlation coefficient distribution using this calculated value. On the other hand, the Rossby deformation radius, which is the range affected by disturbances generated in the ocean, is calculated based on characteristics such as the water depth, Coriolis parameter, and gravitational acceleration of the sea area to be measured. The prediction result of the physical quantity by these numerical ocean models or the calculation result of the Rossby deformation radius is used as the measurement interval. Alternatively, the magnitude relationship between the correlation distance and the Rossby deformation radius is compared, and the smaller one is selected. As a result, it becomes possible to give guidelines for the maximum laying interval and the maximum number of laying of observation devices, and it is possible to efficiently obtain the observation data required for data assimilation processing.

Prior Art Documents

Patent Documents

[0004]

Patent Document 1

Summary of the Invention

Problems to be Solved by the Invention

[0005] In the technology of Patent Document 1 described above, although it is possible to give guidelines for the maximum laying interval and the maximum number of laying of observation devices, it is not possible to present to the user in advance the area where it is possible to obtain observation data suitable for improving the prediction accuracy of the ocean environment. Therefore, the user cannot determine an appropriate observation point when planning ocean observations.

[0006] The present invention has been made in view of the above-described problems, and an object thereof is to support the planning of ocean observations by presenting to the user in advance the area where it is possible to obtain observation data suitable for improving the prediction accuracy of the ocean environment.

Means for Solving the Problems

[0007] The ocean observation planning support program according to the present invention is a program for supporting the planning of ocean observations, and causes a computer to execute a first process of dividing the sea area where the ocean observations are to be performed into a plurality of regions based on the correlation relationship between each point in the sea area, and a second process of displaying the plurality of regions on a screen. , the first process includes a process of obtaining p observation points where the analysis error covariance is minimized within the sea area, a process of obtaining the correlation degree between each of the p observation points and each location within the sea area, and a process of grouping each location within the sea area for each of the observation points based on the correlation degree. Based on the result of the grouping, the sea area is divided into the plurality of regions. In the process of obtaining the p observation points, a process of creating an observation matrix for each of a plurality of grid points set at predetermined intervals for the sea space within the sea area, a process of obtaining the analysis error covariance for each of the plurality of grid points based on the observation matrix, and a process of selecting, as the observation points, the grid points among the plurality of grid points where the analysis error covariance is minimized are repeatedly performed a plurality of times to obtain the p observation points 。 The ocean observation planning support device according to the present invention is a device for supporting the planning of ocean observations, and includes an arithmetic processing unit that performs arithmetic processing, and an input / output unit that receives a user's operation input and presents the result of the arithmetic processing performed by the arithmetic processing unit to the user. The arithmetic processing unit divides the sea area where the ocean observations are to be performed into a plurality of regions based on the correlation relationship between each point in the sea area, and the input / output unit displays the plurality of regions on a screen and presents them to the user. and the arithmetic processing unit executes a process of obtaining p observation points where the analysis error covariance is minimized within the sea area, a process of obtaining the correlation degree between each of the p observation points and each location within the sea area, and a process of grouping each location within the sea area for each of the observation points based on the correlation degree. Based on the result of the grouping, the sea area is divided into the plurality of regions. In the process of obtaining the p observation points, a process of creating an observation matrix for each of a plurality of grid points set at predetermined intervals for the sea space within the sea area, a process of obtaining the analysis error covariance for each of the plurality of grid points based on the observation matrix, and a process of selecting, as the observation points, the grid points among the plurality of grid points where the analysis error covariance is minimized are repeatedly performed a plurality of times to obtain the p observation points 。 The ocean observation planning support method according to the present invention is a method for supporting the planning of ocean observations using a computer, and the computer divides the sea area where the ocean observations are to be performed a process of obtaining p observation points where the analysis error covariance is minimized within the sea area, a process of obtaining the correlation degree between each of the p observation points and each location within the sea area, a process of grouping each location within the sea area for each of the observation points based on the correlation degree, and based on the result of the grouping, the sea area into a plurality of regions a process of displays the plurality of regions on a screen executes a process of obtaining the p observation points. In the process of obtaining the p observation points, a process of creating an observation matrix for each of a plurality of grid points set at predetermined intervals for the sea space within the sea area, a process of obtaining the analysis error covariance for each of the plurality of grid points based on the observation matrix, and a process of selecting, as the observation points, the grid points among the plurality of grid points where the analysis error covariance is minimized are repeatedly performed a plurality of times to obtain the p observation points 。

Effects of the Invention

[0008] According to the present invention, it is possible to support the planning of ocean observations by presenting to the user in advance an area where observation data suitable for improving the prediction accuracy of the ocean environment can be obtained.

Brief Description of the Drawings

[0009] [Fig. 1] It is a block diagram showing a schematic configuration of an ocean observation plan formulation support device according to an embodiment of the present invention.

Figure 2

Figure 3

Figure 4

Figure 5

Figure 6

Modes for Carrying Out the Invention

[0010] Hereinafter, embodiments of the present invention will be described based on the accompanying drawings.

[0011] FIG. 1 is a block diagram showing a schematic configuration of an ocean observation plan formulation support device according to an embodiment of the present invention. The ocean observation plan formulation support device 10 shown in FIG. 1 is a device that supports the user in formulating an ocean observation plan, and includes an arithmetic processing unit 100, an input / output unit 110, and a data storage unit 140.

[0012] The arithmetic processing unit 100 is configured using a computer such as a PC (Personal Computer) or a server, and includes a memory 120 and a CPU (Central Processing Unit) 130. An ocean observation plan formulation support program 121 is deployed in the memory 120. This ocean observation plan formulation support program 121 is composed of instructions that cause the CPU 130 to execute various processes. The CPU 130 performs various arithmetic processes, receives signals from the input / output unit 110, outputs signals to the input / output unit 110, stores and reads data to / from the data storage unit 140, etc., according to the instructions of the ocean observation plan formulation support program 121. In the following description, the processes executed by the ocean observation plan formulation support program 121 are actually executed by the CPU 130 according to the instructions described in the ocean observation plan formulation support program 121.

[0013] The input / output unit 110 is composed of an input unit including a keyboard 111 and a mouse 112, etc., and an output unit including a display 113, etc. The user can input instructions to the arithmetic processing unit 100 by operating the keyboard 111 or the mouse 112, and visually confirm the results of the arithmetic processes performed by the arithmetic processing unit 100 through the screen display performed by the display 113. That is, the input / output unit 110 receives the user's operation input by the keyboard 111 or the mouse 112, and presents the execution results of the ocean observation plan formulation support program 121 by the arithmetic processing unit 100 to the user by the display 113.

[0014] The data storage unit 140 is, for example, an HDD (Hard Disk Drive) or an SSD (Solid State It is configured using a storage device such as a hard disk drive (HDD) or a solid state drive (SSD), and stores various data used when the arithmetic processing unit 100 executes the ocean observation plan formulation support program 121. The data stored in the data storage unit 140 includes, for example, ocean reanalysis data 141, background error covariance data 142, analysis error covariance data 143, segment data 144, map data 145, and the like. The data storage unit 140 can store these data using a data format such as NetCDF (Network Common Data Form). Note that the ocean observation plan formulation support program 121 may be stored in the data storage unit 140, and when executed by the arithmetic processing unit 100, the CPU 130 may read the ocean observation plan formulation support program 121 from the data storage unit 140 and expand it in the memory 120.

[0015] The ocean reanalysis data 141 is two-dimensional or three-dimensional time-series data obtained using ocean simulations regarding the state of the sea area (observation target sea area) to be observed. Ocean simulation means that in the observation target sea area, grids are arranged at predetermined intervals in the sea space expressed by three axes of latitude, longitude, and depth, and each grid point is associated with each point in the sea area, and the numerical model is used to obtain the time variation of various physical quantities such as sea surface height, water temperature, salinity, and flow velocity at each grid point. In the ocean observation plan formulation support apparatus 10 of the present embodiment, data configured by arranging the values at each grid point of each physical quantity obtained by ocean simulation in time series order is stored in the data storage unit 140 as the ocean reanalysis data 141. Note that, as the ocean reanalysis data 141, data created and published by ocean research institutions or the like may be used.

[0016] In the present embodiment, as an example of the ocean reanalysis data 141, a case where the number of grid points is n and the physical quantity targeted for ocean simulation is the water temperature at a depth of 50 m will be described. In this case, assuming that the number of types of physical quantities handled in the ocean observation plan formulation support program 121 is k, then k = 1. Note that each grid point can be uniquely identified by a grid number from 0 to n - 1.

[0017] The background error covariance data 142 is data representing the covariance of the background error of the ocean simulation. The value of the background error covariance data 142 can be obtained in advance from the ocean reanalysis data 141.

[0018] Here, the background error covariance data 142 can be represented by the background error covariance matrix B shown in the following formula (1).

[0019]

Number

[0020] The matrix A in formula (1) can be represented by the following formula (2).

[0021]

Number

[0022] Each element of the matrix A in formula (2) is a temperature anomaly value obtained from the temperature data represented by the ocean reanalysis data 141. For example, δx y0 can be represented by the following formula (3). The same applies to other elements.

[0023]

Number

[0024] In formula (3), x y0 is a state vector representing the water temperature at each grid point in the y0 - th year (1 ≤ y0 ≤ m), and x - represents the average value of the water temperature at each grid point in the y0 - th year. These values can be obtained from the time - series data of the water temperature at each grid point recorded in the ocean reanalysis data 141 over m years.

[0025] In this embodiment, by using the water temperature data over m years as the ocean reanalysis data 141, the ocean reanalysis data 141 is regarded as m ensembles. Also, in this embodiment, since the monthly average value of the water temperature data is handled, the ocean reanalysis data 141 is set to the monthly average value. Here, the state vector x in Equation (3) is n rows by 1 column, the matrix B is n rows by n columns, and the matrix A is n rows by m columns. Note that in order to reduce the amount of data stored in the data storage unit 140, the matrix A may be stored instead of the matrix B as the background error covariance data 142.

[0026] The analysis error covariance data 143 is data representing the covariance of the analysis errors of the ocean simulation. Here, the optimal observation matrix H opt can be expressed by the following Equation (4).

[0027]

Equation

[0028] Equation (4) indicates that the observation matrix H that minimizes the sum of the diagonals (trace) of the analysis error covariance matrix P a is the optimal observation matrix H opt . Here, ||P a || in Equation (4) is the sum of the diagonals of the analysis error covariance matrix P a and is expressed by the following Equation (5).

[0029]

Equation

[0030] The analysis error covariance data 143 can be represented by the analysis error covariance matrix P a shown in the following Equation (6).

[0031]

Equation

[0032] The weight matrix W in Equation (6) is represented by the following Equation (7). In Equation (7), the observation error covariance matrix R represents the error covariance of the observation matrix H.

[0033] [Number]

[0034] Here, if the number of observation points represented by the optimal observation matrix H opt is p, the observation matrix H is p rows by n columns, the observation error covariance matrix R is p rows by p columns, the weight matrix W is n rows by p columns, and the analysis error covariance matrix P a is n rows by n columns. Note that the number of observation points p can be set to any number according to the user's instruction. Also, in order to reduce the amount of data stored in the data storage unit 140, instead of the matrix P a as the analysis error covariance data 143, the sum of the diagonals of the analysis error covariance matrix P a represented by Equation (5) may be stored.

[0035] The segment data 144 is data representing segments corresponding to each of the p observation points represented by the optimal observation matrix H opt . A segment corresponds to each region when the observation target sea area is divided into a plurality of regions by grouping each grid point in the sea area for each observation point. The segment data 144 is data configured to include information such as the segment ID of each segment, the latitude and longitude of the observation points in each segment, the evaluation value of each segment, and the grid points included in each segment. The structure of this segment data 144 will be described later with reference to FIG. 2.

[0036] The map data 145 is data representing a map of the sea area to be observed, and includes, for example, a map image. It is preferable that this map image includes information such as the positions of islands and land existing within or adjacent to the sea area, and the seabed topography of the sea area. In the marine observation plan formulation support device 10 of the present embodiment, a map of the sea area is displayed on the display 113 based on the map data 145, and the execution result of the marine observation plan formulation support program 121 is displayed on the map. Thereby, the execution result can be presented to the user in an easy-to-understand manner. Note that, for the map data 145, a predetermined data format such as WMTS (Web Map Tile Service) can be used.

[0037] FIG. 2 is a diagram showing the data structure of the segment data 144. The segment data 144 is data representing segments corresponding to each of the p observation points as described above, and has the same number p of records as the number of segments. As shown in FIG. 2, each record of the segment data 144 includes information on the segment ID 201, the latitude and longitude 202 of the observation point, the evaluation value 203, and the grid number 204 included in the segment.

[0038] The segment ID 201 stores an ID number for uniquely identifying each segment. The ID numbers of each segment are set by assigning integers starting from, for example, 0 in record order.

[0039] The latitude and longitude 202 of the observation point stores the latitude and longitude of the observation point in each segment. The latitude and longitude of the observation point are determined by an instruction from the user described later or by a segmentation process described later. Note that in the record corresponding to a segment where these processes have not yet been executed and thus the observation point is undetermined, no information is stored in the latitude and longitude 202 of the observation point and it remains blank.

[0040] The evaluation value 203 stores the evaluation values of the observation points in each segment. The evaluation value is an index indicating how suitable the determined observation point for each segment is for observing the physical quantity to be observed (water temperature in this embodiment), and the higher the evaluation value, the more desirable the observation point is. Note that the evaluation value of each observation point is obtained by the segmentation process described later. In the record corresponding to the segment where the segmentation process has not been executed and thus the evaluation value of the observation point has not been obtained, the evaluation value 203 has no information stored and is blank.

[0041] The grid number 204 included in the segment stores the grid numbers representing the grid points included in the area of each segment. This consists of one or more grid numbers.

[0042] The segment data 144 is composed by setting each of the above information for p records. Note that the segment data 144 may also be composed including other information.

[0043] FIG. 3 is a diagram showing an example of a GUI (Graphical User Interface) screen in which the execution result of the ocean observation plan formulation support program 121 is shown on the map of the observation target sea area. The GUI screen 300 shown in FIG. 3 is presented to the user by being displayed on the display 113, and is used to receive instructions from the user through operation inputs of the keyboard 111 and the mouse 112. This GUI screen 300 includes a segment list 310 and a map display section 320.

[0044] The segment list 310 is a list of each segment determined by executing the ocean observation plan formulation support program 121, and includes display columns for the segment ID 311, the latitude and longitude of the observation point 312, and the evaluation value 313. The corresponding data values of each record of the segment data 144 are respectively displayed in these display columns.

[0045] The user can instruct the Ocean Observation Plan Formulation Support Device 10 on the number p of observation points (the number of records in the segment data 144) by selecting the setting icon 314 in the segment list 310 and specifying an arbitrary numerical value through operation input using the keyboard 111 or the mouse 112. The Ocean Observation Plan Formulation Support Device 10 performs segmentation processing of setting p observation points within the observation target sea area and associating each grid point with one of the observation points by executing the Ocean Observation Plan Formulation Support Program 121 according to the number p instructed by the user, and sets segments for each observation point within the observation target sea area. Thereby, the observation target sea area can be grouped by each observation point and divided into a plurality of regions. Note that when there is no instruction of the number p of observation points from the user, the number p of observation points may be set using the default value p0 and the Ocean Observation Plan Formulation Support Program 121 may be executed.

[0046] A map of the ocean including the observation target sea area is displayed on the map display unit 320. The icon 321 on the map in this map display unit 320 indicates the position of the observation point in the segment designated by the user among the segments shown in the segment list 310. Further, in the map display unit 320, in order to present to the user in an easy-to-understand manner which region on the map each segment shown in the segment list 310 corresponds to, the map is divided and displayed in units of segments. That is, in the map displayed on the map display unit 320, a state in which the observation target sea area is divided into a plurality of regions by the segmentation processing is shown.

[0047] Among the regions of each segment shown in the map display unit 320, in the region 322 of the segment designated by the user, the positions of each grid point existing within the segment are drawn. The position of each grid point within the segment can be obtained by referring to the grid number 204 of the record corresponding to the segment in the segment data 144. Specifically, by referring to the grid number 204 of the record corresponding to the segment, the grid number of each grid point can be obtained, and the position of each grid point can be obtained from the grid number.

[0048] The user can specify any position within each segment on the map displayed on the map display unit 320 through operation input using the keyboard 111 or the mouse 112, and set that position as the observation point of the segment. Note that by setting the observation point before executing the segmentation process, any position on the map can be set as the observation point, and the segment including that observation point can be determined by the segmentation process and displayed on the map.

[0049] In the GUI screen 300 of FIG. 3, the state when the user selects the segment with the segment ID value of 0 in the segment list 310 is shown. In the map display unit 320 of this GUI screen 300, the area 322 of the segment is highlighted, and the position of the observation point within the area 322 is indicated by the icon 321. Further, the latitude and longitude 312 of the observation point in the segment in the segment list 310 are shown to be 25.5 degrees and 125.0 degrees, respectively, which represent the position of the observation point indicated by the icon 321. Note that instead of selecting any segment in the segment list 310, the user can also select any segment on the map in the map display unit 320. In that case, the corresponding part of the selected segment is highlighted in the segment list 310.

[0050] When the marine observation plan formulation support program 121 is executed in the arithmetic processing unit 100, the marine observation plan formulation support device 10 of the present embodiment displays a GUI screen 300 as shown in FIG. 3, and allows the user to specify the number p of observation points through the segment list 310 of this GUI screen 300. Also, in the map display unit 320 of the GUI screen 300, a map showing the sea area to be observed is displayed, and the user is allowed to specify observation points on this map. At this time, the user can specify the number p of observation points and the position of each observation point by operating the keyboard 111 or the mouse 112 of the input / output unit 110. These pieces of information specified by the user are received by the input / output unit 110 and output to the arithmetic processing unit 100, and are reflected in the segmentation process performed by the marine observation plan formulation support program 121.

[0051] Figure 4 is a sequence diagram showing the overall flow of the process executed by the ocean observation plan formulation support apparatus 10. The process shown in Figure 4 is performed in the ocean observation plan formulation support apparatus 10 by the CPU 130 expanding and executing the ocean observation plan formulation support program 121 in the memory 120.

[0052] In step 401, the ocean observation plan formulation support program 121 reads the segment data 144 and the map data 145 from the data storage unit 140.

[0053] In step 402, the ocean observation plan formulation support program 121 displays the GUI screen 300 on the display 113 of the input / output unit 110 based on the segment data 144 and the map data 145 read in step 401. Here, the segment list 310 is displayed based on the segment data 144, and a map of the ocean including the observation target sea area is displayed on the map display unit 320 based on the map data 145. On the map, the range of each segment indicated by the segment list 310 and the positions of the observation points of the selected segment are displayed. As a result, the GUI screen 300 as shown in Figure 3 is displayed. However, when the ocean observation plan formulation support program 121 is executed for the first time, since the segment data 144 does not exist, an empty list is displayed in the segment list 310, and a map of the ocean including the observation target sea area is displayed on the map display unit 320.

[0054] In step 403, the ocean observation plan formulation support program 121 receives an input operation from the user instructing the number p of observation points and the positions of each observation point from the input / output unit 110. However, it is not necessary to specify the positions of all the observation points, and it is not necessary to specify any of them. Also, for the number p of observation points, if there is no instruction from the user, it may be set to the default value p0. By the process of this step 403, an integer of 1 or more is set as the number p of observation points, and for the observation points whose positions are specified by the user, their positions are set.

[0055] In step 404, the ocean observation plan formulation support program 121 performs segmentation processing on the observation target sea area. This segmentation processing is composed of an observation point search process for searching the positions of observation points where no position is specified by the user, and a grouping process for grouping each grid point in the observation target sea area for each observation point to set segments. Details of these processes will be described later with reference to FIGS. 5 and 6 respectively.

[0056] In step 405, the ocean observation plan formulation support program 121 creates segment data 144 based on the result of the segmentation processing in step 404, and outputs and stores it in the data storage unit 140.

[0057] After the execution of the process in step 405, the process returns to step 401. That is, the ocean observation plan formulation support program 121 reads the segment data 144 created by the segmentation processing in step 404 and stored in step 405 from the data storage unit 140 in the next process. Thereby, the GUI screen 300 on which the positions of each segment and the observation points obtained by the previous segmentation processing are shown on the map is displayed on the display 113.

[0058] FIG. 5 is a flowchart showing the flow of the observation point search process among the segmentation processes executed in step 404 of FIG. 4.

[0059] In step 501, the ocean observation plan formulation support program 121 performs a loop process of repeating each process from step 502 to step 510 p times. Here, p is the number of observation points instructed by the input operation from the user received in step 403 of FIG. 4. Hereinafter, the index of this loop process will be described as j.

[0060] In step 502, the Ocean Observation Plan Formulation Support Program 121 determines whether the latitude and longitude of the observation point are specified for the segment corresponding to the j-th record in the segment list 310. If the latitude and longitude information corresponding to the position of the observation point specified by the user is stored in the latitude and longitude 312 of the observation point in the j-th record, the process proceeds to step 503. Otherwise, the process proceeds to step 506.

[0061] In step 503, the Ocean Observation Plan Formulation Support Program 121 creates an observation matrix H for the observation point of the j-th record. Here, the grid number corresponding to the observation point is specified from the specified latitude and longitude of the observation point, and the grid number is substituted into the variable I. Then, the value of the element corresponding to the variable I in the observation matrix H is set to 1, and the other elements are set to 0 to create the observation matrix H.

[0062] In step 504, the Ocean Observation Plan Formulation Support Program 121 calculates the diagonal sum of the analysis error covariance matrix P a according to the above formulas (5) to (7) based on the observation matrix H created in step 503, and substitutes the value into the variable amin. After obtaining the value of the variable amin in this way, the process proceeds to step 510.

[0063] In step 505, the Ocean Observation Plan Formulation Support Program 121 performs a loop process of repeating each process from step 506 to step 507 n times. Here, n is the number of grid points. Hereinafter, the index of this loop process will be described as i.

[0064] In step 506, the Ocean Observation Plan Formulation Support Program 121 creates an observation matrix H when the grid point with the grid number i is used as the observation point. Here, the value of the element corresponding to the grid number i in the observation matrix H is set to 1, and the other elements are set to 0 to create the observation matrix H.

[0065] In step 507, based on the observation matrix H created in step 506, the background error covariance data 142, and the analysis error covariance data 143, the ocean observation plan formulation support program 121 calculates the diagonal sum of the analysis error covariance matrix P according to the above-mentioned formulas (5) to (7), and substitutes the value into the i-th element of the array a. a and substitutes the value into the i-th element of the array a.

[0066] In step 508, if the value of i is less than n, the ocean observation plan formulation support program 121 returns to step 505, adds 1 to i, and repeats the loop process of steps 506 to 507. When the value of i reaches n, the loop process of steps 506 to 507 ends and the program proceeds to step 509.

[0067] In step 509, the ocean observation plan formulation support program 121 selects the minimum value among the elements of the array a obtained in the loop process of steps 506 to 507, substitutes the value into the variable amin, and substitutes the value of the loop process index i when the value was obtained in step 507 into the variable I. After obtaining the values of the variable amin and the variable I in this way, the program proceeds to step 510.

[0068] In the observation point search process shown in the flowchart of FIG. 5, the loop process of steps 506 to 507 described above is executed, and then the process of step 509 is executed. As a result, the value of the variable I representing the j-th row element of the optimal observation matrix H shown in the above-mentioned formula (4), and the variable amin representing the diagonal sum of the corresponding analysis error covariance matrix P opt are obtained. a are obtained.

[0069] In step 510, based on the values of the variable amin and the variable I set in step 504 or 509, the ocean observation plan formulation support program 121 creates the j-th record of the segment data 144. Specifically, the value of j is set for the segment ID 201 of the record, and the latitude and longitude 202 of the observation point are set to the latitude and longitude of the grid point whose grid number is the value of the variable I. Also, for the evaluation value 203, the difference between amin in the previous loop process and amin in the current loop process is set.

[0070] In step 511, if the value of j is less than p, the ocean observation plan formulation support program 121 returns to step 501, adds 1 to j, and repeats the loop process of steps 502 to 510. When the value of j reaches p, the loop process of steps 502 to 510 ends, the observation point search process ends, and the program proceeds to the next grouping process.

[0071] By executing the observation point search process described above, the ocean observation plan formulation support program 121 can obtain p observation points at which the analysis error covariance matrix P a becomes the minimum.

[0072] Figure 6 is a flowchart showing the flow of the grouping process in the segmentation process executed in step 404 of Figure 4.

[0073] In step 601, the ocean observation plan formulation support program 121 performs a loop process of repeating each process from step 602 to step 606 n times. Here, n is the number of grid points. Hereinafter, the index of this loop process will be described as i.

[0074] In step 602, the ocean observation plan formulation support program 121 performs a loop process of repeating the process of step 603 p times. Here, p is the number of observation points instructed by the input operation from the user received in step 403 of Figure 4. Hereinafter, the index of this loop process will be described as j.

[0075] In step 603, based on the background error covariance data 142, the ocean observation plan formulation support program 121 obtains the background error correlation coefficient between the j-th observation point and the grid point with grid number i, and substitutes the value into the j-th element of the array a.

[0076] In step 604, if the value of j is less than p, the ocean observation plan formulation support program 121 returns to step 602, adds 1 to j, and repeats the loop process of step 603. When the value of j reaches p, the loop process of step 603 ends and the program proceeds to step 605.

[0077] In the grouping process shown in the flowchart of FIG. 6, by executing the loop process of step 603 described above, the correlation degree between each of the p observation points obtained in the observation point search process and each grid point is obtained.

[0078] In step 605, the ocean observation plan formulation support program 121 substitutes the value of the loop process index j when the maximum value is obtained among the elements of the array a into the variable J. When the value of the variable J is obtained in this way, the program proceeds to step 606.

[0079] In step 606, for the segment corresponding to the value of the variable J set in step 605, the ocean observation plan formulation support program 121 adds the value of the index i in the current loop process as the grid number of the grid points included in the segment. Specifically, in the segment data 144, the record with the segment ID 201 having the same value as the variable J is identified, and the value of the index i is added to the grid number 204 included in the segment of the record.

[0080] In step 607, if the value of i is less than n, the ocean observation plan formulation support program 121 returns to step 601, adds 1 to i, and repeats the loop process of steps 602 to 606. When the value of i reaches n, the loop process of steps 602 to 606 ends, the grouping process ends, and the program proceeds to step 405 in FIG. 4.

[0081] By executing the grouping process described above, the Ocean Observation Plan Formulation Support Program 121 determines the degree of correlation between each of the p observation points and each grid point in the observation target sea area, and based on this degree of correlation, each grid point in the observation target sea area can be grouped by associating it with any one of the segments set for each observation point. As a result, the area of each segment can be specified from the distribution range of each grid point in the group for each segment, and the observation target sea area can be divided into a plurality of areas.

[0082] As described above, in the segmentation process in step 404 of FIG. 4, the observation point search process and the grouping process described in FIGS. 5 and 6 are respectively executed. The result of this segmentation process is stored in the data storage unit 140 as segment data 144 in step 405. Also, in steps 401 and 402 that are subsequently executed, it is read from the data storage unit 140 and used for the display of the segment list 310 and the map display unit 320 on the GUI screen 300.

[0083] According to one embodiment of the present invention described above, the following operational effects are achieved.

[0084] (1) The Ocean Observation Plan Formulation Support Program 121, which is a program for supporting the formulation of an ocean observation plan, performs a first process (step 404: segmentation process) of dividing an observation target sea area, which is a sea area where ocean observations are to be carried out, into a plurality of areas (segments) based on the correlation relationship between each point in the sea area, and a second process (step 402: display process of the GUI screen 300) of displaying the plurality of areas obtained by dividing the observation target sea area in the first process on the screen, and causes the arithmetic processing unit 100 to execute them. By doing so, it is possible to present to the user in advance the areas where observation data suitable for improving the prediction accuracy of the ocean environment can be obtained, and support the formulation of an ocean observation plan.

[0085] (2) The segmentation process is the analysis error covariance matrix P within the observation target sea area aThe process of obtaining p observation points where it is minimized (Fig. 5: Observation point search process), the process of obtaining the background error correlation coefficient representing the correlation degree between each of the p observation points and each point in the observation target sea area (step 603), and the process of grouping each point in the observation target sea area for each observation point based on the obtained background error correlation coefficient (Fig. 6: Grouping process). The marine observation plan formulation support program 121 divides the observation target sea area into a plurality of segments based on the grouping result by this grouping process. By doing so, the segmentation process based on the correlation relationship between each point in the observation target sea area can be realized.

[0086] (3) In step 403, the marine observation plan formulation support program 121 causes the arithmetic processing unit 100 to execute a process of having the user specify an observation point on a map showing the observation target sea area and a process of having the user specify the number p of observation points. By doing so, it becomes possible to set the position of the observation point and the number p according to the user's instruction.

[0087] (4) The marine observation plan formulation support device 10 that supports the formulation of a marine observation plan includes an arithmetic processing unit 100 that performs arithmetic processing, and an input / output unit 110 (keyboard 111, mouse 112, and display 113) that receives the user's operation input and presents the result of the arithmetic processing performed by the arithmetic processing unit 100 to the user. The arithmetic processing unit 100 divides the observation target sea area, which is the sea area where marine observations are to be made, into a plurality of regions (segments) based on the correlation relationship between each point in the sea area (step 404). The input / output unit 110 presents the plurality of regions obtained by dividing the observation target sea area by the arithmetic processing unit 100 on the display 113 to the user by screen display (step 402). By doing so, it is possible to present to the user in advance the areas where observation data suitable for improving the prediction accuracy of the marine environment can be obtained, and support the formulation of a marine observation plan.

[0088] (5) The arithmetic processing unit 100 is the analysis error covariance matrix P within the observation target sea area aThe process of obtaining p observation points where it is minimized (Figure 5: Observation Point Search Process), the process of obtaining background error correlation coefficients representing the correlation degree between each of the p observation points and each point in the sea area to be observed (Step 603), and the process of grouping each point in the sea area to be observed for each observation point based on the obtained background error correlation coefficients (Figure 6: Grouping Process) are executed. Based on the grouping result by this grouping process, the ocean observation plan formulation support device 10 divides the sea area to be observed into a plurality of segments. By doing so, it is possible to realize segment division based on the correlation relationship between each point in the sea area to be observed.

[0089] (6) In step 403, the arithmetic processing unit 100 can set the position specified on the map indicating the sea area to be observed as an observation point by the input operation of the user received by the keyboard 111 or the mouse 112, or set the number specified by the input operation of the user as the number p of observation points. By doing so, it is possible to set the position and the number p of observation points according to the instruction of the user.

[0090] Note that the above-described embodiments and various modifications are merely examples, and the present invention is not limited to these contents as long as the features of the invention are not impaired. Also, although various embodiments and modifications have been described above, the present invention is not limited to these contents. Other aspects conceivable within the scope of the technical idea of the present invention are also included in the scope of the present invention.

Explanation of Reference Numerals

[0091] 10 ··· Ocean Observation Plan Formulation Support Device 100 ··· Arithmetic Processing Unit 110 ··· Input / Output Unit 111 ··· Keyboard 112 ··· Mouse 113 ··· Display 120 ··· Memory 121 ··· Ocean Observation Plan Formulation Support Program 130 ··· CPU 140 ··· Data Storage Unit 141···Ocean reanalysis data 142···Background error covariance data 143···Analysis error covariance data 144···Segment data 145···Map data 300···GUI screen 310···Segment list 320···Map display section

Claims

1. A program for assisting in the planning of ocean observations, comprising: a first process of dividing a sea area where the ocean observations are to be performed into a plurality of regions based on the correlation relationship between each point in the sea area; a second process of displaying the plurality of regions on a screen, and causing a computer to execute the processes; The first process includes: a process of obtaining p observation points where the analysis error covariance is minimized within the sea area; a process of obtaining the degree of correlation between each of the p observation points and each point in the sea area; a process of grouping each point in the sea area for each of the observation points based on the degree of correlation; dividing the sea area into the plurality of regions based on the result of the grouping; In the process of obtaining the p observation points: a process of creating an observation matrix for each of a plurality of grid points set at predetermined intervals in the airspace within the sea area; a process of obtaining the analysis error covariance for each of the plurality of grid points based on the observation matrix; a process of repeatedly selecting, a plurality of times, the grid point with the minimum analysis error covariance among the plurality of grid points as the observation point to obtain the p observation points, an ocean observation plan formulation support program.

2. In the ocean observation plan formulation support program according to Claim 1, causing the computer to execute a process of allowing a user to specify any position on a map showing the sea area; In the process of obtaining the p observation points, including, among the plurality of grid points, the grid point corresponding to the position specified by the user as one of the p observation points, an ocean observation plan formulation support program.

3. In the ocean observation plan formulation support program according to Claim 1, causing the computer to execute a process of allowing a user to specify the number p of the observation points, an ocean observation plan formulation support program.

4. An apparatus for assisting in the planning of ocean observations, comprising: an arithmetic processing unit that performs arithmetic processing; an input / output unit that receives a user's operation input and presents the result of the arithmetic processing performed by the arithmetic processing unit to the user; The arithmetic processing unit divides a sea area where the ocean observations are to be performed into a plurality of regions based on the correlation relationship between each point in the sea area; The input / output unit displays the plurality of regions on a screen and presents them to the user; The arithmetic processing unit: a process of obtaining p observation points where the analysis error covariance is minimized within the sea area; a process of obtaining the degree of correlation between each of the p observation points and each point in the sea area; Performing a process of grouping each point in the sea area based on the correlation degree for each of the observation points; Dividing the sea area into the plurality of regions based on the result of the grouping; In the process of obtaining the p observation points, Performing a process of creating an observation matrix for each of a plurality of grid points set at predetermined intervals in the sea space within the sea area; Performing a process of obtaining the analysis error covariance for each of the plurality of grid points based on the observation matrix; An apparatus for supporting the formulation of an ocean observation plan, which repeatedly performs a process of selecting, as the observation points, grid points among the plurality of grid points where the analysis error covariance is the minimum, a plurality of times to obtain the p observation points.

5. In the apparatus for supporting the formulation of an ocean observation plan according to claim 4, In the process of obtaining the p observation points, the arithmetic processing unit includes, among the plurality of grid points, grid points corresponding to positions specified on the map showing the sea area by the input operation of the user received by the input / output unit in the p observation points. An apparatus for supporting the formulation of an ocean observation plan.

6. In the apparatus for supporting the formulation of an ocean observation plan according to claim 4, The arithmetic processing unit sets, as the number p of the observation points, the number specified by the input operation of the user received by the input / output unit. An apparatus for supporting the formulation of an ocean observation plan.

7. A method for supporting the formulation of an ocean observation plan using a computer, By the computer, Performing a process of obtaining p observation points where the analysis error covariance is the minimum in the sea area where the ocean observation is to be performed; Performing a process of obtaining the correlation degree between each of the p observation points and each point in the sea area; Performing a process of grouping each point in the sea area for each of the observation points based on the correlation degree; Performing a process of dividing the sea area into a plurality of regions based on the result of the grouping; Performing a process of displaying the plurality of regions on the screen; In the process of obtaining the p observation points, Performing a process of creating an observation matrix for each of a plurality of grid points set at predetermined intervals in the sea space within the sea area; Performing a process of obtaining the analysis error covariance for each of the plurality of grid points based on the observation matrix; A method for supporting the formulation of an ocean observation plan, which repeatedly performs a process of selecting, as the observation points, grid points among the plurality of grid points where the analysis error covariance is the minimum, a plurality of times to obtain the p observation points.

8. In the method for supporting the formulation of an ocean observation plan according to claim 7, Causing the user to specify any position on the map showing the sea area; In the process of obtaining the p observation points, a method for supporting the formulation of an ocean observation plan, which includes, among the plurality of grid points, the grid points corresponding to the positions specified by the user among the p observation points.

9. In the method for supporting the formulation of an ocean observation plan according to Claim 7, A method for supporting the formulation of an ocean observation plan, which allows the user to specify the number p of the observation points.

Images