CN109300089B - 一种基于二分法的测井图像均衡处理的方法 - Google Patents

一种基于二分法的测井图像均衡处理的方法 Download PDF

Info

Publication number
CN109300089B
CN109300089B CN201810948463.9A CN201810948463A CN109300089B CN 109300089 B CN109300089 B CN 109300089B CN 201810948463 A CN201810948463 A CN 201810948463A CN 109300089 B CN109300089 B CN 109300089B
Authority
CN
China
Prior art keywords
autocutp
data
img
partmp
range
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.)
Active
Application number
CN201810948463.9A
Other languages
English (en)
Other versions
CN109300089A (zh
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.)
China National Petroleum Corp
China Petroleum Logging Co Ltd
Original Assignee
China National Petroleum Corp
China Petroleum Logging Co Ltd
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 China National Petroleum Corp, China Petroleum Logging Co Ltd filed Critical China National Petroleum Corp
Priority to CN201810948463.9A priority Critical patent/CN109300089B/zh
Publication of CN109300089A publication Critical patent/CN109300089A/zh
Application granted granted Critical
Publication of CN109300089B publication Critical patent/CN109300089B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T5/00Image enhancement or restoration
    • G06T5/90Dynamic range modification of images or parts thereof
    • G06T5/92Dynamic range modification of images or parts thereof based on global image properties
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06TIMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T7/00Image analysis
    • G06T7/10Segmentation; Edge detection

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Image Analysis (AREA)
  • Complex Calculations (AREA)

Abstract

本发明公开了一种基于二分法的测井图像均衡处理的方法,首先进行参数设定;然后对测井图像进行初始化分割,再进行自动分割;判断自动分割结束,根据自动分割好的数据范围进行数据值的映射,输出映射结果。本发明能够对数据值的分布范围不均匀即分割范围的数据个数所占总数百分比大于1/256的分割范围,利用2分法,进行无限次的自动分割。以求达到显示效果均匀分布即所有的显示灰度尽可能的得到有效利用,从而显著提高显示的分辨率和对比度。本发明处理流程操作简便,应用到测井图像的数据处理中,验证其具有较好的处理效果。

Description

一种基于二分法的测井图像均衡处理的方法
技术领域
本发明属于石油开发和工程领域,是一种测井图像均衡处理改进处理的方法。
背景技术
随着测井技术的发展,成像仪器在裸眼井、套管井都有广泛的应用,在井旁构造缝洞识别等有着重要的应用,套管井中在套损检测也有着一定的应用。但是现有的传统动态均衡处理方法,存在一些问题,不能够自动分割,使一些显示像素灰度处于“闲置状态”,导致处理结果的对比度和分辨率不理想。
发明内容
本发明的目的是针对目前传统图像动态均衡处理存在的“数据值的分布范围不均匀”,导致处理结果对比度和分辨率不理想问题,对传统图像动态均衡处理进行了进一步的改进。本发明处理流程具备创新性,操作简便,已经应用到了测井图像的数据处理中,得到了较好的处理效果。
本发明的目的是通过下述技术方案来实现的。
一种基于二分法的测井图像均衡处理的方法,包括以下步骤:
步骤1,首先进行参数设定:
通过标示Flag的两种方式设定切割范围,当Flag=0时,程序自动统计计算最大值、最小值;当Flag=1时,通过人工手动输入方式设定切割范围;
步骤2,对测井图像进行初始化分割:
根据灰度图像显示级别256级,将测井图像值的范围初始化分割为256份,统计每个分割范围内的数据点个数,及相应的所占数据总个数的百分比;
步骤3,进行自动分割:
利用二分法对步骤2分割为256份测井图像值分布不均匀的数据范围进行进一步分割,对于分割后的不同范围,分别统计分割范围点的数据点个数,和所占数据总个数的百分比;
步骤4判断自动分割结束:
当所有的分割份数的所占总数百分比全部小于平均1/256,则程序自动分割结束;如果不满足,则重复步骤3,进行无限次的自动分割,则根据切分范围的总个数大于设定值时,自动分割强制退出;
步骤5,根据自动分割好的数据范围进行数据值的映射,输出映射结果。
对于上述技术方案,本发明还有进一步优选的方案:
进一步,所述步骤1中,步骤1中,当输入标示Flag=0时,程序自动统计计算成像数据的最大值MaxV和最小值MinV;
当输入标示Flag=1时,程序将输入成像数据的最大值、成像数据最小值参数,赋值给MaxV,MinV;同时对测井图像数据进行处理:
IMG(i,j)<MinV,IMG(i,j)=MinV
IMG(i,j)>MaxV,IMG(i,j)=MaxV
i=0,1,...N-1,j=0,1,....M-1
其中,IMG(i,j)表示超声图像数据IMG的第i行、第j列的数据值,N、M为整数,分别为超声图像总行数、总列数,i、j均为整数,i表示测井图像数据的第i行,第j列。
进一步,所述步骤2中,初次分割:将成像数据的最小值MinV到最大值MaxV的范围平均分为256份;
(1)初次分割:定义结构体变量structSecPar,其成员分别为当前分割范围的最小值fMin、最大值fMax、数据值个数iNum及当前数据值个数占与超声图像数据总个数的比值fPercent,定义vector容器类型的变量AutoCutP,里面存放的变量类型structSecPar,将MinV-MaxV分成256份,初始化AutoCutP,里面存放256个structSecPar结构体变量:
AutoCutP(i).fMin=(MaxV-MinV)*i/256+MinV;
AutoCutP(i).fMax=(MaxV-MinV)*(i+1)/256+MinV;
AutoCutP(i).iNum=0;
AutoCutP(i).fPercent=0;
其中,i=0,1,...255其中,i为整数,表示第i个数据范围;
(2)统计每个分割段中的超声图像个数:
如果IMG(i,j)≥AutoCutP[k].fMin并且IMG[i][j]<AutoCutP[k].fMax
则AutoCutP[k].iNum+=1;
i=0,1,...N-1,j=0,1,....M-1,k=0,1,...255
其中,i、j均为整数,IMG(i,j)表示超声图像数据IMG的第i行、第j列的数据值,k为整数,表示容器AutoCutP中的第k个structSecPar类型的变量;
(3)计算每个分割段中的超声图像个数占超声图像数据总个数的比值:
AutoCutP[k].fPercent=AutoCutP[k].iNum/N/M;
k=0,1,...255
其中,N、M为整数,分别为超声图像总行数、总列数,k为整数,AutoCutP[k]表示容器AutoCutP中的第k个structSecPar类型的变量。
进一步,所述步骤3中,测井图像值分布不均匀是指判断每个分隔段范围的所占总数百分比是否超出平均值1/256。
进一步,所述步骤3中,利用二分法对数据范围进行进一步分割方法具体步骤如下:
(1)如果AutoCutP[k].fPercent>1/256
k=0,1...Q,表示容器AutoCutP中的第k个structSecPar类型的变量,Q为整数,为容器AutoCutP中的structSecPar变量总个数,Q值是动态变化的;
则利用2分法,将AutoCutP[k].fMin-AutoCutP[k].fMax进一步分割成两份:
建立临时structSecPar结构体类型的变量ParTmp;
ParTmp.fMin=(AutoCutP(k).fMin+AutoCutP(k).fMax)/2
ParTmp.fMax=AutoCutP(k).fMax
ParTmp.iNum=0
(2)统计ParTmp.fMin-fMax范围的数据个数:
如果IMG(i,j)≥ParTmp.fMin并且IMG(i,j)<ParTmp.fMax
则ParTmp.iNum=ParTmp.iNum+1;
i=0,1,...N-1,j=0,1,....M-1,k=0,1,...255
则ParTmp.fPercent=ParTmp.iNum/N/M;
(3)AutoCutP(k).fMax=(AutoCutP(k).fMin+AutoCutP(k).fMax)/2
AutoCutP(k).iNum=AutoCutP(k).iNum-ParTmp.iNum
AutoCutP(k).fPercent=AutoCutP(k).iNum/N/M;
(4)将ParTmp变量插入AutoCutP中的第k个位置,单次自动分割完毕;
(5)重复步骤3,则可以进行无限次的自动分割。
进一步,所述步骤4中,不满足所有的分割份数的所占总数百分比全部小于平均1/256,进行无限次的自动分割,是由于出现单个数据值所占总数百分比大于1/256,则自动切分是无法将其所占总数百分比降低。
进一步,所述步骤4中,设定值为256份的4-5倍。
进一步,步骤5中,数据值映射是通过以下方法来实现的:
定义新的浮点数组,pMap,初始值设为0,即
pMap[i]=0i=0,1...Q,
其中,AutoCutP(j).表示容器AutoCutP中的第j个structSecPar类型的变量;
对频率进行累加:
Figure BDA0001770870540000041
其中,j为整数,i为整数,pMap[i]表示pMap第i个变量,Q为AutoCutP中的变量总个数,AutoCutP(j).表示容器AutoCutP中的第j个structSecPar类型的变量。
对数据值进行映射:如果IMG(i,j)≥AutoCutP(k).fMin
并且IMG(i,j)<AutoCutP(k).fMax
则映射结果IMGNew(i,j)=pMap(k)*255
i=0,1,...N-1,j=0,1,....N-1,k=0,1,....Q。
本发明的有益效果在于:
本发明能够对数据值的分布范围不均匀即分割范围的数据个数所占总数百分比大于1/256的分割范围,利用2分法,进行无限次的自动分割。以求达到显示效果均匀分布即所有的显示灰度尽可能的得到有效利用,与传统动态均衡相比,能够显著提高显示的分辨率和对比度。本发明处理流程操作简便,应用到测井图像的数据处理中,验证其具有较好的处理效果。
附图说明
图1为本发明流程示意图;
图2为采用本发明方法进行测井图像处理效果图。
具体实施方式
下面结合具体相关流量测井数据,对本发明的具体的实施方式进行说明。
如图1所示,一种基于二分法的测井图像均衡处理的方法,包括以下步骤:
步骤1,首先进行参数设定:
通过标示Flag的两种方式设定切割范围,当Flag=0时,程序自动统计计算最大值、最小值;当Flag=1时,通过人工手动输入方式设定切割范围,本实施例中,设定Flag=0,程序自动统计计算最大值。
步骤1中,当输入标示Flag=0时,程序自动统计计算成像数据的最大值MaxV和最小值MinV;
当输入标示Flag=1时,程序将输入的成像数据最大值、成像数据最小值参数,赋值给MaxVa,MinV;同时对测井图像数据进行处理:
IMG(i,j)<MinV,IMG(i,j)=MinV
IMG(i,j)>MaxV,IMG(i,j)=MaxV
i=0,1,...N-1,j=0,1,....N-1
其中,IMG(i,j)表示图像数据,i表示测井图像数据的第i行,j表示测井图像数据的第j列,N、M为整数,分别为超声图像总行数、总列数,i、j均为整数。
步骤2,对测井图像进行初始化分割:
根据灰度图像显示级别256级,将测井图像值的范围初始化分割为256份,统计每个分割范围内的数据点个数,及相应的所占数据总个数的百分比,程序中设定为256级。
步骤2中,初次分割:将成像数据的最小值MinV到最大值MaxV的范围平均分为256份,定义结构体数据structSecPar,其成员变量为当前切分范围的最小值fMin,最大值fMax,统计的属于当前范围的测井图像数据个数iNum及当前统计的属于当前范围的测井图像数据个数与测井图像总个数的比值fPercent。定义vector容器类型的变量AutoCutP,里面存放的变量类型structSecPar,将MinV-MaxV分成256份,初始化AutoCutP,里面存放256个structSecPar结构体变量。应用C语言容器类型vector,创建容器实例AutoCutP,该容器中存放着structSecPar类型的变量。切分方法如下:
初次分割:定义vector<structSecPar>vecAutoCutPar[256],利用将MinValue-MaxValue分成256份:
AutoCutP(i).fMin=(MaxV-MinV)*i/256+MinV;
AutoCutP(i).fMax=(MaxV-MinV)*(i+1)/256+MinV;
AutoCutP(i).iNum=0;
AutoCutP(i).fPercent=0;
其中,i=0,1,...255其中,i为整数,表示第i个数据范围;
(2)统计每个分割段中的超声图像个数:
如果IMG(i,j)≥AutoCutP[k].fMin并且IMG[i][j]<AutoCutP[k].fMax
则AutoCutP[k].iNum+=1;
i=0,1,...N-1,j=0,1,....M-1,k=0,1,...255
其中,i,j,k均为整数,IMG(i,j)表示超声图像数据IMG的第i行、第j列的数据值,k为整数,表示容器AutoCutP中的第k个structSecPar类型的变量;
(3)计算每个分割段中的超声图像个数占超声图像数据总个数的比值:
AutoCutP[k].fPercent=AutoCutP[k].iNum/N/M;
k=0,1,...255
其中,N、M为整数,分别为超声图像总行数、总列数,k为整数,AutoCutP[k]表示容器AutoCutP中的第k个structSecPar类型的变量。
步骤3,进行自动分割:
利用二分法对步骤2分割为256份测井图像值分布不均匀的数据范围进行进一步分割,对于分割后的不同范围,分别统计分割范围点的数据点个数,和所占数据总个数的百分比。测井图像值分布不均匀是指判断每个分隔段范围的所占总数百分比是否超出平均值1/256。
步骤3中,利用二分法对数据范围进行进一步分割方法,具体步骤如下:
(1)如果AutoCutP[k].fPercent>1/256
k=0,1...Q,其中,Q为AutoCutP中的变量总个数。
则利用2分法,将AutoCutP[k].fMin-AutoCutP[k].fMax进一步分割成两份:
建立临时structSecPar类型的变量ParTmp;
ParTmp.fMin=(AutoCutP(k).fMin+AutoCutP(k).fMax)/2
ParTmp.fMax=AutoCutP(k).fMax
ParTmp.iNum=0
(2)统计ParTmp.fMin-fMax范围的数据个数:
如果IMG(i,j)≥ParTmp.fMin并且IMG(i,j)<ParTmp.fMax
则ParTmp.iNum+=1;
i=0,1,..N-1,j=0,1,....N-1,k=0,1,...255
则ParTmp.fPercent=ParTmp.iNum/N/M;
(3)AutoCutP(k).fMax=(AutoCutP(k).fMin+AutoCutP(k).fMax)/2
AutoCutP(k).iNum=AutoCutP(k).iNum-ParTmp.iNum
AutoCutP(k).fPercent=AutoCutP(k).iNum/N/M;
(4)利用容器vector可动态插入数据特点;将ParTmp变量插入AutoCutP中的第k个位置,单次自动分割完毕;
(5)重复步骤3,则可以进行无限次的自动分割。
步骤4判断自动分割结束:
当所有的分割份数的所占总数百分比全部小于平均1/256,则程序自动分割结束;如果不满足,则重复步骤3,进行无限次的自动分割,则根据切分范围的总个数大于设定值(设定值为256份的4-5倍)时,自动分割强制退出,程序中设定切分范围的总个数大于256*10时,自动分割结束。进行无限次的自动分割,是由于出现单个数据值所占总数百分比大于1/256,则自动切分是无法将其所占总数百分比降低。
步骤5,根据自动分割好的数据范围进行数据值的映射,输出映射结果,即最终的处理结果如图2。
步骤5中,数据值映射是通过以下方法来实现的:
定义新的浮点数组,pMap,初始值设为0,即
pMap[i]=0i=0,1...Q,
其中,i为整数,pMap[i]表示pMap第i个变量,Q为AutoCutP中的变量总个数;
对频率进行累加:
Figure BDA0001770870540000081
其中,j为整数,i为整数,pMap[i]表示pMap第i个变量,Q为AutoCutP中的变量总个数,AutoCutP(j).表示容器AutoCutP中的第j个structSecPar类型的变量;
对数据值进行映射:如果IMG(i,j)≥AutoCutP(k).fMin
并且IMG(i,j)<AutoCutP(k).fMax
则映射结果IMGNew(i,j)=pMap(k)*255
i=0,1,...N-1,j=0,1,....N-1,k=0,1,....Q。
本发明的优点可通过以下处理结果进一步说明:
从处理效果图2来看,传统的图像动态均衡方法细节方面有些模糊,显示不够细腻,原因就是像素值分布不均匀导致的,没有对>1/256的数据值范围没有进行自动分割。而改进后的算法处理效果在细节显示方法有了很大的提升,对比度和分辨率效果有了很大的优化。
以上对本发明实施例所提供的实施方案进行了详细介绍,本文中应用了具体个例对本发明实施例的原理以及实施方式进行了阐述,以上实施例的说明只适用于帮助理解本发明实施例的原理;同时,对于本领域的一般技术人员,依据本发明实施例,在具体实施方式以及应用范围上均会有改变之处,综上所述,本说明书内容不应理解为对本发明的限制。

Claims (7)

1.一种基于2分法的测井图像均衡处理的方法,其特征在于,包括以下步骤:
步骤1,首先进行参数设定:
通过标示Flag的两种方式设定切割范围,当Flag=0时,程序自动统计计算最大值、最小值;当Flag=1时,通过人工手动输入方式设定切割范围;
步骤2,对测井图像进行初始化分割:
根据灰度图像显示级别256级,将测井图像值的范围初始化分割为256份,统计每个分割范围内的数据点个数,及相应的所占数据总个数的百分比;
步骤3,进行自动分割:
利用2分法对步骤2分割为256份测井图像值分布不均匀的数据范围进行进一步分割,对于分割后的不同范围,分别统计分割范围内的数据点个数,和所占数据总个数的百分比;
步骤4,判断自动分割结束:
当所有的分割份数的所占总数百分比全部小于平均1/256,则程序自动分割结束;如果不满足,则重复步骤3,进行无限次的自动分割,则根据切分范围的总个数大于设定值时,自动分割强制退出;
步骤5,根据自动分割好的数据范围进行数据值映射,输出映射结果。
2.根据权利要求1所述的方法,其特征在于,所述步骤1中,步骤1中,当输入标示Flag=0时,程序自动统计计算成像数据的最大值MaxV和最小值MinV;
当输入标示Flag=1时,程序将输入成像数据的最大值、成像数据最小值参数,赋值给MaxV,MinV;同时对测井图像数据进行处理:
IMG(i,j)<MinV,IMG(i,j)=MinV
IMG(i,j)>MaxV,IMG(i,j)=MaxV
i=0,1,...N-1,j=0,1,....M-1
其中,IMG(i,j)表示超声图像数据第i行,第j列的数据值,N、M为整数,分别为超声图像总行数、总列数,i、j均为整数,表示超声图像数据IMG的第i行,第j列。
3.根据权利要求1所述的方法,其特征在于,所述步骤2中,初始化分割:将成像数据的最小值MinV到最大值MaxV的范围平均分为256份;
(1)初次分割:定义结构体变量structSecPar,其成员分别为当前分割范围的最小值fMin、最大值fMax、数据值个数iNum及当前数据值个数与超声图像数据总个数的比值fPercent,定义vector容器类型的变量AutoCutP,里面存放的变量类型structSecPar,将MinV-MaxV分成256份,初始化AutoCutP,里面存放256个structSecPar结构体变量:
AutoCutP(i).fMin=(MaxV-MinV)*i/256+MinV;
AutoCutP(i).fMax=(MaxV-MinV)*(i+1)/256+MinV;
AutoCutP(i).iNum=0;
AutoCutP(i).fPercent=0;
其中,i=0,1,...255其中,i为整数,表示第i个数据范围;
(2)统计每个分割段中的超声图像数据点个数:
如果IMG(i,j)≥AutoCutP[k].fMin并且IMG[i][j]<AutoCutP[k].fMax
则AutoCutP[k].iNum+=1;
i=0,1,...N-1,j=0,1,....M-1,k=0,1,...255
其中,i、j均为整数,IMG(i,j)表示超声图像数据IMG的第i行、第j列的数据值,k为整数,表示容器AutoCutP中的第k个structSecPar类型的变量;
(3)计算每个分割段中的超声图像数据点个数占超声图像数据总个数的比值:
AutoCutP[k].fPercent=AutoCutP[k].iNum/N/M;
k=0,1,...255
其中,N、M为整数,分别为超声图像总行数、总列数,k为整数,AutoCutP[k]表示容器AutoCutP中的第k个structSecPar类型的变量。
4.根据权利要求1所述的方法,其特征在于,所述步骤3中,测井图像值分布不均匀是指存在分隔段范围的所占总数百分比是否超出平均值1/256。
5.根据权利要求3所述的方法,其特征在于,所述步骤3中,利用2分法对数据范围进行进一步分割方法,具体步骤如下:
(1)如果AutoCutP[k].fPercent>1/256
k=0,1...Q,表示容器AutoCutP中的第k个structSecPar类型的变量,Q为整数,为容器AutoCutP中的structSecPar变量总个数,Q值是动态变化的;
则利用2分法,将AutoCutP[k].fMin-AutoCutP[k].fMax进一步分割成两份:
建立临时structSecPar结构体类型的变量ParTmp;
ParTmp.fMin=(AutoCutP(k).fMin+AutoCutP(k).fMax)/2
ParTmp.fMax=AutoCutP(k).fMax
ParTmp.iNum=0
(2)统计ParTmp.fMin-fMax范围的数据个数:
如果IMG(i,j)≥ParTmp.fMin并且IMG(i,j)<ParTmp.fMax
则ParTmp.iNum=ParTmp.iNum+1;
i=0,1,...N-1,j=0,1,....M-1,k=0,1,...255
则ParTmp.fPercent=ParTmp.iNum/N/M;
(3)AutoCutP(k).fMax=(AutoCutP(k).fMin+AutoCutP(k).fMax)/2
AutoCutP(k).iNum=AutoCutP(k).iNum-ParTmp.iNum
AutoCutP(k).fPercent=AutoCutP(k).iNum/N/M;
(4)将ParTmp变量插入AutoCutP中的第k个位置,单次自动分割完毕;
(5)重复步骤3,则可以进行无限次的自动分割。
6.根据权利要求1所述的方法,其特征在于,所述步骤4中,设定值为256份的4-5倍。
7.根据权利要求3所述的方法,其特征在于,步骤5中,数据值映射是通过以下方法来实现的:
定义新的浮点数组,pMap,初始值设为0,即
pMap[i]=0 i=0,1...Q,
其中,i为整数,pMap[i]表示pMap第i个变量,Q为AutoCutP中的变量总个数;
对频率进行累加:
Figure FDA0002621430330000031
其中,AutoCutP(j).表示容器AutoCutP中的第j个structSecPar类型的变量;
对数据值进行映射:如果IMG(i,j)≥AutoCutP(k).fMin
并且IMG(i,j)<AutoCutP(k).fMax
则映射结果IMGNew(i,j)=pMap(k)*255
i=0,1,...N-1,j=0,1,....N-1,k=0,1,....Q。
CN201810948463.9A 2018-08-20 2018-08-20 一种基于二分法的测井图像均衡处理的方法 Active CN109300089B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810948463.9A CN109300089B (zh) 2018-08-20 2018-08-20 一种基于二分法的测井图像均衡处理的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810948463.9A CN109300089B (zh) 2018-08-20 2018-08-20 一种基于二分法的测井图像均衡处理的方法

Publications (2)

Publication Number Publication Date
CN109300089A CN109300089A (zh) 2019-02-01
CN109300089B true CN109300089B (zh) 2020-11-06

Family

ID=65165318

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810948463.9A Active CN109300089B (zh) 2018-08-20 2018-08-20 一种基于二分法的测井图像均衡处理的方法

Country Status (1)

Country Link
CN (1) CN109300089B (zh)

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108363110A (zh) * 2018-01-16 2018-08-03 中石化石油工程技术服务有限公司 成像测井计算页岩储层矿物含量及脆性指数的谱分析方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8948596B2 (en) * 2011-07-01 2015-02-03 CetusView Technologies, LLC Neighborhood node mapping methods and apparatus for ingress mitigation in cable communication systems

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108363110A (zh) * 2018-01-16 2018-08-03 中石化石油工程技术服务有限公司 成像测井计算页岩储层矿物含量及脆性指数的谱分析方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
《Exponential Contrast Maximization of Intensity Images》;Gyorgy Kovacs 等;《7th International Symposium on Image and Signal Processing and Analysis》;20110906;全文 *
《线阵CIS图像采集系统自动增益调节设计》;黄永林 等;《自动化仪表》;20110228;全文 *

Also Published As

Publication number Publication date
CN109300089A (zh) 2019-02-01

Similar Documents

Publication Publication Date Title
CN108564549A (zh) 一种基于多尺度稠密连接网络的图像去雾方法
CN106529543B (zh) 一种动态计算多色级二值化自适应阈值的方法及其系统
CN110084248A (zh) 一种orb特征均匀化提取方法
WO2022100510A1 (zh) 一种图像失真评测的方法、装置及计算机设备
WO2022077863A1 (zh) 视觉定位方法及相关模型的训练方法及相关装置、设备
CN108921868A (zh) 一种改进的Otsu阈值分割方法
CN110049332A (zh) 一种图像压缩方法、图像压缩装置及电子设备
WO2014186363A1 (en) Noise estimation for images
CN115272303A (zh) 基于高斯模糊的纺织物缺陷程度评估方法、装置及系统
CN105931236B (zh) 面向图像分割的模糊c-均值聚类初始聚类中心自动选取方法
CN109191482B (zh) 一种基于区域自适应光谱角阈值的图像合并分割方法
CN109300089B (zh) 一种基于二分法的测井图像均衡处理的方法
CN105828061B (zh) 一种基于视觉掩蔽效应的虚拟视点质量评价方法
Das et al. Histogram equalization techniques for contrast enhancement: a review
US20110052059A1 (en) Generating image histogram by parallel processing
CN110176021B (zh) 结合亮度校正的显著性信息的水平集图像分割方法及系统
CN115660994B (zh) 一种基于区域最小二乘估计的图像增强方法
CN110288604B (zh) 基于K-means的图像分割方法和装置
CN114998186A (zh) 基于图像处理的铜始极片表面结疤缺陷检测方法及系统
CN112419407A (zh) 基于云团边缘识别的云团位移矢量计算方法及装置
CN105260999A (zh) 一种直方图均衡插值的图像增强方法及系统
CN106651898B (zh) 一种基于蒙版插值的快速Grabcut算法
CN101105861A (zh) 可适性影像锐化方法
CN110619636A (zh) 基于rgb-d的变尺度图像分割方法
CN117828737B (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