JPWO2003051034A1 - 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 PDF

Info

Publication number
JPWO2003051034A1
JPWO2003051034A1 JP2002586705A JP2002586705A JPWO2003051034A1 JP WO2003051034 A1 JPWO2003051034 A1 JP WO2003051034A1 JP 2002586705 A JP2002586705 A JP 2002586705A JP 2002586705 A JP2002586705 A JP 2002586705A JP WO2003051034 A1 JPWO2003051034 A1 JP WO2003051034A1
Authority
JP
Japan
Prior art keywords
watermarked
point cloud
cloud data
data
fourier 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.)
Granted
Application number
JP2002586705A
Other languages
Japanese (ja)
Other versions
JP3935078B2 (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.)
Kokusai Kogyo Co Ltd
Original Assignee
Kokusai Kogyo 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 Kokusai Kogyo Co Ltd filed Critical Kokusai Kogyo Co Ltd
Publication of JPWO2003051034A1 publication Critical patent/JPWO2003051034A1/en
Application granted granted Critical
Publication of JP3935078B2 publication Critical patent/JP3935078B2/en
Anticipated expiration legal-status Critical
Expired - Fee Related legal-status Critical Current

Links

Images

Classifications

    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04N—PICTORIAL COMMUNICATION, e.g. TELEVISION
    • H04N1/00—Scanning, transmission or reproduction of documents or the like, e.g. facsimile transmission; Details thereof
    • H04N1/32—Circuits or arrangements for control or supervision between transmitter and receiver or between image input and image output device, e.g. between a still-image camera and its memory or between a still-image camera and a printer device
    • H04N1/32101—Display, printing, storage or transmission of additional information, e.g. ID code, date and time or title
    • H04N1/32144—Display, printing, storage or transmission of additional information, e.g. ID code, date and time or title embedded in the image data, i.e. enclosed or integrated in the image, e.g. watermark, super-imposed logo or stamp
    • H04N1/32149—Methods relating to embedding, encoding, decoding, detection or retrieval operations
    • H04N1/32154—Transform domain methods
    • H04N1/32187—Transform domain methods with selective or adaptive application of the additional information, e.g. in selected frequency coefficients
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01C—MEASURING DISTANCES, LEVELS OR BEARINGS; SURVEYING; NAVIGATION; GYROSCOPIC INSTRUMENTS; PHOTOGRAMMETRY OR VIDEOGRAMMETRY
    • G01C11/00—Photogrammetry or videogrammetry, e.g. stereogrammetry; Photographic surveying
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04N—PICTORIAL COMMUNICATION, e.g. TELEVISION
    • H04N1/00—Scanning, transmission or reproduction of documents or the like, e.g. facsimile transmission; Details thereof
    • H04N1/32—Circuits or arrangements for control or supervision between transmitter and receiver or between image input and image output device, e.g. between a still-image camera and its memory or between a still-image camera and a printer device
    • H04N1/32101—Display, printing, storage or transmission of additional information, e.g. ID code, date and time or title
    • H04N1/32144—Display, printing, storage or transmission of additional information, e.g. ID code, date and time or title embedded in the image data, i.e. enclosed or integrated in the image, e.g. watermark, super-imposed logo or stamp
    • H04N1/32149—Methods relating to embedding, encoding, decoding, detection or retrieval operations
    • H04N1/32154—Transform domain methods
    • H04N1/3216—Transform domain methods using Fourier transforms
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04N—PICTORIAL COMMUNICATION, e.g. TELEVISION
    • H04N2201/00—Indexing scheme relating to scanning, transmission or reproduction of documents or the like, and to details thereof
    • H04N2201/32—Circuits or arrangements for control or supervision between transmitter and receiver or between image input and image output device, e.g. between a still-image camera and its memory or between a still-image camera and a printer device
    • H04N2201/3201—Display, printing, storage or transmission of additional information, e.g. ID code, date and time or title
    • H04N2201/3225—Display, printing, storage or transmission of additional information, e.g. ID code, date and time or title of data relating to an image, a page or a document
    • H04N2201/3233—Display, printing, storage or transmission of additional information, e.g. ID code, date and time or title of data relating to an image, a page or a document of authentication information, e.g. digital signature, watermark
    • H—ELECTRICITY
    • H04—ELECTRIC COMMUNICATION TECHNIQUE
    • H04N—PICTORIAL COMMUNICATION, e.g. TELEVISION
    • H04N2201/00—Indexing scheme relating to scanning, transmission or reproduction of documents or the like, and to details thereof
    • H04N2201/32—Circuits or arrangements for control or supervision between transmitter and receiver or between image input and image output device, e.g. between a still-image camera and its memory or between a still-image camera and a printer device
    • H04N2201/3201—Display, printing, storage or transmission of additional information, e.g. ID code, date and time or title
    • H04N2201/3269—Display, printing, storage or transmission of additional information, e.g. ID code, date and time or title of machine readable codes or marks, e.g. bar codes or glyphs
    • H04N2201/327—Display, printing, storage or transmission of additional information, e.g. ID code, date and time or title of machine readable codes or marks, e.g. bar codes or glyphs which are undetectable to the naked eye, e.g. embedded codes

Abstract

３次元計測に応じて得られた点群データに電子透かしデータとして埋め込む際、点群データが規定されたｘ−ｙ面領域を予め規定された小領域に分割して、小領域毎に含まれる点群を生成して、点群毎にそのｘ，ｙ座標値を点群の重心が原点となるようにオフセットしてオフセット点群に変換する。オフセット点群毎に離散フーリエ変換を行ってフーリエ係数列を求め、電子透かしデータに応じてフーリエ係数列を変更して透かし入りフーリエ係数列とする。透かし入りフーリエ係数列を逆離散フーリエ変換して透かし入り複素数列を生成し、透かし入り複素数列について埋め込みに起因する座標値誤差のトレランスを満たす最適透かし埋め込み強度を求める。最適埋め込み強度に基づいてフーリエ係数列を再度変更して透かし入りフーリエ係数列とする。そして、透かし入りフーリエ係数列を逆オフセットして透かし入り点群データとする。When embedding as digital watermark data in the point cloud data obtained according to the three-dimensional measurement, the xy plane area in which the point cloud data is defined is divided into predetermined small areas and included in each small area. A point cloud is generated, and the x and y coordinate values are offset for each point cloud so that the center of gravity of the point cloud is the origin, and converted to an offset point cloud. A discrete Fourier transform is performed for each offset point group to obtain a Fourier coefficient sequence, and the Fourier coefficient sequence is changed according to the digital watermark data to obtain a watermarked Fourier coefficient sequence. The watermarked Fourier coefficient sequence is subjected to inverse discrete Fourier transform to generate a watermarked complex number sequence, and an optimum watermark embedding strength satisfying the tolerance of the coordinate value error caused by embedding is obtained for the watermarked complex number sequence. Based on the optimum embedding strength, the Fourier coefficient sequence is changed again to obtain a watermarked Fourier coefficient sequence. Then, the watermarked Fourier coefficient sequence is inversely offset to obtain watermarked point cloud data.

Description

技術分野
本発明は、３次元計測に基づいて得られたオリジナルランダム点データ群（３次元点群データ）を情報管理するための方法及びその装置に関し、特に、レーザ３次元計測によって地球表面を計測して得られた３次元点群データの不正使用を防止するために、電子情報を電子透かしデータとして埋め込む際に用いられる方法及び装置に関する。
背景技術
一般に、偽造又は不正使用を防止するため、印刷物上に多数の微細要素からなるパターンを埋め込んで、そのパターンで形成される情報を利用して、不正使用の検出を行うことが行われている。さらに、ベクトル形式（ベクター型）で記述された地図データに対する不正使用を防止するため、地図データに不正利用防止用の情報（埋め込み情報）を埋め込むことが行われており、ベクター型で記述された２次元又は３次元の地図データにおいては、面を構成するポリゴンに対して埋め込み情報を埋め込んでいる。例えば、地図データを三角形ポリゴンの集合で構成して、これら三角形ポリゴンをさらに４つの三角形に分割し、その間に構成される三角形（どの三角形ポリゴンの頂点も含まない三角形）に電子透かしデータを埋め込みデータとして埋め込む方法が知られている。そして、この方法を用いると、ベクター型で記述された地図データに影響を与えることなく、電子透かしデータを取り除くことが困難となる。
一方、特開２００１−１６０８９７公報には、ベクター型で記述された地図データに不正利用防止用の埋め込み情報を埋め込むことが記載されている。ここでは、オブジェクトの構成点の座標列を示す地図図形情報及びオブジェクトの種別を管理するレイヤ情報に、埋め込み基準レイヤと埋め込み参照レイヤを設定したものを入力して、同一の意味のある領域に存在するオブジェクトペアであって、それらが構成する領域内に他のオブジェクトが無いものを埋め込み基準オブジェクトペアとして選別する。そして、埋め込み基準オブジェクトペアに対して、埋め込み箇所を発見されにくいものを選択して、その存在位置及び／又は形状特徴に基づいて新たに埋め込むオブジェクト情報を算出し、このオブジェクト情報に合わせて埋め込み密度に応じて新たなオブジェクトを埋め込んでいる。
ところで、近年、所謂レーザ３次元計測によって地球表面（地表面）を計測して、地表面のデータをオリジナルランダム点データ群として得て、これらオリジナルランダム点群データから地図データを得ることが行われている。例えば、航空機から地上に向けてパルスレーザビームを照射して、地表面の空間座標を得る。この際、航空機の空間位置は、地上に設けられたＧＰＳ基準局と航空機に搭載されたＧＰＳ受信機によって求め、航空機の姿勢は三軸ジャイロスコープによって求める。
なお、１パルス毎の地表座標は、前述のようにして求めた航空機の空間位置及び姿勢に応じて、レーザミラー射角度及びその斜距離によってレーザビーム反射点のｘ，ｙ，及びｚ座標として求められることになる。
上述のようにして得られた地表座標は単に地表面のランダム点群データであるので、このランダム点群データを処理して地図データを生成するようにしている。
上述のようなオリジナル点群データは、単に空間的に離散している点群データであるので、相互に関係を持っておらず、しかも属性も備えていない。つまり、オリジナル点群データは単にｘ，ｙ，及びｚの座標値を示しているだけである。従って、上述のようなベクター型地図データに対して電子透かしデータを埋め込む手法を用いることができず、不正にオリジナルランダム点群データを用いて地図データを生成されてしまうことになる。
しかも、精度誤差の範囲で、オリジナルランダム点群データに対して乱数に応じて電子透かしデータを埋め込んだ場合には、オリジナルランダム点群データを再現することが困難となってしまう。また、オリジナルランダム点群データを部分的に間引いて、電子透かしデータを埋め込んだ場合には、不正使用を防止することが難しいという課題がある。
発明の開示
本発明の目的は、オリジナルランダム点群データの不正使用を防止することのできる電情報埋め込み方法、電子情報の埋め込みに用いられる装置、及びプログラムを提供することにある。
本発明の他の目的は、電子情報が埋め込まれたオリジナルランダム点群データから電子情報を抽出するための電子情報抽出方法、電子情報の抽出に用いられる装置、及びプログラムを提供することにある。
本発明によれば、レーザ３次元計測に応じて得られたオリジナルランダム点群データに電子情報を電子透かしデータとして埋め込む際に用いられる電子情報埋め込み方法であって、前記オリジナルランダム点群データについて離散フーリエ変換を行ってフーリエ係数列を求める第１のステップと、前記電子透かしデータに応じて前記フーリエ係数列を変更して透かし入りフーリエ係数列とする第２のステップと、前記透かし入りフーリエ係数列を逆離散フーリエ変換して該逆離散フーリエ変換に応じて透かし入り点群データを得る第３のステップとを有することを特徴とする電子情報埋め込み方法が得られる。
ここでは、前記第１のステップは、前記オリジナルランダム点群データが規定されたｘ−ｙ面領域を予め規定された小領域に分割して該小領域毎に含まれる点群を生成する第４のステップと、前記点群毎にそのｘ，ｙ座標値を前記点群の重心が原点となるようにオフセットしてオフセット点群に変換する第５のステップと、前記オフセット点群毎に前記離散フーリエ変換を行ってフーリエ係数列を求める第６のステップとを有している。
また、前記第３のステップは、前記透かし入りフーリエ係数列を逆離散フーリエ変換して透かし入り複素数列を生成する第７のステップと、該透かし入り複素数列について埋め込みに起因する座標値誤差のトレランスを満たす最適透かし埋め込み強度を求める第８のステップと、前記最適埋め込み強度に基づいて前記フーリエ係数列を再度変更して透かし入りフーリエ係数列とする第９のステップと、該透かし入りフーリエ係数列を逆オフセットして前記透かし入り点群データとする第１０のステップとを有している。
さらに、本発明では、前述のようにして得られた透かし入り点群データから前記電子透かしデータを抽出する際、前記オリジナルランダム点群データと前記透かし入り点群データとを対応付けてそれぞれ離散型フーリエ変換を行って第１及び第２のフーリエ係数列を求める第１１のステップと、前記第１及び前記第２のフーリエ係数列を比較して前記電子透かしデータを抽出する第１２のステップとを有することを特徴とする電子情報抽出方法が得られる。
ここでは、前記第１１のステップは、前記オリジナルランダム点群データが規定されたｘ−ｙ面領域を予め規定された小領域に分割して該小領域毎に含まれる小領域点群を生成する第１３のステップと、前記小領域点群毎に該小領域点群の頂点と最短距離にある頂点を最短距離頂点として前記透かし入り点群データから探索して前記オリジナルランダム点群データと前記透かし入り点群データとの対応付けを行う第１４のステップとを有している。
例えば、前記第１４のステップは、前記透かし入り点群データに対する２−ｄ木を生成する第１５のステップと、前記小領域点群の頂点位置と埋め込みトレランスとによって規定される問い合わせ領域を設定して前記２−ｄ木内から前記問い合わせ領域に含まれる前記透かし入り点群データを探索して前記オリジナルランダム点群データと前記透かし入り点群データとの対応付けを行う第１６のステップとを有している。
また、本発明によれば、レーザ３次元計測に応じて得られたオリジナルランダム点群データに電子情報を電子透かしデータとして埋め込む際に用いられる電子情報埋め込み装置であって、前記オリジナルランダム点群データについて離散フーリエ変換を行ってフーリエ係数列を求める離散フーリエ変換手段と、前記電子透かしデータに応じて前記フーリエ係数列を変更して透かし入りフーリエ係数列とする変更手段と、前記透かし入りフーリエ係数列を逆離散フーリエ変換して該逆離散フーリエ変換に応じて透かし入り点群データを得る透かし入り点群データ生成手段とを有することを特徴とする電子情報埋め込み装置が得られる。
ここでは、前記離散フーリエ変換手段は、前記オリジナルランダム点群データが規定されたｘ−ｙ面領域を予め規定された小領域に分割して該小領域毎に含まれる点群を生成する分割手段と、前記点群毎にそのｘ，ｙ座標値を前記点群の重心が原点となるようにオフセットしてオフセット点群に変換するオフセット手段と、前記オフセット点群毎に前記離散フーリエ変換を行ってフーリエ係数列を求めるフーリエ係数生成手段とを有している。
さらに、前記透かし入り点群データ生成手段は、前記透かし入りフーリエ係数列を逆離散フーリエ変換して透かし入り複素数列を生成する複素数列生成手段と、該透かし入り複素数列について埋め込みに起因する座標値誤差のトレランスを満たす最適透かし埋め込み強度を求める透かし埋め込み強度算出手段と、前記最適埋め込み強度に基づいて前記フーリエ係数列を再度変更して透かし入りフーリエ係数列とする再変更手段と、該透かし入りフーリエ係数列を逆オフセットして前記透かし入り点群データとする逆オフセット手段とを有している。
さらに、本発明によれば、上述のようにして得られた透かし入り点群データから前記電子透かしデータを抽出する際、前記オリジナルランダム点群データと前記透かし入り点群データとを対応付けてそれぞれ離散型フーリエ変換を行って第１及び第２のフーリエ係数列を求めるフーリエ係数生成手段と、前記第１及び前記第２のフーリエ係数列を比較して前記電子透かしデータを抽出する抽出手段とを有することを特徴とする電子情報抽出装置が得られる。
ここでは、前記フーリエ係数生成手段は、前記オリジナルランダム点群データが規定されたｘ−ｙ面領域を予め規定された小領域に分割して該小領域毎に含まれる小領域点群を生成する分割手段と、前記小領域点群毎に該小領域点群の頂点と最短距離にある頂点を最短距離頂点として前記透かし入り点群データから探索して前記オリジナルランダム点群データと前記透かし入り点群データとの対応付けを行う対応付け手段とを有している。
例えば、前記対応付け手段は、前記透かし入り点群データに対する２−ｄ木を生成する２−ｄ木生成手段と、前記小領域点群の頂点位置と埋め込みトレランスとによって規定される問い合わせ領域を設定して前記２−ｄ木内から前記問い合わせ領域に含まれる前記透かし入り点群データを探索して前記オリジナルランダム点群データと前記透かし入り点群データとの対応付けを行う頂点間対応付け手段とを有している。
加えて、本発明によれば、レーザ３次元計測に応じて得られたオリジナルランダム点群データに電子情報を電子透かしデータとして埋め込む際に電子計算機で用いられる電子情報埋め込みプログラムであって、前記オリジナルランダム点群データについて離散フーリエ変換を行ってフーリエ係数列を求める第１の手続きと、前記電子透かしデータに応じて前記フーリエ係数列を変更して透かし入りフーリエ係数列とする第２の手続きと、前記透かし入りフーリエ係数列を逆離散フーリエ変換して該逆離散フーリエ変換に応じて透かし入り点群データを得る第３の手続きとを有することを特徴とする電子情報埋め込みプログラムが得られる。
ここでは、前記第１の手続きは、前記オリジナルランダム点群データが規定されたｘ−ｙ面領域を予め規定された小領域に分割して該小領域毎に含まれる点群を生成する第４の手続きと、前記点群毎にそのｘ，ｙ座標値を前記点群の重心が原点となるようにオフセットしてオフセット点群に変換する第５の手続きと、前記オフセット点群毎に前記離散フーリエ変換を行ってフーリエ係数列を求める第６の手続きとを有している。
また、前記第３の手続きは、前記透かし入りフーリエ係数列を逆離散フーリエ変換して透かし入り複素数列を生成する第７の手続きと、該透かし入り複素数列について埋め込みに起因する座標値誤差のトレランスを満たす最適透かし埋め込み強度を求める第８の手続きと、前記最適埋め込み強度に基づいて前記フーリエ係数列を再度変更して透かし入りフーリエ係数列とする第９の手続きと、該透かし入りフーリエ係数列を逆オフセットして前記透かし入り点群データとする第１０の手続きとを有している。
さらに、本発明によれば、前述のようにして得られた透かし入り点群データから前記電子透かしデータを抽出する際、前記オリジナルランダム点群データと前記透かし入り点群データとを対応付けてそれぞれ離散型フーリエ変換を行って第１及び第２のフーリエ係数列を求める第１１の手続きと、前記第１及び前記第２のフーリエ係数列を比較して前記電子透かしデータを抽出する第１２の手続きとを有することを特徴とする電子情報抽出プログラムが得られる。
ここでは、前記第１１の手続きは、前記オリジナルランダム点群データが規定されたｘ−ｙ面領域を予め規定された小領域に分割して該小領域毎に含まれる小領域点群を生成する第１３の手続きと、前記小領域点群毎に該小領域点群の頂点と最短距離にある頂点を最短距離頂点として前記透かし入り点群データから探索して前記オリジナルランダム点群データと前記透かし入り点群データとの対応付けを行う第１４の手続きとを有している。
例えば、前記第１４の手続きは、前記透かし入り点群データに対する２−ｄ木を生成する第１５の手続きと、前記小領域点群の頂点位置と埋め込みトレランスとによって規定される問い合わせ領域を設定して前記２−ｄ木内から前記問い合わせ領域に含まれる前記透かし入り点群データを探索して前記オリジナルランダム点群データと前記透かし入り点群データとの対応付けを行う第１６の手続きとを有している。
発明を実施するための最良の形態
以下本発明について図面を参照して説明する。
第１図を参照して、地球表面（例えば、地表面）に関するオリジナルランダム点群データを得る際には、まず、地球表面の３次元計測が行われる。地球表面の３次元計測に当たっては、例えば、レーザ３次元計測が用いられる。３次元レーザ計測を行う際には、航空機１１をオリジナルランダム点群データを収集すべき地域（地区）上を飛行させる。航空機１１にはレーザスキャナ（航空機搭載型レーザスキャナ：図示せず）、ＧＰＳ受信機（図示せず）、及びジャイロスコープ（ＩＭＵ：図示せず）が搭載されている。
第２図も参照して、レーザ３次元計測（ステップＳ１）では、航空機搭載型レーザスキャナから地上（地表面）１２に向けてパルスレーザビームを照射して、地表面１２からの反射パルスをパルスリターンとして受け、パルスレーザビームの照射時刻とパルスリターンを受けた時刻との時間差に応じて地表面１２までの距離（ｚ）を計測して距離データを得る。このパルスレーザビームの照射は予め規定された時間間隔で行われる結果、地表面は離散的な点として規定されることになる。この際、レーザスキャナのスキャンミラーの回転は１軸（航空機１１の飛行方向に直交する）であり、飛行方向がもう一つの走査軸となって、面的にデータが取得されることになる。そして、ＩＭＵは１パルス毎にレーザ照射角度をミラー回転角及びその傾き角度に基づいて求める。
さらに、別に地上において、レーザスキャナ（地上据え置き型レーザスキャナ）を用いて地表面のデータを取得する。このレーザスキャナではスキャンミラーの回転軸は２軸で放射状にしかも面的にデータが得られることになる。
一方、ＧＰＳ電波に応じて基準点測量が行われる（ステップＳ２）。つまり、航空機１１の空間位置は、地上に設けられたＧＰＳ基準局１３、ＧＰＳ受信機、及びＩＭＵによって計測される。その後、レーザ３次元計測によって得られたデータ（スキャンポイントデータ）と基準点測量とによって得られたデータとに基づいてスキャンポイントデータを測地座標データへ変換してオリジナルランダム点群データを得る（ステップＳ３）。つまり、１パルス毎に地表面１２におけるレーザビーム反射点のｘ，ｙ，ｚ（高さ：Ｈ）座標を求めて、オリジナルランダム点群データとする。
そして、上述のオリジナルランダム点群データに電子透かしデータを埋め込み（ステップＳ４）、透かし情報（データ）埋め込み済み３次元空間座標点群として出力する（ステップＳ５）。なお、ステップＳ３〜Ｓ５はコンピュータ（図示せず）によって行われる。つまり、コンピュータはステップＳ３を行う測地座標変換部、ステップＳ４を行う電子透かし埋め込み部、及びステップＳ５を行う透かし入り点群出力部を少なくとも有しており、測地座標変換部では、スキャンポイントデータと基準点測量データとを読み込んでオリジナルランダム点群データを生成する。
ここで、第３図を参照して、電子透かし埋め込み部の動作について説明する。
図示の例では、電子透かし埋め込み部は、小領域分割部、座標値オフセット部、離散フーリエ変換部、フーリエ係数列変更部、逆離散フーリエ変換部、埋め込み強度自動調整部、及び逆オフセット部を有している。
電子透かしデータ埋め込み処理の際には、前述の測地座標変換部からオリジナルランダム点群データ（以下単にオリジナル点群と呼ぶ）Ｖが小領域分割部に与えられるとともに別に入力装置（図示せず）等から小領域ｘ，ｙ方向サイズＤｘ，Ｄｙが小領域分割部に与えられる。
ここでは、オリジナル点群Ｖは

Figure 2003051034
小領域分割部では、まず、オリジナル点群の小領域への分割を行う（ステップＡ１）。第４図に示すように、オリジナル点群Ｖが定義されているｘ−ｙ面内の矩形領域を、Ｄｘ，Ｄｙのサイズ（寸法）を有する小領域Ａａに均等に分割して（等分割）、小領域Ａａ毎に含まれる点群Ｖａを生成する。そして、この点群Ｖａは座標オフセット部に与えられる。
Figure 2003051034
なお、第４図において、その下辺及び右辺のように、ｘ方向及びｙ方向へ分割した際、その最終端では、Ｄｘ，Ｄｙよりも小さいサイズの小領域（この領域に含まれる点群をＶＬａとする）が形成されることもある。
座標オフセット部では、小領域内の点群のｘ，ｙ座標値のオフセットを行う（ステップＡ２）。ここでは、各小領域の点群Ｖａのｘ，ｙ座標値を、Ｖａ
Figure 2003051034
上述のオフセット処理によって、オリジナル点群Ｖのｘ，ｙ座標値の絶対値が大きい場合に、埋め込み処理で用いられる浮動小数点演算における丸め誤差の影響を減少させることができる。
Figure 2003051034
部に与えられ、一方、小領域重心点μａｘ，μａｙは逆オフセット部に与えられる。
Figure 2003051034
れる点のｘ，ｙ座標値から、式（４）で示す複素数列｛ｃｋ｝を生成する。
Figure 2003051034
次に、複素数列｛ｃｋ｝を式（５）に基づいて、離散フーリエ変換部は、離散フーリエ変換して，フーリエ係数列｛Ｃｌ｝を求める。
Figure 2003051034
そして、このフーリエ係数列｛Ｃｌ｝はフーリエ係数列変更部に与えられる。フーリエ係数列変更部には電子透かしデータビット列Ｂ及び埋め込み強度初期値Ｐｉｎｉｔが与えられており、フーリエ係数列変更部では、電子透かしデータビット列Ｂ及び埋め込み強度初期値Ｐｉｎｉｔに応じてフーリエ係数列｛Ｃｌ｝を変更する（ステップＡ４）。
ここで、電子透かしデータビット列Ｂは
Figure 2003051034
つまり、電子透かしデータビット列Ｂ＝｛ｂｍ∈｛０，１｝｜１≦ｍ≦Ｎｂ｝及び埋め込み強度初期値ｐｉｎｉｔを用いて、フーリエ係数列変更部は、式（６）に従ってフーリエ係数列｛Ｃｌ｝を変更して、変更後フーリエ係数列｛Ｃ’ｌ｝を生成する。
Figure 2003051034
なお、直流成分を表す係数Ｃ０は透かしデータによって変更しない。従って、この小領域Ａａ内に埋め込み可能な透かしビット長さは、｜Ｖａ｜−１（ビット）となる。また、埋め込み強度ｐの値は、１回目の埋め込みの際には、ｐ←ｐｉｎｉｔと定め、２回目の埋め込みの際には、ｐ←ｐｏｐｔと定める。
この際、透かしビット長Ｎｂと小領域への埋め込み可能ビット長｜Ｖａ｜−１との関係に基づいて次の規則を適用される。
ｉ）Ｎｂ＜｜Ｖａ｜−１の場合
透かしビットを反復して埋め込み，Ｃ１，Ｃ２，．．．．，Ｃ｜Ｖａ｜−１のフーリエ係数全ての値を前述の式（６）で変更する。
ｉ）Ｎｂ＞｜Ｖａ｜−１の場合
透かしビットの先頭から｜Ｖａ｜−１ビット分のみを埋め込み，それ以降のビットについては埋め込みを行わない。
上述の変更後フーリエ係数列｛Ｃ’ｌ｝、つまり、透かし入りフーリエ係数列｛Ｃ’ｌ｝は逆離散フーリエ変換部に与えられ、逆離散フーリエ変換部では、変更後フーリエ係数列｛Ｃ’ｌ｝に対して逆離散フーリエ変換（ＩＤＦＴ）を行う（ステップＡ５）。つまり、逆離散フーリエ変換部では、透かし入りフーリエ係数列｛Ｃ’ｌ｝に対して、式（７）の逆フーリエ変換を行って、電子透かしによって変更された複素数列｛ｃ’ｋ｝を生成する。
Figure 2003051034
この複素数列｛ｃ’ｋ｝は、埋め込み強度自動調整部に与えられる。埋め
Figure 2003051034
τが与えられており、埋め込み強度自動調整部では、埋め込みトレランスτを考慮して、埋め込み強度を自動調整する（ステップＡ６）。
Figure 2003051034
つまり、埋め込み強度自動調整部では、変更後複素数列｛ｃ’ｋ｝に基づ
Figure 2003051034
Figure 2003051034
となる。
ここで，埋め込みによって生ずるｘ，ｙ方向の頂点座標値誤差のトレランス（許容値）をτ，このトレランスを満たす最適な埋め込み強度をｐｏｐｔとすると、式（９）に示す比例関係が成り立つ．
Figure 2003051034
従って，式（９）から、埋め込み誤差のトレランスτを満たす最適な透かし埋め込み強度ｐｏｐｔは、
Figure 2003051034
となる。そして、この最適な透かし埋め込み強度ｐｏｐｔは、フーリエ係数列変更部にフィードバックされる。
フーリエ係数列変更部では、最適埋め込み強度ｐｏｐｔを用いてフーリエ係数を変更し、逆フーリエ変換部ではその変更結果に応じて逆フーリエ変換を再度実行する。つまり、前述のようにして求めた最適埋め込み強度ｐｏｐｔを用いて、前述のステップＡ４及びＡ５を再度実行して、電子透かし入りオフ
Figure 2003051034
ように、逆オフセット部には、小領域重心点μａｘ，μａｙが与えられており、逆オフセット部は、ｘ，ｙ座標値の逆オフセットを行って透かし入り点群データＶ’ａを生成する（ステップＡ７）。ここでは、前述の式（２）の逆操作を行って、つまり、逆オフセットを式（１１）に基づいて行って、透かし入
Figure 2003051034
上述のようにして、電子透かし埋め込み部では、各小領域毎に処理を行い、最終的に透かし入り点群Ｖ’を生成する。この透かし入り点群は、透かし入り点群出力部によって透かし情報埋め込み済み空間３次元座標点群としてファイルに出力される。
第５図を参照して、透かし入り点群データＶ’から透かしデータを抽出する場合について説明する。ここでは、コンピュータは透かしデータ抽出部を有しており、透かしデータ抽出部には点群小領域分割部、頂点間対応付け部、離散フーリエ変換部、及び透かしビット列抽出部が備えられている。
いま、オリジナル点群Ｖを、
Figure 2003051034
この際、前述のＤｘ，Ｄｙ（分割小領域のｘ，ｙ方向サイズ）及びτ（ｘ，ｙ方向の埋め込みトレランス）は電子透かしデータ埋め込みの際と同一の値を用いる。
まず、オリジナル点群Ｖ及び小領域ｘ，ｙ方向サイズＤｘ，Ｄｙが点群小領域分割部に与えられ、点群小領域分割部はオリジナル点群Ｖを小領域に分割する（ステップＡ８）。
前述の透かし埋め込み時に指定した分割小領域のｘ，ｙ方向サイズＤｘ，Ｄｙの値を用いて、埋め込み時と同様にしてオリジナル点群Ｖを小領域に分割して、領域毎に点群Ｖａを求める。
なお、透かし入り点群Ｖ’は分割を行うことなく、そのまま保存される。
小領域内点群Ｖは頂点間対応付け部に与えられる。この頂点間対応付け部には、前述の埋め込みトレランスτ及び透かし入り点群Ｖ’が与えられており、頂点間対応付け部では、小領域内点群Ｖと透かし入り点群Ｖ’との頂点間の対応づけを行って小領域内透かし入り点群Ｖ’ａを出力する（ステップＡ９）。
例えば、頂点間対応付け部では、各小領域群Ｖａ毎に点群Ｖａ内の頂点ｖｉ（∈Ｖａ）と最短距離にある頂点を透かし入り点群Ｖ’から探索して、この探索の結果得られた頂点を、頂点ｖｉに対応付けられた透かし入り頂点ｖ’ｉ（∈Ｖ’）として採用する。
上述の対応付けにおいて、点群Ｖ’は極めて多数の頂点を含む関係上、総当たりで最短頂点を探索すると、点数の二乗のオーダーで探索時間が必要となる。このため、図示のように、２−ｄ木生成部を備えて、前処理として透かし入り点群Ｖ’に対する２−ｄ木を生成するようにしてもよい（ステップＡ１０）。そして、透かし入り点群Ｖ’の代わりに、透かし入り点群の２−ｄ木を頂点間対応付け部に与えるようにしてもよい。
そして、オリジナル点群内の頂点ｖｉ（∈Ｖａ）の位置と埋め込みトレランスτとによって定まる微小な問い合わせ領域を設定して、２−ｄ木内から問い合わせ領域に含まれる透かし入り頂点群を高速に探索して、これらの中だけで総当たりによって最短頂点を探索する。これによって、オリジナル点群Ｖｉ（∈Ｖａ）と透かし入り頂点ｖ’ｉ（∈Ｖ’）との対応付け処理を効率化することができる。
このようにして、小領域に含まれるオリジナル点群Ｖａに透かし入り点群Ｖ’ａを対応付けた後、離散フーリエ変換部では、オリジナル点群Ｖａ及び透かし入り点群Ｖ’ａに対してＤＦＴを実行する（ステップＡ１１）。オリジナル点群及びこれに対応づけられた透かし入り点群に、埋め込み時と同様にして、式（４）及び（５）で示すＤＦＴを行って、それぞれフーリエ係数列｛Ｃｌ｝及び｛Ｃ’ｌ｝を求める。
そして、透かしビット列抽出部では、フーリエ係数列｛Ｃｌ｝及び透かし入りフーリエ係数列｛Ｃ’ｌ｝に応じて電子透かしデータの抽出を行う（ステップＡ１２）。つまり、上述のようにして得られたフーリエ係数列｛Ｃｌ｝及び｛Ｃ’ｌ｝を比較して、埋め込まれた透かしデータのビット列
Figure 2003051034
式（１２）に基づいて抽出して、別に管理していたオリジナルの透かしビット列と比較する。この比較結果に基づいて透かし入り点群データがオリジナル点群データから生成されたものであるかを判定する。
産業上の利用可能性
このようにして、本発明では、オリジナルランダム点群データに電子透かしデータを埋め込んで、オリジナルランダム点群データの不正使用を防止することができる。さらに、上述のようにして、電子透かしデータを、透かし入り点群データから抽出して、別に管理しているオリジナルの電子透かしデータと比較するようにしたから、その結果に応じて透かし入り点群データがオリジナル点群データから生成されたものであるかを判定することができる。つまり、オリジナルランダム点群データの同一性を判定することもできる。
【図面の簡単な説明】
第１図は３次元レーザ計測による地表面の計測を説明するための図である。
第２図は３次元レーザ計測によって得られたデータから透かし入り点群を得る際の手法を説明するためのフローチャートである。
第３図は電子透かし埋め込みについて説明するためのフローチャートである。
第４図はオリジナルランダム点群データの一例を示す図である。
第５図は透かしデータの抽出を説明するためのフローチャートである。Technical field
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.
Background art
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 (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.
Disclosure of the invention
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, there is provided an electronic information embedding method used when embedding electronic information as electronic watermark data in original random point cloud data obtained in accordance with laser three-dimensional measurement, wherein the original random point cloud data is discrete. A first step of performing a Fourier transform to obtain a Fourier coefficient sequence; a second step of changing the Fourier coefficient sequence according to the digital watermark data to form a watermarked Fourier coefficient sequence; and the watermarked Fourier coefficient sequence And a third step of obtaining watermarked point cloud data in accordance with the inverse discrete Fourier transform.
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.
According to the present invention, there is also provided an electronic information embedding device used when embedding electronic information as electronic watermark data in original random point cloud data obtained in accordance with laser three-dimensional measurement, the original random point cloud data. Discrete Fourier transform means for obtaining a Fourier coefficient sequence by performing discrete Fourier transform on, a changing means for changing the Fourier coefficient sequence according to the digital watermark data into a watermarked Fourier coefficient sequence, and the watermarked Fourier coefficient sequence There is obtained an electronic information embedding apparatus characterized by having a watermarked point cloud data generation means for obtaining a watermarked point cloud data in accordance with the inverse discrete Fourier transform.
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, there is provided an electronic information embedding program used in an electronic computer for embedding electronic information as electronic watermark data in original random point cloud data obtained according to laser three-dimensional measurement, A first procedure for performing a discrete Fourier transform on random point cloud data to obtain a Fourier coefficient sequence; a second procedure for changing the Fourier coefficient sequence according to the digital watermark data to form a watermarked Fourier coefficient sequence; And a third procedure for obtaining watermarked point cloud data according to the inverse discrete Fourier transform by performing the inverse discrete Fourier transform on the watermarked Fourier coefficient sequence.
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.
BEST MODE FOR CARRYING OUT THE INVENTION
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 S1), 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 pulsed. Received as a return, 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.
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. are 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. To small area x, y-direction size Dx, DyIs given to the small area dividing unit.
Here, the original point cloud V is
Figure 2003051034
In the small area dividing unit, first, the original point cloud is divided into small areas (step A1). As shown in FIG. 4, a rectangular area in the xy plane where the original point group V is defined is represented by Dx, DyPoint group V included in each small area Aa by equally dividing the same into small areas Aa having the same size (dimension) (equal division).aIs generated. And this point cloud VaIs given to the coordinate offset part.
Figure 2003051034
In FIG. 4, when divided in the x direction and the y direction as in the lower side and the right side, at the final end, Dx, DySmaller area (points contained in this area are represented by VLaMay be formed.
In the coordinate offset unit, the x and y coordinate values of the point group in the small area are offset (step A2). Here, the point group V of each small regionaX and y coordinate values of Va
Figure 2003051034
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.
Figure 2003051034
While the small area centroid point μax, ΜayIs given to the reverse offset part.
Figure 2003051034
Complex number sequence {ck} Is generated.
Figure 2003051034
Next, a complex sequence {ck} Based on equation (5), the discrete Fourier transform unit performs discrete Fourier transform to obtain a Fourier coefficient sequence {Cl}.
Figure 2003051034
And this Fourier coefficient sequence {Cl} Is given to the Fourier coefficient sequence changing unit. The Fourier coefficient sequence changing unit includes a digital watermark data bit sequence B and an embedding strength initial value P.initIn the Fourier coefficient sequence changing unit, the digital watermark data bit sequence B and the embedding strength initial value PinitFourier coefficient sequence {Cl} Is changed (step A4).
Here, the digital watermark data bit string B is
Figure 2003051034
That is, the digital watermark data bit string B = {bm∈ {0, 1} | 1 ≦ m ≦ Nb} And embedding strength initial value pinitThe Fourier coefficient sequence changing unit uses the Fourier coefficient sequence {C according to the equation (6).l} And the Fourier coefficient sequence {C ′ after the changel} Is generated.
Figure 2003051034
Note that the coefficient C representing the DC component0Is not changed by watermark data. Therefore, the watermark bit length that can be embedded in the small area Aa is | Va| -1 (bit). The embedding strength p is set to p ← p at the first embedding.initP ← p for the second embeddingoptIt is determined.
At this time, the watermark bit length NbAnd bit length that can be embedded in small area | VaThe following rule is applied based on the relationship with | -1.
i) Nb<| VaIn case of -1
Embed watermark bit repeatedly, C1, C2,. . . . , C| Va | -1The values of all the Fourier coefficients are changed by the above equation (6).
i) Nb> | VaIn case of -1
| V from the beginning of the watermark bita| Only one bit is embedded, and the subsequent bits are not embedded.
The modified Fourier coefficient sequence {C ′l}, That is, a watermarked Fourier coefficient sequence {C ′l} Is given to the inverse discrete Fourier transform unit, and in the inverse discrete Fourier transform unit, the post-change Fourier coefficient sequence {C ′l} Is subjected to inverse discrete Fourier transform (IDFT) (step A5). That is, in the inverse discrete Fourier transform unit, the watermarked Fourier coefficient sequence {C ′l} Is subjected to the inverse Fourier transform of equation (7), and the complex number sequence {c ′k} Is generated.
Figure 2003051034
This complex number sequence {c ’k} Is given to the embedding strength automatic adjustment unit. Padding
Figure 2003051034
τ is given, and the embedding strength automatic adjustment unit automatically adjusts the embedding strength in consideration of the embedding tolerance τ (step A6).
Figure 2003051034
That is, the embedding strength automatic adjustment unit changes the complex number sequence {c ′k} Based on
Figure 2003051034
Figure 2003051034
It becomes.
Here, the tolerance (allowable value) of the error in the vertex coordinate values in the x and y directions caused by embedding is τ, and the optimum embedding strength satisfying this tolerance is p.optThen, the proportional relationship shown in Equation (9) holds.
Figure 2003051034
Therefore, from equation (9), the optimum watermark embedding strength p satisfying the embedding error tolerance τoptIs
Figure 2003051034
It becomes. And this optimal watermark embedding strength poptIs fed back to the Fourier coefficient sequence changing unit.
In the Fourier coefficient sequence changing unit, the optimum embedding strength poptAnd the Fourier coefficient is changed, and the inverse Fourier transform unit executes the inverse Fourier transform again according to the change result. That is, the optimum embedding strength p obtained as described above.opt, The above steps A4 and A5 are executed again, and the digital watermark is turned off.
Figure 2003051034
In the reverse offset part, the small area centroid point μax, ΜayThe inverse offset unit performs the inverse offset of the x and y coordinate values and performs watermarked point cloud data V ′.aIs generated (step A7). Here, the reverse operation of the above-described equation (2) is performed, that is, the reverse offset is performed based on the equation (11), and watermarking is performed.
Figure 2003051034
As described above, the digital watermark embedding unit performs processing for each small area, 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 the 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
Figure 2003051034
At this time, the aforementioned Dx, DyThe same values as those used when embedding digital watermark data are used for (size in the x and y directions of the divided small areas) and τ (embedding tolerance in the x and y directions).
First, the original point group V and the small region x, y direction size Dx, DyIs 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).
X, y direction size D of the divided small area specified at the time of embedding the watermarkx, DyIs used to divide the original point cloud V into small regions in the same manner as in embedding, and the point cloud V for each region.aAsk for.
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 given the above-described embedding tolerance τ and the watermarked point group V ′. In the inter-vertex association unit, the vertices of the small-area in-point group V and the watermarked point group V ′ Point group V ′ with a small area watermarkaIs output (step A9).
For example, in the inter-vertex association unit, each small region group VaPoint cloud V every timeaVertex vi(∈Va) And the vertex at the shortest distance are searched from the watermarked point group V ′, and the vertex obtained as a result of this search is obtained as the vertex v.iWatermarked vertex v 'associated withiIt is adopted as (εV ′).
In the above association, since the point group V 'includes a very large number of vertices, searching for the shortest vertex in the round robin requires a search time in the order of the square of the number of points. For this reason, as shown in the drawing, 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 intervertex correspondence unit.
And the vertex v in the original point cloudi(∈Va) And a small inquiry area determined by the embedding tolerance τ, and a watermarked vertex group included in the inquiry area is searched from the 2-d tree at high speed. Explore. As a result, the original point cloud Vi(∈Va) And watermarked vertex v ’iThe association process with (∈V ′) can be made efficient.
In this way, the original point cloud V included in the small areaaWatermarked point cloud V ’aAre associated with each other, the discrete Fourier transform unit then converts the original point group VaAnd watermarked point cloud V 'aA DFT is executed for (step A11). Similar to the case of embedding, the original point cloud and the watermarked point cloud associated therewith are subjected to DFT represented by equations (4) and (5), respectively, and Fourier coefficient sequences {Cl} And {C 'l}.
Then, in the watermark bit string extraction unit, the Fourier coefficient string {Cl} And a watermarked Fourier coefficient sequence {C 'l}, Digital watermark data is extracted (step A12). That is, the Fourier coefficient sequence {Cl} And {C 'l} And a bit string of embedded watermark data
Figure 2003051034
Extraction is performed based on Expression (12), and the result is 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.
Industrial applicability
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.
[Brief description of the drawings]
FIG. 1 is a diagram for explaining ground surface measurement by three-dimensional laser measurement.
FIG. 2 is a flowchart for explaining a technique for obtaining a watermarked point group from data obtained by three-dimensional laser measurement.
FIG. 3 is a flowchart for explaining digital watermark embedding.
FIG. 4 is a diagram showing an example of original random point cloud data.
FIG. 5 is a flowchart for explaining the extraction of watermark data.

Claims (18)

３次元計測に応じて得られたオリジナルランダム点群データに電子情報を電子透かしデータとして埋め込む際に用いられる電子情報埋め込み方法であって、前記オリジナルランダム点群データについて離散フーリエ変換を行ってフーリエ係数列を求める第１のステップと、前記電子透かしデータに応じて前記フーリエ係数列を変更して透かし入りフーリエ係数列とする第２のステップと、前記透かし入りフーリエ係数列を逆離散フーリエ変換して該逆離散フーリエ変換に応じて透かし入り点群データを得る第３のステップとを有することを特徴とする電子情報埋め込み方法。An electronic information embedding method used when embedding electronic information as electronic watermark data in original random point cloud data obtained in accordance with three-dimensional measurement, wherein a Fourier coefficient is obtained by performing discrete Fourier transform on the original random point cloud data A first step of obtaining a sequence; a second step of changing the Fourier coefficient sequence according to the digital watermark data to form a watermarked Fourier coefficient sequence; and performing an inverse discrete Fourier transform on the watermarked Fourier coefficient sequence And a third step of obtaining watermarked point cloud data in accordance with the inverse discrete Fourier transform. 前記第１のステップは、前記オリジナルランダム点群データが規定されたｘ−ｙ面領域を予め規定された小領域に分割して該小領域毎に含まれる点群を生成する第４のステップと、前記点群毎にそのｘ，ｙ座標値を前記点群の重心が原点となるようにオフセットしてオフセット点群に変換する第５のステップと、前記オフセット点群毎に前記離散フーリエ変換を行ってフーリエ係数列を求める第６のステップとを有することを特徴とする請求の範囲第１項に記載の電子情報埋め込み方法。The first step is a fourth step of generating a point group included in each small region by dividing the xy plane region in which the original random point cloud data is defined into predetermined small regions. 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 the offset to the offset point group; and the discrete Fourier transform for each of the offset point groups. The electronic information embedding method according to claim 1, further comprising: a sixth step of performing a Fourier coefficient sequence. 前記第３のステップは、前記透かし入りフーリエ係数列を逆離散フーリエ変換して透かし入り複素数列を生成する第７のステップと、該透かし入り複素数列について埋め込みに起因する座標値誤差のトレランスを満たす最適透かし埋め込み強度を求める第８のステップと、前記最適埋め込み強度に基づいて前記フーリエ係数列を再度変更して透かし入りフーリエ係数列とする第９のステップと、該透かし入りフーリエ係数列を逆オフセットして前記透かし入り点群データとする第１０のステップとを有することを特徴とする請求の範囲第２項に記載の電子情報埋め込み方法。In the third step, a seventh step of generating a watermarked complex number sequence by performing inverse discrete Fourier transform on the watermarked Fourier coefficient sequence, and satisfying a tolerance of a coordinate value error caused by embedding of the watermarked complex number sequence An eighth step of obtaining an optimum watermark embedding strength, a ninth step of changing the Fourier coefficient sequence again based on the optimum embedding strength to obtain a watermarked Fourier coefficient sequence, and an inverse offset of the watermarked Fourier coefficient sequence The electronic information embedding method according to claim 2, further comprising: a tenth step of setting the watermarked point cloud data. 請求の範囲第１項に記載された電子埋め込み方法で得られた透かし入り点群データから前記電子透かしデータを抽出する抽出方法であって、前記オリジナルランダム点群データと前記透かし入り点群データとを対応付けてそれぞれ離散型フーリエ変換を行って第１及び第２のフーリエ係数列を求める第１１のステップと、前記第１及び前記第２のフーリエ係数列を比較して前記電子透かしデータを抽出する第１２のステップとを有することを特徴とする電子情報抽出方法。An extraction method for extracting the digital watermark data from the watermarked point cloud data obtained by the electronic embedding method according to claim 1, wherein the original random point cloud data, the watermarked point cloud data, The digital watermark data is extracted by comparing the first and second Fourier coefficient sequences with the eleventh step of obtaining the first and second Fourier coefficient sequences by performing discrete Fourier transforms in association with each other. An electronic information extraction method comprising: a twelfth step. 前記第１１のステップは、前記オリジナルランダム点群データが規定されたｘ−ｙ面領域を予め規定された小領域に分割して該小領域毎に含まれる小領域点群を生成する第１３のステップと、前記小領域点群毎に該小領域点群の頂点と最短距離にある頂点を最短距離頂点として前記透かし入り点群データから探索して前記オリジナルランダム点群データと前記透かし入り点群データとの対応付けを行う第１４のステップとを有することを特徴とする請求の範囲第４項に記載の電子情報抽出方法。In the eleventh step, a xy plane region in which the original random point cloud data is defined is divided into predetermined small regions, and a small region point group included in each small region is generated. Searching for the original random point cloud data and the watermarked point cloud by searching from the watermarked point cloud data with the shortest distance vertex being the vertex that is the shortest distance from the vertex of the small area point cloud for each small area point cloud; The electronic information extraction method according to claim 4, further comprising a fourteenth step of associating with data. 前記第１４のステップは、前記透かし入り点群データに対する２−ｄ木を生成する第１５のステップと、前記小領域点群の頂点位置と埋め込みトレランスとによって規定される問い合わせ領域を設定して前記２−ｄ木内から前記問い合わせ領域に含まれる前記透かし入り点群データを探索して前記オリジナルランダム点群データと前記透かし入り点群データとの対応付けを行う第１６のステップとを有することを特徴とする請求の範囲第５項に記載の電子情報抽出方法。In the fourteenth step, a fifteenth step of generating a 2-d tree for the watermarked point cloud data, an inquiry region defined by a vertex position of the small region point cloud and an embedding tolerance are set, and A sixteenth step of searching for the watermarked point cloud data included in the inquiry area from a 2-d tree and associating the original random point cloud data with the watermarked point cloud data. The electronic information extraction method according to claim 5. ３次元計測に応じて得られたオリジナルランダム点群データに電子情報を電子透かしデータとして埋め込む際に用いられる電子情報埋め込み装置であって、前記オリジナルランダム点群データについて離散フーリエ変換を行ってフーリエ係数列を求める離散フーリエ変換手段と、前記電子透かしデータに応じて前記フーリエ係数列を変更して透かし入りフーリエ係数列とする変更手段と、前記透かし入りフーリエ係数列を逆離散フーリエ変換して該逆離散フーリエ変換に応じて透かし入り点群データを得る透かし入り点群データ生成手段とを有することを特徴とする電子情報埋め込み装置。An electronic information embedding device used when embedding electronic information as electronic watermark data in original random point cloud data obtained according to three-dimensional measurement, and performing Fourier transform on the original random point cloud data to perform a Fourier coefficient Discrete Fourier transform means for obtaining a sequence, change means for changing the Fourier coefficient sequence in accordance with the digital watermark data into a watermarked Fourier coefficient sequence, and inverse discrete Fourier transform of the watermarked Fourier coefficient sequence An electronic information embedding apparatus comprising: watermarked point cloud data generating means for obtaining watermarked point cloud data in accordance with discrete Fourier transform. 前記離散フーリエ変換手段は、前記オリジナルランダム点群データが規定されたｘ−ｙ面領域を予め規定された小領域に分割して該小領域毎に含まれる点群を生成する分割手段と、前記点群毎にそのｘ，ｙ座標値を前記点群の重心が原点となるようにオフセットしてオフセット点群に変換するオフセット手段と、前記オフセット点群毎に前記離散フーリエ変換を行ってフーリエ係数列を求めるフーリエ係数生成手段とを有することを特徴とする請求の範囲第７項に記載の電子情報埋め込み装置。The discrete Fourier transform means divides the xy plane region in which the original random point cloud data is defined into predetermined small regions, and generates a point group included in each small region; Offset means for 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 the offset point group into an offset point group, and performing the discrete Fourier transform for each of the offset point groups to perform Fourier coefficients 8. The electronic information embedding apparatus according to claim 7, further comprising Fourier coefficient generation means for obtaining a sequence. 前記透かし入り点群データ生成手段は、前記透かし入りフーリエ係数列を逆離散フーリエ変換して透かし入り複素数列を生成する複素数列生成手段と、該透かし入り複素数列について埋め込みに起因する座標値誤差のトレランスを満たす最適透かし埋め込み強度を求める透かし埋め込み強度算出手段と、前記最適埋め込み強度に基づいて前記フーリエ係数列を再度変更して透かし入りフーリエ係数列とする再変更手段と、該透かし入りフーリエ係数列を逆オフセットして前記透かし入り点群データとする逆オフセット手段とを有することを特徴とする請求の範囲第８項に記載の電子情報埋め込み装置。The watermarked point cloud data generation means includes a complex number sequence generation means for generating a watermarked complex number sequence by performing an inverse discrete Fourier transform on the watermarked Fourier coefficient sequence, and a coordinate value error caused by embedding the watermarked complex number sequence. Watermark embedding strength calculating means for obtaining an optimum watermark embedding strength satisfying tolerance, re-changing means for changing the Fourier coefficient sequence again based on the optimum embedding strength into a watermarked Fourier coefficient sequence, and the watermarked Fourier coefficient sequence 9. The electronic information embedding apparatus according to claim 8, further comprising reverse offset means for reverse offsetting the received point cloud data into the watermarked point cloud data. 請求の範囲第７項に記載された電子埋め込み装置で得られた透かし入り点群データから前記電子透かしデータを抽出する抽出装置であって、前記オリジナルランダム点群データと前記透かし入り点群データとを対応付けてそれぞれ離散型フーリエ変換を行って第１及び第２のフーリエ係数列を求めるフーリエ係数生成手段と、前記第１及び前記第２のフーリエ係数列を比較して前記電子透かしデータを抽出する抽出手段とを有することを特徴とする電子情報抽出装置。An extraction device for extracting the digital watermark data from the watermarked point cloud data obtained by the electronic embedding device according to claim 7, wherein the original random point cloud data, the watermarked point cloud data, The digital watermark data is extracted by comparing the Fourier coefficient generation means for obtaining the first and second Fourier coefficient sequences by performing discrete Fourier transform in association with each other and the first and second Fourier coefficient sequences An electronic information extracting apparatus comprising: an extracting means for performing 前記フーリエ係数生成手段は、前記オリジナルランダム点群データが規定されたｘ−ｙ面領域を予め規定された小領域に分割して該小領域毎に含まれる小領域点群を生成する分割手段と、前記小領域点群毎に該小領域点群の頂点と最短距離にある頂点を最短距離頂点として前記透かし入り点群データから探索して前記オリジナルランダム点群データと前記透かし入り点群データとの対応付けを行う対応付け手段とを有することを特徴とする請求の範囲第１０項に記載の電子情報抽出装置。The Fourier coefficient generating means is a dividing means for dividing the xy plane area in which the original random point cloud data is defined into predetermined small areas and generating small area point groups included in the small areas. For each of the small area point groups, the original random point group data and the watermarked point group data are searched from the watermarked point group data using the vertex at the shortest distance from the vertex of the small area point group as the shortest distance vertex. The electronic information extracting apparatus according to claim 10, further comprising: an association unit configured to perform the association. 前記対応付け手段は、前記透かし入り点群データに対する２−ｄ木を生成する２−ｄ木生成手段と、前記小領域点群の頂点位置と埋め込みトレランスとによって規定される問い合わせ領域を設定して前記２−ｄ木内から前記問い合わせ領域に含まれる前記透かし入り点群データを探索して前記オリジナルランダム点群データと前記透かし入り点群データとの対応付けを行う頂点間対応付け手段とを有することを特徴とする請求の範囲第１１項に記載の電子情報抽出装置。The association means sets a query area defined by a 2-d tree generation 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. Vertex-to-vertex association means for searching 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 The electronic information extracting device according to claim 11, wherein ３次元計測に応じて得られたオリジナルランダム点群データに電子情報を電子透かしデータとして埋め込む際に電子計算機で用いられる電子情報埋め込みプログラムであって、前記オリジナルランダム点群データについて離散フーリエ変換を行ってフーリエ係数列を求める第１の手続きと、前記電子透かしデータに応じて前記フーリエ係数列を変更して透かし入りフーリエ係数列とする第２の手続きと、前記透かし入りフーリエ係数列を逆離散フーリエ変換して該逆離散フーリエ変換に応じて透かし入り点群データを得る第３の手続きとを有することを特徴とする電子情報埋め込みプログラム。An electronic information embedding program used by an electronic computer when embedding electronic information as electronic watermark data in original random point cloud data obtained according to three-dimensional measurement, and performing a discrete Fourier transform on the original random point cloud data A first procedure for obtaining a Fourier coefficient sequence, a second procedure for changing the Fourier coefficient sequence according to the digital watermark data to obtain a watermarked Fourier coefficient sequence, and converting the watermarked Fourier coefficient sequence into an inverse discrete Fourier transform. And a third procedure for obtaining watermarked point cloud data according to the inverse discrete Fourier transform. 前記第１の手続きは、前記オリジナルランダム点群データが規定されたｘ−ｙ面領域を予め規定された小領域に分割して該小領域毎に含まれる点群を生成する第４の手続きと、前記点群毎にそのｘ，ｙ座標値を前記点群の重心が原点となるようにオフセットしてオフセット点群に変換する第５の手続きと、前記オフセット点群毎に前記離散フーリエ変換を行ってフーリエ係数列を求める第６の手続きとを有することを特徴とする請求の範囲第１３項に記載の電子情報埋め込みプログラム。The first procedure includes a fourth procedure for dividing the xy plane area in which the original random point cloud data is defined into predetermined small areas and generating a point group included in each small area. , A fifth procedure 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 them into offset point groups; and the discrete Fourier transform for each offset point group The electronic information embedding program according to claim 13, further comprising: a sixth procedure that performs a Fourier coefficient sequence. 前記第３の手続きは、前記透かし入りフーリエ係数列を逆離散フーリエ変換して透かし入り複素数列を生成する第７の手続きと、該透かし入り複素数列について埋め込みに起因する座標値誤差のトレランスを満たす最適透かし埋め込み強度を求める第８の手続きと、前記最適埋め込み強度に基づいて前記フーリエ係数列を再度変更して透かし入りフーリエ係数列とする第９の手続きと、該透かし入りフーリエ係数列を逆オフセットして前記透かし入り点群データとする第１０の手続きとを有することを特徴とする請求の範囲第１４項に記載の電子情報埋め込みプログラム。The third procedure satisfies 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; a ninth procedure for changing the Fourier coefficient sequence again based on the optimum embedding strength to obtain a watermarked Fourier coefficient sequence; and inversely offsetting the watermarked Fourier coefficient sequence 15. The electronic information embedding program according to claim 14, further comprising a tenth procedure for setting the watermarked point cloud data. 電子計算機上で動作し、請求の範囲第１項，第７項，又は第１３項において得られた透かし入り点群データから前記電子透かしデータを抽出するためのプログラムであって、前記オリジナルランダム点群データと前記透かし入り点群データとを対応付けてそれぞれ離散型フーリエ変換を行って第１及び第２のフーリエ係数列を求める第１１の手続きと、前記第１及び前記第２のフーリエ係数列を比較して前記電子透かしデータを抽出する第１２の手続きとを有することを特徴とする電子情報抽出プログラム。A program for operating the electronic computer to extract the electronic watermark data from the watermarked point cloud data obtained in claim 1, claim 7 or claim 13, wherein the original random points An eleventh procedure for obtaining first and second Fourier coefficient sequences by performing discrete Fourier transform in association with group data and the watermarked point group data, and the first and second Fourier coefficient sequences; And a twelfth procedure for extracting the digital watermark data by comparing the electronic watermark data. 前記第１１の手続きは、前記オリジナルランダム点群データが規定されたｘ−ｙ面領域を予め規定された小領域に分割して該小領域毎に含まれる小領域点群を生成する第１３の手続きと、前記小領域点群毎に該小領域点群の頂点と最短距離にある頂点を最短距離頂点として前記透かし入り点群データから探索して前記オリジナルランダム点群データと前記透かし入り点群データとの対応付けを行う第１４の手続きとを有することを特徴とする請求の範囲第１６項に記載の電子情報抽出プログラム。The eleventh procedure is to generate a small area 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. The original random point cloud data and the watermarked point cloud by searching from the watermarked point cloud data with the shortest distance vertex being the vertex that is the shortest distance from the vertex of the small area point cloud for each small area point cloud The electronic information extraction program according to claim 16, further comprising a fourteenth procedure for associating with data. 前記第１４の手続きは、前記透かし入り点群データに対する２−ｄ木を生成する第１５の手続きと、前記小領域点群の頂点位置と埋め込みトレランスとによって規定される問い合わせ領域を設定して前記２−ｄ木内から前記問い合わせ領域に含まれる前記透かし入り点群データを探索して前記オリジナルランダム点群データと前記透かし入り点群データとの対応付けを行う第１６の手続きとを有することを特徴とする請求の範囲第１７項に記載の電子情報抽出プログラム。The fourteenth procedure sets a query region defined by a fifteenth procedure for generating a 2-d tree for the watermarked point group data, a vertex position of the small region point group, and an embedding tolerance, and A sixteenth procedure for searching for the watermarked point cloud data included in the inquiry area from a 2-d tree and associating the original random point cloud data with the watermarked point cloud data. The electronic information extraction program according to claim 17.
JP2002586705A 2001-12-12 2001-12-12 Aviation Electronic Map Information Creation Method and Aviation Electronic Map Information Creation Device Expired - Fee Related JP3935078B2 (en)

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
PCT/JP2001/010904 WO2003051034A1 (en) 2001-12-12 2001-12-12 Electronic information burying method and extracting method, electronic information burying apparatus and extracting apparatus, and programs therefor

Related Child Applications (1)

Application Number Title Priority Date Filing Date
JP2007040094A Division 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

Publications (2)

Publication Number Publication Date
JPWO2003051034A1 true JPWO2003051034A1 (en) 2005-04-21
JP3935078B2 JP3935078B2 (en) 2007-06-20

Family

ID=11738026

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2002586705A Expired - Fee Related JP3935078B2 (en) 2001-12-12 2001-12-12 Aviation Electronic Map Information Creation Method and Aviation Electronic Map Information Creation Device

Country Status (6)

Country Link
US (2) US7370363B2 (en)
EP (1) EP1458177A4 (en)
JP (1) JP3935078B2 (en)
AU (1) AU2002222622A1 (en)
CA (1) CA2438716A1 (en)
WO (1) WO2003051034A1 (en)

Families Citing this family (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7370363B2 (en) * 2001-12-12 2008-05-06 Kokusai Kogyo Co., Ltd. Electronic information embedding method and extracting method, electronic information burying apparatus and extracting apparatus, and programs therefor
TWI288892B (en) * 2005-12-28 2007-10-21 Inst Information Industry Content protection method for vector graph format
US8144925B2 (en) * 2010-08-24 2012-03-27 Seiko Epson Corporation Mapping based message encoding for fast reliable visible watermarking
CN107977515B (en) * 2017-12-01 2021-06-11 中国直升机设计研究所 Helicopter gravity center envelope design method
CN110108257B (en) * 2019-05-20 2021-05-14 苏州中飞遥感技术服务有限公司 Display method based on matching of aerial images of unmanned aerial vehicle and static planning map
KR102363501B1 (en) * 2019-12-10 2022-02-16 주식회사 라이드플럭스 Method, apparatus and computer program for generating earth surface data from 3-dimensional point cloud data
CN113032743B (en) * 2021-03-19 2023-01-06 华南理工大学 Reversible deformation method of tetrahedral mesh model based on information embedding and application
CN113284036B (en) * 2021-06-15 2022-07-26 苏州科技大学 Vector map reversible watermarking algorithm based on virtual grid interval mapping

Family Cites Families (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE2745679C2 (en) * 1977-10-11 1981-12-17 Hoechst Ag, 6000 Frankfurt Pigment preparations, their manufacture and use
JPS63179972A (en) * 1987-01-21 1988-07-23 Miyoshi Kasei:Kk Polyolefin resin treated pigment and cosmetic
US5298355A (en) * 1992-11-13 1994-03-29 Eastman Kodak Company Toner composition with semi-crystalline polyester wax and method of preparation
US6614914B1 (en) * 1995-05-08 2003-09-02 Digimarc Corporation Watermark embedder and reader
JPH10320583A (en) * 1997-05-22 1998-12-04 Meidensha Corp Method for compressing and expanding data
US6359998B1 (en) * 1998-04-23 2002-03-19 3Com Corporation Method and apparatus for wavelet-based digital watermarking
JP4424769B2 (en) * 1998-09-10 2010-03-03 興和株式会社 Digital watermark embedding device
KR100351485B1 (en) * 1998-10-08 2002-09-05 마츠시타 덴끼 산교 가부시키가이샤 Data processor and data recorded medium
JP3397157B2 (en) * 1999-01-13 2003-04-14 日本電気株式会社 Digital watermark insertion system
US7164778B1 (en) * 1999-01-25 2007-01-16 Nippon Telegraph And Telephone Corporation Digital watermark embedding method, digital watermark embedding apparatus, and storage medium storing a digital watermark embedding program
DE19941061A1 (en) * 1999-08-28 2001-03-01 Clariant Gmbh Process for the production of pigment granules
JP3690726B2 (en) * 2000-04-13 2005-08-31 インターナショナル・ビジネス・マシーンズ・コーポレーション Data processing apparatus, image processing apparatus, and methods thereof
DE10054344A1 (en) * 2000-11-02 2002-05-29 Clariant Gmbh Use of coated pigment granules in electrophotographic toners and developers, powder coatings and ink-jet inks
US6950519B2 (en) * 2001-03-05 2005-09-27 Digimarc Corporation Geographically watermarked imagery and methods
US7370363B2 (en) * 2001-12-12 2008-05-06 Kokusai Kogyo Co., Ltd. Electronic information embedding method and extracting method, electronic information burying apparatus and extracting apparatus, and programs therefor

Also Published As

Publication number Publication date
US20090003602A1 (en) 2009-01-01
US7891006B2 (en) 2011-02-15
WO2003051034A1 (en) 2003-06-19
JP3935078B2 (en) 2007-06-20
AU2002222622A1 (en) 2003-06-23
EP1458177A1 (en) 2004-09-15
US20050018846A1 (en) 2005-01-27
US7370363B2 (en) 2008-05-06
EP1458177A4 (en) 2008-07-16
CA2438716A1 (en) 2003-06-19

Similar Documents

Publication Publication Date Title
Lee et al. Vector watermarking scheme for GIS vector map management
JP2019534510A5 (en)
CN110226186A (en) A kind of method and apparatus for indicating map element and vehicle/robot method and apparatus are positioned based on this
US7891006B2 (en) Electronic information embedding method and extracting method, electronic information burying apparatus and extracting apparatus, and programs therefor
JP2011129082A (en) Three-dimensional object recognition device and three-dimensional object recognition method
JP2006323608A (en) Apparatus and method for creating model of group of three-dimensional structure and system for creating three-dimensional model
JP2006214893A (en) Computer software program for measuring three-dimensional shape of object using object measuring method and computer system
JP6945832B2 (en) 3D shape search method and 3D shape search system
JP4427656B2 (en) Survey data processing method
Benedens Robust watermarking and affine registration of 3d meshes
JP4175551B2 (en) Electronic information embedding method and extraction method, electronic information embedding device and extraction device, and program thereof
JP6958993B2 (en) Civil engineering work performance evaluation system, workmanship evaluation method, and program
EP4036859A1 (en) A system and method for providing improved geocoded reference data to a 3d map representation
JP2019040445A (en) Estimation device and program
Bliss et al. Qualitative assessment of surface roughness of limestone specimens from the orientation of unwrapped triangulated point clouds
CN113124816A (en) Antenna work parameter generation method and device, storage medium and computer equipment
Liu et al. A novel approach to automatic registration of point clouds
KR102635566B1 (en) Apparatus for generating map of navigation system and method thereof
WO2000003346A1 (en) Characterization of objects using laser supplemented photogrammetry
KR20120138606A (en) A system for generating road layer using continuity analysis of ground data
JP6618499B2 (en) 3D map separating apparatus and 3D map separating program
Xie et al. Urban 3D building model applied to true orthoimage generation
JP2001109874A (en) Method and system for processing map data
JP2024046092A (en) Development drawing device and development drawing method
Agugiaroa et al. Definition of a transition surface with the purpose of integration between a laser scanner 3D model and a low resolution DTM

Legal Events

Date Code Title Description
A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20061222

A521 Written amendment

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20070220

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: 20070316

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20070316

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: 20100330

Year of fee payment: 3

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20110330

Year of fee payment: 4

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20120330

Year of fee payment: 5

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20120330

Year of fee payment: 5

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20130330

Year of fee payment: 6

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20130330

Year of fee payment: 6

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20140330

Year of fee payment: 7

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