CN104407384B - Seismic section display method - Google Patents
Seismic section display method Download PDFInfo
- Publication number
- CN104407384B CN104407384B CN201410743516.5A CN201410743516A CN104407384B CN 104407384 B CN104407384 B CN 104407384B CN 201410743516 A CN201410743516 A CN 201410743516A CN 104407384 B CN104407384 B CN 104407384B
- Authority
- CN
- China
- Prior art keywords
- pixel
- texture
- seismic profile
- screen
- seismic
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Landscapes
- Image Generation (AREA)
Abstract
The invention provides a seismic section display method. The seismic section display method comprises the following steps: for each pixel of a screen, mapping the covering range of the pixel from a pixel space to a gather space to obtain a seismic trace and a sampling point, on which the pixel covers; correspondingly generating ripple images positioned in each pixel according to the seismic trace and the sampling point, which the pixel covers; obtaining the brightness contribution value of each ripple image in each pixel to the pixel, and accumulating the brightness contribution values of all the ripple images in each pixel to the pixel to obtain the brightness value of the pixel; caching the brightness values of each pixel in a video memory; drawing the brightness values of each pixel cached in the video memory on a screen display. Aliasing phenomena caused by trace decimation are avoided.
Description
Technical field
The present invention relates to image processing of seismic section technical field, especially relate to a kind of method that seismic profile shows.
Background technology
The visualization of geological data especially stacked section data is very important supplementary meanss in seismic data process,
After the geological data of field acquisition passes through and processes, need to be shown in the form of images by means of computer technology
Come, thus for data-processing operator, geological analysis personnel provide a kind of visual basis for estimation, improve oil exploration precision and
Efficiency.
Geological data is all discrete sampling, and is arranged according to road, a series of dispersion number actually in space
Value, its visual principal mode is to be drawn out data in the form of wave molding, variable density.When being drawn in wave molding mode
When, traditional method for drafting is to couple together discrete sampled point straight line, forms a wavy curve, in order to explicit
Earth's surface shows the positive and negative of amplitude of wave form, needs some direction to seismic wave (usually positive) to be filled with, so in earthquake
More visible reflection line-ups are just had at reflecting section.When user's downscaled images, because active computer screen is differentiated
The restriction of rate, after reducing to a certain extent, the distance between adjacent two wave moldings already less than a pixel, even if pressing
Drawn according to a pixel, also resulting in whole screen is whole blacking.General solution is that geological data is entered
Hang Chou road is processed, and extracts a seismic channel every n bar seismic channel and is drawn as representative.Although this method realizes letter
Single, but alias can be produced, and this phenomenon can explain personnel to the overall assurance of geological condition by severe jamming.
Content of the invention
It is an object of the invention to provide a kind of method that seismic profile shows, reduce to go out during display to solve seismic profile
Existing fake frequency phenomenon.
For reaching above-mentioned purpose, the invention provides a kind of method that seismic profile shows, comprise the following steps:
(s1), for each pixel of screen, its coverage is mapped to road collection space from pixel space, obtains each
Seismic channel and sampled point that pixel is covered；
(s2) seismic channel, being covered according to each pixel and sampled point is corresponding generates the moire pattern being located in each pixel
Picture；
(s3), obtain each moire pattern picture in each pixel brightness contribution value to this pixel, and by each picture described
Each moire pattern picture in element the brightness contribution value of this pixel is added up after as this pixel brightness value；
(s4), the brightness value of each pixel is cached in video memory；
(s5), the brightness value of each pixel of caching in described video memory is drawn to screen display.
The method that the seismic profile of the present invention shows, step (s1) to (s4) is by graphic process unit (graphics
Processing unit, gpu) texture generate thread complete；Step (s5) is completed by gpu texture rendering thread.
The method that the seismic profile of the present invention shows, when receiving image scaling instruction, described gpu texture rendering thread will
The currently displaying seismic profile of described screen carries out corresponding scaling and processes the interim seismic profile of generation output display, and simultaneously
Trigger described gpu texture rendering thread, when described gpu texture rendering thread generates seismic profile, described screen is currently shown
The interim seismic profile showing is substituted for the seismic profile that described gpu texture rendering thread generates.
The method that the seismic profile of the present invention shows, described screen is covered by the texture cell of several fixed sizes,
Each texture cell described is configured with an independent texture and generates thread, and each independent texture generates thread parallel and generates
The seismic profile texture of this texture cell.
The method that the seismic profile of the present invention shows, the edge of adjacent texture cell is overlapped.
The method that the seismic profile of the present invention shows, overlapped 2 pixels in border of adjacent texture cell.
The method that the seismic profile of the present invention shows, described by described video memory in the brightness value of each pixel of caching paint
Make to screen display, adjacent texture cell reduces by a pixel at border overlay.
The method that the seismic profile of the present invention shows, described each pixel for screen, by its coverage from pixel
Space reflection, to road collection space, specifically includes:
Each pixel is transformed into world coordinates from screen coordinate: w=p-1*m-1* s, wherein p-1For projection matrix against square
Battle array, m-1For model matrix inverse matrix, s is the screen coordinate of pixel, and w is the corresponding world coordinates of pixel.
The method that the seismic profile of the present invention shows, each moire pattern in each pixel of described acquisition is as to this pixel
Brightness contribution value, and as this picture after each moire pattern picture in each pixel described is added up to the brightness contribution value of this pixel
The brightness value of element, specifically includes:
The world coordinates w of each pixel is evenly dividing into some grid g, puts each grid giInitial value be 0；
Calculate the moire pattern picture intersecting with the world coordinates w of each pixel, and the coverage by these moire pattern pictures
Each grid g with this pixeliIt is test for intersection, and the grid g that will be cappediValue put 1；
Calculate the number of the grid that its value in each pixel is 1, and this number is set to the brightness value of current pixel
The method that the seismic profile of the present invention shows, the world coordinates w of each pixel is evenly dividing into 16 × 16 lattice
Sub- g.
The method that the seismic profile of the present invention shows, each texture cell covers 512 × 512 or 1024 × 1024 pictures
Element.
The method that the seismic profile of the embodiment of the present invention shows, when user reduce earthquake wave molding indication range make adjacent
The distance between twice are less than during a pixel it is not necessary to take out display, but the drafting of wave molding are gradually transformed into similar
In the drafting effect of variable density, so as to effectively avoid due to taking out the fake frequency phenomenon that display leads to.
Brief description
Accompanying drawing described herein is used for providing a further understanding of the present invention, constitutes the part of the application, not
Constitute limitation of the invention.In the accompanying drawings:
Fig. 1 is the method flow diagram that shows of seismic profile of the embodiment of the present invention；
Fig. 2 is that in the embodiment of the present invention, texture cell divides schematic diagram；
Fig. 3 is the schematic diagram that in the embodiment of the present invention, pixel is intersected with seismic channel；
Fig. 4 a-4d is seismic profile variable density display schematic diagram in the embodiment of the present invention.
Specific embodiment
For making the object, technical solutions and advantages of the present invention become more apparent, with reference to embodiment and accompanying drawing, to this
Invention is described in further details.Here, the schematic description and description of the present invention is used for explaining the present invention, but do not make
For limitation of the invention.
Below in conjunction with the accompanying drawings, the specific embodiment of the present invention is described in further detail.
With reference to shown in Fig. 1, the method that the seismic profile of the embodiment of the present invention shows comprises the following steps:
Step s1, each pixel for screen, its coverage is mapped to road collection space from pixel space, obtains every
Seismic channel and sampled point that individual pixel is covered；Specifically include:
Each pixel is transformed into world coordinates from screen coordinate: w=p-1*m-1* s, wherein p-1For projection matrix against square
Battle array, m-1For model matrix inverse matrix, s is the screen coordinate of pixel, and w is the corresponding world coordinates of pixel.
Step s2, the seismic channel being covered according to each pixel and sampled point is corresponding generates the ripple being located in each pixel
Image.For example the discrete sampled point straight line of each seismic channel is coupled together, form a wavy curve, in order to explicitly
Representing the positive and negative of amplitude of wave form, needing some direction to seismic wave (usually positive) to be filled with, thus forming ripple
Image, so just has more visible reflection line-ups at seismic reflection section.
Step s3, obtain each moire pattern picture in each pixel brightness contribution value to this pixel, and by each pixel
Each interior moire pattern picture the brightness contribution value of this pixel is added up after as this pixel brightness value；Specifically include:
The world coordinates w of each pixel is evenly dividing into some grid g (such as 16 × 16), puts each grid gi's
Initial value is 0；
Calculate the moire pattern picture (as shown in Figure 3) intersecting with the world coordinates w of each pixel, and by these moire pattern pictures
Coverage and this pixel each grid giIt is test for intersection, and the grid g that will be cappediValue put 1；
Calculate the number of the grid that its value in each pixel is 1, and this number is set to the brightness value of current pixel
Step s4, the brightness value of each pixel is cached in video memory.
Step s5, the brightness value of each pixel of caching in video memory is drawn to screen display.
Because the generation cost of texture image will be far longer than texture rendering to screen, in order to improve the response to user
Speed, in the method that the seismic profile of the embodiment of the present invention shows, will create two classes and draw threads: gpu texture generate thread and
Gpu texture rendering thread, and corresponding graphic plotting environment, step s1～s4 generates thread by gpu texture and completes, step
Rapid s5 is completed by gpu texture rendering thread.Gpu texture generates thread to be responsible for generating texture image in running background, and gpu texture is painted
Thread processed is responsible for interaction and the response to user of user, and texture image is plotted to screen.Instruct when receiving image scaling
When, seismic profile currently displaying for screen is carried out corresponding scaling and processes the interim seismic profile of generation by gpu texture rendering thread
And output display (is so direct to the response of user and quickly, so directly amplifying the interim seismic profile of generation can produce
Raw mosaic phenomenon), and trigger gpu texture rendering thread simultaneously, when gpu texture rendering thread generates seismic profile, will shield
The currently displaying interim seismic profile of curtain replaces the seismic profile that gpu texture rendering thread generates, so that the present invention is implemented
Example quickly can respond to user, be provided that the drafting preferable image of effect simultaneously again, improve Consumer's Experience.
Because usual video card can be constrained to the maximum magnitude of texture image, and the explanation of seismic profile generally requires
Multiple screen splicings show together, need very big texture, therefore in the embodiment of the present invention, in order to solve multi-screen
The super large being formed draws problem, need to be by the texture cell cover screen of several fixed sizes, and typically each texture cell covers
512 × 512 or 1024 × 1024 pixels, too big texture cell can lead to the waste drawn, and reduces system whole efficiency.
For example shown in Fig. 2, four texture cells a, b, c, d can be screened into, in figure solid line is screen, dotted line is to generate texture
Scope, chain-dotted line is the scope drawing texture.Each texture cell is configured with an independent texture and generates thread, and each is only
Vertical texture generates the seismic profile texture that thread parallel generates this texture cell, and it is permissible that multiple textures generate thread parallel process
The efficiency of texture generation is greatly improved.
During being drawn to screen display due to the brightness value of each pixel of caching in by video memory, the generation meeting of texture
Carry out smothing filtering in the form of bilinear interpolation, can produce between texture cell border and not connect between texture and pixel
Continuous, human eye to this discontinuously very sensitive, in the method that the seismic profile of the therefore embodiment of the present invention shows, adjacent texture
Need to carry out overlapped between unit.Such as between the adjacent texture cell of each two overlapping 2 pixels, texture rendering process
In, the adjacent texture cell of each two draws 1 pixel less in boundary, can effectively seamlessly transit at such Texture Boundaries and
Merge.The size of one texture cell of hypothesis is n*n pixel, needs to generate n*n pixel in texture generating process
The texture of size, but only need to draw the stricture of vagina in the individual pixel coverage of (n-1) * (n-1) during by texture toward screen drawing
Reason.
As can be seen here, the method that the seismic profile of the embodiment of the present invention shows proposes one kind and is based on texture partitioning caching mechanism
Multi-threaded parallel drafting framework, first wave molding is plotted in texture cache, then again texture is drawn directly into screen,
Such benefit is: after one time texture generates, can repeatedly draw, thus reducing the calculation cost of pixel, improve and painting
Efficiency processed.And, when user reduce earthquake wave molding indication range make the distance between adjacent twice be less than a pixel when,
Do not need to take out display, but the drafting of wave molding is gradually transformed into drafting effect similar to variable density (as Fig. 4 a-4d institute
The variable density gradual change shown), so as to effectively avoid due to taking out the fake frequency phenomenon that display leads to, and the embodiment of the present invention
Parallel piecemeal draw and can deal with large screen display, especially Multi-screen display system, thus eliminating video memory texture size limit
Make the constraint bringing.
Those skilled in the art are it will also be appreciated that the embodiment of the present invention various illustrative components, blocks, unit and the step listed
Suddenly can be by hardware, software or both be implemented in combination in.To realize depending on specific to then passing through hardware or software
Application and whole system design requirement.Those skilled in the art can be for every kind of specific application, it is possible to use various
Method realizes described function, but this realization is understood not to the scope beyond embodiment of the present invention protection.
Various illustrative logical block described in the embodiment of the present invention, or unit can pass through general processor,
Digital signal processor, special IC (asic), field programmable gate array or other programmable logic device, discrete gate
Or transistor logic, discrete hardware components, or the design of any of the above described combination are realizing or to operate described function.General place
Reason device can be microprocessor, and alternatively, this general processor can also be any traditional processor, controller, microcontroller
Device or state machine.Processor can also be realized by the combination of computing device, such as digital signal processor and microprocessor,
Multi-microprocessor, one or more microprocessors one Digital Signal Processor Core of joint, or any other like configuration
To realize.
The step of the method described in the embodiment of the present invention or algorithm can be directly embedded into hardware, computing device soft
Part module or the combination of both.Software module can be stored in ram memorizer, flash memory, rom memorizer, eprom storage
Other any form of storage media in device, eeprom memorizer, depositor, hard disk, moveable magnetic disc, cd-rom or this area
In.Exemplarily, storage medium can be connected with processor, so that processor can read information from storage medium, and
Write information can be deposited to storage medium.Alternatively, storage medium can also be integrated in processor.Processor and storage medium can
To be arranged in asic, asic can be arranged in user terminal.Alternatively, processor and storage medium can also be arranged at use
In different part in the terminal of family.
In one or more exemplary designs, above-mentioned functions described by the embodiment of the present invention can be in hardware, soft
The combination in any of part, firmware or this three is realizing.If realized in software, these functions can store and computer-readable
On medium, or it is transmitted on the medium of computer-readable with one or more instructions or code form.Computer readable medium includes electricity
Brain stores medium and is easy to so that allowing computer program transfer to the telecommunication media in other places from a place.Storage medium is permissible
It is any general or special computer can be with the useable medium of access.For example, such computer readable media can include but
It is not limited to ram, rom, eeprom, cd-rom or other optical disc storage, disk storage or other magnetic storage device, or other
What can be used for carrying or storing to instruct or data structure and other can be by general or special computer or general or special handlings
Device reads the medium of the program code of form.Additionally, any connection can be properly termed computer readable medium, example
As if software is to pass through a coaxial cable, fiber optic cables, double from a web-site, server or other remote resource
Twisted wire, Digital Subscriber Line (dsl) or with the wireless way for transmittings such as example infrared, wireless and microwave be also contained in defined
In computer readable medium.Described video disc (disk) and disk (disc) include Zip disk, radium-shine disk, CD, dvd, floppy disk
And Blu-ray Disc, disk is generally with magnetic duplication data, and video disc generally carries out optical reproduction data with laser.Combinations of the above
Can also be included in computer readable medium.
Particular embodiments described above, has carried out detailed further to the purpose of the present invention, technical scheme and beneficial effect
Describe in detail bright, be should be understood that the specific embodiment that the foregoing is only the present invention, the guarantor being not intended to limit the present invention
Shield scope, all any modification, equivalent substitution and improvement within the spirit and principles in the present invention, done etc., should be included in this
Within the protection domain of invention.
Claims (9)
1. a kind of method that seismic profile shows is it is characterised in that comprise the following steps:
(s1), for each pixel of screen, its coverage is mapped to road collection space from pixel space, obtains each pixel
The seismic channel being covered and sampled point；
(s2) seismic channel, being covered according to each pixel and sampled point is corresponding generates the moire pattern picture being located in each pixel；
(s3), obtain each moire pattern picture in each pixel brightness contribution value to this pixel, and by each pixel described
Each moire pattern picture the brightness contribution value of this pixel is added up after as this pixel brightness value；
(s4), the brightness value of each pixel is cached in video memory；
(s5), the brightness value of each pixel of caching in described video memory is drawn to screen display；Wherein:
Described each pixel for screen, its coverage is mapped to road collection space from pixel space, specifically includes:
Each pixel is transformed into world coordinates from screen coordinate: w=p-1*m-1* s, wherein p-1For projection matrix inverse matrix, m-1
For model matrix inverse matrix, s is the screen coordinate of pixel, and w is the corresponding world coordinates of pixel；
The brightness contribution value to this pixel for described each moire pattern picture obtaining in each pixel, and by each pixel described
Each moire pattern picture the brightness contribution value of this pixel is added up after as this pixel brightness value, specifically include:
The world coordinates w of each pixel is evenly dividing into some grid g, puts each grid giInitial value be 0；
Calculate the moire pattern picture intersecting with the world coordinates w of each pixel, and by the coverage of these moire pattern pictures with this
Each grid g of pixeliIt is test for intersection, and the grid g that will be cappediValue put 1；
Calculate the number of the grid that its value in each pixel is 1, and this number is set to the brightness value of current pixel
2. the method that seismic profile according to claim 1 shows is it is characterised in that step (s1) to (s4) is by figure
Reason device gpu texture generates thread and completes；Step (s5) is completed by gpu texture rendering thread.
3. seismic profile according to claim 2 shows method it is characterised in that when receiving image scaling instruction,
Seismic profile currently displaying for described screen is carried out corresponding scaling and processes the interim earthquake of generation by described gpu texture rendering thread
Section output display, and trigger described gpu texture rendering thread simultaneously, treat that described gpu texture rendering thread generates earthquake and cuts open
During face, interim seismic profile currently displaying for described screen is substituted for the seismic profile that described gpu texture rendering thread generates.
4. the method that the seismic profile according to Claims 2 or 3 shows is it is characterised in that described screen is solid by several
The texture cell of sizing is covered, and each texture cell described is configured with an independent texture and generates thread, and each is only
Vertical texture generates the seismic profile texture that thread parallel generates this texture cell.
5. the method that seismic profile according to claim 4 shows is it is characterised in that the edge of adjacent texture cell is mutual
Overlapping.
6. the method that seismic profile according to claim 5 shows is it is characterised in that the border of adjacent texture cell is mutual
Overlapping 2 pixels.
7. seismic profile according to claim 6 shows method it is characterised in that described by described video memory in cache
The brightness value of each pixel draw to screen display, adjacent texture cell reduces by a pixel at border overlay.
8. seismic profile according to claim 1 shows method is it is characterised in that world coordinates w by each pixel
It is evenly dividing into 16 × 16 grid g.
9. seismic profile according to claim 4 shows method it is characterised in that each texture cell cover 512 ×
512 or 1024 × 1024 pixels.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410743516.5A CN104407384B (en) | 2014-12-08 | 2014-12-08 | Seismic section display method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410743516.5A CN104407384B (en) | 2014-12-08 | 2014-12-08 | Seismic section display method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN104407384A CN104407384A (en) | 2015-03-11 |
| CN104407384B true CN104407384B (en) | 2017-01-25 |
Family
ID=52645028
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201410743516.5A Active CN104407384B (en) | 2014-12-08 | 2014-12-08 | Seismic section display method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN104407384B (en) |
Families Citing this family (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104751494A (en) * | 2015-04-27 | 2015-07-01 | 国家电网公司 | Processing and drawing method for OpenGL huge textures |
| CN104849755A (en) * | 2015-05-04 | 2015-08-19 | 中国石油天然气集团公司 | Parallel earthquake data rendering method |
| CN108898643B (en) * | 2018-06-15 | 2022-03-11 | 广东数相智能科技有限公司 | Image generation method, device and computer readable storage medium |
| CN113325475A (en) * | 2020-02-28 | 2021-08-31 | 中国石油天然气集团有限公司 | Seismic profile display method and device |
| CN114063155B (en) * | 2020-07-30 | 2024-04-09 | 中国石油化工股份有限公司 | Method, system, storage medium and electronic device for optimizing seismic profile analysis |
Family Cites Families (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO1998012667A2 (en) * | 1996-08-29 | 1998-03-26 | Johnson Steven A | Wavefield imaging using inverse scattering techniques |
| ITMI20060505A1 (en) * | 2006-03-21 | 2007-09-22 | Eni Spa | METRODO TO DISPLAY AND APPEAR IMAGES AND VOLUMES OF DATA OF PHYSICAL SIZES |
| MX2009002237A (en) * | 2006-09-01 | 2009-10-12 | Landmark Graphics Corp | Systems and methods for imaging waveform volumes. |
| CA2674846C (en) * | 2007-01-05 | 2017-01-03 | Landmark Graphics Corporation, A Halliburton Company | Systems and methods for visualizing multiple volumetric data sets in real time |
| CN102750727B (en) * | 2012-07-05 | 2014-11-19 | 武汉大学 | Access memory method for realizing shear wave data three-dimensional visualization by aiming at parallel volume rendering |
| CN102831577B (en) * | 2012-08-29 | 2016-08-24 | 电子科技大学 | The quick Zoom method of two-dimension earthquake image based on GPU |
-
2014
- 2014-12-08 CN CN201410743516.5A patent/CN104407384B/en active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN104407384A (en) | 2015-03-11 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN104407384B (en) | Seismic section display method | |
| US11954759B2 (en) | Tile-based graphics | |
| CN103945208B (en) | A kind of parallel synchronous zooming engine for multiple views bore hole 3D display and method | |
| CN103308942B (en) | A kind of method and system of visual geological data | |
| CN102902848B (en) | Three-dimensional visualization simulation method directed at boundary element analysis | |
| CN107154063A (en) | The shape method to set up and device in image shows region | |
| CN103605506B (en) | The display packing of a kind of graphical user interface resource and terminal | |
| US9565414B2 (en) | Efficient stereo to multiview rendering using interleaved rendering | |
| CN104104888A (en) | Parallel multi-core FPGA digital image real-time zooming processing method and device | |
| KR20200040883A (en) | Multi-spatial rendering with configurable transformation parameters | |
| US9519996B2 (en) | Virtual view generating method and apparatus | |
| KR102442488B1 (en) | Graphics processing systems and graphics processors | |
| US9519982B2 (en) | Rasterisation in graphics processing systems | |
| CN103116897B (en) | A kind of Three-Dimensional Dynamic data compression based on image space and smoothing method | |
| CN104849755A (en) | Parallel earthquake data rendering method | |
| CN108133454A (en) | Model space geometric image switching method, device, system and interactive device | |
| US20100141658A1 (en) | Two-dimensional shadows showing three-dimensional depth | |
| US20170140569A1 (en) | System and method for optimized sparse volume rendering | |
| US8872827B2 (en) | Shadow softening graphics processing unit and method | |
| CN108062339B (en) | Processing method and device of visual chart | |
| GB2247596A (en) | Anti-aliasing system | |
| Shih et al. | Out-of-core visualization of time-varying hybrid-grid volume data | |
| CN104463941A (en) | Volume rendering method and device | |
| CN107146193A (en) | A kind of GPU parallel calculating methods based on double video cards applied to image procossing | |
| KR102589969B1 (en) | Graphics processing unit, graphics processing system and graphics processing method of performing interpolation in deferred shading |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| C14 | Grant of patent or utility model | ||
| GR01 | Patent grant |