CN103308942A - Method and system for visualizing seismic data - Google Patents

Abstract

The invention discloses a method and system for visualizing seismic data, and belongs to the computer field. The method is characterized by: transforming the seismic data to obtain volume data through a CPU and blocking the volume data to obtain subblocks; filtering the subblocks to obtain visible subblocks; and transmitting the visible subblocks to a GPU to carry out an image rendering. Therefore, computation burden of the GPU is greatly reduced; the seismic data can be 3D visualized through a common GPU in a common computer; and the method and the system are convenient to use, and cost is saved.

Description

A kind of method and system of visual earthquake data

Technical field

The present invention relates to computer realm, particularly a kind of method and system of visual earthquake data.

Background technology

At geophysics field, seismic prospecting is one of important method in the geophysical survey.It is with artificial method earthquake-wave-exciting, and with seismic prospecting instrument the vibrations of the earth is recorded on the tape, with computing machine the data that the field obtains is handled then, obtains information and the oil gas information of relevant underground structure and lithology.Wherein, the geological data that obtains is carried out three-dimensional picture show it is an indispensable part during the three-dimensional data three-dimensional visualization is explained, can not only express and show the feature of all kinds of structural planes and their space distribution rule clearly like this, see tectonic structure, mutual relationship and distribution intuitively, thereby be conducive to pinpoint the problems, problem analysis, for next step geologic interpretation lays the foundation.

Because the geological data amount that obtains of exploration is huge, the calculated amount that produces when carrying out three-dimensional visualization is huge, and 3-D seismics in the past is visual all can only be by disposing powerful supercomputer or professional graphics workstation goes to finish.

In realizing process of the present invention, the inventor finds that there is following problem at least in prior art:

Supercomputer or professional graphics workstation cost all compare expensive, use also inconvenient.

Summary of the invention

In order to solve problems of the prior art, the embodiment of the invention provides a kind of method and system of visual earthquake data.Described technical scheme is as follows:

A kind of method of visual earthquake data, described method comprises:

By CPU (Central Processing Unit, central processing unit) 3D seismic data is converted to volume data, and described volume data is decomposed into sub-piece;

Described CPU carries out filtering with described sub-piece and obtains visible sub-piece;

GPU (Graphic Processing Unit, graphic process unit) carries out obtaining after secondary decomposes the multi-resolution models of described 3D seismic data with described visible sub-piece, wherein, comprises discrete approximation signal and detail signal in the described multi-resolution models;

Transmit the described discrete approximation signal under the J class resolution ratio, the order that is begun to successively decrease according to resolution levels by the J class resolution ratio transmits described detail signal step by step again, and wherein, J is integer;

According to the direction of visual lines vector, the order that arrives according to described discrete approximation signal and detail signal carries out 3 D rendering.

Further, described 3D seismic data is converted to volume data, comprises:

To be converted to a plurality of voxels at the 3D seismic data that a plurality of sampled points obtain, and described a plurality of voxel combinations will be obtained described volume data;

Wherein, comprise 3D seismic data, rgb color value and the transparence information that described each voxel corresponding sampling points obtains in each voxel at least.

Further, described volume data is decomposed into sub-piece, comprises:

Adopt the structure of Octree that described volume data is carried out 2 * 2 * 2 decomposition, obtain 8 sub-pieces, and when satisfying first when pre-conditioned, the sub-piece in lasting described 8 sub-pieces that decomposition is obtained decomposes and obtains littler sub-piece;

Wherein, described first pre-conditioned the comprising:

The consistance of the voxel in the described sub-piece less than the size of user-defined critical value, described sub-piece greater than the size of the size of the video memory corresponding with described GPU, the described sub-piece size greater than described user-defined boy's piece.

Further, describedly described sub-piece carried out filtering obtain visible sub-piece, comprising:

Create the voxel count table;

Obtain the ratio of the visible voxel in all sub-pieces according to the transparent transmission function calculation;

Ratio judgement by the visible voxel in described all sub-pieces obtains described visible sub-piece;

Wherein, described voxel count table is used for statistics through the quantity of the visible voxel behind the described transparent transmission function category, and described transparent transmission function obtains according to described direction of visual lines vector.

Further, described GPU carries out described visible sub-piece to obtain after secondary decomposes before the multiresolution module of described 3D seismic data, and described method also comprises:

Described visible sub-piece is transferred in the video memory of described GPU correspondence.

Further, described GPU carries out obtaining after secondary decomposes the multi-resolution models of described 3D seismic data with described visible sub-piece, comprising:

Each visible sub-piece is carried out 3 D wavelet decompose the level that the order of successively decreasing with resolution levels is created described visible sub-piece correspondence;

Revise the data value at two adjacent described visible sub-block boundary places, with the continuity of the level transition that keeps each resolution levels correspondence;

Carry out the opacity correction by with the resolution of different stage same described visible sub-piece being carried out resampling antithetical phrase piece;

Level to the multiresolution correspondence of described visible sub-piece carries out adaptively selected;

Calculate the multi-resolution models of described 3D seismic data;

Wherein, comprise discrete approximation signal and detail signal in the described multi-resolution models.

Further, the data value at two described visible sub-block boundary places that described modification is adjacent to keep the continuity of level transition, comprising:

At two adjacent described visible sub-block boundary places, copy the data value of current level to last one high-level level, or after the data value of current level carried out interpolation, replace the data value of last one high-level level boundary.

Further, described multi-resolution models is specially:

$f(x,y,z)=\underset{i,j,k\∈z}{\mathrm{\Σ}}{\left(A_m^{d}f\right)}_{i,j,k}{\mathrm{\φ}}_{m,i,j,k}(x,y,z)+\underset{m=1}{\overset{J}{\mathrm{\Σ}}}\underset{n=1}{\overset{7}{\mathrm{\Σ}}}\underset{i,j,k\∈z}{\mathrm{\Σ}}{\left(D_m^{n}f\right)}_{i,j,k}{\mathrm{\ψ}}_{m,i,j,k}^n(x,y,z)$

Wherein, f (x, y z) are described volume data, and integer J is the progression that the wavelet multiresolution rate is analyzed,

Be the discrete level and smooth approximation signal of described volume data under the J class resolution ratio,

Be the discrete detail signal of described volume data under the J class resolution ratio.

Further, described according to the direction of visual lines vector, the order that arrives according to described discrete approximation signal and detail signal carries out 3 D rendering, comprising:

Be orthogonal on the picture plane of described direction of visual lines vector, drawing out the 3-D view of default resolution with the discrete approximation signal under the described J resolution;

Begin the order details of use signal that successively decreases according to the resolution described 3-D view of refinement successively by J resolution again.

Further, after described direction of visual lines vector changed, described method also comprised:

According to the direction of visual lines vector after the described change, the order that arrives according to described discrete approximation signal and detail signal carries out 3 D rendering.

A kind of system of visual earthquake data, described system comprises: central processing unit and image processor;

Wherein, described central processing unit comprises:

Modular converter is used for 3D seismic data is converted to volume data;

Decomposing module is used for described volume data is decomposed into sub-piece;

Filtering module is used for that described sub-piece is carried out filtering and obtains visible sub-piece;

Described graphic process unit comprises:

Model computation module, the multi-resolution models for described visible sub-piece being carried out obtain after secondary decomposes described 3D seismic data wherein, comprises discrete approximation signal and detail signal in the described multi-resolution models;

Transport module is used for the described discrete approximation signal under the transmission J class resolution ratio, and the order that is begun to successively decrease according to resolution levels by the J class resolution ratio transmits described detail signal step by step again, and wherein, J is integer;

Drafting module is used for according to the direction of visual lines vector, and the order that arrives according to described discrete approximation signal and detail signal carries out 3 D rendering.

Further, described modular converter, concrete being used for will be converted to a plurality of voxels at the 3D seismic data that a plurality of sampled points obtain, and described a plurality of voxel combinations are obtained described volume data;

Wherein, comprise 3D seismic data, rgb color value and the transparence information that described each voxel corresponding sampling points obtains in each voxel at least.

Further, described decomposing module, the concrete structure that is used for the employing Octree is carried out 2 * 2 * 2 decomposition to described volume data, obtains 8 sub-pieces, and when satisfying first when pre-conditioned, the sub-piece in described 8 the sub-pieces that continue decomposition is obtained decomposes and obtains littler sub-piece;

Wherein, described first pre-conditioned the comprising:

The consistance of the voxel in the described sub-piece less than the size of user-defined critical value, described sub-piece greater than with the size of the size of video memory, the described sub-piece size greater than described user-defined boy's piece.

Further, described filtering module comprises:

Creating unit is used for creating the voxel count table;

The ratio computing unit is for the ratio that obtains the visible voxel of all sub-pieces according to the transparent transmission function calculation;

Judging unit, the ratio that is used for the visible voxel by described all sub-pieces are judged and are obtained described visible sub-piece;

Wherein, described voxel count table is used for statistics through the quantity of the visible voxel behind the described transparent transmission function category, and described transparent transmission function obtains according to described direction of visual lines vector.

Further, described central processing unit also comprises:

Memory module is used for before described model computation module carries out described visible sub-piece to obtain the multiresolution module of described 3D seismic data after secondary decomposes described visible sub-piece being transferred in the described video memory.

Further, described model computation module comprises:

The 3 D wavelet resolving cell is used for that each visible sub-piece is carried out 3 D wavelet and decomposes the level that the order of successively decreasing with resolution levels is created described visible sub-piece correspondence;

Revise the unit, be used for revising the data value at two adjacent described visible sub-block boundary places, with the continuity of the level transition that keeps each resolution levels correspondence;

Correcting unit is used for carrying out the opacity correction by with the resolution of different stage same described visible sub-piece being carried out resampling antithetical phrase piece;

Selected cell is used for carrying out adaptively selected to the level of the multiresolution correspondence of described visible sub-piece;

The model computing unit is for the multi-resolution models that calculates described 3D seismic data;

Wherein, comprise discrete approximation signal and detail signal in the described multi-resolution models.

Further, described modification unit, copies the data value of current level to last one high-level level at concrete being used at two adjacent described visible sub-block boundary places, or after the data value of current level carried out interpolation, replace the data value of last one high-level level boundary.

Further, described multi-resolution models is specially:

$f(x,y,z)=\underset{i,j,k\∈z}{\mathrm{\Σ}}{\left(A_m^{d}f\right)}_{i,j,k}{\mathrm{\φ}}_{m,i,j,k}(x,y,z)+\underset{m=1}{\overset{J}{\mathrm{\Σ}}}\underset{n=1}{\overset{7}{\mathrm{\Σ}}}\underset{i,j,k\∈z}{\mathrm{\Σ}}{\left(D_m^{n}f\right)}_{i,j,k}{\mathrm{\ψ}}_{m,i,j,k}^n(x,y,z)$

Wherein, f (x, y z) are described volume data, and integer J is the progression that the wavelet multiresolution rate is analyzed, Be the discrete level and smooth approximation signal of described volume data under the J class resolution ratio,

Be the discrete detail signal of described volume data under the J class resolution ratio.

Further, described drafting module comprises:

First drawing unit is used for drawing out the 3-D view of default resolution with the discrete approximation signal under the described J resolution on the picture plane that is orthogonal to described direction of visual lines vector;

Second drawing unit is for the order details of use signal that is begun by J resolution to successively decrease according to the resolution described 3-D view of refinement successively.

Further, described image processor also comprises:

Redraw the molding piece, be used for after described direction of visual lines vector changes, according to the direction of visual lines vector after the described change, the order that arrives according to described discrete approximation signal and detail signal carries out 3 D rendering.

The beneficial effect that the technical scheme that the embodiment of the invention provides is brought is: by using among the CPU geological data is converted to volume data, and volume data carried out piecemeal, filtration obtains visible sub-piece, transfer to GPU and carry out image and draw in that visible sub-piece is passed, can reduce the calculated amount of GPU greatly, make the three-dimensional visualization that on common computer, also can carry out geological data by common GPU, easy to use, and save cost.

Description of drawings

In order to be illustrated more clearly in the technical scheme in the embodiment of the invention, the accompanying drawing of required use is done to introduce simply in will describing embodiment below, apparently, accompanying drawing in describing below only is some embodiments of the present invention, for those of ordinary skills, under the prerequisite of not paying creative work, can also obtain other accompanying drawing according to these accompanying drawings.

Fig. 1 is the schematic flow sheet of the method for a kind of visual earthquake data that provide in the embodiment of the invention 1;

Fig. 2 is the schematic flow sheet of the method for a kind of visual earthquake data that provide in the embodiment of the invention 2;

Fig. 3 is the synoptic diagram that the one-dimensional data of 4 levels of usefulness of providing in the embodiment of the invention 2 is expressed three wavelet decomposition;

Fig. 4 is the adaptive synoptic diagram of the adjacent sub-blocks with 4 levels that provides in the embodiment of the invention 2;

Fig. 5 be provide in the embodiment of the invention 2 same sub-piece is carried out the synoptic diagram of resampling with different resolution;

Fig. 6 is the structural representation of the system of a kind of visual earthquake data that provide in the embodiment of the invention 3;

Fig. 7 is the structural representation of the central processing unit that provides in the embodiment of the invention 3;

Fig. 8 is the structural representation of the image processor that provides in the embodiment of the invention 3;

Fig. 9 is the structural representation of the filtering module that provides in the embodiment of the invention 3;

Figure 10 is second kind of structural representation of the central processing unit that provides in the embodiment of the invention 3;

Figure 11 is the structural representation of the model computation module that provides in the embodiment of the invention 3;

Figure 12 is the structural representation of the drafting module that provides in the embodiment of the invention 3;

Figure 13 is second kind of structural representation of the image processor that provides in the embodiment of the invention 3.

Embodiment

For making the purpose, technical solutions and advantages of the present invention clearer, embodiment of the present invention is described further in detail below in conjunction with accompanying drawing.

Embodiment 1

As shown in Figure 1, present embodiment provides a kind of method of visual earthquake data, and this method comprises:

101, CPU is converted to volume data with 3D seismic data, and volume data is decomposed into sub-piece;

102, CPU carries out filtering with sub-piece and obtains visible sub-piece;

103, GPU carries out obtaining after secondary decomposes the multi-resolution models of 3D seismic data with visible sub-piece, wherein, comprises discrete approximation signal and detail signal in the multi-resolution models;

104, the discrete approximation signal under the transmission J class resolution ratio, the order that is begun to successively decrease according to resolution levels by the J class resolution ratio transmits detail signal step by step again;

Wherein, J is integer.

105, according to the direction of visual lines vector, the order that arrives according to discrete approximation signal and detail signal carries out 3 D rendering.

The method of a kind of visual earthquake data that present embodiment provides, by using among the CPU geological data is converted to volume data, and volume data carried out piecemeal, filtration obtains visible sub-piece, transfer to GPU and carry out image and draw in that visible sub-piece is passed, can reduce the calculated amount of GPU greatly, make the three-dimensional visualization that on common computer, also can carry out geological data by common GPU, easy to use, and save cost.

Embodiment 2

Present embodiment 2 provides a kind of method of visual earthquake data, is the improvement of carrying out on the basis of embodiment 1.

Need to prove, in recent years, GPU (Graphics Processing Unit, graphic process unit) technology has obtained develop rapidly, as a kind of dedicated graphics rendering hardware for PC, GPU has significant meaning to the development of three-dimensional real-time rendering technology.GPU is exactly speed with respect to the main advantage of CPU, and the speed advantage of GPU is mainly derived from its unique hardware systems design.The method of the visual earthquake data that present embodiment provides can make to show CPU and GPU parallel processing 3D seismic data thereby reach the three dimensional stress that uses common PC can carry out geological data.

As shown in Figure 2, the method for a kind of visual earthquake data that present embodiment 2 provides specifically comprises the steps:

201, will be converted to a voxel at the 3D seismic data that each sampled point obtains, obtain volume data;

Wherein, for 3D seismic data is carried out volume drawing, need convert the 3D seismic data that each sampled point obtains to a voxel, each voxel has value (for example amplitude), a RGB (Red Green Blue, RGB) color-values and the transparence information of the 3D seismic data that a corresponding sampling points obtains.

Wherein, transparence information is used for the darkness variable of nominal data transparency.

When a zone is explored, can select a plurality of places in this zone to carry out data acquisition, each place is a seismic trace, by seismic prospecting instrument in this place perpendicular to the vibrations of the prone a plurality of sampled points record of level the earth and store on the tape, each seismic trace is converted into a voxel road, and a plurality of voxels road is volume data.

202, above-mentioned volume data is decomposed obtain sub-piece, and it is standby to deposit sub-blocks of data in disk;

Need to prove that LOD (Level Of Details, level of detail) model is the effective ways of describing the magnanimity 3-D data volume, also is simultaneously to accelerate the technology that three-dimensional picture generates to guarantee the three-dimensional picture real-time rendering.The cardinal principle of LOD model is the different minutias of selecting to draw this object for observer's importance according to observed object, carries out real-time rendering.Can under the condition of different levels and different visions (visual angle and scope), adopt the LOD model of different fine degree to represent same object, to improve the display speed of three-dimensional scenic, realize real-time demonstration and the interactive operation of mass data.In geological data three-dimensional visualization process, different with the general object dimensional graphic plotting pursuit sense of reality, because the restriction of screen size and resolution, the 3-D display of geological data does not often need fully to show all sampling points.The amplitude that can reflect underground structure and stratum characteristic from the extracting data of magnanimity forms the data under the different resolution condition, to reduce unnecessary graphics calculations amount, handle by the piecemeal to geological data, the voxel of a seismic data volume is reduced, to reach the purpose of level of detail management.The general block algorithm that adopts based on Octree of geological data Processing Algorithm based on LOD, this basic idea is to utilize the threshold value of a distance to control the degree of depth of Octree recursive operation, when this threshold ratio is big, obtain less earthquake voxel, otherwise then obtain more earthquake voxel.

Wherein, volume data is continued to carry out adaptive three dimensions decompose, make whole volume data be divided into a plurality of fritters, abbreviate sub-piece as, and it is standby to deposit each sub-blocks of data in disk.

Particularly, from the initial body data, adopt the structure of Octree volume data to be carried out 2 * 2 * 2 decomposition, obtaining little sub-piece, and make the corresponding corresponding volume data of each node of sub-piece, and, from the Octree root, whole volume data is divided into 8 sub-pieces, and when satisfying following default lasting decomposition condition, the antithetical phrase piece continues segmentation.

Default lasting decomposition condition comprises:

In the group piece consistance of voxel less than the defined critical value of user;

Wherein, the consistance described in the sub-piece refers to the similarity of voxel included in the sub-piece.

2. the size of sub-piece is greater than the size (for example, in the present embodiment, the video memory of the video card that uses is 4GB) of video memory;

3. the size of sub-piece is greater than the size of the defined boy's piece of user.

Need to prove that for fear of the jump phenomena between the adjacent sub-blocks, when volume data being carried out the self-adaptation decomposition, the data on the sub-block boundary have to be included in a plurality of sub-pieces, to guarantee the continuity in data.

203, all sub-pieces are carried out filtering and obtain visible sub-piece;

Wherein, this step specifically comprises:

203-1, establishment voxel count table;

The voxel count table is added up the quantity through visible voxel behind the transparent transmission function category.The voxel count table is used to create the visible sub-piece in filter back, jumps with the space of the empty volume elements of effective realization.A plain count table of given sub-block is constructed as follows:

Particularly, create the array of an one dimension, voxel counter by a sequence is formed, the index of this array is corresponding with the tonal range of voxel, travel through all voxels in the sub-piece, record the gray-scale value of all voxels in this sub-piece, and be increased in the corresponding voxel counter, at last all voxel counters are carried out comprehensively in a voxel count table.

203-2, obtain the ratio of the visible voxel in all sub-pieces according to current transparent transmission function calculation;

Particularly, constructed after the voxel count table, to working as previous given transparent transmission function, we can pass through formula (1) and calculate R1 and R2, obtain the visible voxel ratio in the sub-piece.

${\mathrm{<figure-callout\; id="1"\; label="R"\; filenames="HDA0000142663210000032.png"\; state="\{\{state\}\}">R</figure-callout>}}_1=\frac{\mathrm{VCT}\left[{\mathrm{OTF}}_{\mathrm{MAX}}\right]}{\mathrm{VCT}\left[{\mathrm{VCT}}_{\mathrm{MAX}}\right]}$ and ${\mathrm{<figure-callout\; id="2"\; label="R"\; filenames="HDA0000142663210000031.png,HDA0000142663210000032.png"\; state="\{\{state\}\}">R</figure-callout>}}_2=\frac{\mathrm{VCT}\left[{\mathrm{OTF}}_{\mathrm{MAX}}\right]}{\mathrm{VCT}\left[{\mathrm{VCT}}_{\mathrm{MAX}}\right]}---\left(1\right)$

Wherein, OTF _MAXAnd OTF _MINRepresent the minimum and maximum gray-scale value in the transparent transmission function, VCT _MAXThe maximal value of expression voxel count table.

203-3, the ratio by the visible voxel in each sub-piece are judged and are obtained visible sub-piece.

Particularly, the R that in step 203-2, calculates ₁Greater than α or R ₂Less than β, then should add the visible sub-piece in filter back by sub-piece.Wherein, α and β are user-defined 2 threshold values, carry out filtering for the antithetical phrase piece and obtain visible sub-piece.

For example, the gray-scale value in all voxels is greater than OTF _MAX, the R in the formula (1) ₁Be 0.Same, the gray-scale value in all voxels is less than OTF _MIN, R ₂Be 1.As a result, if R ₁Be 0 or R ₂Be 1, think that then specified sub-piece is sky.

To comprising the big volume data of many empty volume elements, we need determine which height piece is the dead zone according to the result of transparent transmission function category.And the sub-piece of the non-NULL that has only comprises the voxel that a small amount of needs are drawn.We utilize the voxel count table, together with user-defined 2 threshold alpha and β, these are almost empty sub-piece carry out filtering.We claim to be the visible sub-piece in filter back by the sub-piece that still remains behind the filtering operation.Each height piece is calculated the visibility information that is not subjected to the viewpoint constraint, and the step of the visible sub-piece in structure filter back is as follows:

1. utilize user-defined transparent transmission minimum of a function value (OTF _MIN) and maximal value (OTF _MAX) test all sub-pieces.

2. to the sub-piece of current accessed, if the R in the formula (1) ₁Greater than α or R ₂Less than β, should add the visible sub-piece in filter back by sub-piece.

3. when the transparent transmission function, upgrade filtered visible sub-piece according to the transparent tansfer function after the conversion.

204, with the data transmission of visible sub-piece in the video memory of video card, and carry out secondary and decompose being passed to data in the video memory according to the number of cores of video card;

In the present embodiment, the data of the visible sub-piece in the disk are read in internal memory, and be sent to video card by internal memory.To the number of cores among the GPU of the data based video card in the video card, it is carried out adaptive secondary decompose, for example, comprise 192 kernels among the video card GTX 460 that selects for use in the present embodiment, therefore the data in the video memory can be divided into 192 parts, distribute to 192 kernels and do parallel computation.

205, in the kernel of described video card, described secondary is decomposed the data that obtain and carry out the multi-resolution models that parallel computation obtains 3D seismic data;

Need to prove that the function that moves at GPU is called nuclear (Kernel), it is characterized in that operating all elements on a plurality of streams and be not only and operate independently element that executive routine is when calling the nuclear program with asynchronous system.Serial part in the executive routine is carried out at CPU, and nuclear is then only carried out at GPU as parallel section, and at this moment, the procedure division that GPU is gone up operation becomes grid (Grid), piece (Block), three thread ranks of thread (Thread) to implement respectively.According to the framework of GPU 3D seismic data being done self-adaptation decomposes, wherein, whole geological data is corresponding with grid (Grid) among the GPU, and sub-piece geological data is corresponding with piece (Block) among the GPU, and the slice of data in the sub-piece geological data is corresponding with thread (Thread) among the GPU.

Particularly, in the kernel of described video card, described secondary is decomposed the data that obtain and carries out the multi-resolution models that parallel computation obtains 3D seismic data, specifically comprise following process:

205-1, each sub-piece is carried out 3 D wavelet decompose the level of creating self correspondence with thicker resolution;

Particularly, before drawing, whole volume data is subdivided into a plurality of little sub-pieces.Each height piece is created the level of self with thicker resolution by the 3 D wavelet decomposition.The one-dimensional data that 4 levels are arranged is as shown in Figure 3 expressed, and at each level, the size of one-dimensional data reduces with 1/2 times speed, and the data on the identical layer second son block boundary have to be included in a plurality of sub-pieces, to guarantee the continuity of data interpolating.Wherein, raw data represents that with the 0th level the interval width of k layer data is 2 times of k-1 layer data interval width, and the quantity of data has only original 1/2.

The data value of 205-2, modification adjacent sub-blocks boundary is to keep the continuity of level transition;

Need to prove that if the schichtenaufbau of volume data as mentioned above, at same level adjacent block boundary, the distortion phenomenon that interpolation is brought just can not produce, this is because the voxel of boundary is shared by adjacent sub-piece.But this can not solve the distortion phenomenon of different levels adjacent block boundary.Therefore, in order to satisfy successional requirement, we have designed the method for transition between different levels, revise the data value of adjacent sub-blocks boundary.

Concrete method is, at the adjacent sub-blocks boundary, directly copies the data value of low level to last layer, or after the data value of low level carried out interpolation, replaces the data value of last layer time boundary.The self-adaptation of the implementation procedure adjacent sub-blocks with 4 levels as shown in Figure 4 is denoted as the data value of black in the level 0,1,2, be modified by copy or interpolation operation, to guarantee the continuity between the different levels adjacent sub-blocks.

205-3, carry out opacity and proofread and correct by with different resolution same sub-piece being carried out resampling antithetical phrase piece;

When drawing sub-piece under different levels, the opacity of sub-piece is different.Existing rendering algorithm depends on along radiation direction each pixel is sampled, but under the situation of multiresolution, volume data is resampled with resolution in a different manner, in order to preserve the optical property between the different resolution piece, need the raw data resampling, produce new sub-piece, transition function is also wanted therefore to make amendment.

Describe with different resolution same sub-piece being carried out resampling by example shown in Figure 5, wherein, each sampled point B _iWith color value c _iWith opacity value α _iCorresponding, then color value A and A ' are:

A＝α ₀c ₀+(1-α ₀)α ₁c ₁+(1-α ₀)(1-α ₁)c ₂ (2)

A′＝α′ ₀c ₀+(1-α′ ₀)c′ ₂ (3)

Here c ₂Be by sampling B ₀, B ₁, the input color that B2 obtains, C ' ₂Be by sampling B ₀, the input color that B2 obtains.Calculate opacity D and D ' that all accumulations are got up, obtain

D＝α ₀+(1-α ₀)α ₁+(1-α ₀)(1-α ₁)D ₂ (4)

D′＝α′ ₀+(1-α′ ₀)D′ ₂ (5)

Suppose to pile up opacity value and the idol sampling of getting up and equate that it is represented as D ₂=D ' ₂And D=D ', then

α ₀+(1-α ₀)α ₁+(1-α ₀)(1-α ₁)D ₂＝α′ ₀+(1-α′ ₀)D′ ₂ (6)

Be α ' ₀Separate this formula, we can obtain

α′ ₀＝1-(1-α ₀)(1-α ₁) (7)

By hypothesis α ₁=α ₀+ ε (wherein, ε is a very little number) can obtain following equation:

α′ ₀＝1-(1-α ₀) ²+O(ε) (8)

Thereby we revise transition function and are:

α′＝1-(1-α ₀) ² (9)

The level of the multiresolution correspondence of 205-4, antithetical phrase piece carries out adaptively selected;

Particularly, before the expression of multi-resolution models was configured, the size of each height piece was determined, and adopted following standard to carry out the selection of level:

Maximum opaque value: the maximum opaque value in sub-piece is determined by the maximum voxel of opaque value in this sub-piece.Its basic thought is exactly should carry out meticulous drafting to opaque big zone.

And the distance between the viewpoint: distance refers to the center of sub-piece and the distance between the viewpoint here.To from the nearer sub-piece of viewpoint, they are even more important, and we should carry out meticulous drafting.For other sub-piece, they from viewpoint away from, should draw with lower resolution.

View field: view field is determined by the border of sub-piece that mainly to those sub-pieces with big view field, we should compare meticulous drafting.

Stare distance: this parameter is to very useful based on the drafting of staring, and staring distance is the distance at the view field center of gaze area center and sub-piece.To those with stare nearer zone, center, we guarantee higher picture quality, and to those away from the zone of staring the center, lower resolution is drawn.

205-6, by calculating the multi-resolution models of 3D seismic data;

Wherein, calculate the following formula of the concrete use of multi-resolution models of 3D seismic data:

$f(x,y,z)=\underset{i,j,k\&Element;z}{\mathrm{\&Sigma;}}{\left(A_m^{d}f\right)}_{i,j,k}{\mathrm{\&phi;}}_{m,i,j,k}(x,y,z)+\underset{m=1}{\overset{J}{\mathrm{\&Sigma;}}}\underset{n=1}{\overset{7}{\mathrm{\&Sigma;}}}\underset{i,j,k\&Element;z}{\mathrm{\&Sigma;}}{\left(D_m^{n}f\right)}_{i,j,k}{\mathrm{\&psi;}}_{m,i,j,k}^n(x,y,z)---\left(10\right)$

In the formula, integer J is the progression that the wavelet multiresolution rate is analyzed, coefficient

Be that volume data dispersing under the J class resolution ratio smoothly approached coefficient

Be the discrete detail signal of volume data under the J class resolution ratio.

206, establishment is based on the volume drawing equation of wavelet field;

Particularly, the ultimate principle of volume drawing is to consider light intensity and the transparency of each voxel in the volume data, and directly the intensity signal fusion with all voxels projects on the picture plane.In the present embodiment, set t and be the vector of looking of expression direction of visual lines, u, v be with it quadrature look like two orthogonal vectors in the plane, then the optical model of volume drawing can be with representing along the curvilinear integral of direction of visual lines t, and is specific as follows:

I _t(u，v，d)＝I _backexp[-∫ _t(0，d)τ(s)dt]+∫ _t(0，d)I(s)exp[-∫ _t(0，d)τ(s)dt]dt (10)

s＝uu+vv+tt (11)

$T({\mathrm{<figure-callout\; id="1"\; label="t"\; filenames="HDA0000142663210000032.png"\; state="\{\{state\}\}">t</figure-callout>}}_1,t_2)=\exp [-{\&Integral;}_{t({\mathrm{<figure-callout\; id="1"\; label="t"\; filenames="HDA0000142663210000032.png"\; state="\{\{state\}\}">t</figure-callout>}}_1,t_2)}\mathrm{\&tau;}\left(s\right)\mathrm{dt}]---\left(12\right)$

α(t ₁，t ₂)＝1-T(t ₁，t ₂)(13)

In formula (10)-(13), I _t(d) expression is merged the comprehensive light intensity that projects on the picture plane along vector t direction for u, v, and d is the orthogonal distance of picture plane and background, I _BackBe the background light intensity, τ (s) is the absorption coefficient of light of s place voxel, and I (s) is the light intensity of s place voxel, T (t ₁, t ₂) be that direction of visual lines t goes up t1 to the comprehensive transparency at t2 place, α (t ₁, t ₂) then be corresponding opacity, symbol t (t ₁, t ₂) represent along the path of integration of vector t direction from t1 to t2.

With volume data f (wavelet field multiresolution expression formula substitution formula (10) z) can obtain the volume drawing equation of following wavelet field for x, y:

$I_t(u,v,d)=\underset{i,j,k\&Element;z}{\mathrm{\&Sigma;}}{\left(A_m^{d}f\right)}_{i,j,k}{\&Integral;}_{t(0,d)}{\mathrm{\&phi;}}_{m,i,j,k}\left(s\right)\&times;\exp [-{\&Integral;}_{t(t,d)}\mathrm{\&tau;}\left(s\right)\mathrm{dt}]$

$+\underset{m=1}{\overset{J}{\mathrm{\&Sigma;}}}\underset{n=1}{\overset{7}{\mathrm{\&Sigma;}}}\underset{i,j,k\&Element;z}{\mathrm{\&Sigma;}}{\left(D_m^{n}f\right)}_{i,j,k}{\&Integral;}_{t(0,d)}{\mathrm{\&psi;}}_{m,i,j,k}^n\left(s\right)\&times;\exp [-{\&Integral;}_{t(t,d)}\mathrm{\&tau;}\left(s\right)\mathrm{dt}]\mathrm{dt}---\left(15\right)$

207, the discrete approximation signal under the transmission volume data J resolution, again by resolution (J, J-1 ..., 2,1) order transmit the detail signal of volume data step by step;

208, according to the direction of visual lines vector, according to the precedence that the multiresolution volume data arrives, utilize the volume drawing equation based on wavelet field that volume data is carried out by thick and smart 3 D rendering.

Wherein, this step specifically comprises:

208-1, being orthogonal on the picture plane of looking vector, draw out the 3-D view of low resolution with the discrete approximation signal under the J resolution;

208-2, with detail signal by resolution (J, J-1 ..., 2,1) refined image progressively, generate successively resolution (J, J-1 ..., 2,1) 3-D view (smoothly approaching).

Wherein, with detail signal by resolution (J, J-1 ..., 2,1) and refined image progressively, specifically realize by applying three-dimensional small echo Mal lat reconstruction filter group.

Further, when changing the direction of visual lines vector, repeated execution of steps 208.

The method of a kind of visual earthquake data that present embodiment provides, by using among the CPU geological data is converted to volume data, and volume data carried out piecemeal, filtration obtains visible sub-piece, transfer to GPU and carry out image and draw in that visible sub-piece is passed, can reduce the calculated amount of GPU greatly, make the three-dimensional visualization that on common computer, also can carry out geological data by common GPU, easy to use, and save cost.

Embodiment 3

As shown in Figure 6, present embodiment provides a kind of system of visual earthquake data, and this system comprises: central processing unit 3 and image processor 4;

Wherein, as shown in Figure 7, central processing unit 3 comprises:

Modular converter 301 is used for 3D seismic data is converted to volume data;

Decomposing module 302 is used for volume data is decomposed into sub-piece;

Filtering module 303 is used for that sub-piece is carried out filtering and obtains visible sub-piece;

As shown in Figure 8, graphic process unit 4 comprises:

Model computation module 401, the multi-resolution models for visible sub-piece being carried out obtain after secondary decomposes 3D seismic data wherein, comprises discrete approximation signal and detail signal in the multi-resolution models;

Transport module 402 is used for the discrete approximation signal under the transmission J class resolution ratio, and the order that is begun to successively decrease according to resolution levels by the J class resolution ratio transmits detail signal step by step again, and wherein, J is integer;

Drafting module 403 is used for according to the direction of visual lines vector, and the order that arrives according to discrete approximation signal and detail signal carries out 3 D rendering.

Further, modular converter 301, concrete being used for will be converted to a plurality of voxels at the 3D seismic data that a plurality of sampled points obtain, and a plurality of voxel combinations are obtained volume data;

Wherein, comprise 3D seismic data, rgb color value and the transparence information that each voxel corresponding sampling points obtains in each voxel at least.

Further, decomposing module 302, the concrete structure that is used for the employing Octree is carried out 2 * 2 * 2 decomposition to volume data, obtains 8 sub-pieces, and when satisfying first when pre-conditioned, the sub-piece in 8 sub-pieces that continue decomposition is obtained decomposes and obtains littler sub-piece;

Wherein, first pre-conditioned the comprising:

The consistance of the voxel in the sub-piece less than the size of user-defined critical value, sub-piece greater than with the size of the size of video memory, the sub-piece size greater than user-defined boy's piece.

Further, as shown in Figure 9, filtering module 303 comprises:

Creating unit 3031 is used for creating the voxel count table;

Ratio computing unit 3032 is for the ratio that obtains the visible voxel of all sub-pieces according to the transparent transmission function calculation;

Judging unit 3033, the ratio that is used for the visible voxel by all sub-pieces are judged and are obtained visible sub-piece;

Wherein, the voxel count table is used for the quantity of the visible voxel behind the statistics process transparent transmission function category, and the transparent transmission function obtains according to the direction of visual lines vector.

Further, as shown in figure 10, central processing unit 3 also comprises:

Memory module 304 is used for before model computation module 401 carries out visible sub-piece to obtain the multiresolution module of 3D seismic data after secondary decomposes visible sub-piece being transferred in the video memory.

Further, as shown in figure 11, model computation module 401 comprises:

3 D wavelet resolving cell 4011 is used for that each visible sub-piece is carried out 3 D wavelet and decomposes the level that the order of successively decreasing with resolution levels is created visible sub-piece correspondence;

Revise unit 4012, be used for revising the data value at two adjacent visible sub-block boundary places, with the continuity of the level transition that keeps each resolution levels correspondence;

Correcting unit 4013 is used for carrying out the opacity correction by with the resolution of different stage same visible sub-piece being carried out resampling antithetical phrase piece;

Selected cell 4014 is used for carrying out adaptively selected to the level of the multiresolution correspondence of visible sub-piece;

Model computing unit 4015 is for the multi-resolution models that calculates 3D seismic data;

Wherein, comprise discrete approximation signal and detail signal in the above-mentioned multi-resolution models.

Further, revise unit 4012, concrete being used at two adjacent visible sub-block boundary places, copy the data value of current level to last one high-level level, or after the data value of current level carried out interpolation, replace the data value of last one high-level level boundary.

Further, above-mentioned multi-resolution models is specially:

$f(x,y,z)=\underset{i,j,k\&Element;z}{\mathrm{\&Sigma;}}{\left(A_m^{d}f\right)}_{i,j,k}{\mathrm{\&phi;}}_{m,i,j,k}(x,y,z)+\underset{m=1}{\overset{J}{\mathrm{\&Sigma;}}}\underset{n=1}{\overset{7}{\mathrm{\&Sigma;}}}\underset{i,j,k\&Element;z}{\mathrm{\&Sigma;}}{\left(D_m^{n}f\right)}_{i,j,k}{\mathrm{\&psi;}}_{m,i,j,k}^n(x,y,z)$

Wherein, f (x, y z) are volume data, and integer J is the progression that the wavelet multiresolution rate is analyzed,

Be the discrete level and smooth approximation signal of volume data under the J class resolution ratio,

Be the discrete detail signal of volume data under the J class resolution ratio.

Further, as shown in figure 12, drafting module 403 comprises:

First drawing unit 4031 is used for drawing out the 3-D view of default resolution with the discrete approximation signal under the J resolution on the picture plane that is orthogonal to the direction of visual lines vector;

Second drawing unit 4032 is for the order details of use signal that is begun by J resolution to successively decrease according to resolution refinement 3-D view successively.

Further, as shown in figure 13, image processor 4 also comprises:

Redraw molding 404, be used for after the direction of visual lines vector changes, according to the direction of visual lines vector after changing, the order that arrives according to discrete approximation signal and detail signal carries out 3 D rendering.

The system of a kind of visual earthquake data that present embodiment provides, by using among the CPU geological data is converted to volume data, and volume data carried out piecemeal, filtration obtains visible sub-piece, transfer to GPU and carry out image and draw in that visible sub-piece is passed, can reduce the calculated amount of GPU greatly, make the three-dimensional visualization that on common computer, also can carry out geological data by common GPU, easy to use, and save cost.

The all or part of step that one of ordinary skill in the art will appreciate that realization above-described embodiment can be finished by hardware, also can instruct relevant hardware to finish by program, described program can be stored in a kind of computer-readable recording medium, the above-mentioned storage medium of mentioning can be ROM (read-only memory), disk or CD etc.

The above only is preferred embodiment of the present invention, and is in order to limit the present invention, within the spirit and principles in the present invention not all, any modification of doing, is equal to replacement, improvement etc., all should be included within protection scope of the present invention.

Claims (20)

1. the method for visual earthquake data is characterized in that, described method comprises:

By central processor CPU 3D seismic data is converted to volume data, and described volume data is decomposed into sub-piece;

Described CPU carries out filtering with described sub-piece and obtains visible sub-piece;

Graphic process unit GPU carries out obtaining after secondary decomposes the multi-resolution models of described 3D seismic data with described visible sub-piece, wherein, comprises discrete approximation signal and detail signal in the described multi-resolution models;

Transmit the described discrete approximation signal under the J class resolution ratio, the order that is begun to successively decrease according to resolution levels by the J class resolution ratio transmits described detail signal step by step again, and wherein, J is integer;

According to the direction of visual lines vector, the order that arrives according to described discrete approximation signal and detail signal carries out 3 D rendering.

2. method according to claim 1 is characterized in that, described 3D seismic data is converted to volume data, comprising:

To be converted to a plurality of voxels at the 3D seismic data that a plurality of sampled points obtain, and described a plurality of voxel combinations will be obtained described volume data;

Wherein, comprise 3D seismic data, rgb color value and the transparence information that described each voxel corresponding sampling points obtains in each voxel at least.

3. method according to claim 1 is characterized in that, described volume data is decomposed into sub-piece, comprising:

Adopt the structure of Octree that described volume data is carried out 2 * 2 * 2 decomposition, obtain 8 sub-pieces, and when satisfying first when pre-conditioned, the sub-piece in lasting described 8 sub-pieces that decomposition is obtained decomposes and obtains littler sub-piece;

Wherein, described first pre-conditioned the comprising:

The consistance of the voxel in the described sub-piece less than the size of user-defined critical value, described sub-piece greater than the size of the size of the video memory corresponding with described GPU, the described sub-piece size greater than described user-defined boy's piece.

4. method according to claim 1 is characterized in that, describedly described sub-piece is carried out filtering obtains visible sub-piece, comprising:

Create the voxel count table;

Obtain the ratio of the visible voxel in all sub-pieces according to the transparent transmission function calculation;

Ratio judgement by the visible voxel in described all sub-pieces obtains described visible sub-piece;

Wherein, described voxel count table is used for statistics through the quantity of the visible voxel behind the described transparent transmission function category, and described transparent transmission function obtains according to described direction of visual lines vector.

5. method according to claim 1 is characterized in that, described GPU carries out described visible sub-piece to obtain after secondary decomposes before the multiresolution module of described 3D seismic data, and described method also comprises:

Described visible sub-piece is transferred in the video memory of described GPU correspondence.

6. method according to claim 1 is characterized in that, described GPU carries out obtaining after secondary decomposes the multi-resolution models of described 3D seismic data with described visible sub-piece, comprising:

Each visible sub-piece is carried out 3 D wavelet decompose the level that the order of successively decreasing with resolution levels is created described visible sub-piece correspondence;

Revise the data value at two adjacent described visible sub-block boundary places, with the continuity of the level transition that keeps each resolution levels correspondence;

Carry out the opacity correction by with the resolution of different stage same described visible sub-piece being carried out resampling antithetical phrase piece;

Level to the multiresolution correspondence of described visible sub-piece carries out adaptively selected;

Calculate the multi-resolution models of described 3D seismic data;

Wherein, comprise discrete approximation signal and detail signal in the described multi-resolution models.

7. method according to claim 6 is characterized in that, the data value at two described visible sub-block boundary places that described modification is adjacent to keep the continuity of level transition, comprising:

At two adjacent described visible sub-block boundary places, copy the data value of current level to last one high-level level, or after the data value of current level carried out interpolation, replace the data value of last one high-level level boundary.

8. method according to claim 6 is characterized in that, described multi-resolution models is specially:

$f(x,y,z)=\underset{i,j,k\&Element;z}{\mathrm{\&Sigma;}}{\left(A_m^{d}f\right)}_{i,j,k}{\mathrm{\&phi;}}_{m,i,j,k}(x,y,z)+\underset{m=1}{\overset{J}{\mathrm{\&Sigma;}}}\underset{n=1}{\overset{7}{\mathrm{\&Sigma;}}}\underset{i,j,k\&Element;z}{\mathrm{\&Sigma;}}{\left(D_m^{n}f\right)}_{i,j,k}{\mathrm{\&psi;}}_{m,i,j,k}^n(x,y,z)$

Wherein, f (x, y z) are described volume data, and integer J is the progression that the wavelet multiresolution rate is analyzed,

Be the discrete level and smooth approximation signal of described volume data under the J class resolution ratio,

Be the discrete detail signal of described volume data under the J class resolution ratio.

9. method according to claim 1 is characterized in that, and is described according to the direction of visual lines vector, and the order that arrives according to described discrete approximation signal and detail signal carries out 3 D rendering, comprising:

Be orthogonal on the picture plane of described direction of visual lines vector, drawing out the 3-D view of default resolution with the discrete approximation signal under the described J resolution;

Begin the order details of use signal that successively decreases according to the resolution described 3-D view of refinement successively by J resolution again.

10. method according to claim 1 is characterized in that, after described direction of visual lines vector changed, described method also comprised:

According to the direction of visual lines vector after the described change, the order that arrives according to described discrete approximation signal and detail signal carries out 3 D rendering.

11. the system of visual earthquake data is characterized in that, described system comprises: central processing unit and image processor;

Wherein, described central processing unit comprises:

Modular converter is used for 3D seismic data is converted to volume data;

Decomposing module is used for described volume data is decomposed into sub-piece;

Filtering module is used for that described sub-piece is carried out filtering and obtains visible sub-piece;

Described graphic process unit comprises:

Model computation module, the multi-resolution models for described visible sub-piece being carried out obtain after secondary decomposes described 3D seismic data wherein, comprises discrete approximation signal and detail signal in the described multi-resolution models;

Transport module is used for the described discrete approximation signal under the transmission J class resolution ratio, and the order that is begun to successively decrease according to resolution levels by the J class resolution ratio transmits described detail signal step by step again, and wherein, J is integer;

Drafting module is used for according to the direction of visual lines vector, and the order that arrives according to described discrete approximation signal and detail signal carries out 3 D rendering.

12. system according to claim 11 is characterized in that, described modular converter, and concrete being used for will be converted to a plurality of voxels at the 3D seismic data that a plurality of sampled points obtain, and described a plurality of voxel combinations are obtained described volume data;

Wherein, comprise 3D seismic data, rgb color value and the transparence information that described each voxel corresponding sampling points obtains in each voxel at least.

13. system according to claim 11, it is characterized in that, described decomposing module, the concrete structure that is used for the employing Octree is carried out 2 * 2 * 2 decomposition to described volume data, obtain 8 sub-pieces, and when satisfying first when pre-conditioned, the sub-piece in described 8 the sub-pieces that continue decomposition is obtained decomposes and obtains littler sub-piece;

Wherein, described first pre-conditioned the comprising:

The consistance of the voxel in the described sub-piece less than the size of user-defined critical value, described sub-piece greater than with the size of the size of video memory, the described sub-piece size greater than described user-defined boy's piece.

14. system according to claim 11 is characterized in that, described filtering module comprises:

Creating unit is used for creating the voxel count table;

The ratio computing unit is for the ratio that obtains the visible voxel of all sub-pieces according to the transparent transmission function calculation;

Judging unit, the ratio that is used for the visible voxel by described all sub-pieces are judged and are obtained described visible sub-piece;

Wherein, described voxel count table is used for statistics through the quantity of the visible voxel behind the described transparent transmission function category, and described transparent transmission function obtains according to described direction of visual lines vector.

15. system according to claim 11 is characterized in that, described central processing unit also comprises:

Memory module is used for before described model computation module carries out described visible sub-piece to obtain the multiresolution module of described 3D seismic data after secondary decomposes described visible sub-piece being transferred in the described video memory.

16. system according to claim 11 is characterized in that, described model computation module comprises:

The 3 D wavelet resolving cell is used for that each visible sub-piece is carried out 3 D wavelet and decomposes the level that the order of successively decreasing with resolution levels is created described visible sub-piece correspondence;

Revise the unit, be used for revising the data value at two adjacent described visible sub-block boundary places, with the continuity of the level transition that keeps each resolution levels correspondence;

Correcting unit is used for carrying out the opacity correction by with the resolution of different stage same described visible sub-piece being carried out resampling antithetical phrase piece;

Selected cell is used for carrying out adaptively selected to the level of the multiresolution correspondence of described visible sub-piece;

The model computing unit is for the multi-resolution models that calculates described 3D seismic data;

Wherein, comprise discrete approximation signal and detail signal in the described multi-resolution models.

17. system according to claim 16, it is characterized in that, described modification unit, concrete being used at two adjacent described visible sub-block boundary places, copy the data value of current level to last one high-level level, or after the data value of current level carried out interpolation, replace the data value of last one high-level level boundary.

18. system according to claim 16 is characterized in that, described multi-resolution models is specially:

$f(x,y,z)=\underset{i,j,k\&Element;z}{\mathrm{\&Sigma;}}{\left(A_m^{d}f\right)}_{i,j,k}{\mathrm{\&phi;}}_{m,i,j,k}(x,y,z)+\underset{m=1}{\overset{J}{\mathrm{\&Sigma;}}}\underset{n=1}{\overset{7}{\mathrm{\&Sigma;}}}\underset{i,j,k\&Element;z}{\mathrm{\&Sigma;}}{\left(D_m^{n}f\right)}_{i,j,k}{\mathrm{\&psi;}}_{m,i,j,k}^n(x,y,z)$

Wherein, f (x, y z) are described volume data, and integer J is the progression that the wavelet multiresolution rate is analyzed,

Be the discrete level and smooth approximation signal of described volume data under the J class resolution ratio,

Be the discrete detail signal of described volume data under the J class resolution ratio.

19. system according to claim 11 is characterized in that, described drafting module comprises:

First drawing unit is used for drawing out the 3-D view of default resolution with the discrete approximation signal under the described J resolution on the picture plane that is orthogonal to described direction of visual lines vector;

Second drawing unit is for the order details of use signal that is begun by J resolution to successively decrease according to the resolution described 3-D view of refinement successively.

20. system according to claim 11 is characterized in that, described image processor also comprises:

Redraw the molding piece, be used for after described direction of visual lines vector changes, according to the direction of visual lines vector after the described change, the order that arrives according to described discrete approximation signal and detail signal carries out 3 D rendering.

Images