JP5380237B2 - Ultrasound image diagnostic apparatus and analysis region setting program - Google Patents

Description

  The present invention generates an ultrasonic image of a subject by transmitting and receiving an ultrasonic beam by an ultrasonic probe, and has a function of setting the position and size of an analysis region for analyzing the motion of a specific part in the subject. The present invention relates to a sound image diagnostic apparatus and an analysis region setting program.

  In the ultrasonic diagnostic imaging apparatus, for example, measurement using the Doppler effect is performed in order to perform functional diagnosis of the heart. What is measured in this way is the blood flow from the heart based on the Doppler signal output from the ultrasound probe, the backflow from the valve, the expansion function, and the like. For example, in the expansion function, functional diagnosis is performed by measuring the elongation of the myocardium according to the blood flow from the left atrium of the heart based on the Doppler waveform of the cardiac muscle (referred to as TDI-PW). In the expansion function, functional diagnosis is also performed by measuring a blood inflow waveform pattern from the mitral valve based on a Doppler waveform (referred to as E / A measurement) near the mitral valve of the heart. As techniques for measuring the function of the heart, there are, for example, Patent Documents 1 and 2. For example, Non-Patent Document 1 discloses a technique related to measuring the expansion / contraction motion of the heart in the left ventricle of the heart using a Doppler waveform and measuring the expansion function of the heart.

JP 2002-140690 A Japanese Patent Laid-Open No. 10-99334

Sohn, D.D. W. , Et al: J. Am. Coll. Cardiol., No.30: pp474-480.1997

  In order to perform measurement using the Doppler effect in order to perform functional diagnosis of the heart or the like as described above, each Doppler waveform such as TDI-PW or E / A measurement must be measured. In order to measure these Doppler waveforms, a sampling mark indicating the measurement position is set on the ultrasonic image. This sampling mark is formed in, for example, a rectangular shape, and its position and width can be adjusted. Based on the Doppler signal in the sampling mark area, for example, the measurement of the myocardial elongation and the measurement of the blood inflow waveform pattern from the mitral valve are performed.

  However, the position and width of the sampling mark must be reset every measurement of each Doppler waveform, and the sampling mark must be angle-corrected in order to accurately measure the blood flow velocity for each measurement. For this reason, it takes time and effort to set the position and width of the sampling mark and angle correction for each measurement.

  An object of the present invention is to provide an ultrasonic diagnostic imaging apparatus and an analysis area setting program that can automatically improve the operability and diagnostic efficiency by automatically setting the position and width of a sampling mark and angle correction for each measurement. There is to do.

  An ultrasonic diagnostic imaging apparatus of the present invention corresponding to claim 1 is an ultrasonic diagnostic imaging apparatus that generates an ultrasonic image of a subject by transmitting and receiving an ultrasonic beam by an ultrasonic probe, and the subject is detected from an ultrasonic image. The tracking processing unit that tracks the movement of the specific part in the system, and the position and size information of the analysis region for analyzing the function of the specific part based on the tracking result of the tracking processing part, and the movement speed of the specific part And an analysis region setting processing unit for setting information on angle correction of the analysis region with respect to the central axis of the ultrasonic beam for calculating.

  An analysis region setting program according to the present invention corresponding to claim 8 is an analysis region for analyzing a function of a specific part of a subject on an ultrasonic image of the subject generated by transmission / reception of an ultrasonic beam by an ultrasonic probe. An analysis region setting program for setting the position and size of a tracking region, which allows a computer to track the movement of a specific part of a subject from an ultrasound image, and the position and size of the analysis region based on the tracking result In addition to setting information, a setting function for setting information on angle correction of the analysis region with respect to the central axis of the ultrasonic beam for calculating the motion speed of the specific part is realized.

  According to the present invention, an ultrasonic diagnostic imaging apparatus and an analysis region setting program capable of automatically setting the position and width of a sampling mark and angle correction for each measurement and improving operability and diagnostic efficiency are provided. it can.

1 is a configuration diagram showing an embodiment of an ultrasonic diagnostic imaging apparatus according to the present invention. The figure which shows the setting of the sampling mark on the ultrasonic image by the apparatus. The figure which shows the electrocardiogram signal used for the apparatus. The schematic diagram which shows the reference image data of the annulus memorize | stored in the annulus tracking process part in the same apparatus. The figure which shows the shape of the valve ring produced | generated by the valve ring tracking process part in the same apparatus. The figure for demonstrating the setting of the information of the position and magnitude | size of a sampling mark by the sampling mark control process part in the same apparatus, and the setting of the information of angle correction of a sampling mark. The figure for demonstrating the setting of the width | variety of the sampling mark by the sampling mark control process part in the apparatus. The figure for demonstrating the angle correction | amendment by the sampling mark control process part in the same apparatus.

Hereinafter, an embodiment of the present invention will be described with reference to the drawings.
FIG. 1 shows a configuration diagram of an ultrasonic diagnostic imaging apparatus. The ultrasonic probe 1 transmits an ultrasonic beam to a subject Q such as a human body, receives a reflected wave from the subject Q, and outputs a Doppler signal.
The transmission / reception circuit 2 applies an electrical pulse to each transducer of the ultrasonic probe 1, scans an ultrasonic beam from the ultrasonic probe 1 to, for example, the subject Q, and transmits the ultrasonic beam from the ultrasonic probe 1. The echo signal from the subject Q when scanning is received and the Doppler signal is output.
The signal processing circuit 3 outputs an ultrasonic signal for displaying an ultrasonic image on the monitor 4 from the output signal of the transmission / reception circuit 2, and includes a preprocess circuit 5, a digital scan converter circuit 6, and a post process circuit 7. And have. The preprocess circuit 5 amplifies the output signal of the transmission / reception circuit 2. The digital scan converter circuit 6 receives the output signal of the preprocess circuit 5 and converts it into ultrasonic two-dimensional image data for projecting and displaying three-dimensional image data acquired by scanning the ultrasonic probe 1 on a two-dimensional surface. Convert. The post-process circuit 7 performs balance, γ correction, color mapping, and the like on the ultrasonic two-dimensional image data converted by the digital scan converter circuit 6 to generate ultrasonic data US.

  The electrocardiograph 8 outputs a biological signal E representing the electrical phenomenon of the heart of the subject Q. The biological circuit 9 generates an electrocardiogram signal ECG from the biological signal E output from the electrocardiograph 8. A peak named R wave appears in the waveform of the electrocardiogram signal ECG.

  The CPU 11 includes a program memory 12. The program memory 12 stores an analysis area setting program in advance. This analysis area setting program is an analysis area for analyzing the function of a specific part such as a heart in the subject Q on the ultrasonic image of the subject Q generated by transmission and reception of the ultrasonic beam of the ultrasonic probe 1. For example, the processing for setting the position and size of the sampling mark S as shown in FIG. 2 is realized.

This analysis region setting program tracks the movement of the mitral valve M in the left ventricle LV of the heart as shown in FIG. 2, for example, as shown in FIG. Are set, and information on the angle correction of the sampling mark S with respect to the central axis of the ultrasonic beam for calculating the motion speed of the specific part is set.
In the analysis region setting program, the movement tracking of a specific part tracks the position of the annulus Ma from an ultrasonic image including the annulus Ma of the heart as a subject.
The processing for setting information such as the position of the sampling mark S in the analysis region setting program is performed by setting the sampling mark S formed in a rectangular shape on the ultrasonic image and with respect to the origin position where the ultrasonic beam is transmitted. The midpoint between the first position where the subject Q is closest and the second position where the subject Q is farthest is set as the center position of the sampling mark S, and then the first position connecting the first position and the second position. The angle formed by the straight line 1 and the second straight line connecting the center position of the sampling mark S and the origin position where the ultrasonic beam is transmitted is set as angle correction information, and then the first and second positions Setting the length of the sampling mark S in the longitudinal direction based on the center position of the sampling mark S and the angle correction information is realized.
The sampling mark S is formed in a rectangular shape. As shown in FIG. 2, for example, when the XY coordinates having the origin G as the upper left coordinate (0, 0) are set on the ultrasonic image as shown in FIG. The The coordinates (X0, Y0) are the origin of the ultrasonic probe 1.

  The CPU 11 executes an analysis area setting program stored in the program memory 12 to thereby analyze the position of an analysis area (hereinafter referred to as a sampling mark) for analyzing the function of a specific part in the subject Q, such as the heart. And an analysis region setting processing unit 10 for setting the size and the like.

Specifically, the analysis region setting processing unit 10 has functions of an annulus tracking processing unit 13 and a sampling mark control processing unit (analysis region control processing unit) 14. The CPU 11 has functions of an input processing unit 15 and a hardware control processing unit 16.
A console 17 having a plurality of operation buttons and the like is connected to the input processing unit 15. The console 17 includes, for example, an operation button for giving an instruction to measure the expansion / contraction motion of the heart muscle by the Doppler effect. The input processing unit 15 issues a command for setting the sampling mark S to the annulus tracking processing unit 13 when receiving an instruction to measure the expansion / contraction motion of the heart's myocardium by the Doppler effect from the operation button.

  The annulus tracking processing unit 13 inputs the ultrasonic two-dimensional image data converted by the digital scan converter circuit 6, and also inputs an electrocardiogram signal ECG as shown in FIG. The position of the heart annulus Ma is tracked on the ultrasonic two-dimensional image data from the time when the R wave appears on the ECG until a predetermined period elapses. Here, the timing at which the R wave appears in the electrocardiogram signal ECG corresponds to the end diastole of the heart. The period from the end-diastolic timing of the heart to the elapse of a predetermined period is a time K that is a half of the systolic phase of the heart or one heartbeat from the end-diastolic timing at which the R wave appears.

  Specifically, the annulus tracking processing unit 13 stores reference image data P of a heart annulus Ma used for pattern matching, for example. The annulus tracking processing unit 13 performs pattern matching between the ultrasonic two-dimensional image data and the reference image data P in a preset search region, and detects the position of the annulus Ma on the ultrasonic two-dimensional image data. Since the annulus Ma of the mitral valve is formed in an annular shape, the position of the annulus Ma is detected by forming a preset search region in a donut shape, for example. And the position of the annulus Ma is detected in these divided areas. Thereby, the shape of the annulus Ma as shown in FIG. 5 is generated.

  Next, after generating the shape of the annulus Ma, the annulus tracking processing unit 13 calculates values such as the amount of movement, speed, and acceleration due to the heartbeat for each position of the annulus Ma detected by pattern matching. To do. Thereby, the tracking result of the moving position of the annulus Ma is obtained.

The sampling mark control processing unit 14 sets, for example, information on the position and size of the sampling mark S for analyzing the function of the heart based on the tracking result of the annulus tracking processing unit 13, and sets the movement speed of the specific part. Information on angle correction of the sampling mark S with respect to the central axis (raster direction) H of the ultrasonic beam for calculation is set.
Specifically, the sampling mark control processing unit 14 sets the sampling mark S on the ultrasonic image as shown in FIG. 6, and the ultrasonic beam is compared with the origin coordinates (X0, Y0) for transmitting the ultrasonic beam. On the central axis, the coordinates (X1, Y1) of the first position where the annulus Ma tracking position is closest and the coordinates (X2, Y2) of the second position which are farthest are detected, and these first positions are detected. The coordinates (Xm, Ym) of the midpoint between the coordinates (X1, Y1) and the coordinates (X2, Y2) of the second position are set as the center position of the sampling mark S.
Xm = (X1 + X2) / 2 (1)
Ym = (Y1 + Y2) / 2 (2)
The sampling mark control processing unit 14 includes the first straight line F1 connecting the coordinates (X1, Y1) of the first position and the coordinates (X2, Y2) of the second position, and the coordinates of the center position of the sampling mark S ( Xm, Ym) and the angle θ formed by the second straight line F2 connecting the origin coordinates (X0, Y0) for transmitting the ultrasonic beam are set as angle correction information.
The sampling mark control processing unit 14 includes coordinates (X1, Y1) of the first position, coordinates (X2, Y2) of the second position, coordinates (Xm, Ym) of the center position of the sampling mark S, and angle correction information θ. Based on the above, the length in the longitudinal direction of the sampling mark S (also referred to as the width of the sampling mark S) is set. For example, as shown in FIG. 7, the width of the sampling mark S includes an upper side position Up and a lower side position Dp. The position Up of the upper side is the position closest to the origin coordinate (X0, Y0) side where the ultrasonic beam is transmitted among the tracking results of the annulus Ma. The position Dp of the lower side is a position farthest from the coordinates (X0, Y0) of the origin where the ultrasonic beam is transmitted most among the annulus Ma tracking results.

In this case, the sampling mark control processing unit 14 determines the movement speed V ′ of the annulus Ma obtained based on the Doppler signal output from the ultrasonic probe 1, the angle correction information θ, and the true movement speed of the annulus Ma. Relationship with V,
V = V ′ / cos θ (3)
In consideration of the above, the upper side position Up ′ and the lower side position Dp ′ of the sampling mark S are set. At this time, the sampling mark control processing unit 14 generates an angle correction mark AH having an upper side position Up ′ and a lower side position Dp ′.
That is, the sampling mark control processing unit 14 calculates L1 · cos θ with respect to the distance L1 between the coordinates (X1, Y1) of the first position and the coordinates (Xm, Ym) of the center position of the sampling mark S. Thus, the position Up ′ of the upper side portion of the sampling mark S from the center position of the sampling mark S toward the origin for transmitting the ultrasonic beam is obtained.
At the same time, the sampling mark control processing unit 14 calculates L2 · cos θ with respect to the distance L2 between the coordinates (X2, Y2) of the second position and the coordinates (Xm, Ym) of the center position of the sampling mark S. Thus, the position Dp ′ of the lower side portion on the raster H of the ultrasonic beam is obtained from the center position of the sampling mark S.

  The mixing circuit 18 includes the ultrasonic data US generated by the signal processing circuit 3, the center position (Xm, Ym) and the width (upper side position Up ′ of the sampling mark S generated by the analysis region setting processing unit 10, The position Dp ′) of the lower side and the angle correction information θ are mixed and sent to the monitor 4.

Next, the setting operation of the sampling mark by the apparatus configured as described above will be described.
First, the ultrasonic probe 1 is applied to the subject Q corresponding to the apex of the heart. The ultrasonic probe 1 transmits an ultrasonic beam to a subject Q such as a human body by the operation of the transmission / reception circuit 2, receives a reflected wave from the subject Q, and outputs a Doppler signal.

  The preprocess circuit 5 amplifies the output signal of the transmission / reception circuit 2. The digital scan converter circuit 6 receives the output signal of the preprocess circuit 5 and converts it into ultrasonic two-dimensional image data for projecting and displaying three-dimensional image data acquired by scanning the ultrasonic probe 1 on a two-dimensional surface. Convert. The post-process circuit 7 performs balance, γ correction, color mapping, and the like on the ultrasonic two-dimensional image data converted by the digital scan converter circuit 6 to generate ultrasonic data US. Thereby, the ultrasonic data US is sent to the monitor 4 through the mixing circuit 18. On the monitor 4, for example, a four-chamber tomographic image of the heart in the subject Q as shown in FIG. 2 is displayed.

When an operation button for instructing to measure the expansion / contraction motion of the heart's myocardium on the console 17 by the Doppler effect is operated, the input processing unit 15 issues a command for setting the sampling mark S to the annulus tracking processing unit. Depart to 13.
When the annulus tracking processing unit 13 receives a command for setting the sampling mark S, the annulus tracking processing unit 13 inputs the ultrasonic two-dimensional image data converted by the digital scan converter circuit 6 and generates a diagram generated by the biological circuit 9. The ECG signal ECG as shown in FIG. 3 is input, and the position of the heart annulus Ma is tracked on the ultrasound two-dimensional image data from the timing when the R wave appears in the ECG signal ECG until a predetermined period elapses. Specifically, the annulus tracking processing unit 13 performs pattern matching between the ultrasonic two-dimensional image data and the reference image data P in a preset search region, and determines the position of the annulus Ma on the ultrasonic two-dimensional image data. To detect. Thereby, the shape of the annulus Ma as shown in FIG. 5 is generated.

  Next, after generating the shape of the annulus Ma, the annulus tracking processing unit 13 calculates values such as the amount of movement, speed, and acceleration due to the heartbeat for each position of the annulus Ma detected by pattern matching. To do. Thereby, the tracking result of the moving position of the annulus Ma is obtained.

Next, the sampling mark control processing unit 14 sets, for example, information on the position and size of the sampling mark S for analyzing the function of the heart based on the tracking result of the annulus tracking processing unit 13, and Information on the angle correction of the sampling mark S with respect to the central axis (raster direction) H of the ultrasonic beam for calculating the motion speed is set.
Specifically, the sampling mark control processing unit 14 sets the sampling mark S on the ultrasonic image as shown in FIG. 6, and performs the annulus Ma tracking with respect to the origin coordinates (X0, Y0) for transmitting the ultrasonic beam. The first position (X1, Y1) at which the position is closest and the second position (X2, Y2) at which the position is farthest are detected, and the first position (X1, Y1) and the second position (X2) are detected. , Y2) is set at the center position of the sampling mark S (Xm, Ym).

Next, the sampling mark control processing unit 14 includes the first straight line F1 connecting the first position (X1, Y1) and the second position (X2, Y2), and the center position (Xm, Ym) of the sampling mark S. ) And the second straight line F2 connecting the origin coordinates (X0, Y0) for transmitting the ultrasonic beam is set as angle correction information.
Next, the sampling mark control processing unit 14 is based on the first position (X1, Y1), the second position (X2, Y2), the center position (Xm, Ym) of the sampling mark S, and the angle correction information θ. To set the width of the sampling mark S. For example, as shown in FIG. 7, the width of the sampling mark S is determined by the position Up of the upper side closest to the origin coordinate (X0, Y0) side where the ultrasonic beam is transmitted among the tracking results of the annulus Ma, and the annulus Ma. It consists of a position Dp on the lower side farthest from the origin coordinates (X0, Y0) where the ultrasonic beam is transmitted most among the tracking results.

In this case, the sampling mark control processing unit 14 determines the movement speed V ′ of the annulus Ma obtained based on the Doppler signal output from the ultrasonic probe 1, the angle correction information θ, and the true movement speed of the annulus Ma. Relationship with V,
V = V ′ / cos θ (3)
In consideration of the above, the upper side position Up ′ and the lower side position Dp ′ of the sampling mark S are set.

That is, the sampling mark control processing unit 14 calculates L1 · cos θ with respect to the distance L1 between the coordinates (X1, Y1) of the first position and the coordinates (Xm, Ym) of the center position of the sampling mark S. Thus, the position Up ′ of the upper side portion of the sampling mark S from the center position of the sampling mark S toward the origin for transmitting the ultrasonic beam is obtained.
At the same time, the sampling mark control processing unit 14 calculates L2 · cos θ with respect to the distance L2 between the coordinates (X2, Y2) of the second position and the coordinates (Xm, Ym) of the center position of the sampling mark S. Thus, the position Dp ′ of the lower side portion on the raster H of the ultrasonic beam is obtained from the center position of the sampling mark S.

  The mixing circuit 18 includes the ultrasonic data US generated by the signal processing circuit 3, the center position (Xm, Ym) and the width (upper side position Up ′ of the sampling mark S generated by the analysis region setting processing unit 10, The position Dp ′) of the lower side and the angle correction information θ are mixed and sent to the monitor 4. Accordingly, the monitor 4 is set with a width by the ultrasonic two-dimensional image data, and the upper side position Up ′ and the lower side position Dp ′ set on the ultrasonic two-dimensional image data, and the angle is set. A sampling mark S whose angle is corrected by the correction information θ is displayed.

  As described above, according to the above-described embodiment, the position of the annulus Ma of the heart is tracked on the ultrasonic two-dimensional image data, and the sampling mark S for analyzing the function of the heart, for example, is analyzed based on the tracking result. Information on the position and size is automatically set, and information on the angle correction of the sampling mark S with respect to the central axis (raster direction) H of the ultrasonic beam for calculating the motion speed of the specific part is automatically set. . Thereby, in functional diagnosis of the heart or the like, each Doppler waveform such as TDI-PW or E / A measurement is measured, but the position and width of the sampling mark S are not reset every time the Doppler waveform is measured. In addition, it is possible to improve the operability during functional diagnosis of the heart and the like and improve the diagnostic efficiency without correcting the angle of the sampling mark S in order to accurately measure the blood flow velocity in the heart chamber for each measurement. Can do. Based on the Doppler signal in the region of the sampling mark S, for example, the measurement of the myocardial elongation according to the blood flow from the left atrium of the heart, the measurement of the blood inflow waveform pattern from the mitral valve, and the like can be performed efficiently. In addition, the burden of diagnosis on a subject such as a patient can be reduced.

Note that the present invention is not limited to the above-described embodiment as it is, and can be embodied by modifying the constituent elements without departing from the scope of the invention in the implementation stage. In addition, various inventions can be formed by appropriately combining a plurality of components disclosed in the embodiment. For example, some components may be deleted from all the components shown in the embodiment. Furthermore, constituent elements over different embodiments may be appropriately combined.
For example, in the above-described embodiment, the position of the heart annulus Ma is tracked. However, the present invention is not limited to this, and the position of the aortic valve may be tracked.

  1: Ultrasonic probe, Q: Subject, 2: Transmission / reception circuit, 3: Signal processing circuit, 4: Monitor, 5: Preprocess circuit, 6: Digital scan converter circuit, 7: Post process circuit, 8: Electrocardiograph 10: analysis region setting processing unit, 11: CPU, 12: program memory, 13: annulus tracking processing unit, 14: sampling mark control processing unit (analysis region control processing unit), 15: input processing unit 1, 16: Hardware control processing unit, 17: console.

Claims (9)

In an ultrasonic diagnostic imaging apparatus that generates an ultrasonic image of a subject by transmitting and receiving an ultrasonic beam with an ultrasonic probe,
A tracking processor that tracks the movement of a specific part of the subject from the ultrasound image;
Based on the tracking result of the tracking processing unit, information on the position and size of the analysis region for analyzing the function of the specific part is set, and the ultrasonic beam for calculating the motion speed of the specific part is set. An analysis region setting processing unit for setting information of angle correction of the analysis region with respect to the central axis;
An ultrasonic diagnostic imaging apparatus comprising:   The ultrasonic image diagnostic apparatus according to claim 1, wherein the tracking processing unit tracks the position of the annulus from the ultrasonic image including a heart annulus as the subject.   The ultrasonic image diagnostic apparatus according to claim 2, wherein the tracking processing unit tracks the position of the annulus within a predetermined period based on the ultrasonic image and the biological signal of the subject.   The analysis area control processing unit sets a sampling mark formed in a rectangular shape as the analysis area on the ultrasonic image, and is on a central axis of the ultrasonic beam with respect to an origin position where the ultrasonic beam is transmitted. 2. The ultrasonic diagnostic imaging apparatus according to claim 1, wherein a midpoint between the first position where the subject is closest and the second position where the subject is farthest is set as a center position of the sampling mark.   The analysis area control processing unit is configured to connect a first straight line connecting the first position and the second position, a second position connecting a center position of the sampling mark and the origin position where the ultrasonic beam is transmitted. The ultrasonic image diagnostic apparatus according to claim 4, wherein an angle formed with a straight line is set as the angle correction information.   The analysis area control processing unit sets a length in a longitudinal direction of the sampling mark based on the first and second positions, a center position of the sampling mark, and the angle correction information. Item 6. The ultrasonic diagnostic imaging apparatus according to Item 5. The analysis area control processing unit is configured to calculate the specific part motion speed V ′ obtained based on the Doppler signal output from the ultrasonic probe, the angle correction information θ, and the true motion speed V of the subject. connection of,
V = V '/ cosθ
The ultrasonic diagnostic imaging apparatus according to claim 6, wherein the length of the sampling mark in the longitudinal direction is set in consideration of the above. An analysis region setting program for setting the position and size of an analysis region for analyzing the function of a specific part of the subject on an ultrasound image of the subject generated by transmission / reception of an ultrasonic beam by an ultrasonic probe. And
On the computer,
A tracking function for tracking movement of a specific part in the subject from the ultrasound image;
Information on the position and size of the analysis region is set based on the tracking result, and information on angle correction of the analysis region with respect to the central axis of the ultrasonic beam for calculating the motion speed of the specific part is set. Setting function
An analysis region setting program characterized by realizing the above. The tracking function tracks the position of the annulus from the ultrasound image including a heart annulus as the subject,
The setting function sets a sampling mark formed in a rectangular shape as the analysis region on the ultrasound image, and a first position where the subject is closest to an origin position where the ultrasound beam is transmitted And the midpoint of the farthest second position as the center position of the sampling mark,
An angle formed by a first straight line connecting the first position and the second position and a second straight line connecting the center position of the sampling mark and the origin position where the ultrasonic beam is transmitted is Set as angle correction information,
Setting the length of the sampling mark in the longitudinal direction based on the first and second positions, the center position of the sampling mark, and the angle correction information;
9. The analysis area setting program according to claim 8, wherein the computer realizes the above.

Images