Ultrasound matrix inspection

0001This invention relates to methods and devices for carrying out ultrasound inspection, and for pipe inspections.

0002<patcit id="pcit0001" dnum="US20110087444"><text>US App. Pub. 2011/0087444 to Volker</text></patcit> (hereinafter the `444 publication) is directed to a "pig" for crawling through the bore of a pipe and performing ultrasound inspection of the inner pipe surface. The reference discloses an algorithm for imaging the pipe surface based on backscatter signals. The '444 publication involves Fermat's principle to determine sound paths with the shortest travel time. The modeling involves first building a grid and determining travel time for each point in the grid. The `444 reference requires scanning a pipe from the inside, where the primary information to be ascertained is 3D information about the inner surface of the pipe. This does not solve than the problem of accurately modeling the inner surface of a pipe using a scanning apparatus positioned on the outer surface.

0003<patcit id="pcit0002" dnum="US7685878B"><text>US 7,685,878 to Brandstrom</text></patcit> (hereinafter the `878 patent) relates to a device for rotating a pair of ultrasound transducers around a pipe circumference for pipe weld inspection. It allows the cables and other apparatus extending away from the transducers to remain stationary, extending away in only a single direction. `878 teaches an apparatus which can be mounted on the pipe at the position adjacent the weld and which carries the transducers and rotates those transducers around the pipe, bearing in mind that effective access to the pipe is generally only available from one side of the pipe.

0004Two transducers are rotated around a circumferential location on a cylindrical body for structural testing of the body, carried on a mounting and drive apparatus including a magnetic attachment which can be manually brought up to a pipe from one side only for fixed connection to the pipe on that side at a position axially spaced from a weld. A collar shaped support for the pair of transducers is formed of a row of separate segments which wrap around the pipe from the one side and is rotated around the axis of the pipe to carry the transducer around the circumferential weld. The segments carry rollers to roll on the surface and are held against the pipe by magnets. The transducers are carried on the support in fixed angular position to track their position but in a manner which allows slight axial or radial movement relative to the pipe.

0005<patcit id="pcit0003" dnum="US7412890B"><text>US 7,412,890 to Johnson </text></patcit>(hereinafter the '890 patent) relates to a method and apparatus for detecting cracks in pipe welds comprising flooding a volume adjacent to the outer pipe surface with water, then using phased array ultrasound to scan the pipe surface. The apparatus has a rectangular cavity that has its open bottom surface pressed against the pipe surface and is flooded with water. The ultrasound array is positioned at the top of the cavity. Phased-array data collection methods are used.

0006<patcit id="pcit0004" dnum="US5515298A"><text>US 5,515,298 to Bicz </text></patcit>(hereinafter the `298 Patent) relates to an apparatus for performing ultrasound scanning of a fingerprint or other object placed on a concave surface. The apparatus projects ultrasound from an array of transducers through an array of pinholes (one per transducer) and against the concave interior of the surface on which the fingerprint rests. The transducers then derive characteristics of the fingerprint from the reflection and scattering of the spherical waveform produced by the pinhole. The apparatus appears to depend on the known structure of the convexo-concave lens structure of the support on which the fingerprint rests.

0007<patcit id="pcit0005" dnum="US6896171B"><text>US 6,896,171 to Den Boer et al</text></patcit> (hereinafter the '171 Patent) relates to an apparatus for performing EMAT (electromagnetic acoustic transducer) scanning of a freshly-made pipe weld while still hot. The apparatus may include an array of EMAT transmitter and receiver coils positioned on a ring structure around the outer surface of the pipe. No post-processing algorithm details are disclosed. The apparatus is described as being able to detect the presence of weld defects, and gives some information as to their size, but neither images, precise locations, nor are any further details of defects discussed in the description.

0008<patcit id="pcit0006" dnum="US20090158850"><text>US App. Pub. No. 2009/0158850 to Alleyne et al</text></patcit> (hereinafter the '850 publication) relates to a method and apparatus for inspecting pipes wherein the pig apparatus is inserted into the bore of the pipe. Ultrasound transducers are pressed against the inner walls of the pipe and use guided waves (e.g. Lamb waves) of ultrasound within the material of the pipe wall itself to detect defects. Data collection and processing appears to be based on a full matrix capture technique from which different wave modes may be extracted, although a phased-array data collection technique may also be used.

0009<patcit id="pcit0007" dnum="US20090078742"><text>US App. Pub. No. 2009/0078742 to Pasquali et al.</text></patcit> (hereinafter the '742 publication) relates to a method and apparatus for inspecting multi-walled pipes, such as those used for undersea transport of hot or cold fluids. The method involves placing an ultrasound probe against the inner pipe surface and scanning at various intervals as the probe rotates around the inner circumference of the pipe wall. The apparatus is a probe positioned at the end of a rotatable arm, which positions the probe within the pipe and then rotates it about the circumference of the inner wall. The `742 publication also discloses methods of positioning the probe at various angles relative to the pipe surface. However, it appears to only teach the use of probes that are displaced from the weld in the pipe's axial direction, and are angled forward or backward toward the location of the pipe weld.

0010Additional prior art references include <patcit id="pcit0008" dnum="US7762136B"><text>US 7762136 to Ume, Ifeanyi C. et al.</text></patcit>, which teaches ultrasound systems and methods for measuring weld penetration depth in real time and off line, <patcit id="pcit0009" dnum="US7694569B"><text>US 7694569 to McGrath, Matthew et al.</text></patcit> which teaches a phased array ultrasonic water wedge apparatus, <patcit id="pcit0010" dnum="US7694564B"><text>US 7694564 to Brignac, Jacques L. et al.</text></patcit> which teaches a boiler tube inspection probe with centering mechanism and method of operating the same, <patcit id="pcit0011" dnum="US6935178B"><text>US 6935178 to Prause, Reinhard</text></patcit> which teaches a device for inspecting pipes using ultrasound, <patcit id="pcit0012" dnum="US6734604B"><text>US 6734604 to Butler, John V. et al.</text></patcit> which teaches a multimode synthesized beam transduction apparatus, <patcit id="pcit0013" dnum="US4872130A"><text>US 4872130 to Pagano, Dominick A.</text></patcit>, which teaches an automated in-line pipe inspection system <patcit id="pcit0014" dnum="JP2004028937B"><text>JP 2004028937 to Furukawa, T. et al.</text></patcit>, which teaches a method for measuring the shape of a welded pipe. <patcit id="pcit0015" dnum="US6409669B1"><text>US 6409669 B1</text></patcit> discloses an ultrasound transducer assembly including an acoustic mirror, an ultrasound transducer positioned to direct a scanned ultrasound at the acoustic mirror, and an actuating device.

0011Example embodiments described in this document relate to methods and devices for performing ultrasound inspection of objects using full matrix data capture techniques.

0012In a first aspect, the application is directed to a device for performing ultrasound scanning of a conduit, comprising a body having a cuff adapted to fit on the object, a carrier adapted to traverse a circumference of the cuff, an ultrasound probe mounted on the carrier and positioned to scan the object, a carrier motor coupled to the carrier for driving the carrier around the circumference of the cuff, a data processing system/controller configured to provide control information, one or more data connections providing control information for the carrier motor and the ultrasound probe and receiving scanning data from the ultrasound probe, an adjustable reflector mounted to the carrier; and a reflector motor for controlling an angle of the adjustable reflector in a plane substantially normal to a longitudinal axis of the object, wherein the ultrasound probe is positioned to scan the object via reflection of ultrasound signals off of the adjustable reflector and the one or more data connections provide control information for the reflector motor.

0013In another aspect, the cuff forms a liquid-resistant seal around the circumference of the conduit, and the device further comprises a liquid feed for receiving a liquid scanning medium and filling the volume defined between the interior of the cuff and the exterior of the conduit with the liquid scanning medium

0014In a further aspect, the device further comprises a power connection for receiving electrical power for the carrier motor.

0015In a further aspect, the cuff is configurable between an open configuration allowing it to be fitted around the conduit and a closed configuration encircling the conduit.

0016In a further aspect, the device further comprises an adjustable reflector mounted to the carrier and a reflector motor for controlling an angle of the adjustable reflector in a plane substantially normal to a longitudinal axis of the object, wherein the ultrasound probe is positioned to scan the object via reflection of ultrasound signals off of the adjustable reflector, and the one or more data connections provide control information for the reflector motor.

0017In a further aspect, the device further comprises a power connection for receiving electrical power for the reflector motor.

0018In a further aspect, the conduit is a cylinder.

0019In a further aspect, the ultrasound probe is an array of ultrasound transceivers.

0020In a further aspect, the cuff comprises a knuckle which releasably secures a first half of said cuff to a second half of said cuff.

0021In a further aspect, the cuff comprises a first half of said cuff detachable from a second half of said cuff.

0022In a further aspect, the application is directed to a method for performing ultrasound scanning of a conduit, comprising providing an ultrasound array having a plurality of ultrasound elements arrayed substantially parallel to a longitudinal axis of the conduit, positioning the ultrasound array to project ultrasound signals toward an external surface of the object at a first point about the circumference of the conduit, performing a full-matrix-capture scan of the first point about the circumference of the conduit, comprising: transmitting an ultrasound signal from a first ultrasound element in the ultrasound array; sensing and recording ultrasound signals received by each other ultrasound element in the ultrasound array; and repeating the steps of transmitting, sensing and recording, wherein the step of transmitting is performed in turn by each ultrasound element in the ultrasound array other than the first ultrasound element; repositioning the ultrasound array at a second point about the circumference of the conduit, performing a full-matrix-capture scan of the second point about the circumference of the conduit, and repeating the steps of repositioning and performing a full-matrix-capture scan.

0023In a further aspect, the method further comprises, before performing each full-matrix-capture scan, transmitting at least one ultrasound signal from at least one ultrasound element in the ultrasound array, sensing at least one ultrasound signal received by at least one ultrasound element in the ultrasound array, evaluating at a processor the quality of the at least one sensed signal, and adjusting a scanning angle of the ultrasound array based on the outcome of the evaluation.

0024in a further aspect, the ultrasound array projects ultrasound signals toward the external surface of the object by reflecting the ultrasound signals off of an adjustable reflector, and adjusting the scanning angle of the ultrasound array comprises adjusting the angle of the adjustable reflector.

0025in a further aspect, the application is directed to a method of modeling the near and far surfaces of an object within a scanning plane passing through the near and far surfaces of the object, comprising providing a set of full-matrix-capture ultrasound scanning data corresponding to a scanning area within the scanning plane, the full-matrix-capture ultrasound scanning data captured using an ultrasound array transmitting and sensing ultrasound signals through a scanning medium situated between the ultrasound array and the near surface of the object and performing the steps of: transmitting an ultrasound signal from a first ultrasound element in the ultrasound array; sensing and recording ultrasound signals received by each other ultrasound element in the ultrasound array; and repeating the steps of transmitting, sensing and recording, wherein the step of transmitting is performed by each ultrasound element in the ultrasound array other than the first ultrasound element;constructing a first intensity map of the scanning area, comprising a plurality of points within the scanning area having associated intensity values, by calculating travel times of ultrasound signals through the scanning medium based on the full-matrix-capture ultrasound scanning data; filtering the first intensity map to model the boundary of the near surface within the scanning area; using the modeled boundary of the near surface as a lens in constructing a second intensity map, comprising a plurality of points within the scanning area having associated intensity values, by the application of Fermat's Principle, to compute ultrasound signal travel times through both the scanning medium and the object based on the full-matrix-capture ultrasound scanning data; and filtering the second intensity map to model the boundary of the far surface within the scanning area.

0026In a further aspect, constructing a first intensity map of the scanning area comprises calculating an intensity <i>I</i> at a plurality of points <i>r</i> within the scanning area where <i>I</i> is defined as the sum of the amplitude of the data-set of analytic time-domain signals from ultrasound array transmitter element <i>i</i> to ultrasound array receiver element <i>j</i> at time <i>t</i> for all <i>i</i> and <i>j,</i> where <i>t</i> is defined for each <i>i,j</i> pair as being the time it takes for sound to travel through the scanning medium.

0027In a further aspect, constructing a first intensity map of the scanning area comprises calculating an intensity I at a plurality of points <i>r</i> within the scanning area defined by the equation <maths id="math0001"><img file="EP3553561B1_D0001.tif" /></maths><maths id="math0002"><img file="EP3553561B1_D0002.tif" /></maths><maths id="math0003"><math display="block"><mi>I</mi><mfenced><mi>r</mi></mfenced><mo>=</mo><mfenced open="|" close="|"><mstyle displaystyle="true"><msub><mo>∑</mo><mrow><mi>i</mi><mo>,</mo><mi>j</mi></mrow></msub><mrow><msub><mi>g</mi><mrow><mfenced><mi>i</mi></mfenced><mi>j</mi></mrow></msub><mfenced><mi>t</mi><mo>=</mo><mfrac><mrow><mfenced open="|" close="|"><msub><mi mathvariant="bold">e</mi><mfenced><mi>i</mi></mfenced></msub><mo>−</mo><mi mathvariant="bold">r</mi></mfenced><mo>+</mo><mfenced open="|" close="|"><msub><mi mathvariant="bold">e</mi><mi>j</mi></msub><mo>−</mo><mi mathvariant="bold">r</mi></mfenced></mrow><mi>c</mi></mfrac></mfenced></mrow></mstyle></mfenced></math><img file="EP3553561B1_D0003.tif" /></maths> wherein <i>g_(i)j(t)</i> is the amplitude of the data-set of analytic time-domain signals from ultrasound array transmitter element <i>i</i> to ultrasound array receiver element <i>j</i> at time <i>t,</i><b>r</b> is the vector defining point <i>r</i> relative to a coordinate origin, e_(i) is a vector defining the position of ultrasound array transmitter element <i>i</i> relative to the coordinate origin, e_j is a vector defining the position of ultrasound array receiver element <i>j</i> relative to the coordinate origin, and <i>c</i> is the speed of sound traveling through the scanning medium.

0028In a further aspect, constructing a first intensity map of the scanning area comprises calculating an intensity <i>I</i> at a plurality of points <i>r</i> within the scanning area, each point intensity being calculated at a plurality of apertures defined by a fixed plurality of ultrasound array elements and the highest intensity of point <i>r</i> calculated for a single aperture being used to represent the intensity of point <i>r</i> in the intensity map.

0029In a further aspect, constructing a first intensity map of the scanning area comprises calculating an intensity <i>I</i> at a plurality of points <i>r</i> within the scanning area defined by the equation <maths id="math0004"><img file="EP3553561B1_D0004.tif" /></maths><maths id="math0005"><img file="EP3553561B1_D0005.tif" /></maths><maths id="math0006"><math display="block"><mi>I</mi><mfenced><mi>r</mi></mfenced><mo>=</mo><munder><mi>max</mi><mrow><mi>a</mi><mo>∈</mo><mi>A</mi></mrow></munder><mfenced open="{" close="}"><mi>I</mi><mfenced><mi mathvariant="bold">r</mi><mi>a</mi></mfenced></mfenced></math><img file="EP3553561B1_D0006.tif" /></maths> wherein <maths id="math0007"><math display="block"><mi>I</mi><mfenced><mi>r</mi><mi>a</mi></mfenced><mo>=</mo><mfenced open="|" close="|"><mstyle displaystyle="true"><msub><mo>∑</mo><mrow><mi>i</mi><mo>,</mo><mi>j</mi><mo>∈</mo><mi>a</mi></mrow></msub><mrow><msub><mi>g</mi><mrow><mfenced><mi>i</mi></mfenced><mi>j</mi></mrow></msub><mfenced><mi>t</mi><mo>=</mo><mfrac><mrow><mfenced open="|" close="|"><msub><mi mathvariant="bold">e</mi><mfenced><mi>i</mi></mfenced></msub><mo>−</mo><mi mathvariant="bold">r</mi></mfenced><mo>+</mo><mfenced open="|" close="|"><msub><mi mathvariant="bold">e</mi><mi>j</mi></msub><mo>−</mo><mi mathvariant="bold">r</mi></mfenced></mrow><mi>c</mi></mfrac></mfenced></mrow></mstyle></mfenced></math><img file="EP3553561B1_D0007.tif" /></maths> and wherein <i>g_(i)j(t)</i> is the amplitude of the data-set of analytic time-domain signals from ultrasound array transmitter element <i>i</i> to ultrasound array receiver element <i>j</i> at time <i>t,</i><b>r</b> is the vector defining point <i>r</i> relative to a coordinate origin, e_(i) is a vector defining the position of ultrasound array transmitter element <i>i</i> relative to the coordinate origin, e_j is a vector defining the position of ultrasound array receiver element <i>j</i> relative to the coordinate origin, <i>c</i> is the speed of sound traveling through the scanning medium, <i>a</i> is an aperture defined by a fixed plurality of adjacent ultrasound elements in the ultrasound array, and <i>A</i> is a set comprising a plurality of such apertures.

0030In a further aspect, using the modeled boundary of the near surface as a lens in constructing a second intensity map comprises calculating an intensity <i>I</i> at a plurality of points <i>r</i> within the scanning area defined by the equation <maths id="math0008"><img file="EP3553561B1_D0008.tif" /></maths><maths id="math0009"><img file="EP3553561B1_D0009.tif" /></maths><maths id="math0010"><math display="block"><mi>I</mi><mfenced><mi>r</mi></mfenced><mo>=</mo><munder><mi>max</mi><mrow><mi>a</mi><mo>∈</mo><mi>A</mi></mrow></munder><mfenced open="{" close="}"><mi>I</mi><mfenced><mi mathvariant="bold">r</mi><mi>a</mi></mfenced></mfenced></math><img file="EP3553561B1_D0010.tif" /></maths> wherein <maths id="math0011"><math display="block"><mi>I</mi><mfenced><mi mathvariant="bold">r</mi><mi>a</mi></mfenced><mo>=</mo><mfenced open="|" close="|"><mstyle displaystyle="true"><munder><mo>∑</mo><mrow><mi>i</mi><mo>,</mo><mi>j</mi><mo>∈</mo><mi>a</mi></mrow></munder><mstyle displaystyle="true"><munder><mo>∑</mo><mrow><mi>t</mi><mo>′</mo><mo>∈</mo><msubsup><mi>T</mi><mi mathvariant="italic">ij</mi><mi>K</mi></msubsup><mfenced><mi mathvariant="bold">r</mi></mfenced></mrow></munder><mrow><msub><mi>g</mi><mrow><mfenced><mi>i</mi></mfenced><mi>j</mi></mrow></msub><mfenced><mi>t</mi><mo>′</mo></mfenced></mrow></mstyle></mstyle></mfenced></math><img file="EP3553561B1_D0011.tif" /></maths> and wherein <maths id="math0012"><math display="block"><msubsup><mi>T</mi><mi mathvariant="italic">ir</mi><mi>K</mi></msubsup><mfenced><mi mathvariant="bold-italic">r</mi></mfenced><mo>=</mo><mfenced open="{" close="}"><msub><mi>t</mi><mi mathvariant="italic">ir</mi></msub><mo>+</mo><mrow><msub><mi>t</mi><mi mathvariant="italic">jr</mi></msub><mo>|</mo></mrow><msub><mi>t</mi><mi mathvariant="italic">ir</mi></msub><mo>∈</mo><msubsup><mi>T</mi><mi mathvariant="italic">ir</mi><mi>K</mi></msubsup><mo>,</mo><msub><mi>t</mi><mi mathvariant="italic">jr</mi></msub><mo>∈</mo><msubsup><mi>T</mi><mi mathvariant="italic">jr</mi><mi>K</mi></msubsup></mfenced></math><img file="EP3553561B1_D0012.tif" /></maths> and wherein <maths id="math0013"><math display="inline"><msubsup><mi>T</mi><mi mathvariant="italic">ir</mi><mi>K</mi></msubsup></math><img file="EP3553561B1_D0013.tif" /></maths> is the set of all first-order and multiple-order derivatives of the time it takes sound to travel from ultrasound array transmitter element <i>i</i> to a point K on the boundary of the near surface, <maths id="math0014"><math display="inline"><msubsup><mi>T</mi><mi mathvariant="italic">jr</mi><mi>K</mi></msubsup></math><img file="EP3553561B1_D0014.tif" /></maths> is the set of all first-order and multiple-order derivatives of the time it takes sound to travel from ultrasound array receiver element <i>j</i> to a point K on the boundary of the near surface, <i>g_(i)j(t)</i> is the amplitude of the data-set of analytic time-domain signals from ultrasound array transmitter element <i>i</i> to ultrasound array receiver element <i>j</i> at time <i>t,</i><b>r</b> is the vector defining point <i>r</i> relative to a coordinate origin, e_(i) is a vector defining the position of ultrasound array transmitter element <i>i</i> relative to the coordinate origin, e_j is a vector defining the position of ultrasound array receiver element <i>j</i> relative to the coordinate origin, <i>c</i> is the speed of sound traveling through the scanning medium, <i>a</i> is an aperture defined by a fixed plurality of adjacent ultrasound elements in the ultrasound array, and <i>A</i> is a set comprising a plurality of such apertures.

0031In a further aspect, the method further comprises, before constructing a first intensity map, filtering the full-matrix-capture ultrasound scanning data to remove noise.

0032In a further aspect, filtering the first intensity map comprises passing the intensity map through an edge-detection filter and using the output as a model of the boundary of the near surface within the scanning area, and filtering the second intensity map comprises passing the intensity map through an edge-detection filter and using the output as a model of the boundary of the far surface within the scanning area.

0033In a further aspect, filtering the first intensity map and filtering the second intensity map each further comprise dilation of the detected edges produced by the edge-detection filter.

0034In a further aspect, filtering the first intensity map and filtering the second intensity map each further comprise thinning the dilated edges.

0035In a further aspect, filtering the first intensity map and filtering the second intensity map each further comprise selecting a single component from each vertical slice of the intensity map and removing all other components in that slice in order to maximize the continuity and length of the remaining components.

0036In a further aspect, the application is directed to a method of modeling the near and far surfaces of an object, comprising applying the methods above to a plurality of sets of full-matrix-capture ultrasound scanning data corresponding to a plurality of scanning planes passing through the near and far surfaces of the object, and modeling the near and far surfaces of the object based on the modeled boundaries within each scanning plane and the relative locations of each scanning plane.

0037In a further aspect, the plurality of scanning planes are parallel to and adjacent to each other.

0038In a further aspect, the object is substantially cylindrical, and the plurality of scanning planes all pass through the longitudinal axis of the object.

0039In a further aspect, the application is directed to a device for performing ultrasound scanning of an object, comprising a body adapted to fit on the object; an ultrasound probe mounted on the body and positioned to scan the body; one or more data connections providing control information for the carrier motor and the ultrasound probe and receiving scanning data from the ultrasound probe; an adjustable reflector mounted to the carrier; and a reflector motor for controlling an angle of the adjustable reflector in a plane substantially normal to a longitudinal axis of the object, wherein: the ultrasound probe is positioned to scan the object via reflection of ultrasound signals off of the adjustable reflector; and the one or more data connections provide control information for the reflector motor.

0040In a further aspect, the body forms a liquid-resistant seal around the circumference of the object, and the device further comprises a liquid feed for receiving a liquid scanning medium and filling the volume defined between the interior of the body and the exterior of the object with the liquid scanning medium.

0041In a further aspect, the device further comprises a power connection for receiving electrical power for the carrier motor.

0042In a further aspect, the device further comprises a power connection for receiving electrical power for the reflector motor.

0043In a further aspect, the ultrasound probe is an array of ultrasound transceivers.

0044Other example embodiments of the present disclosure will be apparent to those of ordinary skill in the art from a review of the following detailed description in conjunction with the drawings. <ul id="ul0001" list-style="none" compact="compact"><li><figref idref="f0001">Figure 1</figref> is a perspective view of an ultrasound probe manipulator according to an example embodiment, operationally fitted to a pipe;</li><li><figref idref="f0001">Figure 2</figref> is a perspective view of the example probe manipulator of <figref idref="f0001">Figure 1</figref>;</li><li><figref idref="f0002">Figure 3</figref> is a side view of the example probe manipulator of <figref idref="f0001">Figure 1</figref> operationally fitted to a pipe;</li><li><figref idref="f0002">Figure 4</figref> is an isometric view of the example probe manipulator of <figref idref="f0001">Figure 1</figref> operationally fitted to a pipe;</li><li><figref idref="f0003">Figure 5</figref> is a perspective view of a hinged probe manipulator according to an example embodiment, showing the manipulator in an open configuration;</li><li><figref idref="f0003">Figure 6</figref> is an isometric view of a hinged probe manipulator according to an example embodiment, showing the manipulator disassembled into two halves and a hinged knuckle;</li><li><figref idref="f0004">Figure 7</figref> is a side view of a linear multi-element ultrasound probe array;</li><li><figref idref="f0004">Figure 8</figref> is a diagram of an N by N matrix of ultrasound scan data generated by the Full Matrix Capture technique using an N-element ultrasound probe array;</li><li><figref idref="f0004">Figure 9</figref> is an isometric representation of an N by N by M matrix of A-scan data as in <figref idref="f0004">Figure 8</figref>, showing the third dimension of M samples recorded over time for each A-scan;</li><li><figref idref="f0005">Figure 10</figref> is a side view of a linear ultrasound probe array scanning a volume through an intermediate medium;</li><li><figref idref="f0006">Figure 11</figref> is a graph of the curve T_ir(K) from point u₁ to u₂ on the surface of the object being scanned;</li><li><figref idref="f0007">Figure 12</figref> is an intensity map in the scanning plane of the outer surface of a pipe wall;</li><li><figref idref="f0007">Figure 13</figref> is an intensity map in the scanning plane of the inner surface of a pipe wall;</li><li><figref idref="f0007">Figure 14</figref> is a graph of the outer surface and inner surface of a pipe wall derived from intensity maps in the scanning plane;</li><li><figref idref="f0008">Figure 15</figref> is a graph of an the real part and the envelope of an A-scan signal before wave packet normalization;</li><li><figref idref="f0008">Figure 16</figref> is a graph of an the real part and the envelope of an A-scan signal after wave packet normalization;</li><li><figref idref="f0009">Figure 17</figref> is a flow chart showing the steps carried out in modeling of outer and inner surfaces of a scanned object according to an example embodiment;</li><li><figref idref="f0010">Figure 18</figref> is a flow chart showing the steps involved in pre-processing of captured ultrasound data according to an example embodiment;</li><li><figref idref="f0011">Figure 19</figref> is a cross-sectional view of an X-Z scanning plane through a pipe wall showing how an intensity map is built in an example embodiment;</li><li><figref idref="f0012">Figure 20</figref> is a flow chart showing the steps involved in building an intensity map using the Shifting Aperture Focus Method according to an example embodiment;</li><li><figref idref="f0013">Figure 21</figref> is a flow chart showing the steps involved in boundary recognition of an intensity map according to an example embodiment;</li><li><figref idref="f0014">Figure 22</figref> is an example black-and-white image used in conjunction with <figref idref="f0014">Figures 23 through 25</figref> to illustrate the effect of erosion and dilation operations;</li><li><figref idref="f0014">Figure 23</figref> is an example 5 x 5 matrix, centered around the origin, used as the structuring element in conjunction with <figref idref="f0014">Figures 22, 24 and 25</figref> to illustrate the effect of erosion and dilation operations;</li><li><figref idref="f0014">Figure 24</figref> is a black-and-white image showing the effect of an erosion operation on <figref idref="f0014">Figure 22</figref> using the structuring element of <figref idref="f0014">Figure 23</figref>;</li><li><figref idref="f0014">Figure 25</figref> is a black-and-white image showing the effect of a dilation operation on <figref idref="f0014">Figure 22</figref> using the structuring element of <figref idref="f0014">Figure 23</figref>;</li><li><figref idref="f0015">Figure 26(a)</figref> is an example black-and-white image used in conjunction with <figref idref="f0015">Figure 26(b)</figref> to illustrate the effect of canny edge detection;</li><li><figref idref="f0015">Figure 26(b)</figref> is a black-and-white image showing the effect of canny edge detection on <figref idref="f0015">Figure 26(a)</figref>;</li><li><figref idref="f0015">Figure 27(a)</figref> is an example black-and-white image used in conjunction with <figref idref="f0015">Figure 27(b)</figref> to illustrate the effect of a thinning algorithm;</li><li><figref idref="f0015">Figure 27(b)</figref> is a black-and-white image showing the effect of a thinning algorithm on <figref idref="f0015">Figure 27(a)</figref>;</li><li><figref idref="f0016">Figure 28(a)</figref> is an example black-and-white image showing a junction;</li><li><figref idref="f0016">Figure 28(b)</figref> is an enlarged view of the junction of <figref idref="f0016">Figure 28(a)</figref>, showing the rectangular path traced around the junction to count light-dark cycles for junction detection;</li><li><figref idref="f0016">Figure 29</figref> is a graph of an example A-scan time-domain signal showing the travel time from transmitter excitation to the leading edge of a received wave packet;</li><li><figref idref="f0017">Figure 30</figref> is an example intensity map of a scanned area;</li><li><figref idref="f0017">Figure 31</figref> is a plot of the edges of the intensity map <figref idref="f0017">Figure 30</figref> plotting the pixels of maximum intensity in the vertical direction;</li><li><figref idref="f0017">Figure 32</figref> is a plot of boundary edges near and above the high intensity pixels of <figref idref="f0017">Figure 31</figref>;</li><li><figref idref="f0018">Figure 33</figref> is a plot of an example output of a boundary detection algorithm without error correction;</li><li><figref idref="f0018">Figure 34</figref> is a plot of the boundary of <figref idref="f0018">Figure 33</figref> after a dilation operation is performed with a rectangular structuring element;</li><li><figref idref="f0018">Figure 35</figref> is a plot of the boundary of <figref idref="f0018">Figure 34</figref> with a thinning algorithm applied to it;</li><li><figref idref="f0018">Figure 36</figref> is a plot of the boundary of <figref idref="f0018">Figure 35</figref> with erroneous pixels removed;</li><li><figref idref="f0019">Figure 37</figref> is an intensity map image of an inner pipe surface (ID);</li><li><figref idref="f0019">Figure 38</figref> is the intensity map of <figref idref="f0019">Figure 37</figref> after applying edge detection, dilation and thinning algorithms, with a junction area of the boundary circled;</li><li><figref idref="f0019">Figure 39</figref> is an enlarged view of the circled junction area of <figref idref="f0019">Figure 38</figref>;</li><li><figref idref="f0019">Figure 40</figref> is the enlarged view of <figref idref="f0019">Figure 39</figref> with the junction removed;</li><li><figref idref="f0020">Figure 41</figref> is the enlarged view of <figref idref="f0019">Figure 40</figref> with bottom portions of connected components removed;</li><li><figref idref="f0020">Figure 42</figref> is the enlarged view of <figref idref="f0020">Figure 41</figref> with small connected components removed;</li><li><figref idref="f0020">Figure 43</figref> is a pipe surface intensity map showing the edge output from a Canny edge detector and a graph of the maximum intensity pixels approximating the shape of the true boundary as a flat plate;</li><li><figref idref="f0021">Figure 44</figref> is an pipe surface intensity map showing interpolation of the boundary with areas between two boundary segments extracted from the edge detection process connected with straight lines;</li><li><figref idref="f0022">Figure 45</figref> is a flow chart of the sequential operation of an example algorithm for boundary recognition and definition;</li><li><figref idref="f0023">Figure 46</figref> is a side view of an ultrasound probe carrier showing an adjustable reflector, for use with the example manipulator of <figref idref="f0001 f0002 f0003">Figures 1 to 6</figref>;</li><li><figref idref="f0023">Figure 47</figref> is a side view of the ultrasound probe carrier of <figref idref="f0023">Figure 46</figref>, shown from the opposite side;</li><li><figref idref="f0024">Figure 48</figref> is a cross-sectional side view of a pipe wall showing uneven thinning;</li><li><figref idref="f0024">Figure 49</figref> is a cross-sectional side view of a pipe wall at a fitting to fitting weld showing a zone inspected for thickness;</li><li><figref idref="f0024">Figure 50</figref> is a cross-sectional side view of a pipe wall at a Grayloc^™ fitting to fitting weld showing a zone inspected for thickness;</li><li><figref idref="f0025">Figure 51</figref> is a network diagram showing remote data acquisition and analysis computers according to an example embodiment;</li><li><figref idref="f0026">Figure 52</figref> is a network diagram showing the relation of local hardware components used at the inspection site according to an example embodiment;</li><li><figref idref="f0027">Figure 53</figref> is a diagram of the equipment configuration at the remote acquisition site according to an example embodiment;</li><li><figref idref="f0028">Figure 54</figref> is an example Distance Amplitude Curve (DAC);</li><li><figref idref="f0029">Figure 55</figref> is a 2" reference block specimen used for calibrating a 2" manipulator according to an example embodiment;</li><li><figref idref="f0030">Figure 56</figref> is a 2" reference block specimen used for calibrating a 2" manipulator in a 6 degree configuration according to an example embodiment;</li><li><figref idref="f0031">Figure 57</figref> is a 2.5" reference block specimen used for calibrating a 2.5" manipulator according to an example embodiment;</li><li><figref idref="f0032">Figure 58</figref> is a perspective view of a calibration block setup showing an example manipulator in an open configuration prior to calibration;</li><li><figref idref="f0033">Figure 59</figref> is a flow chart showing the steps performed during calibration of an example manipulator using an example calibration block according to an example embodiment;</li><li><figref idref="f0034">Figure 60</figref> is a flow chart showing the steps performed during inspection of a pipe diameter according to an example embodiment;</li><li><figref idref="f0035">Figure 61</figref> is a flow chart of steps performed by a data analyst in the process of inspection data analysis according to an example embodiment;</li><li><figref idref="f0036">Figure 62</figref> is a flow chart of steps performed by a data analyst and the various steps applied by the analysis algorithms on the gateway server, blade server, and local analysis PC in the process of inspection data analysis according to an example embodiment;</li><li><figref idref="f0037">Figure 63</figref> is an example Main Results window in the NEOVISION^™ application used by a data analyst in some embodiments;</li><li><figref idref="f0038">Figure 64</figref> is an example 3D window in the NEOVISION^™ application used by a data analyst in some embodiments;</li><li><figref idref="f0039">Figure 65</figref> is an example 3D Overview window in the NEOVISION^™ application used by a data analyst in some embodiments;</li><li><figref idref="f0040">Figure 66</figref> is an example 3D Pop-up window in the NEOVISION^™ application used by a data analyst in some embodiments;</li><li><figref idref="f0041">Figure 67</figref> is an example data flow for IFM OD boundary preparation;</li><li><figref idref="f0042">Figure 68</figref> is a detail view of the samples of interest in a time-domain ultrasound signal showing the time period during which data samples are to be retained according to an example embodiment;</li><li><figref idref="f0043">Figure 69</figref> is a flow chart showing the entire boundary recognition data flow according to an example embodiment;</li><li><figref idref="f0044">Figure 70</figref> is an example data flow for boundary definition;</li><li><figref idref="f0045">Figure 71</figref> is a data flow chart for the Interior Focusing Method process;</li><li><figref idref="f0046">Figure 72</figref> is a side cross-sectional view through an example pipe weld showing example definitions of thickness;</li><li><figref idref="f0046">Figure 73</figref> is a diagram showing the relationship of Cartesian and spherical coordinate systems mapped onto an example linear ultrasound array having length L;</li><li><figref idref="f0047">Figure 74</figref> is a graph of ultrasound beam directivity showing amplitude increasing from 0 to 1 on the vertical axis and angle from -90 degrees to 90 degrees on the horizontal axis for an example embodiment having <i>a</i> = 0.23 mm, <i>f</i>= 7.5 MHz, and <i>c</i> = 1480 m/s;</li><li><figref idref="f0047">Figure 75</figref> is a spectral graph of element directivity showing angle decreasing from 90 to -90 on the vertical axis and ultrasound transducer element size increasing from 0 mm to 0.3 mm on the horizontal axis, and having a legend for colour values on the right side of the figure ranging from amplitude 0 to amplitude 1;</li><li><figref idref="f0048">Figure 76(a)</figref> is a graph of array directivity where the steering angle <i>θ_s</i> equals 0 degrees, showing angle <i>θ</i> decreasing along the vertical axis from 90 degrees to -90 degrees and element size increasing from 0 mm to 0.28 mm along the horizontal axis, and having a legend for colour values on the right side of the figure ranging from amplitude 0 to amplitude 1;</li><li><figref idref="f0048">Figure 76(b)</figref> is a graph of array directivity where the steering angle <i>θ_s</i> equals 30 degrees, showing angle <i>θ</i> decreasing along the vertical axis from 90 degrees to -90 degrees and element size increasing from 0 mm to 0.28 mm along the horizontal axis, and having a legend for colour values on the right side of the figure ranging from amplitude 0 to amplitude 1;</li><li><figref idref="f0048">Figure 77(a)</figref> is the pressure field from a 0.25mm width, 5mm elevation transducer radiating a 7.5 MHz continuous sine wave intersecting on example planar space <i>P₁</i>, showing Z-axis distance increasing from 0 m to 0.015 m along the vertical axis and Y-axis distance increasing from 0 m to 0.005 m along the vertical axis, and having a legend for colour values on the right side of the figure ranging from intensity -0.06 to intensity 0.06;</li><li><figref idref="f0048">Figure 77(b)</figref> is the pressure field from a 0.25mm width, 5mm elevation transducer radiating a 7.5 MHz continuous sine wave intersecting on example planar space <i>P₂</i>, showing Z-axis distance increasing from 0 m to 0.015 m along the vertical axis and Y-axis distance increasing from 0 m to 0.005 m along the vertical axis, and having a legend for colour values on the right side of the figure ranging from intensity -0.06 to intensity 0.06;</li><li><figref idref="f0049">Figure 78</figref> is an example intensity map showing the areas of the intensity map used to generate the Actual Quality Index;</li><li><figref idref="f0050">Figure 79</figref> is a side view of an example probe manipulator having a temperature sensor; and</li><li><figref idref="f0051">Figure 80</figref> is an isometric partial view of the example probe manipulator of <figref idref="f0050">Figure 79</figref>.</li></ul>

0045Example embodiments relate to ultrasound imaging devices and methods for capture and post-processing of ultrasound inspection data. In particular, the described example embodiments relate to devices and methods for inspecting pipe welds using a mechanical cuff that fits around a pipe in the weld region and rotates an ultrasound transceiver array around the circumference of the pipe as the array performs multiple transmit-receive cycles of the pipe volume via the Full Matrix Capture data acquisition technique. All data from the transmit-receive cycles is retained. The data is then post-processed using a two-step algorithm. First, the outer surface of the pipe is modeled by constructing an intensity map of the surface and filtering this map to detect the boundary of the outer surface. Second, the model of the outer surface constructed during the first step is used as a lens in modeling the inner surface of the pipe, using Fermat's principle. The inner surface is modeled the same way as the outer surface: an intensity map is built, then filtered to detect the boundary.

0046The mechanical cuff has a cylindrical outer structure having watertight seals on either end for sealing against a pipe surface. It receives a stream of water via a tube and fills the volume between the structure and the pipe surface with water while in operation in order to facilitate ultrasound scanning.

0047The cuff also has an inner rotating ring having on its inner surface a linear array of ultrasound transceiver crystals with the longitudinal axis of the array aligned along the length of the cylindrical structure, normal to the rotational direction of the inner ring around the circumference of the pipe. The inner ring is automatically rotated around the pipe surface in operation while the outer structure of the cuff remains stationary.

0048Data is acquired by rotating the inner ring around the circumference of the pipe while performing multiple transmit-receive cycles with the ultrasound array for each frame. Each frame uses the Full Matrix Capture technique: a single element is pulsed, with each element in the array measuring the response at that position and storing the resulting time-domain signal (A-scan). This process is then repeated, pulsing each element in turn and recording the response at each element, resulting in a total data corpus of (N x N) A-scans for an array having N elements. The stored time period of each A-scan is determined by monitoring for a signal spike past a set threshold (at time t), then retroactively recording all signal data beginning at a set interval before the spike (at time t-C).

0049In situations where the cuff isn't perfectly normal to the pipe surface at all points around the circumference, it may be preferable to vary the angle of the array to the pipe surface. For this purpose, the inner ring structure incorporates an adjustable reflector or mirror for reflecting ultrasound waves between the transceiver array and the pipe surface at varying angles. The mirror may be adjusted automatically by a local or remote processor or controller module that receives probe data and automatically optimizes the signal quality by adjusting the mirror angle.

0050The post-processing algorithm may feature a number of refinements over the broad outline set out above. Multiple wave modes may be used to improve the reach and resolution of each probe. The outer surface may be modeled as multiple surfaces to further improve resolution of the inner surface where the outer surface is highly irregular. In addition, data from multiple adjacent "slices" of the pipe or other volume may be combined and overlaid to improve the continuity of the surface model, or data from two slices of the same area taken at different times may be overlaid to detect changes in the surfaces over time.

0051While the invention has been described as a pipe inspection tool and technique, the general principles and algorithms are applicable to ultrasound imaging in a number of different contexts and applications.

0052With reference to the drawings, <figref idref="f0001">Figure 1</figref> shows an example embodiment comprising an ultrasound probe manipulator 100. The manipulator 100 comprises a cuff 106 that is fitted around the circumference of a pipe 2 during the scanning process. The center of the cuff 106 is aligned with the longitudinal axis 4 of the pipe 2. The manipulator 100 uses a linear array of ultrasound probe elements, mounted on a carrier (not shown) that traverses the circumference of the cuff 106 by means of a motor 128, to scan the slice of pipe encompassed by the cuff 106.

0053In operation, the cuff 106 is fitted around the pipe 2, with a watertight seal 104 extending from the cuff 106 to the pipe surface. The interior volume defined by the inner surface of the cuff 106, the seal 104, and the outer pipe surface is then filled with water or another fluid suitable for service as an ultrasound scanning medium. In some embodiments, the water is pumped into the interior volume by a hose 110 (shown in <figref idref="f0003">Figure 6</figref>) incorporated into the manipulator 100. The hose 110 is connected to an external water source and/or pump, and feeds into the interior volume of the cuff 106 via a hose intake 132 (shown in <figref idref="f0001">Figure 2</figref>).

0054One or more data connections connect the manipulator 100 to one or more external data processing systems and/or controllers. These external systems may control the operation of the manipulator 100 and/or collect and process the data gathered by the scanning operation of the manipulator 100. <figref idref="f0001">Figure 1</figref> shows a motor connector 130 used to supply power and control data to the motor 128 operative to drive the carrier 102 around the cuff 106. Probe data connectors 108 serve to communicate ultrasound probe control data and data collected by the probe between the probe array and the external data processing systems and/or controllers. In other embodiments, some or all of these functions may take place within the manipulator 100 itself, for example by means of an embedded controller and/or data storage and processing unit. In some embodiments, the motors used by the manipulator 100 may include their own power sources.

0055<figref idref="f0001">Figure 2</figref> shows a similar embodiment to <figref idref="f0001">Figure 1</figref> in isolation instead of fitted to a pipe. The inner surface 112 of the cuff 106 is visible, as is the outer surface 114. The hose intake 132 is also visible in this figure.

0056<figref idref="f0002">Figure 3</figref> is a side view of the manipulator 100 fitted to a curved pipe, showing the longitudinal axis 4 of the pipe portion being scanned. <figref idref="f0002">Figure 4</figref> is an isometric view of the manipulator 100 fitted to a straight pipe, showing the longitudinal axis 4 of the pipe.

0057The manipulator 100 may in some embodiments be fitted or removed from a pipe or other scanning subject by means of a hinged design that allows the cuff 106 to be opened. <figref idref="f0003">Figure 5</figref> shows an example embodiment comprising a hinged manipulator, with a hinged portion 116 allowing the cuff to be opened, and a coupling portion 118 allowing the ends of the cuff to be coupled together into the closed operational position by coupling means such as a latch. <figref idref="f0003">Figure 6</figref> shows the structure of the hinged portion 116, which uses a hinged knuckle 120 to create a double hinge between the two halves of the cuff 106 rather than a simple, single-hinged clamshell design. The hinged knuckle 120 attached to a first half of the cuff 106 at a first connecting point 122, and to the second half of the cuff 106 at a second connecting point 124. Use of a double hinge allows the manipulator 100 to be more easily placed around a pipe circumference due to the greater degrees of freedom afforded. The coupling portion 118 is shown in the example embodiment of <figref idref="f0003">Figure 6</figref> as a latch 128.

0058In operation, the manipulator 100 uses a linear array of ultrasound probe elements, such as resonator crystals, to scan the volume encompassed by the cuff 106. <figref idref="f0004">Figure 7</figref> shows an example linear ultrasound probe array 200 having n elements 202. In operation, the linear array 200 attached to the carrier 102 is aligned parallel to the longitudinal axis 4 of the pipe 2 being scanned. The pipe 2 is scanned by the full array 200 using the Full Matrix Capture technique described below, then the carrier is moved about the circumference of the pipe 2 by a motor included in the manipulator 100, afterwhich the scanning process is repeated for the new circumferential coordinates of the carrier's new position. By performing a number of such scans at regularly-spaced intervals about the circumference of the pipe slice encompassed by the cuff 106, a model of the entire pipe circumference can be built using the scan data.

0059The Full Matrix Capture (FMC) technique used in some embodiments is a known refinement of the phased-array data capture technique widely used for ultrasound scanning. FMC generally requires capturing a larger volume of data than a comparable phased-array scan, but allows more information to be extracted from a single scan. In Full Matrix Capture, a single element 202 of the ultrasound array 200 is pulsed, transmitting ultrasound energy into the medium being scanned. Each 202 element of the array 200 is used as a receiver for this energy, detecting ultrasound vibrations at its coordinates over the time period following this pulse. This detected vibration is recorded and stored for post-processing. Once the data has been recorded for all n elements 202, a second element 202 is pulsed, and the recording process is repeated for all receiving elements 202. This process then repeats again, with each of the n elements 202 being pulsed in turn and data recorded for each receiving element 202, resulting in an n by n matrix of recorded data: each receiving element 202 records scan data from the pulse from each transmitting element 202. This matrix is illustrated by <figref idref="f0004">Figure 8</figref>, showing a matrix of <b>n</b> transmitting elements 206 by <b>n</b> receiving elements 204.

0060In some embodiments, the data from each receiving element 202 is recorded as a series of digital samples taken over time. <figref idref="f0004">Figure 9</figref> shows a three-dimensional matrix of such scan data from a single transmit-receive cycle as described above. The data signal 214 resulting from the pulse of transmitter <b>i</b> 210 captured by receiver <b>j</b> 212 is shown as a series of <b>m</b> samples taken over the time dimension 208, resulting in a total three-dimensional matrix of samples <b>n</b> by <b>n</b> by <b>m</b> in size.

0061In an example embodiment using the manipulator 100 of <figref idref="f0001 f0002 f0003">Figures 1 to 6</figref>, the movement of the carrier 102 and the operation of the ultrasound array 200 is controlled by an external controller connected to the manipulator 100 by the data connections 108. Data recorded by the array 200 is sent to an external data recorder and processor via the data connections 108, where it is stored and processed as further described below. The controller and data processor may also be in communication with each other, and the recorded data may be used by the controller to calibrate or optimize the operation of the carrier 102 and/or array 200 during scanning.

0062A single transmit-receive cycle as described above results in <b>n</b> times <b>n</b> A-scans (i.e., time-domain signals received at a receiving element 202). A single A-scan is generally created by a receiving element by monitoring for vibrations above a set threshold, then recording sensed vibrations and for a set period of time after this threshold is crossed. In some embodiments, a buffer is used to store the sensed data prior to recording, and the recording period is set to include buffered data for a predetermined period before the threshold is crossed, thereby capturing a period beginning shortly before the threshold is crossed and lasting for a set period of time. <figref idref="f0042">Figure 68</figref> shows an illustration of this process. The sensed data 6800 is sampled until a peak 6804 is detected that exceeds a predetermined threshold 6806, thereby signalling the beginning of the period of interest of the signal. However, there may exist data of interest prior to the peak 6804, such as the initial oscillations of the signal 6800 beginning at point 6802. In order to capture these initial data points that precede the peak 6804, buffered data points are retained extending back for a predetermined period before the peak 6804, such as back to an earlier point 6808 which is early enough to capture any initial signal perturbations of interest.

0063Processing of the captured data may be done concurrently with the scan or afterward. Techniques for processing the captured data are described below in accordance with example embodiments. These techniques may involve application of the Shifting Aperture Focusing Method (SFM), the Interior Focus Method (IFM), and boundary detection and recognition to determine the structure of a scanned object, such as the inner and outer surface contours of a pipe wall. These techniques may allow the detection of subtle variations in pipe thickness, defects in pipe walls, and other structural details of arbitrary inner and outer surfaces of a pipe. Some of the mathematical principles applied by various embodiments are hereby described to more fully explain their operation.

0064<figref idref="f0009">Figure 17</figref> is a flowchart showing the operations involved in modeling the outer diameter (OD) and inner diameter (ID) surfaces of a scanned portion of a pipe wall or other object according to an example embodiment. The Full Matrix Capture (FMC) data from a transmit-receive cycle of the probe array 200 is collected at step 1702. At step 1704, the raw FMC data is pre-processed. The OD boundary is modeled in steps 1706 through 1710, then this OD boundary definition 1712 is used to determine the ID boundary 1724 in conjunction with the raw FMC data in steps 1716 through 1722. These steps are described more fully below. In some embodiments, the raw FMC data used includes data from multiple transmit-receive cycles, which is used to improve the modeling of adjacent radial positions of the pipe wall.

0065<figref idref="f0005">Figure 10</figref> shows the vector notation used in describing the processing of A-scan data. The figure depicts a plane defined by the linear ultrasound array 200 and the longitudinal axis 4 of the pipe 2 being scanned. The area of the plane being scanned in this example is shown by image area 14.

0066<b>k</b> is a vector valued function such that <b>k</b><i>(u)</i>=<<i>x,y</i>>=<<i>f(u),g(u)</i>> is a vector in 2-space, <i>R2,</i> where <i>u</i> ∈ <i>[u1,u2].</i> The point <i>k(u)</i> is defined as <i>k(u)</i>=<i>(x,y)</i>=<i>(f(u),g(u)).</i> The vector <b>k</b><i>(u)</i>, represented in boldface, is referred to distinctly from the point <i>k(u).</i> Throughout this description, for any vector <b>v</b>=<x,y>, the point <i>(x,y)</i> is denoted by <i>v,</i> in italics and not boldface.

0067<i>K</i> 10 is the curve defined by the set of all points <i>k(u), u</i> ∈ <i>[u₁,u₂]. K</i> is piecewise-smooth and non self-intersecting. The endpoints of <i>K</i> are <i>k(u₁) and k(u₂)</i> which are denoted <i>k₁</i> 22 and <i>k₂</i> 24 respectively.

0068<i>K</i> separates a first medium 6 from a second medium 8 where ultrasound probe element <i>i</i> 210 lies in the first medium 6 and the endpoint of vector <b>r</b> 16 lies in the second medium 8. For example, when scanning a pipe circumference, the first medium 6 would be composed of water pumped into the interior volume between the cuff 106 and the outer surface of the pipe 2, while the second medium 8 would be the metal of the pipe wall itself. The speed of sound in the first medium 6 and the second medium 8 are denoted as <i>c1</i> and <i>c2</i> respectively.

0069With reference again to <figref idref="f0005">Figure 10</figref>, the figure shows vector <b>k_o,</b> defined by <b>ko</b>=<b>k</b>(uo)=<f(uo),g(uo)>, where the point ko ∈ <i>K</i> (point k_o being denoted by numeral 18, surface <i>K</i> by numeral 10). Also, ultrasound element <i>i</i> 210 is shown with position vector <b>e(i)</b> from coordinate origin 20. The travel time from ultrasound element <i>i</i> 210 to r 16 through <i>k</i>o 18, denoted by <maths id="math0015"><math display="inline"><msup><msub><mi>t</mi><mi mathvariant="italic">ir</mi></msub><mi mathvariant="italic">ko</mi></msup></math><img file="EP3553561B1_D0015.tif" /></maths>, is given by the following equation: <maths id="math0016" num="(Equation 1)"><math display="block"><msubsup><mi>t</mi><mi mathvariant="italic">ir</mi><msub><mi>k</mi><mn>0</mn></msub></msubsup><mo>=</mo><mfrac><mfenced open="|" close="|"><msub><mi mathvariant="bold">k</mi><mn>0</mn></msub><mo>−</mo><msub><mi mathvariant="bold">e</mi><mfenced><mi mathvariant="normal">i</mi></mfenced></msub></mfenced><msub><mi>c</mi><mn>1</mn></msub></mfrac><mo>+</mo><mfrac><mfenced open="|" close="|"><mi mathvariant="bold">r</mi><mo>−</mo><msub><mi mathvariant="bold">k</mi><mn>0</mn></msub></mfenced><msub><mi>c</mi><mn>2</mn></msub></mfrac></math><img file="EP3553561B1_D0016.tif" /></maths>

0070The travel time from <i>r</i> 16 to <i>i</i> 210 through <i>k_o</i> 18 (denoted <maths id="math0017"><math display="inline"><msup><msub><mi>t</mi><mi mathvariant="italic">ri</mi></msub><mi mathvariant="italic">ko</mi></msup></math><img file="EP3553561B1_D0017.tif" /></maths>) is equal to <maths id="math0018"><math display="inline"><msup><msub><mi>t</mi><mi mathvariant="italic">ir</mi></msub><mi mathvariant="italic">ko</mi></msup></math><img file="EP3553561B1_D0018.tif" /></maths>. The times from <i>i</i> 210 to <i>r</i> 16 through all <i>k</i> ∈ <i>K, k</i> varying parametrically from <i>k₁</i> to <i>k₂</i>, are given by the curve, <i>T_ir(K)</i> 26, illustrated in <figref idref="f0006">Figure 11</figref>.

0071In ultrasonic applications, absolute distance measurements may be directly calculated from the travel times of acoustic pulses, and hence may be sensitive to the speed of ultrasonic sound in the materials under study. In some embodiments described below pertaining to pipe weld inspection, thickness is defined as the shortest distance from a point on the outer surface to the inner surface. <figref idref="f0046">Figure 72</figref> illustrates this definition of wall thickness. L1 and L2 are the thicknesses at two different locations on the outer surface of a feeder weld. The inner and outer surface profiles are essential pieces of information used to determine the wall thickness with respect to any location on the outer or inner surface.

0072The ultrasonic distance measurement is determined by the sound velocity of the material and the time taken for a sound wave to travel between start and end points. A distance measurement L can be written as <i>L</i> = <i><o ostyle="single">V</o>τ</i> where <i>τ</i> is the time taken by the wave to travel from start and end points. <i><o ostyle="single">V</o></i> is the average sound velocity along the path. Since the outer surface profile may be measured using immersion techniques and water may be used as a couplant in some embodiments, the sound velocity of water would need to be considered, as it would affect the distant measurement accuracy.

0073The sound velocity of water as a function of its temperature is <i>V_fw(T)</i> = 1405.03+ 4.624<i>T</i> - 0.0383<i>T</i>² where <i>V_fw</i> denotes the fresh water sound velocity in meters per second unit and <i>T</i> denotes the temperature in degrees Celsius.

0074Element directivity is a potentially important factor in the design of the probe array. Generally speaking, element directivity can be thought of as the variance of the amplitude pressure field across different points on the inspection volume. Both the Cartesian (x-y-z) and spherical (<i>ϕ</i> - <i>θ</i> - <i>r</i>) coordinates are standard when discussing directivity (see <figref idref="f0046">Figure 73</figref>). Convenience dictates that the Cartesian coordinates are of use when discussing near field element directivity while spherical coordinates are of use when discussing far field directivity. In the far field, for a rectangular element with <i>L</i> much longer than <i>a</i>, and for a given distance, <i>r,</i> the directivity in the x-z plane is well approximated via a function of only <i>θ</i>, excitation frequency <i>f</i>, and <i>a</i>.

0075Element directivity is ultimately derived from the wave equation, given by the following: <maths id="math0019"><math display="block"><mfrac><mrow><msup><mo>∂</mo><mn>2</mn></msup><mi>p</mi></mrow><mrow><mo>∂</mo><msup><mi>x</mi><mn>2</mn></msup></mrow></mfrac><mo>+</mo><mfrac><mrow><msup><mo>∂</mo><mn>2</mn></msup><mi>p</mi></mrow><mrow><mo>∂</mo><msup><mi>y</mi><mn>2</mn></msup></mrow></mfrac><mo>+</mo><mfrac><mrow><msup><mo>∂</mo><mn>2</mn></msup><mi>p</mi></mrow><mrow><mo>∂</mo><msup><mi>z</mi><mn>2</mn></msup></mrow></mfrac><mo>−</mo><mfrac><mn>1</mn><msup><mi>c</mi><mn>2</mn></msup></mfrac><mfrac><mrow><msup><mo>∂</mo><mn>2</mn></msup><mi>p</mi></mrow><mrow><mo>∂</mo><msup><mi>t</mi><mn>2</mn></msup></mrow></mfrac><mo>=</mo><mn>0</mn></math><img file="EP3553561B1_D0019.tif" /></maths> where <i>p</i> is pressure and <i>t</i> is time.

0076The pressure field at a given point in an inspection volume can be derived numerically. A transducer can be treated as a piston radiating sound waves in water where the transducer generates an infinite number of plane waves, all traveling in the positive z-direction but with different x and y component directions. As such, the pressure field at a point <i>q</i> = <i>(x; y; z)</i> is represented in the form of a 2D integral given by: <maths id="math0020"><math display="block"><mi>p</mi><mfenced><mi>q</mi><mi>ω</mi></mfenced><mo>=</mo><msup><mfenced><mfrac><mn>1</mn><mrow><mn>2</mn><mi>π</mi></mrow></mfrac></mfenced><mn>2</mn></msup><mstyle displaystyle="true"><mrow><msubsup><mo>∫</mo><mrow><mo>−</mo><mi>∞</mi></mrow><mi>∞</mi></msubsup><mstyle displaystyle="true"><mrow><msubsup><mo>∫</mo><mrow><mo>−</mo><mi>∞</mi></mrow><mi>∞</mi></msubsup><mrow><mi>P</mi><mfenced><msub><mi>k</mi><mi>x</mi></msub><msub><mi>k</mi><mi>y</mi></msub></mfenced><msup><mi>e</mi><mrow><mi>i</mi><mfenced><msub><mi>k</mi><mi>x</mi></msub><mi>x</mi><mo>+</mo><msub><mi>k</mi><mi>y</mi></msub><mi>y</mi><mo>+</mo><msub><mi>k</mi><mi>z</mi></msub><mi>z</mi></mfenced></mrow></msup><msub><mi mathvariant="italic">dk</mi><mi>x</mi></msub><msub><mi mathvariant="italic">dk</mi><mi>y</mi></msub></mrow></mrow></mstyle></mrow></mstyle></math><img file="EP3553561B1_D0020.tif" /></maths> where <i><o ostyle="single">k</o></i> = (<i>k_x</i>, <i>k_y</i>, <i>k_z</i>) is the wave vector with magnitude <maths id="math0021"><math display="inline"><mi>k</mi><mo>=</mo><mfrac><mrow><mn>2</mn><mi>π</mi></mrow><mi>λ</mi></mfrac></math><img file="EP3553561B1_D0021.tif" /></maths>.

0077After some derivation, a solution to the above equation for a rectangular element with length <i>l_x</i> in the <i>x</i> direction and length <i>l_y</i> in the <i>y</i> direction can be derived where: <maths id="math0022"><math display="block"><mi>p</mi><mfenced><mi>q</mi><mi>ω</mi></mfenced><mo>=</mo><msup><mfenced><mfrac><mn>1</mn><mrow><mn>2</mn><mi>π</mi></mrow></mfrac></mfenced><mn>2</mn></msup><mstyle displaystyle="true"><mrow><msub><mo>∬</mo><mrow><msubsup><mi>k</mi><mi>x</mi><mn>2</mn></msubsup><mo>+</mo><msubsup><mi>k</mi><mi>y</mi><mn>2</mn></msubsup><mo>≤</mo><msup><mi>k</mi><mn>2</mn></msup></mrow></msub><mfrac><mrow><mi mathvariant="italic">iωρV</mi><mfenced><msub><mi>k</mi><mi>x</mi></msub><msub><mi>k</mi><mi>y</mi></msub></mfenced></mrow><msub><mi mathvariant="italic">ik</mi><mn>2</mn></msub></mfrac></mrow></mstyle><msup><mi>e</mi><mrow><mi>i</mi><mfenced><msub><mi>k</mi><mi>x</mi></msub><mi>x</mi><mo>+</mo><msub><mi>k</mi><mi>y</mi></msub><mi>y</mi><mo>+</mo><msub><mi>k</mi><mi>z</mi></msub><mi>z</mi></mfenced></mrow></msup><msub><mi mathvariant="italic">dk</mi><mi>x</mi></msub><msub><mi mathvariant="italic">dk</mi><mi>y</mi></msub></math><img file="EP3553561B1_D0022.tif" /></maths> where <i>V(k_x,k_y)</i> is the 2D spatial Fourier transform of the velocity of the field of the transducer: <maths id="math0023"><math display="block"><mi>V</mi><mfenced><msub><mi>k</mi><mi>x</mi></msub><msub><mi>k</mi><mi>y</mi></msub></mfenced><mo>=</mo><mstyle displaystyle="true"><mrow><msubsup><mo>∫</mo><mrow><mo>−</mo><mi>∞</mi></mrow><mi>∞</mi></msubsup><mstyle displaystyle="true"><mrow><msubsup><mo>∫</mo><mrow><mo>−</mo><mi>∞</mi></mrow><mi>∞</mi></msubsup><mrow><msub><mi>v</mi><mi>x</mi></msub><mfenced><mi>x</mi><mo>,</mo><mi>y</mi><mo>,</mo><mi>z</mi><mo>=</mo><mn>0</mn><mo>,</mo><mi>ω</mi></mfenced><msup><mi>e</mi><mrow><mo>−</mo><mi>i</mi><mfenced><msub><mi>k</mi><mi>x</mi></msub><mi>x</mi><mo>+</mo><msub><mi>k</mi><mi>y</mi></msub><mi>y</mi></mfenced></mrow></msup><mi mathvariant="italic">dxdy</mi></mrow></mrow></mstyle></mrow></mstyle></math><img file="EP3553561B1_D0023.tif" /></maths>

0078In the far field of an element (when the ultrasound transducer element dimension is very small with respect to inspection distance), for a given radius <i>r</i> in the x-z plane, the directivity can be well approximated by a function varying only in <i>θ</i>, element width <i>a,</i> and frequency <i>f,</i> when the transducer length <i>L</i> is much bigger than its width, <i>a.</i> Far field directivity may be relevant in some embodiments involving weld inspection because the weld will often be in the far field (in the axial direction, x-z plane) of the ultrasound array transducer elements. The formula for far field directivity of the element in the x-z plane is given below: <maths id="math0024"><math display="block"><msub><mi>D</mi><mi>F</mi></msub><mfenced><mi>f</mi><mi>θ</mi></mfenced><mo>=</mo><mi mathvariant="italic">sinc</mi><mfenced><mfrac><mrow><mi mathvariant="italic">πaf</mi><mspace width="1ex" /><mi>sin</mi><mspace width="1ex" /><mi>θ</mi></mrow><mi>c</mi></mfrac></mfenced></math><img file="EP3553561B1_D0024.tif" /></maths> It is worth noting that both <i>a</i> and <i>f</i> increase directivity at <i>a</i> given angle <i>θ</i>. <figref idref="f0047">Figure 74</figref> illustrates directivity in an example embodiment having <i>θ</i> = -90 degrees to 90 degrees, <i>a</i> = 0.23 mm, <i>f</i> = 7.5 MHz, and <i>c</i> = 1480 m/s.

0079Element width a dictates element directivity as is illustrated in the equation above. Smaller element widths will radiate sound omnidirectionally (in all directions). Larger element widths will focus sound in the direction of their surface normals. This is illustrated by <figref idref="f0047">Figure 75</figref>, which shows element directivity as a function of element size a.

0080Analogously, with respect to focusing arrays, arrays with larger element widths will generally focus better than arrays with smaller widths in the direction of the array normal. The direction of the array normal corresponds to a steering angle of zero. On the other hand, arrays with smaller element widths will generally focus better (than arrays with larger element widths) in directions away from the array surface normal.

0081Both claims can be confirmed by examination of <figref idref="f0048">Figure 76(a) and Figure 76(b)</figref>. Both figures simulate the same array with: <ul id="ul0002" list-style="bullet" compact="compact"><li>Number of Elements = 10</li><li>Frequency = 7.5MHz</li><li>Pitch (spacing between centre of elements) = 0.28mm</li></ul>

0082While the element width a is varied between 0 and the pitch (0.28mm), <figref idref="f0048">Figure 76(a)</figref> simulates the array directivity where the steering angle <i>θ_s</i> equals 0 degrees while <figref idref="f0048">Figure 76(b)</figref> simulates the array directivity where the steering angle <i>θ_s</i> equals 30 degrees. The intensity of the main lobe at different element widths is visible by examining the line <i>θ_s</i> = <i>θ</i> for both figures, while the other horizontal lines represent the undesirable effects of grating. For <i>θ_s</i> = 0 (<figref idref="f0048">Figure 76(a)</figref>), larger element widths attenuate the effects of grating while preserving the desirable main lobe intensity. Conversely, for <i>θ_s</i> = 30 (<figref idref="f0048">Figure 76(b)</figref>), smaller element widths preserve the main lobe intensity while grating lobe intensity is unvaried by element size.

0083Element elevation is represented as <i>L</i> in <figref idref="f0046">Figure 73</figref>. In the far field, a large <i>L</i> (with respect to <i>a</i>) will focus energy in the direction normal to the surface (the <i>z</i> direction) while smaller values for <i>L</i> will radiate energy with larger components in <i>y</i>. In the near field, however, which is generally of interest when inspections are performed at distances comparable to <i>L</i>, energy radiated from the transducer will be projected uniformly from the surface of the transducer in the <i>z</i> direction. This is consistent with expectations, as at these inspection distances, inspection points will be comparable to the elevational focal length of the transducer.

0084Recall the coordinate system of <figref idref="f0046">Figure 73</figref>. <figref idref="f0048">Figure 77(a) and Figure 77(b)</figref> reveal the pressure field from a 0.25mm width, 5mm elevation transducer radiating a 7.5 MHz continuous sine wave intersecting on two planar spaces <i>P₁</i> = <i>{x</i> = 0, <i>0</i> ≤ <i>y</i> ≤ <i>5, 0</i> ≤ <i>z</i> ≤ <i>20}</i> and <i>P₂</i> = <i>{x</i> = <i>1</i>, <i>0</i> ≤ <i>y</i> ≤ <i>5, 0</i> ≤ <i>z</i> ≤ <i>20}</i> (all dimensions are in millimetres). The first space <i>P₁</i> is a plane spanned by vectors in the <i>y</i> and z direction and bounded at <i>y</i> = 0,5, and <i>z =</i> 0,20. Space <i>P₂</i> is spanned and limited in the same manner as <i>P₁</i>, but is spaced 1mm away from <i>P₁</i>. Both spaces are shown in <figref idref="f0047">Figure 75</figref>. <figref idref="f0048">Figure 77(a) and Figure 77(b)</figref> illustrate that the intensity of the pressure field radiated by the transducer is concentrated in the y &lt; 2.5 mm area for both cases simulated. The symmetry of the radiated field is exploited in the simulations given. If the radiated pressure field was simulated at planar spaces <maths id="math0025"><math display="inline"><msubsup><mi>P</mi><mn>1</mn><mo>′</mo></msubsup><mo>=</mo><mfenced open="{" close="}"><mi>x</mi><mo>=</mo><mn>0</mn><mo>,</mo><mo>−</mo><mn>5</mn><mo>≤</mo><mi>y</mi><mo>≤</mo><mn>0</mn><mo>,</mo><mo>−</mo><mn>20</mn><mo>≤</mo><mi>z</mi><mo>≤</mo><mn>0</mn></mfenced></math><img file="EP3553561B1_D0025.tif" /></maths> and <maths id="math0026"><math display="inline"><msubsup><mi>P</mi><mn>2</mn><mo>′</mo></msubsup><mo>=</mo><mfenced open="{" close="}"><mi>x</mi><mo>=</mo><mn>1</mn><mo>,</mo><mo>−</mo><mn>5</mn><mo>≤</mo><mi>y</mi><mo>≤</mo><mn>0</mn><mo>,</mo><mo>−</mo><mn>20</mn><mo>≤</mo><mi>z</mi><mo>≤</mo><mn>0</mn></mfenced></math><img file="EP3553561B1_D0026.tif" /></maths> which are reflections of planar spaces <i>P₁</i> and <i>P₂</i> about the z-x plane, the intensity field simulated on these spaces would be a reflection of those on <i>P₁</i> and <i>P₂</i> about the z-x plane.

0085In analog-to-digital conversion, the magnitude of each data sample is converted into an approximated value with finite precision. Quantization is a non-linear process. The smallest quantization level is the resolution. It is determined by the full-scale input amplitude of an analog/digital (A/D) converter and the total number of quantization levels which are usually evenly spaced. The resolution is often expressed by the number of quantization level. A 10-bit A/D converter has 1024 quantization levels. A 12-bit A/D converter has 4096 quantization levels. The resolution of a 12-bit AID converter is four times smaller than that of a 10-bit A/D converter. Since the quantization error is expected to be smaller than the smallest quantization level, a higher resolution A/D converter is in general preferred.

0086If the quantization process rounds the input data value to the nearest quantization level and the errors obey the even-statistical distribution over the quantization intervals, the mean value of the quantization error is clearly zero. The variance of the quantization error is given by <maths id="math0027"><math display="inline"><msubsup><mi>σ</mi><mi>e</mi><mn>2</mn></msubsup><mo>=</mo><mfrac><msup><mi mathvariant="normal">Δ</mi><mn>2</mn></msup><mn>12</mn></mfrac></math><img file="EP3553561B1_D0027.tif" /></maths> where Δ is the quantization interval.

0087The root-mean-squared (rms) value of the quantization error is the standard deviation <maths id="math0028"><math display="inline"><msub><mi>σ</mi><mi>e</mi></msub><mo>=</mo><mfrac><mi mathvariant="normal">Δ</mi><msqrt><mn>12</mn></msqrt></mfrac><mo>≈</mo><mn>0.29</mn><mi mathvariant="normal">Δ</mi></math><img file="EP3553561B1_D0028.tif" /></maths>.

0088If we define a SNR to be the ratio of signal variance to the noise variance, the SNR of a B+1 bit A/D converter can be expressed in some embodiments as: <maths id="math0029"><math display="block"><mi mathvariant="italic">SNR</mi><mo>=</mo><mn>6.02</mn><mi>B</mi><mo>+</mo><mn>10.8</mn><mo>−</mo><mn>20</mn><msub><mi>log</mi><mn>10</mn></msub><mfenced><mfrac><msub><mi>X</mi><mi>m</mi></msub><msub><mi>σ</mi><mi>x</mi></msub></mfrac></mfenced></math><img file="EP3553561B1_D0029.tif" /></maths> where <i>X_m</i> is the full scale level of the A/D converter, and <i>σ_x</i> is the standard deviation of the signal.

0089The SNR limits for 8-bit, 10-bit and 12-bit in this example embodiment are 50 dB, 62 dB and 74 dB respectively. It is worth noting that the optimum SNR can generally only be achieved when the input signal is carefully adjusted to the full-scale amplitude of the A/D converter.

0090In some embodiments, preprocessing consists of several operations which condition raw data for analysis via the SFM and IFM subroutines. In an example embodiment shown in <figref idref="f0010">Figure 18</figref>, these operations are as follows: upsampling the full matrix capture raw data to a sampling frequency of 100MHz 1804, subtracting the DC offset from the full matrix capture raw data 1806, filtering the full matrix capture data set to remove unwanted noise 1808 using digital software filter coefficients 1810, and calculating 1812 the analytic time-domain full matrix capture data-set 1814 from the acquired RF full matrix capture data set 1802.

0091In this example embodiment, full matrix capture (FMC) raw data 1802 is collected by the acquisition system at a sampling frequency of 50MHz. Analyzing the raw data collected at this frequency may deliver results with insufficient accuracy. For this reason, the raw data is upsampled to 100MHz at step 1804. The acquisition system in this embodiment is sensitive to frequencies less than 25MHz, and data is collected at twice this rate, so due to the Nyquist Sampling Theorem, the raw data can be perfectly reconstructed at any sampling rate above 50MHz. The FMC raw data is upsampled from 50MHz to 100MHz.

0092The acquisition system stores raw data via a 12-bit quantization scheme where only positive values are stored. In this embodiment, FMC waveforms are centered about 2¹²/2=2048. The analysis algorithms require that waveforms be centered about zero. In some embodiments, the theoretical DC offset value of 2048 may not be exactly accurate: there may be a DC offset inherent in the hardware controlling the ultrasound probe array. Experimentation may show the right DC offset value to use to get a zero value as closely as possible to reality; in some embodiments using specific hardware, the value is 2058. Therefore a DC offset of 2058 is subtracted from each FMC waveform. The exact value of the DC offset may be a user configurable parameter which can be changed to account for other deviations from the theoretical value. The DC offset is subtracted at step 1806.

0093Unwanted frequency content in the full matrix capture data will sometimes be present due to various noise contributions. These frequencies can be attenuated at step 1808 through the utilization of a digital software filter. Software filtering coefficients may be specified in the filtering process, such as parameters derived from the Filterbuilder program in the Matlab^™ software application, at step 1810.

0094Assigning intensities to points in the inspection medium requires the full matrix data-set of analytic time-domain signals. The full matrix data-set output of the acquisition system 1802 contains the RF data-set (the real part of the analytic time-domain signals). To compute the analytic time-domain signals 1814 from the RF data-set at step 1812, the Hilbert transform of the RF data-set is calculated, multiplied by the imaginary number, i, and added to the RF data-set.

0095Detailed descriptions of functions used in an example embodiment of the system are set out in Table A1 at the end of the Description.

0096Once pre-processing 1704 is done, the Shifting Aperture Focus Method (SFM) is applied to the pre-processed data at step 1706. The Shifting Aperture Focus Method is an algorithm whose purpose is to output an intensity map given full matrix capture raw data. The operation of the SFM to determine the OD intensity map is shown in the flowchart of <figref idref="f0012">Figure 20</figref>. The pre-processed FMC data 2004 may or may not be normalized first, depending on the embodiment or on user-defined parameters for processing. The decision to normalize is made at step 2010. If the data 2004 is to be normalized, normalization occurs at step 2012 as further described below, based on OD normalization parameters 2006 either pre-defined or set by a user. The normalized or non-normalized data is then used to calculate inspection coordinate travel times at step 2014 as further described below, a process which may take into account intensity coordinates calculated at step 2016 based on OD imaging parameters 2002. At step 2018, the intensities at the current inspection coordinates are calculated. At step 2020, the inspection coordinates and their respective intensities are stored. At step 2022, the algorithm may focus further around high intensity coordinates; if it does, the intensity coordinates are calculated at step 2016, creating an iterative loop for steps 2014 to 2022. When the algorithm has iterated through this process one or more times, it stops re-focusing and outputs an OD intensity map 2008. These various steps are described below.

0097The Shifting Aperture Focusing Method is a variation on the Total Focusing Method (TFM). The Total Focusing Method is a known technique for imaging in a single medium where sound travels at speed <i>c.</i>

0098In the Shifting Aperture Focusing Method, first, an intensity function is computed for each point <i>r</i> in the scanning area by summing a function <i>g_(i)j(t)</i> for each transmitting element <i>i</i> and each receiving element <i>j</i> within a fixed-width aperture <i>a</i> spanning a set number of adjacent elements in the array. (in some embodiments, the width of the defined aperture may be user-configurable.) Ultrasound elements <i>i</i> 210 and <i>j</i> 212 belong to aperture <i>a. g_(i)j(t)</i> is the amplitude of the data-set of analytic time-domain signals from transmitter <i>i</i> 210 to receiver <i>j</i> 212 at time <i>t</i> (note that <i>g_(i)j(t)</i> is defined for every i and <i>j,</i> since the full matrix of ultrasonic transmit-receive array data is acquired). <i>i</i> 210 is enclosed in parentheses to represent it as the transmitter, while <i>j</i> 212 is left without parenthesis to represent it as the receiver. The intensity at <i>r</i> is defined as: <maths id="math0030" num="(Equation 2)"><math display="block"><mi>l</mi><mfenced><mi>r</mi><mi>a</mi></mfenced><mo>=</mo><mfenced open="|" close="|"><mstyle displaystyle="true"><msub><mo>∑</mo><mrow><mi>i</mi><mo>,</mo><mi>j</mi><mo>∈</mo><mi>a</mi></mrow></msub><mrow><msub><mi>g</mi><mrow><mfenced><mi>i</mi></mfenced><mi>j</mi></mrow></msub><mfenced><mi>t</mi><mo>=</mo><mfrac><mrow><mfenced open="|" close="|"><msub><mi>e</mi><mfenced><mi>i</mi></mfenced></msub><mo>−</mo><mi>r</mi></mfenced><mo>+</mo><mfenced open="|" close="|"><msub><mi>e</mi><mi>j</mi></msub><mo>−</mo><mi>r</mi></mfenced></mrow><mi>c</mi></mfrac></mfenced></mrow></mstyle></mfenced></math><img file="EP3553561B1_D0030.tif" /></maths> where <b>r</b> is the vector defining point <i>r</i> relative to a coordinate origin, e_(i) is a vector defining the position of transmitter <i>i</i> relative to the coordinate origin, e_j is a vector defining the position of receiver <i>j</i> relative to the coordinate origin, and <i>c</i> is the speed of light.

0099Total focusing is achieved by calculating the above for every point in the imaging area 14 (<b>r</b> is varied). By varying <b>r</b> over a set of points sufficiently close to each other such that <i>I</i>(<b>r</b>,<i>a</i>) does not vary significantly between adjacent values of <b>r,</b> an image of the inspection medium can be formed with image pixels of intensity <i>I</i>(<b>r</b>,<i>a</i>) at positions <i>r</i>. This image is called an intensity map.

0100The next step in the Shifting Aperture Focusing Method is to shift the aperture <i>a</i> along the array and perform the same calculation again for the new aperture. After intensities have been computed for each aperture, the highest such intensity value is used to represent the intensity of point <i>r</i> in the intensity map. Thus, SFM addresses the problem of how to relate intensities from different apertures in assigning an intensity value to <i>r.</i> In many instances, a surface contains reflectors that vary with respect to which apertures they reflect sound best to. For example, one reflector may return sound well to apertures <i>a</i>₁ through <i>a</i>₆, while another reflector may return sound well to only aperture <i>a</i>₃. In imaging a surface, however, reflectors should be imaged with equal intensity, irrespective of how many apertures individual reflectors reflect sound well to. This will lead to increased intensity levelness in imaging the surface. To this end, <i>I(r)</i> is defined as the maximum intensity at <i>r</i> of the set of computed intensities at <i>r</i> with respect to apertures <i>a</i> ∈ A. <i>I(r)</i> is defined as: <maths id="math0031" num="(Equation 3)"><math display="block"><mi>I</mi><mfenced><mi>r</mi></mfenced><mo>=</mo><munder><mi>max</mi><mrow><mi>a</mi><mo>∈</mo><mi>A</mi></mrow></munder><mspace width="1ex" /><mfenced open="{" close="}"><mi>I</mi><mo>′</mo><mfenced><mi>r</mi><mi>a</mi></mfenced></mfenced></math><img file="EP3553561B1_D0031.tif" /></maths>

0101Implementation of the SFM routine can become very computationally intensive if there are many coordinates for which corresponding intensities are evaluated. Limiting the number of coordinates under consideration, while focusing at the appropriate density to meet inspection specifications, may require careful implementation of a focusing strategy. The strategy employed in some embodiments is to first calculate intensities of coordinates on a course grid lying in the area of inspection.

0102Detailed descriptions of OD Imaging Parameters used in an example embodiment of the system are set out in Table A2 at the end of the Description. Detailed descriptions of SFM Functions used in an example embodiment of the system are set out in Table A3 at the end of the Description.

0103The Total Focusing Method and Shifting Aperture Focusing Method (as well as a further variation, the Interior Focusing Method, described below) may use Beam Normalization. Beam-steering at angles away from the direction normal to an array may be optimized when individual elements are omnidirectional. A correction factor may be introduced in imaging to emulate beam-spread omnidirectionality. A potential problem with this method is that it presumes signals found in <i>g(i)_j(t)</i> are located in the direction of <b>r,</b> when assigning a value to <i>I</i>(<b>r</b>,a). While this may not be a problem of great concern when attempting to image small objects, when imaging surfaces this may lead to amplified grating. Amplified grating causes a reduction in signal to noise ratio in areas of the image where the true surface and grating overlap. The method presented here for normalizing beam-spread is to normalize wave-packets found in the real and imaginary parts of <i>g(i)_j(t)</i> such that the envelope of <i>g(i)_j(t),</i> |<i>g(i)_j(t)</i>|, has peak(s) equal to an arbitrary constant. For simplicity, in this example this constant is equal to 1. Thus, if <i>g(i)_j(t)</i> contains wave packets <i>W</i>=<i>{w₁,w₂,...,w_n}</i> and |<i>g(i)_j(t)</i> |has peaks <i>P</i>=<i>{p₁,p₂,...,p_n},</i> to normalize the peaks of |<i>g(i)_j(t)</i>| to 1, the wave packets <i>W</i> are scaled by <i>{p₁^-1,p₂^-1,...,p_n^-1}.</i> Let <i>g'(i)_j(t)</i> denote <i>g(i)_j(t)</i> with normalized wave packets. <figref idref="f0008">Figure 15 and Figure 16</figref> illustrate this concept of analytic time domain signal normalization. The A-scan data 1500 in <figref idref="f0008">Figure 15</figref> has a real part 1502 that exhibits a first peak 1506 and a second peak 1508. The envelope 1504 of the signal is higher for the first peak 1506 than for the second peak 1508. After applying wave packet normalization as described above, the normalized A-scan 1600 is shown in <figref idref="f0008">Figure 16</figref>. The real part 1602 of the A-scan has had the envelope 1604 of its first peak 1606 and its second peak 1608 normalized to the same constant value.

0104At step 2022, depending on preset or user defined parameters, the SFM subroutine either ends, outputting the intensity map 2008 to the boundary detection subroutine, or proceeds to define new coordinates for which to calculate corresponding intensities at step 2016. If the latter course is taken, newly defined coordinates will be positioned around coordinates with high intensities already assigned to them. The cutoff intensity for coordinates of which newly defined coordinates focus around may be defined by a predefined vector. Once the new coordinates have been defined, potentially depending on user defined parameters, the SFM subroutine exits, or proceeds to further focus around coordinates of high intensity. The process of identifying high intensity coordinates and then refocusing around them can be executed an arbitrary number of times, which may be specified by a user in some embodiments.

0105One cycle of the refocusing process (steps 2014 to 2022) is illustrated in <figref idref="f0011">Figure 19</figref>. Coordinates spaced coarsely apart are represented as either white circles 1902 or black circles 1904. The white circles 1902 represent those coordinates whose respective intensities are below the cutoff intensity for which new coordinates are defined. Conversely, the black circles 1904 represent those coordinates whose respective intensities exceed the cutoff intensity for which new coordinates are defined. The gray circles 1906 represent newly defined coordinates around high intensity coordinates. If the gray circles 1906 are defined on iteration i+1 of the coordinate definition and intensity assigning process, dx(i)/dx(i+1) = dz(i)/dz(i+1) = 4 for the example illustrated by <figref idref="f0011">Figure 19</figref> with x-axis 1908 and z-axis 1910.

0106An example intensity map of the OD (outer surface) of a pipe 2 is shown in <figref idref="f0007">Figure 12</figref>. The OD intensity map 50 is mapped within the same plane as <figref idref="f0006">Figure 11</figref>, with the depth 28 of the scan as the vertical axis and the axial position 30 along the longitudinal axis 4 as the horizontal axis. The high-intensity OD regions 32 in the OD intensity map 50 indicate the shape of the outer pipe surface in the axial direction at the radial position of the probe array 200 during the current transmit-receive cycle.

0107The OD intensity map 50 can be further processed, either on its own or in conjunction with neighbouring intensity maps from separate transmit-receive cycles, to build a model of the outer pipe surface. In the example embodiment shown in <figref idref="f0009">Figure 17</figref>, boundary recognition is performed at step 1708 followed by boundary definition at step 1710.

0108Boundary recognition is executed as follows. Given an intensity map (stored in sparse coordinates in some embodiments) of the OD coordinates or an intensity map of the ID coordinates along with corresponding boundary recognition parameters, the boundary recognition algorithm will output the surface boundary.

0109Boundary recognition and definition are intended to define the true boundary of the OD/ID (not edges of aberrations in the image), and the extraction of the boundary in the form of a set of coordinates. The algorithms and tolls used to accomplish this task are common challenges in the field of computer vision, and any of a number of algorithms could be employed to recognize and define the boundary of a surface given an intensity map.

0110In some embodiments, the tools used for these tasks are a robust edge detection algorithm, various morphology operations, and association of high intensity regions of the intensity map with potential boundaries. These tools are described in detail below.

0111In general, an image is given by <i>A</i>=<i>f(m,n),</i> where pixel <i>(m,n)</i> is located at row <i>m</i> and column <i>n</i> of the image. Both <i>m</i> and <i>n</i> are elements of the set of integers, <i>Z. f</i> is a function outputting a real number (an element of the set <i>R</i>) such that <i>f:Z²</i>→<i>R.</i> If we restrict <i>f(m,n)</i> to values of either 0 or 1, then the image <i>f(m,n)</i> is referred to as a binary image.

0112The following portion of the specification defines some basic set operations on images, which form the groundwork for higher level operations defined in subsequent subsections. We first introduce the translation and reflection operations. For an image <i>A</i>, the translation of <i>A</i> by <i>x</i>=<i>(p,q)</i> is represented by <i>A</i> and is defined as <maths id="math0032" num="(Equation 4)"><math display="block"><msub><mi>A</mi><mi>x</mi></msub><mo>=</mo><mfenced open="{" close="}"><mi>f</mi><mfenced><mi>m</mi><mo>+</mo><mi>p</mi><mo>,</mo><mi>n</mi><mo>+</mo><mi>q</mi></mfenced><mrow><mo>|</mo><mrow><mi>f</mi><mfenced><mi>m</mi><mi>n</mi></mfenced><mi mathvariant="normal">ϵ</mi><mi>A</mi></mrow></mrow></mfenced></math><img file="EP3553561B1_D0032.tif" /></maths>

0113The reflection operation, denoted <i>Â</i>, is given as <maths id="math0033" num="(Equation 5)"><math display="block"><mover accent="true"><mi mathvariant="normal">A</mi><mo>^</mo></mover><mo>=</mo><mfenced open="{" close="}"><mi>f</mi><mfenced><mo>−</mo><mi>p</mi><mo>,</mo><mo>−</mo><mi>q</mi></mfenced><mrow><mo>|</mo><mrow><mi>f</mi><mfenced><mi>p</mi><mi>q</mi></mfenced><mspace width="1ex" /><mi mathvariant="normal">ϵ</mi><mspace width="1ex" /><mi>A</mi></mrow></mrow></mfenced></math><img file="EP3553561B1_D0033.tif" /></maths>

0114Furthermore, introduced are the fundamental union and intersection operators, represented by ∪ and ∩ respectively. The union of two images <i>A</i>=<i>f(m,n)</i> and <i>B</i>=<i>g(m,n)</i> is defined as <maths id="math0034" num="(Equation 6)"><math display="block"><mi>A</mi><mo>∪</mo><mi>B</mi><mo>=</mo><mfenced open="{" close="}"><mi mathvariant="italic">max</mi><mfenced><mi>f</mi><mfenced><mi>m</mi><mi>n</mi></mfenced><mo>,</mo><mi>g</mi><mfenced><mi>m</mi><mi>n</mi></mfenced></mfenced><mrow><mo>|</mo><mrow><mi>f</mi><mfenced><mi>m</mi><mi>n</mi></mfenced><mspace width="1ex" /><mi mathvariant="normal">ϵ</mi><mspace width="1ex" /><mi>A</mi><mo>,</mo><mi>g</mi><mfenced><mi>m</mi><mi>n</mi></mfenced><mspace width="1ex" /><mi mathvariant="normal">ϵ</mi><mspace width="1ex" /><mi>B</mi></mrow></mrow></mfenced></math><img file="EP3553561B1_D0034.tif" /></maths> whereas the intersection of the two images <i>A</i> and <i>B</i> is defined as <maths id="math0035" num="(Equation 7)"><math display="block"><mi>A</mi><mo>⋂</mo><mi>B</mi><mo>=</mo><mfenced open="{" close="}"><mi mathvariant="italic">min</mi><mfenced><mi>f</mi><mfenced><mi>m</mi><mi>n</mi></mfenced><mo>,</mo><mi>g</mi><mfenced><mi>m</mi><mi>n</mi></mfenced></mfenced><mrow><mo>|</mo><mrow><mi>f</mi><mfenced><mi>m</mi><mi>n</mi></mfenced><mspace width="1ex" /><mi mathvariant="normal">ϵ</mi><mspace width="1ex" /><mi>A</mi><mo>,</mo><mi>g</mi><mfenced><mi>m</mi><mi>n</mi></mfenced><mspace width="1ex" /><mi mathvariant="normal">ϵ</mi><mspace width="1ex" /><mi>B</mi></mrow></mrow></mfenced></math><img file="EP3553561B1_D0035.tif" /></maths>

0115Regions of adjacent images which share the same values are identified as connected components. While there are different measures of connectivity, the present described example embodiment considers pixels which are 8-connected. If any two pixels are adjacent to one another (including diagonally adjacent), then they are considered 8-connected.

0116Different algorithms exist for labelling connected components of an image. In general, given an image <i>A</i>, a new image <i>B</i> can be defined such that the values of its pixels are the labels of the connected components in image <i>A.</i> For a given image of interest, <i>A</i>, a structuring set <i>8,</i> called the 'structuring element', is an image which, along with either a dilation or erosion operation, is used to modify the image of interest.

0117The dilation operation is defined as <maths id="math0036" num="(Equation 8)"><math display="block"><mi>A</mi><mo>⊕</mo><mi>B</mi><mo>=</mo><mfenced open="{" close="}"><mi>x</mi><mrow><mo>|</mo><msub><mover accent="true"><mi>B</mi><mo>^</mo></mover><mi>x</mi></msub></mrow><mo>⋂</mo><mi>A</mi><mo>≠</mo><mn>0</mn><mo>,</mo><mi>x</mi><mo>=</mo><mfenced><mi>m</mi><mi>n</mi></mfenced><mo>,</mo><mi>m</mi><mo>,</mo><mi>n</mi><mo>∈</mo><mi>Z</mi></mfenced></math><img file="EP3553561B1_D0036.tif" /></maths>

0118Effectively, the dilation operation enlarges image <i>A</i> by reflecting <i>B</i> about its origin and then shifting it by <i>x</i>. Erosion on the other hand is defined as <maths id="math0037" num="(Equation 9)"><math display="block"><mi>A</mi><mo>⊖</mo><mi>B</mi><mo>=</mo><mfenced open="{" close="}"><mi>x</mi><mrow><mo>|</mo><msub><mi>B</mi><mi>x</mi></msub></mrow><mo>⋂</mo><mi>A</mi><mo>=</mo><mi>B</mi><mo>,</mo><mi>x</mi><mo>=</mo><mfenced><mi>m</mi><mi>n</mi></mfenced><mo>,</mo><mi>m</mi><mo>,</mo><mi>n</mi><mo>∈</mo><mi>Z</mi></mfenced></math><img file="EP3553561B1_D0037.tif" /></maths>

0119Erosion preserves 1's in <i>A</i> which when <i>B</i> is translated by <i>x</i> and is intersected with <i>A</i>, equals <i>B.</i> Erosion has the effect of trimming boundary 1's from an image given an image <i>B</i> which is centered on the origin.

0120<figref idref="f0014">Figures 24 and 25</figref> illustrate the dilation and erosion operations, respectively, performed on the image given in <figref idref="f0014">Figure 22</figref>. The structuring element used in the dilation and erosion operations is a 5 x 5 matrix, centered on the origin, given in <figref idref="f0014">Figure 23</figref>.

0121Edge detection can be performed using an edge detection algorithm as known in the art. One of the most popular, robust, and versatile edge detectors used is the Canny Edge detector. It is described in detail below.

0122The first step in Canny edge detection is to smooth the input image. The purpose of this smoothing is to reduce noise and unwanted details and textures. The smoothing process is performed via convolution of the image with a two dimensional Gaussian function. Where the image is given by <i>f(m,n),</i> and the Gaussian is given by <i>Gσ(m,n),</i> the convolution is given by <maths id="math0038" num="(Equation 10)"><math display="block"><mi>g</mi><mfenced><mi>m</mi><mi>n</mi></mfenced><mo>=</mo><mi>f</mi><mfenced><mi>m</mi><mi>n</mi></mfenced><mo>*</mo><mi mathvariant="italic">Gσ</mi><mfenced><mi>m</mi><mi>n</mi></mfenced></math><img file="EP3553561B1_D0038.tif" /></maths>

0123The particular two dimensional Gaussian used here is of the form <maths id="math0039" num="(Equation 11)"><math display="block"><mi mathvariant="italic">Gσ</mi><mfenced><mi>m</mi><mi>n</mi></mfenced><mo>=</mo><msup><mi mathvariant="italic">Ae</mi><mrow><mo>−</mo><mfenced><mfrac><mi>m</mi><mrow><mn>2</mn><msubsup><mi>σ</mi><mi>m</mi><mn>2</mn></msubsup></mrow></mfrac><mo>+</mo><mfrac><mi>n</mi><mrow><mn>2</mn><msubsup><mi>σ</mi><mi>n</mi><mn>2</mn></msubsup></mrow></mfrac></mfenced></mrow></msup></math><img file="EP3553561B1_D0039.tif" /></maths> where <i>σ_m</i> and <i>σ_n</i> are variances of the Gaussian in the vertical and horizontal directions respectively. In traditional implementation of Canny edge detection, <i>σ_m</i>=<i>σ_n</i>, yielding a simpler form of Equation 11. However, in the present described example embodiment, spacing between vertical pixels and horizontal pixels is not necessarily the same. For example, pixels adjacent to each other in the vertical direction may represent locations 0.01mm apart, while pixels adjacent to each other in the horizontal direction may represent locations 0.015mm apart. Therefore, if a Gaussian is used in the convolution process which is circular as opposed to elliptical in shape (not stretched in either the vertical or horizontal directions) with respect to actual locations of pixels, <i>σ_m</i> will not be equal to <i>σ_n</i> if spacing between vertical pixels and horizontal pixels is different. If the spacing between pixels in the vertical direction is a factor if c times as much as spacing between pixels in the horizontal direction, then <i>σ_m</i> = <i>σ_n</i>/<i>c.</i>

0124The second step in the edge detection process is to calculate the gradient of the convoluted image, <i>g(m,n).</i> Since <i>g(m,n)</i> is a discrete function, and generally non-analytic, operators have been developed to approximate the gradient of <i>g(m,n),</i> V <i>g(m,n).</i> The details of these operators are not described here. Utilization of any of the above mentioned operators produce <i>g_m(m,n)</i> and <i>g_n(m,n)</i> which are two images containing vertical and horizontal gradient approximations. Thus the magnitude and direction of ∇ <i>g(m,n)</i> may be represented as <i>M(m,n)</i> and <i>θ(m,n),</i> where <maths id="math0040" num="(Equation 12)"><math display="block"><mi>M</mi><mfenced><mi>m</mi><mi>n</mi></mfenced><mo>=</mo><msqrt><mrow><msub><mi>g</mi><mi>m</mi></msub><mfenced><mi>m</mi><mi>n</mi></mfenced><mo>+</mo><msub><mi>g</mi><mi>n</mi></msub><mfenced><mi>m</mi><mi>n</mi></mfenced></mrow></msqrt></math><img file="EP3553561B1_D0040.tif" /></maths><maths id="math0041" num="(Equation 13)"><math display="block"><mi>θ</mi><mfenced><mi>m</mi><mi>n</mi></mfenced><mo>=</mo><mi mathvariant="italic">arctan</mi><mfenced><mfrac><mrow><msub><mi>g</mi><mi>m</mi></msub><mfenced><mi>m</mi><mi>n</mi></mfenced></mrow><mrow><msub><mi>g</mi><mi>n</mi></msub><mfenced><mi>m</mi><mi>n</mi></mfenced></mrow></mfrac></mfenced></math><img file="EP3553561B1_D0041.tif" /></maths>

0125Step three is to calculate those pixels in <i>M(m,n)</i> which are local maxima. To do this, for every pixel <i>(m,n),</i> two most neighbouring pixels perpendicular to the direction <i>θ(m,n)</i> are considered. If the value of <i>M</i> at these pixels are both less than that of <i>M(m,n),</i> then <i>(m,n)</i> is considered to be a local maximum. The set of all local maxima in <i>M(m,n)</i> is denoted as <i>M̃</i>(m,n).

0126Step four is to define two binary images as functions of two respective real numbers, <i>τ</i>₁ and <i>τ</i>₂, respectively, where 0<<i>τ</i>₁<<i>τ</i>₂<<i>τ'</i>. <i>τ'</i> is equal to the maximum value of <i>M(m,n). M̃τ</i>¹<i>(m,n)</i> and <i>M̃τ</i>²<i>(m,n)</i> are defined as being images containing values of <i>M̃(m,n)</i> greater than <i>τ</i>₁ and <i>τ</i>₂ respectively. Thus, <maths id="math0042"><math display="block"><mover accent="true"><mi>M</mi><mo>˜</mo></mover><msup><mi>τ</mi><mn>1</mn></msup><mfenced><mi>m</mi><mi>n</mi></mfenced><mo>=</mo><mfenced open="[" close="]"><mfenced><mn>1</mn><mspace width="1ex" /><mi>if</mi><mspace width="1ex" /><mover accent="true"><mi>M</mi><mo>˜</mo></mover><mfenced><mi>m</mi><mi>n</mi></mfenced><mo>></mo><msub><mi>τ</mi><mi>i</mi></msub></mfenced><mo>,</mo><mfenced><mn>0</mn><mspace width="1ex" /><mi>otherwise</mi></mfenced><mo>,</mo><mspace width="1ex" /><mi>i</mi><mo>=</mo><mn>1,2</mn></mfenced></math><img file="EP3553561B1_D0042.tif" /></maths>

0127The last step (step five) is an iterative process. <i>E_j</i> is a Boolean image containing edges, output at iteration <i>j</i> of the process. <i>F_j</i> is the intersection of those pixels adjacent (in an 8-connected fashion) to pixels with value 1 in <i>E_j</i> and pixels with value 1 in <i>M̃τ</i>¹<i>(m,n).</i>

0128Set <i>j</i>=1 <i>E₁</i> as <i>M̃τ</i>²<i>(m,n).</i> The iterative computation begins. If the union of <i>E_j</i> and <i>F_j</i> is equal to <i>E_j,</i> then the algorithm terminates, outputting <i>E_j</i> as the canny detected edge. If not <i>E_j</i> is defined as the union of <i>E_j</i> and <i>F_j, and j</i> is incremented. The algorithm then loops, testing the equality of <i>E_j</i> and <i>E_j</i> U <i>F_j,</i> until they are equal. The algorithm outputs <i>E_j</i> when it is equal to the union of <i>E_j</i> U <i>F_j.</i>

0129This can be represented in pseudo-code as follows: <img file="EP3553561B1_D0043.tif" />

0130<figref idref="f0015">Figure 26(a)</figref> shows an image of a set of coins, while <figref idref="f0015">Figure 26(b)</figref> shows this image processed with Canny edge detection.

0131A thinning algorithm may also be used in boundary detection. Reducing connected components in an image to their thin-line representation has a number of useful applications such as data compression, simple structural analysis, and elimination of contour distortion. One iterative thinning used in an example embodiment is presented here. The ∧, ∨, and ^- symbols represent and, or, and not operations respectively.

0132Given an image <i>A</i>=<i>f(m,n),</i> the first step in the thinning algorithm is to divide it into two subfields in a checkerboard pattern. The algorithm then uses a parallel approach where the following two sub-iterations are executed in parallel repeatedly, until they have no effect on the image.

0133The first sub-iteration is to delete pixel <i>p</i> in the checkerboard if and only if conditions G₁, G₂, and G₃ are all satisfied.

0134The second sub-iteration is to delete pixel <i>p</i> in the checkerboard if and only if conditions G₁, G₂, and G₃' are all satisfied. Conditions G₁, G₂, G₃, and G₃' are listed as follows: <ul id="ul0003" list-style="none"><li>Condition G₁ is: <maths id="math0043"><math display="block"><msub><mi>X</mi><mi>H</mi></msub><mfenced><mi>p</mi></mfenced><mo>=</mo><mn>1</mn></math><img file="EP3553561B1_D0044.tif" /></maths> where <maths id="math0044"><math display="block"><msub><mi>X</mi><mi>H</mi></msub><mfenced><mi>p</mi></mfenced><mo>=</mo><mstyle displaystyle="true"><msubsup><mo>∑</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mn>4</mn></msubsup><msub><mi>b</mi><mi>i</mi></msub></mstyle></math><img file="EP3553561B1_D0045.tif" /></maths> and <maths id="math0045"><math display="block"><msub><mi>b</mi><mi>i</mi></msub><mo>=</mo><mfenced open="[" close="]"><mo>(</mo><mrow><mn>1</mn><mspace width="1ex" /><mi>if</mi><mspace width="1ex" /><msub><mi>x</mi><mrow><mn>2</mn><mi>i</mi><mo>−</mo><mn>1</mn></mrow></msub><mo>=</mo><mn>0</mn><mo>∧</mo><mfenced><msub><mi>x</mi><mrow><mn>2</mn><mi>i</mi><mo>−</mo><mn>1</mn></mrow></msub><mo>=</mo><mn>1</mn><mo>∨</mo><msub><mi>x</mi><mrow><mn>2</mn><mi>i</mi><mo>+</mo><mn>1</mn></mrow></msub><mo>=</mo><mn>1</mn></mfenced><mo>,</mo><mspace width="1ex" /><mfenced><mi mathvariant="normal">0</mi><mspace width="1ex" /><mi>otherwise</mi></mfenced></mrow></mfenced></math><img file="EP3553561B1_D0046.tif" /></maths> x₁, x₂, ..., x₈ are the values of the eight neighbours of <i>p</i> respectively beginning with the east neighbour and numbered counter-clockwise until the south-east neighbour is reached.</li><li>Condition G₂ is: <maths id="math0046"><math display="block"><mn>2</mn><mo>≤</mo><mi>min</mi><mfenced open="{" close="}"><msub><mi mathvariant="normal">n</mi><mn>1</mn></msub><mfenced><mi mathvariant="normal">p</mi></mfenced><mo>,</mo><msub><mi mathvariant="normal">n</mi><mn>2</mn></msub><mfenced><mi mathvariant="normal">p</mi></mfenced></mfenced><mo>≤</mo><mn>3</mn></math><img file="EP3553561B1_D0047.tif" /></maths> where <maths id="math0047"><math display="block"><msub><mi mathvariant="normal">n</mi><mn>1</mn></msub><mfenced><mi mathvariant="normal">p</mi></mfenced><mo>=</mo><mstyle displaystyle="true"><msubsup><mo>∑</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mn>4</mn></msubsup><mrow><msub><mi>x</mi><mrow><mn>2</mn><mi>k</mi><mo>−</mo><mn>1</mn></mrow></msub><mo>∨</mo><msub><mi>x</mi><mrow><mn>2</mn><mi>k</mi></mrow></msub></mrow></mstyle></math><img file="EP3553561B1_D0048.tif" /></maths><maths id="math0048"><math display="block"><msub><mi mathvariant="normal">n</mi><mn>2</mn></msub><mfenced><mi mathvariant="normal">p</mi></mfenced><mo>=</mo><mstyle displaystyle="true"><msubsup><mo>∑</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mn>4</mn></msubsup><mrow><msub><mi>x</mi><mrow><mn>2</mn><mi>k</mi></mrow></msub><mo>∨</mo><msub><mi>x</mi><mrow><mn>2</mn><mi>k</mi><mo>+</mo><mn>1</mn></mrow></msub></mrow></mstyle></math><img file="EP3553561B1_D0049.tif" /></maths></li><li>Condition G₃ is: <maths id="math0049"><math display="block"><mfenced><msub><mi>x</mi><mn>2</mn></msub><mo>∨</mo><msub><mi>x</mi><mn>3</mn></msub><mo>∨</mo><mover accent="true"><msub><mi>x</mi><mn>8</mn></msub><mo>‾</mo></mover></mfenced><mo>∧</mo><msub><mi>x</mi><mn>1</mn></msub><mo>=</mo><mn>0</mn></math><img file="EP3553561B1_D0050.tif" /></maths></li><li>Condition G₃' is: <maths id="math0050"><math display="block"><mfenced><msub><mi>x</mi><mn>6</mn></msub><mo>∨</mo><msub><mi>x</mi><mn>7</mn></msub><mo>∨</mo><mover accent="true"><msub><mi>x</mi><mn>4</mn></msub><mo>‾</mo></mover></mfenced><mo>∧</mo><msub><mi>x</mi><mn>5</mn></msub><mo>=</mo><mn>0</mn></math><img file="EP3553561B1_D0051.tif" /></maths></li></ul>

0135<figref idref="f0015">Figure 27(a)</figref> shows a binary image without thinning. <figref idref="f0015">Figure 27(b)</figref> shows the same image after a thinning algorithm has been applied.

0136An algorithm may also be applied to test if a pixel is a junction point of lines in a binary image. Given a 3x3 neighbourhood of a pixel <i>p</i>, <i>p</i> is a junction of lines if and only if when traversing the perimeter of <i>p</i> the number of transitions between 0 and 1 is either 6 or 8. <figref idref="f0016">Figure 28(a)</figref> gives a junction of lines and <figref idref="f0016">Figure 28(b)</figref> zooms in on this junction and illustrates its 3 x 3 neighbourhood. It is clear from <figref idref="f0016">Figure 28(a)</figref> that the junction has 6 transitions between 0 and 1. Thus, the criterion for junction identification is satisfied.

0137The algorithm employed to extract OD or ID boundary coordinates from an OD or ID intensity map in example embodiments may employ all the tools introduced in the previous subsection and tailor their use to the specific domain of ultrasound scanning.

0138The first step to extraction of the true OD/ID boundary from the OD/ID intensity map involves identifying canny edge detected boundaries near and above (on the side near to the array 200) maximum intensity pixels in the vertical direction of the intensity map. Given an intensity map, the true boundary may be extracted. The first question that arises is where the precise coordinates of the boundary lie, given the neighbourhood of and around high intensity pixels in an intensity map. It can be observed from a given intensity map (such as <figref idref="f0017">Figure 30</figref>) that there is a relatively wide region of high intensity pixilation where the true boundary could lie. To answer this question, the method of calculation of distances from A-Scans is extended to intensity maps. To calculate the distance an ultrasonic wave has travelled in a medium from an A-Scan, the time from transmitter excitation to the leading edge of the received wave packet is multiplied by the speed of sound in the medium. The travel time (in samples at a particular digitization frequency) from transmitter excitation to the leading edge of a received wave packet is shown by numeral 2902 in <figref idref="f0016">Figure 29</figref>. The formation of intensity maps involves the mapping of values of analytic time-domain signals (A-Scans along with their Hilbert transforms) to points in the intensity map, per Equation 2. The leading edges of wave packets of A-Scans are mapped to the leading edge of intensity map boundaries where the leading edge of intensity map boundaries are defined as the edges on the near side of the probe array 200. Thus, the OD/ID boundary is derived from leading edges in the intensity map boundaries.

0139The second question that arises is how to programmatically identify the true OD/ID boundary with the potential of imaging aberrations (features appearing due to mechanisms such as grating lobes or multiple back-wall reflections) present in the intensity maps. In <figref idref="f0017">Figure 30</figref>, of the edges in the image, the true boundary is of the edges nearest to the regions of high intensity, contrasting the edges of aberrations due to grating or back wall reflections. Taking any vertical slice of the image, the true boundary will intersect this vertical slice just above where it will intersect the maximum intensity pixel in that slice. These observations motivate the first step of the boundary detection algorithm, isolation of the true boundary edges through identification of edges near and above the highest intensity pixels in the vertical direction. For a given intensity map such as the one given in <figref idref="f0017">Figure 30</figref>, the edges of the intensity map and the pixels of maximum intensity in the vertical direction are plotted in <figref idref="f0017">Figure 31</figref>. The true boundary edges are those that are near and above the high intensity pixels and are plotted in <figref idref="f0017">Figure 32</figref>.

0140Identification of the OD/ID boundary solely through the methods employed described above may yield results which can be improved upon. Additional techniques may be used in some embodiments to eliminate error in the boundary defined via the methods described earlier above.

0141Aberrations can have higher intensity content than the true boundary, which can lead to false identification of the true OD/ID boundary. If these aberrations are small in size, however, operations of dilation and thinning, along with comparison of connected component sizes, can be employed to remove erroneous edges from the boundary definition extracted from the algorithm previously described. Elimination of erroneous edges at this stage in the boundary detection algorithm may be performed via three steps: dilation of the candidate boundary, thinning of the dilated boundary, and trimming of small connected components.

0142<figref idref="f0018">Figure 33 to Figure 36</figref> give the boundary at various stages in the execution of an example algorithm for eliminating erroneous edges. The original boundary output from the algorithm given above is depicted in <figref idref="f0018">Figure 33</figref>. Due to high intensity content in the aberration below the true boundary, the output boundary contains error. This is evident in that a small region of the output boundary appears above the aberration. <figref idref="f0018">Figure 34</figref> shows the boundary after the dilation operation is performed with a rectangular structuring element. The small gap in the location where the boundary should lie (<figref idref="f0018">Figure 33</figref>) has been removed. Thinning of the boundary reduces it back to its desired width. The thinned boundary is depicted in <figref idref="f0018">Figure 35</figref>. Finally, a comparison of boundary pixels is performed. If there is more than one boundary pixel intersecting any vertical slice of the intensity map, only the pixel belonging to the connected component of largest size (of the pixels under consideration) is preserved. The rest are removed. This serves to remove the erroneous boundary pixels. <figref idref="f0018">Figure 36</figref> gives the boundary with erroneous pixels removed.

0143The dilation and thinning operations can introduce line junctions as previously defined. At this point in the boundary detection algorithm, junctions may be removed, as set out above, until no junctions exist in the defined boundary. This serves the purpose of preparing the boundary for removal of more erroneous edges, performed via the algorithms described below.

0144<figref idref="f0019">Figure 37</figref> gives an ID intensity map image. The edge detection, dilation and thinning algorithms give the boundary given in <figref idref="f0019">Figure 38</figref>. The junction area 3802 of the boundary encircled in <figref idref="f0019">Figure 38</figref> reveals the area on and around the boundary to be removed as erroneous. This junction area 3802 is shown magnified in <figref idref="f0019">Figure 39. Figure 40</figref> illustrates the junction area 3802 of the boundary with the junction removed.

0145The next phase of the edge detection process is a removal of bottom portions of connected components. Removal of bottom portions of connected components consists of removal of lower connected component pixels. A lower connected component pixel is defined here as a pixel in some connected component having the same x-component as some other pixel(s) in the same connected component, but with greater z-component values (recall pixels with greater z-component values appear lower on intensity maps. <figref idref="f0020">Figure 41</figref> shows the boundary in <figref idref="f0019">Figure 40</figref> with bottom portions of connected components removed.

0146Small connected components of erroneous edges can still remain in the intensity map. Removal of these edges increases the accuracy of the defined boundary. For every vertical strip of the intensity map, if more than one connected component intersects this slice, then only the largest one is preserved in the intensity map. <figref idref="f0020">Figure 42</figref> illustrates the edge boundary with small connected components removed. This image is the most accurate of any so far in the identification of the true edge boundary of the intensity map depicted in <figref idref="f0019">Figure 37</figref>.

0147Intensity maps are formed from ultrasonic raw data. The decrease in usable ultrasonic raw data near the ends of boundaries will result in a thinning of high intensity content. Using the Canny edge detection process to determine the true boundary in these areas will result in poor determination of the true OD/ID boundary as the edges output curve around the high intensity content, where in fact the true boundary may be flat. <figref idref="f0020">Figure 43</figref> illustrates the edge output from the Canny edge detector 4302 in red. Indeed the curvature defined in the area 4304 where the intensity tapers off is very high and does not do a good job defining the true boundary. A better way to define the boundary in areas where the intensity content tapers off is to use the curvature of the maximum intensity pixels (in the vertical dimension) to approximate the curvature of the true boundary. <figref idref="f0020">Figure 43</figref> provides a good example where the maximum intensity pixels 4306 (in black) approximate the shape of the true boundary. The true boundary in this case is a flat plate.

0148The curvature of the boundary near the end of detected edges is thus approximated by the curvature of the maximum intensity pixels in the vertical strips of the intensity map containing both the detected edges and pixels of maximum intensity. This approximation is performed only if, when dilated, the pixels of maximum intensity are all connected to each other (i.e. they belong to the same connected component). The coordinates of the maximum intensity pixels are then translated vertically such that the pixel of maximum intensity near the side of the detected edge is connected horizontally adjacent to it.

0149In areas between the definitions of the true boundary, where no true boundary edges are defined, the true boundary can be interpolated such that true boundary edges are connected together. This is the last step of the boundary definition process. <figref idref="f0021">Figure 44</figref> illustrates how the boundary can be interpolated. Areas between the boundary segments 4402 extracted from the edge detection process are connected with straight lines 4404.

0150The sequential operation of an example algorithm for boundary recognition and definition is shown in the flow chart of <figref idref="f0022">Figure 45</figref>. The intensity map 4502 is subjected to a series of image processing algorithms, including edge detection 4504, dilation of the edge 4506, thinning of the dilated edge 4508, trimming of erroneous pixels 4510, edge junction removal 4512, removal of bottom portions of connected components 4514, removal of small connected components 4516, and horizontal end approximation 4518, finally outputting a total boundary 4520 that includes the calculated and interpolated result of the previous operations.

0151A comprehensive flow chart showing the entire boundary recognition data flow is given in <figref idref="f0043">Figure 69</figref>. Detailed descriptions of boundary recognition parameters used in an example embodiment of the system are set out in Table A4 at the end of the Description. Detailed descriptions of boundary recognition functions used in an example embodiment of the system are set out in Table A5 at the end of the Description.

0152The Interior Focus Method (IFM), like the Shifting Aperture Focus Method (SFM), is an algorithm whose purpose is to output an intensity map (an image where intensities are assigned to inspection coordinates of interest) given full matrix capture raw data and a smooth approximation of the OD boundary surface.

0153Once the OD intensity map has been determined and an OD surface model determined from the high-intensity OD regions 32, Fermat's Principle can be applied to model details of the inner pipe surface as well. Fermat's Principle is one of the tools by which total focusing is extended beyond the media interface (e.g. curve <i>K</i> 10 in <figref idref="f0005">Figure 10</figref>), a technique which will be referred to as the Interior Focusing Method (IFM).

0154<b>Theorem 1</b> (the Modern Version of Fermat's Principle) states that the path that sound takes from a point source emitter at point p, to another point q, is such that the time taken to traverse the path is a stationary value (i.e. a minimum, maximum, or inflection point).

0155Extending the intensity equation shown above (Equation 2) to the area beyond a media interface involves substituting the times <i>t</i> in the equation to those equalling travel path times whose paths experience refraction at the media interface. Fermat's Principle (Theorem 1) can be applied in solving for these times.

0156Since <i>K</i> is piecewise-smooth, <i>T_ir(K)</i> has a finite number of stationary values. These stationary values are denoted as <maths id="math0051"><math display="inline"><msup><msub><mi>t</mi><mi mathvariant="italic">ir</mi></msub><mrow><mi>K</mi><mo>′</mo></mrow></msup><mo>,</mo><mspace width="1ex" /><msup><msub><mi>t</mi><mi mathvariant="italic">ir</mi></msub><mrow><mi>K</mi><mo>"</mo></mrow></msup><mo>,</mo><mspace width="1ex" /><msup><msub><mi>t</mi><mi mathvariant="italic">ir</mi></msub><mrow><mi>K</mi><mo>‴</mo></mrow></msup><mo>,</mo><mspace width="1ex" /><mo>…</mo></math><img file="EP3553561B1_D0052.tif" /></maths> and the set of these stationary values as <maths id="math0052"><math display="inline"><msup><msub><mi>T</mi><mi mathvariant="italic">ir</mi></msub><mi>K</mi></msup></math><img file="EP3553561B1_D0053.tif" /></maths> where: <maths id="math0053" num="(Equation 4)"><math display="block"><msubsup><mi>T</mi><mi mathvariant="italic">ir</mi><mi>K</mi></msubsup><mo>=</mo><mfenced open="{" close="}"><msubsup><mi>t</mi><mi mathvariant="italic">ir</mi><mrow><mi>K</mi><mo>′</mo></mrow></msubsup><mo>,</mo><msubsup><mi>t</mi><mi mathvariant="italic">ir</mi><mrow><mi>K</mi><mo>"</mo></mrow></msubsup><mo>,</mo><msubsup><mi>t</mi><mi mathvariant="italic">ir</mi><mrow><mi>K</mi><mo>‴</mo></mrow></msubsup><mo>,</mo><mo>…</mo></mfenced></math><img file="EP3553561B1_D0054.tif" /></maths>

0157Sets <maths id="math0054"><math display="inline"><msup><msub><mi>T</mi><mi mathvariant="italic">ir</mi></msub><mi>K</mi></msup></math><img file="EP3553561B1_D0055.tif" /></maths> and <maths id="math0055"><math display="inline"><msup><msub><mi>T</mi><mi mathvariant="italic">jr</mi></msub><mi>K</mi></msup></math><img file="EP3553561B1_D0056.tif" /></maths> are defined where i denotes element i 210 with position vector <i>e_(i)</i> and <i>j</i> denotes element <i>j</i> 212 with position vector <i>e_j</i>. The set of summations of all possible combinations between elements of Ti_r^K and elements of T_jr^K are defined as T_ij^K (<b>r</b>) where: <maths id="math0056" num="(Equation 5)"><math display="block"><msubsup><mi>T</mi><mi mathvariant="italic">ir</mi><mi>K</mi></msubsup><mfenced><mi mathvariant="bold-italic">r</mi></mfenced><mo>=</mo><mfenced open="{" close="}"><msub><mi>t</mi><mi mathvariant="italic">ir</mi></msub><mo>+</mo><mrow><msub><mi>t</mi><mi mathvariant="italic">jr</mi></msub><mo>|</mo></mrow><msub><mi>t</mi><mi mathvariant="italic">ir</mi></msub><mo>∈</mo><msubsup><mi>T</mi><mi mathvariant="italic">ir</mi><mi>K</mi></msubsup><mo>,</mo><msub><mi>t</mi><mi mathvariant="italic">jr</mi></msub><mo>∈</mo><msubsup><mi>T</mi><mi mathvariant="italic">jr</mi><mi>K</mi></msubsup></mfenced></math><img file="EP3553561B1_D0057.tif" /></maths>

0158By Fermat's Principle (Theorem 1), the above set contains travel path times which can be substituted in Equation 2, to extend the total focusing method to imaging areas beyond the interface 10 between media 6, 8. If the transmitter and receiver pair i 210 and j 212 are in the first medium 6, <i>r</i> is in the second medium 8, and <i>K</i> 10 is the curve which divides the first medium 6 and the second medium 8, the intensity of the image at <i>r</i> can be written: <maths id="math0057" num="(Equation 6)"><math display="block"><mi>I</mi><mfenced><mi mathvariant="bold">r</mi><mi>a</mi></mfenced><mo>=</mo><mfenced open="|" close="|"><mstyle displaystyle="true"><munder><mo>∑</mo><mrow><mi>i</mi><mo>,</mo><mi>j</mi><mo>∈</mo><mi>a</mi></mrow></munder><mstyle displaystyle="true"><munder><mo>∑</mo><mrow><mi>t</mi><mo>′</mo><mo>∈</mo><msubsup><mi>T</mi><mi mathvariant="italic">ij</mi><mi>K</mi></msubsup><mfenced><mi mathvariant="bold">r</mi></mfenced></mrow></munder><mrow><msub><mi>g</mi><mrow><mfenced><mi>i</mi></mfenced><mi>j</mi></mrow></msub><mfenced><mi>t</mi><mo>′</mo></mfenced></mrow></mstyle></mstyle></mfenced></math><img file="EP3553561B1_D0058.tif" /></maths>

0159These results can be extended to three dimensional interfaces, as Fermat's Principle holds in three dimensions. Furthermore, the results may be extended for focusing through multiple media interfaces; <i>n</i> media interfaces will imply sound paths that will traverse <i>n</i>+<i>1</i> media. In this case Fermat's Principle can again be utilized to solve for the true travel paths. Some embodiments may employ these techniques to perform more complex analysis, such as scanning volumes containing more than two media or scanning in three dimensions at once rather than in a plane.

0160IFM is applied to create an intensity map of the ID scanning plane, such as the ID intensity map shown in <figref idref="f0007">Figure 13</figref>. After the ID intensity map 52 has been created, the steps of boundary recognition 1720 and boundary definition 1722 are applied using the same algorithms described above to create the ID boundary definition 1724.

0161Implementation of the IFM subroutine can become very computationally intensive if there are many coordinates for which corresponding intensities are evaluated. This is similar to the case for implementation of the SFM subroutine, the difference being the IFM routine is even more computationally taxing. The approach used to reduce the number of necessary computations in the IFM may be identical to the strategy employed in the SFM routine. Limiting the number of coordinates under consideration, while focusing at the appropriate density to meet inspection specifications, imposes the following focusing strategy. First, intensities of coordinates on a course grid may be calculated. The strategy that may be employed is to first calculate intensities of coordinates on a course grid.

0162At this point, depending on user defined parameters, the IFM subroutine may end, outputting the intensity map to the boundary detection subroutine, or proceed to define new coordinates for which to assign intensities. If the latter course is taken, newly defined coordinates will be positioned around coordinates with high intensities already assigned to them. The cutoff intensity for coordinates of which newly defined coordinates focus around is defined in some embodiments by the vector zoom Percentage (see Table A10). Once the new coordinates have been defined, again depending on user defined parameters, the SFM subroutine may exits, or may proceed to further focus around coordinates of high intensity. The process of identifying high intensity coordinates and then refocusing around them can be executed an arbitrary number of times, and may be specified by the user.

0163One cycle of the refocusing process is illustrated in <figref idref="f0011">Figure 19</figref>. Coordinates spaced coarsely apart are represented as either white or black circles. The white circles represent those coordinates whose respective intensities are below the cutoff intensity for which new coordinates are defined. Conversely, the black circles represent those coordinates whose respective intensities exceed the cutoff intensity for which new coordinates are defined. The gray circles represent newly defined coordinates around high intensity coordinates. If the gray circles are defined on iteration i+1 of the coordinate definition and intensity assigning process, dx(i)/dx(i+1)=dz(i)/dz(i+1)=4 for the example illustrated by <figref idref="f0011">Figure 19</figref>.

0164<figref idref="f0045">Figure 71</figref> provides a full data flow for the IFM process. Detailed descriptions of ID imaging parameters used in an example embodiment of the system are set out in Table A10 at the end of the Description. Detailed descriptions of ID imaging functions used in an example embodiment of the system are set out in Table A11 at the end of the Description.

0165The boundary output from the boundary recognition procedure described above may be susceptible to noise. Two main factors contribute to this noise. The first is due to the imaging of aberrations such as grating in the intensity map calculated via SFM. The second is due to the quantization of the intensity map grid. Errors (noise) can arise in a definition of a quantity due to a quantization of its state space.

0166To minimize the effects of noise in determining true boundary (either OD or ID) coordinates, standard filtering techniques may be employed. The boundary outputted from the boundary recognition algorithm may in some embodiments be filtered via two successive processes. The first filter used to eliminate noise from the input data may be a median filter of variable window size. This has the effect of removing high frequency noise from the input coordinates. The next filter used may be a Savitzky-Golay filter. This filter reduces quantization noise by approximating the input coordinates by an unweighted linear least-squares fit using a polynomial of a given degree. Utilization of polynomials of higher degrees makes it possible to smooth heavily while retaining data features of interest.

0167To initially eliminate noise from the calculated boundary, a median filter may be employed. Median filters implement a sliding window to a sequential data sequence, replacing the center value in the window with the median value of all the points within the window. Such a median filter may be chosen for its combination of good noise elimination coupled with its edge-preserving features. For both the OD algorithms and the ID algorithms, the window width of the algorithm may be configured by a user in some embodiments. In some embodiments, the window width may be set to 5mm by default.

0168Savitzky-Golay smoothing filters may also be used due to their superior performance in smoothing out a noisy signal whose frequency range (without noise) is large, as the corroded regions under a pipe weld cap may be seen to be with their potentially sharp edges. Savitsky-Golay filters perform a least-squares fit of a windowed set of consecutive data points to a polynomial and take the calculated central point of the fitted polynomial curve as the new smoothed data point. A set of convolution integers can be derived and used as weighting coefficients for the smoothing operation, a computationally-efficient process. This methodology is exactly equivalent to fitting the data to a specified polynomial. The smoothed data point (y_k)_s by the Savitsky-Golay algorithm is: <maths id="math0058"><math display="block"><msub><mfenced><msub><mi>y</mi><mi>k</mi></msub></mfenced><mi>s</mi></msub><mo>=</mo><mfrac><mstyle displaystyle="true"><msubsup><mo>∑</mo><mrow><mi>i</mi><mo>=</mo><mo>−</mo><mi>n</mi></mrow><mi>n</mi></msubsup><mrow><msub><mi>A</mi><mi>i</mi></msub><msub><mi>y</mi><mrow><mi>k</mi><mo>+</mo><mi>i</mi></mrow></msub></mrow></mstyle><mstyle displaystyle="true"><msubsup><mo>∑</mo><mrow><mi>i</mi><mo>=</mo><mo>−</mo><mi>n</mi></mrow><mi>n</mi></msubsup><msub><mi>A</mi><mi>i</mi></msub></mstyle></mfrac></math><img file="EP3553561B1_D0059.tif" /></maths> where A_i are the convolution integers, and the data window goes from -n to n.

0169In some embodiments, the user-adjustable parameters for the Savitsky-Golay filters used in the OD/ID smoothing are: <ol id="ol0001" compact="compact"><li>1. Smooth Width: The smoothing window width expressed in mm. The algorithm uses the appropriate "X Spacing" parameter to determine the window width in elements: <maths id="math0059"><math display="block"><msub><mi mathvariant="italic">Number</mi><mi mathvariant="italic">elements</mi></msub><mo>=</mo><mfrac><mrow><mi mathvariant="italic">Smooth</mi><mspace width="1ex" /><mi mathvariant="italic">Width</mi></mrow><mrow><mi>X</mi><mspace width="1ex" /><mi mathvariant="italic">Spacing</mi></mrow></mfrac></math><img file="EP3553561B1_D0060.tif" /></maths></li><li>2. Smooth Order: The polynomial order to be used in the Savitsky-Golay algorithm. It is generally smaller than the smoothing window number of elements. A typical value for Smooth Order is 4.</li></ol>

0170An example data flow for boundary definition by the above processes is shown in <figref idref="f0044">Figure 70</figref>. Detailed descriptions of boundary definition parameters used in an example embodiment of the system are set out in Table A6 at the end of the Description. Detailed descriptions of boundary definition functions used in an example embodiment of the system are set out in Table A7 at the end of the Description.

0171The interior focus method (IFM) calculates travel times from probe array elements to points under the outer diameter via Fermat's Principle. Due to a change in the refractive index at the OD, refraction of ultrasonic waves will occur at this interface. Very small, high frequency variations in the definition of the OD surface will invite a large number of travel path solutions between probe array elements and inspection points under the OD, intersecting small regions in the defined OD surface. These travel path solutions are due to either true high frequency variations in the OD surface, high frequency noise introduced into the OD definition, or a combination of both. The signal processing techniques that are described below serve to both to eliminate erroneous travel time solutions caused by high frequency noise content in the input OD boundary and to reduce IFM computation time.

0172While smoothing the defined OD surface can eliminate spurious travel path solutions due to high frequency noise, it is not obvious what effect this will have on travel path solutions induced by true high frequency variations in the OD surface. Eliminating high frequency content of the defined OD surface may reduce the number of travel paths with solutions very close to each other, but will generally preserve effectively distinct solutions. For example, a series of travel path solutions 4.52µs, 4.55µs, 4.48µs, 6.21µs, 6.18µs, 6.23µs could feasibly be reduced to 4.53µs and 6.22 µs. Programmatically, utilizing the latter travel time solutions may produce similar intensity maps output out of the IFM subroutine as those that can be produced utilizing the former travel time solutions, since travel times of interest are retained. Features of interest in the intensity map will be preserved, with the added benefit of reduced computation time of the IFM subroutine, due to a fewer number of computations involving travel time solutions being executed.

0173Thus, the smoothing operations may serve to eliminate the effects of noise with a beneficial side effect. Smoothing of the OD boundary serves firstly to eliminate erroneous travel time solutions due to high frequency noise, while reducing the number of travel time solutions to those that are effectively distinct. The smoothing methods employed on the boundary of the OD may include a median filter followed by a Savitzky-Golay filter. They may be the same filters used for the purpose of boundary definition. These filters have the effect of eliminating high frequency content while preserving data features of interest. Since the utilization of these filters for preparation of the OD boundary as input into IFM has objectives different than the utilization for boundary definition, filtering parameters used may be different. Larger window sizes may be used for the Savitzky-Golay filter in IFM OD boundary preparation than those used in OD boundary preparation, since preservation of high frequency content in the OD boundary is not desired.

0174Finally, the smoothed data may be downsampled. This has the primary effect of reducing computation time necessary to calculate travel times via the IFM subroutine, while retaining travel time solutions of interest.

0175An example data flow for IFM OD boundary preparation by the above processes is shown in <figref idref="f0041">Figure 67</figref>. Detailed descriptions of IFM OD boundary preparation parameters used in an example embodiment of the system are set out in Table A8 at the end of the Description. Detailed descriptions of IFM OD boundary preparation functions used in an example embodiment of the system are set out in Table A9 at the end of the Description.

0176Some example embodiments may be able to adjust the angle of projection of the ultrasound probe array 200 relative to the surface being scanned. For example, when scanning a pipe wall, some embodiments of the manipulator 100 may make use of an adjustable mirror 154 that reflects ultrasound projected from and reflected to the probe array 200. <figref idref="f0023">Figures 46 and 47</figref> show an example embodiment of the carrier 102 having an adjustable mirror 154. The ultrasound probe array 200 projects ultrasound from a linear array of elements arranged on its active face 152. Ultrasound waves are projected toward the mirror 154, whose angle can be adjusted by a mirror adjustment assembly 156 including a motor. By projecting ultrasound from one or more elements and sensing with one or more elements at different angles of reflection, it is possible to optimize the angle of the mirror to ensure that the ultrasound waves are being projected normal to the pipe surface. In some embodiments, this optimizing procedure may be performed prior to each frame of FMC data collection. The mirror adjustment assembly 156 may be controlled by the external controller via the probe data connectors 108.

0177In embodiments having a mirror 154, the transducer is positioned to fire tangentially instead of radially to the inspection surface. The mirror placement is designed to optimize the range of alignment accommodation while minimizing the added UT beam transit distance to the inspection surface.

0178The manipulator 100 may also be controlled by software to adapt the change in inspection step size given possible changes in mirror angle. The control of the mirror may be achieved in some such embodiments by pulsing a group of elements 202 and summing their response. The amplitude of the response is monitored after a series pre-defined step increments. Should the response level drop below a user defined threshold, the module controlling the manipulator 100 will pulse the user defined control aperture while driving the mirror 154 through a defined range of angles at a defined step size. The system sums the responses at each step and then drives the mirror to the angle corresponding to the peak amplitude.

0179The user interface used in such an embodiment may have a GUI where all relevant control variables are entered. The user may have the capability to change the tracking parameters during a data acquisition scan as to effect an optimal control region on the inspection surface.

0180Further details on optimizing the scanning angle using the adjustable mirror are provided below in the description of the data acquisition procedure as applied to pipe weld inspection.

0181Once the OD and ID surface profiles have been modeled, they can be used to determine structural characteristics of the pipe wall or other scanned object. For example, the OD profile 36 and ID profile 38 for a single transmit-receive cycle are combined in the graph shown in <figref idref="f0007">Figure 14</figref>. By juxtaposing these two profiles, information about the thinnest point 40 in the scanned portion of the pipe wall can be assessed and made visible to a user.

0182Thinning of feeder pipe material adjacent to and within the weld area is a concern for power generating stations and other industries. Thinning may be found around the entire circumference of the joint with some highly localized areas directly under the weld cap. The location of these areas may not be predictable due to the association with the weld root condition. <figref idref="f0024">Figure 48</figref> shows an example of a pipe wall exhibiting thinning.

0183A procedure is described below which provides a method and requirements for performing thickness measurements on heat transport feeder piping welds at nuclear power generating stations, using a Micropulse^™ ultrasonic inspection and data acquisition system. This example embodiment is illustrative, and the methods and devices described may be implemented in various embodiments and applied various contexts.

0184The inspection procedure described below is intended to detect localized and broadly based wall thinning in the regions adjacent to circumferential welded joints in feeder piping. This procedure applies to various pipe sizes.

0185As shown in <figref idref="f0024">Figure 49</figref>, the region of inspection is 20 mm on either side of the weld centre line. However, geometric limitations may constrain the ability to achieve this coverage. This inspection records all ultrasonic signals transmitted from each element to all possible receiving elements. Each element in the array is fired in sequence. Thus, for a 128 element linear array, each of the 128 elements that has fired has 128 receptions for a total of 16,384 separate A-scans. This process is repeated at each location as specified in the inspection sequence.

0186The data recorded is then submitted for analysis. The data analysis software produces OD and ID profiles of each inspection data set. The data sets are acquired at discrete circumferential locations. Together the series of profiles around the feeder constitutes the results of the inspection. A separate table records the minimum distance between the OD and ID profiles for each circumferential location, the axial location of the minimum in the slice and the circumferential position of the slice.

0187Special circumstances and exceptions to the procedure may exist. For example, hot spots or high local fields may create a situation where operators cannot access the work area. In addition, in the case of Grayloc^™ fitting to pipe welds the inspection coverage may be modified to collect a zone spanning 10 mm on the Grayloc^™ side of the weld centre line to 20 mm on the fitting side of the weld centre line, as shown in <figref idref="f0024">Figure 50</figref>. (in the context of this embodiment, a Grayloc^™ fitting refers to a conical pipe fitting with its narrow end matching the diameter of a straight pipe section.)

0188The following definitions are used within the context of pipe weld inspection as described below. <i>FMC</i>: Full Matrix Capture. Ultrasonic data collection strategy in which each element in the transducer is individually pulsed while all elements receive. This is repeated for each element in the transducer until all elements have been fired. This strategy creates a data array of <i>n</i> by <i>n</i> where <i>n</i> is the number of elements in the transducer. As a consequence the data files for a FMC inspection is significantly larger than for the equivalent conventional (e.g. phased-array) technique at the same resolution.

0189<i>Home Position</i>: Circumferential location of the inner rotating ring such that the manipulator may be safely unlocked from its closed position. Home position is indicated in both the VIM and the software as identified by the Hall sensor on the manipulator.

0190<i>Main diagonal</i>: A group of send-receive elements in the data collected using FMC where each transmitting element is its own receiver. The main diagonal view of the FMC data set is identical to the conventional linear electronic B scan. The main diagonal view is the default view of the FMC data B scan.

0191<i>Matrix</i>: Data structure created when using the FMC data collection strategy. If the columns of the matrix are assigned to identify the transmitting element, then the rows of the matrix correspond to the receiving elements. Each element of the array then corresponds to an A scan related to that transmitter receiver pair. For example: a combination of transmitting on element 17, receiving on element 32 would produce an A scan that would be located under the 17th column on the 32nd row of the FMC data matrix.

0192<i>Start Position</i>: Circumferential location, with respect to the manipulator, where the scan is initiated. The start location may correspond to the home location or can be offset from the home location.

0193<i>TFM</i>: Total Focus Method. Generic name for a variety of automated data analysis strategies that use the data created via the FMC method. TFM relies on summing up the amplitude values in a range of time indices in A scans from various transmitter-receiver combinations. Where valid surfaces exist, the amplitudes constructively interfere to image the surface. Where no such surface exists, the amplitudes destructively interfere forming no image. TFM is also described as being equivalent to focused phased array throughout the entire inspection volume.

0194The following abbreviations and acronyms may be used within the context of pipe weld inspection as described below. <tables id="tabl0001" num="0001"><table frame="all"><tgroup cols="2"><colspec colnum="1" colname="col1" colwidth="33mm" /><colspec colnum="2" colname="col2" colwidth="133mm" /><tbody><row><entry>A Scan</entry><entry valign="middle">Time-Amplitude plot for a specific Tx-Rx pair</entry></row><row><entry>B Scan</entry><entry valign="middle">Collection of A scans across the inspection array, in this case transverse to the weld, where amplitude is represented in gray or colour scale</entry></row><row><entry>DP</entry><entry valign="middle">Digitization Point - point along the time axis of the A scan</entry></row><row><entry>feeder</entry><entry valign="middle">Pipe carrying heavy water coolant to or from the individual fuel channels</entry></row><row><entry>FMC</entry><entry valign="middle">Full Matrix Capture</entry></row><row><entry>Inspection Array</entry><entry valign="middle">Multi element transducer used for FMC data collection</entry></row><row><entry>ID</entry><entry valign="bottom">Inside Diameter</entry></row><row><entry>NEOVISION</entry><entry valign="middle">Custom UT data acquisition and analysis application for FMC data sets</entry></row><row valign="bottom"><entry>NPS</entry><entry>Nominal Pipe Size</entry></row><row><entry>OD</entry><entry valign="middle">Outside Diameter</entry></row><row><entry>Tmin</entry><entry valign="middle">Minimum wall thickness</entry></row><row><entry>TFM</entry><entry valign="middle">Total Focus Method</entry></row><row><entry>TSSA</entry><entry valign="middle">Technical Standards &amp; Safety Authority</entry></row><row><entry>UT</entry><entry valign="middle">Ultrasonic Testing</entry></row><row><entry>WPIT</entry><entry valign="middle">Weld Profile Inspection Tool</entry></row></tbody></tgroup></table></tables>

0195FMC inspection superimposes a probe trajectory of a cylindrical geometry on the weld configuration. Depending upon the nature of the joint and the placement of the manipulator over the joint, some distortion of the OD and ID signals can occur. Areas where this may occur are the cheek areas of straight to bend geometries or Grayloc^™ to bend geometries. A possible remedy for this is to attempt re-positioning the manipulator over the joint with the intent of optimizing the signals in the regions where distortion is experienced.

0196There are two types of obstructions that limit inspection of welds. The first type is adjacent structures which block installation of the tool in the desired inspection area. In this event the inspection area may be inaccessible. A second type of obstruction exists where the manipulator is prevented from achieving full rotation due to intrusion into the manipulator path by adjacent structures. The data acquisition software may have an option where, should an obstruction be detected, the manipulator will drive to the scan origin and then acquire data in the opposite direction until either full coverage is obtained or the manipulator again contacts the obstruction. In the latter case the portion of the circumference where no data has been obtained may be noted as access restricted.

0197<figref idref="f0025">Figure 51</figref> shows an example configuration of the various data systems used in collecting and processing the inspection data. The inspection system as a whole is divided into two separate spheres; the data acquisition aspect and the data analysis aspect. The UT (Ultrasound Testing) data generated during data acquisition in the vault 5002 is collected, reviewed and recorded at the acquisition site 5004. After the recording process, the data is exported to the analysis site 5006 where the analysis blade system 5008 resides. The export process will take a few minutes (2-3) due to the file size and the network speed.

0198Analysts may be able to remotely access the analysis site 5006 such that data need not be transferred to the analyst machines 5010. The analysts evaluate the scan data to set processing parameters to appropriate values prior to submitting a job. Once the job has been processed, the analysis operator will be able to review the profiles and associated intensity maps.

0199A simplified diagram of the equipment configuration in the vault 5002 is given in <figref idref="f0026">Figure 52</figref>.

0200A diagram of the equipment configuration at the trailer/remote acquisition site 5004 is given in <figref idref="f0027">Figure 53</figref>. If 12 fibre MT has not been installed between the acquisition site 5004 and the acquisition hardware in the vault 5002, this fibre may be run to the from the acquisition site 5004 to the penetration.

0201The following is a list of equipment used for data acquisition under this example: <tables id="tabl0002" num="0002"><table frame="none"><tgroup cols="3" colsep="0"><colspec colnum="1" colname="col1" colwidth="19mm" /><colspec colnum="2" colname="col2" colwidth="50mm" /><colspec colnum="3" colname="col3" colwidth="98mm" /><thead valign="top"><row><entry align="center"><b>Quantity</b></entry><entry><b>Item</b></entry><entry><b>Description</b></entry></row></thead><tbody><row rowsep="0"><entry>3</entry><entry>Micropulse^™ 5PA</entry><entry>Data Acquisition Instrument - 128 channels with Gigabit Ethernet interface, Firmware version P1 4.5.1.29 and P2 4.5.0.10 or latest qualified release</entry></row><row rowsep="0"><entry>3</entry><entry>Vault Interface Module</entry><entry>Interface unit for manipulators, Micropulse^™ and peripherals</entry></row><row rowsep="0"><entry>2</entry><entry>Acquisition computer</entry><entry>4 GByte memory, Core2 Duo 23 GHz or above, 32 bit Windows XP, 500 GByte HDD, Gigabit Ethernet interface, Dual 19" monitors</entry></row><row rowsep="0"><entry>3</entry><entry>Neovision^™ software</entry><entry>Data acquisition software 1.0 R0 or most recent qualified release</entry></row><row rowsep="0"><entry>3</entry><entry>Manipulator 2"</entry><entry /></row><row rowsep="0"><entry>3</entry><entry>Manipulator 2.5"</entry><entry /></row><row rowsep="0"><entry>3 sets</entry><entry>Mounting Brackets for 0 degree configuration</entry><entry /></row><row rowsep="0"><entry>3 sets</entry><entry>Mounting Brackets for 6 degree configuration</entry><entry>Provided by Kinectrics^™</entry></row></tbody></tgroup></table></tables><tables id="tabl0003" num="0003"><table frame="none"><tgroup cols="3" colsep="0"><colspec colnum="1" colname="col1" colwidth="26mm" /><colspec colnum="2" colname="col2" colwidth="52mm" /><colspec colnum="3" colname="col3" colwidth="89mm" /><thead valign="top"><row><entry align="center"><b>Quantity</b></entry><entry><b>Item</b></entry><entry><b>Description</b></entry></row></thead><tbody><row rowsep="0"><entry>3</entry><entry>Array Transducer</entry><entry>Custom 128 elements, 7.5 MHz, 0.27 mm pitch, with 8 m cable and Hypertronix^™ connection</entry></row><row rowsep="0"><entry>3</entry><entry>Temperature sensor</entry><entry /></row><row rowsep="0"><entry>2</entry><entry>Reference specimen 2" 0 degree configuration</entry><entry /></row><row rowsep="0"><entry>2</entry><entry>Reference specimen 2.5"</entry><entry /></row><row rowsep="0"><entry>2</entry><entry>Reference specimen 2" 6 degree configuration</entry><entry /></row><row rowsep="0"><entry>4</entry><entry>Reference Block stand</entry><entry /></row><row rowsep="0"><entry /><entry>Fibre Optic Cables</entry><entry>12 fibre with MT connections for station side and vault side (lengths required vary depending upon configuration)</entry></row><row rowsep="0"><entry>2</entry><entry>Network Switch</entry><entry>Should support Gigabit standard</entry></row><row rowsep="0"><entry>3</entry><entry>CAT 5e or CAT 6 Ethernet cable</entry><entry>Use of CAT 5 cable is not generally recommended</entry></row><row rowsep="0"><entry>1</entry><entry>Ethernet addressable power bar</entry><entry /></row><row rowsep="0"><entry>As required</entry><entry>Programmable headsets</entry><entry>Communications sets to be set to a dedicated channel</entry></row><row rowsep="0"><entry>1</entry><entry>Couplant Pump</entry><entry>Submersible type with tubing</entry></row><row rowsep="0"><entry>1</entry><entry>Water container</entry><entry>Large capacity - 60 litre suggested</entry></row><row rowsep="0"><entry>1</entry><entry>Flow Control Module</entry><entry /></row><row rowsep="0"><entry>As required</entry><entry>Trailer/Remote vault camera monitor</entry><entry>As per current practice</entry></row><row rowsep="0"><entry>As required</entry><entry>Vault camera control unit</entry><entry>As per current practice</entry></row><row rowsep="0"><entry>As required</entry><entry>Containment side camera controller unit</entry><entry>As per current practice</entry></row><row rowsep="0"><entry>As required</entry><entry>PTZ camera</entry><entry>As per current practice</entry></row><row rowsep="0"><entry>1</entry><entry>Penetration insert</entry><entry>Universal or SG penetration (should support MTP fibre)</entry></row><row rowsep="0"><entry>1</entry><entry>Stick Light</entry><entry /></row><row rowsep="0"><entry>2</entry><entry>GFCI</entry><entry>15 Amp Fusible Link to U-ground</entry></row><row rowsep="0"><entry>As required</entry><entry>Adapter</entry><entry>Twist lock to U ground</entry></row></tbody></tgroup></table></tables>

0202Demineralised water for use as couplant is generally allowed to sit for a minimum of 48 hours, preferably 60 hours prior to use. This is to permit gases dissolved in the water to condense out. Heating the water to 50 C and then allowing it to cool will accelerate this process. If the water is not allowed a period to de-aerate it is very likely bubbles will nucleate on both the feeder and transducer surfaces. The bubbles will continue to grow and may impair the quality of the scan data. Swabbing a thin layer of gel couplant on the probe face may reduce the tendency for bubbles to form, however the couplant will dissolve in time with exposure to water.

0203As an alternate to having the water sit for an extended duration another method is available. Effective deaeration may be performed by spraying the water into a suitable container using a common garden nozzle. The nozzle should be set to a fine pattern. The pressure drop at the nozzle causes the air to separate from the water. Properly done, the water should be ready for immediate use.

0204A catenary cable containing breathing air hose(s), 120 VAC power and communication fibre may be prepared and tied off at the platform with the other end tied off at a suitable point on the vault floor. The fibre may be connected to the penetration, the air hose(s) to a dedicated breathing air header, and the 120 VAC to a GFI and designated receptacle.

0205On the platform the fibre may be terminated at the switch and the power at the addressable power bar. The platform equipment may be mounted on a suitable cart or common instrument rack.

0206Conventional UT data collection instruments may benefit from periodic calibration when employing techniques that utilize signal amplitude as a means to detect, discriminate and size indications. Conventional normal beam techniques as applied for feeder thickness measurement rely upon the accuracy and stability of the instrument time base. The time base of digital instrumentation is dependent upon the system clock, devices that are known for their long term stability. As such there may be no need to perform annual calibration of the instrument since signal amplitude is not a factor. The embodiments discussed in this application similarly do not rely upon signal amplitude for measurement; therefore periodic instrument calibration is not mandatory.

0207Calibration may be advisable when significant differences are noted when connecting a transducer to two separate instruments given the same set-up parameters. The data collection instrument may be calibrated at such times when replacing components, upgrading hardware or repairs to the instrument are necessary.

0208Transducers may be characterized upon receipt from the manufacturer or when first placed into service. Transducers may then be characterized at yearly intervals following being placed into service. The characterization activities may be conducted by qualified maintenance technicians. The characterization tasks are considered to be part of the system maintenance activities and are not part of the calibration activities during the inspection. The manufacturer's report for the transducer may be referenced during the characterization. Correct processing of inspection data depends upon accurate measurements of probe parameters. These parameters include: <ol id="ol0002" compact="compact"><li>(1) Element functionality check The NEOVISION^™ application has a feature which introduces a series of linearly varying delays into each channel that when viewed appear like a saw tooth pattern. This feature is the Probe Element Test and is found under the UT Calibration tab. A quick check of the main diagonal B scan response will highlight any missing elements. The operator may use the cursors to identify the channel number of the missing element(s). The channel numbers of any missing element(s) may be recorded in the maintenance record.</li><li>(2) Transducer delay The transducer can be evaluated on the reference block; however, slight misalignment can introduce errors in the delay measurement. Given the clearance dimensions in the various components of the manipulator, it is difficult to ensure precise location of the probe with respect to the reference block surface. For the purposes of this inspection, the transducer delay is measured separately on a metrology fixture. This value is recorded and is subsequently confirmed to be within an acceptable range once the transducer is mounted in the manipulator and measured on the reference block.</li><li>(3) Element delay In addition to the overall transducer delay, the delay of each element may also be evaluated. Transducer elements have been found to vary in delay significantly across the face of the transducer, in some cases up to 1/2 cycle. Significant variation will introduce error into the subsequent processing of inspection data.</li><li>(4) Frequency spectra test The frequency spectra test acquires the individual responses of the elements from a defined target and then performs the Fourier transform on a segment of the A-scan. The resulting frequency spectrum is compared to the original as provided in the manufacturers report. This value is recorded and is subsequently confirmed to be within an acceptable range once the transducer is mounted in the manipulator 100 and measured on the reference block during maintenance.</li><li>(5) Pulse duration test The pulse duration test is associated with the frequency spectra test. This particular test involves measuring the duration of the reflected waveform at both 6dB and 20dB drop points. This value is recorded and is subsequently confirmed to be within an acceptable range once the transducer is mounted in the manipulator 100 and measured on the reference block during maintenance in the rubber area.</li><li>(6) Beam orientation This test verifies the direction (angle) of the emitted beam with respect to the passive and active planes of the transducer. The beam is generally to be within 0.5 degrees of the vector normal to the transducer surface. Transducers found to be outside of this specification may not be suitable for FMC inspection.</li></ol>

0209This series of tests and steps are not associated with the transducer alone; however, they do affect subsequent data analysis processes. It is intended these checks be performed at a maintenance facility or in a properly equipped rubber area.

(1) Water DAC curve

0210The water DAC (Distance Amplitude Curve) is created by scanning over the water column steps of the reference tube while the DAC curve is inactive. The gain is adjusted such that the interface echo does not exceed 80% FSH (Full Screen Height) on any of the steps. The data file is saved and then analyzed.

0211The maintenance technician using the cursors on the A-scan records the time of arrival and peak amplitude of the interface signal. The steps are measured in series from nearest to furthest. The amount of gain needed to raise the interface signal to 80% FSH is calculated and entered along with the time of arrival of the peak in the DAC table.

0212A series of points is created and a best fit curve is plotted through the data series. Any significant outliers to the curve may be investigated for anomalies (potential air bubbles). See <figref idref="f0028">Figure 54</figref> for an example of a DAC curve. The Water DAC is then activated and a final scan of the reference tube is conducted to ensure all interface signals are at 80% FSH regardless of water column distance. If there is any deviation beyond 1.25 dB (+/- 5% FSH) the curve should generally be recalculated.

0213The Water DAC curve in all other acquisition computers may be updated with the most recent curve. Note that while the Water DAC curve is unique to each transducer, any DAC curve demonstrated to provide the appropriate responses for other transducers may be applied.

(2) Transducer Gain Trim

0214Elements within a transducer can vary in sensitivity +/- 4dB or more from each other. Variation in instrumentation, cables and connectors can also compound the element to element sensitivity variation. The consistency of the TFM result may be improved when these variations are corrected.

0215The following test may be performed at a static position rather than as part of the calibration scan. The acquisition operator may use the calibration feature in the NEOVISION^™ application to adjust and confirm the overall gain and the gain of the individual channels. The calibration feature is found under the UT Calibration tab in Set-up mode. <ol id="ol0003" compact="compact"><li>(a) The Acquisition Operator then prepares the gain trim calibration by setting the calibration gate to encompass the interface signal. The Acquisition Operator ensures there are no bubbles on either the transducer face or on the reference tube as these will invalidate the gain trim adjustment.</li><li>(b) The Acquisition Operator starts the calibration. NEOVISION^™ will automatically adjust the gain and then after a short period (approximately 5-10 seconds) return a message indicating the range of gain needed to adjust the individual elements. In some cases the gain trim utility will report error(s) if elements are outside the nominal gain trim range. This may be the case for missing/dead elements and/or conditions where air bubbles exist either on the transducer/mirror or reference tube surface. The Acquisition Operator may investigate the error messages to correct any conditions beyond the expected result.</li><li>(c) The Acquisition Operator may evaluate the reported range to compare it to previous gain trim adjustment.</li><li>(d) The Acquisition Operator should generally note any differences in individual element response from the previous adjustment. The gain trim adjustment may in some embodiments be within 2.0 dB of the previous evaluation for the same transducer. Differences may arise from the presence of very small air bubbles or transducer aging. The differences may be resolved before proceeding. If the gain trim is within acceptable limits the Acquisition Operator may update the Probe File and apply the gain trim settings for the current inspection. This value is recorded and is subsequently confirmed to be within an acceptable range once the transducer is mounted in the manipulator and measured on the reference block.</li></ol>

(3) Metal DAC curve

0216The metal DAC is created by scanning over the metal path steps while using the Water DAC but no secondary DAC. The gain is adjusted from the baseline gain of the Water DAC such that the first backwall echo does not exceed 40% FSH. This the difference between the Water DAC baseline gain and the gain required to being the first backwall reflection to 40% is the gain offset for the Metal DAC. The file is then saved and analyzed.

0217The maintenance technician using the cursors on the A-scan generally records the time of arrival and peak amplitude of the first backwall signal. The steps are measured in series from thinnest to thickest. The amount of gain needed to raise the backwall signal to 40% FSH is calculated and entered along with the time of arrival of the peak in the DAC table.

0218A series of points is created and a best fit curve is plotted through the data series. Any significant outliers to the curve may be investigated for anomalies. The Metal DAC is then activated and a final scan of the reference tube is conducted to ensure the all interface signals are at 40% FSH regardless of metal thickness. In some embodiments, if there is any deviation beyond 1.25 dB (+/- 5% FSH) the curve should generally be recalculated.

0219The metal DAC curve is saved. The Metal DAC curve in all other acquisition computers may be updated with the most recent curve. Note that while the Metal DAC curve is unique to each transducer, any DAC curve demonstrated to provide the appropriate responses for other transducers may be applied.

(4) Temperature transducer check

0220The temperature transducer is checked for functionality and is not part of the transducer characterization. This can be conducted by measuring the temperatures in an ice water bath as well as boiling water. Alternatively calometric devices may be used to test and calibrate the temperature transducers. Although the temperature transducer is calibrated for a limited range, repeatable demonstration of the range will indicate the transducer is operating normally.

0221After calibrations defined in the previous paragraphs are performed, the resulting updated transducer file may be transferred to the Acquisition computer. The new file may be used when the tool containing the probe (e.g. the manipulator 100) is applied.

0222System calibration is defined as verifying the value for operating variables and confirmation of the continued performance of the FMC data acquisition system as a whole. The process of converting the UT data into the final result is time-consuming, thus analysis is not part of the calibration procedure. Instead it is generally sufficient to demonstrate the relevant UT parameters have not deviated from the recommended values over the course of the inspection.

0223System calibration encompasses checking: <ol id="ol0004" compact="compact"><li>1. Encoder verification</li><li>2. Element Functionality Check</li><li>3. Transducer - manipulator Alignment / Probe Axis Function Check</li><li>4. Metal path measurement verification</li><li>5. Transducer delay check</li><li>6. Water path attenuation check</li><li>7. Temperature sensor verification</li><li>8. Water path measurement check</li></ol>

0224Should the calibration values fall outside of the prescribed or previous applicable values, the calibration may have failed. The Analysis Operator may address possible causes but not vary parameters (e.g. eliminate air bubbles but not change velocity) and re-attempt the calibration scan. If the calibration results are acceptable the inspection can proceed. If the calibration has failed then all inspection files obtained following the most recent valid calibration are generally discarded. The feeder welds in question may be rescanned.

0225The calibration equipment consists of the following items: <ol id="ol0005" compact="compact"><li>1. 2" reference block specimen (see <figref idref="f0029">Figure 55</figref>)</li><li>2. 2" reference block specimen for 6 degree configuration (see <figref idref="f0030">Figure 56</figref>)</li><li>3. 2.5" reference block specimen (see <figref idref="f0031">Figure 57</figref>)</li><li>4. reference block stand</li><li>5. transducer/manipulator combination connected to UT instrument</li></ol>

0226The calibration equipment is configured with the manipulator 100 installed on the reference block as depicted in <figref idref="f0032">Figure 58</figref>, with the reference block specimen depending upon the manipulator size to be calibrated. The arrow 5802 in <figref idref="f0032">Figure 58</figref> indicates the approximate start point of the scan.

0227The platform operator clasps the manipulator on the appropriate size reference tube ensuring the manipulator outer ring is engaged in the stand locator pins. In some embodiments there is a separate reference tube for the conventional and 6 degree manipulator configurations. Calibration of each manipulator configuration is conducted on the appropriate reference tube.

0228The stabilizing arms are then latched over the manipulator. The platform operator indicates to the acquisition operator when the system is ready for the calibration verification. A flow chart for the sequence of calibration steps is given in <figref idref="f0033">Figure 59</figref>.

0229The acquisition operator starts the couplant pump and the scan under the set-up window of NEOVISION^™. The acquisition operator observes signals in both A and B scan windows to confirm the manipulator is full of couplant. Should the manipulator not fill with couplant or be incapable of maintaining water column for the duration of the calibration exercise, the cause may be identified and fixed before continuing with the inspection.

0230Verification of the proper behaviour and values for the encoder is easiest to perform by driving the manipulator from the home position one full revolution to arrive back at the home position. These locations are easily observed by lining up the parting line inner ring halves with the parting line of the stationary ring halves. The value should generally match the circumference to within +/- 0.3 mm. In example embodiments with a 2" manipulator and a 2.5" manipulator, the distance is: 195 mm for the 2" manipulator and 245 mm for the 21/2" manipulator.

0231The following steps may be used by an operator: <ol id="ol0006" compact="compact"><li>(a) Approach the home position from the negative side by driving in the positive direction.</li><li>(b) When the split in the inner rings align with the split in the outer ring, set the position to 0.0 under the 'Motors/Relays'. If this position is overshot, the manipulator should not be driven in the negative direction, as doing so may build manipulator lash into the calibration. If the position is overshot, repeat the approach from the negative side.</li><li>(c) From this zero position, drive the manipulator in the forward direction, 195 mm for the 2" manipulator, and 245 mm for the 2.5" manipulator.</li><li>(d) The manipulator should generally perform a 360 degree rotation and stop so that the split in the inner rotating ring aligns with those of the outer rings. Any offset from this location may represent an error in the calibration or mechanical operation of the manipulator.</li><li>(e) The test is conducted in the opposite direction by returning the manipulator to the home position.</li></ol> The final value in some embodiments should be 0 mm within +/- 0.3 mm.

0232Failure to meet these requirements in either direction may be due to several possible causes: incorrect settings, electrical problems or mechanical problems. If the encoder is beyond these values in either direction the manipulator may be removed from service for further evaluation.

0233Measurement of encoder accuracy can be affected by backlash in the manipulator drive train. The backlash varies from manipulator to manipulator and can increase in time due to wear in the gear train and manipulator bearing surfaces.

0234The NEOVISION^™ application has a feature which introduces a series of linearly varying delays into each channel that when viewed appear like a saw tooth pattern. This test can be conducted in a static position while in set-up mode on the reference sample. Alternatively this test may also be conducted as part of the calibration record by noting the presence/absence of the element responses over a calibration feature on the reference sample.

0235A quick check of the main diagonal B scan response will highlight any missing elements. The operator may use the cursors to identify the channel number of the missing element(s). The channel numbers of any missing element(s) may be recorded in the maintenance record.

0236If the transducer and instrument combination fail to meet the requirement identified in the "Element Functionality" factor set out in Table 3 below, the transducer and/or instrument may be removed from service for repair.

0237The Acquisition Operator then loads the UT and scan settings for the Verification scan. See Tables 1 and 2 below. The Platform Operator may ensure proper cable routing for the calibration scan. When ready, the Acquisition Operator starts the calibration scan. <tables id="tabl0004" num="0004"><table frame="all"><title><b>Table 1</b> - <b>Example UT Parameters for 0 degree calibration</b></title><tgroup cols="3"><colspec colnum="1" colname="col1" colwidth="33mm" /><colspec colnum="2" colname="col2" colwidth="43mm" /><colspec colnum="3" colname="col3" colwidth="86mm" /><thead valign="top"><row><entry><b>Parameter</b></entry><entry><b>Value</b></entry><entry><b>Comments</b></entry></row></thead><tbody><row><entry>Averaging</entry><entry>1</entry><entry /></row><row><entry>Active aperture</entry><entry>128</entry><entry /></row><row><entry>Pulser mode</entry><entry>Echo Trigger</entry><entry /></row><row><entry>Digitizing Frequency</entry><entry>100 MHz</entry><entry>Used for measurement accuracy</entry></row><row><entry>Zero Offset</entry><entry>Specific to probe - nominally</entry><entry /></row><row><entry>Range</entry><entry>1250</entry><entry /></row><row><entry>Gate Start</entry><entry>245 DP</entry><entry /></row><row><entry>Gate Length</entry><entry>2500 DP</entry><entry>For acquiring 17 mm water column and angle</entry></row><row><entry>Gain</entry><entry>∼32-35 dB</entry><entry>Function of probe and system</entry></row><row><entry>Threshold</entry><entry>∼13% FSH</entry><entry>Adjusted as required to preclude triggering on Initial Pulse</entry></row></tbody></tgroup></table></tables><tables id="tabl0005" num="0005"><table frame="all"><title><b>Table 2 - Example UT Parameters for 6 degree calibration</b></title><tgroup cols="3"><colspec colnum="1" colname="col1" colwidth="33mm" /><colspec colnum="2" colname="col2" colwidth="43mm" /><colspec colnum="3" colname="col3" colwidth="86mm" /><thead valign="top"><row><entry><b>Parameter</b></entry><entry><b>Value</b></entry><entry><b>Comments</b></entry></row></thead><tbody><row><entry>Averaging</entry><entry>1</entry><entry /></row><row><entry>Active aperture</entry><entry>128</entry><entry /></row><row><entry>Pulser mode</entry><entry>Echo Trigger</entry><entry /></row><row><entry>Digitizing Frequency</entry><entry>100 MHz</entry><entry>Used for measurement accuracy</entry></row><row><entry>Zero Offset</entry><entry>Specific to probe - nominally</entry><entry /></row><row><entry>Range</entry><entry>1500</entry><entry /></row><row><entry>Gate Start</entry><entry>350 DP</entry><entry /></row><row><entry>Gate Length</entry><entry>3000 DP</entry><entry>For acquiring 17 mm water column and angle</entry></row><row><entry>Gain</entry><entry>∼32-35 dB</entry><entry>Function of probe and system</entry></row><row><entry>Threshold</entry><entry>∼13% FSH</entry><entry>Adjusted as required to preclude triggering on Initial Pulse</entry></row></tbody></tgroup></table></tables>

0238The Acquisition Operator may evaluate the calibration scan once complete as per the items listed in the table below.

0239Transducer - manipulator alignment may be evaluated by scanning across the 5 flat regions of the reference block (such as the example blocks shown in <figref idref="f0029 f0030 f0031">Figures 55-57</figref>). These regions in some embodiments are inclined with respect to the reference block axis. When correctly set, the signal amplitude will rise from the 2 degree surface to the 1 degree surface, peak on the surface parallel to the transducer face and then fall to the 1 degree and again to the 2 degree surfaces. The Acquisition Operator may drive the manipulator from the home position forward around the block by a distance of 145 mm for an example 2" manipulator and 185 mm for an example 2.5" manipulator. The Acquisition Operator may then set the motion increment to 0.1 mm and increment relative to the current position while observing the interface signal. The Operator may continue to increment the position until the signal is observed to peak and then fall. At this point the Operator may set the current position of the manipulator; to 149.1 for the 2" manipulator and to 189.6 mm for the 2.5" manipulator. Driving in the opposite direction may result in manipulator backlash being incorporated into the settings.

0240If the peak is found to occur on other surfaces, the transducer and manipulator may not be aligned and should be adjusted. The manipulator may be removed from service to address this issue.

0241Once the peak is identified and the manipulator circumferential position set as per manipulator size, the Acquisition Operator may drive the manipulator to a point beyond zero: -3 mm is suggested. This may be required to remove lash from the manipulator drive train. The Acquisition Operator may drive the manipulator back to the 0 position (not necessarily home position). The Acquisition Operator is then ready to perform the remaining steps of the calibration.

0242The Acquisition Operator should generally select the WPIT Verification scan, not the WPIT Calibration scan. The Verification scan confirms the function of the DAC curve whilst the Calibration scan builds the DAC curve.

0243The verification scan runs a series of positions around the Reference Block, the metal path series, followed by the water column series, followed by the angle verification and the angled surface.

0244This check may be varied or omitted in embodiments incorporating a reflector or mirror 154 that dynamically changes the scanning angle of the transducer. For example, in some embodiments having a mirror 154, transducer-manipulator alignment is evaluated by scanning across the continuously variable region of the reference block. This region is effectively a cam surface that passes through a range of angles. The control system will drive the mirror to an angle to optimize the interface signal. The angles reported should correspond to the angle of the surface at the various points that are checked.

0245To perform this function manually, the Acquisition Operator may drive the probe axis to peak the interface amplitude at various locations along the span of the continuous surface. The probe axis angle is reported and compared to the expected angle. The difference between the reported angle and the expected angle should not exceed +/- 1 degree in some embodiments.

0246The accuracy of the metal path measurement may be verified by checking various steps in the reference block. In some embodiments, a 2 mm, 8 mm, and 14 mm step are checked. Measurements specific to each example reference block shown in Figured 55-57 are found in Tables A12 to A15 at the end of the Description. <ol id="ol0007" compact="compact"><li>(a) This measurement is conducted on the A scan window of the matrix main diagonal channel selection corresponding to the linear electronic B scan. Any channel from the B scan may be selected for this measurement.</li><li>(b) For this measurement the time axis is set to half path and the velocity is set to carbon steel under the Display Options tab.</li><li>(c) The measurement is made by zooming in on the region of first two backwall echoes in the A scan such that the individual data samples can be distinguished.</li><li>(d) Using both pairs of cursors, the measurement is made by placing one cursor on the first rising zero cross of the first echo and the other cursor on the corresponding zero cross of the second echo. The measurement is the difference between metal path values of these two cursors. This value is directly obtained from the INFO button on the A scan window margin.</li><li>(e) The values for these steps are recorded in the calibration record.</li></ol>

0247If the value(s) fall outside the acceptable range, the cause of the discrepancy should be identified and remedied before the inspection continues. All data obtained prior to identification of the discrepancy may no longer be considered valid.

0248Due to the difficulty in precisely locating the manipulator on the reference sample the transducer delay can vary merely due to positional inaccuracy. Transducer delay is therefore generally performed only as an indication.

0249For this test, the interface signal and second interface multiple should be present in the scan data. This can be achieved on either of the 2 mm or 5 mm water column steps.

0250The transducer delay is evaluated by subtracting the measured water column distance from the design water column distance for that step. The delta will be the change in probe delay from that currently set. The most accurate measurement will be obtained from the first minimum of the respective signals as the subsequent distortion increases the error in the measurement. Note this distance is generally evaluated in digitization points instead of water velocity based half path.

0251The Acquisition Operator may use other steps or channels to confirm a reading that appears to be outside expected limits.

0252The water path attenuation curve may be verified by checking the amplitude of the interface signal.

0253The calibration software (WPIT verification) used in some embodiments has a feature that reports the range of amplitude height variation for a given element on the various water column steps. The amplitude of the interface signal should generally remain within 50 to 99% FSH for any main diagonal channel on all steps of the reference block. Deviation from this amount in the absence of other possible factors (air bubbles) suggests the water path attenuation curve is not accurate.

0254If the reported level is outside the range for an individual element, the Acquisition Operator may perform a manual assessment using another transducer element response. The Acquisition Operator may set the cursor level in the A scan window to 80% FSH. The Acquisition operator may select the desired channel and examine the amplitude response from the range of water column steps. If the responses fall within 50 to 99% FSH, the attenuation check is generally acceptable. If deviation is found to exist on several elements of the transducer, or if the range of deviation is found to be excessive, resolve the issue.

0255Given the robustness of these type of devices it is usually sufficient to verify the temperature reading is both stable and within expected bounds. Nominal bounds anticipated for inspection are generally 15 to 45 degrees Celsius.

0256The Acquisition Operator may save the calibration scan and record the results of the calibration in an FMC Feeder Calibration Record. In preparation for the supplementary scan the manipulator is driven to the home position if not already there.

0257Any anomalies may be resolved by appropriate steps (e.g. component replacement, eliminating air bubbles), and re-calibrated before proceeding with inspection.

0258A supplementary scan may be required to perform a water path measurement check. The scan generally runs from 60 mm to 100 mm on the Reference block in steps of 0.5 mm. This scan may be conducted for both the 0 degree and 6 degree configurations, however only the 0 degree result may be evaluated at this time.

0259In some embodiments, water path measurement is confirmed by checking the measurement deltas between the 2, 5, 8, 11, 14 and 17 mm water column steps in the reference block. This measurement is conducted on the A scan window of the matrix main diagonal channel selection corresponding to the linear electronic B scan. Any channel from the B scan may be selected for this measurement. Measurements specific to each reference block may be kept on record and used as a reference in the calibration process.

0260For this measurement the time axis may be set to half path and the velocity may be set to water under the Display Options tab.

0261The measurement is made by zooming in on the region of interface echoes in the A scan such that the individual data samples can be distinguished. Using both pairs of cursors, the measurement is made by placing one cursor on the first rising zero cross of the interface echo and then indexing to the image of the next step, placing the other cursor on the corresponding zero cross of the interface echo. The measurement is the difference between water path values of these two cursors.

0262This value is directly obtained from the INFO button on the A scan window margin. Subsequent water column step measurements are then made in the same manner.

0263Table 3 below identifies the parameter to be measured and applicable target value to be obtained. <tables id="tabl0006" num="0006"><table frame="all"><title><b>Table 3 - List of example calibration criteria and values</b></title><tgroup cols="4"><colspec colnum="1" colname="col1" colwidth="32mm" /><colspec colnum="2" colname="col2" colwidth="37mm" /><colspec colnum="3" colname="col3" colwidth="45mm" /><colspec colnum="4" colname="col4" colwidth="53mm" /><thead valign="middle"><row><entry><b>Factor</b></entry><entry><b>Value</b></entry><entry><b>Range</b></entry><entry><b>Comments</b></entry></row></thead><tbody><row><entry>Element Functionality</entry><entry>Not applicable</entry><entry>No more than 13 missing, no more than 3 adjacent missing</entry><entry>Reject transducer if greater than 13 or 3 adjacent elements</entry></row><row><entry>Probe delay</entry><entry>32 DP</entry><entry>2 DP</entry><entry>Note - when at 100 MHz sampling rate</entry></row><row><entry>Water column DAC</entry><entry>80 % FSH</entry><entry>99% - 50% FSH</entry><entry /></row><row><entry>Metal path measurement</entry><entry>Block specific, Refer to Tables A12 to A15</entry><entry>+/- 0.03 mm</entry><entry /></row><row><entry>Transducer - manipulator alignment</entry><entry>0</entry><entry>+/- 1</entry><entry /></row><row><entry>Temperature sensor function</entry><entry>Not applicable</entry><entry>15 - 45 °C</entry><entry>Observed to be functioning and rational</entry></row><row><entry>Encoder check</entry><entry>0 to 195 (2") or 245 (2.5")</entry><entry>+/- 0.3 mm</entry><entry>Check for backlash or excess play - replace manipulator if necessary</entry></row><row><entry>Water column measurement</entry><entry>Block specific, Refer to Tables A12 to A15</entry><entry>+/- 0.01 mm</entry><entry>Measure the difference between the steps of the Reference Block</entry></row></tbody></tgroup></table></tables>

0264The FMC Feeder Calibration Record may generally contain the following information: <ul id="ul0004" list-style="bullet"><li>Date and Time of calibration file</li><li>Name of the calibration file</li><li>Transducer serial number, manipulator serial number</li><li>Micropulse serial number, Vault Interface Module serial number</li><li>Reference tube serial number</li><li>Name of Acquisition operator performing calibration</li><li>Name of Platform Operator performing calibration</li><li>Data table recording the results of the calibration</li></ul>

0265A calibration scan of the reference block may be conducted at the following intervals: <ul id="ul0005" list-style="bullet"><li>start of a shift</li><li>4 hour intervals</li><li>change of acquisition operator</li></ul>

0266In addition to the above, calibration scans may be performed whenever a transducer, manipulator or instrument is changed. The calibration check is a confirmation that the inspection system remains in calibration and is fully operational. The calibration check file may be recorded with the inspection data.

0267The ultrasonic equipment settings may be available in a pre-defined configuration (set-up file). The Inspection Supervisor responsible for the initial set-up of the equipment for the inspection campaign may ensure the most recent version of all set-up files are loaded into the NEOVISION^™ software.

0268Nominal values for the settings are given in Table 4, below. <tables id="tabl0007" num="0007"><table frame="all"><title><b>Table 4</b> - <b>Nominal values for Ultrasonic testing</b></title><tgroup cols="4"><colspec colnum="1" colname="col1" colwidth="25mm" /><colspec colnum="2" colname="col2" colwidth="20mm" /><colspec colnum="3" colname="col3" colwidth="66mm" /><colspec colnum="4" colname="col4" colwidth="56mm" /><thead valign="top"><row><entry><b>Parameter</b></entry><entry><b>Essential</b></entry><entry><b>Value</b></entry><entry><b>Comments</b></entry></row></thead><tbody><row><entry>Probe Element Spacing</entry><entry>Yes</entry><entry>0.195 mm</entry><entry>Fixed per probe specification</entry></row><row><entry>Probe Element Pitch</entry><entry>Yes</entry><entry>0.27 mm</entry><entry>Fixed per probe specification</entry></row><row><entry>Probe Element Elevation</entry><entry>Yes</entry><entry>5.0 mm</entry><entry>Fixed per probe specification</entry></row><row><entry>Probe Element Count</entry><entry>Yes</entry><entry>128 elements - no greater than 13 elements non responsive with and no greater than 3 elements adjacent</entry><entry>Fixed per probe specification</entry></row><row><entry>Probe Centre Frequency</entry><entry>Yes</entry><entry>As measured, within +/-10% of nominal</entry><entry>Fixed per probe specification</entry></row><row><entry>Transducer element gain trim</entry><entry>Yes</entry><entry>Elements should be within +/- 2.0 dB of previous</entry><entry>As determined via procedure</entry></row><row><entry>Zero Offset</entry><entry>Yes</entry><entry>Should be within +/- 2DP of known delay, however current designs are ∼ 32 DP @ 100 MHz or 16 DP @ 50 MHz</entry><entry>As measured via probe characteristics</entry></row><row><entry>Pulse width</entry><entry>Yes</entry><entry>40 to 60 ns</entry><entry>Optimized to transducer</entry></row><row><entry>Pulse Voltage</entry><entry>Yes</entry><entry>70V</entry><entry>Exceeding 100V will accelerate transducer aging</entry></row><row><entry>Delay</entry><entry>Yes</entry><entry>∼ 30 DP but not less than 20 DP</entry><entry>Variable - Feeder and weld cap geometry influence the range of this variable</entry></row><row><entry>Range</entry><entry>Yes</entry><entry>∼ 500 to 850 DP</entry><entry>Function of thickness &amp; water column range</entry></row><row><entry>Digitization Rate</entry><entry>Yes</entry><entry>50 MHz</entry><entry>Optimized for data transmission rates</entry></row><row><entry>Digitization Resolution</entry><entry>Yes</entry><entry>12 Bit</entry><entry /></row><row><entry>Averaging</entry><entry>Yes</entry><entry>1</entry><entry /></row><row><entry>Data Collection Mode</entry><entry>Yes</entry><entry>FMC</entry><entry /></row><row><entry>Active Aperture</entry><entry>Yes</entry><entry>128</entry><entry /></row><row><entry>Gain</entry><entry>Yes</entry><entry>Variable - nominal values in a range of 40 to 47 dB</entry><entry>See procedure for details</entry></row><row><entry>Filter</entry><entry>Yes</entry><entry>None</entry><entry>Using the filter influences subsequent processing</entry></row><row><entry>Water DAC</entry><entry>Yes</entry><entry>On - as determined by procedure</entry><entry /></row><row><entry>Metal DAC</entry><entry>No</entry><entry>Off</entry><entry>As per Neovision^™ 1.0 analysis software</entry></row><row><entry>Interface Gate</entry><entry>Yes</entry><entry>On</entry><entry /></row><row><entry>Set interface to Zero</entry><entry>Yes</entry><entry>Off (unchecked)</entry><entry /></row><row><entry>Interface Gate Start</entry><entry>Yes</entry><entry>∼ 225 - 350 DP but not less than 200 DP</entry><entry>Function of scan type, i.e., 0 degree, 6 degree in case of 2"</entry></row><row><entry>Interface Gate Range</entry><entry>Yes</entry><entry>∼ 1000 DP</entry><entry>Function of feeder geometry</entry></row><row><entry>Interface Gate Amplitude</entry><entry>Yes</entry><entry>∼ 11% FSH</entry><entry /></row><row><entry>Pre - trigger acquisition</entry><entry>Yes</entry><entry>On - 30 DP</entry><entry /></row><row><entry>Gate 1</entry><entry>No</entry><entry>Off</entry><entry>Input for data quality indicator</entry></row><row><entry>Gate 1 Start</entry><entry>No</entry><entry /><entry>Input for data quality indicator</entry></row><row><entry>Gate 1 length</entry><entry>No</entry><entry /><entry>Input for data quality indicator</entry></row><row><entry>Gate 1 Amplitude</entry><entry>No</entry><entry>5% FSH</entry><entry>Input for data quality indicator</entry></row><row><entry>Outer Diameter Apertu re</entry><entry>No</entry><entry>18</entry><entry>Input for data quality indicator</entry></row><row><entry>Minimum Percentage</entry><entry>No</entry><entry>50%</entry><entry>Input for data quality indicator</entry></row><row><entry>Preferred Percentage</entry><entry>No</entry><entry>80%</entry><entry>Input for data quality indicator</entry></row><row><entry>Inner Diameter Apertu re</entry><entry>No</entry><entry>32</entry><entry>Input for data quality indicator</entry></row><row><entry>Minimum Percentage</entry><entry>No</entry><entry>50%</entry><entry>Input for data quality indicator</entry></row><row><entry>Preferred Percentage</entry><entry>No</entry><entry>80%</entry><entry>Input for data quality indicator</entry></row></tbody></tgroup></table></tables>

0269Table 4 above identifies the nominal Ultrasonic parameters generally used in this inspection. The Acquisition Operator may confirm the appropriate parameters have been loaded following the calibration run.

0270The temperature of the couplant during acquisition can change by an appreciable margin, with a typical range of 18 to 40 degrees and, in some instances, beyond this range. There is a corresponding change in sound velocity within the couplant. Unlike other inspection methods, velocity change in the couplant can have a significant effect on the results. Therefore the temperature at which the inspection is conducted should generally be monitored. Acceptable temperature range for inspection in some embodiments is 15 to 45 C. If temperatures are encountered outside of this range, efforts to bring the temperature within the specified range may be made.

0271A temperature sensor may be included in the hardware in some embodiments. <figref idref="f0050">Figures 79</figref> and <figref idref="f0051">80</figref> show an example temperature sensor 7902 connected to the manipulator 100 by means of a temperature sensor attachment clip 7904. The Neovision^™ application may record the temperature for each acquisition position.

0272Some embodiments may include a Data Validity Indicator software feature as an aid to the Acquisition Operator in evaluating the data set. The Indicator displays the presence of signals in the interface gate and gate 1 for the OD and ID surfaces respectively.

0273Input parameters are: <ul id="ul0006" list-style="bullet" compact="compact"><li>Interface gate</li><li>Interface gate amplitude</li><li>Gate 1</li><li>Gate 1 amplitude</li><li>OD aperture size, minimum percentage, preferred percentage</li><li>ID aperture size, minimum percentage, preferred percentage</li></ul> Input parameters may be adjusted in set-up mode.

0274A separate inspection sequence may be defined for the 2" and 2.5" manipulators. The settings for the manipulators are given in Tables 5 and 6 below. <tables id="tabl0008" num="0008"><table frame="all"><title><b>Table 5</b> - <b>Inspection sequence parameters - 2.5" manipulator</b></title><tgroup cols="3"><colspec colnum="1" colname="col1" colwidth="34mm" /><colspec colnum="2" colname="col2" colwidth="26mm" /><colspec colnum="3" colname="col3" colwidth="74mm" /><thead valign="top"><row><entry><b>Parameter</b></entry><entry><b>Value</b></entry><entry><b>Comments</b></entry></row></thead><tbody><row><entry>Start</entry><entry>0</entry><entry>Start position may differ from home position</entry></row><row><entry>End</entry><entry>245</entry><entry>Based on pipe diameter</entry></row><row><entry>Scan Interval</entry><entry>0.5</entry><entry>Procedure requirement</entry></row><row><entry>Speed</entry><entry>10</entry><entry /></row><row><entry>Counts/mm</entry><entry>91.80</entry><entry>Calculated value</entry></row><row><entry>Manipulator radius</entry><entry>47.45</entry><entry>Design value</entry></row><row><entry>Reduction ratio</entry><entry>3748.6</entry><entry>Calculated value</entry></row><row><entry>Maximum current</entry><entry>500</entry><entry>Related to maximum surge current - do not exceed</entry></row><row><entry>Stall current</entry><entry>450</entry><entry>Related to maximum continuous current</entry></row><row><entry>Proportional Gain</entry><entry>2000</entry><entry>Determined via testing</entry></row><row><entry>Integral Gain</entry><entry>80</entry><entry>Determined via testing</entry></row><row><entry>Derivative Gain</entry><entry>1000</entry><entry>Determined via testing</entry></row></tbody></tgroup></table></tables><tables id="tabl0009" num="0009"><table frame="all"><title><b>Table 6 - Inspection sequence parameters - 2" manipulator</b></title><tgroup cols="3"><colspec colnum="1" colname="col1" colwidth="33mm" /><colspec colnum="2" colname="col2" colwidth="26mm" /><colspec colnum="3" colname="col3" colwidth="73mm" /><thead valign="top"><row><entry><b>Parameter</b></entry><entry><b>Value</b></entry><entry><b>Comments</b></entry></row></thead><tbody><row><entry>Start</entry><entry>0</entry><entry>Start position may differ from home position</entry></row><row><entry>End</entry><entry>195</entry><entry>Based on pipe diameter</entry></row><row><entry>Scan Interval</entry><entry>0.5</entry><entry>Procedure requirement</entry></row><row><entry>Speed</entry><entry>10</entry><entry /></row><row><entry>Counts/mm</entry><entry>65.64</entry><entry>Calculated value</entry></row><row><entry>Manipulator radius</entry><entry>39.3</entry><entry>Design value</entry></row><row><entry>Reduction ratio</entry><entry>2133.3</entry><entry>Calculated value</entry></row><row><entry>Maximum current</entry><entry>500</entry><entry>Related to maximum surge current - do not exceed</entry></row><row><entry>Stall current</entry><entry>450</entry><entry>Related to maximum continuous current</entry></row><row><entry>Proportional Gain</entry><entry>2000</entry><entry>Determined via testing</entry></row><row><entry>Integral Gain</entry><entry>80</entry><entry>Determined via testing</entry></row><row><entry>Derivative Gain</entry><entry>1000</entry><entry>Determined via testing</entry></row></tbody></tgroup></table></tables>

0275Due to offset between the feeder and the manipulator as well as surface variation of the feeder, the manipulator in some example embodiments can only acquire UT signals in specific orientations. A compromise may be effected wherein the manipulator is configured in a 0 degree mode to acquire intrados and extrados areas, a 6 degree forward configuration to acquire the left cheek and a 6 degree reverse to acquire the right cheek. The UT signals may be degraded when the manipulator is scanning over areas other than the intended region.

0276The start location for the 0 degree scan is in some embodiments defined as the location where the transducer is directly over the right cheek. The manipulator is oriented with the cables facing toward the header (away from the Grayloc^™) and the scan proceeds in the positive direction through the intrados of the bend to the left cheek, then to the extrados of the bend terminating on the right cheek.

0277The start location for the 6 degree forward scan is in some embodiments defined as the location where the transducer is slightly to the right of the fitting intrados. The orientation of the cables is the same as the 0 degree configuration, facing toward the header, away from the Grayloc^™. The scan proceeds from the intrados over the left cheek to just past the extrados centre line. The scan should generally acquire both intrados and extrados and not be less than 100 mm in length.

0278The start location for the 6 degree reverse scan is in some embodiments defined as the location where the transducer is slightly to the left of the fitting intrados. The orientation of the cables is opposite the 0 degree configuration, facing away from the header, towards the Grayloc^™. The scan proceeds from the intrados over the right cheek to just past the extrados centre line. The scan must acquire both intrados and extrados and not be less than 100 mm in length.

0279The Acquisition Operator may identify the scan type, start and length in the Inspection Record.

0280In some embodiments having an adjustable reflector or mirror 154, the dual axis configuration may improve the capabilities of the system. By articulating the UT beam via a mirror 154, the tool may be capable of gathering the required data in one scan for most configurations. In the case of compound joints 2 scans may be required, one that tracks on the first fitting, and a second scan that track on the 2nd fitting.

0281The following steps may be carried out generally in all scan types: <ol id="ol0008" compact="compact"><li>(1) The platform operator may install the manipulator on the feeder, ensuring that the manipulator halves meet and are securely fastened. If direct access to the inspection area is not possible the Platform Operator can install the manipulator at an accessible area of the feeder and reposition the manipulator over the area to be inspected. The drive enclosure on the manipulator should generally be at the top of the feeder for horizontal feeder runs, or with the gear train at the bottom in case of vertical runs. The platform operator may also check the seals on the manipulator to ensure the direction they are flexed in is uniform.</li><li>(2) The platform operator may loop the cable over an adjacent end fitting to relieve strain on the manipulator, ensuring that there is sufficient transducer cable slack available to permit rotation of the manipulator in the intended direction.</li><li>(3) The platform operator and the acquisition operator may re-confirm the identity of the feeder to be inspected after the manipulator has been installed on the feeder.</li><li>(4) Once manipulator placement is confirmed, the acquisition operator may start the couplant pump to fill the manipulator. The platform operator may check the manipulator seals for excessive leakage. Adjustment of the seals may be required to reduce leakage. The platform operator may be required to adjust the flow rate to optimize couplant consumption via the valve on the couplant flow solenoid assembly.</li><li>(5) Coincident with checking the seals, the Platform Operator may also check and confirm that both halves of the manipulator are securely latched. Excessive leakage and poor motion control may be encountered if the halves are not fully engaged.</li></ol>

0282A flowchart for the inspection process is given in <figref idref="f0034">Figure 60</figref>. The Acquisition Operator has the option of using different manipulator scanning modes depending upon the circumstance encountered. The following steps may be carried out generally in all scan types: <ol id="ol0009" compact="compact"><li>(1) The Acquisition Operator may start the static scan once the couplant pump has been started and couplant is confirmed to be leaking from the manipulator. The B scan window should be populated with linear B scan image, the A scan window should be populated with an A scan from the B scan window. If the B scan is present without the A scan, the Acquisition Operator may re-establish the link of the A scan selection to the B scan window. If neither window is populated, the Acquisition Operator and Inspection Supervisor may investigate and remedy the cause before proceeding. Once the A and B scans have been established the Acquisition Operator should make a cursory assessment of the scan quality to identify and correct any obvious deficiencies.</li><li>(2) The Acquisition Operator may jog the transducer to the desired start location, if not using the Home position as the start position. In the alternative the Platform Operator can perform this task using the local control on the VIM. When the transducer has been positioned at the scan start position the Acquisition Operator may set Start position to zero. Again, the Platform Operator may ensure sufficient cable slack and proper cable feed for the intended direction of travel.</li><li>(3) If at any time the Platform Operator sees the possibility of interference or restricted access for the inspection, the Platform Operator may notify the Acquisition Operator. The Acquisition Operator may choose to start the scan at a more advantageous position if clearance is an issue.</li><li>(4) The Acquisition Operator may choose to perform an optional test run. The purpose of the test run is to check the position of weld, signal amplitude over weld cap, manipulator alignment, presence of air bubbles or pockets, or other conditions/artefacts. The suggested value for circumferential resolution is 10 to 15 mm in some embodiments. Any adverse conditions noted in the test run can be corrected before the scan is conducted.</li><li>(5) The Acquisition Operator may initiate the scan from the start position. The platform operator may monitor the progress of the manipulator, paying close attention to proper cable feed. The acquisition operator may monitor the main diagonal B scan and corresponding A scans for data quality. In regions of reduced clearance the transducer may contact adjacent obstructions. It is usually preferable for the manipulator to stall out rather than the transducer pushing the manipulator aside. In the latter case, erroneous data may be recorded.</li><li>(6) The manipulator control may return the manipulator to the start position upon completion of the scan provided a stall condition has not been detected. If a stall condition has been detected, the behaviour of the manipulator may be as per the scan option chosen by the Acquisition Operator.</li><li>(7) The Acquisition Operator may review the completed file for factors affecting scan acceptability. The Acquisition Operator may shut off the couplant pump while reviewing the file. If the scan is not acceptable the Acquisition and Platform Operators may attempt to correct the deficiency and rescan the weld. There are no specified limits to the number of attempts at acquiring acceptable data. The Acquisition Operator may save data that does not meet a set of quality criteria if a number of scan attempts have been made. Notations describing the limitations may be made in the inspection record.</li><li>(8) If the Acquisition Operator deems the file to be acceptable, the Acquisition Operator may save the file and the Inspection Record.</li><li>(9) The Acquisition Operator may drive the manipulator to home position if it not there already.</li><li>(10) The Acquisition Operator may direct the Platform Operator to either detach the manipulator or move it to the next location if on the same feeder.</li><li>(11) The Acquisition and Platform Operators may perform the calibration at specified intervals as noted above. Data obtained subsequent to the most recent valid calibration scan may be discarded and the affected welds may be re-inspected.</li></ol>

0283The extent of the inspection may be as per scope defined above and in <figref idref="f0024">Figures 49-50</figref>. Note the inspection zone defined for these geometries is based on a known position of the weld centreline. The manipulator may in some embodiments be built to track a cylindrical profile and may not necessarily track the weld centreline on complex geometries. On these geometries, the weld centreline may follow a sinusoidal trajectory in the scan data. Users may elect to reposition the manipulator for a second pass if it is deemed necessary to obtain data from other areas of the feeder.

0284With respect to the NEOVISION^™ application, all signals following the interface signal may be recorded.

0285The Acquisition Operator may review the data file with respect to the following criteria to determine scan acceptability. Given the number of factors contributing to signal quality, the Acquisition Operator may monitor the data as it is acquired and spend time after acquisition performing any detailed review of selected areas.

0286The following criteria apply to the area of interest for the scan in question, i.e., extrados/intrados in the case of the 0 degree configuration, left cheek for the 6 degree forward scan, right cheek for the 6 degree reverse scan according to the example embodiments described above. Adverse conditions encountered in areas not relevant to the scan are generally not considered in data quality evaluation. <ol id="ol0010" compact="compact"><li>(1) Skew &amp; Offset: Skew conditions may occur when the transducer is not aligned with the inspection surface in the axial direction. Offset is a similar condition, but may instead occur when the transducer, while parallel to the pipe axis, is a lateral distance sideways from the pipe axis. Both of these conditions may lead to a significant drop in OD signal amplitude as well as splitting of the OD signal. Under these conditions the ID signal may be completely lost. These conditions may be, to a degree, inherent when inspecting tight elbows or joint to joint configurations, but should not generally exceed 40% of the total circumference. Repositioning the manipulator in the axial and/or radial position with respect to the pipe surface may reduce the amount of distortion in the scan.</li><li>(2) Air bubbles: Air bubbles can occur in 3 places: on the inspection surface, on the transducer surface or oscillating in the water column. Air bubbles in general have a negative effect on data quality. Of the three conditions, air bubbles on the pipe surface have the least effect due to their localized nature and their tendency to be small in size. The larger of these bubbles may yield a response that causes acquisition to occur prematurely. Surface air bubbles are excessive when they result in an observable decrease in OD or ID signal amplitude. Air bubbles on the pipe may be reduced by a combination proper cleaning of the pipe surface and use of a surfactant during the cleaning process. Alternately, moving the manipulator axially along the pipe such that the seals scrub the surface can dislodge the majority of the bubbles. Air bubbles located on the transducer/mirror face are more significant insofar as they constantly dampen the pulse from specific elements for the entire duration of the inspection. This type of air bubble is detected by an abnormal persistent drop in OD/ID signals particular to a few elements. This behaviour can be confirmed on the reference block. This type of air bubble may be unacceptable in some embodiments when it causes a drop of 2 dB or more in the OD/ID signals for the length of the scan. Air bubbles on the transducer/mirror face may be mitigated by gently swabbing the surface of the transducer with an approved UT gel couplant. Air bubbles in the water column are evident as erratic triggering of the data acquisition interface. The most common cause of this type of air bubble is air entrained in the couplant supply. Inappropriate triggering of the data acquisition may lead to problems in the data analysis routines and could, under some conditions, give rise to erroneous results. Given the size of the water column, air bubbles in the water column should not generally occupy more than a few individual scan lines of the entire data file in some embodiments. If air bubbles in the water column are found in any appreciable quantity (more than 5% of scan locations in some embodiments), the source of the air bubbles should be identified and eliminated. Rescanning the affected weld may be conducted and noted in the inspection record.</li><li>(3) Air pockets: Air pockets may occur at the upper extremity of the manipulator and may be caused when the leak rate exceeds the couplant supply rate. An air pocket is observed when there is a partial to complete loss of UT signal across the transducer on consecutive scan lines. If the loss of signal is isolated to few channels or is short in duration (~5 frames) the acquisition operator may accept the scan file. Otherwise it may be advisable to address the cause of the excessive leak rate and/or increase the flow to the manipulator and rescan the feeder weld.</li><li>(4) Signal range limits: Factors such as position of the manipulator on the feeder surface, weld cap height and feeder geometry may have a significant influence on the range of water path distance in the scan. Excessive variation in the water path length may have an adverse effect on the scan quality. The most problematic situation may occur when the transducer is too close to the inspection surface (less than 200 DP in some embodiments). In this case the interface signal may merge into the run down from the initial pulse of the transmitting element, potentially leading to a false trigger. Furthermore the interface signal at this point may be adjacent the near zone for the element, in which case the signal amplitude is highly variable. Both of these conditions may introduce errors in the analysis process. The most significant effect, however, may be that the second water column interface will be captured in the data. This feature may be erroneously identified as the backwall since the amplitude may be greater than the true ID signals. In extreme cases such as excessive weld bead height, there may be insufficient clearance between the transducers and the inspection surface resulting in damage to the transducer itself. A potentially less problematic condition is where there is excessive distance between the transducer and the inspection surface. This may occur on the opposite side of the feeder from the condition found in the above paragraph. The effect on data analysis is more subtle as the further away the transducer is from the inspection surface, the more restricted the viewing angle to any given feature. One way to correct for range limit related problems is to reposition the manipulator with the intent of bring the manipulator close to the pipe on the long water column side while moving the manipulator away from the pipe on the short column side. The weld may be rescanned after re-positioning the manipulator. Note that problems introduced by excessive weld cap height might not be addressed by re-positioning the manipulator.</li><li>(5) Electrical Noise: Electrical noise in the scan may be identified as random (in few cases, repeating) signals in the one or more A scan traces. Under the assumption the noise is not correlated, i.e. does not form a repeating signal, the noise may be effectively averaged out in the data analysis. Problems may occur if the amplitude of the noise becomes high relative to the amplitude of the ID signals (within -6 dB of ID in some embodiments). The source of the electrical noise could be from a number of areas related to transducer, cable or connector damage or in some instances may be related to the Micropulse^™ unit itself. If electrical noise is identified in scan data, one means to address the problem is to replace the transducer or the Micropulse^™ instrument and rescan the weld.</li><li>(6) Temperature Range: The temperature range for the data may be 15 to 45 C in some embodiments. Temperature outside this range may cause inaccurate results. Acquisition Operators may rescan such that the temperatures are within these bounds.</li><li>(7) Missing/dead elements: Missing or dead elements can arise from the natural aging of the transducer, damage to the cable or connector, or a fault in the Micropulse^™ instrument. The scan is generally acceptable if no greater than 13 dead/missing elements are identified, provided that of the 13 elements, no more than 3 are immediately adjacent each other. If the fault can be traced to the Micropulse^™ it is generally recommended that the Micropulse^™ be removed for repair and the weld rescanned.</li><li>(8) Interface amplitude: The gain may be set to a level that under nominal conditions interface amplitude is saturated by 6dB. The purpose of this setting is to ensure the weld cap interface signal remains at a sufficiently high level to trigger data acquisition. Excessive gain settings can be problematic for subsequent data analysis. If the weld cap signals are insufficient to trigger data acquisition, the gain may be raised, preferably not more than 6 dB in some embodiments, along with increasing the interface gate amplitude. Caution is advised when lowering gate amplitude since unwanted signals may now cause false acquisition.</li><li>(9) ID signal amplitude: Under nominal conditions, transducer normal to surface, stable water column, inspection surface at intermediate range, the pipe ID signals should attain an amplitude of approximately 10% FSH in some embodiments. If the pipe ID signals deviate by more than +/-6 dB from this value the cause may be investigated and corrected if possible (air bubbles, loose surface dirt).</li><li>(10) Inspection Coverage: Inadequate inspection area cover with respect to the desired coverage may require repositioning of the manipulator and rescanning. In some instances the inspection geometry may limit one-pass coverage necessitating two passes to achieve desired coverage. If this is the case two files may be used for the inspection with appropriate notations made in the inspection record. If obstructions limit the ability to position the manipulator to attain the required coverage, appropriate notations maybe made in the inspection record.</li><li>(11) Start/End Correspondence: The beginning and final frame of a full circumference scan may display the same UT feature at the same location to within +/- 2.0 mm in the axial direction and within +/- 50 DP in the radial direction in some embodiments. If the start and end frames do not correspond to within these limits, the manipulator may have been displaced during scanning. If so, the scan may be repeated.</li></ol>

0287The Acquisition Operator may export the data file to the analysis site. The analysis site may have multiple redundant storage devices to retain the data files. Analysis of the data may occur upon completion of the data acquisition phase.

0288The Inspection Record in some embodiments is a document that logs the parameters for the inspection conducted. There may be a separate inspection record for each feeder inspected. Note there may be multiple welds associated with any given feeder.

0289The following information may be contained in the inspection record: Station, Unit, Date, Time, Face, Feeder, Weld, File Name, Transducer, Manipulator, Micropulse Serial, Vault Interface, ID, Serial Number, Number, Module SN, Acquisition, Platform, Reference Tube ID, Calibration File, Operator, Operator, Resolution, Data Quality, Scan start, Scan length, Scan Type, Restrictions, a sketch of the actual feeder identifying the location and nature of any scanning restrictions, and any other relevant information about the inspection.

0290The data acquisition method for inspection of feeder weld profiles in the present described examples is the application of the Full Matrix Capture (FMC) inspection technique. This method uses multiple independent elements to transmit and receive sound, and thus acquires a very large data array. The array is so large as to effectively preclude manual analysis of the data.

0291The basis for this analysis method is sequential back propagation of time indices in the source A scans. This technique is sometimes generically referred to as the Total Focus Method (TFM). The result of the analysis method is Computed Ultrasonic Tomography (CUT).

0292The following section defines the procedure by which FMC data may be analyzed using the NEOVISION^™ software application in an example embodiment. It provides the procedures for using the NEOVISION^™ software application as an analysis tool in the inspection of a feeder weld area. The data to be analysed is obtained from the inspection of feeder weld joints. This procedure is applicable to FMC data sets obtained from: pipe to pipe, pipe to fitting, fitting to fitting, and fitting to Grayloc^™ welds. See <figref idref="f0024">Figure 49</figref> as an example of a section of a fitting to fitting weld. The intended coverage is 20 mm on either side of the weld centreline. However geometric limitations may constrain the region of access.

0293This example embodiment accommodates the analysis of data files obtained using both 0 and 6 degree offsets as described in the example acquisition procedure above.

0294FMC UT data from the inspection campaign is provided prior to analysis. A multi processor blade server system may be used to support the Neovision^™ analysis software application. This system in some embodiments is linked by an appropriate network to both the acquisition and analysis sites. The network connections should be able to support a data transfer rate 1 GBit/s.

0295Table 7 below lists example equipment used under this procedure in an example embodiment. <tables id="tabl0010" num="0010"><table frame="all"><title><b>Table 7</b> - <b>Hardware and software required for data analysis</b></title><tgroup cols="3"><colspec colnum="1" colname="col1" colwidth="28mm" /><colspec colnum="2" colname="col2" colwidth="26mm" /><colspec colnum="3" colname="col3" colwidth="113mm" /><thead valign="top"><row><entry><b>Quantity</b></entry><entry><b>Item</b></entry><entry><b>Description</b></entry></row></thead><tbody><row><entry>1 per user</entry><entry>Desktop PC (remote)</entry><entry>4 GByte memory, Core2 Duo 2.3 GHz or above, 32 bit Windows XP^™, 500 GByte HDD, 1 GBit Ethernet interface, Dual 19" monitors</entry></row><row><entry>1</entry><entry>Client PC</entry><entry>Configured as required to support the Desktop PC</entry></row><row><entry>1</entry><entry>Storage Device</entry><entry>Network based local device for storage of inspection data and processed results. Suggested size 5 TB minimum</entry></row><row><entry>1</entry><entry>Analysis Blade Server</entry><entry>Analysis system consisting of 2 servers and 28 dual processor quad core blades</entry></row><row><entry>1</entry><entry>Linux^™ OS</entry><entry>SUSE</entry></row><row><entry>1</entry><entry>Matlab^™</entry><entry>Version 2008b</entry></row><row><entry>1</entry><entry>Windows^™ 64 bit XP Server</entry><entry>SP3</entry></row><row><entry>1</entry><entry>NEOVISION^™</entry><entry>Version 1.0 or most recent qualified release</entry></row></tbody></tgroup></table></tables>

0296Verification of the analysis results is achieved by having a second qualified analyst independently perform analysis on the same data set. In this strategy the initial work is called the primary analysis while the following work is termed the secondary analysis.

0297Any qualified analyst may perform either primary or secondary analysis on any data set. Generally, the same analyst does not perform both primary and secondary analysis on the same data set.

0298The value and location of the minimum thickness for both primary and secondary analysis are compared. If the results are within 0.06 mm for the same location in some embodiments the results of the lesser of the two are reported. If the results do not agree within these bounds a third analyst, Lead Analyst, may perform independent evaluation of the data.

0299A Results Coordinator may be responsible to ensure the analysis is conducted per this procedure.

0300Analysis of FMC data sets is a serial chain of processing steps, the results of any given step are built upon the results of the preceding steps. Thus any inaccuracy or loss of data may significantly impair all subsequent steps. The reliability of the results is a function of the input data quality.

0301The following two sections discuss the Analysis workflow and Analysis data flow. The workflow describes the steps taken by the analyst in order to generate the final results. See the flow chart presented in <figref idref="f0035">Figure 61</figref>. The data flow describes the processing performed by the software. The data flow is summarized here to provide background information when determining the processing parameters. The data analysis algorithms are discussed in the previous sections above.

0302The first step in the workflow is initiated in the Data Acquisition phase. The Acquisition Operator uploads a completed scan file that meets the acceptability criteria to the Gateway server. A completed Inspection Record generally also accompanies the data file.

0303The Analyst may then retrieve the UT data file and the corresponding Inspection Record form and initiate an Analysis Record form for the job. The Analyst may then load the UT data file and review the data quality as per acceptance criteria. Any deviations or exceptions from the acceptance criteria should generally already be identified in the Inspection Record for the file. If so, the exceptions may be noted in the Analysis Record. If exceptions to data quality are identified in the UT data that are not reflected in the Inspection Record, the Analyst may initiate a request for rescan of the weld.

0304The Analyst may then review the file recording the UT parameters that influence the analysis settings. The processing parameters may be adjusted as per the guidelines described above. Processing parameters may be recorded in the Analysis record. The job is then submitted for processing.

0305When the processing has been completed the Analyst may retrieve the Results file. The Analyst may initiate a Results Record form. The Analyst may evaluate the results against a number of criteria for overall validity. The data may be re-processed with different series of parameters to address areas that do not meet the validity criteria.

0306When the Analyst has completed evaluation and no further gains in results are obtained, the Analyst may record the minimum thickness information in the Results Record. The Analyst may also produce trend file(s) as specified below.

0307The analysis data flow is depicted in <figref idref="f0036">Figure 62</figref>. In the example embodiment described currently, the process is initiated with the transfer of data from the acquisition site to the Gateway server. The job list will update as the Gateway receives the data, preprocessing of the data begins immediately. Pre-processing includes such steps as data file conversion and digital filtering.

0308The files for the individual scan frames are downloaded to available processors located in the Blade servers. The following steps are performed in sequential order. <ul id="ul0007" list-style="bullet" compact="compact"><li>Up-sampling the data set to equivalent of 100 MHz sampling rate.</li><li>Compensation for DC offset and digital filtering of UT data</li><li>OD Intensity Map formation. The OD Intensity map is formed by evaluating the time indices to each defined point in the inspection volume for each element as per the parameters set. The time indices are summed for each transmitter-receiver pair and the amplitude of the A scan at that time index is summed into the intensity map. The summation process is repeated for all transmitter-receiver combinations that are valid for the specified point in the inspection volume.</li><li>OD Surface Recognition. The OD surface is defined using the parameters provided. The X,Z coordinate pairs are generated by the algorithm. Interpolation may be conducted on low level or missing portions of the Intensity Map. The OD Surface is smoothed as an input to the next processing step. A separate smoothing value is applied for output purposes. Note the OD Surface X,Z pairs form part of the Results output.</li><li>ID Intensity Map formation. The ID Intensity Map is calculated applying the parameters modified. The ID Intensity Map is calculated in a similar manner as the OD Intensity Map however the X,Z coordinates of the OD surface are used to determine the appropriate time indices for the various transmitter-receiver pairs. Options such as OD signal suppression and signal normalization are implemented at this stage prior to forming the ID Intensity Map.</li><li>ID Surface Recognition. The ID surface is defined using the parameters provided. The X,Z coordinate pairs are generated by the algorithm. Interpolation may be conducted on low level or missing portions of the Intensity Map. The X,Z coordinate pairs are smoothed for output purposes. Note the ID Surface X,Z coordinate pairs form part of the Results output.</li></ul>

0309The processed OD and ID coordinates are transferred back to the Gateway server as part of the Results file. This process is repeated for all individual frame files in the Gateway queue. When all available FMC frame files have been processed the job status on the Gateway server is updated to "post-processing". A short interval after the job status is updated to "complete".

0310The Main results window is illustrated in <figref idref="f0037">Figure 63</figref>. The Analyst can use this window to step through the results for each one of the frames. The contents of the Main results window are listed below: <ol id="ol0011" compact="compact"><li>(1) OD Intensity map is plotted in the upper left hand of the Results window</li><li>(2) ID Intensity map is plotted in the middle of the left hand of the Results window</li><li>(3) The combined OD/ID profiles (X,Z coordinates) are plotted in the lower left hand of the Results window.</li><li>(4) Acquisition, Analysis and File information are listed in the lower right hand of the Results window.</li><li>(5) A listing of all TMIN values and locations for each frame, located in the mid right hand side of the Results window.</li><li>(6) The global TMIN is found in the upper right side of the Results view.</li><li>(7) Controls for scrolling through and sorting the Results data are found in the upper right hand.</li><li>(8) Options for generating various 3D views are located in the upper right hand of the Results view.</li></ol>

0311The 3D Generation windows are a group that provide a 3D reconstruction of the inspected volume. The reconstruction is achieved by plotting the OD and ID profiles either in a linear space or in a radial format. These views are exceptionally useful when assessing OD surface or ID surface related features or evaluating thinning trends. The usefulness becomes compromised however when an appreciable number of the scan profiles have significant interpolation or incorrectly identify Intensity map artefacts.

0312The 3D window provides the base 3D development and 3D radial views. An example of this window is given in <figref idref="f0038">Figure 64</figref>. Display options located on the right side of the window permit the Analyst to vary the degree of detail in the view, the content of the view, surface filtering, pan, rotation, zoom and the colour scale applied.

0313The Overview window provides the Analyst with the 3D radial reconstruction of the inspection volume. In addition, this view also provides the Analyst with axial and radial profiles as plotted with respect to the inspection volume. This is useful when evaluating thinning trends on either surface. Again the plots may be panned, rotated, zoomed with different surfaces and detail levels presented. See <figref idref="f0039">Figure 65</figref> as an example of the Overview window.

0314The Pop-up window provides the Analyst with a dedicated window for any of the selected view types. The purpose of this window is to optimize the view for evaluation or reporting purposes. See <figref idref="f0040">Figure 66</figref> for an example of the Pop-up window.

0315The following is a brief synopsis of the Analysis method as per <figref idref="f0036">Figure 62</figref>. Further detail for each step is provided in the paragraph identified for that step.

0316Data analysis can begin with notification that a file has status has changed from "Job Pre-processing" to "Job Ready for Processing". The Analyst may be assigned the role of primary or secondary analysis on a given file provided they do not perform the alternate role on the same file. The Analyst may initiate an Analysis Record form for the file to be analyzed.

Step 1: Calibration verification

0317If the Analyst has not yet verified the calibration scan and record for the current session, the Analyst may retrieve both the calibration data file and record. Both of these records may be checked for metrics as outlined in Table 8 below. If the calibration is not valid, the Analyst may notify the Lead Analyst to resolve the issue with the Work Package Owner and Acquisition FLM. <tables id="tabl0011" num="0011"><table frame="all"><title><b>Table 8</b> - <b>List of Calibration variables and criteria</b></title><tgroup cols="4"><colspec colnum="1" colname="col1" colwidth="37mm" /><colspec colnum="2" colname="col2" colwidth="29mm" /><colspec colnum="3" colname="col3" colwidth="47mm" /><colspec colnum="4" colname="col4" colwidth="54mm" /><thead valign="top"><row><entry><b>Factor</b></entry><entry><b>Value</b></entry><entry><b>Range</b></entry><entry><b>Comments</b></entry></row></thead><tbody><row><entry>Relationship to acquisition setup</entry><entry>Not applicable</entry><entry>Not applicable</entry><entry>Must be same transducer, equipment and personnel as used for data acquisition</entry></row><row><entry>Validity Interval</entry><entry /><entry>Data acquired within 4 hours of Calibration activity</entry><entry /></row><row><entry>Element Functionality</entry><entry>Not applicable</entry><entry>No more than 13 missing, no more than 3 adjacent missing</entry><entry>Reject transducer if outside spec</entry></row><row><entry>Probe delay</entry><entry>32 DP</entry><entry>+/- 2 DP</entry><entry>Note - when at 100 MHz sampling rate</entry></row><row><entry>Water column DAC</entry><entry>80 % FSH</entry><entry>50% to 99% FSH</entry><entry /></row><row><entry>Water column measurement</entry><entry>As per block, see Appendix B</entry><entry>+/- 0.01 mm</entry><entry>Note - when at 100 MHz sampling rate</entry></row><row><entry>Metal path measurement</entry><entry>As per block, see Appendix B</entry><entry>+/- 0.03 mm</entry><entry>Note - when at 100 MHz sampling rate</entry></row><row><entry>Transducer - manipulator alignment</entry><entry>Not applicable</entry><entry>+/- 1 degree</entry><entry>Peak to occur on the zero degree flat</entry></row><row><entry>Temperature sensor function</entry><entry>Not applicable</entry><entry>15-45 C</entry><entry>Observed to be functioning and rational</entry></row><row><entry>Encoder check</entry><entry>0, 195/245</entry><entry>+/- 0.3 mm</entry><entry>Value dependent upon manipulator in use</entry></row></tbody></tgroup></table></tables>

Step 2: Inspection data quality and record verification

0318The Analyst may review the UT data file and corresponding Inspection Record as per the criteria provided. Any exceptions to data quality may be noted in the Inspection Record. If exceptions are not identified in the Inspection Record, the Analyst may request a rescan. <tables id="tabl0012" num="0012"><table frame="all"><title><b>Table 9 - Summary of scan data quality factors</b></title><tgroup cols="4"><colspec colnum="1" colname="col1" colwidth="18mm" /><colspec colnum="2" colname="col2" colwidth="29mm" /><colspec colnum="3" colname="col3" colwidth="27mm" /><colspec colnum="4" colname="col4" colwidth="93mm" /><thead valign="top"><row><entry><b>Number</b></entry><entry><b>Issue</b></entry><entry><b>Condition</b></entry><entry><b>Criteria</b></entry></row></thead><tbody><row><entry>1</entry><entry>Transducer - Feeder Misalignment</entry><entry>Skew</entry><entry>No more than a cumulative total of 40% of the scan missing</entry></row><row><entry>1</entry><entry>Transducer - Feeder Misalignment</entry><entry>Offset</entry><entry>No more than a cumulative total of 40% of the scan missing</entry></row><row><entry>2</entry><entry>Air Bubbles</entry><entry>Pipe surface:</entry><entry>No appreciable reduction in interface or ID signal amplitude</entry></row><row><entry>2</entry><entry>Air Bubbles</entry><entry>Transducer / Mirror surface</entry><entry>Not acceptable</entry></row><row><entry>2</entry><entry>Air Bubbles</entry><entry>Water Column</entry><entry>Not greater than 5 scan lines mistriggered</entry></row><row><entry>3</entry><entry>Air Pockets</entry><entry>Horizontal</entry><entry>Air pocket less than 5 consecutive scan lines</entry></row><row><entry>3</entry><entry>Air Pockets</entry><entry>Vertical</entry><entry>Less than 10 elements and no greater than 40 scan lines</entry></row><row><entry>4</entry><entry>Signal Range Limits</entry><entry>Near</entry><entry>Interface signal not closer than 200 DP</entry></row><row><entry>4</entry><entry>Signal Range Limits</entry><entry>Far</entry><entry>Sufficiently close to trigger Interface Gate, should not exceed 900 DP</entry></row><row><entry>5</entry><entry>Electrical Noise</entry><entry /><entry>Not to exceed 5 % FSH if uncorrelated, to 2% FSH if correlated</entry></row><row><entry>6</entry><entry>Temperature Reading</entry><entry /><entry>Within range of 15 to 45 C and rational</entry></row><row><entry>7</entry><entry>Missing - Dead Elements</entry><entry /><entry>No more than 13 total, and not greater than 3 adjacent</entry></row><row><entry>8</entry><entry>Interface Amplitude</entry><entry /><entry>Approximately 6 dB beyond 100% FSH on straight pipe sections or sufficient to reliably trigger interface gate on weld cap</entry></row><row><entry>9</entry><entry>ID Signal Amplitude</entry><entry>Under pipe material</entry><entry>Observable in Linear Electronic B scan</entry></row><row><entry>10</entry><entry>Inspection Coverage</entry><entry>Axial</entry><entry>Transducer centred over weld - subject to fitting geometry and obstructions</entry></row><row><entry>10</entry><entry>Inspection Coverage</entry><entry>Circumferential</entry><entry>Full circumference scanned - subject to obstructions</entry></row><row><entry>11</entry><entry>Start/End Correspondence</entry><entry>Axial</entry><entry>Start &amp; End UT features in same axial position +/- 1.0 mm</entry></row><row><entry>11</entry><entry>Start/End Correspondence</entry><entry>Radial</entry><entry>Start &amp; End UT features in same radial position +/- 50 DP</entry></row></tbody></tgroup></table></tables>

Step 3: Inspection data characteristics measurement

0319The Analyst may repeat the review of the UT data measuring specific characteristics of the feeder. These characteristics are then used as a basis for modifying the analysis parameters as required. Some specific characteristics that can be measured are set out in Table 10 below. <tables id="tabl0013" num="0013"><table frame="all"><title><b>Table 10</b> - <b>UT Data characteristics related to processing parameters.</b></title><tgroup cols="5"><colspec colnum="1" colname="col1" colwidth="26mm" /><colspec colnum="2" colname="col2" colwidth="20mm" /><colspec colnum="3" colname="col3" colwidth="14mm" colsep="0" /><colspec colnum="4" colname="col4" colwidth="25mm" /><colspec colnum="5" colname="col5" colwidth="83mm" /><thead valign="top"><row><entry><b>Condition</b></entry><entry /><entry namest="col3" nameend="col4" align="center"><b>Processing Parameter(s)</b></entry><entry><b>Comment</b></entry></row></thead><tbody><row rowsep="0"><entry morerows="2" rowsep="1">Water path to interface</entry><entry morerows="2" rowsep="1">Maximum</entry><entry /><entry>X Min, X Max</entry><entry>Increasing water column reduces the effective aperture of the transducer.</entry></row><row rowsep="0"><entry /><entry>Aperture</entry><entry morerows="1" rowsep="1">Excessive Height can introduce ID artefacts if the 2^nd interface is recorded in the UT data. Insufficient Height and X Max can clip ID Intensity map images in areas of thick cross section</entry></row><row><entry /><entry>Height</entry></row><row rowsep="0"><entry morerows="3" rowsep="1">Water path to interface</entry><entry morerows="3" rowsep="1">Minimum</entry><entry /><entry>X Min, X Max</entry><entry morerows="3" rowsep="1">While reducing the water column height allows signals to be observed across a wider aperture, the 2^nd water column interface will limit the maximum thickness that can be imaged. The wider aperture is only effective up to the point where Directivity Cutoff limits the number of channels used.</entry></row><row rowsep="0"><entry /><entry>Aperture</entry></row><row rowsep="0"><entry /><entry>Directivity Cutoff</entry></row><row><entry /><entry>Height</entry></row><row rowsep="0"><entry morerows="2" rowsep="1">Interface signal amplitude</entry><entry morerows="2" rowsep="1" /><entry /><entry>Normalization</entry><entry morerows="1">Low interface signals amplitude may be compensated for by using Normalization - with the potential of introducing artefacts</entry></row><row rowsep="0"><entry /><entry>OD suppression</entry></row><row><entry /><entry /><entry>OD suppression will improve S/N ratio on ID only if minimal/correct interpolation of OD is achieved</entry></row><row><entry>ID signal amplitude in Linear Electronic B scan</entry><entry>Under weld cap</entry><entry /><entry>Normalization</entry><entry>Reducing the threshold levels used for normalization may enhance the ID Intensity map</entry></row><row rowsep="0"><entry morerows="1" rowsep="1">ID signal amplitude in Linear Electronic B scan</entry><entry morerows="1" rowsep="1">Under parent material</entry><entry /><entry>Directivity Cutoff</entry><entry morerows="1" rowsep="1">Changing both the Directivity Cutoff and Max Tx-Rx Difference may aid in the formation of ID Intensity map.</entry></row><row><entry /><entry>Maximum Transmitter- Receiver Difference</entry></row><row rowsep="0"><entry morerows="1" rowsep="1">Weld cap width</entry><entry /><entry /><entry>Maximum Transmitter-Receiver Difference</entry><entry>The aspect ratio (width to height) is an indication of what may be observed under the weld cap, the higher the</entry></row><row><entry /><entry /><entry>Normalization</entry><entry>ratio, the better the reconstruction of the ID</entry></row><row><entry>Weld cap height</entry><entry /><entry namest="col3" nameend="col4" align="left">NA</entry><entry>Only relevant with respect to weld cap width, see above</entry></row><row><entry>Weld cap interface</entry><entry /><entry /><entry>Normalization</entry><entry>Indication of ripple, high modulation may defeat ID reconstruction in one or</entry></row></tbody></tgroup></table></tables>

Step 4: Parameter variation and processing

0320The default parameters can be applied for the first pass in the Analysis process. Specific Inspection data characteristics may warrant changing some processing parameters. The parameters used for each pass may be saved to a unique file associated with the Results file.

Step 5: Results review and verification

0321The Analyst may retrieve and review the results. Depending upon the findings of the Results review, the Analyst may choose to re-process the data one or more times for the entire file or subsets of the file. The Analyst may save each set of processing parameters under a separate file name for each round of processing.

Step 6: Reporting results and documentation

0322The Analyst may record specific results information in the Results record. The Analyst may also create a Trend file. Supplementary outputs may be included as well.

0323In some embodiments, the software will provide the user with additional information to aid the analyst in determining the quality of acquired UT data. In some embodiments, this analysis may take place during Step 5: Results review and verification.

0324In an example embodiment, two indices are provided to the analyst to aid in the analysis. The first index is referred to as the "Theoretical Quality Index". This is based on the geometry of the detected OD surface, and indicates the theoretical best possible quality level achievable for ID UT. The second index is referred to as the "Actual Quality Index". This is based on the ID intensity map from the collected UT data.

0325The Actual Quality Index is generated using the collected UT data. First, the software takes the average of the whole intensity map, ignoring some small areas based on the detected surface. Specifically, the ignored areas are those that extend downwards in the Y axis from the detected surface on the intensity map. By ignoring these areas when calculating the average, the calculation removes all parts of the intensity map that contribute to the detected surfaces. The lower the average intensity of the rest of the intensity map is, the better the "Actual Quality Index" value will be. To arrive at the Actual Quality Index value, this calculated average is compared to one or both of two numbers: <ol id="ol0012" compact="compact"><li>1) the maximum intensity found anywhere in the intensity map (which will likely be found in one of the areas ignored by the averaging), and/or</li><li>2) the maximum possible intensity, based on theoretical UT values that could be achieved based on OD geometry.</li></ol>

0326<figref idref="f0049">Figure 78</figref> shows an example implementation of generating the Actual Quality Index from an intensity map. The detected boundary 7802 is used to identify those areas 7804 extending downwardly therefrom on the Y axis. These areas 7804 are ignored, and only the remaining areas 7806 are used to generate an average intensity value. The lower this average intensity is, the better the Actual Quality Index.

0327Software used to perform the various procedures is described in detail below as part of an example embodiment. The tables A1 to A10 at the end of the Description provide example functions and parameters to be used within the context of this example embodiment. The sections below describing the software used in this example embodiment are intended to be illustrative only, as a possible implementation of the system described in more general terms above.

0328The Feeder Weld Area Thickness Measurement Tool (WPIT) contains a software suite to acquire, analyze and display feeder weld area profiles to determine minimum feeder thicknesses caused by Flow Assisted Corrosion.

0329Neovision^™ is the primary user interface software, responsible for acquisition of the ultrasonic data, control of the probe manipulator, setup of the analysis parameters, communication between itself and the Grid Middleware program in Gateway mode to submit analysis job and retrieve combined result, and initiation of the Results Display program to display the analyzed result. This software can only be executed under the Windows^™ operating system.

0330The Grid Middleware program is a java program and it has 2 modes: Gateway mode and Agent mode.

0331In Gateway mode, it communicates with Neovision^™ and receives an analysis job request from it. It breaks down the Neovision^™ ultrasonic data file into individual FMC data files by invoking the Exporter program. It controls the flow of data to the Grid Middleware Agent, by sending FMC data out and receiving FMC result files back. It then merges all the FMC result files into one combined result file by invoking the Merger program. The combined result file is then sent back to Neovision^™ based on user request.

0332In Agent mode, it receives FMC data files from the Gateway and starts up a set of Analysis programs so that the data can be analyzed. It reports back to the Gateway on its status and the statuses of the Analysis programs. When the FMC result files are created by the Analysis programs, it sends the result files back to the Gateway. It is designed so that each server needs one instance of the Agent and the Agent controls and invokes a set of Analysis programs as defined and limited by configuration settings. The Grid Middleware program is written in Java and it can be run on Windows^™ operating system or Linux^™ operating system. For Gateway mode, it must be executed on Windows operating system since it needs to invoke the Exporter program, which only works on Windows^™ OS. A typical deployment includes a 64-bits Windows^™ server which runs the Grid Middleware in Gateway mode, and a farm of 64-bits Linux^™ Servers which run the Grid Middleware in Agent mode. For a small-scale deployment, the Grid middleware will be executed on the same machine as the Neovision^™ program and it will contain one instance of the Grid Middleware running in Gateway mode and one instance of the Grid Middleware running in Agent mode.

0333The Exporter program exports the Neovision^™ data file into a set of FMC data files for the Analysis program. The FMC data file contains only one set of FMC signals and it is saved as MATLAB^™ data file. These FMC data files will be used as input file for the Analysis program. This program can only be executed under Windows^™ OS. This program is invoked by the Grid Middleware program in Gateway mode.

0334The Analysis program is a MATLAB^™ stand-alone executable program which can be run on Windows^™ operating system or Linux^™ operating system. With a single set of FMC data, it calculates 2D intensity map for the OD and determines the OD boundary definition. When the OD boundary is defined, it then calculates the 2D intensity map for the ID and determines the ID boundary definition. Using the OD and ID boundary definitions, it calculates the minimum thickness (TMin) for this pair of OD and ID boundaries. It then converts the reportable results into a standardized format and saves them. The Analysis program outputs a single result file which contains all the reportable results. It is a highly memory intensive and CPU intensive program. This program is run by the Grid Middleware program in Agent mode.

0335The Merger program takes output analysis result files from the Analysis program and combines all the reportable results to create a single combined result output file. By combining the OD and ID boundaries from each profile, the Merger program can reconstruct a 3D representation of the scan and recreate the OD and ID surfaces. The Merger program is written in MATLAB^™ and it is a MATLAB^™ stand-alone executable program. It is run by the Middleware in Gateway mode.

0336The Result Display program displays the combined result output file so that the user can inspect and interrogate the result. It displays the 2D Intensity Maps, the OD and ID boundaries, the TMin for each profile and the absolute TMin for the whole data file. It also displays the 3D representation of the scan to give user a better understanding of the data. The Result Display program provides a simple user interface so that user can filter out unreasonable TMin values and profiles. Users can create the thickness trend data file and export the filtered data to other CAD (Computer-Aided Design) programs for further investigation such as stress analysis. It runs on the same machine that runs Neovision^™. It is a MATLAB stand-alone executable program and it is invoked from Neovision^™.

0337Here are the data flow steps when analyzing the Neovision data file: <ol id="ol0013" compact="compact"><li>(1) After the Neovision^™ data file is acquired, either the acquisition operator or the analyst can submit the valid data file into Neovision^™ and Neovision^™ will send the data file to Grid Middleware Gateway mode and create a job for this data file.</li><li>(2) The analyst selects the job and inspects the data file. The analyst enters the analysis parameters after reviewing the data file and starts the job. Neovision^™ sends a request to the Grid Middleware and the Grid Middleware Gateway agent will start the analysis by calling the Exporter program</li><li>(3) After the Exporter program converts the Neovision^™ data file into individual FMC data files in MATLAB^™ format, it sends them back to Gateway.</li><li>(4) When the Gateway finds an available Agent machine, it will send FMC .MAT files to the Grid Middleware Agent.</li><li>(5) The Agent creates a node and invokes the Analysis program with the FMC .MAT file as one of the input parameters.</li><li>(6) After the Analysis program completes its calculation, it sends the Result FMC .MAT file back to the Agent.</li><li>(7) The Agent sends the Result FMC .MAT file to the Gateway and informs the Gateway that it has an available node for more analysis computation.</li><li>(8) When Gateway receives all Result FMC .MAT files for the Neovision^™ data file, it will invoke the Merger program and send it all the Result FMC .MAT files.</li><li>(9) The Merger program will extract reportable data and merge individual results into a meaningful order and output a Combined Result file back to the Gateway.</li><li>(10) The Gateway informs Neovision^™ that the result is ready. The user then retrieves the Combined Result file from the Gateway via Neovision^™ and saves a local copy of the Combined Result file.</li><li>(11) After Neovision^™ receives the Combined Result file, it invokes the Results Display program to display the data. The user can review, inspect, modify and export the data in the Result Display program.</li></ol>

0338The grid middleware component can be run in two modes which behave very differently. One is the gateway and the other is the agent. A functioning analysis system will have one instance of Neovision^™, one gateway, and one or more agents. The gateway receives messages from both Neovision^™ and the agent(s) during analysis. Generally, the gateway does not initiate communications, though there are exceptions to this.

0339Analysis of the FMC data is implemented as described above.

0340Where the acquisition software delivers multiple FMC data-sets to the analysis software, the analysis software analyzes each FMC data-set received independently and prescriptively. Thus, a description of the analysis of a single FMC data-set is sufficient to describe the analysis of all FMC data-sets, in that the analysis process for all datasets is identical. Upon completion of the analysis of the FMC data-sets, the output (OD and ID boundary coordinates) are input into the Graphical Output Process.

0341The Profile-Based Result is based on the analysis result. Data is remapped and axes values are changed as needed. Thickness values for each profile are created. The scan position and other profile based information are created in this section.

0342It takes a small amount of time (about 3 seconds at maximum resolution) for each profile to create the Profile-Based Result. However, to create 460 profiles together, even using MATLAB^™ Vectorization technique, it will still take a significant amount of time in the whole process. In order to shorten the total execution time, the Profile-Based Result is created at the end of the analysis program. So the Profile-Based Result is created in parallel without any significant delay on the execution time.

0343Merge Result is a stand-alone MATLAB^™ executable. It merges individual Profile-Based Results in a scan and creates a result file. If the result file already exists, it will merge the specified profiles of Profile-Based Results into the result file.

0344One of the main tasks of Merge Result is to reconstruct the 3D representation of the OD and ID. Since MATLAB^™ uses a 2D matrix to represent any surface and uses the matrix row and column indices to maintain the logical relative positions between data points, Merge Result will merge the OD and ID based on their relative position by first determining the global spacing resolution and the limits. The result from each profile is then remapped, and valid set of data points is inserted into the final OD and ID matrices.

0345The final profiles for the OD and ID are saved into matrices based on their profile ordering. A sorting vector is calculated based on profile's position and their relative positions so that it can be used when displaying the OD and ID in 3D views.

0346All Profile-Based Results are saved into the result file as different structures. All the useful information is stored in the merged result file, so there is no need to keep the output file from the analysis command line program. Invalid, partial (OD only), and empty profiles are merged into the result as null value when the required data does not exist.

0347Merge Result gathers each profile's minimum thickness and saves them into a matrix so that it is easier for the Display Program to access the information.

0348One of the limitations of the Merge Result is that the Profile-Based Result from different profiles must have the same spacing resolution. Merge Result will skip the profile if the spacing settings are not the same.

0349Also, the global spacing resolution and the limits are determined when the result file is being created by iterating the profiles until the first valid profile that has the information is found. Even if newer Profile-Based Results are merged in which have larger ranges, the Merge Result will only merge based on the limits defined previously. If the user wants to re-create a different spacing resolution and different range limit of data, the result file must be removed first so that the global spacing resolutions and limits can be re-determined. Within the FP6 analysis steps, this requires the user to create a new job so that there will not have any existing result file and resubmit the scan for analysis again.

0350This is the main GUI for the Display Result program. All other display windows are created based on this window. The user can invoke the rest of the display windows from this window and the user can also modify TMin and export data from this window.

0351It has an initialize function (initialize) to setup data by reading the result file and cache the data into memory. It loads global setting from a settings file, which defines most of the display properties and the trend output properties. It loads an analysis parameter mapping file so that when the program displays the internal analysis parameters, the parameter names can be the same or similar to the one used in the analysis parameter window in NeoVision^™.

0352This window displays Profile-Based Results by accessing the individual structure for each profile saved in the result file. They are the OD intensity map, ID intensity map, and the OD &amp; ID profiles. It also displays the TMin information and the parameter values for the scan settings, the analysis settings and the profile information.

0353When the user launches the sub-display windows such as popup, overview or 3D, the program stores the handles from those figures (handlePopup, handleOverview, and handle3D) so that it could close all sub-windows when it is terminated. It also uses this list of handles to transmit the Modified TMin information to the sub-windows.

0354A custom data cursor is added into the program. All windows use the same custom function, CursorUpdateText, to display the data cursor.

0355There is no sortable 'uitable' control in MATLAB^™ so the software has to re-create the feature manually for the TMin lists. Based on the button the user clicks, different sorting parameter values are set and the sortUpdateListGUI amd DisplayTMin are called to perform the sorting and display.

0356The Popup window is based on the main window and this is the only window that can have more than 1 instance within the Analysis Display program.

0357The main usage of this window is to allow the maximum screen resolution by displaying a single figure only. The user can view a larger version of any of the figures shown in other display windows.

0358Similar to other sub-windows, the Popup window can be invoked directly and a stand-alone executable program can be created.

0359The Overview window is based on the main window and it is a singleton window.

0360The main usage of this window is to provide a high level overview of the scan. Along with the 3D view of the scan, it provides the cross section view and the axial view of the scan so that user can inspect the result. Refer to the 3D window for more details and description of the 3D View.

0361Similar to the popup window and the 3D window, a median filter can be applied to the display data in the circumferential direction. This should reduce noise and data error so that user can review the data effectively.

0362Similar to the popup window and the 3D window, four levels of data sampling are defined. The default level is medium and it should have the best detail and performance combination.

0363The viewing plane of the cross section view and the axial view are also drawn on the 3D view to enhance the visual representation of the data.

0364Similar to other sub-windows, the Overview window can be invoked directly and a stand-alone executable program can be created.

0365The 3D window is based on the main window and it is a singleton window. The main usage of this window is to provide the 3D representation of the scan. It contains the developmental view and 3D view. The developmental view displays all the profiles of data based on their profile ordering and displays them side by side. The 3D view is the 3D reconstruction of the scan.

0366Both views can be displayed as 'surf or 'mesh' objects in MATLAB^™ drawing. For 'surf', color can be used to represent the thickness of the feeder, or the surface can be drawn monotonically to review the detail of the surface.

0367Similar to the popup window and the Overview window, a median filter can be applied to the display data in the circumferential direction. This should reduce noise and data error so that the user can review the data effectively.

0368Similar to the popup window and the Overview window, four levels of data sampling are defined. The default level is medium and it should have the best detail and performance combination.

0369Similar to other sub-windows, the Overview window can be invoked directly and a stand-alone executable program can be created.

0370The functionality of Modify TMin is handled by the LoadModifiedTMin, SaveModifiedTMin, DisplayTMin, TMinUitableCellEditCallback and menuui_ExportToFigures Callback functions.

0371LoadModifiedTMin reads the MATLAB^™ data file and loads modified TMin information into the program.

0372SaveModifiedTMin writes modified TMin information into data file.

0373DisplayTMin will recreate the TMin lists based on sorting mode and ignore list, and display the TMin lists on the main window.

0374TMinUitableCellEditCallback is invoked when the user modifies the local TMin list by selecting the ignore checkbox to remove or re-add TMins from the list. This function will update the internal definition of the TMin list. Note that when the user removes a TMin from the list, it removes that profile of the scan from the display, not just the TMin definition. The profile removal will affect the display and the export result as those profiles will be omitted.

0375menuui_ExportToFigures_Callback is invoked when the user selects the 'Export Modified TMins To Other Window(s)' option from the menu. The modified TMin list will be sent to all sub-windows and the data on those windows will be updated. It is a two part process where the user will need to refresh individual sub-windows to display the updated TMin information.

0376The Export 3D Point Cloud exports the OD and ID surface as point cloud values.

0377The output filenames will be selected by the user and the OD and ID surfaces will be saved into separate output file with file extension .xyz, which is a defined point cloud file extension for SolidWorks^™.

0378The output file format contain a simple point definition per line and the point is defined as 'x y z', where x, y, z are the coordination of the point.

0379Export Trend Result exports the trend information of the scan to an Excel^™ file.

0380Any existing export file will be deleted.

0381The trend information contains the channel name, the global TMin list, and the TMin Trend.

0382The global TMin list contains the following fields: <ol id="ol0014" compact="compact"><li>1. Profile Number</li><li>2. TMin</li><li>3. Cric Position (mm)</li><li>4. OD Axial Position (X in mm)</li><li>5. OD Depth (Z in mm)</li><li>6. IF Axial Position (X in mm)</li><li>7. ID Depth (z in mm)</li></ol>

0383The TMin Trend contains equally spaced thickness values from the scan. Each line represents a scan profile with the following fields: <ol id="ol0015" compact="compact"><li>1. Profile Number</li><li>2. Cric Position (mm)</li><li>3. Axial Positions (mm)</li></ol>

0384Axial Positions are equally spaced throughout the profile. When the thickness is not available for the particular sampling point, it will be marked as 'NaN'. When the thickness for a particular sampling point is only available using interpolated values, it will be marked as 'IV'. If a particular sampling point is ignored based on the modified TMin selection, it will be marked as 'Ign'. Currently the ignored selection only applies to the profile so the whole row will be marked as 'Ign'. All trend information will be saved into a new tab in the Excel^™ file and the tab name is the channel name.

0385The present invention is defined by the appended claims. <tables id="tabl0014" num="0014"><table frame="all"><title><b>Table A1 - Preprocessing Functions</b></title><tgroup cols="7"><colspec colnum="1" colname="col1" colwidth="19mm" /><colspec colnum="2" colname="col2" colwidth="25mm" /><colspec colnum="3" colname="col3" colwidth="15mm" /><colspec colnum="4" colname="col4" colwidth="42mm" /><colspec colnum="5" colname="col5" colwidth="20mm" /><colspec colnum="6" colname="col6" colwidth="17mm" /><colspec colnum="7" colname="col7" colwidth="31mm" /><thead valign="top"><row><entry><b>Function Name</b></entry><entry><b>Input Data</b></entry><entry><b>Input Data Class</b></entry><entry><b>Input Data Description</b></entry><entry><b>Output Data</b></entry><entry><b>Output Data Class</b></entry><entry><b>Output Description</b></entry></row></thead><tbody><row rowsep="0"><entry morerows="2" rowsep="1">filtfiithd</entry><entry morerows="1">1. Three dimension al Matrix (real).</entry><entry>1. Single</entry><entry>1. Full matrix RF data-set.</entry><entry morerows="2" rowsep="1">Three dimensio nal matrix.</entry><entry morerows="2" rowsep="1">Single</entry><entry morerows="2" rowsep="1">Filtered full matrix RF data-set.</entry></row><row rowsep="0"><entry>2. Filter object</entry><entry morerows="1" rowsep="1">2. Filter object which defines frequency attenuating properties of the digital filter.</entry></row><row><entry>2. Digital software filter.</entry><entry /></row><row><entry>ndimhilbe rt</entry><entry>Three dimension al Matrix.</entry><entry>Single</entry><entry>Full matrix RF data-set.</entry><entry>Three dimensio nal matrix</entry><entry>Single (comple x)</entry><entry>Full matrix analytic time-domain data-set.</entry></row></tbody></tgroup></table></tables><tables id="tabl0015" num="0015"><table frame="all"><title><b>Table A2 - OD Imaging Parameters</b></title><tgroup cols="3"><colspec colnum="1" colname="col1" colwidth="32mm" /><colspec colnum="2" colname="col2" colwidth="17mm" /><colspec colnum="3" colname="col3" colwidth="117mm" /><thead valign="top"><row><entry><b>Variable</b></entry><entry><b>Class</b></entry><entry><b>Description</b></entry></row></thead><tbody><row><entry>blockSize</entry><entry>Single</entry><entry>A Scalar defining the size of the physical aperture used in SFM imaging. This scalar can range from 1 to the size of the array.</entry></row><row rowsep="0"><entry>directivityCutoff</entry><entry>Single</entry><entry>A real number between 0 and 1 (inclusive) used to determine intensity map coordinates considered (with respect to the SFM algorithm) to reflect ultrasonic waves to/from individual array elements.</entry></row><row><entry /><entry /><entry>If the calculated directivity from an element to a coordinate in the intensity map is lower than directivityCutoff, then that element is programmatically not considered to reflect sound to/from the coordinate in question. To calculate the directivity of an element to a specified location, a formula which approximates directivity is used (see Theory Manual [R - 2]).</entry></row><row><entry>xMin</entry><entry>Single</entry><entry>A real number defining the lower limit of x-components of intensity map coordinates (in metres). xMin is less than xMax. X-components of intensity map coordinates will range from xMin to xMax.</entry></row><row><entry>xMax</entry><entry>Single</entry><entry>A real number defining the upper limit of x-components of intensity map coordinates (in metres). xMax is greater than xMin. xMin is less than xMax. X-components of intensity map coordinates will range from xMin to xMax.</entry></row><row rowsep="0"><entry>xSpacing</entry><entry>Single</entry><entry>A vector of length n specifying x-spacings at which to iteratively perform computation of intensities at intensity map coordinates. A total of n iterations of intensity computations are performed (assuming intensity computations do not need to be further broken up to satisfy memory requirements).</entry></row><row><entry /><entry /><entry>Coarser spacings are specified toward the lesser indices of the vector while finer spacings are specified toward the finer indices of the vector. Additionally for every i between 1 and n-1, dx(i+1)/dx(i) is a positive integer.</entry></row><row rowsep="0"><entry>zMin</entry><entry>Single</entry><entry>A real number defining the lower limit of z-components of intensity map coordinates (in metres). zMin is greater or equal than zero.</entry></row><row><entry /><entry /><entry>z-components of intensity map coordinates will range from x to x+height, where x is the maximum of zMin and the distance corresponding to the minimum delay in lookupTable.Offset minus minDistanceBuffer.</entry></row><row rowsep="0"><entry>zSpacing</entry><entry>Single</entry><entry>A vector of length n specifying z-spacings at which to iteratively perform computation of intensities at intensity map coordinates. A total of n iterations of intensity computations are performed (assuming intensity computations do not need to be further broken up to satisfy memory requirements).</entry></row><row><entry /><entry /><entry>Coarser spacings are specified toward the lesser indices of the vector while finer spacings are specified toward the finer indices of the vector. Additionally for every i between 1 and n-1, dz(i+1)/dz(i) is a positive integer.</entry></row><row><entry>height</entry><entry>Single</entry><entry>A real number representing a distance (in metres).</entry></row><row rowsep="0"><entry>minDistanceBuffer</entry><entry>Single</entry><entry>A real number representing a distance (in metres).</entry></row><row><entry /><entry /><entry>z-components of intensity map coordinates will range from x to x+height, where x is the maximum of zMin and the distance corresponding to the minimum delay in lookupTable.Offset minus minDistanceBuffer.</entry></row><row rowsep="0"><entry>zoomPercentage</entry><entry>Single</entry><entry>A vector of length n-1 implicitly specifying coordinates in the intensity map to focus around upon iterative computation of coordinate intensities.</entry></row><row><entry /><entry /><entry>For i between 1 and n-1, zoomPercentage(i) contains a real number between 0 and 1 inclusive. On iteration i+1 of coordinate intensity computation, where Imax is the maximum intensity calculated in the intensity map, the coordinates in which to focus around are those of intensity greater than zoomPercentage(i) * Imax.</entry></row><row><entry>runNorm</entry><entry>Boolean</entry><entry>This determines whether or not normalization will be run on the input data. If this value is false, then the values of normThresh, normWidth and normRunUp are unimportant as they will not be used.</entry></row><row rowsep="0"><entry>norm Thresh</entry><entry>Single</entry><entry>A real number greater than 0 and less than or equal to 100.</entry></row><row rowsep="0"><entry /><entry /><entry>This number specifies the cutoff envelope height of wave packets to be normalized in A-Scan data. Where we let 'y' be the maximum envelope value of all FMC data, wavepacket information with peak height greater or equal to threshold*y/100 will be normalized. All other information in the A-Scans will be set to 0</entry></row><row><entry /><entry /><entry>Recall the envelope of hilbertized data is simply its magnitude.</entry></row><row rowsep="0"><entry>normWidth</entry><entry>Single</entry><entry>A positive scalar greater than 1 and less than or equal to any nonzero length of AScan contained in the FMC data.</entry></row><row><entry /><entry /><entry>This number specifies the distance (in digitization points) between distinct wave packets in an A-Scan. If the distance between envelope peaks greater than the threshold variable exceeds the variable width, then envelope peaks are considered to belong to distinct wave packets. Otherwise, the envelope peaks are considered to belong to the same wave packets.</entry></row><row rowsep="0"><entry>normRunUp</entry><entry>Single</entry><entry>A positive scalar greater than 1 and less than or equal to any nonzero length of AScan contained in the FMC data.</entry></row><row><entry /><entry /><entry>This number specifies the number (in digitization points) of A-Scan values ahead of envelope peaks with values greater than the threshold variable to normalize along with the envelope peaks themselves.</entry></row><row rowsep="0"><entry>gratingThreshold</entry><entry>Single</entry><entry>A real number greater than or equal 0 and less than or equal to 100.</entry></row><row><entry /><entry /><entry>This value is used to help reduce the grating present in the final intensity map. It works by removing low intensity contributions from each aperture. In each aperture, this threshold is applied as follows. The highest intensity location in that threshold is the standard to which 100% is measured. Every other location in the intensity map must be at least gratingThreshold percentage of this value in order to be included in the final summation, otherwise it is zeroed out.</entry></row><row><entry>total Memory</entry><entry>Single</entry><entry>This value is the amount of memory available to this node, in bytes. It is used to fragment the computation of certain values, if necessary. This allows the program to fit within memory constraints. If space allows computation to occur without fragmenting, no fragmenting of computation occurs.</entry></row><row><entry>tempOverride</entry><entry>Boolean</entry><entry>This value determines whether or not the user input temperature value will override the value stored in the slice data file. If this value is false, the value for tempValue is unimportant as it will not be used.</entry></row><row><entry>tempValue</entry><entry>Single</entry><entry>This is the user input value for the temperature. This value will be used for the temperature in calculations if the tempt3verride value is true. Otherwise, this value will be ignored.</entry></row></tbody></tgroup></table></tables><tables id="tabl0016" num="0016"><table frame="all"><title><b>Table A3 - SFM Functions</b></title><tgroup cols="5"><colspec colnum="1" colname="col1" colwidth="28mm" /><colspec colnum="2" colname="col2" colwidth="21mm" /><colspec colnum="3" colname="col3" colwidth="45mm" /><colspec colnum="4" colname="col4" colwidth="21mm" /><colspec colnum="5" colname="col5" colwidth="53mm" /><thead valign="top"><row><entry><b>Function Name</b></entry><entry><b>Input Data / Input Data Class</b></entry><entry><b>Input Data Description</b></entry><entry><b>Output Data / Output Data Class</b></entry><entry><b>Output Description</b></entry></row></thead><tbody><row rowsep="0"><entry morerows="4" rowsep="1">normalize FMClocal 3</entry><entry>1. 2D matrix/single (complex).</entry><entry>1. Filtered and hilbertized FMC data set</entry><entry morerows="4" rowsep="1">1. 2D matrix/single (complex).</entry><entry morerows="4" rowsep="1">1. FMC data set with normalized peaks. All unnormalized data is set to zero.</entry></row><row rowsep="0"><entry>2. scalar/single</entry><entry>2. cutoff threshold used for normalization.</entry></row><row rowsep="0"><entry>3. positive integer/single</entry><entry morerows="1">3. minimum distance (in samples) between wave packet peaks exceeding cutoff threshold such that the wave packets are considered distinct.</entry></row><row rowsep="0"><entry morerows="1" rowsep="1">4. positive integer/single</entry></row><row><entry>4. distance (in samples) to normalize wavepacket information before detected peaks exceeding cutoff thresholds.</entry></row><row rowsep="0"><entry morerows="3" rowsep="1">preparera wdataforo d</entry><entry>1. 2D matrix/ single (complex).</entry><entry>1. Filtered and possibly normalized (if user has set the analysis to run normalization) FMC data set.</entry><entry morerows="3" rowsep="1">1. 2D matrix/ single (complex).</entry><entry morerows="3" rowsep="1">1. FMC data set ordered for correct interpretation in sfmzoom2 algorithm. A-Scans belonging to specific transmitter/receiv er pairs are expected at specific columns in the FMC data set.</entry></row><row rowsep="0"><entry>2. positive integer/single</entry><entry>2. number of probearray elements</entry></row><row rowsep="0"><entry>3. positive integer/single.</entry><entry>3. size of the aperture</entry></row><row><entry>4. 2D matrix/single.</entry><entry>4. 2D matrix whose entries indicate the column in the FMC data set for which the A-Scan is stored, for each transmitter/receiver pair. The transmitter and receivers correspond to the rows and columns of the matrix, respectively.</entry></row><row><entry>pre_comp utearrayc oords</entry><entry>N/A</entry><entry>N/A</entry><entry>1. column matrix/single.</entry><entry>1. Entry x in the column matrix is the x-component position of probe array element x.</entry></row><row rowsep="0"><entry morerows="3" rowsep="1">o_comput einitcoord s3</entry><entry morerows="3" rowsep="1">2D matrix/double</entry><entry morerows="3" rowsep="1">2D matrix whose entries correspond to delay (in samples) between transmission and recording of ultrasonic data for transmitter/receiver pairs. The row and column of each entry corresponds to its transmitter and receiver pair respectively.</entry><entry>1. column matrix/sing le</entry><entry morerows="1">1. x-coordinate for first pass of intensity map calculation.</entry></row><row rowsep="0"><entry>2. column matrix/sing le</entry></row><row rowsep="0"><entry morerows="1" rowsep="1">3. column matrix/sing le</entry><entry>2. z-coordinatefor passof intensity map calculation.</entry></row><row><entry>3. profile number for first pass of intensity map calculation. For each (x,z) pair, the profile number indicates which FMC data set this pair belongs to. Note: this variable array is used for parallel processing of OD intensity maps belonging to different data sets. If only one data set is being processed at a time, every entry in the profile array is set to 1.</entry></row><row rowsep="0"><entry morerows="5" rowsep="1">o_comput ecurrtinec oords4</entry><entry>1. column matrix/single</entry><entry>1. x-coordinates of calculated intensity map coordinates.</entry><entry>1. column matrix/sing le.</entry><entry morerows="3">1. x-coordinate of neighbourhoo d of intensity map coordinates exceeding specified intensity cutoff.</entry></row><row rowsep="0"><entry>2. column matrix/single</entry><entry>2. z-coordinates of calculated intensitymap coordinates.</entry><entry>2. column matrix/single.</entry></row><row rowsep="0"><entry>3. column matrix/single</entry><entry morerows="1">3. profile numbers for calculated intensity map coordinates.</entry><entry morerows="3" rowsep="1">3. column matrix single.</entry></row><row rowsep="0"><entry morerows="2" rowsep="1">4. scalar/double</entry></row><row rowsep="0"><entry morerows="1" rowsep="1">4. loop iteration</entry><entry>2. z-coordinate of neighbourhoo d of intensity map coordinates exceeding specified intensity cutoff.</entry></row><row><entry>3. profile number of neighbourhoo d of intensity map coordinates exceeding specified intensity cutoff.</entry></row><row rowsep="0"><entry morerows="5" rowsep="1">o currcomputeindi cesandvis ibility</entry><entry>1. column matrix/single</entry><entry>1. probe array x- coordinates</entry><entry>1. 2D matrix/sing le</entry><entry morerows="4">1. Entries of matrix yield travel times (in samples) for a given probe array element and intensity map coordinate. The row of the matrix corresponds to the (x,z) coordinate index and the column of the matrix corresponds to the probe array element index.</entry></row><row rowsep="0"><entry>2. column matrix/single</entry><entry morerows="1">2. x-coordinates used for imminent intensity map computation</entry><entry morerows="4" rowsep="1">2. 2D matrix/logi cal</entry></row><row rowsep="0"><entry>3. column matrix/single</entry></row><row rowsep="0"><entry>4. column matrix/single</entry><entry>3. z-coordinates used for imminent intensity map computation.</entry></row><row rowsep="0"><entry /><entry>4. coordinate profiles used for imminent intensity map computation.</entry></row><row><entry /><entry /><entry>2. Entries of matrix yield 'true' if an intensity map coordinate falls within the the specified user defined element directivity, and 'false' otherwise. The row of the matrix corresponds to the (x,z) coordinate index and the column of the matrix corresponds to the probe array element index.</entry></row><row rowsep="0"><entry>sfmzoom 2</entry><entry>1. 3D Matrix/single (complex).</entry><entry>1. Filtered and hilbertized FMC data.</entry><entry>1. 1D Matrix/single</entry><entry morerows="1">1. Computed intensities for coordinates stored in indices identical to those of input item 2.</entry></row><row rowsep="0"><entry /><entry>2. 1D column matrix/single</entry><entry morerows="1">2. Column matrix indicating profile number for each row index in input items 2 and 3.</entry><entry /></row><row rowsep="0"><entry /><entry>3. 2D matrix/single</entry><entry /><entry /></row><row rowsep="0"><entry /><entry>4. 2D matrix/logical</entry><entry morerows="2">3. Entries of matrix yield travel times (in samples) for a given probe array element and intensity map coordinate. The row of the matrix corresponds to the (x,z) coordinate index and the column of the matrix corresponds to the probe array element index.</entry><entry /><entry /></row><row rowsep="0"><entry /><entry>5. 3D matrix/single</entry><entry /><entry /></row><row rowsep="0"><entry /><entry>6. scalar/single</entry><entry /><entry /></row><row rowsep="0"><entry /><entry /><entry>4. Entries of matrix yield 'true' if an intensity map coordinate falls within the the specified user defined element directivity, and 'false' otherwise. The row of the matrix corresponds to the (x,z) coordinate index and the column of the matrix corresponds to the probe array element index.</entry><entry /><entry /></row><row><entry rowsep="0" /><entry rowsep="0" /><entry rowsep="0">5. Entry [i,j,k] of the matrix gives the samples between pulse transmission from element j and receiver reception by element k</entry><entry rowsep="0" /><entry rowsep="0" /></row><row><entry /><entry /><entry>6. Grating threshold that values from each aperture should pass before being added to the main intensity map</entry><entry /><entry /></row><row rowsep="0"><entry morerows="6" rowsep="1">o_comput ecoordina tesabove cutoff</entry><entry>1. 1D column matrix/single.</entry><entry>1. x-coordinates used in intensity map computation</entry><entry>1. 1D column matrix/single.</entry><entry morerows="1">1. newly defined x-coordinate s.</entry></row><row rowsep="0"><entry>2. 1D column matrix/single.</entry><entry>2. z-coordinates used in intensity map computation.</entry><entry>2. 1D column matrix/single.</entry></row><row rowsep="0"><entry>3. 1D column matrix/single.</entry><entry morerows="2">3. coordinate profiles used in imminent intensity map computation.</entry><entry morerows="4" rowsep="1">3. 1D column matrix/sing le.</entry><entry>2. newly defined z-coordinate s.</entry></row><row rowsep="0"><entry>4. 1D column matrix/single.</entry><entry morerows="3" rowsep="1">3. profiles for given newly deifned coordinate s.</entry></row><row rowsep="0"><entry morerows="2" rowsep="1">5. 1D column matrix</entry></row><row rowsep="0"><entry>4. calculated intensities for each coordinate/profile combination.</entry></row><row><entry>5. Entry i in column matrix gives the intensity cutoff for profile i such that new coordinates around coordinates with intensities exceeding the intensity in entry i will have corresponding intensities calculated.</entry></row></tbody></tgroup></table></tables><tables id="tabl0017" num="0017"><table frame="all"><title><b>Table A4</b> - <b>Boundary Recognition Parameters</b></title><tgroup cols="3"><colspec colnum="1" colname="col1" colwidth="36mm" /><colspec colnum="2" colname="col2" colwidth="16mm" /><colspec colnum="3" colname="col3" colwidth="114mm" /><thead valign="top"><row><entry><b>Variable</b></entry><entry><b>Class</b></entry><entry><b>Description</b></entry></row></thead><tbody><row><entry>sigmaMetres</entry><entry>Double</entry><entry>A scalar which defines the size of the convolution window to implement in Canny edge detection.</entry></row><row><entry>thresholdParams</entry><entry>Double</entry><entry>A two element vector where thresholdParams(1) defines the low cutoff threshold and thresholdParams(2) defines the high cutoff threshold in Canny edge detection.</entry></row><row><entry>SEDimsZX</entry><entry>Double</entry><entry>A two element vector (containing only positive integers) defining rectangular structuring element parameters used to dilate the edge detected boundary. SEDimsZX(1) contains the rectangular structuring element height (z-direction) while DEDimsZX(2) contains the rectangular structuring element width(x-direction).</entry></row><row><entry>maxSEDims</entry><entry>Double</entry><entry>A two element vector (containing only positive integers) defining rectangular structuring element parameters used to dilate the pixels of maximum intensity found in vertical slices of the input intensity map. maxSEDims (1) contains the rectangular structuring element height (z-direction) while maxSEDims (2) contains the rectangular structuring element width(x-direction).</entry></row><row><entry>edgeMaxDistance</entry><entry>Double</entry><entry>A scalar defining the maximum allowable distance (in the z-direction) between Canny edge detected points and points of points of maximum intensity (in the z-direction) in the intensity map.</entry></row><row><entry>lengthTrim</entry><entry>Double</entry><entry>A scalar defining the length (in the x-direction) of edge detected connected components ends to trim and replace with translated maximum intensity coordinates found in the same vertical slices as the trimmed edges.</entry></row><row><entry>minAcceptableLength</entry><entry>Double</entry><entry>A scalar defining the minimum acceptable length (in the x-direction) of the calculated boundary. If the length of the calculated boundary is less than the minimum acceptable length, no boundary is returned.</entry></row><row><entry>maxDisjointPieces</entry><entry>Double</entry><entry>A scalar defining the maximum acceptable number of horizontally connected components in the calculated boundary. If the number of horizontally connected components exceeds maxDisjointPieces, no boundary is returned.</entry></row></tbody></tgroup></table></tables><tables id="tabl0018" num="0018"><table frame="all"><title><b>Table A5 - Boundary Recognition Functions</b></title><tgroup cols="6"><colspec colnum="1" colname="col1" colwidth="20mm" /><colspec colnum="2" colname="col2" colwidth="24mm" /><colspec colnum="3" colname="col3" colwidth="68mm" /><colspec colnum="4" colname="col4" colwidth="35mm" /><colspec colnum="5" colname="col5" colwidth="53mm" colsep="0" /><colspec colnum="6" colname="col6" colwidth="43mm" /><thead valign="top"><row><entry><b>Function Name</b></entry><entry><b>Input Data / Input Data Class</b></entry><entry><b>Input Data Description</b></entry><entry><b>Output Data / Output Data Class</b></entry><entry><b>Output Description</b></entry><entry /></row></thead><tbody><row rowsep="0"><entry>makematr ixoutofsp arse</entry><entry>1. column matrix/ single</entry><entry morerows="4" rowsep="1">Items 1-2 are (x-z) coordinates of the intensity map. Item 3 is the stored intensity for each (x,z) coordinate. Thus, items 1-3 have equal length. Items 4-5 are the minimum spacings between x and z coordinates respectively.</entry><entry morerows="4" rowsep="1">Two dimensional matrix/double</entry><entry namest="col5" nameend="col6" morerows="4" rowsep="1" align="left">Matrix containing intensity values. intensity values and their locations are stored in items 1-3. If no intensity can be mapped to a matrix entry, NaN is stored in that entry.</entry></row><row rowsep="0"><entry /><entry>2. column matrix/ single</entry></row><row rowsep="0"><entry /><entry>3. column matrix/ single</entry></row><row rowsep="0"><entry /><entry>4. Scalar variable/ single</entry></row><row><entry /><entry>5. Scalar variable/ single</entry></row><row><entry>fillnan</entry><entry>Two dimensional matrix/double</entry><entry>Matrix containing intensity values of mapped sparse coordinate triples, (x,z,i).</entry><entry>Two dimensional matrix/double</entry><entry>Output matrix is identical to input matrix with sparse NaN entries filled with intensity values of neighbouring entry intensities. Large blocks of NaN entries are filled with the value 0.</entry><entry /></row><row rowsep="0"><entry>leading_e dge</entry><entry>1. Two dimensional matrix/ double</entry><entry>Item 1 contains intensity map in 2D matrix format.</entry><entry>1. Two dimension al matrix/logical</entry><entry>1. Binary image containing detected edges of input item 1.</entry><entry /></row><row rowsep="0"><entry /><entry>2. scalar variable/ double</entry><entry morerows="1" rowsep="1">Items 2 and 3 contain Canny edge detection low and high threshold values respectively.</entry><entry morerows="1" rowsep="1">2. Two dimension al matrix/dou</entry><entry namest="col5" nameend="col6" morerows="1" rowsep="1" align="left">2. smoothed version of input item 1.</entry></row><row><entry /><entry>3. scalar</entry></row><row rowsep="0"><entry /><entry>variable/ double</entry><entry morerows="2" rowsep="1">Items 4 and 5 contain Canny edge detection horizontal and vertical sigma values respectively.</entry><entry>ble</entry><entry /><entry /></row><row rowsep="0"><entry /><entry>4. scalar variable/ double</entry><entry /><entry /><entry /></row><row><entry /><entry>5. scalar variable/ double</entry><entry /><entry /><entry /></row><row rowsep="0"><entry morerows="4" rowsep="1">computep ointsbelo wedgewit hindistanc e</entry><entry>1. row matrix/single</entry><entry morerows="3">Items 1 and 2 contain x and z coordinates (respectively) of maximum intensity pixels of vertical slices of smoothed map output from leading_edge function.</entry><entry>1. row matrix/sing le</entry><entry namest="col5" nameend="col6" morerows="4" rowsep="1" align="left">Items 1 and 2 contain respectively the x and z coordinates of maximum intensity pixels lying directly below (having identical x values but greater z values) the edges outputted from function leading_edge.</entry></row><row rowsep="0"><entry>2. row matrix/single</entry><entry morerows="3" rowsep="1">2. row matrix/dou ble</entry></row><row rowsep="0"><entry>3. row matrix/single</entry></row><row rowsep="0"><entry>4. row</entry></row><row><entry>5. scalar variable/ double</entry><entry>Items 3 and 4 contain x and z coordinates (respectively) of edges output from leading_edge function.</entry></row><row rowsep="0"><entry morerows="3" rowsep="1">computep ointsabov emaxline</entry><entry>1. row matrix/single</entry><entry morerows="1">Items 1 and 2 contain respectively the x and z coordinates of edges output from the leading_edge function.</entry><entry>1. row matrix/sing le</entry><entry namest="col5" nameend="col6" morerows="3" rowsep="1" align="left">Items 1 and 2 are respectively the x and z coordinates of edges outputted from function leading edge above (having identical x values but smaller z values) input items 3 and 4.</entry></row><row rowsep="0"><entry>2. row matrix/single</entry><entry morerows="2" rowsep="1">2. row matrix/dou ble</entry></row><row rowsep="0"><entry>3. row matrix/single</entry><entry morerows="1" rowsep="1">Items 3 and 4 contain respectively x and z coordinates of maximum intensity pixels outputted from function computepointsbelow edgewithindistance.</entry></row><row><entry>4. row matrix/ double</entry></row><row rowsep="0"><entry morerows="5" rowsep="1">computep ointswithi ndistance</entry><entry>1. row matrix/single</entry><entry morerows="2">Items 1 and 2 contain respectively the x and z edge coordinates output from function computepointsabovemaxline</entry><entry>1. row matrix/sing le</entry><entry namest="col5" nameend="col6" morerows="5" rowsep="1" align="left">Items 1 and 2 are respectively the x and z coordinates of edges within the vertical distance of maximum intensity coordinates output from function computepointsbelowed gewithindistance</entry></row><row rowsep="0"><entry>2. row matrix/single</entry><entry morerows="4" rowsep="1">2. row matrix/dou ble</entry></row><row rowsep="0"><entry>3. row matrix/single</entry></row><row rowsep="0"><entry>4. row matrix/doubl e</entry><entry morerows="1">Items 3 and 4 contain respectively the x and z coordinates the maximum intensity pixels output from function computepointsbelow edgewithindistance</entry></row><row rowsep="0"><entry>5. scalar variable/ double</entry></row><row><entry /><entry>Item 5 is the distance (in metres) specifying the maximum distance between input coordinates (1 and 2) and input coordinates (3 and 4)</entry></row><row rowsep="0"><entry>coords2m atrix</entry><entry>1. Two element vector double</entry><entry morerows="2">Item 1 contains the height (first element) and width (second element) of the intensity map matrix output from function makematrixoutofspar se</entry><entry morerows="5" rowsep="1">Two dimensional matrix/logical</entry><entry namest="col5" nameend="col6" morerows="5" rowsep="1" align="left">A binary image containing true in entries mapped from edge locations in input items 4 and 5. False is contained in all other entries.</entry></row><row rowsep="0"><entry /><entry>2. Two element vector/single</entry></row><row rowsep="0"><entry /><entry>3. Two element vector/single</entry></row><row rowsep="0"><entry /><entry>4. row vector/single</entry><entry morerows="1">Item 2 (resp. 3) contains the lower and higher x-axis (resp. z-axis) limits of the intensity map in its first and second elements.</entry></row><row rowsep="0"><entry /><entry>5. row vector/ double</entry></row><row><entry /><entry /><entry>Item 4 (resp. 5) contains the x-axis (resp. z-axis) values of the edge coordinates output from function.</entry></row><row rowsep="0"><entry>imdilate</entry><entry>1. Two dimensional matrix/logical</entry><entry colsep="0" /><entry>1. binary image output from function coords2matrix</entry><entry colsep="1">Two dimensional matrix/logical</entry><entry>Input two dimensional binary matrix post dilation operation.</entry></row><row><entry /><entry>2. rectangular structuring element/ structuring element.</entry><entry colsep="0" /><entry>2. rectangular structuring element used in morphological dilation operation of input item 1.</entry><entry colsep="1" /><entry /></row><row><entry>bwlabeln</entry><entry>Two dimensional matrix/logical</entry><entry colsep="0" /><entry>Binary image output from function imdilate</entry><entry colsep="1">Two dimensional matrix/double</entry><entry>Input matrix with false entries set to 0 and individual connected components enumerated by the positive integers 1,2,3,...</entry></row><row rowsep="0"><entry>bwmorph</entry><entry>1. Two dimensional matrix/logical</entry><entry namest="col3" nameend="col4" align="left">1. Binary image output from function imdilate.</entry><entry colsep="1">Two dimensional matrix/logical</entry><entry>input binary image thinned.</entry></row><row rowsep="0"><entry /><entry>2. row vector/string</entry><entry namest="col3" nameend="col4" morerows="1" align="left">2. string defining morphology operation. String is input as constant, 'thin'.</entry><entry colsep="1" /><entry /></row><row rowsep="0"><entry /><entry morerows="1" rowsep="1">3. scalar variable/ double</entry><entry colsep="1" /><entry /></row><row><entry /><entry namest="col3" nameend="col4" align="left">3. scalar defining the number of thinning operations to perform on input binary image (input item 1). This scalar is set to Inf such that thinning operations will continue to be successively applied until they are no longer effective.</entry><entry colsep="1" /><entry /></row><row><entry>trimcomp onsameli ne</entry><entry>Two dimensional matrix/logical</entry><entry>Binary image output from bwmorph</entry><entry>Two dimensional matrix/logical</entry><entry>Connected components of input binary image trimmed such that smaller connected components with pixels on the same vertical slices as larger connected components have these pixels set to false.</entry><entry /></row><row><entry>removeall connectio nsbetwee nlines</entry><entry>Two dimensional matrix/logical</entry><entry>Binary image output from trimcomponsameline</entry><entry>Two dimensional matrix/logical</entry><entry>Input binary image with all junctions between lines removed.</entry><entry /></row><row><entry>removebo ttomcomp onent</entry><entry>Two dimensional matrix/logical</entry><entry>Binary image output from removeallconnection sbetweenlines</entry><entry>Two dimensional matrix/logical</entry><entry>Connected components of input binary image trimmed such that only their top (lowest corresponding z value) pixels with respect to any column in the matrix are preserved.</entry><entry /></row><row><entry>removeco mponsam evertline</entry><entry>Two dimensional matrix/logical</entry><entry>Binary image output from removebottomcompo nent</entry><entry>Two dimensional matrix/logical</entry><entry>Connected components of input binary image totally removed such that smaller connected components with pixels on the same vertical slices as larger connected components are set to false.</entry><entry /></row><row><entry>replaceci gartips</entry><entry>1. Two dimensional</entry><entry>1. Binary image output from function removecomponsame</entry><entry>Two dimensional</entry><entry>Binary image containing input item 1 with ends of</entry><entry /></row><row rowsep="0"><entry /><entry>matrix/logical</entry><entry morerows="1">vertline containing edges output from function bwlabeln belonging to a particular label</entry><entry>matrix/logical</entry><entry namest="col5" nameend="col6" morerows="5" rowsep="1" align="left">horizontally connected components replaced by vertically translated points in input item 2.</entry></row><row rowsep="0"><entry /><entry>2. Two dimensional matrix/logical</entry><entry /></row><row rowsep="0"><entry /><entry>3. scalar/ double</entry><entry morerows="1">2. Binary image with maximum intensity coordinates output from function computepointsbelow edgewithindistance</entry><entry /></row><row rowsep="0"><entry /><entry morerows="2" rowsep="1">4. rectangular structuring element/stru cturing element.</entry><entry /></row><row rowsep="0"><entry /><entry>3. Calculated number of pixels to trim off the horizontal ends of input item 1.</entry><entry /></row><row><entry /><entry>4. rectangular structuring element used in morphological dilation operation of input item 2.</entry><entry /></row><row rowsep="0"><entry morerows="2" rowsep="1">matrix2co ordswithn an</entry><entry>1. Two dimensional matrix/logical</entry><entry morerows="1">Item 1 contains a binary image outputted from function replacecigartps</entry><entry>1. column matrix/single</entry><entry namest="col5" nameend="col6" morerows="2" rowsep="1" align="left">Items 1 and 2 contain respectively the x-axis values and z-axis values of the calculated edge binary image output from function replacecigartips</entry></row><row rowsep="0"><entry>2. Two element vector/single</entry><entry morerows="1" rowsep="1">2. column matrix/single</entry></row><row><entry>3. Two element vector/single</entry><entry>Item 2 (resp. 3) contains the lower and higher x-axis (resp. z-axis) limits of the intensity map in its first and second elements.</entry></row><row rowsep="0"><entry morerows="1" rowsep="1">Matrix2co ordsinterp only</entry><entry>1. Two dimensional matrix/logical</entry><entry morerows="1" rowsep="1">Item 1 contains a binary image outputted from function replacecigartps</entry><entry>1. column matrix/single</entry><entry namest="col5" nameend="col6" morerows="1" rowsep="1" align="left">Items 1 and 2 contain respectively the x-axis values and z-axis values of linearly interpolated coordinates of empty</entry></row><row><entry>2. Two element</entry><entry>2. column matrix/single</entry></row><row rowsep="0"><entry /><entry>vector/single</entry><entry morerows="1" rowsep="1">Item 2 (resp. 3) contains the lower and higher x-axis (resp. z-axis) limits of the intensity map in its first and second elements.</entry><entry colsep="0" /><entry colsep="1" /><entry morerows="1" rowsep="1">regions of the calculated edge binary image output from function replacecigartips</entry></row><row><entry /><entry>3. Two element vector/single</entry><entry colsep="0" /><entry colsep="1" /></row></tbody></tgroup></table></tables><tables id="tabl0019" num="0019"><table frame="all"><title><b>Table A6 - Boundary Definition Parameters</b></title><tgroup cols="3"><colspec colnum="1" colname="col1" colwidth="24mm" /><colspec colnum="2" colname="col2" colwidth="15mm" /><colspec colnum="3" colname="col3" colwidth="128mm" /><thead valign="top"><row><entry><b>Variable</b></entry><entry><b>Class</b></entry><entry><b>Description</b></entry></row></thead><tbody><row><entry>width</entry><entry>Single</entry><entry>This determines the width of the median filter used in the MedianFilter function. Width determines the number of data points used when choosing a median to replace a given data point with.</entry></row><row><entry>smoothWidth</entry><entry>Single</entry><entry>This determines the width of the Savitzky-Golay filter applied in the PerformSmoothing function.</entry></row><row><entry>smoothOrder</entry><entry>Single</entry><entry>This determines the order of the polynomial used to approximate the surface in the Savitzky-Golay filter applied in the PerformSmoothing function.</entry></row></tbody></tgroup></table></tables><tables id="tabl0020" num="0020"><table frame="all"><title><b>Table A7 - Boundary Definition Functions</b></title><tgroup cols="5"><colspec colnum="1" colname="col1" colwidth="23mm" /><colspec colnum="2" colname="col2" colwidth="27mm" /><colspec colnum="3" colname="col3" colwidth="52mm" /><colspec colnum="4" colname="col4" colwidth="31mm" /><colspec colnum="5" colname="col5" colwidth="34mm" /><thead valign="top"><row><entry><b>Function Name</b></entry><entry><b>Input Data / Input Data Class</b></entry><entry><b>Input Data Description</b></entry><entry><b>Output Data / Output Data Class</b></entry><entry><b>Output Description</b></entry></row></thead><tbody><row rowsep="0"><entry morerows="3" rowsep="1">MedianFilter</entry><entry>1. 1D column matrix/single</entry><entry>1. x-components of surface to be filtered</entry><entry>1. 1D column matrix/single</entry><entry>1. x-components of filtered surface</entry></row><row rowsep="0"><entry>2. 1D column matrix/single</entry><entry morerows="1">2. y/z-components of surface to be filtered</entry><entry morerows="2" rowsep="1">2. 1D column matrix/single</entry><entry morerows="2" rowsep="1">2. y/z-components of filtered surface</entry></row><row rowsep="0"><entry morerows="1" rowsep="1">3. positive integer/single</entry></row><row><entry>3. Width of median filter to be used</entry></row><row rowsep="0"><entry morerows="4" rowsep="1">Perform Smoothing</entry><entry>1. 1D column matrix/single</entry><entry>1. x-components of surface to be filtered</entry><entry>1. 1D column matrix/single</entry><entry>1. x-components of smoothed surface</entry></row><row rowsep="0"><entry>2. 1D column matrix/single</entry><entry morerows="1">2. y/z-components of surface to be filtered</entry><entry morerows="3" rowsep="1">2. 1D column matrix/single</entry><entry morerows="3" rowsep="1">2. y/z-components of smoothed surface</entry></row><row rowsep="0"><entry>3. positive integer/single</entry></row><row rowsep="0"><entry>4. positive integer/single</entry><entry>3. Width of smoothing window to be used</entry></row><row><entry /><entry>4. Order of polynomial used in smoothing approximation</entry></row></tbody></tgroup></table></tables><tables id="tabl0021" num="0021"><table frame="all"><title>Table A8 - IFM OD boundary Preparation Parameters</title><tgroup cols="3"><colspec colnum="1" colname="col1" colwidth="40mm" /><colspec colnum="2" colname="col2" colwidth="15mm" /><colspec colnum="3" colname="col3" colwidth="112mm" /><thead valign="top"><row><entry><b>Variable</b></entry><entry><b>Class</b></entry><entry><b>Description</b></entry></row></thead><tbody><row><entry>width</entry><entry>Single</entry><entry>This determines the width of the median filter used in the MedianFilter function. Width determines the number of data points used when choosing a median to replace a given data point with.</entry></row><row><entry>smoothWidth</entry><entry>Single</entry><entry>This determines the width of the Savitzky-Golay filter applied in the PerformSmoothing function.</entry></row><row><entry>smoothOrder</entry><entry>Single</entry><entry>This determines the order of the polynomial used to approximate the surface in the Savitzky-Golay filter applied in the PerformSmoothing function.</entry></row><row><entry>surfaceSamplinglnterval</entry><entry>Single</entry><entry>This number is a positive integer. The number 1 would indicate that every point on the given OD surface should be used. The number 2 would indicate to use every other point, the number 3 every third, et cetera.</entry></row></tbody></tgroup></table></tables><tables id="tabl0022" num="0022"><table frame="all"><title><b>Table A9 - IFM OD boundary Preparation Functions</b></title><tgroup cols="5"><colspec colnum="1" colname="col1" colwidth="23mm" /><colspec colnum="2" colname="col2" colwidth="27mm" /><colspec colnum="3" colname="col3" colwidth="52mm" /><colspec colnum="4" colname="col4" colwidth="31mm" /><colspec colnum="5" colname="col5" colwidth="34mm" /><thead valign="top"><row><entry><b>Function Name</b></entry><entry><b>Input Data / Input Data Class</b></entry><entry><b>Input Data Description</b></entry><entry><b>Output Data / Output Data Class</b></entry><entry><b>Output Description</b></entry></row></thead><tbody><row rowsep="0"><entry morerows="2" rowsep="1">MedianFilter</entry><entry>1. 1D column matrix/single</entry><entry>1. x-components of surface to be filtered</entry><entry>1. 1D column matrix/single</entry><entry>1. x-components of filtered surface</entry></row><row rowsep="0"><entry>2. 1D column matrix/single</entry><entry>2. y/z-components of surface to be filtered</entry><entry morerows="1" rowsep="1">2. 1D column matrix/single</entry><entry morerows="1" rowsep="1">2. y/z-components of filtered surface</entry></row><row><entry>3. positive integer/single</entry><entry>3. Width of median filter to be used</entry></row><row rowsep="0"><entry morerows="4" rowsep="1">Perform Smoothing</entry><entry>1. 1D column matrix/single</entry><entry>1. x-components of surface to be filtered</entry><entry>1. 1D column matrix/single</entry><entry morerows="1">1. x-components of smoothed surface</entry></row><row rowsep="0"><entry>2. 1D column matrix/single</entry><entry morerows="1">2. y/z-components of surface to be filtered</entry><entry morerows="3" rowsep="1">2. 1D column matrix/single</entry></row><row rowsep="0"><entry>3. positive integer/single</entry><entry morerows="2" rowsep="1">2. y/z-components of smoothed surface</entry></row><row rowsep="0"><entry morerows="1" rowsep="1">4. positive integer/single</entry><entry>3. Width of smoothing window to be used</entry></row><row><entry>4. Order of polynomial used in smoothing approximation</entry></row></tbody></tgroup></table></tables><tables id="tabl0023" num="0023"><table frame="all"><title><b>Table A10- ID Imaging Parameters</b></title><tgroup cols="3"><colspec colnum="1" colname="col1" colwidth="33mm" /><colspec colnum="2" colname="col2" colwidth="15mm" /><colspec colnum="3" colname="col3" colwidth="119mm" /><thead valign="top"><row><entry><b>Variable</b></entry><entry><b>Class</b></entry><entry><b>Description</b></entry></row></thead><tbody><row rowsep="0"><entry>directivityCutoff</entry><entry>Single</entry><entry>A real number between 0 and 1 (inclusive) used to determine OD boundary coordinates considered (with respect to the SFM algorithm) to refract ultrasonic waves from individual array elements.</entry></row><row><entry /><entry /><entry>If the calculated directivity from an element to an OD boundary coordinate is lower than directivityCutoff, then that OD boundary coordinate is programmatically not considered to refract sound from the element. To calculate the directivity of an element to a specified location, a formula which approximates directivity is used (see Theory Manual [R - 1]).</entry></row><row><entry>speedOfMedium2</entry><entry>Single</entry><entry>A real number defining the speed at which sound travels in the second inspection material found between the OD and ID surfaces.</entry></row><row rowsep="0"><entry>xMin</entry><entry>Single</entry><entry>A real number defining the absolute lower limit of x-components of intensity map coordinates (in millimetres). xMin is less than xMax.</entry></row><row><entry /><entry /><entry>Where y is the minimum x-component of the OD boundary, the actual lower limit of the x-components of intensity map coordinates is the maximum of the following two scalars, xMin and y-horzBuffer.</entry></row><row rowsep="0"><entry>xMax</entry><entry>Single</entry><entry>A real number defining the absolute upper limit of x-components of intensity map coordinates (in millimetres). xMin is less than xMax.</entry></row><row><entry /><entry /><entry>Where y' is the maximum x-component of the OD boundary, the actual upper limit of the x-components of intensity map coordinates is the minimum of the following two scalars, xMax and y'+horzBuffer.</entry></row><row rowsep="0"><entry>horzBuffer</entry><entry>Single</entry><entry>A real number bounding x-component coordinates of the IFM intensity map.</entry></row><row rowsep="0"><entry /><entry /><entry>Where y is the minimum x-component of the OD boundary, the actual lower limit of the x-components of intensity map coordinates is the maximum of the following two scalars, xMin and y-horzBuffer.</entry></row><row><entry /><entry /><entry>Where v' is the maximum x-component of the OD boundary, the actual upper limit of the x-components of intensity map coordinates is the minimum of the following two scalars, xMax and y'+horzBuffer.</entry></row><row rowsep="0"><entry>xSpacing</entry><entry>Single</entry><entry>A vector of length n specifying x-spacings at which to iteratively perform computation of intensities at intensity map coordinates. A total of n iterations of intensity computations are performed (assuming intensity computations do not need to be further broken up to satisfy memory requirements).</entry></row><row><entry /><entry /><entry>Coarser spacings are specified toward the lesser indices of the vector while finer spacings are specified toward the finer indices of the vector. Additionally for every i between 1 and n-1, dx(i+1)/dx(i) is a positive integer.</entry></row><row rowsep="0"><entry>zSpacing</entry><entry>Single</entry><entry>A vector of length n specifying z-spacings at which to iteratively perform computation of intensities at intensity map coordinates. A total of n iterations of intensity computations are performed (assuming intensity computations do not need to be further broken up to satisfy memory requirements).</entry></row><row><entry /><entry /><entry>Coarser spacings are specified toward the lesser indices of the vector while finer spacings are specified toward the finer indices of the vector. Additionally for every i between 1 and n-1, dz(i+1)/dz(i) is a positive integer.</entry></row><row><entry>firstDiagIndex</entry><entry>Single</entry><entry>A positive scalar greater than 0 and less than lastdiagindex. This number indicates the first element in the probe array that should be used for the IFM intensity-map computation.</entry></row><row><entry>lastdiagindex</entry><entry>single</entry><entry>A positive scalar greater than firstdiagindex and less than the number of elements in the probe array. This number indicates the last element in the probe array that should be used for the IFM intensity-map computation.</entry></row><row rowsep="0"><entry>zoomPercentages</entry><entry>Single</entry><entry>A vector of length n-1 implicitly specifying coordinates in the intensity map to focus around upon iterative computation of coordinate intensities.</entry></row><row><entry /><entry /><entry>For i between 1 and n-1, zoomPercentage(i) contains a real number between 0 and 1 inclusive. On iteration i+1 of coordinate intensity computation, where Imax is the maximum intensity calculated in the intensity map, the coordinates in which to focus around are those of intensity greater than zoomPercentage(i) * Imax.</entry></row><row rowsep="0"><entry>maxTransmitter</entry><entry morerows="1" rowsep="1">single</entry><entry morerows="1" rowsep="1">A positive scalar defining the maximum absolute difference between a transmitter and receiver used in combination for IFM intensity map computation.</entry></row><row><entry>ReceiverDifference</entry></row><row><entry>criticalPoints</entry><entry>Char</entry><entry>One of the following three strings: 'mins', 'maxes', or 'minsandmaxes.' If criticalPoints is equal to 'mins' (resp. 'maxes'), then only local minimums (resp. maximums) of travel path times are considered to be true travel paths of refracted ultrasonic waves. If criticalPoints is equal to 'minsandmaxes', then both local minimums and local maximums of travel path times are considered to be true travel paths of refracted ultrasonics waves.</entry></row><row><entry>doNormalization</entry><entry>logical</entry><entry>If doNormalization is set to true, then normalization of all A-Scans for IFM computation is executed. Otherwise, if doNormalization is set to false, no normalization is calculated.</entry></row><row rowsep="0"><entry>theshold</entry><entry>Single</entry><entry>A real number greater than 0 and less than or equal to 100.</entry></row><row rowsep="0"><entry /><entry /><entry>This number specifies the cutoff envelope height of wave packets to be normalized in A-Scan data. Where we let 'y' be the maximum envelope value of all FMC data, wavepacket information with peak height greater or equal to threshold*y/100 will be normalized. All other information in the A-Scans will be set to 0.</entry></row><row><entry /><entry /><entry>Recall the envelope of hilbertized data is simply its magnitude.</entry></row><row rowsep="0"><entry>width</entry><entry>Single</entry><entry>A positive scalar greater than 1 and less than or equal to any nonzero length of AScan contained in the FMC data.</entry></row><row><entry /><entry /><entry>This number specifies the distance (in digitization points) between distinct wave packets in an A-Scan. If the distance between envelope peaks greater than the threshold variable exceeds the variable width, then envelope peaks are considered to belong to distinct wave packets. Otherwise, the envelope peaks are considered to belong to the same wave packets.</entry></row><row rowsep="0"><entry>Runup</entry><entry>single</entry><entry>A positive scalar greater than 1 and less than or equal to any nonzero length of AScan contained in the FMC data.</entry></row><row><entry /><entry /><entry>This number specifies the number (in digitization points) of A-Scan values ahead of envelope peaks with values greater than the threshold variable to normalize along with the envelope peaks themselves.</entry></row></tbody></tgroup></table></tables><tables id="tabl0024" num="0024"><table frame="all"><title><b>Table A11 - ID Imaging Functions</b></title><tgroup cols="5"><colspec colnum="1" colname="col1" colwidth="28mm" /><colspec colnum="2" colname="col2" colwidth="20mm" /><colspec colnum="3" colname="col3" colwidth="45mm" /><colspec colnum="4" colname="col4" colwidth="20mm" /><colspec colnum="5" colname="col5" colwidth="56mm" /><thead valign="top"><row><entry><b>Function Name</b></entry><entry><b>Input Data / Input Data Class</b></entry><entry><b>Input Data Description</b></entry><entry><b>Output Data / Output Data Class</b></entry><entry><b>Output Description</b></entry></row></thead><tbody><row rowsep="0"><entry morerows="3" rowsep="1">O_compu tesurface directivityf actor</entry><entry>1. 1D column matrix/single</entry><entry>1. x-coordinates of probe array elements.</entry><entry morerows="3" rowsep="1">1. 2D matrix/logical</entry><entry morerows="3" rowsep="1">Entry (j,k) contains 'true' if the directivity from element j to surface coordinate with index k is greater than or equal to the 'surfaceDirectivityCutoff', defined globally. It contains 'false' otherwise.</entry></row><row rowsep="0"><entry>2. 1D column matrix/single</entry><entry>2. x-coordinates of surface.</entry></row><row rowsep="0"><entry>3. 1D column matrix/single</entry><entry>3. z-coordinates of surface.</entry></row><row><entry>4. scalar/single.</entry><entry>4. wavelength of excitation pulse in water.</entry></row><row rowsep="0"><entry morerows="3" rowsep="1">O_compu tesurface probearra ytimes</entry><entry>1. 1D column matrix/single</entry><entry>1. x-coordinates of surface.</entry><entry>2D matrix/single</entry><entry morerows="3" rowsep="1">Entry (j,k) contains the time needed for a sound pulse to travel from element k in the probe array to surface coordinate with index j.</entry></row><row rowsep="0"><entry>2. 1D column matrix/single</entry><entry>2. z-coordinates of surface.</entry><entry /></row><row rowsep="0"><entry>3. 1D column matrix/single</entry><entry>3. x-coordinates of probe array elements.</entry><entry /></row><row><entry>4. scalar/single.</entry><entry>4. speed of sound in water.</entry><entry /></row><row rowsep="0"><entry morerows="2" rowsep="1">o_comput eidcoordi nates</entry><entry>1. 1D column matrix/single</entry><entry morerows="1">1. x-component of surface coordinates.</entry><entry morerows="1">1. 1D column matrix/ single</entry><entry morerows="1">1. x-component of coordinates of initial intensity computation.</entry></row><row rowsep="0"><entry morerows="1" rowsep="1">2 1D column matrix/single</entry></row><row><entry>2. z-component of surface coordinates.</entry><entry>2. 1D column matrix/ single.</entry><entry>2. z-component of coordinates of initial intensity computation.</entry></row><row rowsep="0"><entry morerows="1" rowsep="1">o_comput ecurrfinec oords4_id</entry><entry>1. 1D column matrix/single</entry><entry morerows="1" rowsep="1">1. x-component of coordinates of previous loop iteration's</entry><entry morerows="1" rowsep="1">1. 1D column matrix/</entry><entry morerows="1" rowsep="1">1. x-coordinate of neighbourhood of intensity map coordinates</entry></row><row><entry>2. 1D column</entry></row><row rowsep="0"><entry /><entry>matrix/single</entry><entry>intensity computation.</entry><entry>single.</entry><entry>exceeding specified intensity cutoff.</entry></row><row rowsep="0"><entry /><entry>3. 1D column matrix/single</entry><entry>2. z-component of coordinates of previous loop</entry><entry morerows="4" rowsep="1">2. ID column matrix/ single.</entry><entry morerows="4" rowsep="1">2. z-coordinate of neighbourhood of intensity map coordinates exceeding specified intensity cutoff.</entry></row><row rowsep="0"><entry /><entry>4. 1D column matrix/single</entry><entry>iteration's intensity computation.</entry></row><row rowsep="0"><entry /><entry>5. positive integer/ double.</entry><entry>3. x-component of surface coordinates.</entry></row><row rowsep="0"><entry /><entry /><entry>4. z-component of surface coordinates.</entry></row><row><entry /><entry /><entry>5. iteration of intensity computation loop.</entry></row><row rowsep="0"><entry morerows="4">O_remov eidcoords outsidebo unds3</entry><entry>1. 1D column matrix/single</entry><entry morerows="2">1. x-component of candidate coordinates of current loop iteration's intensity computation.</entry><entry morerows="1">1. 1D column matrix/ single</entry><entry morerows="2">1. x-component of coordinates of current loop iteration's intensity computation within defined intensity map boundaries.</entry></row><row rowsep="0"><entry>2. 1D column matrix/single</entry></row><row rowsep="0"><entry morerows="1">3. 1D column matrix/single</entry><entry morerows="4" rowsep="1">2. 1D column matrix/ single.</entry></row><row rowsep="0"><entry morerows="1">2. z-component of candidate coordinates of current loop iteration's intensity computation.</entry><entry morerows="3" rowsep="1">2. z-component of coordinates of current loop iteration's intensity computation within defined intensity map boundaries.</entry></row><row rowsep="0"><entry morerows="2" rowsep="1">4. 1D column matrix/single</entry></row><row rowsep="0"><entry /><entry>3. x-component of surface coordinates.</entry></row><row><entry /><entry>4. z-component of surface coordinates.</entry></row><row rowsep="0"><entry morerows="7" rowsep="1">O_computecoordvi sibilityfro msurface 10</entry><entry>1. 1D column matrix/single</entry><entry morerows="2">1. x-component of coordinates of current loop iteration's intensity computation.</entry><entry morerows="1">1. 1D column matrix/ single</entry><entry morerows="1">1. sorted x-component of coordinates of current loop iteration's intensity computation.</entry></row><row rowsep="0"><entry morerows="1">2. 1D column matrix/single</entry></row><row rowsep="0"><entry morerows="1">2. 1D column matrix/ single.</entry><entry morerows="2">2. sorted z-component of coordinates of current loop iteration's intensity computation.</entry></row><row rowsep="0"><entry>3. 1D column matrix/single</entry><entry morerows="2">2. z-component of coordinates of current loop iteration's intensity computation.</entry></row><row rowsep="0"><entry morerows="3" rowsep="1">4. 1D column matrix/single</entry><entry morerows="3" rowsep="1">3. 2D matrix/ logical</entry></row><row rowsep="0"><entry morerows="2" rowsep="1">3. Entry (j,k) is 'true' if there is a straight line from surface coordinate k to intensity map coordinate j which does not intersect the input surface (except at coordinate j). Entry (j,k) is 'false' otherwise.</entry></row><row rowsep="0"><entry>3. x-component of surface coordinates.</entry></row><row><entry>4. z-component of surface coordinates.</entry></row><row rowsep="0"><entry morerows="5" rowsep="1">O_compu tesurface coordtime s</entry><entry>1. 1D column matrix/single</entry><entry morerows="1">1. x-component of coordinates of current loop iteration's intensity computation.</entry><entry morerows="5" rowsep="1">1. 2D matrix/ single.</entry><entry morerows="5" rowsep="1">1. Entry (j,k) is the travel time of sound in water from element k to surface coordinate with index j.</entry></row><row rowsep="0"><entry>2. 1D column matrix/single</entry></row><row rowsep="0"><entry>3. 1D column matrix/single</entry><entry morerows="1">2. z-component of coordinates of current loop iteration's intensity computation.</entry></row><row rowsep="0"><entry morerows="2" rowsep="1">4. 1D column matrix/single</entry></row><row rowsep="0"><entry>3. x-component of surface coordinates.</entry></row><row><entry>4. z-component of surface coordinates.</entry></row><row><entry /><entry /><entry /><entry /><entry /></row><row rowsep="0"><entry morerows="9" rowsep="1">O_compu tetravelti mes4</entry><entry>1. 2D matrix/single</entry><entry morerows="2">1. Entry (j,k) is the travel time of sound in water from element k to surface coordinate with index j.</entry><entry>1. Structure with fields:</entry><entry morerows="2">With respect to index i of the structure (corresponding to element i of the linear probe array),</entry></row><row rowsep="0"><entry>2. 2D matrix/single</entry><entry morerows="1">a) column matrix/unsign ed int (16-bit).</entry></row><row rowsep="0"><entry>3. 2D</entry></row><row rowsep="0"><entry>matrix/logical</entry><entry morerows="2">2. Entry (j,k) contains the time needed for a sound pulse to travel from element k in the probe array to surface coordinate with index j.</entry><entry>b) column matrix/single</entry><entry>a) entry j contains theindex of a coordinate</entry></row><row rowsep="0"><entry morerows="5" rowsep="1">4. 2D matrix/logical</entry><entry>c) column matrix/unsign ed integer (8 bit)</entry><entry>b) entry j contains thetravel time (in samples) from element i to the coordinate with index j in field a).</entry></row><row rowsep="0"><entry morerows="1">d) column matrix/unsign ed integer (8-bit)</entry><entry morerows="1">c) entry j contains the number travel-time solutions from element I to the coordinate with index j in field a).</entry></row><row rowsep="0"><entry morerows="2">3. Entry (j,k) contains 'true' if the directivity from element j to surface coordinate with index k is greater than or equal to the 'surfaceDirectivity Cutoff', defined globally. It contains 'false'</entry></row><row rowsep="0"><entry>e) scalar /unsigned integer (8-bit)</entry><entry>d) contains the set of number of solutions from element i to all inspection coordinates.</entry></row><row rowsep="0"><entry /><entry>e) contains the length of field d).</entry></row><row><entry>4. Entry (j,k) is 'true' if there is a straight line from surface coordinate k to intensity map coordinate j which does not intersect the input surface (except at coordinate j). Entry (j,k) is 'false' otherwise.</entry><entry /><entry /></row><row rowsep="0"><entry>O_ifmzoo m8</entry><entry>1 . 2x1 element matrix/doubl e</entry><entry>1. Size of coordinate matrices</entry><entry morerows="12" rowsep="1">1. column matrix/single</entry><entry morerows="12" rowsep="1">1. computed intensities for coordinates with matching indices as computed intensities.</entry></row><row rowsep="0"><entry /><entry>2. 2D matrix/ complex single</entry><entry>2. FMC data-set where A-Scans are contained in</entry></row><row rowsep="0"><entry /><entry>3. 2D matrix/single</entry><entry>the vertical dimension.</entry></row><row rowsep="0"><entry /><entry>4. 2D matrix/single</entry><entry morerows="1">3. Entry (i,j) contains the column index of the A-Scan belonging to transmitter i and receiver j.</entry></row><row rowsep="0"><entry /><entry>5. Structure with fields:</entry></row><row rowsep="0"><entry /><entry>a) column matrix/unsigned int (16-bit).</entry><entry morerows="3">4. Entry (i,j) contains the offset (in samples) between transmission of ultrasonic wave from element i and recording of A-Scan data from receiver j.</entry></row><row rowsep="0"><entry /><entry>b) column matrix/single</entry></row><row rowsep="0"><entry /><entry>c) column matrix/unsigned integer (8 bit)</entry></row><row rowsep="0"><entry /><entry morerows="1">d) column matrix/unsigned integer (8-bit)</entry></row><row rowsep="0"><entry /><entry morerows="1">5. With respect to index i of the structure (corresponding to element i of the linear probe array),</entry></row><row rowsep="0"><entry /><entry morerows="2" rowsep="1">e) scalar /unsigned integer (8-bit)</entry></row><row rowsep="0"><entry /><entry>a) entry j contains the index of a coordinate</entry></row><row><entry /><entry>b) entry j contains the travel time (in samples) from</entry></row><row rowsep="0"><entry /><entry /><entry>element i to the coordinate with index j in field a).</entry><entry /><entry /></row><row rowsep="0"><entry /><entry /><entry>c) entry j contains the number travel-time solutions from element I to the coordinate with index j in field a).</entry><entry /><entry /></row><row rowsep="0"><entry /><entry /><entry>d) contains the set of number of solutions from element i to all inspection coordinates.</entry><entry /><entry /></row><row><entry /><entry /><entry>e) contains the length of field d).</entry><entry /><entry /></row><row rowsep="0"><entry morerows="4" rowsep="1">O_compu tecoordin atesabov ecutoff_id</entry><entry>1. column matrix/single</entry><entry>1. coordinate x-components.</entry><entry>1. column matrix! single.</entry><entry morerows="1">1. coordinate x-components with intensity values exceeding cutoff intensity value (input item 4).</entry></row><row rowsep="0"><entry morerows="1">2. column matrix/single</entry><entry>2. coordinate z-components.</entry><entry morerows="3" rowsep="1">2. column matrix/ single.</entry></row><row rowsep="0"><entry>3. coordinate intensities.</entry><entry morerows="2" rowsep="1">2. coordinate z-components with intensity values exceeding cutoff intensity value (input item 4).</entry></row><row rowsep="0"><entry>3. column matrix/single</entry><entry morerows="1" rowsep="1">4. cutoff intensity value.</entry></row><row><entry>4. scalar/single.</entry></row></tbody></tgroup></table></tables><tables id="tabl0025" num="0025"><table frame="all"><title><b>Table A12</b> - <b>2" Reference Block (Serial Number 001) Dimensions</b></title><tgroup cols="5"><colspec colnum="1" colname="col1" colwidth="15mm" /><colspec colnum="2" colname="col2" colwidth="16mm" /><colspec colnum="3" colname="col3" colwidth="57mm" /><colspec colnum="4" colname="col4" colwidth="38mm" /><colspec colnum="5" colname="col5" colwidth="21mm" /><thead valign="top"><row><entry namest="col1" nameend="col2" align="left">Location</entry><entry>Measurement converted from Imperial</entry><entry>Rounded Measurement</entry><entry>Step Delta</entry></row><row><entry>Scan</entry><entry>Block</entry><entry /><entry /><entry /></row></thead><tbody><row><entry /><entry>1</entry><entry>2.032</entry><entry>2.03</entry><entry /></row><row><entry /><entry>2</entry><entry>5.0038</entry><entry>5.00</entry><entry /></row><row><entry /><entry>3</entry><entry>8.00608</entry><entry>8.01</entry><entry /></row><row><entry /><entry>4</entry><entry>10.9982</entry><entry>11.00</entry><entry /></row><row><entry /><entry>5</entry><entry>13.9827</entry><entry>13.98</entry><entry /></row><row><entry /><entry>6</entry><entry>16.99006</entry><entry>16.99</entry><entry /></row><row><entry /><entry>7</entry><entry>23.89632</entry><entry>23.90</entry><entry>3.00</entry></row><row><entry /><entry>8</entry><entry>20.89912</entry><entry>20.90</entry><entry>3.00</entry></row><row><entry /><entry>9</entry><entry>17.89684</entry><entry>17.90</entry><entry>3.00</entry></row><row><entry /><entry>10</entry><entry>14.89964</entry><entry>14.90</entry><entry>3.00</entry></row><row><entry /><entry>11</entry><entry>11.89736</entry><entry>11.90</entry><entry>2.97</entry></row><row><entry /><entry>12</entry><entry>8.9281</entry><entry>8.93</entry><entry /></row></tbody></tgroup></table></tables><tables id="tabl0026" num="0026"><table frame="all"><title><b>Table A13 - 2" Reference Block (Serial Number 002) Dimensions</b></title><tgroup cols="5"><colspec colnum="1" colname="col1" colwidth="15mm" /><colspec colnum="2" colname="col2" colwidth="16mm" /><colspec colnum="3" colname="col3" colwidth="57mm" /><colspec colnum="4" colname="col4" colwidth="38mm" /><colspec colnum="5" colname="col5" colwidth="21mm" /><thead valign="top"><row><entry namest="col1" nameend="col2" align="left">Location</entry><entry>Measurement converted from Imperial</entry><entry>Rounded Measurement</entry><entry>Step Delta</entry></row><row><entry>Scan</entry><entry>Block</entry><entry /><entry /><entry /></row></thead><tbody><row><entry /><entry>1</entry><entry>2.032</entry><entry>2.03</entry><entry /></row><row><entry /><entry>2</entry><entry>5.00888</entry><entry>5.01</entry><entry /></row><row><entry /><entry>3</entry><entry>8.0264</entry><entry>8.03</entry><entry /></row><row><entry /><entry>4</entry><entry>11.0109</entry><entry>11.01</entry><entry /></row><row><entry /><entry>5</entry><entry>13.99286</entry><entry>13.99</entry><entry /></row><row><entry /><entry>6</entry><entry>17.00022</entry><entry>17.00</entry><entry /></row><row><entry /><entry>7</entry><entry>23.89378</entry><entry>23.89</entry><entry>2.99</entry></row><row><entry /><entry>8</entry><entry>20.9042</entry><entry>20.90</entry><entry>3.01</entry></row><row><entry /><entry>9</entry><entry>17.8943</entry><entry>17.89</entry><entry>3.00</entry></row><row><entry /><entry>10</entry><entry>148971</entry><entry>14.90</entry><entry>3.00</entry></row><row><entry /><entry>11</entry><entry>11.8999</entry><entry>11.90</entry><entry>2.98</entry></row><row><entry /><entry>12</entry><entry>8.9154</entry><entry>8.92</entry><entry /></row></tbody></tgroup></table></tables><tables id="tabl0027" num="0027"><table frame="all"><title><b>Table A14-2.5" Reference Block (Serial Number 001) Dimensions</b></title><tgroup cols="5"><colspec colnum="1" colname="col1" colwidth="15mm" /><colspec colnum="2" colname="col2" colwidth="16mm" /><colspec colnum="3" colname="col3" colwidth="58mm" /><colspec colnum="4" colname="col4" colwidth="38mm" /><colspec colnum="5" colname="col5" colwidth="19mm" /><thead valign="top"><row><entry namest="col1" nameend="col2" align="left">Location</entry><entry>Measurement converted from Imperial</entry><entry>Rounded Measurement</entry><entry>Step Delta</entry></row><row><entry>Scan</entry><entry>Block</entry><entry /><entry /><entry /></row></thead><tbody><row><entry align="center" /><entry align="center">1</entry><entry align="center">1.99898</entry><entry align="center">2.00</entry><entry align="center" /></row><row><entry align="center" /><entry align="center">2</entry><entry align="center">4.99872</entry><entry align="center">5.00</entry><entry align="center" /></row><row><entry align="center" /><entry align="center">3</entry><entry align="center">8.001</entry><entry align="center">8.00</entry><entry align="center" /></row><row><entry align="center" /><entry align="center">4</entry><entry align="center">10.9982</entry><entry align="center">11.00</entry><entry align="center" /></row><row><entry align="center" /><entry align="center">5</entry><entry align="center">13.9954</entry><entry align="center">14.00</entry><entry align="center" /></row><row><entry align="center" /><entry align="center">6</entry><entry align="center">16.99768</entry><entry align="center">17.00</entry><entry align="center" /></row><row><entry align="center" /><entry align="center">7</entry><entry align="center">24.0919</entry><entry align="center">24.09</entry><entry align="center">3.00</entry></row><row><entry align="center" /><entry align="center">8</entry><entry align="center">21.0947</entry><entry align="center">21.10</entry><entry align="center">3.01</entry></row><row><entry align="center" /><entry align="center">9</entry><entry align="center">18.0848</entry><entry align="center">18.09</entry><entry align="center">2.99</entry></row><row><entry align="center" /><entry align="center">10</entry><entry align="center">15.09522</entry><entry align="center">15.10</entry><entry align="center">3.00</entry></row><row><entry align="center" /><entry align="center">11</entry><entry align="center">12.09802</entry><entry align="center">12.10</entry><entry align="center">3.01</entry></row><row><entry align="center" /><entry align="center">12</entry><entry align="center">9.0932</entry><entry align="center">9.09</entry><entry align="center" /></row></tbody></tgroup></table></tables><tables id="tabl0028" num="0028"><table frame="all"><title><b>Table A15-2.5" Reference Block (Serial Number 002) Dimensions</b></title><tgroup cols="5"><colspec colnum="1" colname="col1" colwidth="15mm" /><colspec colnum="2" colname="col2" colwidth="16mm" /><colspec colnum="3" colname="col3" colwidth="58mm" /><colspec colnum="4" colname="col4" colwidth="38mm" /><colspec colnum="5" colname="col5" colwidth="19mm" /><thead valign="top"><row><entry namest="col1" nameend="col2" align="left">Location</entry><entry morerows="1">Measurement converted from Imperial</entry><entry morerows="1">Rounded Measurement</entry><entry morerows="1">Step Delta</entry></row><row><entry>Scan</entry><entry>Block</entry></row></thead><tbody><row><entry /><entry align="center">1</entry><entry>2.00406</entry><entry>2.00</entry><entry /></row><row><entry /><entry align="center">2</entry><entry>4.99364</entry><entry>4.99</entry><entry /></row><row><entry /><entry align="center">3</entry><entry>7.99592</entry><entry>8.00</entry><entry /></row><row><entry /><entry align="center">4</entry><entry>10.9982</entry><entry>11.00</entry><entry /></row><row><entry /><entry align="center">5</entry><entry>13.9954</entry><entry>14.00</entry><entry /></row><row><entry /><entry align="center">6</entry><entry>16.9926</entry><entry>1699</entry><entry /></row><row><entry /><entry align="center">7</entry><entry>24.09698</entry><entry>24.10</entry><entry>3.00</entry></row><row><entry /><entry align="center">8</entry><entry>21.10232</entry><entry>21.10</entry><entry>3.01</entry></row><row><entry /><entry align="center">9</entry><entry>18.0975</entry><entry>18.10</entry><entry>2.99</entry></row><row><entry /><entry align="center">10</entry><entry>15.10538</entry><entry>15.11</entry><entry>3.01</entry></row><row><entry /><entry align="center">11</entry><entry>12.09294</entry><entry>12.09</entry><entry>3.00</entry></row><row><entry /><entry align="center">12</entry><entry>9.0932</entry><entry>9.09</entry><entry /></row></tbody></tgroup></table></tables>