1 Introduction

Analyzer-based imaging(ABI)is an important X-ray phase-contrast imaging(XPCI)method that has higher contrast and resolution than conventional X-ray absorptioncontrast imaging(XACI),particularly for materials made of light elements.Therefore,the ABI method has great potential in medical diagnosis[1–5],biology[6–9],material science[10,11],and others.Thanks to progress in synchrotron radiation(SR)sources,ABI has been developed and applied extensively.In parallel with the advances in synchrotron-based ABI methods,there has been a continualdevelopmentoflaboratory source technology[12–15]that aims to translate the ABI method into practical applications.ABI has the higher signal-to-noise ratio(SNR)[21]than other well-known XPCI methods,such as propagation-based imaging(PBI)[16,17]and grating interferometry(GI)[18–20],especially when high X-ray energies are employed.This can effectively reduce the absorbed dose to the objects and present a big advantage for those applications in which thick or high-Z materials have to be investigated(such as imaging of whole and large human organs)[13,21].Furthermore,the effects of beam hardening and scattering in ABI can be negligible due to the monochromaticity of the X-rays and the angular selectivity of the analyzer.

Several types of information,i.e.,absorption,refraction,and ultra-small-angle X-ray scattering(USAXS),can be retrieved by ABI using appropriate extraction methods,including diffraction-enhanced imaging (DEI) [22],extended diffraction-enhanced imaging(E-DEI)[23,24],generalized diffraction-enhanced imaging(G-DEI)[25],multiple-image radiography(MIR)[26–28],and Gaussian curve fitting(GCF)[29].Two images are acquired in DEI and E-DEI to extract absorption and refraction information.Image contrast is poor in samples with a subpixel structure—namely with a strong USAXS signal,such as the lung.This is because no USAXS is considered[30].Further,G-DEI,MIR,and GCF can extract absorption,refraction,and USAXS simultaneously.Only three images are collected in G-DEI;however,the extraction results are limited by the width of the rocking curve,which is greatly affected by the noise[30,31].MIR and GCF are both statistical methods that acquire a set of images by rocking the analyzer;thus,the extraction results are more accurate and have better anti-noise performance.As Gaussian fitting is performed for each pixel,the computational time of the GCF increases signi ficantly[29,30].In general,MIR is a better information extraction method,especially when the radiation dose and acquisition time are not critical.The MIR method was separately proposed by Oltulu et al.[26],Wernick et al.[28],and Pagot et al.[27]in 2003,and it has been extensively studied due to its robustness,scattering rejection,and high extraction accuracy.Previously,research and improvements have been carried out on the characteristics of the MIR method.Wernick et al.[32,33]evaluated the effects of the number of acquired images on the accuracy of extraction results using the Crame´r–Rao lower bound(CRLB).Hu et al.[34]extracted three kinds of information using MIR and fused the contrast images in one image,while a comparison of refraction information extraction methods in analyzer-based imaging was carried out[24].Huang et al.[35,36]and Zhu and his coworkers[37,38]discussed the limitations of MIR and proposed some improved information extraction methods.Diemoz et al.[21,30]compared the extraction results of several ABI information extraction methods in different conditions and theoretically compared the SNR and FOM of three XPCI techniques of ABI,PBI,and GI.Li et al.[31,39]improved the MIR method to extract more accurate refraction information and compared the refraction angle extraction results of the DEI,G-DEI,and MIR methods with different sampling strategies.

However,MIR has been the subject of extensive research in recent years,particularly when using SR sources.However,many obstacles remain to use conventional X-ray sources.Poisson noise has a signi ficant effect on extraction results due to the low intensity of the source[32,33].Here,the selection of the sampling strategies becomes more complicated.In addition,previous reports have focused on refraction information[24,31].In fact,absorption,refraction,and USAXS dominate these different cases.Absorption is dominant for samples with high absorption coefficients and/or large dimensions—USAXS is prominent for samples with a subpixel structure,such as lung,powders,and paper.

In this paper,the effects of Poisson noise on extraction results were studied to guide the application of MIR in lowintensity sources(such as X-ray tube sources),so that an optimized radiation doseand imaging time can be employed.The effects of different sampling strategies(including range and interval)on the extraction results are studied while considering the Poisson noise.This provides references for the selection of sampling strategies in practical applications.Furthermore,the mutual effects of the three kinds of information in MIR were analyzed,and its manifestations,causes,and trends are given below.These can be used to optimize information extraction and lay the foundation for correcting MIR extraction results.

现阶段，多数的上市公司为了提升自身利益在信息披露过程中有选择性，没有根据规定要求披露一些会产生影响的信息内容。导致披露的信息缺乏价值，这样不仅仅破坏了信息披露的规则要求，无法保障消费者的合法权益，也制约了公司的持续发展。

2 Principles of the MIR method

Here,P is the projection of the absorption coefficient,Δθ is the refraction angle,and σ2is the variance of USAXS.Here,(I0,ref,θ0,ref,σ20，ref)and(I0,obj,θ0,obj,σ20，obj)correspond to the intrinsic RC and the sample RC,respectively,and can be calculated by Eq.(2).

通过紫菜的单因素试验和正交试验，优选出微波功率、堆积厚度、微波时间等因素水平。采用L9（34）正交优化试验研究结果表明，影响紫菜加工工艺的主要因素为微波功率、微波时间，其次为堆积厚度。优选方案为A3B3C1，即微波时间12 min，堆积厚度5 mm，微波功率300 W。在此微波条件下，紫菜感官评分为93.2分，水分含量为5.0%。紫菜微波冻干工艺技术作为交叉学科发展的产物，在紫菜加工过程中能保持其色、香及水分含量，并在贮藏方面有重要的意义。微波冻干技术易于控制，效率高，在食品加工行业有着广阔的前景。

Fig.1(Color online)Diagram of the ABI setup.Here,(x,y,z)refer to the coordinate system fixed in an imaging system

The sample/intrinsic RC is obtained by acquiring multiple images(≥3)at different analyzer positions(corresponding to θ)with/without sample.Figure 2 shows that the sample RC is reduced,shifted,and broadened compared with the intrinsic RC.

Information on absorption,refraction,and USAXS can be extracted by comparing the sample RC and the intrinsic RC.For each pixel,the zeroth-, first-,and second-order moments are calculated with the following[27]:

Figure 6 shows that ε approaches 1 when the Poisson noise level is less than PN(1000),indicating that the effects of Poisson noise decrease when the detector counts are more than 1000 ph/pixel.For larger theoretical refraction angles(Δθt),acceptable results can still be extracted at higher Poisson noise levels,e.g.,PN(100).As a result,the detector counts should be more than 1000 ph/pixel in practical applications to confirm that the effects of Poisson noise are negligible.Under these circumstances,a photon flux of 2.5×106ph mm-1is required when the pixel unit of 20 μm × 20 μm is exploited.The required photon flux can be reduced with the same Poisson noise level when a larger pixel unit is used—increasing the pixel unit would also decrease the spatial resolution so that some appropriate trade-offs are seen in each speci fic case.For the SR sources,the photon flux is generally up to 107ph s-1 mm-2[28].Thus,Poisson noiselevelsbetterthan PN(10000)are easily implemented.However,the photon flux of conventional,low-intensity,X-ray sources is only about 104ph s-1mm-2[12,13];thus,an additional exposure time is needed for this case.In general,counts of 1000 are an ideal reference value.When the theoretical refraction angles are large(equivalent to FWHM),counts of 100 can reduce the radiation dose and exposure time with acceptable image quality(ε＞0.95;Fig.6a).These results confirm prior conclusions in the literature[28].Images can still be well identi fied with a count of 50.Of note,the dark current should be considered when the photon flux is too small—several temperature control modules can be exploited to reduce its impact[12].

摘 要：结合当前高校图书馆的服务情况，在分析“互联网+”的基本特性的基础上，从多方面论述并探讨了“互联网+”环境下高校图书馆存在的问题，据此提出了有针对性的多方面对策，希望能对今后“互联网+”下的高校图书馆发展所有帮助。

by measuring the intrinsic RC and the sample RC,the information of absorption,refraction,and USAXS can be obtained

中高煤阶煤基质孔隙表面具有很强的甲烷吸附能力，其能够吸附的气体量要远远大于煤层孔隙的气体含量，因此，煤层常为不饱和吸附而处于欠压状态。速度计算过程中需考虑压力和气体含量2个变量的影响。在煤层气储层处于饱和吸附状态时，压力变化主要影响煤基质固体颗粒弹性性质，与气体吸附引起的变化相比，这个变化较小可以忽略。在计算过程中主要考虑储层气体含量的变化，最大吸附气体量与压力关系根据Langmuir方程计算得出，即

Fig.2(Color online)The rocking curves of Si(333)at 59.31 keV with(sample RC)/without(intrinsic RC)objects.FWHM is the full width at half maximum of the intrinsic RC and is about 1 μrad.The Δθpand Δθcare the peak and centroid offsets between the two RCs,respectively.Here, Δθcis smaller than Δθpbecause parts of the sample RC are out of the sampling range

A sketch of the ABI system is shown in Fig.1.It is mainly composed of the monochromator and analyzer crystals.The M1and M2are monochromator crystals,which are mostly fabricated of perfect Si or Ge crystals.The X-rays become monochromatic and are collimated after re flecting from M2according to Bragg’s law.The analyzer is fabricated from the same monochromator crystals and is placed between the sample and detector.The X-rays will be recorded by the detector after being re flected by the analyzer when the X-rays are incident to the analyzer at the Bragg angle(θB)with a narrow range(typically on the order of a few microradians).Normally,the analyzer parallels with the monochromator crystals,and the initial incident angle of the X-ray is equal to θB.The analyzer rocks around the z-axis near θB,and the rocking angle is expressed as θ.Note that clockwise is de fined as the‘positive’ direction.The relationship ofre flectivity(R)and θ is the rocking curve(RC).The shape of the intrinsic RC is determined by the crystals and the X-ray energy.

3 Simulation conditions

The simulation conditions are as follows:The sample used in the simulation experiment is a cylinder—a common structure in biomedical and materials science.The cylinder is made of polystyrene(ρ=1.06 g/cm3),and its density is similar to soft tissue.The X-ray energy is 59.31 keV,i.e.,the Kα1characteristic energy of tungsten,and tungsten is the most applied target materials.Since the higher-energy X-rays have the stronger penetrating ability,the radiation dose delivered to the sample can be signi ficantly reduced.For instance,a 59.31-keV X-ray passes through 5 cm water with a transmittance of 35%,while only 0.65%of an 18 keV X-rays is transmitted[13].The absorption coef ficient of polystyrene is 0.1982 cm-1,and the refractive index decrement δ is ～ 6.7239× 10-8at 59.31 keV.Both the monochromator and the analyzer are Si(333)crystals because of the high angular sensitivity of Si(333)[13,22].Figure 2 shows the intrinsic RC of Si(333)at 59.31 keV.Its full width at half maximum(FWHM)is～1 μrad.The transverse and longitudinal sections of the sample are displayed in Fig.3a,b,respectively.The inside and outside diameters of the cylinder are 5 and 10 mm,respectively,and the height is 10 mm.The pixel unit of the detector is 20 μm × 20 μm.

Fig.3(Color online)a Transverse and b longitudinal section schematic diagram of the model.Different USAXS levels are de fined and indicated as A,B,C,D,and E

The deviation of sample RC produced by refraction will affect the extraction results within a certain and limited sampling range.To study the effects under different theoretical refraction angles,the parameter of kr= Δθt_max/A characterizes the degree of refraction,where A is 1 μrad at a sampling strategy of 1_0.2.The initial P is calculated based on the simulation conditions mentioned in Sect.3 with σ2=0.The extraction results of P with kr=1 and kr=0 are comparable far from the edge of the sample(Δθt is small),but the extraction results of P when kr=1 are much larger than that when kr=0 near the edge of the sample(Δθtis large)(Fig.10a).

以学校考研工作、大学生数学竞赛以及数学建模竞赛为切入点，以学校教学质量工程建设和创新强校工作为抓手，突出应用型、创新型和创业型人才培养，设计公共数学课程分类分级教学改革方案，探索具有地方高校特色的公共数学课程教学模式的改革与实践．每年组织校内大学生数学竞赛、数学建模竞赛活动和非数学类专业学生参加全国、全省大学生数学竞赛和数学建模竞赛等．采取多轨并行教学模式开展培训与竞赛活动，抓好大学生第二课堂活动，实施人才培养自主学习模式．组建数学建模协会，有计划地进行学校数学建模活动的普及性工作和参赛队员的初级培训，数学建模协会干部既是校内数学建模竞赛的组织者、赛前辅导者同时又是竞赛的参与者．

4 Results and discussion

4.1 Effects of Poisson noise

Different PN(n)s were used to investigate theeffects of Poisson noise on the extraction results of absorption,refraction,and USAXS,where n=10xand x=1:0.1:4.The sampling range was set as-FWHM to FWHM,i.e.,-1 to 1 μrad,and an interval of 0.1 μrad was adopted.To compare the results extracted under different Poisson noise levels(test images)and noise free(reference image),we de fine the coefficient of similarity ε as

where fiand giare the gray values of the ith pixel in the reference image and the test image,respectively.The faand gadenote the average gray values of the reference image and the test image.The two images are exactly the same when ε=1.

Fig.4(Color online)The theoretical refraction angles of the sample in the simulation experiment.The refraction angle is determined by the gradient of the phase shift,which is obvious at the edge of the objects

4.3.2 USAXS

where I(θj)refers to the X-ray intensity in a pixel for the angular setting θj,and N is the number of acquired images.The integrated intensity(I0),the center of mass(θ0),and the variance of the RC(σ20)can be expressed as

Fig.5 Extraction results with different Poisson noise levels.a,b Correspond to the results when the theoretical refraction angles are Δθtand 0.1Δθt,respectively.P is the projection of absorption coefficient,Δθ is the extracted refraction angle,and σ2is the variance of the USAXS.The extraction results are shown as the images of P,Δθ,and σ2

4.2 Effects of the sampling strategies

An appropriate sampling strategy(including range and interval)can improve the extraction results and reduce the radiation dose.Here,the sampling strategy is de fined in terms of(-A:s:A)and is abbreviated as A_s,which refers to the sampling range of-A to A μrad with an interval of s μrad.Consequently,the number of acquired images N is(2A/s+1).In our previous study[31],only the refraction information was considered,and the effects of Poisson noise were not considered.Here,the effects of the sampling strategies on extraction results are studied under two Poisson noise levels of PN(10000)and PN(10)while considering all three kinds of information.Three symmetric sampling strategies(1_0.1,2_0.2,and 1_0.2)were de fined by referring to the FWHM.The values of ε were calculated by comparing the theoretical images(reference images)and the extracted images(test images).The extraction results with PN(10000)and PN(10)are shown in Fig.7a,b,respectively.The values of ε are listed in Table 1.

At a low Poisson noise(PN(10000)),better extraction results are obtained using a wider sampling range,2_0.2,and an appropriate decrease in the sampling interval,1_0.2,barely impacts the extraction results compared with 1_0.1.However,the extraction results are more affected when the wider sampling range 2_0.2 is used at a high Poisson noise level(PN(10));image degradation is serious,as shown in Fig.7b.In addition,the ε of 2_0.2 has a more substantial reduction than that of 1_0.1 as shown in Table 1;this agrees with Fig.7.The deviation and broadening of the sample RC is produced by refraction and USAXS,respectively.Incomplete results will be obtained when a limited sampling range is employed(Fig.2).

2016年7月在德国汉堡召开的ICME-13上有两件大事值得一提，一是在会前的ICMI会员代表大会（General Assembly，GA）上，中国数学会推荐的候选人徐斌艳在九选五的差额选举中当选下一届ICMI EC成员，这是中国的第五位EC成员，中国大陆第四位EC成员；二是在ICME-13开幕式上，中国香港大学梁贯成教授正式被授予弗赖登塔尔奖（2013年度）．

The degree of incompleteness increased with enhanced refraction and USAXS and can generally be reduced by using a wider sampling range.However,since the wider sampling range(2_0.2)acquires data at low points of the sample RC,the number of photons incident on the detector is small due to the low re flectivity at these points.This corresponds to the high Poisson noise.The mass center of the sample RC will be strongly affected by these points because they match with large rocking angles θ.However,for the narrow sampling range(1_0.1 and 1_0.2),more photons are incident on the detector due to the relatively high re flectivity at the sampling points.This implies lower Poisson noise,particularly for the small theoretical values of refraction and USAXS.Of course,the effects of Poisson noise will increase if the number of acquired images,N,is lower(1_0.2).However,the effects of Poisson noise are negligible in the case of the Poisson noise level of PN(10000).

Fig.6(Color online)The relations between the coefficient of similarity,ε,and different Poisson noise levels when the theoretical refraction angles are a Δθtand b 0.1Δθt

The extraction results of refraction angle Δθ at the edge of the sample(the 51st column)are calculated via two noise levels,PN(10000)and PN(10),with theoretical refraction angles of Δθtand 0.1Δθt.The results extracted by 2_0.2 are better than that of 1_0.1 and 1_0.2 despite the larger data fluctuation when the theoretical refraction angles are Δθtwith Δθt_maxof 1.5 μrad(Fig.8).As illustrated in Fig.9,when the theoretical refraction angles are 0.1Δθtwith 0.1Δθt_maxof 0.15 μrad—which are smaller than FWHM by one order of magnitude—the values extracted by these three sampling strategies are comparable without USAXS.Meanwhile,1_0.1 has a better SNR than 2_0.2.Of note,the data fluctuation enlarges with the increase in USAXS,while the extracted values are reduced by USAXS.

4.3 Interaction of the three kinds of contrast information

Fig.7 Comparison of the results extracted by different sampling strategies(2_0.2,1_0.1,and 1_0.2)at two levels of Poisson noise a PN(10000)and b PN(10).The extraction results are shown as the images of P,Δθ,and σ2

Deriving from the principle of the MIR method,the contrast information of absorption,refraction,and USAXS is affected by each other.Some manifestations are in Sects.4.1 and 4.2;herein,a more systematic study is performed.To study the intrinsic shortcomings of the MIR method,the effects of Poisson noise are not considered.Section 4.2 and the conclusions from the literature[33]suggest a sampling strategy of 1_0.2,where 11 images should be acquired in total.

4.3.1 Refraction

Fig.9 The results of Δθ extracted at the edge of the sample using different sampling strategies(2_0.2,1_0.1,and 1_0.2)at two different Poisson noises:a PN(10000),b PN(10).The theoretical refraction angles are 0.1Δθt

Table 1 The coefficientofsimilarity with differentsampling strategies and Poisson noise

PN(10000) PN(10)P Δθ σ2 P Δθ σ2 2_0.2 0.3303 0.9277 0.8531 0.2622 0.7256 0.7012 1_0.1 0.1564 0.8130 0.7031 0.1419 0.6506 0.5943 1_0.2 0.1506 0.8060 0.6984 0.1393 0.5686 0.5364

Fig.8 The results of Δθ extracted at the edge of the sample using different sampling strategies(2_0.2,1_0.1,and 1_0.2)at two different Poisson noises:a PN(10000)and b PN(10).The theoretical refraction angles are Δθt

The absorption,refraction,and USAXS are all simultaneously considered in this simulation by assuming the cylinder has only absorption and refraction,while also considering that the levels of USAXS de fined by σ2can be controlled by setting the thickness of paper(no absorption and refraction)behind the sample[28,40].Figure 3b shows five regions(12 mm×2 mm)marked as A,B,C,D,and E that correspond to the five USAXS levels,i.e.,σ2=0,0.5,1,2,and 5 μrad2,respectively.Poisson noise is also added to study the effects of photon noise[28,41,42].Poisson noise level is negatively related to the number of photons incident on the detector—it can be de fined by the average number of photons n irradiating each pixel and can be described in terms of PN(n)when no sample is present.The initial counts of each pixel will fluctuate around n.After n is de fined,Poisson noise can be generated by the noise-generated function of poissrnd(n,p,q)in MATLAB.This function can generate Poisson random numbers with mean n,where scalars p and q are the row and column dimensions of the image;each Poisson random number corresponds to one pixel.The effects of Poisson noise become more complicated when the sample is inserted because of the sample’s absorption,refraction,and USAXS of X-rays.Note that each incident photon is assumed to produce one count[12,21].

This will illuminate the edges of the details(Fig.7a).The cylinder edges are obviously illuminated.The illumination partially enhances the contours of the sample,but the intrinsic image contrast signi ficantly declined.Thus,wrong quantitative information will be extracted,which agrees with the literature[30].The extraction results of Δθ when kr=0.5 and kr=1 are shown in Fig.10b.These are marked as Δθ_0.5and Δθ_1,respectively.At the edge of the sample,Δθ_1＜2Δθ_0.5,indicating that the accuracy of the Δθ declines with increasing Δθt.Besides,as Fig.10c shows,several additional USAXS signals are extracted due to the effect of refraction and presented at the edges of the sample as ‘negative peaks’.This leads to the dark lines in the USAXS image in Fig.7a.

The results in Fig.11 show the dependence of the extraction results on kr,wherein the extraction value corresponds to the 51st pixel in Fig.10.Figure 11 shows that P remains mostly constant until kr～0.5.It then increases with kr.The Δθ increases linearly until kr ～ 0.6.When kr increases further,the Δθ gradually becomes smooth and tends to 1 μrad,which is equal to A(the sampling range of the sampling strategy of 1_0.2).This indicates that the extracted maximum value of Δθ is related to the sampling range.Here,σ2decreases with kr,and the changes in σ2are not obvious if kr＜0.15.These results show that as kr increases(denoting the degree of refraction),the effects on the extraction results due to refraction behave as σ2＞P＞Δθ within a certain and limited sampling range.

Fig.10(Color online)The effects of refraction on the extraction results of a P,b Δθ,and c σ2(the 250th row).Here,krcharacterizes the degree of refraction

Fig.11(Color online)The relationships between refraction and the extraction results of P,Δθ,and σ2

Figure 4 shows the theoretical refraction angles Δθtof the sample in the simulation experiments.The maximum Δθt,i.e.,Δθt_max,is ～ 1.5 μrad.This is observed at the edge of the sample and is comparable to the FWHM(～ 1 μrad).Two different orders of theoretical refraction angles(Δθtand 0.1Δθt)are employed to compare the effects of Poisson noise.The 0.1Δθt_maxis smaller than FWHM by one order of magnitude.The extracted contrast images(P,Δθ,and σ2)with different Poisson noise values(PN(10),PN(100),PN(1000),PN(10000),and noise free)when the theoretical refraction angles are set as Δθtand 0.1Δθtare shown in Fig.5a,b,respectively.The images are degraded and blurred with an increasing Poisson noise.The ε of the images at different levels of Poisson noise is shown in Fig.6.The degree of image degradation caused by Poisson noise is Δθ＞σ2＞P.Obviously,refraction is more sensitive to Poisson noise when the theoretical refraction angles are small(0.1Δθt).This agrees with the results presented in Fig.5.

据了解，本次技能大赛分为初赛、决赛两部分，60%的参赛人员通过初赛考核，可进入决赛环节。在决赛之前，青岛市家庭服务业生态圈统一组织入围参赛者接受为期1个月的培训。通过决赛阶段综合笔试和实操两项成绩，最终评选出“首席家庭服务员”8人、“首席婴幼儿营养烹饪师”“首席母婴护理师”各10名，共28位首席技师，优胜奖50名，由青岛市家庭服务产业生态圈颁发荣誉证书并授予称号，举行优胜获奖者与客户供需见面会，为岛城家庭服务市场输送优质服务人员。

Fig.12(Color online)The effects of USAXS on the extraction results of a P,b Δθ,and c σ2.Here,kscharacterizes the degree of USAXS

The results are affected by the broadening of sample RC caused by USAXS.Analogous to Sect.4.3.1,the parameter ks=2.36σ/A is de fined,where σ is the standard deviation of USAXS.The initial P is calculated based on the simulation conditions mentioned in Sect.3.The theoretical refraction angles are set to 0.1Δθt,and in this case,the effects of refraction are negligible.A comparison of the extraction results when ks=0 and ks=1 is shown in Fig.12.The values of P increased overall because of the effect of USAXS,and the degree of the increase depends on the level of USAXS.This will cause layering,contrast reduction,and image distortion in the absorption image(Fig.7a).In Fig.12b,USAXS will make Δθ smaller than the theoretical one—especially at the edges of the sample.Figure 7 shows that Δθ decreased signi ficantly with strong USAXS.Figure 12c indicates that the extraction result of σ2 (0.06 μrad2) is less than the theoreticalvalue(0.18 μrad2)that is due to the limitations of the sampling range.The effects of refraction can be neglected if the theoretical refraction angles are 0.1Δθtand correspond to kr=0.15(Fig.12).

The extraction results as a function of ksare shown in Fig.13,wherein the extraction value corresponds to the 51st pixel in Fig.12.The results in Fig.13 show that P increases with ks,and the increment is negligible if ks＜0.5.Here,Δθ decreases slowly when ks＜0.4,and it then obviously decreases with increasing ks.It finally becomes constant at 0 when ksis larger than 3.The σ2 changes inversely as Δθ,and the maximum value extracted,approaches 0.18 μrad2due to the limitation of the sampling range.We can infer that the USAXS leads to layering and distortion in the P and Δθ image(Fig.7).The extraction results of σ2are inaccurate if the USAXS is strong.

4.3.3 Absorption

再者，《劳动合同法》第六十八、第六十九条、第七十一条分别规定：“非全日制用工，是指以小时计酬为主，劳动者在同一用人单位一般平均每日工作时间不超过四小时，每周工作时间累计不超过二十四小时的用工形式。”“全日制用工双方当事人可以订立口头协议。从事非全日制用工的劳动者可以与一个或者一个以上用人单位订立劳动合同。”“非全日制用工双方当事人任何一方都可以随时通知对方终止用工。”

Fig.13(Color online)The relationship between USAXS and the extraction results of P,Δθ,and σ2

Ideally,the sample absorption of X-ray does not affect the extraction results of refraction and USAXS because of the information extraction characteristics of the MIR method.However,samples that have strong absorption of X-rays will cause the Poisson noise level to increase due to the large reduction in the number of photons incident on the detector.Thus,additional exposure is needed.In general,phase-contrast imaging techniques are applied in materials made of light elements,where a low contrast is obtained via absorption-contrastimaging.Here,the absorption is weak—especially when using high X-ray energy(59.31 keV).For instance,only 18%of the X-rays can be absorbed by 10-mm polystyrene with an absorption coefficient of 0.1982 cm-1.As a consequence,the effects of absorption on the extraction results can be somewhat neglected.

其中，Pi为单因子污染指数，Cis为重金属浓度实测值，Cin为重金属参比值，Eir为单因子危害系数，Tir为毒性响应系数。

5 Conclusion

In this paper,the characteristics of MIR information extraction were studied via a simulation experiment.The results suggest that Poisson noise has signi ficant effects on information extraction.As the Poisson noise increases,the results obviously degraded,and the image blurs.The effects of the Poisson noise on refraction,USAXS,and absorption are sequentially increased.The effects of Poisson noise are insigni ficant when the Poisson noise level is better than PN(1000).When the theoretical refraction angle is large,the results are slightly affected by Poisson noise.Acceptable results can still be obtained even if the Poisson noise level is as high as PN(100).Here,the requirements of radiation dose and imaging time are reduced.

For SR sources,the effects due to Poisson noise are negligible at high photon flux:The noise level is normally better than(PN(10000)).On the contrary,X-ray tube sources have Poisson noise levels that are high due to their low intensity.Thus,they require additional exposure and better system stability and temperature control.When Poisson noise is negligible(e.g.,PN(10000)),the sampling interval barely impacts the information extraction results.More accurate results can be achieved by increasing the sampling range.The interaction of the three kinds of information can also be reduced if a wide sampling range is employed.However,in the case of high Poisson noise level(e.g.,PN(10)),the results extracted with a narrow sampling range have better SNR.Consequently,the selection of a sampling strategy should be based on the Poisson noise level and the theoretical value to optimize the extraction results and reduce imaging time.

The three kinds of information will affect each other during extraction.Across a certain and limited sampling range,the refraction will cause bright and dark edges in the absorption image and USAXS image,respectively.There are incomplete results for the refraction image.This will decrease image contrast,and the resulting quantitative results will be incorrect.The sampling range and the degree of refraction determine the refraction effects.The maximum extracted value is related to the sampling range.The effects of the extraction results due to refraction can be neglected when the theoretical value of the refraction angle is smaller than the sampling range by an order of magnitude(e.g.,kr=0.15).The USAXS of X-rays will lead to the layering of the absorption image and a decrease in the image contrast,as well as a reduction in the refraction results with incomplete USAXS results.The effects of USAXS can be ignored when the standard deviation of USAXS σ＜0.4(A/2.36),i.e.,ks＜0.4.In general,it is unnecessary to consider the effects of sample absorption on the results of information extraction;however,an appropriate increase in radiation exposure should be performed to reduce the Poisson noise caused by the absorption to X-rays when the sample is made of high-Z elements and/or is large.

References

1.Y.Z.Zhao,E.Brun,P.Coan et al.,High-resolution,low-dose phase contrast X-ray tomography for 3D diagnosis of human breast cancers.Proc.Natl.Acad.Sci.USA 109,18290–18294(2012).https://doi.org/10.1073/pnas.1204460109

2.H.S.Rocha,G.R.Pereira,P.Faria et al.,Diffraction-enhanced imaging microradiography applied in breast samples.Eur.J.Radiol.68,S37–S40(2008).https://doi.org/10.1016/j.ejrad.2008.04.032

3.C.Muehleman,J.Li,Z.Zhong et al.,Multiple-image radiography for human soft tissue.J.Anat.208,115–124(2006).https://doi.org/10.1111/j.1469-7580.2006.00502.x

4.P.Coan,J.Mollenhauer,A.Wagner et al.,Analyzer-based imaging techniquein tomography ofcartilageand metal implants:a study at the ESRF.Eur.J.Radiol.68,S41–S48(2008).https://doi.org/10.1016/j.ejrad.2008.04.036

5.P.Suortti,J.Keyrila¨inen,W.Thomlinson,Analyser-based x-ray imaging for biomedical research.J.Phys.D Appl.Phys.46,494002(2013).https://doi.org/10.1088/0022-3727/46/49/494002

6.M.J.Kitchen,D.M.Paganin,K.Uesugi et al.,X-ray phase,absorption and scatter retrieval using two or more phase contrast images.Opt.Express 18,19994(2010).https://doi.org/10.1364/OE.18.019994

7.C.H.Hu,T.Zhao,L.Zhang et al.,Information extraction and CT reconstruction of liver images based on diffraction enhanced imaging.Prog.Nat.Sci.19,955–962(2009).https://doi.org/10.1016/j.pnsc.2008.06.031

8.M.J.Kitchen,K.M.Pavlov,S.B.Hooper et al.,Simultaneous acquisition of dual analyser-based phase contrast X-ray images for small animal imaging.Eur.J.Radiol.68,S49–S53(2008).https://doi.org/10.1016/j.ejrad.2008.04.028

9.X.Zhang,X.R.Yang,Y.Chen et al.,Visualising liver fibrosis by phase-contrast X-ray imaging in common bile duct ligated mice.Eur.Radiol.23,417–423(2013).https://doi.org/10.1007/s00330-012-2630-z

10.D.Chapman,W.Thomlinson,Z.Zhong et al.,Diffraction enhanced imaging applied to materials science and medicine.Synchrotron Radiat News 11,4–11(1998).https://doi.org/10.1080/08940889808260849

11.T.M.Wang,J.J.Xu,J.Li et al.,In situ study on dendrite growth of metallic alloy by a synchrotron radiation imaging technology.Sci.China Technol.Sci.53,1278–1284(2010).https://doi.org/10.1007/s11431-010-0087-3

12.W.Zhou,K.Majidi,J.G.Brankov,Analyzer-based phase-contrast imaging system using a micro focus X-ray source.Rev.Sci.Instrum.85,085114(2014).https://doi.org/10.1063/1.4890281

13.C.Parham,Z.Zhong,D.M.Connor et al.,Design and implementation of a compact low-dose diffraction enhanced medical imaging system.Acad.Radiol.16,911–917(2009).https://doi.org/10.1016/j.acra.2009.02.007

14.I.Nesch,D.P.Fogarty,T.Tzvetkov et al.,The design and application of an in-laboratory diffraction-enhanced x-ray imaging instrument.Rev.Sci.Instrum.80,093702(2009).https://doi.org/10.1063/1.3213621

15.D.J.Vine,D.M.Paganin,K.M.Pavlov et al.,Analyzer-based phase contrast imaging and phase retrieval using a rotating anode x-ray source.Appl.Phys.Lett.91,254110(2007).https://doi.org/10.1063/1.2825426

16.A.Snigirev,I.Snigireva,V.Kohn et al.,On the possibilities of x-ray phase contrast microimaging by coherent high-energy synchrotron radiation.Rev.Sci.Instrum.66,5486–5492(1995).https://doi.org/10.1063/1.1146073

17.S.W.Wilkins,T.E.Gureyev,D.Gao et al.,Phase-contrast imaging using polychromatic hard X-rays.Nature 384,335–338(1996).https://doi.org/10.1038/384335a0

18.C.David,B.No¨hammer,H.H.Solak et al.,Differential x-ray phase contrast imaging using a shearing interferometer.Appl.Phys.Lett.81,3287(2002).https://doi.org/10.1063/1.1516611

19.T.Weitkamp,A.Diaz,C.David et al.,X-ray phase imaging with a grating interferometer.Opt.Express 13,6296–6304(2005).https://doi.org/10.1364/OPEX.13.006296

20.F.Pfeiffer,T.Weitkamp,O.Bunk et al.,Phase retrieval and differential phase-contrast imaging with low-brilliance X-ray sources.Nat.Phys.2,258–261(2006).https://doi.org/10.1038/nphys265

21.P.C.Diemoz,A.Bravin,M.Langer et al.,Analytical and experimental determination of signal-to-noise ratio and figure of merit in three phase-contrast imaging techniques.Opt.Express 20,27670–27690(2012).https://doi.org/10.1364/OE.20.027670 22.D.Chapman,W.Thomlinson,R.Johnston et al.,Diffraction enhanced x-ray imaging.Phys.Med.Biol.42,2015(1997).https://doi.org/10.1088/0031-9155/42/11/001

23.A.Maksimenko,Nonlinear extension of the x-ray diffraction enhanced imaging.Appl.Phys.Lett.90,154106(2007).https://doi.org/10.1063/1.2721378

24.C.H.Hu,L.Zhang,H.Li et al.,Comparison of refraction information extraction methods in diffraction enhanced imaging.Opt.Express 16,16704–16710(2008).https://doi.org/10.1364/OE.16.016704

25.L.Rigon,F.Arfelli,R.H.Menk,Three-image diffraction enhanced imaging algorithm to extract absorption,refraction,and ultrasmall-angle scattering.Appl.Phys.Lett.90,114102(2007).https://doi.org/10.1063/1.2713147

26.O.Oltulu,Z.Zhong,M.Hasnah et al.,Extraction of extinction,refraction and absorption properties in diffraction enhanced imaging.J.Phys.D Appl.Phys.36,2152–2156(2003).https://doi.org/10.1088/0022-3727/36/17/320

27.E.Pagot,P.Cloetens,S.Fiedler et al.,A method to extract quantitative information in analyzer-based x-ray phase contrast imaging.Appl.Phys.Lett.82,3421(2003).https://doi.org/10.1063/1.1575508

28.M.N.Wernick,O.Wirjadi,D.Chapman et al.,Multiple-image radiography.Phys.Med.Biol.48,3875–3895(2003).https://doi.org/10.1088/0031-9155/48/23/006

29.Y.I.Nesterets,P.Coan,T.E.Gureyev et al.,On qualitative and quantitative analysis in analyser-based imaging.Acta.Crystallogr.Sect.A 62,296–308 (2006).https://doi.org/10.1107/S0108767306017843

30.P.C.Diemoz,P.Coan,C.Glaser et al.,Absorption,refraction and scattering in analyzer-based imaging:comparison of different algorithms.Opt.Express 18,3494–3509(2010).https://doi.org/10.1364/OE.18.003494

31.Y.B.Wang,G.P.Li,X.D.Pan et al.,Simulation of X-ray refraction information extraction using multiple image-collecting strategies.Acta Phys.Sin.63,150–161(2014).https://doi.org/10.7498/aps.63.104206.(in Chinese)

32.B.Marquet,J.G.Brankov,M.N.Wernick,Noise and sampling analysis for multiple-image radiography,in Proc.3rd IEEE International Symposium on Biomedical Imaging:Nano To Macro,Virginia,CA,United States,2006,pp.1232–1235.https://doi.org/10.1109/isbi.2006.1625147

33.K.Majidi,J.G.Brankov,M.N.Wernick,Sampling strategies in multiple-image radiography,in Proc.5th IEEE International Symposium on Biomedical Imaging:Nano To Macro,Pairs,CA,France,2008,pp.688–691.https://doi.org/10.1109/isbi.2008.4541089

34.C.H.Hu,T.Zhao,H.Li et al.,The study of phase information extraction and fusion based on diffraction enhanced imaging.J.Image Graph.13,1622–1628(2008).https://doi.org/10.11834/jig.20080838.(in Chinese)

35.Z.H.Huang,K.J.Kang,Y.G.Yang,Extraction methods of phase information for X-ray diffraction enhanced imaging.Nucl.Instrum.Methods A 579,218–222(2007).https://doi.org/10.1016/j.nima.2007.04.043

36.Z.Q.Chen,F.Ding,Z.F.Huang et al.,Polynomial curve fitting method for refraction-angle extraction in diffraction enhanced imaging.Chin.Phys.C 33,969–974(2009).https://doi.org/10.1088/1674-1137/33/11/008

37.P.Y.Li,K.Zhang,W.X.Huang et al.,Cosine fitting radiography and computed tomography.Chin.Phys.B 24,068704(2015).https://doi.org/10.1088/1674-1056/24/6/068704

38.X.J.Zhao,K.Zhang,Y.L.Hong et al.,A simple method of extracting multiple-information with diffraction enhanced imaging.Acta Phys.Sin.62,124202(2013).https://doi.org/10.7498/aps.62.124202.(in Chinese)

39.Y.B.Wang,G.P.Li,X.D.Pan et al.,An improved multiple-image radiography method extracting refraction information in analyzer-based imaging.Nucl.Instrum.Methods A 770,182–188(2015).https://doi.org/10.1016/j.nima.2014.10.035

40.C.Y.Chou,M.A.Anastasio,J.G.Brankov et al.,An extended diffraction-enhanced imaging method for implementing multipleimage radiography.Phys.Med.Biol.52,1923–1945(2007).https://doi.org/10.1088/0031-9155/52/7/011

41.L.Rigon,F.Arfelli,R.H.Menk,Generalized diffraction enhanced imaging to retrieve absorption,refraction and scattering effects.J.Phys.D Appl.Phys.40,3077–3089(2007).https://doi.org/10.1088/0022-3727/40/10/011

42.J.G.Brankov,M.N.Wernick,Y.Yang et al.,A computed tomography implementation of multiple-image radiography.Med.Phys.33,278–289 (2006).https://doi.org/10.1118/1.2150788