Huihua Kenny Chiang

Prof. Huihua Kenny Chiang

Institute of Biomedical Engineering, National Yang Ming Chiao Tung University, Taipei City 11221, Taiwan
orcid

Research Focus & Metrics

Academic Career Length

22Years (2005-Present)

(Calculated once Scilit publication data is synced: current year − year of first publication + 1)

Representative Publications

View All Publications
1 January 2011 in Journal of Biomedical Optics
Verified
Journal Article
Open access
PDF file
TSTai Yuan Su
Tai Yuan Su
CHChen Kerh Hwa
Chen Kerh Hwa
PLPo Hsuan Liu
Po Hsuan Liu
MWMing Hong Wu
Ming Hong Wu
DCDavid O. Chang
David O. Chang
Huihua Kenny ChiangHuihua Kenny Chiang
Huihua Kenny Chiang
Tear film provides lubrication and constant moisture to maintain vision and comfort of the eye. The outer layer of the tear film is a lipid layer preventing excessive tear evaporation and maintaining lubrication of the eye. A person may experience dry eye symptoms when an imbalance occurs in the tear film system or the tear film cannot remain continuous between blinks.1

1 December 2005 in Optical and Quantum Electronics
Verified
Journal Article
Closed access
PDF file
TWTsuei-Lian Wang
Tsuei-Lian Wang
Huihua Kenny ChiangHuihua Kenny Chiang
Huihua Kenny Chiang
HLHui-Hsin Lu
Hui-Hsin Lu
FPFang-Ying Peng
Fang-Ying Peng

22 September 2021 in Advanced Healthcare Materials
Verified
Journal Article
Closed access
PDF file
YLYi‐Xuan Li
Yi‐Xuan Li
SSShih‐Po Su
Shih‐Po Su
CYChou‐Hsun Yang
Chou‐Hsun Yang
MLMing‐Ho Liu
Ming‐Ho Liu
PLPin‐Ho Lo
Pin‐Ho Lo
YCYang-Hsiang Chan
Yang-Hsiang Chan

16 December 2010 in Journal of Raman Spectroscopy
Verified
Journal Article
Closed access
PDF file
PHPo‐Hsiang Hsu
Po‐Hsiang Hsu
Huihua Kenny ChiangHuihua Kenny Chiang
Huihua Kenny Chiang

16 July 2019 in BioMed Research International
Verified
Journal Article
Open access
WTWei-Nung Teng
Wei-Nung Teng
CTChien-Kun Ting
Chien-Kun Ting
YWYu-Tzu Wang
Yu-Tzu Wang
MHMing-Chih Hou
Ming-Chih Hou
WCWen-Kuei Chang
Wen-Kuei Chang
Chunli LinChunli Lin
Chunli Lin
During sedated endoscopic examinations, upper airway obstruction occurs. Nasal breathing often shifts to oral breathing during open mouth esophagogastroduodenoscopy (EGD). High-flow nasal cannula (HFNC) which delivers humidified 100% oxygen at 30 L min−1 may prevent hypoxemia. A mandibular advancement (MA) bite block with oxygen inlet directed to both mouth and nose may prevent airway obstruction during sedated EGD. The purpose of this study was to evaluate the efficacy of these airway devices versus standard management. One hundred and eighty-nine patients were assessed for eligibility. One hundred and fifty-three were enrolled. This study randomly assigned eligible patients to three arms: the standard bite block and standard nasal cannula, HFNC, and MA bite block groups. EGD was performed after anaesthetic induction. The primary endpoint was the oxygen desaturation area under curve at 90% (AUCDesat). The secondary endpoints were percentage of patients with hypoxic, upper airway obstruction, and apnoeic and rescue events. One hundred and fifty-three patients were enrolled. AUCdesat was significantly lower for HFNC and MA bite blocks versus the standard management (p= 0.019). The HFNC reduced hypoxic events by 18% despite similar airway obstruction and apnoeic events as standard group. The MA bite block reduced hypoxic events by 12% and airway obstructions by 32%. The HFNC and MA groups both showed a 16% and 14% reduction in the number of patients who received rescue intervention, respectively, compared to the standard group. The HFNC and MA bite block may both reduce degree and duration of hypoxemia. HFNC may decrease hypoxemic events while maintaining nasal patency is crucial during sedative EGD. The MA bite block may prevent airway obstruction and decrease the need for rescue intervention.

1 January 2008 in Journal of Biomedical Optics
Verified
Journal Article
Open access
PDF file
BRBin Rao
Bin Rao
LYLingfeng Yu
Lingfeng Yu
Huihua Kenny ChiangHuihua Kenny Chiang
Huihua Kenny Chiang
LZLeandro C. Zacharias
Leandro C. Zacharias
RKRonald Michael Kurtz
Ronald Michael Kurtz
ZCZhongping Chen
Zhongping Chen
A functional Fourier domain optical coherence tomography instrument offering spectral Doppler imaging of in vivo pulsatile human retinal blood flow was constructed. An improved phase-resolved algorithm was developed to correct bulk motion artifacts. Spectral Doppler imaging provides complementary temporal flow information to the spatially distributed flow information of the color Doppler image by providing direct visualization of the Doppler spectrum of the flow whose pattern can be further quantified with various velocity envelope curves and their corresponding flow indices. The coefficient of repeatability on resistance index measurement was assessed by analyzing 14 measurements on two vessels within two normal subjects. Previous studies suggest that retinal haemodynamics play an important role in glaucoma,1 diabetic retinopathy,2 and age-related macular degeneration.3 Therefore, it is desirable to assess retinal blood flow in vivo in order to investigate their role in major eye diseases. Although many flow measurement technologies have been developed, they do not provide adequate information for understanding the relationship between retinal blood flow and important ocular diseases.4 Doppler optical coherence tomography 5, 6, 7, 8, 9, 10, 11 (DOCT), a variation of optical coherence tomography (OCT) that combines coherence gating and laser Doppler effects, is an alternative that can overcome some limitations of the technologies mentioned above.4 Because of the unknown Doppler angle between the blood flow and incident light beam and the absence of a quantification method that can generate interpretable results for clinicians, DOCT is not widely used in eye clinics. Conventional DOCT only generates a snapshot of pulsatile ocular blood flow that is projected along the light beam direction in a cardiac cycle. Most recent developments in quantifying blood flow information of the human eye utilizes 3-D vascular orientation information to estimate Doppler angle and the absolute flow velocity.10, 11 On the other hand, it is noteworthy to quantify the pulsatile flow pattern as an alternative method to investigate retinal flow dynamics. A simple projected ocular blood flow velocity (integration over the whole blood vessel) plot through a cardiac cycle was chosen by White9 to demonstrate the pulsatile flow property, although it is possible to acquire much more hemodynamics information from the same raw data. The M-mode scanning method has been used to acquire temporal flow information in time-domain DOCT systems.12, 13 The short-time fast Fourier transformation method from Doppler ultrasound14 was used to generate Doppler spectrum wave forms13 but without further quantification, which provides the most valuable information for clinicians. The purpose of this paper is to implement the full concept of spectral Doppler imaging, developed by scientists and clinicians in ultrasound medicine, in a Fourier-domain DOCT system and provide an alternative quantification method for an ocular blood flow pattern that may be further investigated for vascular related eye diseases. “Spectral Doppler”14 is a terminology from Doppler ultrasound and should not be related to spectral OCT that uses spectral information as a contrast mechanism of OCT. Spectral Doppler imaging of pulsatile retinal blood flow includes Doppler spectrum visualization, using spectral Doppler wave forms, and a method for quantifying the temporal properties of flow, using various velocity envelopes and their corresponding Doppler-angle-independent indices.14 Continuous color Doppler data are acquired when the light beam performs repeated dense scans over the region of interest. Spectral Doppler analysis on the data shows how the velocity components and longitudinally projected flow-volume-rate change over time for scatters within the imaging volume with spectral Doppler wave forms. Various velocity envelope curves can be derived from spectral Doppler wave forms and used to extract the corresponding pulsatility index (PI), resistance index (RI) and several other indices that can provide interpretable Doppler-angle-independent information needed to quantify the pulsatile nature of ocular blood flow. A Fourier domain functional OCT system was developed for retinal blood flow imaging.8, 9 Briefly, low coherence light with a center wavelength of 890nm and FWHM bandwidth of 150nm was protected from optical feedback using a broadband optical isolator before entering a 2×2 broadband fiber coupler-based interferometer. The light from the reference arm was focused onto a reference mirror with an optical attenuator inserted into the optical path. The sample arm was modified from the patient module of a Zeiss Stratus OCT instrument. The detection arm was connected to a high-performance spectrometer that allows the system bench-top sensitivity of 100dB with 650μw light out of the sample arm fiber and 50μs CCD integration. A 9dB of SNR roll-off from zero imaging depth to 2mm imaging depth was observed. The system has an axial resolution of 3.5μm . In this study, the system speed was set at 16.7K A-lines/s with its CCD A-line integration time set at 50μs and the line period set at 60μs . The maximum longitudinal velocity in retinal tissue that corresponded to a phase difference of π was determined to be 2.69mm∕s according to Vmax=λ∕(4T) , where T is the line period of the CCD camera. The measured phase noise from a mirror was 2.66mrad .9 The velocity measurement error was determined to be <5% by imaging steady-state scattering flow (polystyrene bead solution with mean diameter of 0.3μm and volume concentration of 0.26%) pumped at different velocities. An improved phase-resolved algorithm that uses a Doppler variance image15 to select tissue pixels for histogram analysis16 of bulk phase was developed in this instrument to compensate for the axial eye movement. Most bulk motion presented in Fig. 1a has been corrected by the conventional histogram algorithm in Fig. 1b. Additional motion artifacts resulting from bulk motion phase estimation error by normal histogram algorithm shown in Fig. 1c were corrected in Fig. 1d by the improved phase-resolved algorithm that eliminates the assumption that the vessel has to be small compared to the tissues in the same A-line.16 <!-- named anchor --> Fig. 1 (a) Doppler image without bulk phase correction, (b) Doppler image with bulk phase correction by the conventional histogram method, (c) enlarged Doppler image of a selected area within the red window in (b), (d) doppler image generated by the improved phase-resolved algorithm with more reduced artifacts than normal histogram method in (c), (e) Spectral Doppler scan trace shown on 5×5mm OCT fundus image ( 140slices×512 A-lines), (f) the intensity and color Doppler image for one snapshot (700×256pixels) , (g) spectral Doppler wave forms that show the change of (axial) velocity and flow-volume-rate within a time span of ∼7.9s (the right grayscale bar is used to represent the volume-rate contribution for a given velocity bin) and (h) the maximum velocity envelope curve of eight cardiac cycles from which the quantitative Doppler-angle-independent flow indices were calculated. For the scan protocol of spectral Doppler imaging, repeated color Doppler scans over the broken region of the red line shown in Fig. 1e for a short period of time were performed after a specific vessel was selected from the OCT fundus17 image generated from a 3-D scan immediately before spectral Doppler imaging. Spectral Doppler analysis was performed on the color Doppler images. In order to evaluate variation of measurements within one session, spectral Doppler imaging of 512 snapshots through eight cardiac cycles in 7.9s was performed on one normal subject after the vessel was selected as shown in Fig. 1e. Each snapshot used 256 A-lines on 132.5μm tissue. The imaging speed of 65 snapshots per second is fast enough to capture the dynamic flow during a cardiac cycle. The structure image and velocity image for a typical snapshot are shown in Fig. 1f. The longitudinal velocity sensitivity was 174μm∕s by fitting the velocity profile and calculating the standard deviation for the velocity image across the center of the vessel8 in Fig. 1f. Before spectral Doppler analysis, a threshold was applied to the Doppler phase image and then some morphological operations were performed to get a vessel mask. The vessel center could then be easily determined from the peak positions after projecting the vessel mask to horizontal and vertical directions. A rectangular window centered at the estimated vessel center position can be applied to the velocity image in Fig. 1f. Next, the spectral Doppler wave forms in Fig. 1g could be generated after spectral Doppler analysis was performed for every velocity image snapshot that corresponded to one vertical line of the spectral Doppler wave forms. The velocity range was digitized into 256 velocity bins (v0−v255) . An iteration of all the pixels within the window generates a function n(vi) that represents the number of pixels having velocity of vi . Because each pixel represents a vortex for the vessel lumen area, the product of the velocity of the velocity bin and the number of pixels that fall within the velocity bin produced gray scale amplitude, which is proportional to the longitudinal projection of the flow-volume-rate contributed by each velocity bin. The summation of the gray scale amplitude along the y -axis gives the total longitudinal projection of flow-volume-rate of a vessel at a given time point. The summation of the gray scale amplitude along the y - and x -axes for one cardiac cycle provides the longitudinal projection of total flow volume within that cardiac cycle. Note that spectral Doppler wave forms with other gray scale intensity definitions can be defined according to other physical meanings, such as particle numbers or particle energies. Different color channels may be introduced to represent multiple flow-related parameters in a color-coded plot of spectral Doppler wave forms. In order to obtain more quantitative, interpretable results, different Doppler velocity envelopes, such as maximum velocity envelope, mean velocity envelope, and flow-volume-rate envelope, can be derived from the spectral Doppler wave forms accordingly. Figure 1h shows an example of the maximum velocity envelope derived from the spectral Doppler wave forms in Fig. 1g. For each cardiac cycle in Fig. 1h, symbol S represents the peak systolic maximum velocity, symbol D represents the end diastolic maximum velocity and symbol A represents the temporal average of maximum velocity. The PI and RI are defined as follows to characterize the curve and remove the dependence on Doppler angle: <!-- named anchor --> 1 PI=(S−D)A,RI=(S−D)S. Similar Doppler indices, such as the S∕D ratio and D∕A ratio, can be defined accordingly. Cycle 2 was excluded for statistical calculation due to eye motion. Table 1 summarizes the flow indices measured for other cardiac cycles, their average values, standard deviation values, and coefficients of variance for the maximum velocity envelope curve. The quantitative indices indicated above can be derived accordingly for other envelope curve definitions. <!-- named anchor --> Table 1 Intrasession coefficient of variation of Doppler flow indices for one measurement. Cycle PI RI SD DA 1 1.0753 0.62559 2.6709 0.64352 3 1.0551 0.621 2.6386 0.64393 4 1.1869 0.6805 3.1299 0.55727 5 1.1362 0.65217 2.875 0.606 6 1.4098 0.7561 4.1 0.45478 7 1.0667 0.61789 2.617 0.65965 8 1.0539 0.64 2.7778 0.59284 Mean 1.140557 0.656179 2.972743 0.593999 Std. Dev. 0.128572 0.049147 0.527933 0.070866 Coe. Var. 0.112727 0.0749 0.177591 0.119303 The intersession repeatability18 of RI measurement was assessed from multiple pairs of measurements. One measurement pair was defined as two separate measurements on the same vessel site. We measured two retinal vessels from two normal patients. Each vessel was measured seven times in the same day. There were C72=21 independent measurement pairs out of seven repeated measurements on each vessel. The intersession coefficient of repeatability18 (CoR) of 0.08336 was calculated from the 42 pairs of measurement on two vessels (see Table 2 ) according to the formula: CoR=1.96[∑(d0−d1)2]∕(n−1) . <!-- named anchor --> Table 2 Intersession coefficient of repeatability18 for 42 pairs of measurements on two vessels. Repeat Number 1 2 3 4 5 6 7 RI (Vessel 1) 0.7807 0.7326 0.7031 0.7968 0.7032 0.7805 0.7541 RI...