CN103377456B - A kind of adaptive digital raster map visible watermark method - Google Patents

A kind of adaptive digital raster map visible watermark method Download PDF

Info

Publication number
CN103377456B
CN103377456B CN201210114437.9A CN201210114437A CN103377456B CN 103377456 B CN103377456 B CN 103377456B CN 201210114437 A CN201210114437 A CN 201210114437A CN 103377456 B CN103377456 B CN 103377456B
Authority
CN
China
Prior art keywords
map
visible watermark
watermark
information
coefficient
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
Application number
CN201210114437.9A
Other languages
Chinese (zh)
Other versions
CN103377456A (en
Inventor
朱长青
符浩军
曹曲
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Nanjing Jiyin Information Technology Co., Ltd.
Original Assignee
Changshu Nanjing Normal University Development Research Institute Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Changshu Nanjing Normal University Development Research Institute Co Ltd filed Critical Changshu Nanjing Normal University Development Research Institute Co Ltd
Priority to CN201210114437.9A priority Critical patent/CN103377456B/en
Publication of CN103377456A publication Critical patent/CN103377456A/en
Application granted granted Critical
Publication of CN103377456B publication Critical patent/CN103377456B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Image Processing (AREA)
  • Editing Of Facsimile Originals (AREA)

Abstract

The present invention relates to a kind of adaptive digital raster map visible watermark method: for digital raster map, there is compared with natural image higher brightness and less saturation, spatial domain pixel value change the features such as relatively larger, based on wavelet transformation and human visual system's characteristic, first visible watermark is carried out pretreatment, then at digital raster map, visible watermark embedding after wavelet transformation embedded location selected in map adaptively, the selection of embedded location is determined according to digital raster map data characteristics and human visual system's characteristic.Visible watermark is semi-transparently presented on former digital raster map by the present invention, maintains the detailed information of original map well, provides better copyright protection for digital raster map.

Description

A kind of adaptive digital raster map visible watermark method
Technical field
The present invention relates to the technical field of geography information copyright protection, particularly to a kind of adaptive digital raster map visible watermark method.
Background technology
Digital raster map is the important component part of numerical map, is also the key data source of geospatial information acquisition.Along with going deep into of NSDI construction, the maturation of GIS-Geographic Information System (GIS) technology and extensive use, digital raster map has important using value in country's key project construction, military surveillance, dynamic monitoring important goal, calamity emergency and process, territory supervision etc..
In recent years, along with developing rapidly of digitized and network times, the storage of digital raster map, transmission, duplication all become very convenient fast.The illegal copies of digital raster map and duplication, not only can make the owner of digital raster map suffer tremendous economic loss, more entail dangers to national security, national defense safety etc..This makes data owner be reluctant open easily or issue its product, is reluctant to share, dare not share, seriously hinders the development of geomatics industry.Therefore, the safety of digital raster map and intellectual property how is protected to become extremely urgent problem.
Digital watermark technology is as a kind of brand-new information security technology, and the safeguard protection for digital raster map provides a kind of practicable solution route.Digital watermarking refers to embedding watermark information in digitalized data, watermark information and source data is combined together, becomes the inseparable part of source data.Thus determine copyright owner, proprietary rights certification, tracking abuse, the verity in authenticated digital content source, identification buyer, other additional informations etc. about digital content are provided.Simultaneously for strengthening responsibility, fright illegal act, having good grounds to look into soon etc. has important function.The effect that digital watermark technology plays in military security guarantee, national security has been subjected to the great attention of national government office.
Digital watermarking achieves many application in the safeguard protection in the fields such as image, figure, video, audio frequency, have also been obtained important application in mapping association area in recent years.
Invisible watermark technology about digital raster map has much research, and these researchs have carried out more deep analysis with regard to problems such as the copyright protection of digital raster map, content authentication, tracking uses and proposed corresponding solution;Comparatively speaking, as a kind of more positive copyright protection technology, it is seen that digital watermark also has some to study in digital raster map, its research of present stage and application focus primarily upon on natural image.GUOJie (Visiblewatermarkingtechniquebasedonluminanceandtexturefe atures [J] .InfraredandLaserEngineering, 2003,32 (1): 92-100.) visible watermarking arithmetic of a kind of DCT domain is devised, utilize the local luminance of image and textural characteristics to determine that watermark embeds the factor, so under the premise meeting visible watermark characteristic, strengthen watermark embedment strength, improve watermark robustness;MKankanhalli (Adaptivevisiblewatermarkingofimages [A] .Proc.ofICMCS [C] .Florence, Italy, IEEE, 1999.) each DCT block obtains watermark drawing coefficient to propose to utilize the sensitivity of image content features and brightness to determine, SPMohanty (ADCTdomainvisiblewatermarkingtechniqueforimages [A] .Proc.Int.Conf.onMultimediaandExpo [C] .USA:IEEE, 2000.) this algorithm has been improved, enhance the adaptivity of watermark processing.But the common problem that DCT domain visible watermark processes is that piecemeal stretching watermarking images may result in visual discontinuity;GONGSheng-rong (AVisibleWatermarkingAlgorithmHoldingImageContent [J] .GeomaticsandInformaticaScienceofWuhanUniversity, 2006,31 (9): 757-764.) a kind of wavelet field visible watermarking arithmetic is devised according to the masking by noise of image, and the textural characteristics according to watermarking images improves algorithm, original image content can be kept preferably;HUANGBiao-bing (VisibleWatermarkingSchemeWithHumanVisualSystem [J] .ComputerEngineeringAndApplications, 2006,13:63-65.) according to the local of master image and watermarking images and global characteristics, set up drawing coefficient mathematical model, proposing a kind of visible watermarking arithmetic with picture material retention performance, this algorithm experimental result is better;ZHAOYou-jun (AdaptiveVisibleWatermarkTechniqueinRegionsofInterest [J] .ComputerEngineering, 2007,33 (9): 180-181.) propose the self adaptation visible watermarking technique of a kind of area-of-interest based on image particular content, but its visible watermark embeds region and must user select.
These researchs achieve some achievements, but in current visible watermarking technique research and application, are primarily directed to natural image, and to digital raster map, this aspect such as product own characteristic, visible watermarking arithmetic robustness of manually charting considers also deficiency.The present invention, on the basis that the visible watermarking arithmetic of current natural image is analyzed, in conjunction with digital raster map data characteristic, provides better copyright protection for digital raster map.
Summary of the invention
It is an object of the invention to: compare natural image for digital raster map and have the advantages that higher brightness and less saturation, the change of spatial domain pixel value are relatively larger, watermark information is embedded in an adaptive way digital raster map, so under the premise meeting visible watermark characteristic, farthest embed watermark, improve watermark robustness.Having carried out watermark period expansion before embedding simultaneously, and embedded it in digital raster map vision important area, this makes under the premise not damaging map use quality, it is not easy to remove the visible watermark in map, thus protecting the copyright of map preferably.
In order to realize foregoing invention purpose, the technical solution used in the present invention is: a kind of adaptive digital raster map visible watermark method, there is higher brightness and less saturation for digital raster map compared with natural image, the features such as spatial domain pixel value change is relatively larger, based on wavelet transformation and human visual system's characteristic, first visible watermark is carried out pretreatment, then at digital raster map, visible watermark is embedded adaptively the embedded location selected in map after wavelet transformation, the selection of embedded location is determined according to digital raster map data characteristics and human visual system's characteristic.
Technique scheme farther includes three below step:
Step one: watermark information pretreatment;
Step 2: the selection of visible watermark embedded location;
Step 3: visible watermark is embedded digital raster map based on wavelet transformation.
Described step one, is that visible watermark W is expanded to visible watermark is W ', to ensure that it can cover most of region of original place figure;
Described step 2, is on the basis of step one, according to a definite sequence, whole numerical map is traveled through, and calculates the importance rate index I of chart portfolio block, at all of I(i,j)Middle maximizing, the region corresponding to this maximum is exactly the visible watermark embedded location found;
Described step 3, it is on the basis of step 2, utilize low (height) after original place figure and visible watermark wavelet transformation frequency coefficient, linear adder is adopted to embed criterion, obtain low (height) frequency coefficient after visible watermark embeds original place figure, subsequently by wavelet inverse transformation, obtain the digital raster map containing visible watermark.
The present invention, based on wavelet transformation and human visual system's characteristic, farthest embeds watermark, improves watermark robustness;Make under the premise not damaging map use quality, it is not easy to remove the visible watermark in map, thus protecting the copyright of map better simultaneously.
Accompanying drawing explanation
The watermark that Fig. 1 is the inventive method embeds flow chart.
Fig. 2 is a width original figure grating map of the prior art.
Fig. 3 is visible watermark image.
Fig. 4 is ω1(x, value y).
Fig. 5 is the design sketch after Fig. 2 embeds Fig. 3.
Detailed description of the invention
Below in conjunction with Figure of description, the invention is described in further detail.
Embodiment one
As shown in Figures 2 to 5, select in the present embodiment a width be sized to 808 × 758 256 grades of digital raster maps and a width to be sized to 256 grades of images of 200 × 200 be experimental data, the adaptive feed-forward network of visible watermark is carried out according to the present invention, provide one embodiment of the present of invention, further describe the present invention.
Step one: watermark information pretreatment;
Digital raster map size according to the present embodiment and watermark size, it is determined that watermark extended amplitude E=min{808/200,758/200}, after extension, visible watermark is W ', is sized to m ' × n ', wherein m '=E × m, n '=E × n.
Step 2: the selection of visible watermark embedded location;
On the basis of step one, with m ' × n ' for template size, poor with h for template movement, according to a definite sequence, whole numerical map is traveled through, for finding out digital raster map vision important area.If the region that when traversal starts, template covers is Reg (1,1), then set template successively in the region that traversal covers on the way as Reg (i, j), wherein i=1,2 ..., [(M-m ')/h], j=1,2 ..., [(N-n ')/h], [] represents that logarithm rounds.It is utilized respectively following formula and obtains the importance rate index I of these map blocks.
I ( i , j ) = Σ x = i i + m ′ Σ y = j j + n ′ ω 1 ( x , y ) · ω 2 ( x , y ) · F ( x , y )
Wherein I(i,j)Represent that (i, j) the importance rate index of map block represent that in the figure of original place, (x, y) the characters of ground object value at place, mainly by (x, y) place's brightness value determines, ω to Reg in traversal1(x,y)、ω2(x y) represents (x, y) weight at place, wherein ω in the figure of original place1(x, y) by (x, y) the liftoff figure centre distance in place determines, due to the optic centre of human eye, to be positioned at observed middle parts of objects slightly elevated, and therefore the visual importance in map center region generally should be higher than that other regions.Here map is divided into nine regions, ω1(x, value y) is as shown in Figure 4.ω2(x, y) by (x, y) whether place is that atural object marginal information determines, in general, in map, its importance of marginal information should be higher than that non-edge information.Test shows, works as ω2(x, during the such as following formula of value condition y), algorithm effect is more satisfactory.
After so digital raster map being traveled through, calculate the importance rate index of chart portfolio block, at all of I(i,j)Middle maximizing, the region corresponding to this maximum is exactly the visible watermark embedded location found.
Step 3: visible watermark is embedded digital raster map based on wavelet transformation.
On the basis of step 2, original place figure and visible watermark carrying out low (height) after wavelet transformation frequency coefficient, adopts linear adder to embed criterion, obtain low (height) frequency coefficient after visible watermark embeds original place figure, its formula can be described as:
C′ij(n)=α (i, j) Ci,j(n)+β(i,j)·Wi,j(n)
α in formula (i, j) and β (i j) is the embedding factor of the scale factor of original place figure and visible watermark, C' respectivelyijSmall echo low (height) the frequency coefficient of map, C after (n) expression visible watermark embeddingi,j(n) and Wi,jN () represents original place figure and visible watermark low (height) frequency coefficient after wavelet transformation respectively, α (i, j) with β (i, j) determined by digital raster map data characteristics and wavelet coefficient, and the drawing coefficient of low frequency sub-band and high-frequency sub-band determines that rule is different.
In low frequency sub-band, its drawing coefficient is mainly determined by monochrome information.Different brightness are had different sensitivity by human vision, in general, intermediate light is most sensitive, have the characteristic of non-linear decline to low-light level and high brightness direction, and thus the present invention adopts following formula definition low frequency sub-band drawing coefficient:
α ( i , j ) = w i j · ( D i j - D a v e D max ) 2
β (i, j)=1-α (i, j)
D in formulamaxRepresent the coefficient maximum of low frequency sub-band, DaveRepresent the coefficient meansigma methods of low frequency sub-band, DijRepresent the coefficient of low frequency sub-band, wijRepresenting the quantization weight factor, its effect is can to keep original place figure vision validity as far as possible, makes watermark have good observability simultaneously.Consider that in map, pixel value change in spatial domain is relatively larger, atural object marginal information is relatively more, and atural object marginal area is most important visual information and geography information in map, and human eye is very sensitive to the change at atural object edge, in order to keep the complete of map, it is necessary for keeping the complete of atural object edge.So determining scale factor α (i, j) with embedding factor-beta (i, j) time, atural object marginal area answers special consideration, atural object marginal area is changed minimum, therefore can only visible watermark information low-intensity be embedded in the edge block of original place figure, that is can define scale factor with above formula and embed the factor in map non-edge block, map edge block then needs adjust accordingly, scale factor α (i herein, j) can suitably strengthen, and (i j) should suitably reduce to embed factor-beta.And in map background region, not including meaning cartographic information, human eye is also insensitive to it, therefore (i, j), improves and embeds factor-beta (i can suitably to reduce scale factor α, j), ensureing on the basis of visible watermark visual consistency, prominent visible watermark information.So, by adjusting scale factor and embedding the factor, original place figure visual quality can be kept preferably, also make the observability that watermark has had simultaneously.
The high-frequency sub-band of digital raster map wavelet decomposition reflects the minutia of map.Its approximate Laplace of obedience is distributed, and major part coefficient is close to 0, and only fraction has bigger peak value corresponding to the coefficient of map edge and texture region.High frequency coefficient distinguishes the smooth region on correspondence spatial domain, texture region and marginal area from small to large, and coefficient non-zero number and situation of change represent the texture information of this block, and nonzero coefficient is more many, and index variation is more big, and this block texture is more complicated.In general, flat site at original place figure answers performance visible watermark texture information as much as possible, namely at map flat site, if there is abundant texture information in corresponding visible watermark region, then based on visible watermark texture information, so can highlight visible watermark information, and if visible watermark region is also flat site, then with map area information for main display.Simultaneously, human visual system is insensitive to the change of high texture region ground pixel, and the texture number in the same area is more many, and the sheltering of noise is more big, vision distortion is more little, is necessary for strengthening visible watermark high frequency coefficient for display visible watermark detailed information and embeds the factor.Therefore, in a map area being rich in texture information, information is uniformly distributed in each high frequency coefficient, and this means that the variances sigma in this region(i,j)Less, need to add more watermark information for display.In the region atlas dimension factor-alpha (i, j) and visible watermark embed factor-beta (i, j) can be defined as:
α (i, j)=σ '(i,j)·w(i,j)
β (i, j)=1-α (i, j)
σ ' (i, j)=σ in formula(i,j)max, σ(i,j)Represent the high frequency variance of map block, σmaxRepresent map block high frequency variance maximum, w(i,j)Represent weight factor.
Draw low (height) frequency coefficient after visible watermark embedding original place figure by the above-mentioned steps of the present invention, subsequently by wavelet inverse transformation, obtain the digital raster map containing visible watermark, such as Fig. 5.
5. the feature of this method and technical advantage:
(1) this method is compared natural image for digital raster map and is had the advantages that higher brightness and less saturation, the change of spatial domain pixel value are relatively larger, watermark information is embedded in an adaptive way digital raster map, so under the premise meeting visible watermark characteristic, farthest embed watermark, improve watermark robustness.
(2) this method has carried out watermark period expansion before embedding; and embed it in digital raster map vision important area; this makes under the premise not damaging map use quality, it is not easy to remove the visible watermark in map, thus protecting the copyright of map preferably.
6. the test of the present invention and analysis:
Method proposed by the invention is specific to the water mark method of digital raster map, adopts the method can develop and realize the digital raster copyright protection system based on digital watermark technology.
(1) visual effect
As it is shown in figure 5, visible watermark can more visible being presented in map, and cover most of region of map, also more visible can reflect former cartographic information containing visible watermark map.
(2) Map quality aspect
As shown in Figure 5, the method can keep former digital raster map quality preferably, is presented in map very important visual region with making again visible watermark translucence, and covers map major part region, especially at texture region and marginal area, the fabulous details maintaining original place figure.
(3) algorithm robustness.
Visible watermark covers map major part region, so makes it that map trimming operation is had good robustness;Visible watermark is embedded in digital raster map by algorithm so that cartographic information and visible watermark information organically combine, thus to various map operations such as JPEG compression, add make an uproar, filter, resistivity that sharpening etc. has had;For various geometric attacks such as translation, convergent-divergent, stretching, deletion etc., owing to visible watermark is embedded into map major part region, and organically combining with original place figure, algorithm is effective against geometric attack;Visible watermark adopts adaptive feed-forward network mode, it is determined by this pixel data characteristic at the embedment strength of each pixel, namely the watermark embedment strength of each position is not necessarily identical, therefore the map containing visible watermark carries out reject attack and is difficult to remove visible watermark mark.
Embodiment two
As shown in Figure 1, it is an object of the invention to: compare natural image for digital raster map and have the advantages that higher brightness and less saturation, the change of spatial domain pixel value are relatively larger, watermark information is embedded in an adaptive way digital raster map, so under the premise meeting visible watermark characteristic, farthest embed watermark, improve watermark robustness.Having carried out watermark period expansion before embedding simultaneously, and embedded it in digital raster map vision important area, this makes under the premise not damaging map use quality, it is not easy to remove the visible watermark in map, thus protecting the copyright of map preferably.
In order to realize foregoing invention purpose, the technical solution used in the present invention is: a kind of adaptive digital raster map visible watermark method, there is higher brightness and less saturation for digital raster map compared with natural image, the features such as spatial domain pixel value change is relatively larger, based on wavelet transformation and human visual system's characteristic, first visible watermark is carried out pretreatment, then at digital raster map, visible watermark is embedded adaptively the embedded location selected in map after wavelet transformation, the selection of embedded location is determined according to digital raster map data characteristics and human visual system's characteristic.
Technique scheme farther includes three below part:
Step one: watermark information pretreatment;
Step 2: the selection of visible watermark embedded location;
Step 3: visible watermark is embedded digital raster map based on wavelet transformation.
Described step one, is that visible watermark W is expanded to visible watermark is W ', to ensure that it can cover most of region of original place figure;
Described step 2, is on the basis of step one, according to a definite sequence, whole numerical map is traveled through, and calculates the importance rate index I of chart portfolio block, at all of I(i,j)Middle maximizing, the region corresponding to this maximum is exactly the visible watermark embedded location found;
Described step 3, it is on the basis of step 2, utilize low (height) after original place figure and visible watermark wavelet transformation frequency coefficient, linear adder is adopted to embed criterion, obtain low (height) frequency coefficient after visible watermark embeds original place figure, subsequently by wavelet inverse transformation, obtain the digital raster map containing visible watermark.
The present invention, based on wavelet transformation and human visual system's characteristic, farthest embeds watermark, improves watermark robustness;Make under the premise not damaging map use quality, it is not easy to remove the visible watermark in map, thus protecting the copyright of map better simultaneously.

Claims (1)

1. an adaptive digital raster map visible watermark method, it is characterized in that: for digital raster map, there is compared with natural image higher brightness and less saturation, spatial domain pixel value change relatively larger feature, based on wavelet transformation and human visual system's characteristic, first visible watermark is carried out pretreatment, then at digital raster map, visible watermark embedding after wavelet transformation embedded location selected in map adaptively, the selection of embedded location is determined according to digital raster map data characteristics and human visual system's characteristic;Described method comprises the steps:
Step one: watermark information pretreatment;
According to digital raster map size and watermark size, it is determined that watermark extended amplitude E=min{808/200,758/200}, after extension, visible watermark is W ', is sized to m ' × n ', wherein m '=E × m, n '=E × n;M and n represents the size of the front original watermark of extension;
Step 2: the selection of visible watermark embedded location;
On the basis of step one, with m ' × n ' for template size, poor with h for template movement, according to a definite sequence, whole numerical map is traveled through, for finding out digital raster map vision important area;If the region that when traversal starts, template covers is Reg (1,1), then set template successively in the region that traversal covers on the way as Reg (i, j), wherein i=1,2, ..., [(M-m ')/h], j=1,2, ..., [(N-n ')/h], [] represents that logarithm rounds, and is utilized respectively following formula and obtains the importance rate index I of these map blocks;M and N represents original figure grating map size;
I ( i , j ) = Σ x = i i + m ′ Σ y = j j + n ′ ω 1 ( x , y ) · ω 2 ( x , y ) · F ( x , y )
Wherein, (x y) represents digital raster map (x, y) the attribute correlation at place to F;I(i,j)Represent that (i, j) the importance rate index of map block represent that in the figure of original place, (x, y) the characters of ground object value at place, by (x, y) place's brightness value determines, ω to Reg in traversal1(x,y)、ω2(x y) represents (x, y) weight at place, wherein ω in the figure of original place1(x, y) by (x, y) the liftoff figure centre distance in place determines, due to the optic centre of human eye, to be positioned at observed middle parts of objects slightly elevated, and therefore the visual importance in map center region is higher than other regions, and map is divided into nine regions here;ω2(x, y) by (x, y) whether place is that atural object marginal information determines, in map, its importance of marginal information should be higher than that non-edge information;Test shows, works as ω2(x, during the such as following formula of value condition y), algorithm effect is more satisfactory:
After so digital raster map being traveled through, calculate the importance rate index of chart portfolio block, at all of I(i,j)Middle maximizing, the region corresponding to this maximum is exactly the visible watermark embedded location found;
Step 3: visible watermark is embedded digital raster map based on wavelet transformation;
On the basis of step 2, original place figure and visible watermark carrying out the low high coefficient after wavelet transformation, adopt linear adder to embed criterion, obtain the low high coefficient after visible watermark embeds original place figure, its formula can be described as:
C′ij(n)=α (i, j) Ci,j(n)+β(i,j)·Wi,j(n)
α in formula (i, j) and β (i j) is the embedding factor of the scale factor of original place figure and visible watermark, C ' respectivelyijThe small echo low high coefficient of map, C after (n) expression visible watermark embeddingi,j(n) and Wi,jN () represents original place figure and the visible watermark low high coefficient after wavelet transformation respectively, α (i, j) (i, j) is determined by digital raster map data characteristics and wavelet coefficient, and the drawing coefficient of low frequency sub-band and high-frequency sub-band determines that rule is different with β;
In low frequency sub-band, its drawing coefficient is mainly determined by monochrome information;Different brightness are had different sensitivity by human vision, and intermediate light is most sensitive, have the characteristic of non-linear decline to low-light level and high brightness direction, thus adopt following formula definition low frequency sub-band drawing coefficient:
α ( i , j ) = w i j · ( D i j - D a v e D max ) 2
β (i, j)=1-α (i, j)
D in formulamaxRepresent the coefficient maximum of low frequency sub-band, DaveRepresent the coefficient meansigma methods of low frequency sub-band, DijRepresent the coefficient of low frequency sub-band, wijRepresenting the quantization weight factor, its effect is can to keep original place figure vision validity as far as possible, makes watermark have good observability simultaneously;Consider that in map, pixel value change in spatial domain is relatively larger, atural object marginal information is relatively more, and atural object marginal area is most important visual information and geography information in map, and human eye is very sensitive to the change at atural object edge, in order to keep the complete of map, it is necessary for keeping the complete of atural object edge;So determining scale factor α (i, j) with embedding factor-beta (i, j) time, atural object marginal area answers special consideration, atural object marginal area is changed minimum, therefore can only visible watermark information low-intensity be embedded in the edge block of original place figure, that is can define scale factor with above formula and embed the factor in map non-edge block, map edge block then needs adjust accordingly, scale factor α (i herein, j) can suitably strengthen, and embed factor-beta (i, j) should suitably reduce, and in map background region, do not include meaning cartographic information, human eye is also insensitive to it, therefore scale factor α (i can suitably be reduced, j), improve and embed factor-beta (i, j), ensureing on the basis of visible watermark visual consistency, prominent visible watermark information;So, by adjusting scale factor and embedding the factor, original place figure visual quality can be kept preferably, also make the observability that watermark has had simultaneously;
The high-frequency sub-band of digital raster map wavelet decomposition reflects the minutia of map, and its approximate Laplace of obedience is distributed, and major part coefficient is close to 0, and only fraction has bigger peak value corresponding to the coefficient of map edge and texture region;High frequency coefficient distinguishes the smooth region on correspondence spatial domain, texture region and marginal area from small to large, and coefficient non-zero number and situation of change represent the texture information of this block, and nonzero coefficient is more many, and index variation is more big, and this block texture is more complicated;Flat site at original place figure answers performance visible watermark texture information as much as possible, namely at map flat site, if there is abundant texture information in corresponding visible watermark region, then based on visible watermark texture information, so can highlight visible watermark information, if and visible watermark region is also flat site, then with map area information for main display;Simultaneously, human visual system is insensitive to the change of high texture region ground pixel, and the texture number in the same area is more many, and the sheltering of noise is more big, vision distortion is more little, is necessary for strengthening visible watermark high frequency coefficient for display visible watermark detailed information and embeds the factor;Therefore, in a map area being rich in texture information, information is uniformly distributed in each high frequency coefficient, and this means that the variances sigma in this region(i,j)Less, for display need can add more watermark information, in the region atlas dimension factor-alpha (i, j) and visible watermark embed factor-beta (i, j) is defined as:
α (i, j)=σ '(i,j)·w(i,j)
β (i, j)=1-α (i, j)
σ ' (i, j)=σ in formula(i,j)max, σ(i,j)Represent the high frequency variance of map block, σmaxRepresent map block high frequency variance maximum, w(i,j)Represent weight factor;
Draw the low high coefficient after visible watermark embedding original place figure by above-mentioned steps, subsequently by wavelet inverse transformation, obtain the digital raster map containing visible watermark.
CN201210114437.9A 2012-04-19 2012-04-19 A kind of adaptive digital raster map visible watermark method Active CN103377456B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201210114437.9A CN103377456B (en) 2012-04-19 2012-04-19 A kind of adaptive digital raster map visible watermark method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201210114437.9A CN103377456B (en) 2012-04-19 2012-04-19 A kind of adaptive digital raster map visible watermark method

Publications (2)

Publication Number Publication Date
CN103377456A CN103377456A (en) 2013-10-30
CN103377456B true CN103377456B (en) 2016-06-29

Family

ID=49462526

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201210114437.9A Active CN103377456B (en) 2012-04-19 2012-04-19 A kind of adaptive digital raster map visible watermark method

Country Status (1)

Country Link
CN (1) CN103377456B (en)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US10163254B2 (en) * 2016-06-24 2018-12-25 Adobe Systems Incorporated Rendering of digital images on a substrate
CN107274330B (en) * 2017-05-24 2019-12-17 中国地质大学(武汉) BFA and LSB based grid map copyright protection method
CN110415153A (en) * 2018-04-27 2019-11-05 北京大学 A kind of watermark embedding method, system and a kind of watermark extracting method and system
CN113450243A (en) * 2020-03-24 2021-09-28 北京四维图新科技股份有限公司 Watermark adding method and device
CN115564632B (en) * 2022-10-20 2024-04-26 中科星图数字地球合肥有限公司 Watermark embedding method and device for grid map and electronic equipment

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
基于小波分析的栅格数字地图水印技术研究;符浩军;《中国优秀硕士学位论文全文数据库(电子期刊)信息科技辑》;20120215(第2012年第02期);正文第6章 *
基于小波变换的栅格数字地图水印算法;符浩军 等;《测绘学报》;20090531;第34卷(第3期);91,107,108 *
基于整数小波变换的栅格数字地图数字水印算法;朱长青 等;《武汉大学学报·信息科学版》;20090531;第34卷(第5期);619-625 *

Also Published As

Publication number Publication date
CN103377456A (en) 2013-10-30

Similar Documents

Publication Publication Date Title
CN103377456B (en) A kind of adaptive digital raster map visible watermark method
CN101493928B (en) Digital watermarking embedding, extracting and quantizing step size coordinating factor optimizing method and device
Khan et al. Digital watermarking for image authenticationbased on combined dct, dwt and svd transformation
CN103295204B (en) A kind of image self-adapting enhancement method based on non-down sampling contourlet transform
Li et al. Watermarking 3D mesh by spherical parameterization
CN104103082A (en) Image saliency detection method based on region description and priori knowledge
CN104331450B (en) Video copying detection method based on multi-mode feature and tensor resolution
CN112528969B (en) Face image authenticity detection method and system, computer equipment and storage medium
CN112801846B (en) Watermark embedding and extracting method and device, computer equipment and storage medium
Zhang et al. An image watermarking method based on visual saliency and contourlet transform
US10169908B2 (en) Method, apparatus, storage medium and device for controlled synthesis of inhomogeneous textures
CN106408497A (en) Strong-robustness watermark embedding and extraction method for original remote sensing images
Soppari et al. Development of improved whale optimization-based FCM clustering for image watermarking
Ramakrishnan et al. Svd based robust digital watermarking for still images using wavelet transform
Agarwal et al. Visible watermarking based on importance and just noticeable distortion of image regions
Wan et al. JND-guided perceptually color image watermarking in spatial domain
Jiang et al. Adaptive spread transform QIM watermarking algorithm based on improved perceptual models
Wang et al. A novel attention-guided JND Model for improving robust image watermarking
Wang et al. Multiscale single image dehazing based on adaptive wavelet fusion
CN105023237A (en) Method for improving concealment performance of image digital watermarks
Hoshi et al. A robust watermark algorithm for copyright protection by using 5-level DWT and two logos
Abodena et al. Hybrid technique for robust image watermarking using discrete time fourier transform
Liu et al. Novel robust multiple watermarking against regional attacks of digital images
Kamble et al. DWT-based multiple watermarking for privacy and security of digital images in e-commerce
Jadav Comparison of LSB and Subband DCT Technique for Image Watermarking

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
TR01 Transfer of patent right

Effective date of registration: 20170410

Address after: Qixia District of Nanjing City, Jiangsu province 210034 yuan Xianlin University City Road No. 8

Patentee after: Nanjing Jiyin Information Technology Co., Ltd.

Address before: 215500 Suzhou, Changshu, Southeast of Changshou City Economic Development Zone, Jin Road, No. 8, No.

Patentee before: Changshu Nanjing Normal University Development Research Academy Institute Co., Ltd.

TR01 Transfer of patent right