JP2012122906A - Curvature distribution generation system for cable, computer program of the same and curvature distribution generation method for cable - Google Patents
Curvature distribution generation system for cable, computer program of the same and curvature distribution generation method for cable Download PDFInfo
- Publication number
- JP2012122906A JP2012122906A JP2010275038A JP2010275038A JP2012122906A JP 2012122906 A JP2012122906 A JP 2012122906A JP 2010275038 A JP2010275038 A JP 2010275038A JP 2010275038 A JP2010275038 A JP 2010275038A JP 2012122906 A JP2012122906 A JP 2012122906A
- Authority
- JP
- Japan
- Prior art keywords
- cable
- curvature
- curvature distribution
- points
- distribution generation
- 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
Links
- 238000009826 distribution Methods 0.000 title claims abstract description 129
- 238000000034 method Methods 0.000 title claims description 61
- 238000004590 computer program Methods 0.000 title claims description 5
- 238000009499 grossing Methods 0.000 claims abstract description 91
- 238000012545 processing Methods 0.000 claims abstract description 35
- 238000006073 displacement reaction Methods 0.000 claims description 49
- 230000008569 process Effects 0.000 claims description 43
- 238000005452 bending Methods 0.000 claims description 15
- 238000001514 detection method Methods 0.000 abstract description 17
- 230000006870 function Effects 0.000 description 20
- 239000013598 vector Substances 0.000 description 18
- 238000004364 calculation method Methods 0.000 description 15
- 230000008859 change Effects 0.000 description 10
- 230000010349 pulsation Effects 0.000 description 7
- 238000010586 diagram Methods 0.000 description 6
- 238000006243 chemical reaction Methods 0.000 description 5
- 238000004891 communication Methods 0.000 description 4
- 238000005259 measurement Methods 0.000 description 4
- 230000008054 signal transmission Effects 0.000 description 4
- 238000010521 absorption reaction Methods 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 3
- 229910052751 metal Inorganic materials 0.000 description 3
- 239000002184 metal Substances 0.000 description 3
- 239000012212 insulator Substances 0.000 description 2
- 239000010410 layer Substances 0.000 description 2
- 230000001629 suppression Effects 0.000 description 2
- 238000012935 Averaging Methods 0.000 description 1
- RYGMFSIKBFXOCR-UHFFFAOYSA-N Copper Chemical compound [Cu] RYGMFSIKBFXOCR-UHFFFAOYSA-N 0.000 description 1
- LKJPSUCKSLORMF-UHFFFAOYSA-N Monolinuron Chemical compound CON(C)C(=O)NC1=CC=C(Cl)C=C1 LKJPSUCKSLORMF-UHFFFAOYSA-N 0.000 description 1
- FFBHFFJDDLITSX-UHFFFAOYSA-N benzyl N-[2-hydroxy-4-(3-oxomorpholin-4-yl)phenyl]carbamate Chemical compound OC1=C(NC(=O)OCC2=CC=CC=C2)C=CC(=C1)N1CCOCC1=O FFBHFFJDDLITSX-UHFFFAOYSA-N 0.000 description 1
- 239000011247 coating layer Substances 0.000 description 1
- 229910052802 copper Inorganic materials 0.000 description 1
- 239000010949 copper Substances 0.000 description 1
- 230000004069 differentiation Effects 0.000 description 1
- 239000006185 dispersion Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000002093 peripheral effect Effects 0.000 description 1
- 238000012887 quadratic function Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
Landscapes
- Length-Measuring Devices Using Wave Or Particle Radiation (AREA)
- Electric Cable Arrangement Between Relatively Moving Parts (AREA)
Abstract
Description
本発明は、ケーブルの曲率分布生成システム及びそのコンピュータプログラム、並びにケーブルの曲率分布生成方法に関するものである。 The present invention relates to a cable curvature distribution generation system, a computer program thereof, and a cable curvature distribution generation method.
自動車や産業機器等には、信号送信や電力供給のためのケーブルが用いられている。このようなケーブルには、複数の導電金属製の素線を撚り束ねたものや、素線を撚り束ねたものを芯線として複数本撚り束ね、一本のケーブルを構成したもの等がある。
このようなケーブルには、自動車のドア部等の屈曲可動部分に配索されるものもあり、そのような場合、屈曲可動部分の屈曲動作に伴って複雑に屈曲したり捻れたりし、繰り返し変形することで断線に至ることがある。
Cables for signal transmission and power supply are used in automobiles and industrial equipment. Such cables include those in which a plurality of conductive metal strands are twisted and bundled, and those in which a plurality of strands of strands are twisted and bundled to form a single cable.
Some of these cables are routed to movable bending parts such as automobile doors. In such cases, the cables are bent or twisted in a complicated manner as the bending movable part bends, and repeatedly deformed. Doing so may lead to disconnection.
このようなケーブルの変形を評価するために、屈曲したケーブルの長手方向における曲率分布を求めることが求められている。本発明者は、ケーブルの曲率分布を求める手法について、特許文献1にて提案している。 In order to evaluate such deformation of the cable, it is required to obtain a curvature distribution in the longitudinal direction of the bent cable. The inventor has proposed a technique for obtaining the curvature distribution of a cable in Patent Document 1.
特許文献1記載の手法では、まず、ケーブルがCTスキャナで撮像され、複数の断面画像データが取得される。その断面画像データから、ケーブル断面の中心位置が検出される。そして、複数の中心位置からケーブル長手方向の曲率分布が生成される。 In the method described in Patent Document 1, first, a cable is imaged by a CT scanner, and a plurality of cross-sectional image data is acquired. From the cross-sectional image data, the center position of the cable cross-section is detected. And the curvature distribution of a cable longitudinal direction is produced | generated from a some center position.
本発明者は、さらに研究を進めた結果、生成される曲率分布の精度向上のための新たな着想を得た。つまり、断面画像データからケーブル断面の中心位置を画像処理により検出する際に、当該中心位置の検出誤差が生じ、その検出誤差がノイズとなることがある。本発明者は、このノイズを低減させると、曲率分布の精度を向上できることを見出した。 As a result of further research, the present inventor has obtained a new idea for improving the accuracy of the generated curvature distribution. That is, when the center position of the cable section is detected from the section image data by image processing, a detection error of the center position may occur, and the detection error may be noise. The present inventor has found that the accuracy of the curvature distribution can be improved by reducing this noise.
本発明は、かかる観点からなされたものであり、生成される曲率分布の精度向上を目的とする。 The present invention has been made from this point of view, and aims to improve the accuracy of the generated curvature distribution.
(1)本発明は、ケーブルの曲率分布を生成するシステムであって、ケーブルの長手方向における複数の点それぞれの位置情報を取得する取得部と、前記複数の点の前記位置情報に基づいて前記複数の点の曲率を算出することで、前記ケーブルの長手方向における曲率分布を生成する曲率分布生成部と、前記曲率分布、及び/又は、前記複数の点の前記位置情報の平滑化処理を行う平滑化処理部と、を備えているケーブル曲率分布生成システムである。
上記本発明によれば、ケーブルの複数の点それぞれについて、検出誤差が生じても、その検出誤差によるノイズを、平滑化処理によって抑制することができる。
(1) The present invention is a system that generates a curvature distribution of a cable, the acquisition unit acquiring position information of each of a plurality of points in the longitudinal direction of the cable, and the above-described information based on the position information of the plurality of points. By calculating the curvature of a plurality of points, a curvature distribution generating unit that generates a curvature distribution in the longitudinal direction of the cable, and smoothing processing of the curvature distribution and / or the position information of the plurality of points is performed. A cable curvature distribution generation system including a smoothing processing unit.
According to the present invention, even if a detection error occurs at each of a plurality of points on the cable, noise due to the detection error can be suppressed by the smoothing process.
(2)前記平滑化処理部は、複数のケーブル曲げ方向ごとに、前記曲率分布、及び/又は、前記複数の点の前記位置情報の平滑化処理を行うのが好ましい。この場合、ケーブル曲げ方向ごとの平滑化処理が可能である。 (2) It is preferable that the smoothing processing unit smoothes the curvature distribution and / or the position information of the plurality of points for each of a plurality of cable bending directions. In this case, a smoothing process for each cable bending direction is possible.
(3)前記曲率分布生成部は、複数のケーブル曲げ方向ごとの曲率分布を示す複数の方向成分別曲率分布を生成するのが好ましい。この場合、曲率分布を曲げの方向を考慮して評価することができる。例えば、ケーブルを左右に繰り返し曲げた場合の曲率変化の的確な把握に好適である。例えば、ケーブルを左右に同程度繰り返し曲げた場合を想定する。ここで、左に曲げた場合の曲率をKLとし、右に曲げた場合の曲率をKRとする。この場合、曲げの方向を考慮しない絶対曲率で考えると、KL=KRであり、ケーブルを左に曲げた状態から右に曲げた状態に変化させても、絶対曲率は変化しない(ΔK=KL−KR=0)。
一方、方向成分別曲率分布であれば、KL=−KRであり、左に曲げた状態から右に曲げた状態に変化させると、ΔK=KL−KR=2KLであり、曲率変化の的確な把握が可能となる。
(3) It is preferable that the curvature distribution generation unit generates a plurality of directional component-specific curvature distributions indicating a curvature distribution for each of a plurality of cable bending directions. In this case, the curvature distribution can be evaluated in consideration of the bending direction. For example, it is suitable for accurately grasping the curvature change when the cable is repeatedly bent left and right. For example, assume a case where the cable is repeatedly bent left and right to the same extent. Here, a curvature when bent to the left and K L, the curvature when bent to the right with K R. In this case, when considering the absolute curvature without considering the bending direction, K L = K R , and even if the cable is changed from the bent state to the right state, the absolute curvature does not change (ΔK = K L -K R = 0).
On the other hand, in the case of the directional component-specific curvature distribution, K L = −K R , and when changed from the left bent state to the right bent state, ΔK = K L −K R = 2K L and the curvature It is possible to accurately grasp changes.
(4)前記平滑化処理部は、前記複数の方向成分曲率分布それぞれについて平滑化処理を行うことができる。 (4) The smoothing processing unit can perform a smoothing process on each of the plurality of directional component curvature distributions.
(5)前記複数の点の前記位置情報は、それぞれ、局所座標における座標値として示され、前記複数の点に対応する複数の前記局所座標は、それぞれ、各位置情報が示す点又はその近傍におけるケーブル長手方向が第1の座標軸方向として設定され、他の第2及び第3の座標軸方向が、前記第1の座標軸方向に直交する平面内に設定された3次元座標であり前記曲率分布生成部は、複数のケーブル曲げ方向ごとの方向成分別曲率分布を生成するのが好ましい。 (5) The position information of the plurality of points is indicated as coordinate values in local coordinates, and the plurality of local coordinates corresponding to the plurality of points are respectively in the points indicated by the position information or in the vicinity thereof. The longitudinal direction of the cable is set as the first coordinate axis direction, and the other second and third coordinate axis directions are three-dimensional coordinates set in a plane orthogonal to the first coordinate axis direction, and the curvature distribution generation unit Preferably generates a curvature distribution by direction component for each of a plurality of cable bending directions.
(6)前記位置情報は、前記複数の点における各点の隣接する点からの変位が、前記局所座標を用いて表されたものであり、前記平滑化処理部は、前記局所座標における前記第2座標軸方向及び第3座標軸方向ごとに、前記変位の平滑化処理を行うのが好ましい。 (6) In the position information, a displacement from an adjacent point of each of the plurality of points is expressed using the local coordinates, and the smoothing processing unit is configured to perform the first processing in the local coordinates. The displacement smoothing process is preferably performed for each of the two coordinate axis directions and the third coordinate axis direction.
(7)前記曲率分布生成部は、複数の点それぞれの曲率を、局所座標における座標値の2次微分値として算出するのが好ましい。この場合、複数の点の位置情報から、方向成分ごとの曲率を容易に得ることができるとともに、演算速度が向上し、大量データの処理に適したものとなる。 (7) It is preferable that the curvature distribution generation unit calculates the curvature of each of the plurality of points as a secondary differential value of coordinate values in local coordinates. In this case, the curvature for each direction component can be easily obtained from the position information of a plurality of points, and the calculation speed is improved, which is suitable for processing a large amount of data.
(8)前記位置情報は、前記複数の点における各点の隣接する点からの変位が、前記局所座標を用いて表され、前記平滑化処理部は、前記複数の方向成分曲率分布それぞれについて平滑化処理を行うのが好ましい。 (8) In the position information, a displacement of each of the plurality of points from an adjacent point is expressed using the local coordinates, and the smoothing processing unit smoothes each of the plurality of directional component curvature distributions. It is preferable to perform the conversion treatment.
(9)前記曲率分布生成部は、平滑化処理がなされた複数の方向成分別別曲率分布に基づいて、絶対曲率分布を生成するのが好ましい。 (9) It is preferable that the curvature distribution generation unit generates an absolute curvature distribution based on a plurality of directional component-specific curvature distributions that have been smoothed.
(10)前記平滑化処理部は、曲率分布生成の対象となっているケーブルの長手方向範囲を複数に分割した区間毎に、平滑化処理を行うのが好ましい。 (10) It is preferable that the smoothing processing unit performs the smoothing process for each section obtained by dividing the longitudinal range of the cable that is the object of curvature distribution generation into a plurality of sections.
(11)他の観点からみた本発明は、コンピュータを、(1)項〜(10)項のいずれか1項に記載のケーブル曲率分布生成システムとして機能させるためのコンピュータプログラムである。 (11) The present invention viewed from another viewpoint is a computer program for causing a computer to function as the cable curvature distribution generation system according to any one of items (1) to (10).
(12)さらに他の観点から見た本発明は、ケーブルの曲率分布をコンピュータにより生成する方法であって、ケーブルの長手方向における複数の点それぞれの位置情報をコンピュータが取得するステップと、前記複数の点の前記位置情報に基づいて前記複数の点の曲率を算出することで、前記ケーブルの長手方向における曲率分布を前記コンピュータが生成するステップと、前記曲率分布、及び/又は、前記複数の点の前記位置情報の平滑化処理を行う前記コンピュータがステップと、を含む方法である。なお、平滑化処理を行うステップは、曲率分布を生成するステップの前に行っても良いし、後に行っても良い。 (12) The present invention viewed from still another viewpoint is a method for generating a curvature distribution of a cable by a computer, wherein the computer acquires position information of each of a plurality of points in the longitudinal direction of the cable; Calculating a curvature distribution in the longitudinal direction of the cable by calculating curvatures of the plurality of points based on the position information of the points; and the curvature distribution and / or the plurality of points. The computer performing the smoothing process of the position information. Note that the step of performing the smoothing process may be performed before or after the step of generating the curvature distribution.
本発明によれば、生成される曲率分布の精度を向上させることができる。 According to the present invention, it is possible to improve the accuracy of the generated curvature distribution.
以下、本発明の好ましい実施形態について図面を参照しながら説明する。
[1.ケーブル曲率分布生成システム]
図1は、ケーブル曲率分布生成システムを示すブロック図である。システムは、被測定物であるケーブルの断面画像を取得するためのX線CTスキャナ装置(測定装置)1と、測定装置1にて得られたデータを処理するパーソナルコンピュータ(処理装置)2と、を備えている。
Hereinafter, preferred embodiments of the present invention will be described with reference to the drawings.
[1. Cable curvature distribution generation system]
FIG. 1 is a block diagram showing a cable curvature distribution generation system. The system includes an X-ray CT scanner device (measuring device) 1 for acquiring a cross-sectional image of a cable that is an object to be measured, a personal computer (processing device) 2 that processes data obtained by the measuring device 1, It has.
パーソナルコンピュータ2は、図に示すように、各部の制御を行うCPU21、キーボードやマウス等の入力部22、ディスプレイ等の表示部23、外部機器とのデータの入出力を行うためのLAN等の通信部24、及び、メモリやハードディスク等から構成される記憶部25を備えている。 As shown in the figure, the personal computer 2 includes a CPU 21 that controls each unit, an input unit 22 such as a keyboard and a mouse, a display unit 23 such as a display, and a LAN communication for performing input / output of data with external devices. And a storage unit 25 including a memory and a hard disk.
記憶部25には、オペレーティングシステムの他、曲率分布生成プログラム26等が格納されている。曲率分布生成プログラム26がコンピュータ2によって実行されることで、コンピュータ2が曲率分布生成システムとして機能する。
また、通信部24には、X線CTスキャナ装置1が接続されており、X線CTスキャナ装置1から出力されるデータを受信する。前記曲率分布生成プログラム26は、X線CTスキャナ装置1から出力されるデータに基づいてパーソナルコンピュータ2にケーブルの曲率分布測定に係る処理を実行させる。なお、X線CTスキャナ装置1とパーソナルコンピュータとは、接続されている必要は無い。
In addition to the operating system, the storage unit 25 stores a curvature distribution generation program 26 and the like. When the curvature distribution generation program 26 is executed by the computer 2, the computer 2 functions as a curvature distribution generation system.
The communication unit 24 is connected to the X-ray CT scanner apparatus 1 and receives data output from the X-ray CT scanner apparatus 1. The curvature distribution generation program 26 causes the personal computer 2 to execute processing related to cable curvature distribution measurement based on data output from the X-ray CT scanner device 1. Note that the X-ray CT scanner device 1 and the personal computer need not be connected.
CTスキャナ装置1は、被測定物であるケーブルの長手方向に沿った複数の断面内部画像を生成する。そして、CTスキャナ装置1は、取得した被測定物の断面内部画像データをパーソナルコンピュータ2に対して出力する。そして、パーソナルコンピュータ2は、前記断面画像データを通信部24により受信し、記憶部25に保存する。記憶部25に格納された曲率分布生成プログラム26は、上述したように、X線CTスキャナ装置1から出力されるデータに基づいて、パーソナルコンピュータ2に、ケーブルの曲率分布生成に係る処理を実行させる。 The CT scanner device 1 generates a plurality of cross-section internal images along the longitudinal direction of a cable that is a measurement object. Then, the CT scanner device 1 outputs the acquired cross-sectional internal image data of the measured object to the personal computer 2. The personal computer 2 receives the cross-sectional image data by the communication unit 24 and stores it in the storage unit 25. As described above, the curvature distribution generation program 26 stored in the storage unit 25 causes the personal computer 2 to execute processing related to cable curvature distribution generation based on the data output from the X-ray CT scanner apparatus 1. .
図2は、曲率分布生成プログラム26に基づく曲率分布生成システムのブロック図を示している。このシステムは、ケーブル4の長手方向における複数の点それぞれの位置情報を取得する取得部31と、前記ケーブル4の長手方向における曲率分布を生成する曲率分布生成部32と、第1及び第2の平滑化処理部33,34を備えている。 FIG. 2 shows a block diagram of a curvature distribution generation system based on the curvature distribution generation program 26. This system includes an acquisition unit 31 that acquires position information of each of a plurality of points in the longitudinal direction of the cable 4, a curvature distribution generation unit 32 that generates a curvature distribution in the longitudinal direction of the cable 4, and first and second Smoothing processing units 33 and 34 are provided.
次に、上記構成を有する形状状態測定システムによって実現される、ケーブルの曲率分布の生成方法について説明する。図3(a)は、前記曲率分布生成方法における測定対象となるケーブルの構成を示した概略図であり、図4は、前記曲率分布生成方法を示すフローチャートである。 Next, a method for generating a cable curvature distribution realized by the shape state measurement system having the above configuration will be described. FIG. 3A is a schematic diagram illustrating a configuration of a cable to be measured in the curvature distribution generation method, and FIG. 4 is a flowchart illustrating the curvature distribution generation method.
図3(a)において、ケーブル4は、その内部に螺旋状に撚り束ねられた7本の芯線41と、その外周側を覆う絶縁体からなる外部被覆層42とを有している。また、各芯線41は、7本の銅線等の導電金属線材を螺旋状に撚り束ねた素線43と、その外周側を覆う絶縁体からなる被覆層44とにより構成されている。すなわち、各芯線41は、それぞれ独立して信号伝達や電力供給を行うことができるように構成されており、一本のケーブル4によって、多数の信号伝達や電力供給を行えるように構成されている。
なお、ケーブル4を構成している芯線41、芯線41を構成している素線43においても、それぞれ、信号送信や電力供給のためのケーブルとみなすことができる。
In FIG. 3A, the cable 4 has seven core wires 41 spirally bundled therein and an outer covering layer 42 made of an insulator that covers the outer peripheral side thereof. Each core wire 41 is constituted by a strand 43 in which conductive metal wires such as seven copper wires are spirally twisted and a covering layer 44 made of an insulator covering the outer periphery thereof. That is, each core wire 41 is configured to be able to perform signal transmission and power supply independently, and is configured to be able to perform a large number of signal transmission and power supply by a single cable 4. .
Note that the core wire 41 constituting the cable 4 and the strand 43 constituting the core wire 41 can also be regarded as cables for signal transmission and power supply, respectively.
[2.曲率分布生成方法]
以下、上記構成のケーブル4の曲率分布の生成方法について、図4のフローチャートに基づいて説明する。
[2.1 ケーブルの3次元情報取得]
まず、上記構成のケーブル4を、治具等によって、実機に配索された状態を再現するように屈曲変形させて固定し、それをX線CTスキャナ装置1で撮影する。これにより、ケーブル4の断面内部画像(ケーブル4の3次元情報)を取得する(ステップS1)。
図3(b)にケーブル4の断面内部画像の一例を示す。図3(b)中の断面内部画像gは、図3(a)のように、X線CTスキャナ装置1による撮影面Aをケーブル4と交差させ断面内部を透過撮影したものである。これらの断面内部画像gは、コンピュータに読み込み可能な画像データとして取得される。
[2. Curvature distribution generation method]
Hereinafter, a method for generating the curvature distribution of the cable 4 having the above configuration will be described with reference to the flowchart of FIG.
[2.1 Three-dimensional information acquisition of cable]
First, the cable 4 having the above configuration is bent and deformed so as to reproduce the state of being routed to an actual machine with a jig or the like, and is imaged by the X-ray CT scanner device 1. Thereby, the cross-sectional internal image of the cable 4 (three-dimensional information of the cable 4) is acquired (step S1).
FIG. 3B shows an example of a cross-sectional internal image of the cable 4. A cross-sectional internal image g in FIG. 3B is a cross-sectional image of the X-ray CT scanner device 1 crossed with the cable 4 as shown in FIG. These cross-sectional internal images g are acquired as image data that can be read by a computer.
また、図3(b)の断面内部画像gは、周囲よりも濃く(暗く)現れる図中破線で示した濃部g1と、濃部g1よりも薄い(明るい)淡部g2とで現される濃淡画像として得られるものであり、本実施形態の場合、ケーブル4において7本の素線43が撚り束ねられた芯線41の断面像が7箇所の濃部g1として現れている。断面内部画像gは、X線の吸収率の差によって、導電金属線材からなる素線43と、被覆層42,44を含むその他の部分と、の間で大きくコントラストが生じるものであり、素線43の部分を濃部として表した場合、その他の部分は淡部として現れ、図3(b)のように芯線41の断面内部を濃部g1として認識することができる。 Also, the cross-sectional internal image g in FIG. 3B is expressed by a dark portion g1 indicated by a broken line in the figure that appears darker (darker) than the surroundings, and a light portion g2 that is lighter (brighter) than the dark portion g1. In this embodiment, a cross-sectional image of the core wire 41 in which the seven strands 43 are twisted and bundled appears as seven dark portions g1. The cross-sectional internal image g has a large contrast between the strand 43 made of a conductive metal wire and the other portions including the coating layers 42 and 44 due to the difference in the X-ray absorption rate. When the portion 43 is represented as a dark portion, the other portions appear as light portions, and the inside of the cross section of the core wire 41 can be recognized as the dark portion g1 as shown in FIG.
なお、一般にX線画像は、X線吸収率の高い部分が明るく現れ、吸収率の低い部分が暗く現れるため、ケーブル4の撮影像は、金属である素線43が存在する部分が明るい淡部、他の部分が暗い濃部となって得られるのであるが、本実施形態の断面内部画像gでは、淡部と濃部とを反転させることで、素線43(芯線41)の断面像が濃部、他の部分が淡部となるように示し、素線43の断面像をより明瞭に把握できるようにしている。 In general, in an X-ray image, a portion where the X-ray absorption rate is high appears bright, and a portion where the absorption rate is low appears dark. The other portion is obtained as a dark dark portion, but in the cross-sectional internal image g of the present embodiment, a cross-sectional image of the strand 43 (core wire 41) is obtained by inverting the light portion and the dark portion. The dark portion and other portions are shown as light portions, so that the cross-sectional image of the strands 43 can be grasped more clearly.
なお、複数の断面内部画像gが並ぶ方向(ケーブルの長手方向にほぼ対応する)を、Z軸方向とし、Z軸方向に直交する一の方向をX軸方向とし、Z軸方向及びX軸方向に直交する方向をY軸方向とする。また、以下では、このXYZ3次元座標を、全体座標という。 Note that the direction in which a plurality of cross-sectional internal images g are arranged (corresponding substantially to the longitudinal direction of the cable) is the Z-axis direction, and one direction orthogonal to the Z-axis direction is the X-axis direction. The direction orthogonal to is the Y-axis direction. Hereinafter, the XYZ three-dimensional coordinates are referred to as global coordinates.
X線CTスキャナ装置1は、Z方向において所定の幅寸法に設定された測定範囲内を所定のピッチ間隔(一定の間隔)Lで、撮影する。これにより、ケーブル4のZ方向における複数の断面内部画像gが得られる。図3(b)は、Z方向に連続して撮影された複数の断面内部画像gを模式的に示しており、図中各矢印は、X,Y,Z方向を示している。断面内部画像gは、X−Y平面に対して平行であるとともに、Z方向に沿って一定のピッチ間隔Lで連続的に複数枚撮影される。
ここで、ケーブル4の長手方向と、全体座標におけるZ方向と、は、できるだけ一致しているのが好ましいが、ケーブル4は、屈曲変形しているので、両者は必ずしも一致していない。
なお、ケーブル4の3次元情報は、CTスキャナ装置4によって断面画像として取得するほか、3次元読み取り装置等によってケーブルの複数の点の座標値を読み取っても良い。
The X-ray CT scanner device 1 captures images within a measurement range set to a predetermined width dimension in the Z direction at a predetermined pitch interval (fixed interval) L. Thereby, the some cross-section internal image g in the Z direction of the cable 4 is obtained. FIG. 3B schematically shows a plurality of cross-section internal images g taken continuously in the Z direction, and each arrow in the drawing indicates the X, Y, and Z directions. The cross-section internal image g is parallel to the XY plane, and a plurality of images are continuously photographed at a constant pitch interval L along the Z direction.
Here, it is preferable that the longitudinal direction of the cable 4 and the Z direction in the overall coordinates coincide as much as possible, but since the cable 4 is bent and deformed, they do not necessarily coincide.
Note that the three-dimensional information of the cable 4 may be acquired as a cross-sectional image by the CT scanner device 4, and the coordinate values of a plurality of points on the cable may be read by a three-dimensional reading device or the like.
上記ステップS1にて取得された複数の断面内部画像gは、X線CTスキャナ装置1からコンピュータ2に送信される。そして、コンピュータ2は、送信された複数の断面内部画像gを通信部24により受信する。
その後、コンピュータ2のCPU21は、曲率分布生成プログラム26に基づいて、以下の処理(ステップS2〜S7)を実行する。
The plurality of cross-sectional internal images g acquired in step S1 are transmitted from the X-ray CT scanner apparatus 1 to the computer 2. Then, the computer 2 receives the transmitted plurality of cross-sectional internal images g by the communication unit 24.
Thereafter, the CPU 21 of the computer 2 executes the following processing (steps S2 to S7) based on the curvature distribution generation program 26.
[2.2 ケーブル中心位置検出]
まず、取得部31の中心位置検出部31aは、複数の断面内部画像gそれぞれについて、ケーブル4の中心位置を検出する(ステップS2)。中心位置検出部31aとして機能するCPU21は、コンピュータ2の記憶部25に格納した複数の断面内部画像gそれぞれについて、画像処理等によって濃淡画像から2値化画像に変換させる。そして、画像g中、ケーブル4の部分を示す濃部g1の輪郭を示す境界線g3の2次元データを取得する。
[2.2 Cable center position detection]
First, the center position detection unit 31a of the acquisition unit 31 detects the center position of the cable 4 for each of the plurality of cross-section internal images g (step S2). The CPU 21 functioning as the center position detection unit 31a converts each of the plurality of cross-sectional internal images g stored in the storage unit 25 of the computer 2 from a grayscale image to a binarized image by image processing or the like. And the two-dimensional data of the boundary line g3 which shows the outline of the dark part g1 which shows the part of the cable 4 in the image g is acquired.
次に、CPU21は、ケーブル4の断面形状である境界線g3の2次元データから、ケーブル中心位置g4を特定する。中心位置は、例えば、境界線g3の2次元データから、当該境界線g3で囲まれた領域の面積中心を求めることで、特定することができる。
本実施形態において、ケーブル中心位置g4は、各画像gにおけるケーブル4の位置を示す基準位置となる。なお、各画像gにおけるケーブル4の位置は、ケーブル4の他の位置を基準として定めても良い。
Next, the CPU 21 specifies the cable center position g4 from the two-dimensional data of the boundary line g3 that is the cross-sectional shape of the cable 4. The center position can be specified by, for example, obtaining the area center of the region surrounded by the boundary line g3 from the two-dimensional data of the boundary line g3.
In the present embodiment, the cable center position g4 is a reference position indicating the position of the cable 4 in each image g. Note that the position of the cable 4 in each image g may be determined with reference to another position of the cable 4.
ここで、複数の断面内部画像gは、X−Y平面に平行かつZ方向に沿ってピッチ間隔Lで連続的に撮影されているため、各断面内部画像gに対してはそれぞれZ方向における座標値を割り当てることができる。従って、CPU21は、上記各断面内部画像g上の点である断面中心g4のX−Y座標値に、それぞれ対応するZ座標値を割り当て、断面中心g4におけるX−Y−Zの3方向の3次元座標値を取得できる。
つまり、CPU21は、画像gに対応するケーブル4の長手方向における複数の点iについて、全体座標における3次元座標値(位置情報)を取得できる。複数の点iについて、全体座標における3次元座標値(位置情報)は、記憶部25に格納される。
なお、複数の断面内部画像gのZ方向におけるそれぞれの間隔を予め把握しておけば、複数の断面内部画像gのZ方向の間隔を一定のピッチ間隔とせずとも、各断面中心g4のX−Y座標値に対応するZ座標値を特定することができる。
Here, since the plurality of cross-sectional internal images g are taken continuously at a pitch interval L along the Z direction and parallel to the XY plane, the coordinates in the Z direction are respectively set for the cross-sectional internal images g. A value can be assigned. Therefore, the CPU 21 assigns corresponding Z coordinate values to the XY coordinate values of the cross-sectional center g4 that are points on the cross-sectional internal images g, and 3 in the three directions of XYZ at the cross-sectional center g4. The dimension coordinate value can be acquired.
That is, the CPU 21 can acquire three-dimensional coordinate values (position information) in the overall coordinates for a plurality of points i in the longitudinal direction of the cable 4 corresponding to the image g. For a plurality of points i, the three-dimensional coordinate values (position information) in the overall coordinates are stored in the storage unit 25.
If the intervals in the Z direction of the plurality of cross-section internal images g are grasped in advance, the X− of each cross-section center g4 can be obtained without setting the intervals in the Z direction of the plurality of cross-section internal images g to be constant pitch intervals. A Z coordinate value corresponding to the Y coordinate value can be specified.
上記のように得られたケーブル中心位置g4の座標値には、画像gにおける中心位置検出の際の検出誤差がノイズとして含まれることがある。 In the coordinate value of the cable center position g4 obtained as described above, a detection error in detecting the center position in the image g may be included as noise.
[2.3 全体座標から局所座標への変換]
続いて、取得部31の局所座標変換部31bは、複数の中心位置g4(点i)それぞれの座標値(全体座標での座標値)を、以下に説明する局所座標の座標値(位置情報)に変換する(ステップS3)。
局所座標は、u軸、v軸、w軸を有する3次元直交座標であり、点iを局所座標原点とする局所座標の各軸を、ui軸,vi軸、wi軸と表すこととする。なお、以下では、ui,vi、wiは、局所座標の軸を表すだけでなく、各軸方向の単位ベクトルを表す場合がある。
[2.3 Conversion from global coordinates to local coordinates]
Subsequently, the local coordinate conversion unit 31b of the acquisition unit 31 uses the coordinate values (coordinate values in global coordinates) of the plurality of center positions g4 (points i) as the coordinate values (position information) of local coordinates described below. (Step S3).
Local coordinates, u-axis, v-axis, a three-dimensional orthogonal coordinate having an w-axis, the respective axes of the local coordinates of the point i and the local coordinate origin, u i axis, v i axis, be expressed as w i axis And In the following, u i , v i , and w i may represent not only axes of local coordinates but also unit vectors in the respective axis directions.
図5(a)に示すように、任意の点iの局所座標のwi軸は、点iに隣接する点(一つ手前の点)i−1から、点iに向かう方向に一致するように設定される。すなわち、全体座標からみると、各点iを原点とする局所座標におけるwi軸は、それぞれ、ケーブルの屈曲により、異なる方向を向き得る。 As shown in FIG. 5 (a), w i axis of the local coordinates of an arbitrary point i is a point adjacent to the point i (one front point) from i-1, to match the direction toward the point i Set to That is, when viewed from the global coordinate, w i axis in the local coordinate each point i and origin, respectively, by bending of the cable may orientations different directions.
各局所座標のui軸及びvi軸は、各局所座標におけるwi軸に直交する平面内に設定される。
各局所座標のui軸は、iの変化によっては大きく変化しないか、又は、点iの位置の変化に応じて連続的に変化し、かつ、wi軸と直交するベクトルの方向に設定される。各局所座標のvi軸は、ui,wi軸と直交するベクトルの方向に設定される。
ui軸及びvi軸を設定するには、例えば、図5(b)に示すように、曲率の算出対象となっている全ての点iのwi軸に対して、平行とならない固定の単位ベクトルu0を予め規定しておく。そして、vi=wi×u0のベクトル演算、及び、ui=vi×wiのベクトル演算により、ui軸及びvi軸の方向を持つ単位ベクトルを算出すればよい。なお、上記ベクトル演算式において「×」は、外積を求めるためのベクトル演算記号である。
The u i axis and the v i axis of each local coordinate are set in a plane orthogonal to the w i axis in each local coordinate.
The u i axis of each local coordinate does not change greatly depending on the change of i, or changes continuously according to the change of the position of the point i, and is set in the direction of the vector orthogonal to the w i axis. The The v i axis of each local coordinate is set in the direction of a vector orthogonal to the u i and w i axes.
In order to set the u i axis and the v i axis, for example, as shown in FIG. 5 (b), a fixed value that is not parallel to the w i axis of all points i for which the curvature is to be calculated. A unit vector u 0 is defined in advance. The vector operation of v i = w i × u 0 , and, by vector operation of u i = v i × w i, may be calculated unit vector having a direction u i axis and v i axis. In the above vector arithmetic expression, “x” is a vector arithmetic symbol for obtaining an outer product.
また、ui軸及びvi軸を設定するには、例えば、図5(c)に示すように、ケーブル4の中心位置g4(位置i)からみた特定位置、例えば、特定の素線43の中心位置g5への単位ベクトルをu0(i)として規定しておく。そして、上記と同様に、vi=wi×u0(i)のベクトル演算、及び、ui=vi×wiのベクトル演算により、ui軸及びvi軸の方向を持つ単位ベクトルを算出すればよい。 Further, in order to set the u i axis and the v i axis, for example, as shown in FIG. 5C, a specific position viewed from the center position g4 (position i) of the cable 4, for example, a specific wire 43 A unit vector to the center position g5 is defined as u 0 (i). Then, similarly to the above, a unit vector having directions of the u i axis and the v i axis is obtained by the vector operation of v i = w i × u 0 (i) and the vector operation of u i = v i × w i. May be calculated.
各点iの局所座標の設定の仕方は、上記のものに限られないが、事後的に局所座標値を全体座標値に変換したり、局所座標でのケーブル形状状態を評価したりできるように、所定のルールに基づいて、一意に決定されるものであればよい。なお、例示した局所座標の設定の仕方によれば、ui軸及びvi軸は、wi軸の方向の連続的な変化に応じて、その方向が連続的に変化するため、点iの位置の変位や曲率を算出するのに適切であるとともに、後述の平滑化にも適切なものとなっている。 The method of setting the local coordinates of each point i is not limited to the above, but it is possible to convert the local coordinate values into global coordinate values later and to evaluate the cable shape state at the local coordinates. As long as it is uniquely determined based on a predetermined rule. In addition, according to the setting method of the illustrated local coordinates, the u i axis and the v i axis change continuously in accordance with the continuous change in the direction of the w i axis. It is suitable for calculating the displacement and curvature of the position, and also suitable for smoothing described later.
局所座標変換部31bとして機能するCPU21は、記憶部25に格納されている全体座標での点iの位置情報(座標値)を、下記式に基づく演算により、上記のように設定される局所座標における座標値に変換する。複数(4以上)の点iの局所座標における座標値(位置情報)は、記憶部25に格納される。
ここで、Aは、任意の点iについての全体座標における座標値(ベクトル値)であり、Bは、局所座標における座標値(ベクトル値)であり、oは、前記局所座標の原点の全体座標における座標値(ベクトル値)であり、eu,ev,ewは、前記局所座標軸u,v,wの単位ベクトルである。A,B,o,eu,ev,ewは、それそれ、次のように定義される(図6参照)。
Here, A is a coordinate value (vector value) in global coordinates for an arbitrary point i, B is a coordinate value (vector value) in local coordinates, and o is the global coordinates of the origin of the local coordinates a coordinate value in the (vector value), e u, e v, e w , the local axes u, v, is a unit vector in w. A, B, o, e u , e v and e w are respectively defined as follows (see FIG. 6).
本実施形態では、点i+1の座標値を示す局所座標の原点は、点iが採用される。したがって、局所座標における点i+1の座標値Bは、隣接する点(一つ手前の点)iからみた、点i+1の変位量を示すことになる。より具体的は、座標値Bは、隣接する点iからみた、ui軸方向の変位と、vi軸方向の変位と、を示しており、これらの変位は、点i+1におけるケーブル4の曲げ度合いを示す情報となっている。また、当然ながら、座標値Bは、wi軸方向の変位も示している。
なお、点i+1の座標値を示す局所座標の原点は、点iに限定されるものではなく、点i+1近傍の他の点であってもよい。
In the present embodiment, the point i is adopted as the origin of the local coordinates indicating the coordinate value of the point i + 1. Therefore, the coordinate value B of the point i + 1 in local coordinates indicates the amount of displacement of the point i + 1 as seen from the adjacent point (one point before) i. More specifically, the coordinate value B indicates the displacement in the u i axis direction and the displacement in the v i axis direction as seen from the adjacent point i, and these displacements are the bending of the cable 4 at the point i + 1. Information indicating the degree. Also, of course, the coordinate values B also shows w i axis direction of displacement.
Note that the origin of the local coordinates indicating the coordinate value of the point i + 1 is not limited to the point i, and may be another point in the vicinity of the point i + 1.
本実施形態では、点i+1の座標値を示す局所座標のwi軸の方向は、点i−1から点iに向かう方向としたが、点i−1から点i+1に向かう方向、又は、点iから点i+1に向かう方向などとしても良い。つまり、点i+1の座標値を示す局所座標のw軸の方向は、点i+1又はその付近のケーブル4の長手方向にほぼ沿っていればよい。 In the present embodiment, the direction of the w i axis of the local coordinates representing the coordinates of the point i + 1 is set to the direction from the point i-1 to the point i, the direction from the point i-1 to the point i + 1, or, point The direction from i to the point i + 1 may be used. That is, the w-axis direction of the local coordinates indicating the coordinate value of the point i + 1 only needs to be substantially along the longitudinal direction of the cable 4 at or near the point i + 1.
[2.4 位置データの平滑化]
第1平滑化処理部33は、各点iの局所座標値の平滑化を行うことで、前記ノイズ(検出誤差)の除去/抑制を行う(ステップS4)。まず、第1平滑化処理部33として機能するCPU21は、記憶部25に格納されている複数の点iの局所座標での座標値(位置情報)を用いて、ui軸方向、vi軸方向ごとの変位データを生成する。
ui軸方向の変位データは、ケーブル4上での長手方向位置を示す指標値Li(ケーブル端部からの距離など)が示す位置(点i)におけるケーブル4のui軸方向変位(隣接する点i−1からの変位)を示すものである。つまり、ui軸方向の変位データは、複数の点iの局所座標値におけるui軸の値の集合から構成される。
vi軸方向の変位データは、前記指標値Liが示す位置(点i)におけるケーブル4のvi軸方向変位(隣接する点i−1からの変位)を示すものである。つまり、vi軸方向の変位データは、複数の点iの局所座標値におけるvi軸の値の集合から構成される。
[2.4 Smoothing of position data]
The first smoothing processing unit 33 removes / suppresses the noise (detection error) by smoothing the local coordinate value of each point i (step S4). First, the CPU 21 functioning as the first smoothing processing unit 33 uses the coordinate values (position information) at the local coordinates of a plurality of points i stored in the storage unit 25 to use the u i axis direction and the v i axis. Generate displacement data for each direction.
The displacement data in the u i axis direction is the displacement of the cable 4 in the u i axis direction (adjacent) at the position (point i) indicated by the index value L i (the distance from the cable end or the like) indicating the position in the longitudinal direction on the cable 4. (Displacement from point i-1). That is, the displacement data in the u i axis direction is composed of a set of u i axis values in the local coordinate values of a plurality of points i.
v i axis direction displacement data indicates the v i axial displacement of the cable 4 in the index value L i indicates the position (point i) (displacement from i-1 adjacent points). That is, the displacement data in the v i axis direction is composed of a set of v i axis values in the local coordinate values of a plurality of points i.
CPU21は、これらのui軸方向、vi軸方向ごとの変位データそれぞれに対して、平滑化処理(第1の平滑化処理)を行う。第1の平滑化処理により、中心位置の検出誤差により生じた変位データそれぞれの脈動を抑制する。平滑化処理は、検出誤差が、ほぼサンプリング周波数で発生するという特徴を用いて行う。 The CPU 21 performs a smoothing process (first smoothing process) on each of the displacement data for each of the u i axis direction and the v i axis direction. By the first smoothing process, the pulsation of each displacement data caused by the detection error of the center position is suppressed. The smoothing process is performed using the feature that the detection error occurs at almost the sampling frequency.
第1の平滑化処理は、例えば、変位データを、その変位データの脈動に追従しない関数(1回微分して0にならない関数)で近似し、その近似関数を平滑化された変位データとすることで行われる。
近似関数を多項式とする場合、点iの数Nに対して、N/2次以下の多項式とするのが好ましい。図7は、方向ごとの変位データ(N−1次関数)を、N/2次関数で近似した様子を示している。図7から明らかなように、N/2次関数で近似することで、変位データは、脈動が少なくなっている。
In the first smoothing process, for example, the displacement data is approximated by a function that does not follow the pulsation of the displacement data (a function that is differentiated once and does not become 0), and the approximate function is used as smoothed displacement data. Is done.
When the approximation function is a polynomial, it is preferable to use a polynomial of N / 2 order or less with respect to the number N of points i. FIG. 7 shows a state in which displacement data (N-1 order function) for each direction is approximated by an N / 2 order function. As is apparent from FIG. 7, the displacement data has less pulsation by approximating with an N / 2 quadratic function.
近似関数を三角関数(sin関数,cos関数)とする場合は、基底関数を、sin(αL),cos(αL)として、α≦π(N−1)/Lm(ここで、Lmは、Liの最大値)とする)とするのが好ましい。図8は、方向ごとの変位データ(基底関数におけるα=πN/Lm)を、基底関数におけるα=π(N−1)/Lmとした近似関数で近似した様子を示している。図8においても、変位データの脈動が少なくなっている。 When the approximate function is a trigonometric function (sin function, cos function), the basis functions are sin (αL) and cos (αL), and α ≦ π (N−1) / L m (where L m is preferably, the maximum value) to) of L i. FIG. 8 shows a state in which displacement data for each direction (α = πN / L m in the basis function) is approximated by an approximation function with α = π (N−1) / L m in the basis function. Also in FIG. 8, the pulsation of the displacement data is reduced.
なお、第1の平滑化処理は、近似関数による近似のほか、ローパスフィルタによって脈動を抑制/除去することによっておこなってもよい。 Note that the first smoothing process may be performed by suppressing / removing pulsation by a low-pass filter, in addition to approximation by an approximation function.
また、図8に示すように、複数の点iが属するケーブル範囲(曲率分布生成の対象となっているケーブルの長手方向範囲;Li=0〜Lm)を、複数の区間P1,P2に分割しておき、複数に分割した区間毎に、上述のような平滑化処理を行っても良い。ケーブル4の大部分が直線状で、ケーブル端部だけが部分的に屈曲しているような場合のように、単一の近似関数で近似するのが困難な場合もある。このような場合であっても、第1平滑化処理部33では、複数の区間への分割と、区間毎に異なる近似関数の適用と、を行うことで、適切な平滑化を行うことができる。なお、複数区間への分割は、システムユーザからの分割ポイントの指定を受けて行っても良いし、よい近似が得られない場合に、第1平滑化処理部33が自動的に行ってもよい。 Further, as shown in FIG. 8, a cable range to which a plurality of points i belong (longitudinal range of a cable for which a curvature distribution is generated; L i = 0 to L m ) is assigned to a plurality of sections P1 and P2. The above-described smoothing process may be performed for each of the divided sections. It may be difficult to approximate with a single approximation function, such as when the majority of the cable 4 is straight and only the cable end is partially bent. Even in such a case, the first smoothing processing unit 33 can perform appropriate smoothing by performing division into a plurality of sections and application of different approximation functions for each section. . The division into a plurality of sections may be performed in response to the designation of the division point from the system user, or may be automatically performed by the first smoothing processing unit 33 when a good approximation cannot be obtained. .
なお、複数の区間P1,P2それぞれの近似関数は、曲率の連続性を保持するため、区間P1,P2の境界において1次微分値がほぼ等しくなるものを用いるのが好ましい。 In addition, it is preferable to use the approximate function of each of the plurality of sections P1 and P2 in which the primary differential values are substantially equal at the boundary between the sections P1 and P2 in order to maintain the continuity of curvature.
変位データには、検出誤差等により生じた位置変動誤差が含まれているが、上記のような平滑化処理を行うことで、変位データから算出される曲率の誤差を抑制できる。
平滑化処理された変位データは、記憶部25に格納される。
The displacement data includes a position variation error caused by a detection error or the like. By performing the smoothing process as described above, an error in curvature calculated from the displacement data can be suppressed.
The smoothed displacement data is stored in the storage unit 25.
[2.5 曲げ方向成分別の曲率算出]
[2.5.1 第1の曲率算出方法(外接円近似)]
曲率分布生成部32の曲率算出部32aは、記憶部25に格納された平滑化処理後の変位データから、各方向ui,viごとの曲率分布を生成する(ステップS5)。
曲率算出部32aとして機能するCPU21は、まず、平滑化処理されたui軸方向、vi軸方向ごとの変位データ(位置情報)を、全体座標における各点iごとの座標値(位置情報)に戻す。そして、CPU21は、複数の点iの中の任意の点iの曲率を求める場合、その点iの両隣の点i−1,i+1の座標を、それぞれ、その点iを原点とする局所座標の座標値に変換する。
[2.5 Calculation of curvature by bending direction component]
[2.5.1 First curvature calculation method (circumscribed circle approximation)]
Curvature calculating portion 32a of the curvature distribution generating unit 32, from the displacement data after stored smoothing process in the storage unit 25, generates the curvature distribution for each direction u i, v i (step S5).
CPU21 functioning as curvature calculating section 32a, first, the smoothing treated u i axis direction, v i axis displacement data for each direction (position information), the coordinate values of each point i in the global coordinate (position information) Return to. And when calculating | requiring the curvature of the arbitrary points i in the some point i, the CPU 21 uses the coordinates of the points i−1 and i + 1 on both sides of the point i as local coordinates with the point i as the origin. Convert to coordinate values.
さらに、CPU21は、点i−1,i+1の局所座標値から、外接円の式に基づき、点iの曲率を算出する。曲率の算出は、局所座標のvi−wi面に投影した曲率Ku,ui−wi面に投影した曲率Kvそれぞれついて別々に行われる。つまり、複数の方向成分ごとの曲率Ku,Kvが算出される。なお、vi−wi面に投影した曲率Kuについては、各点についてのui軸の座標値を0として演算すればよい。
点iにおける曲率Ku,Kvの算出は、例えば、外接円の式を用いて行うことができる。
外接円の式による曲率の算出には、図10に示すように、点i−1をA、点iをO,点i+1をBとし、ベクトルOA,OBを算出する。そして、OA,OBの角度を算出する。OA,OBの角度θは、cosθ=OA・OB/((|OA|)|OB|)により求めることが可能である。
この角度θを用いて、外接円の式(R=|AB|/sinθ)により、曲率半径Rを算出できる。そして、この曲率半径Rより、曲率K=1/Rを求めることができる。
Further, the CPU 21 calculates the curvature of the point i from the local coordinate values of the points i−1 and i + 1 based on the circumscribed circle formula. Calculation of the curvature, the local coordinates of v i -w i curvature projected on the plane K u, is performed separately with each u i -w i curvature projected on the plane K v. That is, the curvatures K u and K v for each of the plurality of direction components are calculated. Note that the curvature K u projected on the v i -w i plane may be calculated by setting the coordinate value of the u i axis for each point to zero.
The curvatures K u and K v at the point i can be calculated using, for example, a circumscribed circle formula.
In calculating the curvature by the circumscribed circle formula, as shown in FIG. 10, the point i-1 is A, the point i is O, the point i + 1 is B, and vectors OA and OB are calculated. Then, the angles of OA and OB are calculated. The angle θ between OA and OB can be obtained by cos θ = OA · OB / ((| OA |) | OB |).
Using this angle θ, the radius of curvature R can be calculated by the circumscribed circle equation (R = | AB | / sin θ). From this curvature radius R, the curvature K = 1 / R can be obtained.
[2.5.2 第2の曲率算出方法(ケーブル軸(w軸)からの変化量2次微分換算)]
各軸ui,viまわりの曲率Ku,Kvは、次の式によって求めることもできる。
The curvatures K u and K v around the axes u i and v i can also be obtained by the following equations.
上記式によって曲率Ku,Kvを求めるには、まず、第1の平滑化処理後の変位データを用いて、点iと点i+1間のZ微分値を、次のように求める。
点i−1,点iにおける微分値は、先の定義上、
である。
したがって、点i−1から点i間の距離をsiとし、点iから点i+1間の距離をsi+1とすると、数値解析手法における差分法に基づき、曲率Ku,Kv≒2次微分値は、下記にように算出できる(図11参照)。
It is.
Therefore, assuming that the distance between the point i-1 and the point i is s i and the distance between the point i and the point i + 1 is s i + 1 , the curvatures K u , K v ≈secondary differentiation based on the difference method in the numerical analysis method. The value can be calculated as follows (see FIG. 11).
上記のようにして各点iについて算出された曲率Kuからなる曲率分布(第1の方向成分別別曲率分布)と、各点iについて算出された曲率Kvからなる曲率分布(第2の方向成分別曲率分布)は、記憶部25に格納される。
なお、上記演算方法は、一例であり、局所座標の定義の仕方に応じて演算方法を適宜変更すればよい。例えば、ケーブル長手方向(w軸方向)を、点i−1から点i+1に向かう方向、又は、点iから点i+1に向かう方向とした場合、それに応じて演算方法を変更すればよい。
また、形状に与えられた微分方程式を解くことができる数値解析手法として、上記では、差分法を採用したが、有限要素法などの他の数値解析手法を用いて、軸からの位置変動の2次微分(曲率成分)を求めることもできる。
The curvature distribution (first curvature component-specific curvature distribution) composed of the curvature K u calculated for each point i as described above and the curvature distribution (second curvature distribution) composed of the curvature K v calculated for each point i. The direction component curvature distribution) is stored in the storage unit 25.
In addition, the said calculation method is an example and what is necessary is just to change a calculation method suitably according to the method of the definition of a local coordinate. For example, when the cable longitudinal direction (w-axis direction) is a direction from point i-1 to point i + 1 or a direction from point i to point i + 1, the calculation method may be changed accordingly.
In the above, the difference method is adopted as a numerical analysis method capable of solving the differential equation given to the shape. However, by using another numerical analysis method such as a finite element method, the position variation from the axis 2 A second derivative (curvature component) can also be obtained.
[2.6 曲率分布の平滑化]
第2平滑化処理部34は、曲率分布の平滑化を行うことで、曲率分布に含まれるノイズの除去/抑制を行う(ステップS6)。第2平滑化処理部34による平滑化処理は、平滑化の対象が異なるだけで、第1平滑化処理部34による平滑化処理と同様の処理である。
[2.6 Smoothing of curvature distribution]
The second smoothing processing unit 34 removes / suppresses noise included in the curvature distribution by smoothing the curvature distribution (step S6). The smoothing process by the second smoothing processing unit 34 is the same as the smoothing process by the first smoothing processing unit 34 except that the smoothing target is different.
つまり、第2平滑化処理部33として機能するCPU21は、第1及び第2の方向成分別曲率分布それぞれに対して、ui軸方向、vi軸方向ごとの変位データに対して行った第1の平滑化処理と同様の平滑化処理(第2の平滑化処理)を行う。
なお、第2の平滑化処理においても、複数区間P1,P2に分割した上での平滑化処理が可能である。
平滑化処理された第1及び第2の方向成分別曲率分布は、記憶部25に格納される。
In other words, the CPU 21 functioning as the second smoothing processing unit 33 performs the first change performed on the displacement data for each of the u i axis direction and the v i axis direction for each of the first and second direction component curvature distributions. A smoothing process (second smoothing process) similar to the smoothing process 1 is performed.
In the second smoothing process as well, the smoothing process can be performed after being divided into a plurality of sections P1 and P2.
The smoothed first and second direction component curvature distributions are stored in the storage unit 25.
[2.7 曲率分布出力]
曲率分布生成部32の絶対曲率算出部32は、平滑化処理された第1及び第2の方向成分別曲率分布から、絶対曲率分布を求める。そして、曲率分布生成部32は、平滑化処理された第1・第2の方向成分別曲率分布、及び/又は、絶対曲率分布を出力する(ステップS7)。
[2.7 Curvature distribution output]
The absolute curvature calculation unit 32 of the curvature distribution generation unit 32 obtains an absolute curvature distribution from the smoothed first and second direction component-specific curvature distributions. Then, the curvature distribution generation unit 32 outputs the smoothed first and second direction component-specific curvature distributions and / or absolute curvature distributions (step S7).
絶対曲率分布は、第1及び第2の方向成分別曲率分布のように、曲げ方向成分毎の曲率分布ではなく、方向とは無関係に、各点iにおけるケーブル曲率の分布を示すものである。
絶対曲率算出部32は、各点iにおける絶対曲率Kiを、Ki=Kui2+Kvi2によって算出する。
The absolute curvature distribution is not a curvature distribution for each bending direction component, like the first and second direction component curvature distributions, but indicates the distribution of the cable curvature at each point i regardless of the direction.
The absolute curvature calculator 32 calculates an absolute curvature K i at each point i by K i = K ui2 + K vi2 .
なお、システムは、曲率分布のほか、変位データを出力してもよい。変位データは、局所座標によるものであってもよいが、ケーブル形状の評価をする上で、全体座標が好ましい場合には、下記式に基づいて、局所座標値Bを全体座標値Aに変換して出力すればよい。
[3.考察]
以上のような本実施形態に係る曲率分布生成システムによれば、検出誤差を含む変位データから、そのまま曲率を算出した場合の曲率のばらつきを抑えることができる。つまり、図12に示すように、変位データの脈動は、隣接する点間で曲率を大きく変動させ、曲率分布の精度を低下させる。したがって、検出誤差が大きい場合には、変位データを平滑化することで、曲率のばらつきを抑えることができる。
また、変位データに脈動が含まれたまま曲率を求めると、図13において「平滑化なし」のデータとして示すように、各方向成分の曲率として、極めて大きな正・負の曲率が交互に発生する。このような曲率分布から絶対曲率を求めると、図14において「平滑化なし」のデータとして示すように、絶対曲率に、極めて大きなノイズが発生する。
[3. Discussion]
According to the curvature distribution generation system according to the present embodiment as described above, it is possible to suppress variation in curvature when the curvature is calculated as it is from the displacement data including the detection error. That is, as shown in FIG. 12, the pulsation of the displacement data causes the curvature to fluctuate greatly between adjacent points, thereby reducing the accuracy of the curvature distribution. Therefore, when the detection error is large, the variation in curvature can be suppressed by smoothing the displacement data.
Further, when the curvature is obtained with the pulsation included in the displacement data, extremely large positive and negative curvatures are alternately generated as the curvature of each direction component, as shown as “no smoothing” data in FIG. . When the absolute curvature is obtained from such a curvature distribution, extremely large noise is generated in the absolute curvature, as shown as “no smoothing” data in FIG.
これに対して、方向成分ごとに分離した変位データ又は曲率分布に対して、平滑化処理を行うと、図13において、「第1の平滑化のみ」及び「第2の平滑化のみ」のデータとして示すように、「平滑化なし」のデータに比べて、ノイズの少ない曲率分布が得られる。
また、図14において「平滑化データ」のデータとして示すように、第1又は第2の平滑化処理を行うと、「平滑化なし」のデータに比べて、ノイズを抑制した絶対曲率分布が得られる。
On the other hand, when the smoothing process is performed on the displacement data or the curvature distribution separated for each directional component, data of “first smoothing only” and “second smoothing only” in FIG. As can be seen from the graph, a curvature distribution with less noise is obtained as compared with the data of “no smoothing”.
Further, as shown as “smoothed data” data in FIG. 14, when the first or second smoothing process is performed, an absolute curvature distribution in which noise is suppressed is obtained as compared with the data of “no smoothing”. It is done.
図14において、「絶対曲率算出後平滑化」で示されるデータは、平滑化せずに絶対曲率分布を求めた後に、平滑化した場合のデータである。この場合、「平滑化なし」のデータにおけるノイズを平均化したノイズが残っており、「平滑化なし」よりはノイズが抑制されているものの、第1の平滑化処理又は第2の平滑化処理を行った場合よりも、ノイズ抑制効果が低下した。
したがって、第1の平滑化処理又は第2の平滑化処理のように、複数の方向成分ごとのデータ(変位データ又は曲率分布)に対して平滑化処理を行ったほうが、ノイズ抑制効果が高いことがわかる。
In FIG. 14, data indicated by “smoothing after calculation of absolute curvature” is data when smoothing is performed after obtaining an absolute curvature distribution without smoothing. In this case, the noise obtained by averaging the noise in the “no smoothing” data remains, and the noise is suppressed as compared with “no smoothing”, but the first smoothing process or the second smoothing process is performed. The noise suppression effect was lower than when performing the above.
Therefore, the noise suppression effect is higher when the smoothing process is performed on the data (displacement data or curvature distribution) for each of the plurality of direction components as in the first smoothing process or the second smoothing process. I understand.
本実施形態では、第1の平滑化処理と第2の平滑化処理の双方を実行しているが、いずれか一方だけであってもよい。いずれか一方の平滑化処理を実行する場合、より次数が低い曲率分布の平滑化(第2の平滑化処理)の方を実行するのが好ましい。特に、変位データにおける検出誤差が比較的小さければ、第2の平滑化処理だけでも十分である。 In the present embodiment, both the first smoothing process and the second smoothing process are executed, but only one of them may be performed. When either one of the smoothing processes is executed, it is preferable to execute the smoothing of the curvature distribution having a lower degree (second smoothing process). In particular, if the detection error in the displacement data is relatively small, only the second smoothing process is sufficient.
なお、今回開示された実施の形態はすべての点で例示であって制限的なものではないと考えられるべきである。本発明の範囲は、上記した意味ではなく、特許請求の範囲によって示され、特許請求の範囲と均等の意味、及び範囲内でのすべての変更が含まれることが意図される。 The embodiment disclosed this time should be considered as illustrative in all points and not restrictive. The scope of the present invention is defined by the terms of the claims, rather than the meanings described above, and is intended to include any modifications within the scope and meaning equivalent to the terms of the claims.
2 コンピュータ
26 曲率分布生成コンピュータプログラム
31 取得部
31a 中心位置検出部
31b 局所座標変換部
32 曲率分布生成部
32a 曲率算出部
32b 絶対曲率算出部
33 第1平滑化処理部
34 第2平滑化処理部
g4 ケーブル中心位置(点)
2 Computer 26 Curvature distribution generation computer program 31 Acquisition unit 31a Center position detection unit 31b Local coordinate conversion unit 32 Curvature distribution generation unit 32a Curvature calculation unit 32b Absolute curvature calculation unit 33 First smoothing processing unit 34 Second smoothing processing unit g4 Cable center position (point)
Claims (12)
ケーブルの長手方向における複数の点それぞれの位置情報を取得する取得部と、
前記複数の点の前記位置情報に基づいて前記複数の点の曲率を算出することで、前記ケーブルの長手方向における曲率分布を生成する曲率分布生成部と、
前記曲率分布、及び/又は、前記複数の点の前記位置情報の平滑化処理を行う平滑化処理部と、
を備えているケーブル曲率分布生成システム。 A system for generating a curvature distribution of a cable,
An acquisition unit for acquiring position information of each of a plurality of points in the longitudinal direction of the cable;
A curvature distribution generation unit that generates a curvature distribution in the longitudinal direction of the cable by calculating the curvature of the plurality of points based on the position information of the plurality of points;
A smoothing processing unit that performs a smoothing process on the curvature distribution and / or the position information of the plurality of points;
Cable curvature distribution generation system equipped with.
請求項1記載のケーブル曲率分布生成システム。 The cable curvature distribution generation system according to claim 1, wherein the smoothing processing unit performs smoothing processing of the curvature distribution and / or the position information of the plurality of points for each of a plurality of cable bending directions.
請求項1又は2記載のケーブル曲率分布生成システム。 The cable curvature distribution generation system according to claim 1, wherein the curvature distribution generation unit generates a plurality of direction component-specific curvature distributions indicating a curvature distribution for each of a plurality of cable bending directions.
請求項3記載のケーブル曲率分布生成システム。 The cable curvature distribution generation system according to claim 3, wherein the smoothing processing unit performs a smoothing process on each of the plurality of directional component curvature distributions.
前記複数の点に対応する複数の前記局所座標は、それぞれ、各位置情報が示す点又はその近傍におけるケーブル長手方向が第1の座標軸方向として設定され、他の第2及び第3の座標軸方向が、前記第1の座標軸方向に直交する平面内に設定された3次元座標であり
前記曲率分布生成部は、複数のケーブル曲げ方向ごとの方向成分別曲率分布を生成する
請求項1〜4のいずれか1項に記載のケーブル曲率分布生成システム。 Each of the position information of the plurality of points is indicated as a coordinate value in local coordinates,
In the plurality of local coordinates corresponding to the plurality of points, the cable longitudinal direction at or near the point indicated by each position information is set as the first coordinate axis direction, and the other second and third coordinate axis directions are set. The three-dimensional coordinates set in a plane orthogonal to the first coordinate axis direction, the curvature distribution generation unit generates a curvature distribution by direction component for each of a plurality of cable bending directions. The cable curvature distribution generation system according to claim 1.
前記平滑化処理部は、前記局所座標における前記第2座標軸方向及び第3座標軸方向ごとに、前記変位の平滑化処理を行う
請求項5記載のケーブル曲率分布生成システム。 In the position information, a displacement from an adjacent point of each point in the plurality of points is expressed using the local coordinates,
The cable curvature distribution generation system according to claim 5, wherein the smoothing processing unit performs the displacement smoothing process for each of the second coordinate axis direction and the third coordinate axis direction in the local coordinates.
請求項5又は6記載のケーブル曲率分布生成システム。 The cable curvature distribution generation system according to claim 5 or 6, wherein the curvature distribution generation unit calculates the curvature of each of a plurality of points as a secondary differential value of coordinate values in local coordinates.
前記平滑化処理部は、前記複数の方向成分曲率分布それぞれについて平滑化処理を行う
請求項5〜7のいずれか1項に記載のケーブル曲率分布生成システム。 In the position information, a displacement from an adjacent point of each point in the plurality of points is expressed using the local coordinates,
The cable curvature distribution generation system according to any one of claims 5 to 7, wherein the smoothing processing unit performs a smoothing process on each of the plurality of directional component curvature distributions.
請求項5〜8のいずれか1項に記載のケーブル曲率分布生成システム。 The cable curvature distribution generation system according to any one of claims 5 to 8, wherein the curvature distribution generation unit generates an absolute curvature distribution based on a plurality of directional component-specific curvature distributions that have been smoothed.
請求項1〜9のいずれか1項に記載のケーブル曲率分布生成システム。 The cable curvature according to any one of claims 1 to 9, wherein the smoothing processing unit performs a smoothing process for each section obtained by dividing a longitudinal range of a cable that is a target of curvature distribution generation into a plurality of sections. Distribution generation system.
ケーブルの長手方向における複数の点それぞれの位置情報をコンピュータが取得するステップと、
前記複数の点の前記位置情報に基づいて前記複数の点の曲率を算出することで、前記ケーブルの長手方向における曲率分布を前記コンピュータが生成するステップと、
前記曲率分布、及び/又は、前記複数の点の前記位置情報の平滑化処理を前記コンピュータが行うステップと、
を含む方法。 A method of generating a cable curvature distribution by a computer,
A computer acquiring position information of each of a plurality of points in the longitudinal direction of the cable;
Calculating the curvature distribution in the longitudinal direction of the cable by calculating the curvature of the plurality of points based on the position information of the plurality of points; and
The computer performing a smoothing process on the curvature distribution and / or the position information of the plurality of points;
Including methods.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2010275038A JP5464135B2 (en) | 2010-12-09 | 2010-12-09 | Cable curvature distribution generation system, computer program thereof, and cable curvature distribution generation method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2010275038A JP5464135B2 (en) | 2010-12-09 | 2010-12-09 | Cable curvature distribution generation system, computer program thereof, and cable curvature distribution generation method |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2012122906A true JP2012122906A (en) | 2012-06-28 |
JP5464135B2 JP5464135B2 (en) | 2014-04-09 |
Family
ID=46504483
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2010275038A Active JP5464135B2 (en) | 2010-12-09 | 2010-12-09 | Cable curvature distribution generation system, computer program thereof, and cable curvature distribution generation method |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP5464135B2 (en) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR101402810B1 (en) | 2012-09-17 | 2014-06-03 | 한국전력공사 | Apparatus and method for determining deformation rate of cable in pipe |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH07183349A (en) * | 1993-12-22 | 1995-07-21 | Toshiba Corp | Inspection of semiconductor device |
JPH11178952A (en) * | 1997-12-19 | 1999-07-06 | Mizuno Corp | Method for determining bending point of golf club, golf club shaft with features expressed by the method, and golf club equipped with the golf club shaft |
JP2002500762A (en) * | 1997-05-26 | 2002-01-08 | エステーエヌ アトラス エレクトロニク ゲゼルシャフト ミット ベシュレンクテル ハフツング | Measurement of rail waviness and long undulations by longitudinal light lines. |
JP2005350801A (en) * | 2004-06-09 | 2005-12-22 | Kao Corp | Method for evaluating shape of fiber |
JP2007205802A (en) * | 2006-01-31 | 2007-08-16 | Sumitomo Electric Ind Ltd | Cable shape state measuring method, cable shape state measuring system used in it, cable shape state measuring program, and cable state evaluation method |
JP2009266775A (en) * | 2008-04-30 | 2009-11-12 | Sumitomo Electric Ind Ltd | Bending-state visualizing method, bending-state visualization system, bending-state visualization program, cable-bending life estimation method, and cable-bending life estimation system for conductive wire |
-
2010
- 2010-12-09 JP JP2010275038A patent/JP5464135B2/en active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH07183349A (en) * | 1993-12-22 | 1995-07-21 | Toshiba Corp | Inspection of semiconductor device |
JP2002500762A (en) * | 1997-05-26 | 2002-01-08 | エステーエヌ アトラス エレクトロニク ゲゼルシャフト ミット ベシュレンクテル ハフツング | Measurement of rail waviness and long undulations by longitudinal light lines. |
JPH11178952A (en) * | 1997-12-19 | 1999-07-06 | Mizuno Corp | Method for determining bending point of golf club, golf club shaft with features expressed by the method, and golf club equipped with the golf club shaft |
JP2005350801A (en) * | 2004-06-09 | 2005-12-22 | Kao Corp | Method for evaluating shape of fiber |
JP2007205802A (en) * | 2006-01-31 | 2007-08-16 | Sumitomo Electric Ind Ltd | Cable shape state measuring method, cable shape state measuring system used in it, cable shape state measuring program, and cable state evaluation method |
JP2009266775A (en) * | 2008-04-30 | 2009-11-12 | Sumitomo Electric Ind Ltd | Bending-state visualizing method, bending-state visualization system, bending-state visualization program, cable-bending life estimation method, and cable-bending life estimation system for conductive wire |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR101402810B1 (en) | 2012-09-17 | 2014-06-03 | 한국전력공사 | Apparatus and method for determining deformation rate of cable in pipe |
Also Published As
Publication number | Publication date |
---|---|
JP5464135B2 (en) | 2014-04-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP4285486B2 (en) | Cable shape state measuring method, cable shape state measuring system used therefor, cable shape state measuring program, cable state evaluating method | |
US20110262015A1 (en) | Image processing apparatus, image processing method, and storage medium | |
WO2012114574A1 (en) | Image magnification device and method | |
EP3300025B1 (en) | Image processing device and image processing method | |
CN112017233A (en) | Method and device for measuring conical appearance of reaction force, computer equipment and system | |
JP5405766B2 (en) | Bending state visualization method of conductor wire, method of outputting bending state, bending state visualization system, bending state visualization program, cable bending life prediction method, and cable bending life prediction system | |
JP4535061B2 (en) | Method for determining arrangement period of cable strands, cable shape state measuring method, cable shape state measuring system used therefor, cable shape state measuring program | |
JP5464135B2 (en) | Cable curvature distribution generation system, computer program thereof, and cable curvature distribution generation method | |
JP5919949B2 (en) | Cable analysis system, cable analysis method, and computer program | |
Ye et al. | A Method of Binocular Laser 3-D Scanning Imaging for Reflective Workpieces | |
US11055865B2 (en) | Image acquisition device and method of operating image acquisition device | |
KR102333768B1 (en) | Hand recognition augmented reality-intraction apparatus and method | |
CN116518868A (en) | Deformation measurement method, device, equipment and storage medium based on artificial intelligence | |
CN111079523A (en) | Object detection method, object detection device, computer equipment and storage medium | |
JP5620741B2 (en) | Information processing apparatus, information processing method, and program | |
KR101469851B1 (en) | Multi-objective optimization-based approach for human posture prediction | |
JP5107413B2 (en) | Moving distance calculation device and moving distance calculation method | |
Odehnal | Subdivision algorithms for ruled surfaces | |
CN111156948A (en) | Three-dimensional data coordinate correction method and device for 3C glass panel detection | |
JP2013178113A (en) | Apparatus, method, and program for restoring three-dimensional shape of subject | |
JPH11339052A (en) | Device and method for analyzing shape | |
JP6448459B2 (en) | Rotation state calculation device | |
CN115115700B (en) | Object attitude estimation method and device, electronic equipment and storage medium | |
JP6448458B2 (en) | Rotation state calculation device | |
CN113808050B (en) | Denoising method, device and equipment for 3D point cloud and storage medium |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20130625 |
|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20131008 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20131015 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20131203 |
|
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: 20131224 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20140106 |
|
R150 | Certificate of patent or registration of utility model |
Ref document number: 5464135 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
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 |
|
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 |