JP4175551B2 - Electronic information embedding method and extraction method, electronic information embedding device and extraction device, and program thereof - Google Patents
Electronic information embedding method and extraction method, electronic information embedding device and extraction device, and program thereof Download PDFInfo
- Publication number
- JP4175551B2 JP4175551B2 JP2007040094A JP2007040094A JP4175551B2 JP 4175551 B2 JP4175551 B2 JP 4175551B2 JP 2007040094 A JP2007040094 A JP 2007040094A JP 2007040094 A JP2007040094 A JP 2007040094A JP 4175551 B2 JP4175551 B2 JP 4175551B2
- Authority
- JP
- Japan
- Prior art keywords
- watermarked
- point cloud
- cloud data
- data
- electronic information
- 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.)
- Expired - Fee Related
Links
Description
本発明は、３次元計測に基づいて得られたオリジナルランダム点データ群（３次元点群データ）を情報管理するための方法及びその装置に関し、特に、レーザ３次元計測によって地球表面を計測して得られた３次元点群データの不正使用を防止するために、電子情報を電子透かしデータとして埋め込む際に用いられる方法及び装置に関する。 The present invention relates to a method and apparatus for information management of an original random point data group (three-dimensional point group data) obtained based on three-dimensional measurement, and in particular, by measuring the earth surface by laser three-dimensional measurement. The present invention relates to a method and apparatus used when embedding electronic information as electronic watermark data in order to prevent unauthorized use of the obtained three-dimensional point cloud data.
一般に、偽造又は不正使用を防止するため、印刷物上に多数の微細要素からなるパターンを埋め込んで、そのパターンで形成される情報を利用して、不正使用の検出を行うことが行われている。さらに、ベクトル形式（ベクター型）で記述された地図データに対する不正使用を防止するため、地図データに不正利用防止用の情報（埋め込み情報）を埋め込むことが行われており、ベクター型で記述された２次元又は３次元の地図データにおいては、面を構成するポリゴンに対して埋め込み情報を埋め込んでいる。例えば、地図データを三角形ポリゴンの集合で構成して、これら三角形ポリゴンをさらに４つの三角形に分割し、その間に構成される三角形（どの三角形ポリゴンの頂点も含まない三角形）に電子透かしデータを埋め込みデータとして埋め込む方法が知られている。そして、この方法を用いると、ベクター型で記述された地図データに影響を与えることなく、電子透かしデータを取り除くことが困難となる。 In general, in order to prevent forgery or unauthorized use, a pattern made up of a large number of fine elements is embedded in a printed material, and unauthorized use is detected using information formed by the pattern. Furthermore, in order to prevent unauthorized use of map data described in a vector format (vector type), information for preventing unauthorized use (embedded information) is embedded in the map data. In the two-dimensional or three-dimensional map data, embedding information is embedded in the polygons constituting the surface. For example, map data is composed of a set of triangle polygons, these triangle polygons are further divided into four triangles, and digital watermark data is embedded in the triangles formed between them (the triangles that do not contain any triangle polygon vertices) The method of embedding is known. If this method is used, it becomes difficult to remove the digital watermark data without affecting the map data described in the vector type.
一方、特開２００１−１６０８９７公報には、ベクター型で記述された地図データに不正利用防止用の埋め込み情報を埋め込むことが記載されている。ここでは、オブジェクトの構成点の座標列を示す地図図形情報及びオブジェクトの種別を管理するレイヤ情報に、埋め込み基準レイヤと埋め込み参照レイヤを設定したものを入力して、同一の意味のある領域に存在するオブジェクトペアであって、それらが構成する領域内に他のオブジェクトが無いものを埋め込み基準オブジェクトペアとして選別する。そして、埋め込み基準オブジェクトペアに対して、埋め込み箇所を発見されにくいものを選択して、その存在位置及び／又は形状特徴に基づいて新たに埋め込むオブジェクト情報を算出し、このオブジェクト情報に合わせて埋め込み密度に応じて新たなオブジェクトを埋め込んでいる。
ところで、近年、所謂レーザ３次元計測によって地球表面（地表面）を計測して、地表面のデータをオリジナルランダム点データ群として得て、これらオリジナルランダム点群データから地図データを得ることが行われている。例えば、航空機から地上に向けてパルスレーザビームを照射して、地表面の空間座標を得る。この際、航空機の空間位置は、地上に設けられたＧＰＳ基準局と航空機に搭載されたＧＰＳ受信機によって求め、航空機の姿勢は三軸ジャイロスコープによって求める。
なお、１パルス毎の地表座標は、前述のようにして求めた航空機の空間位置及び姿勢に応じて、レーザミラー射角度及びその斜距離によってレーザビーム反射点のｘ，ｙ，及びｚ座標として求められることになる。
On the other hand, Japanese Patent Laid-Open No. 2001-160897 describes embedding embedded information for preventing unauthorized use in map data described in a vector type. Here, the map graphic information indicating the coordinate sequence of the constituent points of the object and the layer information for managing the type of the object are input with the embedded reference layer and the embedded reference layer set, and exist in the same meaningful area. The object pairs that do not have other objects in the area constituted by them are selected as the embedding reference object pairs. Then, for the embedding reference object pair, an object where it is difficult to find an embedding location is selected, and object information to be newly embedded is calculated based on the existence position and / or shape feature, and the embedding density is matched to this object information. New objects are embedded depending on the situation.
By the way, in recent years, the surface of the earth (the ground surface) is measured by so-called laser three-dimensional measurement, the ground surface data is obtained as an original random point data group, and map data is obtained from the original random point group data. ing. For example, a pulse laser beam is irradiated from the aircraft toward the ground to obtain spatial coordinates of the ground surface. At this time, the spatial position of the aircraft is obtained by a GPS reference station provided on the ground and a GPS receiver mounted on the aircraft, and the attitude of the aircraft is obtained by a three-axis gyroscope.
The ground coordinates for each pulse are obtained as the x, y, and z coordinates of the laser beam reflection point based on the laser mirror firing angle and the oblique distance according to the spatial position and attitude of the aircraft obtained as described above. Will be.
上述のようにして得られた地表座標は単に地表面のランダム点群データであるので、このランダム点群データを処理して地図データを生成するようにしている。
上述のようなオリジナル点群データは、単に空間的に離散している点群データであるので、相互に関係を持っておらず、しかも属性も備えていない。つまり、オリジナル点群データは単にｘ，ｙ，及びｚの座標値を示しているだけである。従って、上述のようなベクター型地図データに対して電子透かしデータを埋め込む手法を用いることができず、不正にオリジナルランダム点群データを用いて地図データを生成されてしまうことになる。
しかも、精度誤差の範囲で、オリジナルランダム点群データに対して乱数に応じて電子透かしデータを埋め込んだ場合には、オリジナルランダム点群データを再現することが困難となってしまう。また、オリジナルランダム点群データを部分的に間引いて、電子透かしデータを埋め込んだ場合には、不正使用を防止することが難しいという課題がある。
Since the surface coordinates obtained as described above are simply random point cloud data on the ground surface, the random point cloud data is processed to generate map data.
Since the original point cloud data as described above is merely point cloud data that is spatially discrete, it does not have a mutual relationship and has no attributes. That is, the original point cloud data merely indicates the coordinate values of x, y, and z. Therefore, the method of embedding digital watermark data in the vector type map data as described above cannot be used, and map data is illegally generated using original random point cloud data.
In addition, when the electronic watermark data is embedded in the original random point group data in accordance with the random number within the accuracy error range, it becomes difficult to reproduce the original random point group data. Further, when the original random point cloud data is partially thinned out and the digital watermark data is embedded, there is a problem that it is difficult to prevent unauthorized use.
本発明の目的は、オリジナルランダム点群データの不正使用を防止することのできる電情報埋め込み方法、電子情報の埋め込みに用いられる装置、及びプログラムを提供することにある。
本発明の他の目的は、電子情報が埋め込まれたオリジナルランダム点群データから電子情報を抽出するための電子情報抽出方法、電子情報の抽出に用いられる装置、及びプログラムを提供することにある。
An object of the present invention is to provide an electronic information embedding method, an apparatus used for embedding electronic information, and a program capable of preventing unauthorized use of original random point cloud data.
Another object of the present invention is to provide an electronic information extraction method for extracting electronic information from original random point cloud data in which electronic information is embedded, an apparatus used for extracting electronic information, and a program.
本発明によれば、複数の頂点座標（ｘ，ｙ，…）で表され、ベクトル型に記述される不規則なデータ集合であるオリジナルランダム点群データに電子情報を電子透かしデータとして埋め込む際に用いられる電子情報埋め込み方法であって、複数の頂点座標の２つが特に規則的配置を与えられていない前記オリジナルランダム点群データについて、前記２つの頂点座標の一方を実数部に他方を虚数部にして１次元的に順序付けて１次元複素離散フーリエ変換を行ってフーリエ係数列を求める第１のステップと、前記電子透かしデータに応じて前記フーリエ係数列を変更して透かし入りフーリエ係数列とする第２のステップと、前記透かし入りフーリエ係数列を逆離散フーリエ変換して該逆離散フーリエ変換に応じて透かし入り点群データを得る第３のステップとを有することを特徴とする電子情報埋め込み方法が得られる。
好ましくは、前記複数の頂点座標が、３次元計測に応じて得られた頂点座標（ｘ，ｙ，ｚ）で表されることを特徴とする。
更に好ましくは、前記１次元複素離散フーリエ変換が１次元複素フーリエ変換を実施した際の離散フーリエ係数Ｃ ｋ の実数部がｘ ｋ 、虚数部がｙ ｋ であるＣ ｋ ＝ｘ ｋ ＋ｉｙ ｋ （ｋは非負整数、ｉは虚数単位）なる複素数列に関する１次元複素離散フーリエ変換を行うことを特徴とする。
According to the present invention, when electronic information is embedded as electronic watermark data in original random point cloud data that is an irregular data set represented by a plurality of vertex coordinates (x, y,...) And described in a vector type. An electronic information embedding method used, wherein two of the plurality of vertex coordinates are not particularly given a regular arrangement, and for the original random point cloud data, one of the two vertex coordinates is a real part and the other is an imaginary part. A first step of performing one-dimensional complex discrete Fourier transform to obtain a Fourier coefficient sequence, and changing the Fourier coefficient sequence according to the digital watermark data to form a watermarked Fourier coefficient sequence. Step 2 and inverse discrete Fourier transform of the watermarked Fourier coefficient sequence to obtain watermarked point cloud data according to the inverse discrete Fourier transform Electronic information embedding method which is characterized by having 3 of the steps is obtained.
Preferably, the plurality of vertex coordinates are represented by vertex coordinates (x, y, z) obtained according to three-dimensional measurement.
More preferably, the real part of the discrete Fourier coefficients C k when 1 dimensional complex discrete Fourier transform is performed one-dimensional complex Fourier transformation is x k, the imaginary part is y k C k = x k + iy k (k 1 is a non-negative integer, i is an imaginary unit), and one-dimensional complex discrete Fourier transform is performed.
ここでは、前記第１のステップは、前記オリジナルランダム点群データが規定されたｘ−ｙ面領域を予め規定された小領域に分割して該小領域毎に含まれる点群を生成する第４のステップと、前記点群毎にそのｘ，ｙ座標値を前記点群の重心が原点となるようにオフセットしてオフセット点群に変換する第５のステップと、前記オフセット点群毎に前記離散フーリエ変換を行ってフーリエ係数列を求める第６のステップとを有している。
また、前記第３のステップは、前記透かし入りフーリエ係数列を逆離散フーリエ変換して透かし入り複素数列を生成する第７のステップと、該透かし入り複素数列について埋め込みに起因する座標値誤差のトレランスを満たす最適透かし埋め込み強度を求める第８のステップと、前記最適埋め込み強度に基づいて前記フーリエ係数列を再度変更して透かし入りフーリエ係数列とする第９のステップと、該透かし入りフーリエ係数列を逆オフセットして前記透かし入り点群データとする第１０のステップとを有している。
さらに、本発明では、前述のようにして得られた透かし入り点群データから前記電子透かしデータを抽出する際、前記オリジナルランダム点群データと前記透かし入り点群データとを対応付けてそれぞれ離散型フーリエ変換を行って第１及び第２のフーリエ係数列を求める第１１のステップと、前記第１及び前記第２のフーリエ係数列を比較して前記電子透かしデータを抽出する第１２のステップとを有することを特徴とする電子情報抽出方法が得られる。
Here, in the first step, the xy plane region in which the original random point cloud data is defined is divided into predetermined small regions, and a point group included in each small region is generated. A fifth step of offsetting the x and y coordinate values of each point group so that the center of gravity of the point group is the origin, and converting them into offset point groups, and the discrete points for each offset point group And a sixth step of obtaining a Fourier coefficient sequence by performing a Fourier transform.
The third step includes a seventh step of generating a watermarked complex number sequence by performing an inverse discrete Fourier transform on the watermarked Fourier coefficient sequence, and a tolerance of a coordinate value error caused by embedding the watermarked complex number sequence. An eighth step for obtaining an optimum watermark embedding strength satisfying the following: a ninth step for changing the Fourier coefficient sequence again based on the optimum embedding strength to obtain a watermarked Fourier coefficient sequence; and And a tenth step of reversely offsetting the watermarked point cloud data.
Further, in the present invention, when the digital watermark data is extracted from the watermarked point cloud data obtained as described above, the original random point cloud data and the watermarked point cloud data are associated with each other in a discrete type. An eleventh step of performing Fourier transform to obtain first and second Fourier coefficient sequences, and a twelfth step of extracting the digital watermark data by comparing the first and second Fourier coefficient sequences. The electronic information extraction method characterized by having it is obtained.
ここでは、前記第１１のステップは、前記オリジナルランダム点群データが規定されたｘ−ｙ面領域を予め規定された小領域に分割して該小領域毎に含まれる小領域点群を生成する第１３のステップと、前記小領域点群毎に該小領域点群の頂点と最短距離にある頂点を最短距離頂点として前記透かし入り点群データから探索して前記オリジナルランダム点群データと前記透かし入り点群データとの対応付けを行う第１４のステップとを有している。
例えば、前記第１４のステップは、前記透かし入り点群データに対する２−ｄ木を生成する第１５のステップと、前記小領域点群の頂点位置と埋め込みトレランスとによって規定される問い合わせ領域を設定して前記２−ｄ木内から前記問い合わせ領域に含まれる前記透かし入り点群データを探索して前記オリジナルランダム点群データと前記透かし入り点群データとの対応付けを行う第１６のステップとを有している。
Here, in the eleventh step, the xy plane area in which the original random point cloud data is defined is divided into predetermined small areas to generate small area point groups included in the small areas. A thirteenth step, searching for the original random point group data and the watermark by searching from the watermarked point group data with the vertex at the shortest distance from the vertex of the small area point group as the shortest distance vertex for each small area point group; A fourteenth step of associating with the entry point group data.
For example, in the fourteenth step, an inquiry region defined by a fifteenth step of generating a 2-d tree for the watermarked point cloud data, a vertex position of the small region point cloud, and an embedding tolerance is set. And searching for the watermarked point cloud data included in the inquiry area from the 2-d tree and associating the original random point cloud data with the watermarked point cloud data. ing.
また、本発明によれば、複数の頂点座標（ｘ，ｙ，…）で表され、ベクトル型に記述される不規則なデータ集合であるオリジナルランダム点群データに電子情報を電子透かしデータとして埋め込む際に用いられる電子情報埋め込み装置であって、複数の頂点座標２つが特に規則的配置を与えられていない前記オリジナルランダム点群データについて前記２つの頂点座標の一方を実数部に他方を虚数部にして１次元的に順序付けて１次元複素離散フーリエ変換を行ってフーリエ係数列を求める離散フーリエ変換手段と、前記電子透かしデータに応じて前記フーリエ係数列を変更して透かし入りフーリエ係数列とする変更手段と、前記透かし入りフーリエ係数列を逆離散フーリエ変換して該逆離散フーリエ変換に応じて透かし入り点群データを得る透かし入り点群データ生成手段とを有することを特徴とする電子情報埋め込み装置が得られる。
好ましくは、前記複数の頂点座標が、３次元計測に応じて得られた頂点座標（ｘ，ｙ，ｚ）で表されることを特徴とする。
更に好ましくは、前記１次元複素離散フーリエ変換が１次元複素フーリエ変換を実施した際の離散フーリエ係数Ｃ ｋ の実数部がｘ ｋ 、虚数部がｙ ｋ であるＣ ｋ ＝ｘ ｋ ＋ｉｙ ｋ （ｋは非負整数、ｉは虚数単位）なる複素数列に関する１次元複素離散フーリエ変換を行うことを特徴とする。
ここでは、前記離散フーリエ変換手段は、前記オリジナルランダム点群データが規定されたｘ−ｙ面領域を予め規定された小領域に分割して該小領域毎に含まれる点群を生成する分割手段と、前記点群毎にそのｘ，ｙ座標値を前記点群の重心が原点となるようにオフセットしてオフセット点群に変換するオフセット手段と、前記オフセット点群毎に前記離散フーリエ変換を行ってフーリエ係数列を求めるフーリエ係数生成手段とを有している。
さらに、前記透かし入り点群データ生成手段は、前記透かし入りフーリエ係数列を逆離散フーリエ変換して透かし入り複素数列を生成する複素数列生成手段と、該透かし入り複素数列について埋め込みに起因する座標値誤差のトレランスを満たす最適透かし埋め込み強度を求める透かし埋め込み強度算出手段と、前記最適埋め込み強度に基づいて前記フーリエ係数列を再度変更して透かし入りフーリエ係数列とする再変更手段と、該透かし入りフーリエ係数列を逆オフセットして前記透かし入り点群データとする逆オフセット手段とを有している。
さらに、本発明によれば、上述のようにして得られた透かし入り点群データから前記電子透かしデータを抽出する際、前記オリジナルランダム点群データと前記透かし入り点群データとを対応付けてそれぞれ離散型フーリエ変換を行って第１及び第２のフーリエ係数列を求めるフーリエ係数生成手段と、前記第１及び前記第２のフーリエ係数列を比較して前記電子透かしデータを抽出する抽出手段とを有することを特徴とする電子情報抽出装置が得られる。
ここでは、前記フーリエ係数生成手段は、前記オリジナルランダム点群データが規定されたｘ−ｙ面領域を予め規定された小領域に分割して該小領域毎に含まれる小領域点群を生成する分割手段と、前記小領域点群毎に該小領域点群の頂点と最短距離にある頂点を最短距離頂点として前記透かし入り点群データから探索して前記オリジナルランダム点群データと前記透かし入り点群データとの対応付けを行う対応付け手段とを有している。
例えば、前記対応付け手段は、前記透かし入り点群データに対する２−ｄ木を生成する２−ｄ木生成手段と、前記小領域点群の頂点位置と埋め込みトレランスとによって規定される問い合わせ領域を設定して前記２−ｄ木内から前記問い合わせ領域に含まれる前記透かし入り点群データを探索して前記オリジナルランダム点群データと前記透かし入り点群データとの対応付けを行う頂点間対応付け手段とを有している。
Further, according to the present invention, electronic information is embedded as electronic watermark data in original random point cloud data that is an irregular data set represented by a plurality of vertex coordinates (x, y,...) And described in a vector type. In the electronic information embedding apparatus used in this case, one of the two vertex coordinates is used as a real part and the other is used as an imaginary part for the original random point cloud data in which a plurality of vertex coordinates are not particularly given a regular arrangement. Discrete Fourier transform means for obtaining a Fourier coefficient sequence by performing one-dimensional complex discrete Fourier transform by ordering one-dimensionally, and changing the Fourier coefficient sequence according to the digital watermark data to change it to a watermarked Fourier coefficient sequence And a discrete Fourier transform of the watermarked Fourier coefficient sequence to obtain watermarked point cloud data according to the inverse discrete Fourier transform. Electronic information embedding apparatus characterized by having a watermarked point group data generating means is obtained.
Preferably, the plurality of vertex coordinates are represented by vertex coordinates (x, y, z) obtained according to three-dimensional measurement.
More preferably, the real part of the discrete Fourier coefficients C k when 1 dimensional complex discrete Fourier transform is performed one-dimensional complex Fourier transformation is x k, the imaginary part is y k C k = x k + iy k (k 1 is a non-negative integer, i is an imaginary unit), and one-dimensional complex discrete Fourier transform is performed.
Here, the discrete Fourier transform means divides the xy plane area in which the original random point cloud data is defined into predetermined small areas, and generates a point group included in each small area. Offset means for offsetting the x and y coordinate values for each point group so that the center of gravity of the point group is the origin, and converting the offset point group to the offset point group, and performing the discrete Fourier transform for each offset point group And Fourier coefficient generation means for obtaining a Fourier coefficient sequence.
Further, the watermarked point cloud data generating means includes a complex sequence generating means for generating a watermarked complex sequence by performing inverse discrete Fourier transform on the watermarked Fourier coefficient sequence, and a coordinate value resulting from embedding the watermarked complex sequence A watermark embedding strength calculating means for obtaining an optimum watermark embedding strength satisfying an error tolerance; a re-changing means for changing the Fourier coefficient sequence again to a watermarked Fourier coefficient sequence based on the optimum embedding strength; and the watermarked Fourier Reverse offset means for reversely offsetting the coefficient sequence to obtain the watermarked point cloud data.
Further, according to the present invention, when the digital watermark data is extracted from the watermarked point cloud data obtained as described above, the original random point cloud data and the watermarked point cloud data are associated with each other. Fourier coefficient generation means for obtaining first and second Fourier coefficient sequences by performing discrete Fourier transform; and extraction means for extracting the digital watermark data by comparing the first and second Fourier coefficient sequences. An electronic information extraction apparatus characterized by having the electronic information extraction apparatus is obtained.
Here, the Fourier coefficient generation means divides the xy plane area in which the original random point cloud data is defined into predetermined small areas, and generates small area point groups included in the small areas. Dividing means, and for each of the small area point groups, search for the original random point group data and the watermarked points by searching from the watermarked point group data with the vertex that is the shortest distance from the vertex of the small area point group as the shortest distance vertex. And an association means for associating with the group data.
For example, the associating means sets a query area defined by a 2-d tree generating means for generating a 2-d tree for the watermarked point group data, a vertex position of the small area point group, and an embedding tolerance. A vertex-to-vertex association unit that searches the 2-d tree for the watermarked point cloud data included in the inquiry area and associates the original random point cloud data with the watermarked point cloud data. Have.
加えて、本発明によれば、複数の頂点座標（ｘ，ｙ，…）で表され、ベクトル型に記述される不規則なデータ集合であるオリジナルランダム点群データに電子情報を電子透かしデータとして埋め込む際に電子計算機で用いられる電子情報埋め込みプログラムであって、複数の頂点座標の２つが特に規則的配置を与えられていない前記オリジナルランダム点群データについて、前記２つの頂点座標の一方を実数部に他方を虚数部にして１次元的に順序付けて、１次元離散フーリエ変換を行ってフーリエ係数列を求める第１の手続きと、前記電子透かしデータに応じて前記フーリエ係数列を変更して透かし入りフーリエ係数列とする第２の手続きと、前記透かし入りフーリエ係数列を逆離散フーリエ変換して該逆離散フーリエ変換に応じて透かし入り点群データを得る第３の手続きとを有することを特徴とする電子情報埋め込みプログラムが得られる。
好ましくは、前記複数の頂点座標が、３次元計測に応じて得られた頂点座標（ｘ，ｙ，ｚ）で表されることを特徴とする。
更に好ましくは、前記１次元複素離散フーリエ変換が１次元複素フーリエ変換を実施した際の離散フーリエ係数Ｃ ｋ の実数部がｘ ｋ 、虚数部がｙ ｋ であるＣ ｋ ＝ｘ ｋ ＋ｉｙ ｋ （ｋは非負整数、ｉは虚数単位）なる複素数列に関する１次元複素離散フーリエ変換を行うことを特徴とする。
ここでは、前記第１の手続きは、前記オリジナルランダム点群データが規定されたｘ−ｙ面領域を予め規定された小領域に分割して該小領域毎に含まれる点群を生成する第４の手続きと、前記点群毎にそのｘ，ｙ座標値を前記点群の重心が原点となるようにオフセットしてオフセット点群に変換する第５の手続きと、前記オフセット点群毎に前記離散フーリエ変換を行ってフーリエ係数列を求める第６の手続きとを有している。
また、前記第３の手続きは、前記透かし入りフーリエ係数列を逆離散フーリエ変換して透かし入り複素数列を生成する第７の手続きと、該透かし入り複素数列について埋め込みに起因する座標値誤差のトレランスを満たす最適透かし埋め込み強度を求める第８の手続きと、前記最適埋め込み強度に基づいて前記フーリエ係数列を再度変更して透かし入りフーリエ係数列とする第９の手続きと、該透かし入りフーリエ係数列を逆オフセットして前記透かし入り点群データとする第１０の手続きとを有している。
さらに、本発明によれば、前述のようにして得られた透かし入り点群データから前記電子透かしデータを抽出する際、前記オリジナルランダム点群データと前記透かし入り点群データとを対応付けてそれぞれ離散型フーリエ変換を行って第１及び第２のフーリエ係数列を求める第１１の手続きと、前記第１及び前記第２のフーリエ係数列を比較して前記電子透かしデータを抽出する第１２の手続きとを有することを特徴とする電子情報抽出プログラムが得られる。
ここでは、前記第１１の手続きは、前記オリジナルランダム点群データが規定されたｘ−ｙ面領域を予め規定された小領域に分割して該小領域毎に含まれる小領域点群を生成する第１３の手続きと、前記小領域点群毎に該小領域点群の頂点と最短距離にある頂点を最短距離頂点として前記透かし入り点群データから探索して前記オリジナルランダム点群データと前記透かし入り点群データとの対応付けを行う第１４の手続きとを有している。
例えば、前記第１４の手続きは、前記透かし入り点群データに対する２−ｄ木を生成する第１５の手続きと、前記小領域点群の頂点位置と埋め込みトレランスとによって規定される問い合わせ領域を設定して前記２−ｄ木内から前記問い合わせ領域に含まれる前記透かし入り点群データを探索して前記オリジナルランダム点群データと前記透かし入り点群データとの対応付けを行う第１６の手続きとを有している。
In addition, according to the present invention, electronic information is added as electronic watermark data to original random point cloud data that is an irregular data set represented by a plurality of vertex coordinates (x, y,...) And described in a vector type. An electronic information embedding program used by an electronic computer for embedding, wherein one of the two vertex coordinates is a real part for the original random point cloud data in which two of the plurality of vertex coordinates are not given a regular arrangement. The other is an imaginary part and is one-dimensionally ordered, and a one-dimensional discrete Fourier transform is performed to obtain a Fourier coefficient sequence, and the Fourier coefficient sequence is changed according to the digital watermark data and watermarked. A second procedure for making a Fourier coefficient sequence, and an inverse discrete Fourier transform of the watermarked Fourier coefficient sequence and a watermark in accordance with the inverse discrete Fourier transform Electronic information embedding program; and a third procedure to obtain to obtain the point cloud data Ri.
Preferably, the plurality of vertex coordinates are represented by vertex coordinates (x, y, z) obtained according to three-dimensional measurement.
More preferably, the real part of the discrete Fourier coefficients C k when 1 dimensional complex discrete Fourier transform is performed one-dimensional complex Fourier transformation is x k, the imaginary part is y k C k = x k + iy k (k 1 is a non-negative integer, i is an imaginary unit), and one-dimensional complex discrete Fourier transform is performed.
Here, the first procedure generates a point group included in each small area by dividing the xy plane area in which the original random point cloud data is defined into predetermined small areas. A fifth procedure for offsetting the x and y coordinate values of each point cloud so that the center of gravity of the point cloud is the origin, and converting the offset into a point cloud, and the discrete for each offset point cloud A sixth procedure for performing a Fourier transform to obtain a Fourier coefficient sequence.
The third procedure includes a seventh procedure for generating a watermarked complex number sequence by performing an inverse discrete Fourier transform on the watermarked Fourier coefficient sequence, and a tolerance of a coordinate value error caused by embedding the watermarked complex number sequence. An eighth procedure for obtaining the optimum watermark embedding strength satisfying the above, a ninth procedure for changing the Fourier coefficient sequence again based on the optimum embedding strength to obtain a watermarked Fourier coefficient sequence, and the watermarked Fourier coefficient sequence: And a tenth procedure for reverse offset to obtain the watermarked point cloud data.
Further, according to the present invention, when the digital watermark data is extracted from the watermarked point cloud data obtained as described above, the original random point cloud data and the watermarked point cloud data are associated with each other. An eleventh procedure for obtaining first and second Fourier coefficient sequences by performing discrete Fourier transform and a twelfth procedure for extracting the digital watermark data by comparing the first and second Fourier coefficient sequences An electronic information extraction program characterized by having:
Here, in the eleventh procedure, the xy plane region in which the original random point cloud data is defined is divided into predetermined small regions to generate small region point groups included in the small regions. A thirteenth procedure, searching for the original random point group data and the watermark by searching from the watermarked point group data with the vertex at the shortest distance from the vertex of the small region point group as the shortest distance vertex for each small region point group; And a fourteenth procedure for associating with the entry point cloud data.
For example, in the fourteenth procedure, an inquiry region defined by a fifteenth procedure for generating a 2-d tree for the watermarked point cloud data, a vertex position of the small region point cloud, and an embedding tolerance is set. And searching for the watermarked point cloud data included in the inquiry area from the 2-d tree and associating the original random point cloud data with the watermarked point cloud data. ing.
このようにして、本発明では、オリジナルランダム点群データに電子透かしデータを埋め込んで、オリジナルランダム点群データの不正使用を防止することができる。さらに、上述のようにして、電子透かしデータを、透かし入り点群データから抽出して、別に管理しているオリジナルの電子透かしデータと比較するようにしたから、その結果に応じて透かし入り点群データがオリジナル点群データから生成されたものであるかを判定することができる。つまり、オリジナルランダム点群データの同一性を判定することもできる。 In this way, in the present invention, digital watermark data can be embedded in original random point cloud data to prevent unauthorized use of the original random point cloud data. Further, as described above, since the digital watermark data is extracted from the watermarked point cloud data and compared with the original digital watermark data managed separately, the watermarked point cloud is selected according to the result. It can be determined whether the data is generated from the original point cloud data. That is, it is possible to determine the identity of the original random point cloud data.
以下本発明について図面を参照して説明する。
図１を参照して、地球表面（例えば、地表面）に関するオリジナルランダム点群データを得る際には、まず、地球表面の３次元計測が行われる。地球表面の３次元計測に当たっては、例えば、レーザ３次元計測が用いられる。３次元レーザ計測を行う際には、航空機１１をオリジナルランダム点群データを収集すべき地域（地区）上を飛行させる。航空機１１にはレーザスキャナ（航空機搭載型レーザスキャナ：図示せず）、ＧＰＳ受信機（図示せず）、及びジャイロスコープ（ＩＭＵ：図示せず）が搭載されている。
図２も参照して、レーザ３次元計測（ステップＳ１）では、航空機搭載型レーザスキャナから地上（地表面）１２に向けてパルスレーザビームを照射して、地表面１２からの反射パルスをパルスリターンとして受け、パルスレーザビームの照射時刻とパルスリターンを受けた時刻との時間差に応じて地表面１２までの距離（ｚ）を計測して距離データを得る。このパルスレーザビームの照射は予め規定された時間間隔で行われる結果、地表面は離散的な点として規定されることになる。この際、レーザスキャナのスキャンミラーの回転は１軸（航空機１１の飛行方向に直交する）であり、飛行方向がもう一つの走査軸となって、面的にデータが取得されることになる。そして、ＩＭＵは１パルス毎にレーザ照射角度をミラー回転角及びその傾き角度に基づいて求める。
さらに、別に地上において、レーザスキャナ（地上据え置き型レーザスキャナ）を用いて地表面のデータを取得する。このレーザスキャナではスキャンミラーの回転軸は２軸で放射状にしかも面的にデータが得られることになる。
一方、ＧＰＳ電波に応じて基準点測量が行われる（ステップＳ２）。つまり、航空機１１の空間位置は、地上に設けられたＧＰＳ基準局１３、ＧＰＳ受信機、及びＩＭＵによって計測される。その後、レーザ３次元計測によって得られたデータ（スキャンポイントデータ）と基準点測量とによって得られたデータとに基づいてスキャンポイントデータを測地座標データへ変換してオリジナルランダム点群データを得る（ステップＳ３）。つまり、１パルス毎に地表面１２におけるレーザビーム反射点のｘ，ｙ，ｚ（高さ：Ｈ）座標を求めて、オリジナルランダム点群データとする。
そして、上述のオリジナルランダム点群データに電子透かしデータを埋め込み（ステップＳ４）、透かし情報（データ）埋め込み済み３次元空間座標点群として出力する（ステップＳ５）。なお、ステップＳ３〜Ｓ５はコンピュータ（図示せず）によって行われる。つまり、コンピュータはステップＳ３を行う測地座標変換部、ステップＳ４を行う電子透かし埋め込み部、及びステップＳ５を行う透かし入り点群出力部を少なくとも有しており、測地座標変換部では、スキャンポイントデータと基準点測量データとを読み込んでオリジナルランダム点群データを生成する。
The present invention will be described below with reference to the drawings.
Referring to FIG. 1, when obtaining original random point cloud data relating to the earth surface (for example, the earth surface), first, three-dimensional measurement of the earth surface is performed. For the three-dimensional measurement of the earth surface, for example, laser three-dimensional measurement is used. When performing three-dimensional laser measurement, the aircraft 11 is allowed to fly over an area (district) where original random point cloud data should be collected. The aircraft 11 is equipped with a laser scanner (aircraft laser scanner: not shown), a GPS receiver (not shown), and a gyroscope (IMU: not shown).
Referring also to FIG. 2, in the laser three-dimensional measurement (step S <b> 1), a pulse laser beam is irradiated from the airborne laser scanner toward the ground (ground surface) 12, and the reflected pulse from the ground surface 12 is pulse-returned. The distance data is obtained by measuring the distance (z) to the ground surface 12 according to the time difference between the irradiation time of the pulse laser beam and the time of receiving the pulse return. The irradiation of the pulse laser beam is performed at a predetermined time interval. As a result, the ground surface is defined as discrete points. At this time, the rotation of the scan mirror of the laser scanner is one axis (perpendicular to the flight direction of the aircraft 11), and the flight direction becomes another scan axis, and data is acquired in a plane. The IMU obtains the laser irradiation angle for each pulse based on the mirror rotation angle and the tilt angle.
Further, on the ground, ground surface data is acquired using a laser scanner (a stationary laser scanner on the ground). In this laser scanner, the rotation axis of the scan mirror is radial with two axes, and data can be obtained in a plane.
On the other hand, reference point surveying is performed according to GPS radio waves (step S2). That is, the spatial position of the aircraft 11 is measured by the GPS reference station 13, the GPS receiver, and the IMU provided on the ground. Thereafter, the original random point cloud data is obtained by converting the scan point data into geodetic coordinate data based on the data (scan point data) obtained by the laser three-dimensional measurement and the data obtained by the reference point survey (step) S3). That is, for each pulse, the x, y, z (height: H) coordinates of the laser beam reflection point on the ground surface 12 are obtained and used as original random point cloud data.
Then, the electronic watermark data is embedded in the original random point cloud data (step S4), and is output as a watermark information (data) embedded 3-dimensional space coordinate point cloud (step S5). Steps S3 to S5 are performed by a computer (not shown). That is, the computer has at least a geodetic coordinate conversion unit that performs step S3, a digital watermark embedding unit that performs step S4, and a watermarked point group output unit that performs step S5. In the geodetic coordinate conversion unit, scan point data and The reference point survey data is read to generate original random point cloud data.
ここで、図３を参照して、電子透かし埋め込み部の動作について説明する。
図示の例では、電子透かし埋め込み部は、小領域分割部、座標値オフセット部、離散フーリエ変換部、フーリエ係数列変更部、逆離散フーリエ変換部、埋め込み強度自動調整部、及び逆オフセット部を有している。
電子透かしデータ埋め込み処理の際には、前述の測地座標変換部からオリジナルランダム点群データ（以下単にオリジナル点群と呼ぶ）Ｖが小領域分割部に与えられるとともに別に入力装置（図示せず）等から小領域ｘ，ｙ方向サイズＤｘ，Ｄｙが小領域分割部に与えられる。
ここでは、オリジナル点群Ｖは
で表されるものとする。
ｖｉ：頂点３次元座標，Ｎ：全頂点数．
●Ｄｘ，Ｄｙ：分割小領域のx, y方向サイズ （標準値Dx=100[m],Dy=100[m]）
小領域分割部では、まず、オリジナル点群の小領域への分割を行う（ステップＡ１）。図４に示すように、オリジナル点群Ｖが定義されているｘ−ｙ面内の矩形領域を、Ｄｘ，Ｄｙのサイズ（寸法）を有する小領域Ａａに均等に分割して（等分割）、小領域Ａａ毎に含まれる点群Ｖａを生成する。そして、この点群Ｖａは座標オフセット部に与えられる。
ＬＡ：全小領域数
なお、図４において、その下辺及び右辺のように、ｘ方向及びｙ方向へ分割した際、その最終端では、Ｄｘ，Ｄｙよりも小さいサイズの小領域（この領域に含まれる点群をＶＬａとする）が形成されることもある。
座標オフセット部では、小領域内の点群のｘ，ｙ座標値のオフセットを行う（ステップＡ２）。ここでは、各小領域の点群Ｖａのｘ，ｙ座標値を、Ｖａの重心が原点となるようにオフセットした点群
へ変換する。つまり、オフセット点群
は、
μａｘ，μａｙ：点群Ｖａの重心点x, y座標値
上述のオフセット処理によって、オリジナル点群Ｖのｘ，ｙ座標値の絶対値が大きい場合
に、埋め込み処理で用いられる浮動小数点演算における丸め誤差の影響を減少させること
ができる。
このオフセット点群
は離散フーリエ変換部及び埋め込み強度自動調整部に与えられ、一方、小領域重心点μａｘ，μａｙは逆オフセット部に与えられる。
離散フーリエ変換部は、オフセット点群
に対して離散フーリエ変換（ＤＦＴ）を行う（ステップＡ３）。ここでは、まず、オフセット点群
に含まれる点のｘ，ｙ座標値から、式（４）で示す複素数列｛ｃｋ｝を生成する。
i:虚数単位
次に、複素数列｛ｃｋ｝を式（５）に基づいて、離散フーリエ変換部は、離散フーリエ変換して，フーリエ係数列｛Ｃｌ｝を求める。
そして、このフーリエ係数列｛Ｃｌ｝はフーリエ係数列変更部に与えられる。フーリエ係数列変更部には電子透かしデータビット列Ｂ及び埋め込み強度初期値Ｐｉｎｉｔが与えられており、フーリエ係数列変更部では、電子透かしデータビット列Ｂ及び埋め込み強度初期値Ｐｉｎｉｔに応じてフーリエ係数列｛Ｃｌ｝を変更する（ステップＡ４）。
ここで、電子透かしデータビット列Ｂは
●Ｂ＝｛ｂｍ∈｛０，１｝｜１＜ｍ＜Ｎｂ｝で表される。
●Ｎｂ：透かしデータビット長
●ｐｉｎｉｔ：透かし埋め込み強度初期値（標準値ｐｉｎｉｔ＝５０００）
つまり、電子透かしデータビット列Ｂ＝｛ｂｍ∈｛０，１｝｜１≦ｍ≦Ｎｂ｝及び埋め込み強度初期値ｐｉｎｉｔを用いて、フーリエ係数列変更部は、式（６）に従ってフーリエ係数列｛Ｃｌ｝を変更して、変更後フーリエ係数列｛Ｃ’ｌ｝を生成する。
なお、直流成分を表す係数Ｃ０は透かしデータによって変更しない。従って、この小領域Ａａ内に埋め込み可能な透かしビット長さは、｜Ｖａ｜−１（ビット）となる。また、埋め込み強度ｐの値は、１回目の埋め込みの際には、ｐ←ｐｉｎｉｔと定め、２回目の埋め込みの際には、ｐ←ｐｏｐｔと定める。
この際、透かしビット長Ｎｂと小領域への埋め込み可能ビット長｜Ｖａ｜−１との関係に基づいて次の規則を適用される。
Here, the operation of the digital watermark embedding unit will be described with reference to FIG.
In the illustrated example, the digital watermark embedding unit includes a small region dividing unit, a coordinate value offset unit, a discrete Fourier transform unit, a Fourier coefficient sequence changing unit, an inverse discrete Fourier transform unit, an embedding strength automatic adjustment unit, and an inverse offset unit. is doing.
In the digital watermark data embedding process, original random point cloud data (hereinafter simply referred to as original point cloud) V is given from the above-described geodetic coordinate conversion unit to the small area dividing unit and an input device (not shown) or the like. The small area x and the y-direction sizes Dx and Dy are given to the small area dividing unit.
Here, the original point cloud V is
It shall be represented by
v i : vertex three-dimensional coordinates, N: total number of vertices.
● D x , D y : Size of divided small area in x and y direction (standard values Dx = 100 [m], Dy = 100 [m])
In the small area dividing unit, first, the original point cloud is divided into small areas (step A1). As shown in FIG. 4, the rectangular area in the xy plane in which the original point group V is defined is equally divided into small areas Aa having the sizes (dimensions) of Dx and Dy (equal division). A point group Va included in each small area Aa is generated. This point group Va is given to the coordinate offset portion.
L A : the number of all small areas In FIG. 4, when divided in the x direction and the y direction, such as the lower side and the right side, at the final end, a small area having a size smaller than D x and D y (this A point group included in the region may be defined as V La ).
In the coordinate offset unit, the x and y coordinate values of the point group in the small area are offset (step A2). Here, a point group in which the x and y coordinate values of the point group Va of each small region are offset so that the center of gravity of Va is the origin.
Convert to In other words, offset point cloud
Is
μ ax , μ ay : barycentric point x, y coordinate value of point group V a
By the offset processing described above, when the absolute values of the x and y coordinate values of the original point group V are large, it is possible to reduce the influence of the rounding error in the floating point calculation used in the embedding processing.
This offset point cloud
Is given to the discrete Fourier transform unit and the embedding strength automatic adjustment unit, while the small area centroid points μ ax and μ ay are given to the inverse offset unit.
The discrete Fourier transform unit uses the offset point group
Discrete Fourier transform (DFT) is performed on (step A3). Here, first, the offset point cloud
From the x and y coordinate values of the points included in the complex number sequence {c k } shown in Expression (4).
i: Imaginary Unit Next, based on the equation (5), the discrete Fourier transform unit performs a discrete Fourier transform on the complex sequence {c k } to obtain a Fourier coefficient sequence {C l }.
This Fourier coefficient sequence {C 1 } is given to the Fourier coefficient sequence changing unit. The Fourier coefficient sequence changing unit is given a digital watermark data bit sequence B and an embedding strength initial value Pinit , and the Fourier coefficient sequence changing unit has a Fourier coefficient sequence according to the digital watermark data bit sequence B and the embedding strength initial value Pinit. {C l } is changed (step A4).
Here, the digital watermark data bit string B is expressed by: B = {b m ε {0, 1} | 1 <m <N b }.
● N b: watermark data bit length ● p init: watermark embedding intensity initial value (standard value p init = 5000)
That is, using the digital watermark data bit sequence B = {b m ε {0,1} | 1 ≦ m ≦ N b } and the embedding strength initial value p init , the Fourier coefficient sequence changing unit performs the Fourier coefficient according to Equation (6). The sequence {C l } is changed to generate a post-change Fourier coefficient sequence {C ′ l }.
Note that the coefficient C 0 representing the DC component is not changed by the watermark data. Therefore, the watermark bit length that can be embedded in the small area Aa is | V a | −1 (bits). Further, the value of the embedding strength p is defined as p ← p init at the first embedding, and p ← p opt at the second embedding.
At this time, implantable bit length of the watermark bit length N b and a small area | V a |, based on the relationship between -1 apply the following rules.
ｉ）Ｎｂ＜｜Ｖａ｜−１の場合
透かしビットを反復して埋め込み，Ｃ１，Ｃ２，．．．．，Ｃ｜Ｖａ｜−１のフーリエ係数全ての値を前述の式（６）で変更する。
ｉ）Ｎｂ＞｜Ｖａ｜−１の場合
透かしビットの先頭から｜Ｖａ｜−１ビット分のみを埋め込み，それ以降のビットについては埋め込みを行わない。
上述の変更後フーリエ係数列｛Ｃ’ｌ｝、つまり、透かし入りフーリエ係数列｛Ｃ’ｌ｝は逆離散フーリエ変換部に与えられ、逆離散フーリエ変換部では、変更後フーリエ係数列｛Ｃ’ｌ｝に対して逆離散フーリエ変換（ＩＤＦＴ）を行う（ステップＡ５）。つまり、逆離散フーリエ変換部では、透かし入りフーリエ係数列｛Ｃ’ｌ｝に対して、式（７）の逆フーリエ変換を行って、電子透かしによって変更された複素数列｛ｃ’ｋ｝を生成する。
この複素数列｛ｃ’ｋ｝は、埋め込み強度自動調整部に与えられる。埋め込み強度自動調整部には、前述のオフセット点群
及び埋め込みトレランスτが与えられており、埋め込み強度自動調整部では、埋め込みトレランスτを考慮して、埋め込み強度を自動調整する（ステップＡ６）。
●τ：ｘ，ｙ方向の埋め込みトレランス (標準値τ＝０．３[ｍ])
i) If N b <| V a | −1, the watermark bits are repeatedly embedded and C 1 , C 2 ,. . . . , C | V a | −1 all the Fourier coefficients are changed by the above-described equation (6).
i) In the case of N b > | V a | −1, only | V a | −1 bits are embedded from the beginning of the watermark bit, and the subsequent bits are not embedded.
The modified Fourier coefficient sequence {C ′ l } described above, that is, the watermarked Fourier coefficient sequence {C ′ l } is given to the inverse discrete Fourier transform unit, and the inverse discrete Fourier transform unit uses the modified Fourier coefficient sequence {C ′ l } is subjected to inverse discrete Fourier transform (IDFT) (step A5). That is, the inverse discrete Fourier transform unit performs the inverse Fourier transform of the equation (7) on the watermarked Fourier coefficient sequence {C ′ l } to generate a complex sequence {c ′ k } modified by the digital watermark. To do.
The complex number sequence {c ′ k } is given to the embedding strength automatic adjustment unit. The embedding strength automatic adjustment unit includes the offset point group described above.
And an embedding tolerance τ, and the embedding strength automatic adjustment unit automatically adjusts the embedding strength in consideration of the embedding tolerance τ (step A6).
● τ: Embedding tolerance in the x and y directions (standard value τ = 0.3 [m])
つまり、埋め込み強度自動調整部では、変更後複素数列｛ｃ’ｋ｝に基づいて透かし入り点群データ
の座標値を知り、透かし埋め込み前の点群
の頂点
と埋め込み後の点群
の頂点
との間において、ｘ方向とｙ方向のうちの最大誤差を
とし，その最大誤差を与える頂点番号
を求める。即ち、
となる。
ここで，埋め込みによって生ずるｘ，ｙ方向の頂点座標値誤差のトレランス（許容値）をτ，このトレランスを満たす最適な埋め込み強度をｐｏｐｔとすると、式（９）に示す比例関係が成り立つ．
従って，式（９）から、埋め込み誤差のトレランスτを満たす最適な透かし埋め込み強度ｐｏｐｔは、
となる。そして、この最適な透かし埋め込み強度ｐｏｐｔは、フーリエ係数列変更部にフィードバックされる。
フーリエ係数列変更部では、最適埋め込み強度ｐｏｐｔを用いてフーリエ係数を変更し、逆フーリエ変換部ではその変更結果に応じて逆フーリエ変換を再度実行する。つまり、前述のようにして求めた最適埋め込み強度ｐｏｐｔを用いて、前述のステップＡ４及びＡ５を再度実行して、電子透かし入りオフセット点群
を生成する。
電子透かし入りオフセット点群
は逆オフセット部に与えられる。前述のように、逆オフセット部には、小領域重心点μａ
ｘ，μａｙが与えられており、逆オフセット部は、ｘ，ｙ座標値の逆オフセットを行って
透かし入り点群データＶ’ａを生成する（ステップＡ７）。ここでは、前述の式（２）の逆操作を行って、つまり、逆オフセットを式（１１）に基づいて行って、透かし入り点群
から透かし入り点群Ｖ’ａを求める。
上述のようにして、電子透かし埋め込み部では、各小領域毎に処理を行い、最終的に透かし入り点群Ｖ’を生成する。この透かし入り点群は、透かし入り点群出力部によって透かし情報埋め込み済み空間３次元座標点群としてファイルに出力される。
図５を参照して、透かし入り点群データＶ’から透かしデータを抽出する場合について説明する。ここでは、コンピュータは透かしデータ抽出部を有しており、透かしデータ抽出部には点群小領域分割部、頂点間対応付け部、離散フーリエ変換部、及び透かしビット列抽出部が備えられている。
いま、オリジナル点群Ｖを、
とし、透かし入り点群Ｖ’を、
とする。
この際、前述のＤｘ，Ｄｙ（分割小領域のｘ，ｙ方向サイズ）及びτ（ｘ，ｙ方向の埋め込みトレランス）は電子透かしデータ埋め込みの際と同一の値を用いる。
That is, the embedding strength automatic adjustment unit performs watermarked point cloud data based on the changed complex number sequence {c ′ k }.
Know the coordinate value of the point cloud before embedding the watermark
Vertex of
And point cloud after embedding
Vertex of
Between the x direction and the y direction
And the vertex number giving the maximum error
Ask for. That is,
It becomes.
Here, if the tolerance (allowable value) of the vertex coordinate value error in the x and y directions caused by embedding is τ, and the optimum embedding strength satisfying this tolerance is p opt , the proportional relationship shown in Equation (9) is established.
Therefore, from equation (9), the optimum watermark embedding strength p opt that satisfies the embedding error tolerance τ is:
It becomes. The optimum watermark embedding strength p opt is fed back to the Fourier coefficient sequence changing unit.
The Fourier coefficient sequence changing unit changes the Fourier coefficient using the optimum embedding strength p opt , and the inverse Fourier transform unit executes the inverse Fourier transform again according to the change result. That is, the above-described steps A4 and A5 are executed again using the optimum embedding strength p opt obtained as described above, and the watermark-added offset point group
Is generated.
Digital watermarked offset point cloud
Is given to the reverse offset part. As described above, the reverse offset portion includes the small area centroid point μ a.
x and μ ay are given, and the inverse offset unit performs inverse offset of the x and y coordinate values to generate watermarked point group data V ′ a (step A7). Here, the inverse operation of the above-described equation (2) is performed, that is, the inverse offset is performed based on the equation (11), and the watermarked point cloud
Is used to obtain a watermarked point group V′a.
As described above, the digital watermark embedding unit performs processing for each small region, and finally generates a watermarked point group V ′. The watermarked point group is output to a file as a watermark information embedded space three-dimensional coordinate point group by the watermarked point group output unit.
A case where watermark data is extracted from watermarked point cloud data V ′ will be described with reference to FIG. Here, the computer has a watermark data extraction unit, and the watermark data extraction unit includes a point group small region division unit, an inter-vertex association unit, a discrete Fourier transform unit, and a watermark bit string extraction unit.
Now, the original point cloud V
And a watermarked point cloud V ′,
And
At this time, D x , D y (size in the x and y directions of the divided small area) and τ (embedding tolerance in the x and y directions) are the same values as those used when embedding the digital watermark data.
まず、オリジナル点群Ｖ及び小領域ｘ，ｙ方向サイズＤｘ，Ｄｙが点群小領域分割部に与えられ、点群小領域分割部はオリジナル点群Ｖを小領域に分割する（ステップＡ８）。
前述の透かし埋め込み時に指定した分割小領域のｘ，ｙ方向サイズＤｘ，Ｄｙの値を用いて、埋め込み時と同様にしてオリジナル点群Ｖを小領域に分割して、領域毎に点群Ｖａを求める。
なお、透かし入り点群Ｖ’は分割を行うことなく、そのまま保存される。
小領域内点群Ｖは頂点間対応付け部に与えられる。この頂点間対応付け部には、前述の埋め込みトレランスτ及び透かし入り点群Ｖ’が与えられており、頂点間対応付け部では、小領域内点群Ｖと透かし入り点群Ｖ’との頂点間の対応づけを行って小領域内透かし入り点群Ｖ’ａを出力する（ステップＡ９）。
例えば、頂点間対応付け部では、各小領域群Ｖａ毎に点群Ｖａ内の頂点ｖｉ（∈Ｖａ）と最短距離にある頂点を透かし入り点群Ｖ’から探索して、この探索の結果得られた頂点を、頂点ｖｉに対応付けられた透かし入り頂点ｖ’ｉ（∈Ｖ’）として採用する。
上述の対応付けにおいて、点群Ｖ’は極めて多数の頂点を含む関係上、総当たりで最短頂点を探索すると、点数の二乗のオーダーで探索時間が必要となる。このため、図示のように、２−ｄ木生成部を備えて、前処理として透かし入り点群Ｖ’に対する２−ｄ木を生成するようにしてもよい（ステップＡ１０）。そして、透かし入り点群Ｖ’の代わりに、透かし入り点群の２−ｄ木を頂点間対応付け部に与えるようにしてもよい。
そして、オリジナル点群内の頂点ｖｉ（∈Ｖａ）の位置と埋め込みトレランスτとによって定まる微小な問い合わせ領域を設定して、２−ｄ木内から問い合わせ領域に含まれる透かし入り頂点群を高速に探索して、これらの中だけで総当たりによって最短頂点を探索する。これによって、オリジナル点群Ｖｉ（∈Ｖａ）と透かし入り頂点ｖ’ｉ（∈Ｖ’）との対応付け処理を効率化することができる。
このようにして、小領域に含まれるオリジナル点群Ｖａに透かし入り点群Ｖ’ａを対応付けた後、離散フーリエ変換部では、オリジナル点群Ｖａ及び透かし入り点群Ｖ’ａに対してＤＦＴを実行する（ステップＡ１１）。オリジナル点群及びこれに対応づけられた透かし入り点群に、埋め込み時と同様にして、式（４）及び（５）で示すＤＦＴを行って、それぞれフーリエ係数列｛Ｃｌ｝及び｛Ｃ’ｌ｝を求める。
そして、透かしビット列抽出部では、フーリエ係数列｛Ｃｌ｝及び透かし入りフーリエ係数列｛Ｃ’ｌ｝に応じて電子透かしデータの抽出を行う（ステップＡ１２）。つまり、上述のようにして得られたフーリエ係数列｛Ｃｌ｝及び｛Ｃ’ｌ｝を比較して、埋め込まれた透かしデータのビット列
を式（１２）に基づいて抽出する。
このようにして、埋め込まれた透かしビット列
を式（１２）に基づいて抽出して、別に管理していたオリジナルの透かしビット列と比較する。この比較結果に基づいて透かし入り点群データがオリジナル点群データから生成されたものであるかを判定する。
First, the original point group V and the small region x, y-direction sizes D x and D y are given to the point group small region dividing unit, and the point group small region dividing unit divides the original point group V into small regions (step A8). ).
The original point cloud V is divided into small areas in the same manner as in the embedding using the values of the x, y direction sizes D x and D y of the divided small areas specified at the time of embedding the watermark described above, and the point cloud for each area. determine the V a.
The watermarked point group V ′ is stored as it is without being divided.
The small area in-point group V is given to the inter-vertex association unit. The inter-vertex association unit is provided with the above-described embedding tolerance τ and the watermarked point group V ′. In the inter-vertex association unit, the vertexes of the small-area in-point group V and the watermarked point group V ′ are provided. association and outputs a small area within the watermarked point group V 'a go between (step A9).
For example, the vertex association unit, in search from the vertex v i (∈Va) and watermarked point group vertices in the shortest distance V 'in the point group V a for each sub-area groups V a, the search the resulting vertex is employed as vertex v i to the associated watermarked vertex v 'i (∈V').
In the above-described association, the point group V ′ includes a very large number of vertices. Therefore, when searching for the shortest vertex in the round robin, the search time is required in the order of the square of the number of points. For this reason, as shown in the figure, a 2-d tree generation unit may be provided to generate a 2-d tree for the watermarked point group V ′ as preprocessing (step A10). Then, instead of the watermarked point group V ′, a 2-d tree of the watermarked point group may be given to the inter-vertex association unit.
Then, a minute inquiry area determined by the position of the vertex v i (∈V a ) in the original point group and the embedding tolerance τ is set, and the watermarked vertex group included in the inquiry area is quickly moved from the 2-d tree. Search and search for the shortest vertex by brute force only in these. As a result, the process of associating the original point group V i (∈V a ) with the watermarked vertex v ′ i (∈V ′) can be made efficient.
In this way, after associating the watermarked point group V ′ a with the original point group V a included in the small region, the discrete Fourier transform unit applies the original point group V a and the watermarked point group V′a. The DFT is executed (step A11). Similar to the case of embedding, the original point cloud and the watermarked point cloud associated with the original point cloud are subjected to DFT represented by equations (4) and (5), respectively, and Fourier coefficient sequences {C 1 } and {C ′ l }.
The watermark bit string extraction unit extracts digital watermark data in accordance with the Fourier coefficient string {C 1 } and the watermarked Fourier coefficient string {C ′ 1 } (step A12). That is, by comparing the Fourier coefficient sequences {C 1 } and {C ′ 1 } obtained as described above, the bit sequence of the embedded watermark data
Is extracted based on the equation (12).
In this way, the embedded watermark bit string
Is extracted based on Expression (12) and compared with the original watermark bit string managed separately. Based on the comparison result, it is determined whether the watermarked point cloud data is generated from the original point cloud data.
本発明によれば、オリジナルランダム点群データの不正使用を防止することのできる電子情報埋め込み方法、電子情報の埋め込みに用いられる装置、及びプログラムを提供することが出来る。
又本発明によれば、電子情報が埋め込まれたオリジナルランダム点群データから電子情報を抽出するための電子情報抽出方法、電子情報の抽出に用いられる装置、及びプログラムを提供出来る。
According to the present invention, it is possible to provide an electronic information embedding method capable of preventing unauthorized use of original random point cloud data, an apparatus used for embedding electronic information, and a program.
In addition, according to the present invention, it is possible to provide an electronic information extraction method for extracting electronic information from original random point cloud data in which electronic information is embedded, an apparatus used for extracting electronic information, and a program.
Claims (24)
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2007040094A JP4175551B2 (en) | 2007-02-20 | 2007-02-20 | Electronic information embedding method and extraction method, electronic information embedding device and extraction device, and program thereof |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2007040094A JP4175551B2 (en) | 2007-02-20 | 2007-02-20 | Electronic information embedding method and extraction method, electronic information embedding device and extraction device, and program thereof |
Related Parent Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2002586705A Division JP3935078B2 (en) | 2001-12-12 | 2001-12-12 | Aviation Electronic Map Information Creation Method and Aviation Electronic Map Information Creation Device |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JP2007189724A JP2007189724A (en) | 2007-07-26 |
| JP4175551B2 true JP4175551B2 (en) | 2008-11-05 |
Family
ID=38344525
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2007040094A Expired - Fee Related JP4175551B2 (en) | 2007-02-20 | 2007-02-20 | Electronic information embedding method and extraction method, electronic information embedding device and extraction device, and program thereof |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP4175551B2 (en) |
-
2007
- 2007-02-20 JP JP2007040094A patent/JP4175551B2/en not_active Expired - Fee Related
Also Published As
| Publication number | Publication date |
|---|---|
| JP2007189724A (en) | 2007-07-26 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Eltner et al. | Analysis of different methods for 3D reconstruction of natural surfaces from parallel‐axes UAV images | |
| Verhoeven et al. | Computer vision‐based orthophoto mapping of complex archaeological sites: The ancient quarry of Pitaranha (Portugal–Spain) | |
| US7891006B2 (en) | Electronic information embedding method and extracting method, electronic information burying apparatus and extracting apparatus, and programs therefor | |
| JP2006214893A (en) | Computer software program for measuring three-dimensional shape of object using object measuring method and computer system | |
| Javadnejad et al. | An assessment of UAS-based photogrammetry for civil integrated management (CIM) modeling of pipes | |
| Postolov et al. | Registration of airborne laser data to surfaces generated by photogrammetric means | |
| Pavan et al. | A global closed-form refinement for consistent TLS data registration | |
| JP6945832B2 (en) | 3D shape search method and 3D shape search system | |
| Arias et al. | ORTHOIMAGE‐BASED DOCUMENTATION OF ARCHAEOLOGICAL STRUCTURES: THE CASE OF A MEDIAEVAL WALL IN PONTEVEDRA, SPAIN | |
| JP4175551B2 (en) | Electronic information embedding method and extraction method, electronic information embedding device and extraction device, and program thereof | |
| Toschi et al. | Improving automated 3D reconstruction methods via vision metrology | |
| EP4036859A1 (en) | A system and method for providing improved geocoded reference data to a 3d map representation | |
| Murtiyos et al. | Multi-sensor 3D recording pipeline for the documentation of Javanese Temples | |
| KR101459004B1 (en) | Method for converting 3D Image based plan to 3D Image based spherical surface | |
| Wang et al. | A novel capacity expansion and recognition acceleration method for dot-dispersing coded targets in photogrammetry | |
| Liu et al. | A novel approach to automatic registration of point clouds | |
| Hämmerle et al. | Terrestrial and unmanned aerial system imagery for deriving photogrammetric three-dimensional point clouds and volume models of mass wasting sites | |
| Ivanova | Adaptive digitalization methods and digital transformation trends for security. | |
| Rozas et al. | Comparative study of building footprint estimation methods from LiDAR point clouds | |
| KR20120138606A (en) | A system for generating road layer using continuity analysis of ground data | |
| JP7122901B2 (en) | breakline extraction program, breakline extraction system | |
| Ponomarev et al. | Automatic structural matching of 3D image data | |
| Zhou et al. | A pyramid image matching algorithm for chang'e-2 CCD images of the moon | |
| AL-Qadri et al. | Accuracy Assessment for Geometric Features Extraction from Multirotor UAV's Images | |
| Rzonca et al. | Lidarometry as a Variant of Integration of Photogrammetric and Laser Scanning Data |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20080425 |
|
| A521 | Written amendment |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20080624 |
|
| TRDD | Decision of grant or rejection written | ||
| A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20080808 |
|
| A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20080814 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110829 Year of fee payment: 3 |
|
| R150 | Certificate of patent or registration of utility model |
Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110829 Year of fee payment: 3 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120829 Year of fee payment: 4 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130829 Year of fee payment: 5 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| LAPS | Cancellation because of no payment of annual fees |