CN107280698A - A kind of mouse heart imaging system and method based on prediction formula ecg-gating - Google Patents

A kind of mouse heart imaging system and method based on prediction formula ecg-gating Download PDF

Info

Publication number
CN107280698A
CN107280698A CN201710367645.2A CN201710367645A CN107280698A CN 107280698 A CN107280698 A CN 107280698A CN 201710367645 A CN201710367645 A CN 201710367645A CN 107280698 A CN107280698 A CN 107280698A
Authority
CN
China
Prior art keywords
mouse
signal transacting
board analysis
detector
dsp signal
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201710367645.2A
Other languages
Chinese (zh)
Other versions
CN107280698B (en
Inventor
朱守平
樊宇
吕高旗
王民
曹旭
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Xidian University
Original Assignee
Xidian University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Xidian University filed Critical Xidian University
Priority to CN201710367645.2A priority Critical patent/CN107280698B/en
Publication of CN107280698A publication Critical patent/CN107280698A/en
Application granted granted Critical
Publication of CN107280698B publication Critical patent/CN107280698B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/54Control of apparatus or devices for radiation diagnosis
    • A61B6/541Control of apparatus or devices for radiation diagnosis involving acquisition triggered by a physiological signal
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/02Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
    • A61B6/03Computed tomography [CT]
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/50Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications
    • A61B6/503Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications for diagnosis of the heart
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/50Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications
    • A61B6/508Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications for non-human patients
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/52Devices using data or image processing specially adapted for radiation diagnosis

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Medical Informatics (AREA)
  • Pathology (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Veterinary Medicine (AREA)
  • Biophysics (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Public Health (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Optics & Photonics (AREA)
  • General Health & Medical Sciences (AREA)
  • Radiology & Medical Imaging (AREA)
  • Biomedical Technology (AREA)
  • Physics & Mathematics (AREA)
  • Molecular Biology (AREA)
  • Surgery (AREA)
  • Animal Behavior & Ethology (AREA)
  • Dentistry (AREA)
  • Oral & Maxillofacial Surgery (AREA)
  • Cardiology (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Physiology (AREA)
  • Apparatus For Radiation Diagnosis (AREA)

Abstract

The invention belongs to image processing technique field, a kind of mouse heart imaging system and method based on prediction formula ecg-gating are disclosed, including:Ecg signal acquiring front end, for obtaining mouse real-time heart beat information during CT IMAQs;Control card, triggering level is produced by obtaining mouse ecg information, and trigger the collection that X-ray detector completes mouse projected image;Micro CT acquisition systems, for being rebuild to 360 data for projection.The present invention is by catching intravital mouse ecg information, extract the R ripple crests of mouse electrocardiosignal, to judge the motion phase of mouse heart, and X-ray detector is triggered using Micro CT imaging techniques to mouse heart position progress data for projection sequence acquisition by R crests, multiple phase images of the mouse in the whole heart motion cycle under each angle are obtained, so as to recover the reconstruction image of multiple phases of the mouse heart in a complete heart beat cycle.

Description

A kind of mouse heart imaging system and method based on prediction formula ecg-gating
Technical field
The invention belongs to image processing technique field, more particularly to a kind of mouse heart based on prediction formula ecg-gating into As system and method.
Background technology
Micro-CT has imaging resolution high as the medical image system for toy compared to conventional medical CT The characteristics of, spatial resolution can reach micron order, while its reconstruction precision is also higher compared to conventional medical CT.Micro-CT Imaging is significant to cardiac imaging, and traditional Micro-CT imagings be present:Due to having used gate skill The timesharing that art carries out data for projection is mutually extracted, and traditional dynamic Micro-CT systems, which have triggering light source every time, can only gather one It is the problem of image, excessively slow this results in picking rate, considerably increase the light exposure of X-ray and reduce efficiency;During collection Between shortening mean the time for exposure reduce so that the contrast of image is relatively low, reduce the signal to noise ratio of reconstructed results.
In summary, the problem of prior art is present be:It is excessively slow to there is picking rate in traditional Micro-CT imagings, significantly Add the light exposure of X-ray and reduce efficiency;The contrast of image is relatively low, reduces the signal to noise ratio of reconstructed results.
The content of the invention
The problem of existing for prior art, the invention provides it is a kind of based on prediction formula ecg-gating mouse heart into As system and method.
The present invention is achieved in that a kind of mouse heart imaging system based on prediction formula ecg-gating, described to be based on The mouse heart imaging system of prediction formula ecg-gating includes:
Front end ecg signal acquiring module, for obtaining mouse real-time heart beat information during CT IMAQs;
DSP Signal transacting board analysis, triggering level is produced by obtaining the ecg information of mouse, and trigger X X-ray detection Xs Device completes the collection of mouse projected image;
Micro-CT medical image acquisition systems, are connected by RS232 with DSP Signal transacting board analysis, for to be scanned Object presses the multiple sections of 360 degree of scannings, for being rebuild to 360 data for projection.
Another object of the present invention is to provide a kind of mouse heart imaging system based on prediction formula ecg-gating Based on prediction formula ecg-gating mouse heart imaging method, it is described based on prediction formula ecg-gating mouse heart imaging side Method includes:Following steps:
Step one, mouse is placed on objective table, between X source and X-ray detector, required imaging region is placed in Within detector available field of view, the front end of bipolar lead is inserted in mouse left fore and right hind musculature, acupuncture needle end It is connected with button-shaped cardiac diagnosis lead joint;
Step 2, bipolar electrode lead access electrocardiogram acquisition module, is fixed on rear side of mouse, effective in X-ray detector Outside visual field, the analog output of electrocardio module accesses the AD conversion input interface of DSP Signal transacting board analysis;
Step 3, the AD conversion module of DSP Signal transacting board analysis receives the input of analog electrocardiogram signal, and records mouse electrocardio The maximum potential value of signal imitation input quantity, the most diastole moment beaten at the time of numerical value is produced as mouse heart, and Measure close to producing one at the time of this maximum by 0 to 1 Transistor-Transistor Logic level change, high level should be continued above 2ms, reach Required triggering the most short of X-ray detector is held time;
Step 4, DSP Signal transacting board analysis by AD conversion module judges mouse electrocardiosignal twice between crest Time, and judge the motion phase that a complete cycle planted agent of mouse heart beating obtains according to the time of integration needed for detector Position number;
Step 5, determines that configuration of the PC to X-ray detector is joined according to the monocycle phase number obtained by step 4 Number, the single frames time for exposure of configuration detector, image sequence totalframes, Binning parameters, detector triggering pattern, detector Scan pattern;
Step 6, the rotating disk for starting Micro-CT systems is enabled, and starts to catch mouse heartbeat signal, at mouse electrocardio R peaks Place sends trigger signal from DSP Signal transacting board analysis to X-ray detector, completes the projection data acquisitions under single angle;
Step 7, big one angle of disc spins is controlled after each angle acquisition is finished by computer, repeat afterwards into Operation in row step 6, until the data for projection of 360 angles is all gathered and finished.
Further, physical relationship is as follows in the step 4:
Monocycle phase number=heart beat cycle/individual projecting integral's time.
Further, X-ray detector is controlled to carry out IMAQ using hardware external trigger mode of operation, by DSP signals Handle board and produce short time high level, and input the Sync In interfaces of X-ray detector;X ray detectors are receiving From after the trigger request of DSP Signal transacting board analysis, IMAQ is proceeded by after one not fixed time period, and will adopt Collection result is put in caching;After an angle completes the projection image sequence collection of whole heart motion cycle, computer Data in caching are taken out, and carry out decoding preservation, and are numbered according to scanning angle and motion state, calculating is stored in In machine.
Further, after data acquisition is finished, according to the SMV document definitions preserved, the number that each angle is gathered According to being split, it is grouped by phase, i.e., the data for projection of same phase in 360 SMV files is proposed, be grouped, Carry out single-phase data reconstruction.
Further, the signal processing method of the DSP Signal transacting board analysis includes:
Step one, initial configuration is carried out to DSP Signal transacting board analysis;
Step 2, carries out the setting of the initial work and related work mode of analog-digital converter;By analog-digital converter Working frequency be set to standard working frequency 25MHz, then will enable analog-digital converter clock;
Step 3, completes system timer Timer0 configuration work;Elect Timer0 working frequency as 150MHz, The minimum time precision of Timer0 timers is 1 microsecond, and timer interruption time interval is set as into 1 millisecond, system meter secondary When device interrupt counter can store the count information of 32 bit lengths;
Step 4, carries out CPU and interrupts enable, CPU will receive the interrupt requests from system timer, DSP signal transactings Board will read serial ports reception area incessantly, determine whether the gate request signal from DSP Signal transacting board analysis;When DSP Signal transacting board analysis is to during the transmission gate request of DSP Signal transacting board analysis, and DSP Signal transacting board analysis end will believe for request Number respond, start global interrupt;
Step 5, DSP Signal transacting board analysis constantly carries out repeating reading to serial ports reception area, determines whether to come from The servomotor that DSP Signal transacting board analysis is sent enables switching operation signal;
Step 6, while waiting DSP Signal transacting board analysis to send servo enable operation signal by serial ports, DSP letters Number processing board global interrupt is in opening;When DSP Signal transacting board analysis programs are run to interrupt function, GPIO is read The value of pin, servomotor in off position, does not then jump out interrupt function, receives the request of DSP Signal transacting board analysis again;When When servomotor is closed, DSP Signal transacting board analysis proceeds by AD conversion, while Timer0 interrupt counters are certainly Increase;
Step 7, after DSP Signal transacting board analysis system resets, will carry out mouse ecg-r wave using initial 3 seconds The judgement of crest digital quantity, DSP Signal transacting board analysis is constantly AD converted;
Step 8, after interrupt function is again introduced into, DSP Signal transacting board analysis judges AD sampled value for R ripple crests When, GPIO01 pins will be driven high at once, continue 3 milliseconds, triggering X-ray detector carries out the collecting work of data for projection.
Further, the decoding of the detector acquired projections data of the DSP Signal transacting board analysis and storage method include:
The first step, the bus to computer is scanned for;
Second step, searches after X-ray detector, is realized by Dexela1512 types detector to X-ray detector Configuration;
3rd step, carries out the configuration of servomotor, calls admittedly high Control card to carry Dynamic Link Library Function and realizes opening Motion controller, sets the card number of current kinetic controller and reads the system clock of motion controller;To the current shape of motor State is judged;
4th step, the configuration of serial ports, computer is configured by RS232 serial ports to x-ray source, to the electricity of X radiographic sources Stream and voltage are configured;
5th step, computer is connected by serial ports with DSP Signal transacting board analysis control panels;When the imitative body setting in motion of dynamic When, computer enters wait state, and serial ports buffering area will be constantly examined, if having the letter from DSP Signal transacting board analysis Number;When DSP Signal transacting board analysis end completes to capture the period of motion of the imitative body of dynamic, the periodic quantity will be write to serial ports, and Read by serial ports of computers;
6th step, parameter configuration is carried out to X-ray detector by computer, and parameter includes projection sequence frame number, individual throwing Shadow time for exposure, detector triggering mode and high sensitivity/high dynamic mode of operation selection;
7th step, the configuration of data for projection frame number in the cycle is once triggered to X-ray detector to be needed to use by DSP signals The mouse electrocardiosignal cycle parameter of process plate card door control unit feedback;Get after mouse electrocardio periodic quantity, will be according to detection The time for exposure of individual data for projection of device carries out the division of motion phase to the imitative body of dynamic, and projection image sequence number is by mouse heart The mechanical movement cycle divided by individual data for projection time for exposure of X-ray detector obtain, projection image sequence number is mouse The motion phase number divided in one entire motion cycle of heartbeat;
8th step, carries out the selection that X-ray detector triggers pattern, and detector is set into external hardware triggering pattern, Triggering wait state is entered by computer control detector, then detector will be completed by DSP Signal transacting board analysis door control unit Triggering work;
9th step, detector will enter the wait state of trigger signal, and computer software will show current X-ray detector The number of image frames captured, and the fast refresh parameter;When detector receives the external trigger from DSP Signal transacting board analysis During signal, detector starts to carry out IMAQ according to the parameter of configuration, and the data for projection collected will pass through Camera Link High-speed transmission line writes Computer Cache;
Tenth step, after each detector is triggered and completes collecting work, checks whether and completes mouse heart 360 degrees of data are gathered, will be by computer software control GTS400 motion controls if not completing the collecting work of whole angles Card, the rotation of Micro-CT deep bids is realized through motor servo driver, and the anglec of rotation is 1 degree;When 360 degree of view data are adopted After collection is completed, EP (end of program).
Another object of the present invention is to provide it is a kind of be provided with the mouse heart based on prediction formula ecg-gating into As the Micro-CT systems of system.
Advantages of the present invention and good effect are:By catching intravital mouse ecg information, the R of mouse electrocardiosignal is extracted Ripple crest, judges the mechanical movement phase of mouse heart, and utilizes Micro-CT imaging techniques to mouse heart under phase Position carries out data for projection sequence acquisition, thus obtains mouse in the whole cardiac mechanical period of motion under each angle Multiple phase images, the pattern that Micro-CT is scanned using " turn-stopping " is operated, when 360 degree of projected images gather completion Afterwards, using it is angled under leggy image sequence carry out image reconstruction so that it is complete at one to recover mouse heart The reconstruction image of multiple phases in heart beat cycle.It make use of the projection image sequence acquisition technique of X-ray detector so that Whole mouse in the heart beat cycle can be disposably collected in a period of motion (i.e. one heart beat cycle) for mouse heart The projected image of motion phase, the dynamic imaging techniques that instead of former Micro-CT are only capable of collection one in a heart beat cycle The method of projected image is opened, system operating efficiency is improved.
Brief description of the drawings
Fig. 1 is the mouse heart imaging system structural representation provided in an embodiment of the present invention based on prediction formula ecg-gating Figure;
In figure:1st, front end ecg signal acquiring module;2nd, DSP Signal transacting board analysis;3rd, Micro-CT medical images are gathered System.
Fig. 2 is the flow chart of program of lower computer design provided in an embodiment of the present invention.
Fig. 3 is the design flow diagram of corresponding host computer acquisition software provided in an embodiment of the present invention.
Fig. 4 is the detector collecting work mode sequential signal of Micro-CT ecg-gatings provided in an embodiment of the present invention Figure.
Fig. 5 is the mouse heart imaging method flow chart provided in an embodiment of the present invention based on prediction formula ecg-gating.
Fig. 6 is the reconstructed results schematic diagram that the data provided in an embodiment of the present invention gathered using gating technology are carried out.
Fig. 7 is provided in an embodiment of the present invention not utilize the reconstructed results schematic diagram that gating technology is gathered.
Embodiment
In order to make the purpose , technical scheme and advantage of the present invention be clearer, with reference to embodiments, to the present invention It is further elaborated.It should be appreciated that the specific embodiments described herein are merely illustrative of the present invention, it is not used to Limit the present invention.
The application principle of the present invention is explained in detail below in conjunction with the accompanying drawings.
As shown in figure 1, the mouse heart imaging system provided in an embodiment of the present invention based on prediction formula ecg-gating includes:
Front end ecg signal acquiring module 1, for obtaining mouse real-time heart beat information during CT IMAQs;
DSP Signal transacting board analysis 2, triggering level is produced by obtaining the ecg information of mouse, and trigger the spy of X rays Survey the collection that device completes mouse projected image;
Micro-CT medical images acquisition system 3, is connected by RS232 with DSP Signal transacting board analysis 2, mainly by rotating disk, X-ray detector and X source are constituted, and whole Micro-CT will use drum-type structure, and it is similar to the CT clinically used, Now detector and x-ray source are then fixed in above a rotatable gantry, are rotated when scanning around testee One week, realize the IMAQ of multi-angle, and to be placed on holding position above the objective table at gantry center constant for sample.And under it Respectively there is a translation stage side, can be for regulation x-ray source and X-ray detector relative position.Toy is placed in rotating disk On the objective table of the heart, deep bid is at the uniform velocity rotated when collection image, realizes the scanning of multi-angle, and toy centre is fixed not It is dynamic.
Mouse heart imaging system provided in an embodiment of the present invention based on prediction formula ecg-gating is divided into host computer with Position machine two parts are realized:
Wherein slave computer (DSP Signal transacting board analysis) programming partial process view as shown in Fig. 2 specific method such as Under:
1. first, program will carry out initial configuration to DSP Signal transacting board analysis, its Program is from ROM to the complete of RAM Move entirely and the configuration work of GPIO pin will be partially completed at this.
2. after the basic initialization for completing DSP Signal transacting board analysis, ADC initial work and related work will be carried out Make the setting of mode.It is the setting of ADC work clock first, it is necessary to which ADC working frequency is set into standard work frequency Rate 25MHz, will then enable ADC clocks.ADC working method is set as sequential sampling, and scaling down processing is not done under 25MHz, Continuous sampling is set and covering sequence is opened, A0 sampling channels are chosen.
3. complete system timer Timer0 configuration work.Elect Timer0 working frequency as 150MHz, Timer0 The minimum time precision of timer is 1 microsecond, and timer interruption time interval is set as into 1 millisecond secondary, system timer Interrupt counter can store the count information of 32 bit lengths.
4. completing after ADC configuration work and the setting of system timer, carry out CPU and interrupt enable, now CPU will connect By the interrupt requests from system timer.Afterwards, DSP Signal transacting board analysis will read serial ports reception area incessantly, judge Whether gate request signal from host computer is had.When host computer sends gate request to DSP Signal transacting board analysis, DSP letters Number processing board end will be responded for request signal, start global interrupt.
5.DSP Signal transacting board analysis constantly carries out repeating reading to serial ports reception area, determines whether to come from host computer The servomotor of transmission enables switching operation signal.If DSP Signal transacting board analysis authenticates to the private clothes electricity that host computer is sent Machine enables open signal, can be by electromagnetism relay by the pin by controlling GPIO02 to electromagnetic relay progress conducting control Device regards a switch as, and when the GPIO02 pins of DSP Signal transacting board analysis export high level, electromagnetic relay is in and led Lead to state, now motor servo driver will be enabled, and computer can be moved by control card to Micro-CT rotating disks Control.After computer control deep bid rotates an angle, send servo to DSP Signal transacting board analysis from computer and enable closing Signal, DSP Signal transacting board analysis upon receipt of the signal drags down GPIO02 pins, so as to realize the closing that servo is enabled Operation.
6. while waiting host computer to send servo enable operation signal by serial ports, DSP Signal transacting board analysis is global Interrupt and be in opening, now DSP Signal transacting board analysis can be interrupted normal through Timer0 enters interrupt function.Work as DSP When Signal transacting board analysis program is run to interrupt function, the value of GPIO pin is read first, it is ensured that now servomotor, which is in, closes Closed state, if servomotor not in off position, jumps out interrupt function, receives the request of host computer again.Work as servomotor When being closed, DSP Signal transacting board analysis proceeds by AD conversion, while Timer0 interrupt counters increase certainly.
7. after DSP Signal transacting board analysis system resets, mouse ecg-r wave crest will be carried out using initial 3 seconds The judgement of digital quantity, now DSP Signal transacting board analysis will be constantly AD converted, and be occurred when AD changes larger numerical value When, the smaller value occurred before covering, and maximum digital quantity is saved be designated as Max.Then, a memory space is utilized The successively ecg signal data that collects is preserved for 5 array A [5] circulations, only when array element is all higher than 0.95Max, and Meet the 3rd elements A [2] of array and be in array element maximum, A [0], A [1] are incremented by, when A [3], A [4] successively decrease, during by this Carve the approximate most diastole moment for regarding mouse heart mechanical movement as., can be according between adjacent peaks according to the judgement of R ripple crests Timer0 interrupts counting and is worth to mouse core hop cycle, and unit is millisecond.DSP Signal transacting board analysis calculates the heartbeat for obtaining mouse After cycle, the information is fed back into computer by RS232 serial ports.
8. after interrupt function is again introduced into, once DSP Signal transacting board analysis judges AD sampled value for R ripple crests When, GPIO01 pins will be driven high at once, continue 3 milliseconds, so as to trigger the collection work that X-ray detector carries out data for projection Make.After being successfully triggered to X-ray detector, DSP Signal transacting board analysis will not wait detector collection to finish, and enter at once In next system break stage, the gathered data of X-ray detector is obtained and storage operation will be completed by computer.
Upper computer software (the meter of mouse heart imaging system provided in an embodiment of the present invention based on prediction formula ecg-gating Calculation machine acquisition software) cooperation DSP Signal transacting board analysis carry out X-ray detector relevant exposure parameter configuration, and detection The setting of device triggering mode, and complete the decoding and storage work of detector acquired projections data, specific Software for Design flow As shown in figure 3, method is as follows:
1. program starts, system scans for the bus to computer, to ensure that X-ray detector is located at bus On, if not finding detector device in computer bus, eject the prompting frame for not finding equipment, and terminator.
2. after X-ray detector is successfully searched, system be able to can now pass through by hanging equipment in bus The function that the software development kit of Dexela1512 type detectors is provided realizes the configuration to X-ray detector, now required to match somebody with somebody The parameter detector put is also imperfect, temporarily hangs up detector and enters state to be configured in bus.
3. following program will carry out the configuration of servomotor.Before program operation, Killmorgen servo electricity need to be ensured Machine driver is in "on" position, and motor enables switch and need to ensure to have opened, and driver passes through admittedly high GTS400 motion controls Board is connected with computer.On the basis of ensureing that hardware setting and connection are correct, carry out the initialization to servomotor and grasp Make, admittedly high Control card will be called to carry Dynamic Link Library Function and realize open movement controller, current kinetic control is then set The card number of device processed and the system clock for reading motion controller.Initialization procedure also needs to sentence the current state of motor It is disconnected, if be kept in motion under motor present case, allow its stop motion.After the completion of initialization operation, program will be returned One Boolean, when returning to TRUE, expression is initialized successfully.
4. following software will realize the configuration of serial ports, PC will pass through serial ports and x-ray source and DSP Signal transacting board analysis Board is communicated.Computer is configured by RS232 serial ports to x-ray source, it is necessary to the electric current and voltage of x-ray source Configured, tube voltage is generally set to 50kV, tube current is set to 300uA, and focal spot size is set to middle focal spot.Meanwhile, The current real-time voltage current parameters of X-ray are fed back to computer by x-ray source information monitoring module, facilitate user to be seen Examine.
5. simultaneously, computer is connected by serial ports with DSP Signal transacting board analysis control panels.The baud rate of serial communication will be set It is set to 9600, with 8 data bit, 1 stop position, no parity.It should be noted that servomotor is transported in upper electricity When dynamic, stronger magnetic interference can be produced, 3 common line RS232 Serial Port Lines will be affected by magnetic fields, so as to influence To the normal serial communication between computer and DSP Signal transacting board analysis.Therefore, be RS422 interfaces by RS232 interface conversions, And by way of external 5V voltage-stabilized power supplies, improve the driving force and antijamming capability of Serial Port Transmission circuit.Due to RS232 Interface conversion with RS422 is entirely that hardware is completed, it is not necessary to corresponding computer software and DSP Signal transacting board analysis journeys Sequence is modified.When dynamically imitative body setting in motion, computer enters wait state, and serial ports buffering area will constantly be examined, and see Whether signal from DSP Signal transacting board analysis is had.When DSP Signal transacting board analysis end completes to catch the period of motion of the imitative body of dynamic When obtaining, the periodic quantity will be write to serial ports, and read by serial ports of computers.
6. next needing to carry out X-ray detector parameter configuration by computer, parameter includes projection sequence frame number, list Open projection exposure time, detector triggering mode and high sensitivity/high dynamic mode of operation selection.
7. the configuration that pair X-ray detector once triggers data for projection frame number in the cycle needs to use by DSP signal transactings The mouse electrocardiosignal cycle parameter of board door control unit feedback.After program gets mouse electrocardio periodic quantity, program is by root The division of motion phase is carried out to the imitative body of dynamic according to the time for exposure of detector individual data for projection, projection image sequence number is by small The mechanical movement cycle of rat heart divided by individual data for projection time for exposure of X-ray detector obtain, and projection image sequence number is The motion phase number divided in the entire motion cycle beaten by mouse heart.
8. needing to carry out the selection that X-ray detector triggers pattern, detector is set to external hardware triggering pattern, by Computer control detector enters triggering wait state, then will complete detector by DSP Signal transacting board analysis door control unit Trigger work.
9. when detector, which is configured, to be completed, detector will enter the wait state of trigger signal, now computer software will Show the number of image frames that current X-ray detector has been captured, and the fast refresh parameter.When detector is received from DSP letters Number processing board outer triggering signal when, detector start according to configuration parameter carry out IMAQ, the projection collected Data will write Computer Cache by Camera Link high-speed transmission lines.When all images two in a triggering cycle enter Data processed, which are all gathered, to be completed and writes after Computer Cache, completes the decoding operate to buffer area initial data by software, and press Decoded raw data are saved in computer specified location by parameter according to detector setup parameter.
10. after each detector is triggered and completes collecting work, program, which will be checked whether, completes mouse heart 360 degrees of data collection, will be by computer software control GTS400 motion controls if not completing the collecting work of whole angles Card, the rotation of Micro-CT deep bids is realized through motor servo driver, and the anglec of rotation is 1 degree.When 360 degree of view data are adopted After collection is completed, EP (end of program).
The hardware type external trigger working method of present invention application X-ray detector, i.e., be responsible for X-ray detector by PC Parameter configuration, but order detector carry out actual exposure operation trigger signal then transfer to third party device to complete, by based on The Control cards of DSP Signal transacting board analysis chips completes the generation of external trigger signal.
As shown in figure 5, the mouse heart imaging method provided in an embodiment of the present invention based on prediction formula ecg-gating includes Following steps:
Step one, mouse is placed on objective table, between X source and X-ray detector, required imaging region is placed in Within detector available field of view, the front end (acupuncture needle) of bipolar lead is inserted in mouse left fore and right hind musculature, silver Pin end is connected with button-shaped cardiac diagnosis lead joint, and contact area is increased with conducting resinl;
Step 2, bipolar electrode lead access electrocardiogram acquisition module, is fixed on rear side of mouse, effective in X-ray detector Outside visual field, the analog output of electrocardio module accesses the AD conversion input interface of DSP Signal transacting board analysis;
Step 3, the AD conversion module of DSP Signal transacting board analysis receives the input of analog electrocardiogram signal, and records mouse electrocardio The maximum potential value of signal imitation input quantity, the most diastole moment beaten at the time of this numerical value is produced as mouse heart, and Produce one at the time of measuring close to this maximum by 0 to 1 TTL level to be changed, high level should be continued above 2ms, with Triggering the most short of X-ray detector is held time needed for reaching;
Step 4, DSP Signal transacting board analysis by AD conversion module judges mouse electrocardiosignal twice between crest Time, and judge the motion phase that a complete cycle planted agent of mouse heart beating obtains according to the time of integration needed for detector Position number, physical relationship is as follows:
Monocycle phase number=heart beat cycle/individual projecting integral's time;
Step 5, determines that configuration of the PC to X-ray detector is joined according to the monocycle phase number obtained by step 4 Number, the single frames time for exposure of configuration detector, image sequence totalframes, Binning parameters, detector triggering pattern, detector Scan pattern;
Step 6, the rotating disk for starting Micro-CT systems is enabled, and starts to catch mouse heartbeat signal, at mouse electrocardio R peaks Place sends trigger signal from DSP Signal transacting board analysis to X-ray detector, completes the projection data acquisitions under single angle;
Step 7, big one angle of disc spins is controlled after each angle acquisition is finished by computer, repeat afterwards into Operation in row step 6, until the data for projection of 360 angles is all gathered and finished.
In the step 5, it is contemplated that mouse heart jumping frequency rate is very fast, in order to be able to which mouse heart exists in accurate distinguish Significant difference under different heartbeat phase states, the time for exposure of detector is unsuitable long, it is proposed that be chosen at current Binning Minimum exposure time under parameter carries out data for projection scanning, Binning parameter settings and detector minimum exposure time;
The sequence totalframes of detector should be consistent with gained monocycle phase number in step 4, but is considered as when rebuilding The significant difference of different motion phasetophase simultaneously shortens the time needed for rebuilding as far as possible, can not exclusively rebuild all motion phases Position, but be separated by a fixed numbers and rebuild, reduce the motion state number finally rebuild and obtained;
For X-ray detector, the software triggering pattern of general Micro-CT systems is different from, is touched using outside hardware Send out mode of operation to control X-ray detector to carry out IMAQ, short time high level is produced by DSP Signal transacting board analysis, and Input the Sync In interfaces of X-ray detector;X-ray detector is receiving the trigger request from DSP Signal transacting board analysis Afterwards, by a not fixed time period (when the time length is not gathered finish remaining by X-ray detector previous frame and scanned Between determine) after proceed by IMAQ, and collection result is put in caching;Whole heart movement is completed in an angle After the projection image sequence collection in cycle, computer takes out the data in caching, and carries out decoding preservation, and according to scanning Angle and motion state are numbered, and are stored in computer.
After data acquisition is finished, according to the SMV document definitions preserved, the data that each angle is gathered are torn open Point, it is grouped by phase, i.e., the data for projection of same phase in 360 SMV files is proposed, be grouped, then carried out Single-phase data reconstruction.
The application effect of the present invention is explained in detail with reference to contrast experiment.
During system performance testing, the dynamic for carrying out Micro-CT using an imitative body of the motion rotated at a high speed is adopted Collection, the movement velocity for imitating body is close with mouse heart, is 500 complete cycles of motion per minute.Fig. 6 is to utilize this patent The reconstructed results that data that gating technology is gathered are carried out, Fig. 7 is the reconstructed results that are gathered without gating technology.It can see Arrive, certain flashy phantom image in the imitative body and successful reconstitution of high-speed motion can be captured using gating technology, and it is unfavorable It is then disorderly and unsystematic with the reconstructed results of gating technology.Because mouse heart motion amplitude is more much smaller than imitative body, therefore the present invention Prediction formula gating technology can be very good adapt to mouse palmic rate and amplitude, obtain preferable collection result.
The foregoing is merely illustrative of the preferred embodiments of the present invention, is not intended to limit the invention, all essences in the present invention Any modifications, equivalent substitutions and improvements made within refreshing and principle etc., should be included in the scope of the protection.

Claims (8)

1. a kind of mouse heart imaging system based on prediction formula ecg-gating, it is characterised in that described to be based on prediction formula electrocardio The mouse heart imaging system of gate includes:
Front end ecg signal acquiring module, for obtaining mouse real-time heart beat information during CT IMAQs;
DSP Signal transacting board analysis, triggering level is produced by obtaining the ecg information of mouse, and it is complete to trigger X-ray detector Into the collection of mouse projected image;
Micro-CT medical image acquisition systems, are connected by RS232 with DSP Signal transacting board analysis, for object to be scanned By 360 degree of scannings data for projection, rebuild for the data for projection to 360 degree.
2. it is a kind of as claimed in claim 1 based on prediction formula ecg-gating mouse heart imaging system based on prediction formula electrocardio The mouse heart imaging method of gate, it is characterised in that the mouse heart imaging method bag based on prediction formula ecg-gating Include:Following steps:
Step one, mouse is placed on objective table, between x-ray source and X-ray detector, required imaging region is placed in X Within ray detector available field of view, the front end of bipolar lead is inserted in mouse left fore and right hind musculature, acupuncture needle End is connected with button-shaped cardiac diagnosis lead joint;
Step 2, bipolar electrode lead access electrocardiogram acquisition module, is fixed on rear side of mouse, in X-ray detector available field of view Outside, the analog output of electrocardio module accesses the AD conversion input interface of DSP Signal transacting board analysis;
Step 3, the AD conversion module of DSP Signal transacting board analysis receives the input of analog electrocardiogram signal, and records mouse electrocardiosignal The maximum potential value (mouse electrocardio R peaks) of simulation input amount, the most diastole beaten at the time of numerical value is produced as mouse heart Moment, and one is produced by 0 to 1 Transistor-Transistor Logic level change at the time of measuring close to this maximum, high level should be continued above 2ms, triggering the most short of X-ray detector is held time needed for reaching;
Step 4, DSP Signal transacting board analysis by AD conversion module judges time of the mouse electrocardiosignal twice between crest, And judge the motion phase that a complete cycle planted agent of mouse heart beating obtains according to the time of integration needed for X-ray detector Position number;
Step 5, determines configuration parameter of the PC to X-ray detector according to the monocycle phase number obtained by step 4, matches somebody with somebody Put single frames time for exposure of detector, image sequence totalframes, Binning parameters, detector triggering pattern, detector scanning mould Formula;
Step 6, the rotating disk for starting x-ray source and Micro-CT systems is enabled, and starts to catch mouse heartbeat signal, in mouse core At electric R peaks trigger signal, the projected image of whole heart motion cycle are sent from DSP Signal transacting board analysis to X-ray detector Sequence acquisition.
Step 7, big one angle of disc spins is controlled after each angle acquisition is finished by computer, repeats step afterwards Operation in rapid six, until 360 degree of data for projection is all gathered and finished.
3. the mouse heart imaging method as claimed in claim 2 based on prediction formula ecg-gating, it is characterised in that the step Physical relationship is as follows in rapid four:
Monocycle phase number=heart beat cycle/individual projecting integral's time.
4. the mouse heart imaging method as claimed in claim 2 based on prediction formula ecg-gating, it is characterised in that using hard Part external trigger mode of operation controls X-ray detector to carry out IMAQ, by DSP Signal transacting board analysis produces short time high electricity It is flat, and input the Sync In interfaces of X-ray detector;X-ray detector is receiving the triggering from DSP Signal transacting board analysis After request, IMAQ is proceeded by after one not fixed time period, and collection result is put in caching;At one Angle is completed after the projection image sequence collection of whole heart motion cycle, and computer takes out the data in caching, goes forward side by side Row decoding is preserved, and is numbered according to scanning angle and motion state, is stored in computer.
5. the mouse heart imaging method as claimed in claim 2 based on prediction formula ecg-gating, it is characterised in that data are adopted After collection is finished, each angle acquired image sequence is split, is grouped by phase, i.e., by 360 degree of projection The projected image of same phase is combined in data, participates in rebuilding.
6. the mouse heart imaging method as claimed in claim 2 based on prediction formula ecg-gating, it is characterised in that described The signal processing method of DSP Signal transacting board analysis includes:
Step one, initial configuration is carried out to DSP Signal transacting board analysis;
Step 2, carries out the setting of ADC initial work and related work mode;ADC working frequency is set to mark Quasi- working frequency 25MHz, will then enable ADC clocks;
Step 3, completes system timer Timer0 configuration work;Elect Timer0 working frequency as 150MHz, Timer0 The minimum time precision of timer is 1 microsecond, and timer interruption time interval is set as into 1 millisecond secondary, system timer Interrupt counter can store the count information of 32 bit lengths;
Step 4, carries out CPU and interrupts enable, CPU will receive the interrupt requests from system timer, DSP Signal transacting board analysis Serial ports reception area will be read incessantly, determine whether the gate request signal from DSP Signal transacting board analysis;When DSP letters Number processing board to DSP Signal transacting board analysis send gate request when, DSP Signal transacting board analysis end will be made for request signal Response, starts global interrupt;
Step 5, DSP Signal transacting board analysis constantly carries out repeating reading to serial ports reception area, determines whether from DSP letters Number processing board send servomotor enable switching operation signal;
Step 6, while waiting DSP Signal transacting board analysis to send servo enable operation signal by serial ports, at DSP signals Manage board global interrupt and be in opening;When DSP Signal transacting board analysis programs are run to interrupt function, GPIO pin is read Value, servomotor in off position, then do not jump out interrupt function, receives the request of DSP Signal transacting board analysis again;Work as servo When motor is closed, DSP Signal transacting board analysis proceeds by AD conversion, while Timer0 interrupt counters increase certainly;
Step 7, after DSP Signal transacting board analysis system resets, will carry out mouse ecg-r wave crest using initial 3 seconds The judgement of digital quantity, DSP Signal transacting board analysis is constantly AD converted;
Step 8, after interrupt function is again introduced into, when DSP Signal transacting board analysis judges AD sampled value for R ripple crests, GPIO01 pins will be driven high at once, continue 3 milliseconds, triggering X-ray detector carries out the collecting work of data for projection.
7. the mouse heart imaging method as claimed in claim 2 based on prediction formula ecg-gating, it is characterised in that described The decoding of the detector acquired projections data of DSP Signal transacting board analysis and storage method include:
The first step, the bus to computer is scanned for;
Second step, searches after X-ray detector, parameter configuration is carried out to X-ray detector;
3rd step, carries out the configuration of servomotor, and open movement controller sets the card number of current kinetic controller and reads fortune The system clock of movement controller;Current state to servomotor is judged;
4th step, the configuration of x-ray source, computer is configured by RS232 serial ports to x-ray source, to the electric current of x-ray source Configured with voltage;
5th step, computer is connected by serial ports with DSP Signal transacting board analysis control panels;Computer enters wait state, serial ports Buffering area will be constantly examined, if having the signal from DSP Signal transacting board analysis;When the completion pair of DSP Signal transacting board analysis end When the mouse core coron cycle captures, the periodic quantity will be write to serial ports, and read by serial ports of computers;
6th step, parameter configuration is carried out to X-ray detector by computer, and parameter includes projection sequence frame number, individual projection and exposed Between light time, detector triggering mode and high sensitivity/high dynamic mode of operation selection;
7th step, the configuration of data for projection frame number in the cycle is once triggered to X-ray detector to be needed to use by DSP signal transactings The mouse electrocardiosignal cycle parameter of board door control unit feedback;Get after mouse electrocardio periodic quantity, will be according to detector list Division of the time for exposure of data for projection to mouse heart motion phase is opened, projection image sequence number is by the mouse heart period of motion Divided by individual data for projection time for exposure of X-ray detector obtains, projection image sequence number is one of mouse heart beating The motion phase number divided in the entire motion cycle;
8th step, carries out the selection that X-ray detector triggers pattern, detector is set into external hardware triggering pattern, by counting Calculation machine control detector enters triggering wait state, then will complete touching for detector by DSP Signal transacting board analysis door control unit Send out work;
9th step, detector will enter the wait state of trigger signal, and computer software will show that current X-ray detector has been caught The number of image frames received, and the fast refresh parameter;When detector receives the outer triggering signal from DSP Signal transacting board analysis When, detector starts to carry out IMAQ according to the parameter of configuration, and the data for projection collected will be by Camera Link at a high speed Transmission line writes Computer Cache;
Tenth step, after each detector is triggered and completes collecting work, checks whether complete mouse heart 360 degree Data acquisition, if not completing the collecting work of whole angles, will realize Micro-CT by computer software control servomotor The rotation of deep bid, the anglec of rotation is 1 degree;After 360 degree of view data gather completion, EP (end of program).
8. a kind of Micro-CT systems for being provided with the mouse heart imaging system based on prediction formula ecg-gating described in claim 1 System.
CN201710367645.2A 2017-05-23 2017-05-23 Prospective electrocardio-gating-based mouse heart imaging system and method Active CN107280698B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710367645.2A CN107280698B (en) 2017-05-23 2017-05-23 Prospective electrocardio-gating-based mouse heart imaging system and method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710367645.2A CN107280698B (en) 2017-05-23 2017-05-23 Prospective electrocardio-gating-based mouse heart imaging system and method

Publications (2)

Publication Number Publication Date
CN107280698A true CN107280698A (en) 2017-10-24
CN107280698B CN107280698B (en) 2020-12-08

Family

ID=60094611

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710367645.2A Active CN107280698B (en) 2017-05-23 2017-05-23 Prospective electrocardio-gating-based mouse heart imaging system and method

Country Status (1)

Country Link
CN (1) CN107280698B (en)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109222964A (en) * 2018-07-20 2019-01-18 广州视源电子科技股份有限公司 atrial fibrillation detection device and storage medium
US11557037B2 (en) 2019-08-20 2023-01-17 Shanghai United Imaging Healthcare Co., Ltd. Systems and methods for scanning data processing
CN117115577A (en) * 2023-10-23 2023-11-24 南京安科医疗科技有限公司 Cardiac CT projection domain optimal phase identification method, equipment and medium

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20090207968A1 (en) * 2006-06-22 2009-08-20 Koninklijke Philips Electronics N. V. Dual x-ray tube gating
CN105142503A (en) * 2013-03-20 2015-12-09 皇家飞利浦有限公司 Neurophysiological monitoring for prospective motion gating in radiological imaging

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20090207968A1 (en) * 2006-06-22 2009-08-20 Koninklijke Philips Electronics N. V. Dual x-ray tube gating
CN105142503A (en) * 2013-03-20 2015-12-09 皇家飞利浦有限公司 Neurophysiological monitoring for prospective motion gating in radiological imaging

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
NANCY L.FORD ET AL.: "Prospective respiratory-gated micro-CT of free breathing rodents", 《MED.PHUS.》 *
XIAOLIAN GUO ET AL.: "4D Micro-CT using Fast Prospective Gating", 《PHYS MED BIOL》 *

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109222964A (en) * 2018-07-20 2019-01-18 广州视源电子科技股份有限公司 atrial fibrillation detection device and storage medium
US11557037B2 (en) 2019-08-20 2023-01-17 Shanghai United Imaging Healthcare Co., Ltd. Systems and methods for scanning data processing
US11972565B2 (en) 2019-08-20 2024-04-30 Shanghai United Imaging Healthcare Co., Ltd. Systems and methods for scanning data processing
CN117115577A (en) * 2023-10-23 2023-11-24 南京安科医疗科技有限公司 Cardiac CT projection domain optimal phase identification method, equipment and medium
CN117115577B (en) * 2023-10-23 2023-12-26 南京安科医疗科技有限公司 Cardiac CT projection domain optimal phase identification method, equipment and medium

Also Published As

Publication number Publication date
CN107280698B (en) 2020-12-08

Similar Documents

Publication Publication Date Title
US7020234B2 (en) Method for producing tomograms of a periodically moving object with the aid of a focus detector combination
CN101283912B (en) X-ray ct apparatus
RU2654611C2 (en) Automatic positioning of standard planes for real-time fetal heart evaluation
US8777856B2 (en) Diagnostic system and method for obtaining an ultrasound image frame
CN107280698A (en) A kind of mouse heart imaging system and method based on prediction formula ecg-gating
JP5764058B2 (en) Ultrasound assessment of cardiac synchrony and viability
CN101049243B (en) X-ray ct device
US7052460B2 (en) System for producing an ultrasound image using line-based image reconstruction
CN102469984A (en) Three dimensional fetal heart imaging by non-ecg physiological gated acquisition
US20050111622A1 (en) Method for production of tomographic section images of a periodically moving object with a number of focus detector combinations
CN108042125B (en) High-speed endoscopic optical coherent blood flow imaging system
CN105592799A (en) Ultrasound systems and methods for automated fetal heartbeat identification
WO2011059584A1 (en) Apparatus and methods for computed tomography imaging
CN101541232A (en) Reconstruction window adaption in ECG-gated computed tomography
CN102462508B (en) System and method for ultrasound imaging
CN101331521A (en) Iterative image reconstruction of a moving object from projection data
JP2024512663A (en) System and method of image enhancement for in-motion multipulse-actuated X-ray source tomosynthesis devices using electrocardiogram gating
CN103247071B (en) A kind of structure three-dimensional blood vessel model method and apparatus
CN103635145A (en) Ultrasound diagnostic device, ultrasound image processing device, and medical image diagnostic device
US20050175141A1 (en) Method for producing tomograms of a periodically moving object with the aid of a focus/detector combination
CN101352352B (en) X-ray ct device
CN106859683A (en) Heart dynamic analog device
JP4071101B2 (en) Ultrasonic diagnostic equipment
JP2001161680A (en) Method and apparatus for controlling quantity of x-ray irradiation during synchronized heart scanning
CN102258383B (en) X ray computer tomography device and image processing method

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant