本發(fā)明涉及圖像解析裝置、圖像解析方法以及程序，特別是適用于對纖維增強復合材料中包含的纖維束取向進行解析的圖像解析裝置、圖像解析方法以及程序。

背景技術：

近年來，正在進行纖維增強復合材料的一種即陶瓷基復合材料(cmc：ceramicmatrixcomposites)的開發(fā)。cmc是通過母材(基質)將陶瓷纖維增強而得到的復合材料，具有輕質且耐熱性優(yōu)良這種的特征。正在研究利用該特征例如將cmc用于飛機發(fā)動機部件，當前正在進行實用化。通過將cmc用作飛機發(fā)動機部件，能夠期待大幅提高燃料經濟性。

以下是cmc的一般的形成工藝。首先，捆扎幾百根左右的陶瓷纖維來制作纖維束，編織該纖維束來制作編織物。關于纖維束的編織方法，例如具有三維編織或平紋編織。三維編織是在xyz方向這三個方向上編織纖維束來制作編織物的方法，平紋編織為在xy方向這兩個方向上編織纖維束來制作編織物的方法。

在制作編織物之后，通過cvi(chemicalvaporinfiltration：化學氣相滲透法)和pip(polymerimpregnationandpyrolysis：浸漬裂解工藝)來形成基質，最后進行機械加工和表面涂層等，由此形成cmc。在此，此時形成的cmc中的纖維束取向對cmc的強度影響大。

即，當纖維束在本來應該為直線的地方蜿蜒，或者整體偏離本來應該配置的基準軸從而發(fā)生偏向，或者在途中斷裂的情況下，cmc的強度降低。與此相對，在纖維束不發(fā)生蜿蜒、偏向或斷裂而在固定方向上使朝向一致地來恰當地進行排列時，強度高且耐熱性優(yōu)良。因此為了確認所形成的cmc的強度是否足夠，重要的是評價纖維束取向。

在專利文獻1中公開了以下取向解析方法：對樹脂成型品的切片圖像進行二值化來獲取二值圖像，對該二值圖像進行傅里葉變換來獲取能譜圖像，將與根據該能譜圖像描繪的橢圓正交的橢圓的主軸方向設為樹脂成型品中包含的填充劑(纖維)的取向方向。

另外，在非專利文獻1中公開了以下技術：通過x射線ct裝置對編織纖維束而得到的編織物進行拍攝來獲取x射線ct圖像，使用特殊的濾波函數對該x射線ct圖像進行計算，由此對構成纖維束的一根一根纖維的取向進行解析。

在先技術文獻

專利文獻1：日本特開2012-2547號公報

非專利文獻1：t.shinohara,j.takayama,s.ohyama,anda.kobayashi,“extractionofyarnpositionalinformationfromathree‐dimensionalctimageoftextilefabricusingyarntracingwithafilamentmodelforstructureanalysis”,textileresearchjournal,vol.80,no.7,pp.623-630(2010)

技術實現要素：

但是，在專利文獻1所記載的技術中，關于切片圖像內包含的填充劑(纖維)的取向，作為解析結果只能得到一個方向。因而，例如在如三維編織、平紋編織那樣纖維束在多個方向上排列的情況下，無法得到各個纖維束取向來作為解析結果。另外，無法解析纖維束是否不會蜿蜒、偏向或斷裂而在固定方向上使朝向一致地來適當地排列。

另外，在非專利文獻1所記載的技術中，需要獲取能夠一根一根地識別構成纖維束的纖維的程度的高分辨率的x射線ct圖像。在該情況下，用于得到x射線ct圖像的拍攝時間變長，無法用于產品檢查，因此并不實用。另外，是針對截面為圓形的纖維有效的技術，無法作為對截面為扁平的纖維束取向進行解析的技術而直接使用。并且存在需要輸入x射線ct圖像中的各纖維的開始點從而操作變得繁雜這種問題。

本發(fā)明是考慮上述問題點而做出的，提出一種從cmc的三維圖像能夠容易地解析纖維束取向的圖像解析裝置、圖像解析方法以及程序。

為了解決該課題，本發(fā)明公開一種圖像解析裝置，其從由x紗線、y紗線以及z紗線的纖維束構成的編織物的三維圖像解析x紗線和y紗線的纖維束的取向，其特征在于，具備：二值化部，其對三維圖像進行二值化；重疊區(qū)域提取部，其從二值化后的圖像提取x紗線和y紗線垂直地立體交叉的重疊區(qū)域；參考方向決定部，其使重疊區(qū)域中包含的各體素的重疊方向平均化，將平均化后的方向決定為參考方向；z紗線去除部，其在相對于參考方向垂直的參考平面上應用方向性距離法，從二值化后的圖像中去除z紗線；以及纖維束取向推定部，其針對去除了z紗線的圖像在參考平面上再次應用方向性距離法，根據應用時計算的方向性距離，推定x紗線和y紗線的纖維束取向。

另外，為了解決該課題，本發(fā)明公開了一種圖像解析方法，其從由x紗線、y紗線以及z紗線的纖維束構成的編織物的三維圖像解析x紗線和y紗線的纖維束取向，其特征在于，計算機具有以下步驟：第一步驟，對三維圖像進行二值化；第二步驟，從二值化后的圖像中提取x紗線和y紗線垂直地立體交叉的重疊區(qū)域；第三步驟，使重疊區(qū)域內包含的各體素的重疊方向平均化，將平均化后的方向決定為參考方向；第四步驟，在相對于參考方向垂直的參考平面上應用方向性距離法，從二值化后的圖像中去除z紗線；以及第五步驟，針對去除了z紗線的圖像在參考平面上再次應用方向性距離法，根據應用時計算的方向性距離，推定x紗線和y紗線的纖維束取向。

另外，為了解決該課題，本發(fā)明公開一種程序，其從由x紗線、y紗線以及z紗線的纖維束構成的編織物的三維圖像解析x紗線和y紗線的纖維束取向，其特征在于，程序使計算機執(zhí)行以下步驟：第一步驟，對三維圖像進行二值化；第二步驟，從二值化后的圖像中提取x紗線和y紗線垂直地立體交叉的重疊區(qū)域；第三步驟，使重疊區(qū)域內包含的各體素的重疊方向平均化，將平均化后的方向決定為參考方向；第四步驟，在相對于參考方向垂直的參考平面上應用方向性距離法，從二值化后的圖像中去除z紗線；以及第五步驟，針對去除了z紗線的圖像在參考平面上再次應用方向性距離法，根據應用時計算的方向性距離，推定x紗線和y紗線的纖維束取向。

另外，為了解決該課題，公開的圖像解析裝置的特征在于，具備：二值化部，其對由x紗線、y紗線以及z紗線的纖維束構成的編織物的三維圖像進行二值化；重疊區(qū)域提取部，其從上述二值化后的圖像中提取上述x紗線和上述y紗線立體交叉的重疊區(qū)域；以及重疊區(qū)域形態(tài)解析部，其解析上述提取的重疊區(qū)域的形態(tài)。

另外，為了解決該課題，公開的圖像解析方法的特征在于，具備以下步驟：對由x紗線、y紗線以及z紗線的纖維束構成的編織物的三維圖像進行二值化的步驟；從上述二值化后的圖像中提取上述x紗線和上述y紗線立體交叉的重疊區(qū)域的步驟；以及解析上述提取的重疊區(qū)域的形態(tài)的步驟。

另外，為了解決該課題，公開的程序使計算機執(zhí)行以下步驟：對由x紗線、y紗線以及z紗線的纖維束構成的編織物的三維圖像進行二值化的步驟；從上述二值化后的圖像中提取上述x紗線和上述y紗線立體交叉的重疊區(qū)域的步驟；以及解析上述提取的重疊區(qū)域的形態(tài)的步驟。

根據公開的本發(fā)明，能夠從cmc的三維圖像容易地解析纖維束取向。

附圖說明

圖1是圖像解析裝置的整體結構圖。

圖2是三維編織的編織物的概念圖。

圖3是圖像解析處理的整體流程圖。

圖4是二值化處理的詳細流程圖。

圖5是通過二值化處理生成的處理圖像。

圖6是重疊區(qū)域提取處理的詳細流程圖。

圖7是通過重疊區(qū)域提取處理生成的處理圖像。

圖8是參考方向決定處理的詳細流程圖。

圖9是應用通常的方向性距離法來推定方向性距離的處理的概要圖。

圖10是通過參考方向決定處理生成的處理圖像。

圖11是z紗線去除處理的詳細流程圖。

圖12是應用了通常的方向性距離法時的概念圖。

圖13是附帶參考的方向性距離法的概念圖。

圖14是應用了附帶參考的方向性距離法時的概念圖。

圖15是對附帶條件的方向性距離進行推定的處理的概要圖。

圖16是表示x紗線或y紗線的區(qū)域與z紗線的區(qū)域的概念圖。

圖17是纖維束取向推定處理的詳細流程圖。

圖18是輸入圖像為模擬圖像時的解析結果。

圖19是輸入圖像為高分辨率的x射線ct圖像時的解析結果。

圖20是輸入圖像為低分辨率的x射線ct圖像時的解析結果。

圖21是輸入圖像為缺陷模擬圖像時的解析結果。

圖22是輸入圖像為曲面模擬圖像時的解析結果。

圖23是圖像解析裝置的整體結構圖。

圖24是圖像解析處理的整體流程圖。

圖25表示支柱的體積的計算例。

圖26表示支柱伸長的方向的計算例。

圖27表示支柱彎曲時的圖像處理。

圖28表示支柱的一部分不規(guī)則地排列的例子。

圖29是通過xy平面表示x紗線和y紗線的纖維束的理想的配置。

圖30示意性地表示圖29的纖維束的配置中的xz平面(或yz平面)上的重疊區(qū)域提取圖像。

圖31示意性地表示圖29的纖維束的配置中的xy平面上的重疊區(qū)域提取圖像。

圖32表示在xy平面上纖維束的一部分在與應該延伸的方向相反的xy平面內的方向上延伸的狀態(tài)。

圖33示意性地表示圖32的纖維束的狀態(tài)下的yz平面的重疊區(qū)域提取圖像。

圖34示意性地表示圖32的纖維束的狀態(tài)下的xy平面的重疊區(qū)域提取圖像。

圖35表示在xy平面上纖維束的一部分在與應該延伸的方向相反的z方向上延伸的狀態(tài)。

圖36示意性地表示圖35的纖維束的狀態(tài)下的yz平面的重疊區(qū)域提取圖像。

圖37表示在xy平面上纖維束的一部分在與應該延伸的方向不同的一同具有xy平面內的方向成分和z方向成分的方向上進行延伸的狀態(tài)。

圖38示意性地表示圖37的纖維束的狀態(tài)下的、yz平面的重疊區(qū)域提取圖像。

圖39表示在xy平面上纖維束作為整體分別在傾斜的方向上延伸的狀態(tài)。

圖40表示在xy平面上纖維束作為整體分別在傾斜的方向上延伸的狀態(tài)。

圖41示意性地表示圖39或40所示的纖維束的狀態(tài)下的、xy平面的重疊區(qū)域提取圖像。

圖42表示在xy平面上不存在纖維束的一部分(例如x紗線)的狀態(tài)的圖。

圖43表示在xy平面上纖維束的一部分的粗細比其它纖維束粗或細的情況。

圖44示意性地表示圖42或圖43的纖維束的狀態(tài)下的、yz平面的重疊區(qū)域提取圖像。

圖45表示在xy平面上纖維束的一部分(例如x紗線)折返的狀態(tài)。

圖46示意性地表示圖45的纖維束的狀態(tài)下的、xy平面的重疊區(qū)域提取圖像。

圖47是圖像解析裝置的整體結構圖。

圖48是圖像解析處理的整體流程圖。

具體實施方式

以下，根據附圖詳細說明本發(fā)明的一個實施方式。

(1)圖像解析裝置1的整體結構

圖1表示本實施方式的圖像解析裝置1的整體結構。圖像解析裝置1為具備cpu(centralprocessingunit：中央處理器)11、輸入部12、存儲部13、顯示部14以及存儲器15而構成的計算機。

cpu11是與存儲器15中存儲的各種程序進行協(xié)作來綜合控制圖像解析裝置1的動作的處理器。輸入部12為接受來自用戶的輸入的接口，例如為鍵盤、鼠標。本實施方式的輸入部12可以是用于輸入構成cmc(ceramicmatrixcomposites：陶瓷基復合材料)的編織物的x射線ct圖像g10的接口。

在此，cmc是指捆扎幾百根左右的陶瓷纖維來制作纖維束，在編織該纖維束而作成編織物之后用碳等覆蓋纖維表面，并經過用于形成基質的cvi(chemicalvaporinfiltration化學氣相滲透)工藝和pip(polymerimpregnationandpyrolysis浸漬裂解)工藝等而形成的纖維增強復合材料。

在制作編織物時的編織方法中存在被稱為三維編織或平紋編織的方法。三維編織為從xyz方向這三個方向編織纖維束來制作編織物的方法，平紋編織為從xy方向這兩個方向上編織纖維束來制作編織物的方法。

圖2表示三維編織的編織物的概念圖。在本實施方式中，通過x射線ct裝置對特別是通過復雜的三維編織而制作出的編織物或由該編織物形成的cmc進行拍攝，根據得到的x射線ct圖像g10來自動解析纖維束取向。

如圖2所示，將由等間隔排列的x紗線的纖維束形成的纖維層以及由等間隔排列的y紗線的纖維束形成的纖維層交替地進行層疊，并通過z紗線的纖維束進行固定以使層疊的很多纖維層不會在cmc形成過程中剝離，從而形成三維編織的編織物。

關于通過該編織物形成的cmc，通常預想向x紗線方向或y紗線方向進行伸縮來進行設計。由此，與x紗線和y紗線大致垂直而編織的z紗線不對直接影響cmc的強度。另一方面，關于z紗線的存在，在對x紗線和y紗線的取向進行解析時，有時造成精度不良。

因此，在本實施方式中，想要從三維編織的編織物的x射線ct圖像g10中去除z紗線，高精度地解析x紗線和y紗線的纖維束取向。

此外，通常，取向是指在固定方向上使朝向一致地進行排列或進行了排列的狀態(tài)的術語，在本實施方式中也以相同的意思來使用。即使在蜿蜒、偏向或斷裂地排列的狀態(tài)下，如果是在固定方向上使朝向一致地進行了排列的狀態(tài)則也稱為取向。

返回至圖1，存儲部13為用于存儲從輸入部12輸入的x射線ct圖像g10和對該x射線ct圖像g10實施了各種圖像處理后的處理圖像的存儲介質。顯示部14為顯示x射線ct圖像g10圖像和處理圖像的lcd(liquidcrystaldisplay：液晶顯示器)等顯示裝置。例如，顯示部14從x射線ct圖像g10中自動解析纖維束取向，將纖維束取向推定圖像g100顯示在顯示畫面上。

存儲器15為用于存儲各種程序的存儲介質，各種程序與cpu11進行協(xié)作來執(zhí)行圖像解析處理。在各種程序中具有二值化部151、重疊區(qū)域提取部152、參考方向決定部153、z紗線去除部154以及纖維束取向推定部155。在后文說明通過這些各種程序而執(zhí)行的圖像解析處理(圖3)。

(2)圖像解析處理的流程圖

圖3表示本實施方式的圖像解析處理p1的整體流程圖。以輸入部12接受了來自用戶的執(zhí)行指示為契機，通過cpu11與存儲在存儲器15中的各種程序的協(xié)作來執(zhí)行該圖像解析處理。以下，為了便于說明，將處理主體作為各種程序來進行說明。

首先，當經由輸入部12輸入了x射線ct圖像g10時(sp1)，二值化部151根據規(guī)定閾值對輸入的x射線ct圖像g10進行二值化，生成在前景中表示x紗線、y紗線以及z紗線的各纖維束的二值圖像(sp2)。

接著，重疊區(qū)域提取部152提取x紗線與y紗線垂直地立體交叉的重疊區(qū)域(sp3)，參考方向決定部153將提取出的重疊區(qū)域的重疊方向決定為參考方向(sp4)。

在此，提取重疊區(qū)域的理由在于，為了對提取出的重疊區(qū)域應用通常的方向性距離法，推定重疊方向。在后文中說明通常的方向性距離法。

另外，將重疊方向決定為參考方向的理由在于，為了在與該參考方向垂直的平面上應用二維的通常的方向性距離法，來推定纖維束取向。

將與參考方向垂直的平面稱為參考平面，將在參考平面上應用二維的通常的方向性距離法的方法稱為附帶參考的方向性距離法。

在參考平面上存在x紗線或y紗線，因此通過應用附帶參考的方向性距離法，能夠高精度地推定x紗線和y紗線的纖維束取向。在后文中說明附帶參考的方向性距離法。

接著，z紗線去除部154對二值圖像應用附帶參考的方向性距離法，推定x紗線、y紗線以及z紗線的纖維束取向。然后，z紗線去除部154根據推定出的取向來去除二值圖像中包含的z紗線(sp5)。

接著，纖維束取向推定部155對去除了z紗線的二值圖像再次應用附帶參考的方向性距離法，推定x紗線和y紗線的纖維束取向(sp6)。

然后，纖維束取向推定部155制作纖維束取向推定圖像g100，并通過顯示部14進行顯示(sp7)，結束該圖像解析處理。

(3)各處理的詳細內容

以下，參照圖4～圖17，此外使用計算式(式1～式3)來詳細說明在圖3中說明的各處理(sp2～sp6)。此外，x射線ct圖像g10的輸入處理(sp1)和纖維束取向推定圖像g100的顯示處理(sp7)為通常的方法，因此省略說明。

(3-1)二值化處理

圖4表示二值化處理的詳細流程圖。二值化部151首先根據規(guī)定閾值對經由輸入部12輸入的x射線ct圖像g10進行二值化，生成二值圖像(sp21)。生成的二值圖像有時由于纖維束內的纖維密度不均勻而產生的黑斑。即，存在將本來應該為背景的體素當作前景體素的地方。

為了將該本來應該為背景的前景體素返回至背景體素，二值化部151執(zhí)行形態(tài)處理中的閉合處理(sp22)。通過執(zhí)行閉合處理，能夠去除二值圖像中產生的黑斑。二值化部151生成去除了黑斑從而清晰的二值圖像，結束本處理。

圖5表示通過二值化處理生成的處理圖像。具體地說，表示閉合處理前的二值圖像g21以及閉合處理后的二值圖像g22。

如圖5所示，在二值圖像g21中存在將本來應該為背景的體素當作前景體素的地方。通過執(zhí)行閉合處理，能夠得到去除了黑斑的閉合處理后的二值圖像g22。

(3-2)重疊區(qū)域提取處理

圖6表示重疊區(qū)域提取處理的詳細流程圖。重疊區(qū)域提取部152為了提取x紗線與y紗線立體交叉的重疊區(qū)域，執(zhí)行形態(tài)處理中的開放處理(sp31)。

通過執(zhí)行開放處理，能夠提取粗略的重疊區(qū)域。開放處理后的二值圖像的形狀發(fā)生變化，存在將本來應該為重疊區(qū)域的前景體素當作不是重疊區(qū)域的背景體素的地方。

為了將該本來為重疊區(qū)域的背景體素返回至前景體素，重疊區(qū)域提取部152執(zhí)行形態(tài)處理中的膨脹處理(sp32)。

接著，重疊區(qū)域提取部152計算膨脹處理后的二值圖像與開放處理前的二值圖像的交集，提取正確的重疊區(qū)域(sp33)。重疊區(qū)域提取部152生成提取了重疊區(qū)域的重疊區(qū)域提取圖像，然后結束本處理。

圖7表示通過重疊區(qū)域提取處理生成的處理圖像。具體地說，表示了開放處理前的二值圖像g3、開放處理后的二值圖像g31、膨脹處理后的二值圖像g32以及重疊區(qū)域提取圖像g33。

如圖7所示，通過計算開放處理前的二值圖像g3與膨脹處理后的二值圖像g32的交集，能夠得到提取了正確的重疊區(qū)域的重疊區(qū)域提取圖像g33。

(3-3)參考方向決定處理

圖8表示了參考方向決定處理的詳細流程圖。參考方向決定部153對重疊區(qū)域提取圖像g33(圖7)應用通常的方向性距離法，針對重疊區(qū)域中的各體素推定重疊方向(sp41)。

圖9表示應用通常的方向性距離法來推定關注體素的方向性距離的處理的概要。在輸入了二值化后的輸入圖像g411的情況下，對關注體素，計算直到背景為止的方向性距離。關于朝向與關注體素接近的體素的方向，在二維的情況下為8個方向，在三維的情況下為26個方向。

在此，說明二維的情況。當把去除了朝向為反向的向量后的方向設為vsi時，在圖像上在4個方向上行進，在到達背景的時間點停止行進。關于－vsi也同樣地在圖像上在4個方向行進，在到達背景的時間點停止行進。

此時，當把在vsi方向上行進的距離設為d(vsi)，把在－vsi方向上行進的距離設為d(－vsi)時，用以下式1表示±vsi方向的方向性距離d(±vsi)。

[式1]

d(±vsi)＝d(vsi)+d(-vsi)……(1)

另外，分別用以下的式2和3定義方向向量pi和方向張量m。在對方向張量m進行固有值分解的情況下，針對最大固有值的固有向量表示纖維束取向。另外，針對第二大固有值的固有向量表示纖維束的寬度方向。在三維的情況下，針對最小固有值的固有向量表示厚度方向。

[式2]

pi＝d(±vsi)·vsi……(2)

[式3]

返回至圖8，關于在步驟sp41中推定出的重疊方向，根據方向性距離法的性質，對于重疊區(qū)域與背景的邊界部分的體素無法正確地計算。為了去除無法正確計算的該邊界部分的體素，參考方向決定部153執(zhí)行形態(tài)處理中的收縮處理，從重疊區(qū)域中提取中心區(qū)域(sp42)。

接著，因為構成了重疊區(qū)域各自獨立的連結成分，因此參考方向決定部153對于提取出的各中心區(qū)域執(zhí)行6附近的標記處理來進行區(qū)域分割，從而將重疊區(qū)域分離(sp43)。

接著，參考方向決定部153使在步驟sp41中推定出的重疊方向中的、中心區(qū)域的各體素的重疊方向平均化，將平均化而得到的方向決定為參考方向(sp44)，結束本處理。

圖10表示通過參考方向決定處理生成的處理圖像。具體地說，表示重疊方向推定圖像g41、中心區(qū)域提取圖像g42、標記圖像g43以及參考方向圖像g44。

如圖10所示，通過對重疊區(qū)域提取圖像g33(圖7)應用通常的方向性距離法，能夠得到重疊方向推定圖像g41。另外，通過對重疊方向推定圖像g41執(zhí)行收縮處理，能夠得到中心區(qū)域提取圖像g42。

另外，通過對中心區(qū)域提取圖像g42執(zhí)行標記處理，能夠得到標記圖像g43。另外，通過使中心區(qū)域的重疊方向平均化，能夠得到參考方向圖像g44。

(3-4)z紗線去除處理

圖11表示z紗線去除處理的詳細流程圖。z紗線去除部154對在步驟sp2中得到的二值圖像g22(圖5)應用附帶參考的方向性距離法，計算方向性距離(sp51)。

附帶參考的方向性距離法是指使計算方向性距離的平面旋轉為與參考方向垂直的平面(參考平面)，在該參考平面上通過二維的通常的方向性距離法來計算方向性距離的方法。

圖12表示應用通常的方向性距離法來計算方向性距離時的概念圖。當在x紗線與y紗線立體交叉的重疊區(qū)域中應用通常的方向性距離法來計算方向性距離時，存在無法計算正確的方向性距離的問題。

具體地說，關于接近纖維束取向的方向，由于x紗線和y紗線的纖維束厚度薄，因此很快到達背景從而距離變短，另一方面，關于接近參考方向的方向，由于重疊區(qū)域的截面粗因此不容易到達背景從而距離變長。

因此，與接近纖維束取向的方向的距離相比，接近參考方向的方向的距離長。結果，存在將接近參考方向的方向的方向性距離作為重疊區(qū)域中的方向性距離來進行計算的問題。因此，在本實施方式中應用附帶參考的方向性距離法來計算方向性距離。

圖13表示附帶參考的方向性距離法的概念圖。附帶參考的方向性距離法是指使成為方向性距離計算對象的平面進行旋轉，在旋轉后的平面上應用二維的方向性距離法來計算方向性距離的方法。

具體地說，使方向±vs0旋轉從而使其與參考方向一致，并按照使±vs0進行旋轉的量，還使與±vs0垂直的其它方向±vs1和±vs2進行旋轉。當把旋轉后的方向分別設為±vs0’、±vs1’以及±vs2’時，x紗線或y紗線存在于通過±vs1’以及±vs2’定義的平面(參考平面)上。

通過在該參考平面上應用二維的方向性距離法來計算方向性距離，能夠不受重疊區(qū)域的影響地推定纖維束取向。

圖14表示應用附帶參考的方向性距離法來計算參考平面上的x紗線或y紗線的方向性距離時的概念圖。在針對位于重疊區(qū)域內的關注體素計算方向性距離的情況下，應用附帶參考的方向性距離法來計算方向性距離。

具體地說，如上所述，使方向±vs0旋轉從而使其與參考方向一致，還使與±vs0垂直的其它方向±vs1、±vs2、±vs3以及±vs4旋轉相同程度。

當把旋轉后的方向分別設為±vs0’、±vs1’、±vs2’、±vs3’以及±vs4’時，在通過±vs1’、±vs2’、±vs3’以及±vs4’定義的參考平面上應用二維的方向性距離法來計算方向性距離。

在應用附帶參考的方向性距離法來計算方向性距離的情況下，上述式1～3表示的關系成立。即，在對方向張量m進行固有值分解的情況下針對最大固有值的固有向量表示纖維束取向。另外，針對第二大固有值的固有向量表示纖維束的寬度方向。在三維的情況下針對最小固有值的固有向量表示厚度方向。

返回至圖11，z紗線去除部154使用在進行附帶參考的方向性距離法的計算時得到的針對第二大固有值的固有向量來計算附帶條件的方向性距離(sp52)。然后，z紗線去除部154根據附帶條件的方向性距離，將x紗線或y紗線的區(qū)域與z紗線的區(qū)域分離(sp53)。

附帶條件的方向性距離是指在以下的條件下計算的方向性距離，該條件為從關注體素開始在圖像上行進，在到達背景的時間點，或者在位于下一個行進的位置的體素的固有向量所示的方向與行進方向所形成的角度大于規(guī)定閾值的時間點停止行進。

圖15表示用于推定附帶條件的方向性距離的處理概要。以位于位置x的關注體素為起點，在關注體素的固有向量v2(x)所示的方向上在圖像上行進，在到達背景的時間點停止行進。將從最后到達的體素至關注體素為止的距離設為xe1。

另一方面，關于反向的向量－v2(x)的方向，在圖像上行進，在位于下一個行進的位置x’的體素的固有向量v2(x’)所示的方向與－v2(x)所示的方向所形成的角度大于規(guī)定閾值的時間點停止行進。將從最后到達的體素至關注體素為止的距離設為xe2。然后，將距離xe1與xe2的絕對值的總和作為附帶條件的方向性距離dc(x)。

圖16表示根據附帶條件的方向性距離將x紗線或y紗線的區(qū)域與z紗線的區(qū)域分離的概念圖。x紗線或y紗線的寬度wxy大于z紗線的截面的對角線距離wz，因此x紗線或y紗線的區(qū)域中的體素p的附帶條件的方向性距離dc(p)大于z紗線的區(qū)域中的體素q的附帶條件的方向性距離dc(q)。

例如，能夠將大于z紗線的截面的對角線距離wz且小于x紗線或y紗線的寬度wxy的值設為閾值，將大于該閾值的附帶條件的方向性距離的體素包含在x紗線或y紗線的區(qū)域中，將小于閾值的附帶條件的方向性距離的體素包含在z紗線的區(qū)域中。

返回至圖11，z紗線去除部154刪除z紗線的區(qū)域中包含的體素，由此從二值圖像g22(圖5)中去除z紗線(sp54)。當去除了z紗線時，有時還將x紗線或y紗線的區(qū)域的一部分去除。

因此，z紗線去除部154對去除z紗線后的圖像執(zhí)行形態(tài)處理中的膨脹處理(sp55)，計算膨脹處理后的圖像與膨脹處理前的圖像的交集來提取x紗線和y紗線(sp56)，結束本處理。

(3-5)纖維束取向推定處理

圖17表示纖維束取向推定處理的詳細流程圖。纖維束取向推定部155對提取了x紗線和y紗線的圖像再次應用附帶參考的方向性距離法，推定x紗線和y紗線的纖維方向(sp61)。

接著，纖維束取向推定部155在推定了纖維方向的體素中僅保留視為精度較好的中心部附近的體素而刪除其它體素，提取中心部附近的體素(sp62)。

接著，纖維束取向推定部155對中心部附近的體素執(zhí)行以下的聚類處理：將相同程度的方向的體素進行連接，使連接的體素群屬于同一簇，使在纖維束的途中向量急劇變化的體素屬于不同的簇。

纖維束的取向不會在數個體素單位的短距離中急劇變化。由此，能夠判斷在纖維束的途中向量急劇變化的體素為噪聲。因此，纖維束取向推定部155刪除所屬的體素少的簇來去除噪聲(sp63)，結束本處理。

(4)解析結果

以下，參照圖18～圖21來說明對各種輸入圖像執(zhí)行上述說明的本實施方式中的圖像解析處理時得到的處理圖像。

圖18表示輸入圖像為模擬圖像g11的情況下執(zhí)行圖像解析處理而得到的處理圖像。模擬圖像g11為設定x紗線和y紗線的纖維束的長度、寬度、厚度、x紗線的纖維層和y紗線的纖維層的層疊數、纖維束彼此的間隔以及z紗線的截面形狀等而制作的三維圖像。即，是各種信息為已知的圖像。

在對該模擬圖像g11執(zhí)行本實施方式的圖像解析處理的情況下，能夠得到三維圖像g110、x紗線截面圖像g111以及y紗線截面圖像g112。三維圖像g110為表示x紗線和y紗線的纖維束取向的三維圖像。

另外，x紗線截面圖像g111為表示x紗線的纖維束取向的二維圖像，y紗線截面圖像g112為表示y紗線的纖維束取向的二維圖像。

對得到這些處理圖像g110～g112時計算出的纖維束的方向與制作模擬圖像g11時設定的纖維束的方向所形成的角度進行計算來作為誤差，結果，誤差的最大值為89.9度，平均值為4.9度。

此外，已知通常方向性距離法本身的平均誤差為4.3度。另外，全部體素中94.7％的體素的誤差在6度以內。根據上述結果，能夠確認本實施方式的圖像解析處理的有效性。

圖19表示輸入圖像為x射線ct圖像g10的情況下執(zhí)行圖像解析處理而得到的處理圖像。x射線ct圖像g10為通過x射線ct裝置對實際的編織物進行拍攝而得到的三維圖像。另外，x射線ct圖像g10與圖20表示的x射線ct圖像g12相比為高分辨率。

在對該x射線ct圖像g10執(zhí)行本實施方式的圖像解析處理的情況下，能夠得到三維圖像g100(纖維束取向推定圖像g100)、x紗線截面圖像g101和y紗線截面圖像g102。三維圖像g100為表示x紗線和y紗線的纖維束取向的三維圖像。

另外，x紗線截面圖像g101為表示x紗線的纖維束取向的二維圖像，y紗線截面圖像g102為表示y紗線的纖維束取向的二維圖像。通過參照這些處理圖像g100～g102，能夠容易地確認x紗線和y紗線的纖維束取向。

另外，計算時間為243.8秒左右。當以往要對相同程度的體素數的圖像進行解析時，需要20倍左右的時間，由此通過執(zhí)行本實施方式的圖像解析處理，能夠謀求縮短計算時間。

圖20表示輸入圖像為x射線ct圖像g12的情況下執(zhí)行圖像解析處理而得到的處理圖像。x射線ct圖像g12與圖19表示的x射線ct圖像g10相比，不同點在于為低分辨率。

在對該x射線ct圖像g12執(zhí)行本實施方式的圖像解析處理的情況下，能夠得到三維圖像g120、x紗線截面圖像g121以及y紗線截面圖像g122。當參照這些處理圖像g120～g122時，雖然一部分纖維束缺失，但是能夠大致容易地確認其纖維束取向。

圖21表示輸入圖像為缺陷模擬圖像g13的情況下執(zhí)行圖像解析處理而得到的處理圖像。缺陷模擬圖像g13與圖18表示的模擬圖像g11相比，不同點在于具有x紗線的纖維束在途中彎曲10度的缺陷。

在對該缺陷模擬圖像g13執(zhí)行本實施方式的圖像解析處理的情況下，能夠得到三維圖像g130和x紗線截面圖像g131。當參照這些處理圖像g130和g131時，雖然存在一些誤差，但是能夠容易地確認x紗線的纖維束從彎曲地點開始彎曲。

圖22表示輸入圖像為曲面模擬圖像g14的情況下執(zhí)行圖像解析處理而得到的處理圖像。曲面模擬圖像g14與圖18表示的模擬圖像g11相比，不同點在于，具有在圓弧上變形的曲面部分。

在對該曲面模擬圖像g14執(zhí)行本實施方式的圖像解析處理的情況下，能夠得到三維圖像g140和x紗線截面圖像g141。當參照這些處理圖像g140和g141時，能夠容易地確認纖維束在圓弧上取向。

(5)本實施方式的效果

如上所述通過本實施方式的圖像解析裝置、圖像解析方法以及程序，針對通過三維編織制作的編織物的x射線ct圖像應用附帶參考的方向性距離法而去除z紗線，對去除了z紗線的圖像再次應用附帶參考的方向性距離法來推定x紗線和y紗線的纖維束取向，因此能夠去除z紗線的影響而高精度地在短時間內推定纖維束取向。另外，對于具有曲面形狀的編織物的x射線ct圖像，也能夠推定纖維束取向。因此，能夠用于實際的產品檢查。

接著，參照圖23～圖46說明作為本發(fā)明的一實施方式的圖像解析裝置2。

(1)圖像解析裝置2的整體結構

如圖23所示，本實施方式的圖像解析裝置2具備與圖1的圖像解析裝置1相同的cpu11、輸入部12、存儲部13、顯示部14以及存儲器15，存儲器15具有二值化部151和重疊區(qū)域提取部152。這些結構與圖1中具有相同的名稱和附圖標記的結構相同，并進行相同的處理，因此省略重復說明。在此，在本實施方式的圖像解析裝置2中，在存儲器15內還具有重疊區(qū)域形態(tài)解析部253。重疊區(qū)域形態(tài)解析部253對通過重疊區(qū)域提取部152的處理而得到的重疊區(qū)域提取圖像g33中的“重疊區(qū)域”的形態(tài)進行解析。在此，在重疊區(qū)域的形態(tài)中不僅包含一個重疊區(qū)域的形態(tài)，還包含將多個散布的重疊區(qū)域進行組合后的形態(tài)。另外，在附圖中，圖像解析裝置2具有顯示部14，但是也可以是不具有不進行顯示地進行解析的顯示部14的結構。

(2)圖像解析處理p2的流程圖

在圖24中表示通過圖像解析裝置2進行的圖像解析處理p2的流程圖。如該流程圖所示，在圖像解析處理p2中，首先，進行x射線ct圖像的輸入(sp1)、二值化處理(sp2)以及重疊區(qū)域提取處理(sp3)。這些處理與上述圖3表示的sp1～sp3的處理相同，因此省略說明。接著，在重疊區(qū)域形態(tài)解析處理(sp14)中，通過重疊區(qū)域形態(tài)解析部253，使用通過重疊區(qū)域提取處理而得到的重疊區(qū)域提取圖像g33(參照圖7)來進行重疊區(qū)域21的形態(tài)解析的處理。在此，在本實施方式中，以后，將重疊區(qū)域提取圖像g33中的各個重疊區(qū)域21稱為“支柱21”。

(3)重疊區(qū)域形態(tài)解析處理的示例

重疊區(qū)域形態(tài)解析部253例如可以通過對三維圖像的支柱21內包含的體素數進行計數等來計算各支柱21的體積。在圖25中例示地表示了在二維圖像中用黑點表示的支柱21內的像素，但是在三維圖像中也能夠同樣地對體素數進行計數。在該情況下，例如也可以將多個支柱21的體積平均作為參照值，通過與特定支柱21的體積進行比較，來判斷特定支柱21的體積是否適當。在此，在體積的參照值中，不僅能夠使用多個支柱21的平均，還能夠使用設計上的基準值這樣的其它值。另外，雖然設為計算支柱21的體積，但是在使用二維圖像進行解析的情況下也可以計算面積。通過這樣的計算體積的處理，能夠提供與周圍的支柱21的體積相比不同的體積等用于檢測纖維束取向異常的信息。另外，根據所提供的信息，例如能夠檢測應該立體交叉的x紗線或y紗線數量不足等纖維束取向異常。

另外，重疊區(qū)域形態(tài)解析部253還可以計算支柱21伸長的方向。關于支柱21伸長的方向，例如圖26所示，即可以使用形態(tài)圖像處理的收縮處理來使伸長方向顯現化(s21的斜線部)、向量化(s22)，也能夠使用其它的圖像處理來求出伸長的方向。另外，例如，重疊區(qū)域形態(tài)解析部253將多個支柱21伸長的方向的平均作為參照值，通過與特定的支柱21伸長的方向進行比較，能夠判斷支柱21伸長的方向是否適當。在此，關于伸長的方向的參照值，不僅能夠使用多個支柱21伸長的方向的平均，還能夠使用設計上的基準值等其它值。另外，也可以與支柱21伸長的方向不同地求出支柱21的中立軸的形狀。例如如圖27的斜線部所示，能夠使用形態(tài)圖像處理的收縮處理來提取中立軸。在該情況下，還能夠檢測支柱21并非向某一固定方向伸長而彎曲的情況等，并且通過與參照形狀進行比較，能夠判斷是否為適當的形狀。通過這樣的對支柱21伸長的方向、中立軸進行計算的處理，例如能夠提供層疊產生偏差的部分等其他的用于檢測纖維束取向異常的信息。另外，根據提供的信息，能夠檢測纖維束取向異常。

另外，重疊區(qū)域形態(tài)解析部253還可以計算多個支柱21的重心位置。理想地是，將相鄰支柱21的重心位置進行連結的線成為x紗線和y紗線延伸的方向，因此在xy平面上，檢測這些紗線不在平滑的線上排列等，求出不規(guī)則的重心位置，由此能夠判斷重心位置是否適當。圖28表示了用斜線表示的支柱21的重心g與周圍支柱21的重心g相比并非平滑地排列的樣子。另外，確定某一區(qū)域，通過對該區(qū)域內包含的重心g的數量進行計數等，求出在區(qū)域內是否存在適當數量的重心g，由此可以判斷是否存在不規(guī)則的配置的支柱21。通過這樣的計算重心位置的處理，能夠提供用于檢測纖維束取向異常的信息。另外，根據所提供的信息，能夠檢測纖維束取向異常。

在上述重疊區(qū)域形態(tài)解析處理的示例中，還記載了包含纖維束取向異常檢測處理的情況，但是也可以在重疊區(qū)域形態(tài)解析部253的處理中不包含這樣的檢測，僅進行計算支柱21的體積、計算支柱21伸長的方向或計算支柱21的重心位置等重疊區(qū)域形態(tài)解析處理，并且，例如通過輸出到顯示部14或向其它裝置發(fā)送計算出的數據等來提供處理后的信息。通過這種處理，能夠提供用于檢測纖維束取向異常的信息。

(4)纖維束取向異常的具體例

以下，使用圖29～46說明通過支柱21的形態(tài)解析檢測出的纖維束取向異常的具體例。在檢測纖維束取向異常的具體例中，還記載了在通過重疊區(qū)域形態(tài)解析部253進行的重疊區(qū)域形態(tài)解析處理中包含檢測纖維束取向異常的情況，但是在重疊區(qū)域形態(tài)解析部253的處理中也可以不包含這樣的判斷而僅提供通過支柱21的形態(tài)解析而得到的與纖維束取向有關的信息。參考faa(federalaviationadministration：美國聯邦航空管理局)發(fā)布的報告“dot/faa/ar-06/10”來記載了圖30、33、36、38、40、41、43、44以及46。

圖29在xy平面上表示x紗線和y紗線的纖維束的理想的配置例子，圖30和31分別示意性地表示圖29那樣的理想的纖維束的配置中的xz平面(或yz平面)以及xy平面上的重疊區(qū)域提取圖像。如圖30所示，在x紗線和y紗線的纖維束的理想的配置中，支柱21在x紗線和y紗線的纖維束的層疊方向(z方向)上不間斷地延伸，如圖31所示，支柱21(或支柱21的重心g)規(guī)則性地排列。

圖32(－y方向)表示在xy平面上纖維束的一部分在與應該延伸的方向不同的xy平面內的方向上延伸的狀態(tài)。圖33和34分別示意性地表示圖32的纖維束的狀態(tài)下的yz平面和xy平面的重疊區(qū)域提取圖像。如這些附圖所示，在圖33的重疊區(qū)域提取圖像中，由于向－y方向延伸的纖維束(斜線部)的影響，具有該纖維素的支柱211的重心位置g向－y方向移動。在圖34的重疊區(qū)域提取圖像中，支柱211的重心位置g從附近的支柱21的重心位置g的規(guī)則性偏離。因而，重疊區(qū)域形態(tài)解析部253例如能夠解析重疊區(qū)域提取圖像中的各個支柱21的重心位置g，檢測纖維束的一部分在與應該延伸的方向不同的xy平面內的方向上延伸的區(qū)域。

圖35表示在xy平面上纖維束的一部分在與應該延伸的方向不同的z方向上延伸的狀態(tài)。圖36示意性地表示圖35的纖維束的狀態(tài)下的yz平面的重疊區(qū)域提取圖像。如該圖所示，在圖36的重疊區(qū)域提取圖像中，由于向z方向延伸的纖維束(斜線部)的影響，支柱212被分割成多個。因而，重疊區(qū)域形態(tài)解析部253能夠使用重疊區(qū)域提取圖像來檢測纖維束的一部分在與應該延伸的方向不同的z方向上延伸的區(qū)域。

圖37表示在xy平面上纖維束的一部分在與應該延伸的方向不同的一同具有xy平面內的方向成分和z方向成分的方向上延伸的狀態(tài)。圖38示意性地表示圖37的纖維束的狀態(tài)下的yz平面的重疊區(qū)域提取圖像。如該圖所示，在圖38的重疊區(qū)域提取圖像中，由于在－z方向和－y方向上延伸的纖維束(斜線部)的影響，支柱213被分割成多個。另外，被分割的支柱213的一方的重心位置g向+y方向移動。因而，重疊區(qū)域形態(tài)解析部253能夠使用重疊區(qū)域提取圖像來檢測纖維束的一部分在與應該延伸的方向不同的一同具有xy平面內的方向成分和z方向成分的方向上延伸的區(qū)域。

圖39和圖40表示在xy平面上各個纖維束作為整體在傾斜的方向上延伸的狀態(tài)。圖41示意性地表示圖39或40所示的纖維束的狀態(tài)下的xy平面的重疊區(qū)域提取圖像。如該圖所示，在圖41的重疊區(qū)域提取圖像中，規(guī)則排列的支柱21的排列中的一部分列即支柱214被錯開配置。因而，重疊區(qū)域形態(tài)解析部253通過使用重疊區(qū)域提取圖像來檢測例如支柱214的重心位置g的排列從規(guī)定的規(guī)則性偏離等，能夠檢測纖維束作為整體在傾斜的方向上延伸的區(qū)域。

圖42表示在xy平面上纖維束的一部分(例如x紗線)不存在的狀態(tài)。圖43表示在xy平面上纖維束的一部分的粗細比其它纖維束粗或細的情況。圖44示意性地表示圖42或圖43的纖維束的狀態(tài)下的yz平面的重疊區(qū)域提取圖像。在圖44的重疊區(qū)域提取圖像中，一部分支柱215被分割成多個，即使將分割后的多個支柱215的體積合并，也小于通常的支柱21的體積。因而，重疊區(qū)域形態(tài)解析部253能夠使用重疊區(qū)域提取圖像來檢測纖維束的一部分不存在的情況以及纖維束的一部分的粗細比其它纖維束粗或細的區(qū)域。

圖45表示在xy平面上纖維束的一部分(例如x紗線)折返的狀態(tài)。圖46示意性地表示圖45的纖維束的狀態(tài)下的xy平面的重疊區(qū)域提取圖像。在圖46的重疊區(qū)域提取圖像中，由于纖維素的一部分折返，折返周圍的支柱216被推開，與其周圍的支柱21的有規(guī)則性的重心位置g相比，支柱216的重心位置g發(fā)生混亂。因而，重疊區(qū)域形態(tài)解析部253能夠使用重疊區(qū)域提取圖像來對這些重心位置g的混亂進行解析，由此檢測纖維束的一部分折返的區(qū)域。

(5)本實施方式的效果

如上所述，在本實施方式的圖像解析裝置2中，二值化部151對由x紗線、y紗線以及z紗線的纖維束構成的編織物的三維圖像進行二值化，重疊區(qū)域提取部152從二值化后的圖像中提取x紗線和y紗線立體交叉的重疊區(qū)域，重疊區(qū)域形態(tài)解析部253對提取出的重疊區(qū)域的形態(tài)進行解析。由此，僅通過處理時間短的圖像處理的組合來進行解析，因此能夠更簡單且在短時間內提供與x紗線和y紗線的纖維束取向有關的信息。

另外，重疊區(qū)域形態(tài)解析部253在進一步檢測纖維束取向異常的情況下，能夠更簡單且在短時間內發(fā)現纖維束取向異常的部位。另外，即使在對于纖維束取向異常需要進行詳細解析的情況下，僅對包含檢測出的取向異常的區(qū)域另行進行詳細解析即可，因此作為整體能夠謀求縮短解析時間。

另外，無論全體纖維束的形狀為平面的情況還是為曲面的情況，都能夠使用本實施方式的針對重疊區(qū)域的形態(tài)的解析。無論全體纖維束為怎樣的形狀都能夠提取重疊區(qū)域，因此能夠應用重疊區(qū)域的形態(tài)的解析處理。另外，支柱21沿全體纖維束的形狀有規(guī)則地配置，因此無論全體纖維束的形狀為平面的情況還是為曲面的情況，都能夠通過解析重疊區(qū)域的形態(tài)來檢測纖維束取向的異常等。

另外，重疊區(qū)域形態(tài)解析部253可以計算重疊區(qū)域的體積。在該情況下，能夠與重疊區(qū)域的“體積的參照值”進行比較，并且能夠將該“體積的參照值”設為多個重疊區(qū)域的體積的平均值。

另外，重疊區(qū)域形態(tài)解析部253可以計算重疊區(qū)域伸長的方向。在該情況下，能夠與重疊區(qū)域的“朝向的參照值”進行比較，能夠將該“朝向的參照值”設為多個上述重疊區(qū)域的方向的平均值。

另外，重疊區(qū)域形態(tài)解析部253可以計算多個重疊區(qū)域的重心位置。在該情況下，能夠檢測重心位置的配置不規(guī)則的區(qū)域。例如求出重心位置排列的參照線，計算從該參照線以怎樣的程度分離或計算在某一區(qū)域內包含的重心位置的數量等，由此可以檢測重心位置的配置不規(guī)則的區(qū)域。

另外，重疊區(qū)域形態(tài)解析部253可以計算重疊區(qū)域的中立軸，并與參照形狀進行比較。在該情況下，在中立軸的計算中能夠使用形態(tài)的收縮處理。

另外，本實施方式的圖像解析方法的特征在于，具備以下步驟：對由x紗線、y紗線以及z紗線的纖維束構成的編織物的三維圖像進行二值化；從上述二值化后的圖像中提取上述x紗線和上述y紗線立體交叉的重疊區(qū)域；以及對上述提取出的重疊區(qū)域的形態(tài)進行解析。

另外，本實施方式的程序使計算機執(zhí)行以下步驟：對由x紗線、y紗線以及z紗線的纖維束構成的編織物的三維圖像進行二值化；從上述二值化后的圖像中提取上述x紗線和上述y紗線立體交叉的重疊區(qū)域；以及對上述提取出的重疊區(qū)域的形態(tài)進行解析。

(6)其它的實施方式

在上述各實施方式中，說明了分別使用圖1所示的圖像解析裝置1的圖像解析處理(圖3)以及圖23所示的圖像解析裝置2的圖像解析處理p2(圖24)的情況，但是能夠使這些處理關聯地進行動作。在圖47中表示了將這些處理進行組合來執(zhí)行的一實施方式的圖像解析裝置3的結構。圖像解析裝置3與上述圖像解析裝置1或2相同，具備cpu11、輸入部12、存儲部13、顯示部14以及存儲器15，存儲器15具有二值化部151、重疊區(qū)域提取部152、參考方向決定部153、z紗線去除部154、纖維束取向推定部155以及重疊區(qū)域形態(tài)解析部253。各結構進行與上述對應的結構的處理相同的處理。

在圖48中表示了該圖像解析裝置3的圖像解析處理p3的流程圖。圖像解析處理p3首先與圖24的圖像解析處理p2同樣地進行x射線ct圖像的輸入(sp1)、二值化處理(sp2)、重疊區(qū)域提取處理(sp3)以及重疊區(qū)域形態(tài)解析處理(sp14)。在此，在重疊區(qū)域形態(tài)解析處理中，不僅解析重疊區(qū)域的形態(tài)，還進行纖維束取向異常檢測。接著，對于包含檢測出的纖維束取向異常的區(qū)域，與圖3的圖像解析處理p1中的對應的處理相同地，進行參考方向決定處理(sp4)、z紗線去除處理(sp5)、纖維束取向推定處理(sp6)。

這樣，將處理時間短的圖像解析處理p2應用于例如成為檢查對象的x射線ct圖像，將處理時間較長且進行詳細解析的圖像解析處理p1應用于通過圖像解析處理p2檢測出的包含纖維束取向異常的區(qū)域，由此能夠更高效且詳細地解析纖維束取向異常。

符號說明

1：圖像解析裝置；11：cpu；12：輸入部；13：存儲部；14：顯示部；15：存儲器；151：二值化部；152：重疊區(qū)域提取部；153：參考方向決定部；154：z紗線去除部；155：纖維束取向推定部；2：圖像解析裝置；21：支柱；253：重疊區(qū)域形態(tài)解析部；3：圖像解析裝置。