CN106558086B - 一种林木应力波断层成像的速度误差校正方法 - Google Patents
一种林木应力波断层成像的速度误差校正方法 Download PDFInfo
- Publication number
- CN106558086B CN106558086B CN201610908293.2A CN201610908293A CN106558086B CN 106558086 B CN106558086 B CN 106558086B CN 201610908293 A CN201610908293 A CN 201610908293A CN 106558086 B CN106558086 B CN 106558086B
- Authority
- CN
- China
- Prior art keywords
- grid cell
- simulation drawing
- wood cross
- section
- section simulation
- 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.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims abstract description 36
- 238000012937 correction Methods 0.000 title claims abstract description 32
- 238000003384 imaging method Methods 0.000 title claims abstract description 32
- 239000002023 wood Substances 0.000 claims abstract description 140
- 238000004088 simulation Methods 0.000 claims abstract description 111
- 230000007246 mechanism Effects 0.000 claims abstract description 9
- 230000011218 segmentation Effects 0.000 claims description 9
- 239000003086 colorant Substances 0.000 claims description 5
- 230000007547 defect Effects 0.000 abstract description 12
- 230000015572 biosynthetic process Effects 0.000 abstract description 2
- 238000009659 non-destructive testing Methods 0.000 description 7
- 239000000463 material Substances 0.000 description 5
- 238000000354 decomposition reaction Methods 0.000 description 4
- 230000008901 benefit Effects 0.000 description 3
- 238000002474 experimental method Methods 0.000 description 3
- 230000002950 deficient Effects 0.000 description 2
- 239000011159 matrix material Substances 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 241001444167 Picus Species 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000018109 developmental process Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000010079 rubber tapping Methods 0.000 description 1
- 230000007480 spreading Effects 0.000 description 1
- 238000003325 tomography Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T11/00—2D [Two Dimensional] image generation
- G06T11/003—Reconstruction from projections, e.g. tomography
- G06T11/005—Specific pre-processing for tomographic reconstruction, e.g. calibration, source positioning, rebinning, scatter correction, retrospective gating
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T11/00—2D [Two Dimensional] image generation
- G06T11/40—Filling a planar surface by adding surface attributes, e.g. colour or texture
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种林木应力波断层成像的速度误差校正方法，包括以下步骤：生成木材横截面模拟图，对木材横截面模拟图进行网格划分；计算木材横截面模拟图中网格单元的传播速度；计算木材横截面模拟图中网格单元的传播速度相对误差值；根据投票机制对木材横截面模拟图中非边缘网格单元的传播速度相对误差值进行校正；根据校正后的木材横截面模拟图中网格单元的传播速度相对误差值对该木材横截面模拟图中网格单元进行填色，形成该木材横截面的断层成像。本发明方法在传统的应力波速度反演成像算法基础上，利用速度误差校正机制对反演结果进行优化，形成的木材横截面断层成像准确率高，可以准确检测出木材内部是否存在缺陷以及识别缺陷的面积及形状。
Description
技术领域
本发明属于木材无损检测技术领域，具体是涉及一种林木应力波断层成像的速度误差校正方法。
背景技术
面对木材资源匮乏的现状，如何提高木材利用率，促进我国林业现代化的可持续发展，已成为林业工作者关注的焦点。木材无损检测对于发现林木内部缺陷，防止腐朽扩散，提高木材利用率具有重要意义。应力波无损检测技术具有低成本、操作方便、无辐射，不受被测木材的大小形状影响等优点，已经在林木无损检测中得到了广泛的应用。虽然现在国际上主流的林木应力波无损检测产品如德国的ARBOTOM和PICUS以及匈牙利的FAKOPP都提供了二维或三维断层成像功能，但成像效果和检测精度仍然有待提高。
发明内容
发明目的：为了克服现有技术中存在的不足，本发明提供一种林木应力波断层成像的速度误差校正方法，利用速度误差校正机制对速度反演结果进行优化，利用该方法可以准确检测出木材内部是否存在缺陷以及识别缺陷的面积及形状，与其他方法相比，本发明方法准确率高。
技术方案：为实现上述目的，本发明的一种林木应力波断层成像的速度误差校正方法，包括以下步骤：
S1生成木材横截面模拟图，对该木材横截面模拟图进行网格划分；
S2计算所述木材横截面模拟图中网格单元的传播速度；
S3计算所述木材横截面模拟图中网格单元的传播速度相对误差值；
S4根据投票机制对所述木材横截面模拟图中非边缘网格单元的传播速度相对误差值进行校正；
S5根据校正后的所述木材横截面模拟图中网格单元的传播速度相对误差值对该木材横截面模拟图中网格单元进行填色，形成该木材横截面的断层成像。
进一步地，对所述木材横截面模拟图进行网格划分基准是以该木材横截面模拟图轮廓为内切圆，对该内切圆的外切正四边形进行网格划分，设该外切正四边形的边长为a，则划分成a*a个网格单元，设所述木材横截面中网格单元的个数为m。
进一步地，所述步骤S2包括以下步骤：
S21利用应力波时间测定仪测量应力波在木材内部的传播时间tj，tj表示第j条应力波传播路径的传播时间，j＝(1,2,…,n)；
S22计算应力波传播路径在所述木材横截面模拟图中网格单元上的分割距离dij，i＝(1,2,…,m)，dij表示第j条应力波传播路径在木材横截面模拟图中第i个网格单元上的分割距离，j＝(1,2,…,n)；
S23通过方程组(一)计算所述木材横截面模拟图中网格单元的慢度pi，pi表示木材横截面模拟图中第i个网格单元的慢度，i＝(1,2,…,m)：
S24通过公式(一)计算所述木材横截面模拟图中网格单元的传播速度：
其中，vi表示木材横截面模拟图中第i个网格单元的传播速度，i＝(1,2,…,m)。
进一步地，所述步骤S22中，基于零点法求交点坐标的方法确定应力波传播路径与所述木材横截面模拟图中网格单元的交点，具体步骤如下：确定应力波传播路径直线的方程式；确定所述木材横截面模拟图中网格线的方程式；将应力波传播路径直线的方程式与所述木材横截面模拟图中网格线的方程式作差构成新的方程式；求解该新的方程式的零点，该零点即为应力波传播路径与所述木材横截面模拟图中网格线的交点，该应力波传播路径对应的发射传感器和接收传感器的位置坐标即为该应力波传播路径与所述木材横截面模拟图边界的两个交点；以该应力波传播路径对应的发射传感器的位置为起点，以该应力波传播路径对应的接收传感器的位置为终点，将该应力波传播路径与所述木材横截面模拟图中网格线的交点按照该交点与起点之间的距离排序，计算排序好的相邻两个交点之间的欧式距离，该欧式距离即为该应力波传播路径在所述木材横截面模拟图中网格单元上的分割距离。
进一步地，所述步骤S3中，根据公式(二)计算所述木材横截面模拟图中网格单元的传播速度相对误差值：
其中，δi表示所述木材横截面模拟图中第i个网格单元的传播速度相对误差值，vi表示木材横截面模拟图中第i个网格单元的传播速度，vref表示预设参考速度， i＝(1,2,…,m)。
进一步地，所述步骤S4包括以下步骤：
S41遍历校正所述木材横截面模拟图中非边缘网格单元的传播速度相对误差值；
S42判断非边缘网格单元的上、下、左、右四个邻居网格单元的传播速度相对误差值，如果四个邻居网格单元中有三个或三个以上的传播速度相对误差值大于或等于10％，则将该非边缘网格单元的传播速度相对误差值校正为10％；如果四个邻居网格单元中有三个或三个以上的传播速度相对误差值小于10％，则将该非边缘网格单元的传播速度相对误差值校正为[0,10％)的随机数；
S43经过一次遍历校正后，判断所述木材横截面模拟图中网格单元的传播速度相对误差值是否发生变化，若发生变化则进入S41，否则进入S5。
进一步地，所述步骤S5中，根据遍历校正后的所述木材横截面模拟图中网格单元的传播速度相对误差值对所述木材横截面模拟图中网格单元进行填色，若所述木材横截面模拟图中网格单元的传播速度相对误差值小于10％，则将对应的网格单元填充为颜色一；若所述木材横截面模拟图中网格单元的传播速度相对误差值大于或等于10％，则将对应的网格单元填充为颜色二，形成该木材横截面的断层成像。
有益效果：本发明与现有技术比较，具有的优点是：本发明方法首先具体是通过LSQR(Least Square QR-factorization)最小平方QR分解法求解速度反演方程组，通过反演计算得到木材横截面上各网格单元的应力波传播速度，然后求出各网格单元的传播速度相对误差值，接着通过速度误差校正机制对传播速度相对误差值进行误差校正，最后根据校正后的网格单元的传播速度相对误差值对木材横截面进行填色，形成该木材横截面的断层成像；本发明方法在传统的应力波速度反演成像算法基础上，利用速度误差校正机制对反演结果进行优化，形成的木材横截面断层成像准确率高，可以准确检测出木材内部是否存在缺陷以及识别缺陷的面积及形状。
附图说明
图1是本发明方法流程图。
图2是应力波木材无损检测实验平台。
图3是采用本发明方法对被测木材进行试验得到的木材横截面的断层成像。
图4是图3中采用的被测木材的同一横截面实际缺陷情况图。
具体实施方式
下面结合附图对本发明作更进一步的说明。
本发明提出了一种林木应力波断层成像的速度误差校正方法，参照图1，包括以下步骤，S1生成木材横截面模拟图，对该木材横截面模拟图进行网格划分；
S2计算所述木材横截面模拟图中网格单元的传播速度；
S3计算所述木材横截面模拟图中网格单元的传播速度相对误差值；
S4根据投票机制对所述木材横截面模拟图中非边缘网格单元的传播速度相对误差值进行校正；
S5根据校正后的所述木材横截面模拟图中网格单元的传播速度相对误差值对该木材横截面模拟图中网格单元进行填色，形成该木材横截面的断层成像。
实施例：
在2016年度江苏省普通高效专业学位研究生实践创新计划项目-林木内部应力波传播仿真及成像算法研究(项目编号为SJLX16_0499)中，通过实验来实现本实施例，本实施例中，采用12个传感器，首先生成木材横截面模拟图，对该模拟图进行网格划分，对木材横截面模拟图进行网格划分基准是以该木材横截面模拟图轮廓为内切圆，对该内切圆的外切正四边形进行网格划分，设该外切正四边形的边长为a，则划分成a*a个网格单元，设所述木材横截面中网格单元的个数为m；如图3所示，图3中在该木材横截面轮廓上等分12等分，将该木材横截面轮廓模拟成正圆(模拟后的正圆上等分点显示不均匀)，将该正圆的外切正四边形进行网格划分，设该外切正四边形的边长为17，则划分成17*17＝289个网格单元，设所述木材横截面中网格单元的个数m小于289个；
参照图2，采用应力波木材无损检测实验平台测试得到应力波在木材内部的传播时间，图2中包括被测木材试样1、固定在木材试样1某个横截面上的一组传感器、用于敲击传感器的脉冲锤2，与传感器连接的应力波时间测试仪3以及与应力波测试仪连接的PC电脑4，PC电脑上设有分析软件系统；首先选择木材试样1的某个横截面，测量其半径及周长，根据其半径及周长将横截面轮廓12等分，分别在12等分点处(等分点数量根据具体实验需要而定)固定传感器，规定1号传感器为横截面水平正方向的传感器，逆时针依次顺序编号1～12号传感器，分别将12个传感器均匀固定在木材试样1横截面上，使用脉冲锤2从1号传感器开始敲击，每个传感器被敲击数次，记录PC电脑显示的度数并取其平均值作为最终数据，得到应力波在木材内部的传播时间，如表1所示，表1中显示了各发送端传感器和各接收端传感器之间的应力波传播时间，得到144 条应力波传播路径的传播时间，其中存在12个时间是发送端传感器和接收端传感器为同一个传感器时的时间，时间为0，将其删除，还剩132条应力波传播路径的传播时间，可以标记为t1～t132，得到时间列向量为n＝132；
表1应力波传播时间矩阵
计算应力波传播路径在所述木材横截面模拟图中网格单元上的分割距离dij， i＝(1,2,…,m)，dij表示第j条应力波传播路径在木材横截面模拟图中第i个网格单元上的分割距离，j＝(1,2,…,n)，基于零点法求交点坐标的方法确定应力波传播路径与所述木材横截面模拟图中网格单元的交点，具体步骤如下：确定应力波传播路径直线的方程式；确定所述木材横截面模拟图中网格线的方程式；将应力波传播路径直线的方程式与所述木材横截面模拟图中网格线的方程式作差构成新的方程式；求解该新的方程式的零点，该零点即为应力波传播路径与所述木材横截面模拟图中网格线的交点，该应力波传播路径对应的发射传感器和接收传感器的位置坐标即为该应力波传播路径与木材横截面模拟图边界的两个交点；以该应力波传播路径对应的发射传感器的位置为起点，以该应力波传播路径对应的接收传感器的位置为终点，将该应力波传播路径与木材横截面模拟图中网格线的交点按照该交点与起点之间的距离排序，计算排序好的相邻两个交点之间的欧式距离，该欧式距离即为该应力波传播路径在木材横截面模拟图中网格单元上的分割距离；上述132条应力波传播路径经过木材横截面模拟图中各网格(有部分网格不经过)，最终计算得到应力波传播路径在木材横截面模拟图中网格单元上的分割距离dij， dij表示第j条应力波传播路径在木材横截面模拟图中第i个网格单元上的分割距离， i＝(1,2,…,m)，j＝(1,2,…,n)，得到距离矩阵本实施例中取n＝132；
通过列方程组(一)计算木材横截面模拟图中网格单元的慢度pi，pi表示木材横截面模拟图中第i个网格单元的慢度，i＝(1,2,…,m)：
通过LSQR(Least Square QR-factorization)最小平方QR分解法求解速度反演方程组，得到木材横截面模拟图中m个网格单元的慢度p1～pm；设方程组(一)为T＝DP，n＝132，m为木材横截面模拟图中网格单元的个数，LSQR(Least Square QR-factorization)最小平方QR分解法中对于方程组AX＝B，调用函数X＝lsqr(A，B)，就可以求得X，因此对于T＝DP方程组(一)通过最小平方QR分解法可以求得木材横截面模拟图中m个网格单元的慢度p1～pm；
通过公式(一)计算木材横截面模拟图中网格单元的传播速度：
其中，vi表示木材横截面模拟图中第i个网格单元的传播速度，i＝(1,2,…,m)，得到木材横截面模拟图中m个网格单元的慢度v1～vm；
计算木材横截面模拟图中网格单元的传播速度相对误差值：根据公式(二)计算木材横截面模拟图中网格单元的传播速度相对误差值：
其中，δi表示所述木材横截面模拟图中第i个网格单元的传播速度相对误差值，vi表示木材横截面模拟图中第i个网格单元的传播速度，vref表示预设参考速度， i＝(1,2,…,m)，最终得到木材横截面模拟图中m个网格单元的传播速度相对误差值δ1～δm；
根据投票机制对木材横截面模拟图中非边缘网格单元的传播速度相对误差值进行校正：遍历校正木材横截面模拟图中所有非边缘网格单元的传播速度相对误差值，判断非边缘网格单元的上、下、左、右四个邻居网格单元的传播速度相对误差值，如果四个邻居网格单元中有三个或三个以上的传播速度相对误差值大于或等于10％，则将该非边缘网格单元的传播速度相对误差值校正为10％；如果四个邻居网格单元中有三个或三个以上的传播速度相对误差值小于10％，则将该非边缘网格单元的传播速度相对误差值校正为[0,10％)的随机数；经过一次遍历校正后，判断木材横截面模拟图中网格单元的传播速度相对误差值是否发生变化，若发生变化则再次进行遍历校正，直至相对误差值不再发生变化；
最后根据校正后的木材横截面模拟图中网格单元的传播速度相对误差值对该木材横截面模拟图中网格单元进行填色，形成该木材横截面的断层成像，若木材横截面模拟图中网格单元的传播速度相对误差值小于10％，则将对应的网格单元填充为颜色一；若木材横截面模拟图中网格单元的传播速度相对误差值大于或等于10％，则将对应的网格单元填充为颜色二，形成该木材横截面的断层成像，本实施例中将颜色一设置为成淡灰色，将颜色二设置为深灰色，得到如图3所示的木材横截面的断层成像，在图3中可以清楚显示木材横截面上的缺陷位置和大小，将通过本发明方法得到的木材横截面的断层成像(图3所示)与采用的被测木材的同一横截面实际缺陷情况图作对比，被测木材的同一横截面实际缺陷情况图如图4所示，可以发现本发明方法得到的缺陷位置和大小和实际缺陷位置和大小基本一致，通过本发明方法可以准确检测出木材内部是否存在缺陷以及识别缺陷的面积及形状。
以上所述仅是本发明的优选实施方式，应当指出：对于本技术领域的技术人员来说，在不脱离本发明原理的前提下，还可以做出若干改进和润饰，这些改进和润饰也应视为本发明的保护范围。
Claims (7)
1.一种林木应力波断层成像的速度误差校正方法，其特征在于：包括以下步骤：
S1 生成木材横截面模拟图，对该木材横截面模拟图进行网格划分；
S2 计算所述木材横截面模拟图中网格单元的传播速度；
S3 计算所述木材横截面模拟图中网格单元的传播速度相对误差值；
S4 根据投票机制对所述木材横截面模拟图中非边缘网格单元的传播速度相对误差值进行校正；
S5 根据校正后的所述木材横截面模拟图中网格单元的传播速度相对误差值对该木材横截面模拟图中网格单元进行填色，形成该木材横截面的断层成像。
2.根据权利要求1所述的林木应力波断层成像的速度误差校正方法，其特征在于：对所述木材横截面模拟图进行网格划分基准是以该木材横截面模拟图轮廓为内切圆，对该内切圆的外切正四边形进行网格划分，设该外切正四边形的边长为a，则划分成a*a个网格单元，设所述木材横截面中网格单元的个数为m。
3.根据权利要求1所述的林木应力波断层成像的速度误差校正方法，其特征在于：所述步骤S2包括以下步骤：
S21 利用应力波时间测定仪测量应力波在木材内部的传播时间tj，tj表示第j条应力波传播路径的传播时间，j＝(1,2,…,n)；
S22 计算应力波传播路径在所述木材横截面模拟图中网格单元上的分割距离dij，dij表示第j条应力波传播路径在木材横截面模拟图中第i个网格单元上的分割距离，i＝(1,2,…,m)，j＝(1,2,…,n)；
S23 通过方程组(一)计算所述木材横截面模拟图中网格单元的慢度pi，pi表示木材横截面模拟图中第i个网格单元的慢度，i＝(1,2,…,m)：
S24 通过公式(一)计算所述木材横截面模拟图中网格单元的传播速度：
其中，vi表示木材横截面模拟图中第i个网格单元的传播速度，i＝(1,2,…,m)。
4.根据权利要求3所述的林木应力波断层成像的速度误差校正方法，其特征在于：所述步骤S22中，基于零点法求交点坐标的方法确定应力波传播路径与所述木材横截面模拟图中网格单元的交点，具体步骤如下：确定应力波传播路径直线的方程式；确定所述木材横截面模拟图中网格线的方程式；将应力波传播路径直线的方程式与所述木材横截面模拟图中网格线的方程式作差构成新的方程式；求解该新的方程式的零点，该零点即为应力波传播路径与所述木材横截面模拟图中网格线的交点，该应力波传播路径对应的发射传感器和接收传感器的位置坐标即为该应力波传播路径与所述木材横截面模拟图边界的两个交点；以该应力波传播路径对应的发射传感器的位置为起点，以该应力波传播路径对应的接收传感器的位置为终点，将该应力波传播路径与所述木材横截面模拟图中网格线的交点按照该交点与起点之间的距离排序，计算排序好的相邻两个交点之间的欧式距离，该欧式距离即为该应力波传播路径在所述木材横截面模拟图中网格单元上的分割距离。
5.根据权利要求1所述的林木应力波断层成像的速度误差校正方法，其特征在于：所述步骤S3中，根据公式(二)计算所述木材横截面模拟图中网格单元的传播速度相对误差值：
其中，δi表示所述木材横截面模拟图中第i个网格单元的传播速度相对误差值，vi表示木材横截面模拟图中第i个网格单元的传播速度，vref表示预设参考速度，i＝(1,2,…,m)。
6.根据权利要求1所述的林木应力波断层成像的速度误差校正方法，其特征在于：所述步骤S4包括以下步骤：
S41 遍历校正所述木材横截面模拟图中非边缘网格单元的传播速度相对误差值；
S42 判断非边缘网格单元的上、下、左、右四个邻居网格单元的传播速度相对误差值，如果四个邻居网格单元中有三个或三个以上的传播速度相对误差值大于或等于10％，则将该非边缘网格单元的传播速度相对误差值校正为10％；如果四个邻居网格单元中有三个或三个以上的传播速度相对误差值小于10％，则将该非边缘网格单元的传播速度相对误差值校正为[0,10％)的随机数；
S43 经过一次遍历校正后，判断所述木材横截面模拟图中网格单元的传播速度相对误差值是否发生变化，若发生变化则进入S41，否则进入S5。
7.根据权利要求1所述的林木应力波断层成像的速度误差校正方法，其特征在于：所述步骤S5中，根据遍历校正后的所述木材横截面模拟图中网格单元的传播速度相对误差值对所述木材横截面模拟图中网格单元进行填色，若所述木材横截面模拟图中网格单元的传播速度相对误差值小于10％，则将对应的网格单元填充为颜色一；若所述木材横截面模拟图中网格单元的传播速度相对误差值大于或等于10％，则将对应的网格单元填充为颜色二，形成该木材横截面的断层成像。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610908293.2A CN106558086B (zh) | 2016-10-18 | 2016-10-18 | 一种林木应力波断层成像的速度误差校正方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610908293.2A CN106558086B (zh) | 2016-10-18 | 2016-10-18 | 一种林木应力波断层成像的速度误差校正方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106558086A CN106558086A (zh) | 2017-04-05 |
| CN106558086B true CN106558086B (zh) | 2019-08-06 |
Family
ID=58443126
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610908293.2A Expired - Fee Related CN106558086B (zh) | 2016-10-18 | 2016-10-18 | 一种林木应力波断层成像的速度误差校正方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106558086B (zh) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107655978B (zh) * | 2017-08-14 | 2020-06-23 | 浙江农林大学 | 基于速度修正插值法的木材径切面缺陷成像方法 |
| CN110196281A (zh) * | 2019-06-25 | 2019-09-03 | 石家庄铁道大学 | 一种信号波传输波速的检测方法、系统及终端设备 |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP0614309B1 (en) * | 1993-03-03 | 2001-11-07 | Seiko Epson Corporation | Image processing apparatus |
| CN103440614A (zh) * | 2013-02-05 | 2013-12-11 | 浙江工业大学 | 基于虚拟图形水印的三维网格数据保护方法 |
| CN104036475A (zh) * | 2013-07-22 | 2014-09-10 | 成都智慧星球科技有限公司 | 适应于任意投影机群及投影屏幕的高鲁棒性几何校正方法 |
-
2016
- 2016-10-18 CN CN201610908293.2A patent/CN106558086B/zh not_active Expired - Fee Related
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP0614309B1 (en) * | 1993-03-03 | 2001-11-07 | Seiko Epson Corporation | Image processing apparatus |
| CN103440614A (zh) * | 2013-02-05 | 2013-12-11 | 浙江工业大学 | 基于虚拟图形水印的三维网格数据保护方法 |
| CN104036475A (zh) * | 2013-07-22 | 2014-09-10 | 成都智慧星球科技有限公司 | 适应于任意投影机群及投影屏幕的高鲁棒性几何校正方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106558086A (zh) | 2017-04-05 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105528288B (zh) | 一种软件测试方法以及装置 | |
| CN101424520B (zh) | 物件曲面的局部轮廓检测方法 | |
| CN103353923B (zh) | 基于空间特征分析的自适应空间插值方法及其系统 | |
| CN108153978B (zh) | 一种基于水动力水质模型的污染源普查结果达标评估方法 | |
| CN110398320A (zh) | 一种易于持续优化的燃气泄漏检测定位方法及系统 | |
| CN103970994B (zh) | 基于数据同化的内陆水体叶绿素浓度多模型协同反演方法 | |
| CN107780923B (zh) | 一种基于泥质校正的含水饱和度模型的建立、仿真方法 | |
| CN110242865A (zh) | 一种易于持续优化的燃气泄漏检测判定方法及系统 | |
| CN103513270A (zh) | 一种基于岩石声学特性的气层识别评价方法及装置 | |
| CN101750030A (zh) | 曲面检测系统及方法 | |
| CN101231167A (zh) | 海洋测线网系统误差的检测与调整方法 | |
| CN106558086B (zh) | 一种林木应力波断层成像的速度误差校正方法 | |
| CN110414415A (zh) | 面向课堂场景的人体行为识别方法 | |
| CN109900789B (zh) | 一种树木纵截面内部缺陷成像方法 | |
| CN109163674B (zh) | 一种面结构光自动化三维测量中传感器测量视点规划方法 | |
| CN112946751B (zh) | 一种预测不同储层平面分布的方法及系统 | |
| CN114492984A (zh) | 粉尘浓度的时空分布预测方法、装置、设备和存储介质 | |
| CN114034334A (zh) | 岩溶管道污染源和流量的识别方法 | |
| CN111122813B (zh) | 基于区域地下水流场方向的水质类别评价方法 | |
| CN107402257A (zh) | 基于路径分组插值法的木材径切面缺陷成像方法 | |
| CN110688767B (zh) | 一种评价岩体裂隙网络模型精度的综合差异度的方法 | |
| CN116932679A (zh) | 一种基于gis技术多模式行业气象服务检验应用方法 | |
| CN112259172B (zh) | 一种存在空气污染物化学反应情况的污染源逆向辨识方法 | |
| CN112257354B (zh) | 一种动态风场条件下的空气污染源逆向定位方法 | |
| CN103175749B (zh) | 一种汽车塑料内外饰件耐刮擦性能的测试方法 |
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 | ||
| CF01 | Termination of patent right due to non-payment of annual fee | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20190806 |