阅读说明：本技术 航空重力测量点的布格重力异常值的计算装置 (Calculation device for bump gravity abnormal value of aviation gravity measurement point ) 是由 罗锋 郭志宏 骆遥 屈进红 孙艳云 王明 李行素 于 2020-06-08 设计创作，主要内容包括：本发明公开了一种航空重力测量点的布格重力异常值的计算装置,包括：获取模块,用于获取航空重力测量点的地理位置信息、所述航空重力测量点的GNSS高程以及重力异常值；投影模块,用于将航空重力测量点的GNSS高程转换至大地水准面高程,并将所述航空重力测量点的地理位置信息投影至地形高程数据网格中,得到高程投影点；构建模块；地形改正值计算模块；布格重力异常计算模块,用于根据该航空重力测量点的重力异常值以及所述地形改正量获得该航空重力测量点的布格重力异常。本发明提供的航空重力测量点的布格重力异常值的计算装置,可以获得航空布格重力异常。(The invention discloses a device for calculating a grid gravity abnormal value of an aviation gravity measurement point, which comprises: the acquisition module is used for acquiring the geographical position information of the aviation gravity measurement point, the GNSS elevation of the aviation gravity measurement point and the gravity abnormal value; the projection module is used for converting the GNSS elevation of the aerial gravity measurement point to the ground level elevation, and projecting the geographic position information of the aerial gravity measurement point to a terrain elevation data grid to obtain an elevation projection point; building a module; a terrain correction value calculation module; and the grid gravity anomaly calculation module is used for obtaining the grid gravity anomaly of the aviation gravity measurement point according to the gravity anomaly value of the aviation gravity measurement point and the terrain correction quantity. The calculation device for the grid gravity abnormal value of the aviation gravity measurement point can obtain the aviation grid gravity abnormality.)

1. An apparatus for calculating a Bouger gravity abnormal value of an aviation gravity measurement point, comprising:

the acquisition module is used for acquiring the geographical position information of the aviation gravity measurement point, the GNSS elevation of the aviation gravity measurement point and the gravity abnormal value;

the projection module is used for converting the GNSS elevation of the aerial gravity measurement point to the ground level elevation, and projecting the geographic position information of the aerial gravity measurement point to a terrain elevation data grid to obtain an elevation projection point;

the construction module is used for constructing a first cuboid in a terrain elevation data grid according to four nodes closest to the elevation projection point as a bottom surface, the height of the first cuboid is an average value of elevation values of the four nodes, the terrain elevation data grid comprises grid data formed by coordinate values of a plurality of nodes, and the terrain elevation data grid comprises the elevation values of the plurality of nodes;

the building module is further used for respectively taking the four nodes as starting points, taking the grid step length in the terrain elevation data grid as the side length, extending in the direction away from the projection point, building a second cuboid, wherein the height of the second cuboid is the average value of the elevation values corresponding to the four nodes on the bottom surface, and repeating the step of building the second cuboid until the coordinate value of the node exceeds a preset range;

the terrain correction value calculation module is used for calculating the terrain correction value of each cuboid according to the height difference from the aviation gravity measurement point to each cuboid, the distance difference from the aviation gravity measurement point to each cuboid, the density of each cuboid and the universal gravitation constant, wherein each cuboid is a first cuboid or a second cuboid;

the terrain correction value calculation module is used for summing all the terrain correction values to obtain the terrain correction value of the aviation gravity measurement point;

and the grid gravity anomaly calculation module is used for obtaining the grid gravity anomaly of the aviation gravity measurement point according to the gravity anomaly value of the aviation gravity measurement point and the terrain correction quantity.

Technical Field

The invention relates to the technical field of aviation gravity measurement, in particular to a device for calculating a grid gravity abnormal value of an aviation gravity measurement point.

Background

With the rapid development of science and technology, the accuracy requirement of gravity measurement in the field of gravity exploration is higher and higher.

And the data obtained by the gravity measurement is gravity observation data. The gravity observation data includes abnormalities caused by geological bodies with uneven underground density, and also includes influences caused by different terrains around each measuring point, different latitudes where the measuring points are located and the like, so the gravity observation data needs to be corrected. The aviation gravity measurement is one kind of gravity measurement, and is especially one kind of gravity measurement method with aviation gravity measurement system installed onto airplane for continuous measurement.

Currently, the treatment of gravity data by a Bouguer gravity anomaly (Bouguer gravity anomaly) has a very important meaning. The ground weave gravity anomaly is typically obtained by ground gravity measurements. The gravity obtained by subtracting the normal gravity value after the latitude correction, the height correction, the middle layer correction and the terrain correction of the observation result of the ground gravimeter is called as the ground grid gravity anomaly.

Based on this, the inventor of the present application finds that the calculation method for acquiring the bogey anomaly from the ground gravity measurement value in the prior art has a certain unsuitability in the calculation for acquiring the aviation bogey anomaly from the aviation gravity anomaly.

The information disclosed in this background section is only for enhancement of understanding of the general background of the invention and should not be taken as an acknowledgement or any form of suggestion that this information forms the prior art already known to a person skilled in the art.

Disclosure of Invention

In order to solve the above problem, an object of an embodiment of the present invention is to provide a device for calculating a grid gravity abnormal value of an aviation gravity measurement point.

In a first aspect, an embodiment of the present invention provides a method for calculating a grid gravity abnormal value of an aviation gravity measurement point, including: acquiring geographic position information of an aviation gravity measuring point, GNSS elevation of the aviation gravity measuring point and a gravity abnormal value; converting the GNSS elevation of the aerial gravity measurement point to the ground level elevation, and projecting the geographical position information of the aerial gravity measurement point to a terrain elevation data grid to obtain an elevation projection point; in a terrain elevation data grid, constructing a first cuboid according to four nodes closest to the elevation projection point as a bottom surface, wherein the height of the first cuboid is an average value of elevation values of the four nodes, the terrain elevation data grid comprises grid data formed by coordinate values of a plurality of nodes, and the terrain elevation data grid comprises the elevation values of the plurality of nodes; respectively taking the four nodes as starting points, taking the grid step length in the terrain elevation data grid as the side length, extending in the direction away from the projection point, constructing a second cuboid, wherein the height of the second cuboid is the average value of the elevation values corresponding to the four nodes on the bottom surface, and repeating the step of constructing the second cuboid until the coordinate value of the node exceeds a preset range; calculating a terrain correction value of each cuboid according to the height difference from the aviation gravity measurement point to each cuboid, the distance difference from the aviation gravity measurement point to each cuboid, the density of each cuboid and a universal gravitation constant, wherein each cuboid is a first cuboid or a second cuboid; summing all terrain correction values to obtain the terrain correction quantity of the aviation gravity measurement point; and acquiring the Booth gravity anomaly of the aviation gravity measurement point according to the gravity anomaly value of the aviation gravity measurement point and the terrain correction quantity.

In one possible implementation, the terrain elevation data grid comprises a local terrain elevation data grid and an area terrain elevation data grid; the preset range of the node coordinate values in the local terrain elevation data grid is 0-20km, and the preset range of the node coordinate values in the regional terrain elevation data grid is 20km-166.7 km.

In a possible implementation manner, the coordinate value of the node exceeding the preset range is: the coordinate value of the node exceeds 166.7 km.

In one possible implementation manner, before constructing, in the terrain elevation data grid, a first rectangular parallelepiped according to four nodes closest to the elevation projection point as a bottom surface, the method further includes: projecting the geographical position information of the measuring points to a rock density distribution data grid to obtain density projection points, wherein the rock density distribution data grid comprises grid data formed by coordinate values of a plurality of nodes, and the rock density distribution data grid comprises density values of the plurality of nodes; the density of the first cuboid is a density value corresponding to the density projection point.

In one possible implementation manner, the density value of the second cuboid is an average value of four nodes projected to the rock density distribution data grid respectively on the bottom surface of the second cuboid.

In one possible implementation mode, the division rule of the local terrain elevation data grids is high precision and a large scale; the division rule of the regional terrain elevation data grids is low-level precision and a small scale.

In one possible implementation, the calculation method further includes: acquiring the geographical position information of the next aviation gravity measurement point, the GNSS elevation and the gravity measurement value of the next aviation gravity measurement point, and repeating the step of calculating the terrain correction amount corresponding to the next aviation gravity measurement point until the calculation of the terrain correction amount of all the aviation gravity measurement points is completed; performing point-by-point low-pass filtering on the terrain correction quantities of all the aviation gravity measurement points, wherein the filtering scale of the low-pass filtering corresponds to the filtering scale of the aviation gravity measurement point gravity anomaly; the obtaining of the grid gravity abnormal value of the aviation gravity measurement point according to the gravity abnormal value of the aviation gravity measurement point and the terrain correction amount comprises: and obtaining the grid gravity abnormal value of the aviation gravity measurement point according to the gravity abnormal value of the aviation gravity measurement point and the filtered terrain correction quantity.

In a second aspect, an embodiment of the present invention further provides a device for calculating a grid gravity abnormal value of an aviation gravity measurement point, including: the acquisition module is used for acquiring the geographical position information of the aviation gravity measurement point, the GNSS elevation of the aviation gravity measurement point and the gravity abnormal value; the projection module is used for converting the GNSS elevation of the aerial gravity measurement point to the ground level elevation, and projecting the geographic position information of the aerial gravity measurement point to a terrain elevation data grid to obtain an elevation projection point; the construction module is used for constructing a first cuboid in a terrain elevation data grid according to four nodes closest to the elevation projection point as a bottom surface, the height of the first cuboid is an average value of elevation values of the four nodes, the terrain elevation data grid comprises grid data formed by coordinate values of a plurality of nodes, and the terrain elevation data grid comprises the elevation values of the plurality of nodes; the building module is further used for respectively taking the four nodes as starting points, taking the grid step length in the terrain elevation data grid as the side length, extending in the direction away from the projection point, building a second cuboid, wherein the height of the second cuboid is the average value of the elevation values corresponding to the four nodes on the bottom surface, and repeating the step of building the second cuboid until the coordinate value of the node exceeds a preset range; the terrain correction value calculation module is used for calculating the terrain correction value of each cuboid according to the height difference from the aviation gravity measurement point to each cuboid, the distance difference from the aviation gravity measurement point to each cuboid, the density of each cuboid and the universal gravitation constant, wherein each cuboid is a first cuboid or a second cuboid; the terrain correction value calculation module is used for summing all the terrain correction values to obtain the terrain correction value of the aviation gravity measurement point; and the grid gravity anomaly calculation module is used for obtaining the grid gravity anomaly of the aviation gravity measurement point according to the gravity anomaly value of the aviation gravity measurement point and the terrain correction quantity.

In a third aspect, an embodiment of the present invention further provides a storage medium, where the storage medium stores computer-executable instructions, and the computer-executable instructions are configured to perform the calculation of the grid gravity abnormal value of the aviation gravity measurement point.

According to the device for calculating the grid distribution gravity abnormal value of the aviation gravity measurement point, the GNSS elevation of the aviation gravity measurement point is converted into the ground level surface elevation, the geographic position information of the aviation gravity measurement point is projected into the terrain elevation data grid, a first cuboid is constructed by taking four nodes closest to the elevation projection point as the bottom surface, and the cuboid is repeatedly constructed by taking the four nodes as the starting points; calculating the terrain correction value of each cuboid, and summing all the terrain correction values to obtain the terrain correction amount of the aviation gravity measurement point; and acquiring the grid distribution gravity anomaly of the aviation gravity measurement point according to the gravity anomaly value of the aviation gravity measurement point and the terrain correction amount, so that the calculation of the aviation grid distribution gravity anomaly can be realized.

Drawings

In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings used in the description of the embodiments or the prior art will be briefly described below, it is obvious that the drawings in the following description are only some embodiments of the present invention, and for those skilled in the art, other drawings can be obtained according to the drawings without creative efforts.

FIG. 1 is a flow chart illustrating a calculation method of a Booth gravity abnormal value of an aviation gravity measurement point according to an embodiment of the present invention;

FIG. 2 is a schematic diagram of a terrain elevation data grid provided by an embodiment of the present invention;

FIG. 3 is a schematic diagram illustrating a cuboid calculation method according to an embodiment of the present invention;

fig. 4 is a schematic structural diagram of a calculation apparatus for a grid gravity abnormal value of an aviation gravity measurement point according to an embodiment of the present invention.

Detailed Description

The following detailed description of the present invention is provided in conjunction with the accompanying drawings, but it should be understood that the scope of the present invention is not limited to the specific embodiments.

In the description of the present invention, it is to be understood that the terms "center", "longitudinal", "lateral", "length", "width", "thickness", "upper", "lower", "front", "rear", "left", "right", "vertical", "horizontal", "top", "bottom", "inner", "outer", "clockwise", "counterclockwise", and the like, indicate orientations and positional relationships based on those shown in the drawings, and are used only for convenience of description and simplicity of description, and do not indicate or imply that the device or element being referred to must have a particular orientation, be constructed and operated in a particular orientation, and thus, should not be considered as limiting the present invention.

Throughout the specification and claims, unless explicitly stated otherwise, the word "comprise", or variations such as "comprises" or "comprising", will be understood to imply the inclusion of a stated element or component but not the exclusion of any other element or component.

Furthermore, the terms "first", "second" and "first" are used for descriptive purposes only and are not to be construed as indicating or implying relative importance or implicitly indicating the number of technical features indicated. Thus, a feature defined as "first" or "second" may explicitly or implicitly include one or more of that feature. In the description of the present invention, "a plurality" means two or more unless specifically defined otherwise.

Because the actual topographic relief change of the aviation gravity measurement is complex, the terrain around the aviation gravity measurement point is divided into a plurality of cuboids, the gravity influence value of the terrain quality of each cuboid on the aviation gravity measurement point is calculated, and finally the terrain influence value of the point is obtained through accumulation and summation. The terrain around the measuring point is an extremely irregular curved surface and cannot be expressed by a function formula of a space coordinate, the digital simulation of the ground terrain is realized by adopting terrain elevation data, and grid data of regular terrain elevation is established. When the aviation gravity terrain correction calculation is carried out, reading elevation data in the local terrain elevation grid data and the regional terrain elevation grid data in a certain range according to grid nodes.

The embodiment of the invention provides a flow chart of a calculation method of a grid gravity abnormal value of an aviation gravity measurement point, which is shown in a figure 1 and comprises the following steps: step 1-step 7.

Step 1, acquiring geographic position information of an aviation gravity measurement point, GNSS elevation of the aviation gravity measurement point and a gravity abnormal value;

and 2, converting the GNSS elevation of the aerial gravity measurement point to the ground level surface elevation, and projecting the geographical position information of the aerial gravity measurement point to a terrain elevation data grid to obtain an elevation projection point.

In one implementation, the terrain elevation data grid comprises a local terrain elevation data grid and an area terrain elevation data grid; the preset range of the node coordinate values in the local terrain elevation data grid is 0-20km, and the preset range of the node coordinate values in the regional terrain elevation data grid is 20km-166.7 km. Referring to fig. 2, which is a schematic diagram of a terrain elevation data grid, EFGH is a local terrain area, IJKL is an area terrain, and a middle black area is a range of an airborne gravity measurement point.

The division rule of the local terrain elevation data grids is high precision and large scale; the division rule of the regional terrain elevation data grids is low-level precision and a small scale.

Step 3, constructing a first cuboid in a terrain elevation data grid according to four nodes closest to the elevation projection point as a bottom surface, wherein the height of the first cuboid is an average value of elevation values of the four nodes, the terrain elevation data grid comprises grid data formed by coordinate values of a plurality of nodes, and the terrain elevation data grid comprises the elevation values of the plurality of nodes;

step 4, respectively taking the four nodes as starting points, taking the grid length in the terrain elevation data grid as side length, extending the grid length in the direction away from the projection point, constructing a second cuboid, wherein the height of the second cuboid is the average value of elevation values corresponding to the four nodes on the bottom surface, and repeating the step of constructing the second cuboid until the coordinate value of the node exceeds a preset range;

the coordinate value of the node exceeds the preset range, and the coordinate value of the node can exceed 166.7 km.

Step 5, calculating a terrain correction value of each cuboid according to the height difference from the aviation gravity measurement point to each cuboid, the distance difference from the aviation gravity measurement point to each cuboid, the density of each cuboid and a universal gravitation constant, wherein each cuboid is a first cuboid or a second cuboid;

step 6, summing all terrain correction values to obtain the terrain correction amount of the aviation gravity measurement point;

and 7, acquiring the lattice gravity anomaly of the aviation gravity measurement point according to the gravity anomaly value of the aviation gravity measurement point and the terrain correction amount.

In the method for calculating the grid distribution gravity abnormal value of the airborne gravity measurement point, the GNSS elevation of the airborne gravity measurement point is converted into the elevation of the ground level surface, the geographical position information of the airborne gravity measurement point is projected into the terrain elevation data grid, a first cuboid is constructed by taking four nodes closest to the elevation projection point as the bottom surface, and the cuboid is repeatedly constructed by taking the four nodes as the starting points; calculating the terrain correction value of each cuboid, and summing all the terrain correction values to obtain the terrain correction amount of the aviation gravity measurement point; and acquiring the grid distribution gravity anomaly of the aviation gravity measurement point according to the gravity anomaly value of the aviation gravity measurement point and the terrain correction amount, so that the calculation of the aviation grid distribution gravity anomaly can be realized.

In one implementation, step 3 further includes, before: projecting the geographical position information of the measuring points to a rock density distribution data grid to obtain density projection points, wherein the rock density distribution data grid comprises grid data formed by coordinate values of a plurality of nodes, and the rock density distribution data grid comprises density values of the plurality of nodes; the density of the first cuboid is a density value corresponding to the density projection point.

Correspondingly, the density value of the second cuboid is an average value of four nodes on the bottom surface of the second cuboid which are respectively projected into the rock density distribution data grid.

In one implementation, step 5 may specifically include:

for each cuboid, calculating an aviation gravity terrain correction value deltag of the cuboid according to a first formula, wherein the first formula comprises:

wherein G is a gravitational constant, ρ is a density of each rectangular solid, and X, Y, and Z are observation point spatial positions, (X1, X2), (Y1, Y2), (Z1, and Z2) are rectangular solid spatial positions. And R is the distance from the observation point to each vertex of the cuboid.

Further, the implementation of step 6 and step 7 is as follows. Referring to fig. 3, which is a schematic diagram of the cuboid calculation method provided in this embodiment, the position of the projection of the aerial gravity observation point to the terrain elevation grid data is a point P (X0, Y0), where a point A, B, C, D is a node closest to the point P in the terrain grid data. Calculating the terrain influence value delta g from the cuboid to the aeronautical gravity measurement point by taking the point A, B, C, D as the side of the cuboid and the height as the average value of the elevation values of the nodes A, B, C, D₁. Expanding towards two directions of the grid, and calculating the terrain influence value delta g from the fitting cuboid corresponding to every four nodes to the aviation gravity measurement point A by analogy_i。

The sum of the terrain correction values of the aviation gravity measurement points is as follows:

the Booth gravity abnormal value of the aviation gravity measurement point is as follows:

G_{cloth grid}＝G_{From sky}-G_Ground (2)

Wherein G is_{From sky}The space gravity abnormal value of the aviation gravity measurement point A is obtained through aviation gravity measurement.

In one implementation, step 7 may be followed by: and acquiring the geographical position information of the next aviation gravity measurement point, the GNSS elevation and the gravity measurement value of the next aviation gravity measurement point, and repeating the steps of calculating the terrain correction amount corresponding to the next aviation gravity measurement point until the calculation of the terrain correction amount of all the aviation gravity measurement points is completed.

Specifically, the calculation steps from step 1 to step 6 are repeated for each aviation gravity measurement point until the terrain correction amount of all aviation gravity measurement points is calculated.

In one implementation, the method may further include: and performing point-by-point low-pass filtering on the terrain correction quantities of all the aviation gravity measurement points, wherein the filtering scale of the low-pass filtering corresponds to the filtering scale of aviation space gravity anomaly.

Accordingly, step 7 may be: and obtaining the grid gravity abnormal value of the aviation gravity measurement point according to the gravity abnormal value of the aviation gravity measurement point and the filtered terrain correction quantity.

The embodiment of the present invention provides a schematic structural diagram of a device for calculating a grid gravity abnormal value of an aviation gravity measurement point, which is shown in fig. 4, and includes: the device comprises an acquisition module 1, a projection module 2, a construction module 3, a terrain correction value calculation module 4, a terrain correction amount calculation module 5 and a grid gravity anomaly calculation module 6.

The acquisition module 1 is used for acquiring the geographical position information of an aviation gravity measurement point, the GNSS elevation of the aviation gravity measurement point and a gravity abnormal value;

the projection module 2 is used for converting the GNSS elevation of the aviation gravity measurement point to the ground level elevation, and projecting the geographic position information of the aviation gravity measurement point to a terrain elevation data grid to obtain an elevation projection point;

the construction module 3 is used for constructing a first cuboid in a terrain elevation data grid according to four nodes closest to the elevation projection point as a bottom surface, wherein the height of the first cuboid is an average value of elevation values of the four nodes, the terrain elevation data grid comprises grid data formed by coordinate values of a plurality of nodes, and the terrain elevation data grid comprises the elevation values of the plurality of nodes;

the building module 3 is further configured to respectively use the four nodes as starting points, use the grid length in the terrain elevation data grid as a side length, extend in a direction away from the projection point, build a second cuboid, use the height of the second cuboid as an average value of elevation values corresponding to the four nodes on the bottom surface, and repeat the step of building the second cuboid until the coordinate value of the node exceeds a preset range;

the terrain correction value calculating module 4 is used for calculating the terrain correction value of each cuboid according to the height difference from the aviation gravity measuring point to each cuboid, the distance difference from the aviation gravity measuring point to each cuboid, the density of each cuboid and the universal gravitation constant, wherein each cuboid is a first cuboid or a second cuboid;

the terrain correction amount calculation module 5 is used for summing all the terrain correction values to obtain the terrain correction amount of the aviation gravity measurement point;

and the lattice gravity anomaly calculation module 6 is used for obtaining the lattice gravity anomaly of the aviation gravity measurement point according to the gravity anomaly value of the aviation gravity measurement point and the terrain correction quantity.

The device for calculating the grid distribution gravity abnormal value of the airborne gravity measurement point provided by this embodiment is configured to convert the GNSS elevation of the airborne gravity measurement point to the elevation of the ground level surface, project the geographical location information of the airborne gravity measurement point into the terrain elevation data grid, construct a first cuboid according to four nodes closest to the elevation projection point as a bottom surface, and repeatedly construct a cuboid with the four nodes as starting points; calculating the terrain correction value of each cuboid, and summing all the terrain correction values to obtain the terrain correction amount of the aviation gravity measurement point; and acquiring the grid distribution gravity anomaly of the aviation gravity measurement point according to the gravity anomaly value of the aviation gravity measurement point and the terrain correction amount, so that the calculation of the aviation grid distribution gravity anomaly can be realized.

An embodiment of the present invention further provides a storage medium storing computer-executable instructions including a program for executing the method for calculating the grid gravity abnormal value of the aviation gravity measurement point, where the computer-executable instructions may execute the method in any of the method embodiments.

The storage medium may be any available medium or data storage device that can be accessed by a computer, including but not limited to magnetic memory (e.g., floppy disks, hard disks, magnetic tape, magneto-optical disks (MOs), etc.), optical memory (e.g., CDs, DVDs, BDs, HVDs, etc.), and semiconductor memory (e.g., ROMs, EPROMs, EEPROMs, nonvolatile memory (NAND FLASH), Solid State Disks (SSDs)), etc.

As will be appreciated by one skilled in the art, embodiments of the present application may be provided as a method, system, or computer program product. Accordingly, the present application may take the form of an entirely hardware embodiment, an entirely software embodiment or an embodiment combining software and hardware aspects.

The present application is described with reference to flowchart illustrations and/or block diagrams of methods, apparatus (systems), and computer program products according to embodiments of the application. It will be understood that each flow and/or block of the flow diagrams and/or block diagrams, and combinations of flows and/or blocks in the flow diagrams and/or block diagrams, can be implemented by computer program instructions. These computer program instructions may be provided to a processor of a general purpose computer, special purpose computer, embedded processor, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, create means for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.

These computer program instructions may also be stored in a computer-readable memory that can direct a computer or other programmable data processing apparatus to function in a particular manner, such that the instructions stored in the computer-readable memory produce an article of manufacture including instruction means which implement the function specified in the flowchart flow or flows and/or block diagram block or blocks.

These computer program instructions may also be loaded onto a computer or other programmable data processing apparatus to cause a series of operational steps to be performed on the computer or other programmable apparatus to produce a computer implemented process such that the instructions which execute on the computer or other programmable apparatus provide steps for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.

The foregoing descriptions of specific exemplary embodiments of the present invention have been presented for purposes of illustration and description. It is not intended to limit the invention to the precise form disclosed, and obviously many modifications and variations are possible in light of the above teaching. The exemplary embodiments were chosen and described in order to explain certain principles of the invention and its practical application to enable one skilled in the art to make and use various exemplary embodiments of the invention and various alternatives and modifications as are suited to the particular use contemplated. It is intended that the scope of the invention be defined by the claims and their equivalents.