CN109300083A - 一种分块处理Wallis匀色方法及装置 - Google Patents
一种分块处理Wallis匀色方法及装置 Download PDFInfo
- Publication number
- CN109300083A CN109300083A CN201811177734.1A CN201811177734A CN109300083A CN 109300083 A CN109300083 A CN 109300083A CN 201811177734 A CN201811177734 A CN 201811177734A CN 109300083 A CN109300083 A CN 109300083A
- Authority
- CN
- China
- Prior art keywords
- image
- standard deviation
- gray average
- pixel
- wallis
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 38
- 238000012545 processing Methods 0.000 title abstract description 14
- 230000009466 transformation Effects 0.000 claims abstract description 17
- 230000009467 reduction Effects 0.000 claims description 17
- 238000004364 calculation method Methods 0.000 claims description 9
- 230000008569 process Effects 0.000 claims description 6
- 238000006243 chemical reaction Methods 0.000 claims description 5
- 238000003860 storage Methods 0.000 claims description 4
- 238000004590 computer program Methods 0.000 claims description 3
- 238000010025 steaming Methods 0.000 description 5
- 238000010586 diagram Methods 0.000 description 4
- 238000004458 analytical method Methods 0.000 description 3
- 230000006870 function Effects 0.000 description 3
- VMXUWOKSQNHOCA-UKTHLTGXSA-N ranitidine Chemical compound [O-][N+](=O)\C=C(/NC)NCCSCC1=CC=C(CN(C)C)O1 VMXUWOKSQNHOCA-UKTHLTGXSA-N 0.000 description 3
- 238000011426 transformation method Methods 0.000 description 3
- 238000004422 calculation algorithm Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000013178 mathematical model Methods 0.000 description 2
- 238000005457 optimization Methods 0.000 description 2
- 230000009471 action Effects 0.000 description 1
- 230000004075 alteration Effects 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000008901 benefit Effects 0.000 description 1
- 238000004040 coloring Methods 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 239000006185 dispersion Substances 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000005286 illumination Methods 0.000 description 1
- 238000005192 partition Methods 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 230000007480 spreading Effects 0.000 description 1
- 238000003892 spreading Methods 0.000 description 1
- 230000001360 synchronised effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T3/00—Geometric image transformations in the plane of the image
- G06T3/40—Scaling of whole images or parts thereof, e.g. expanding or contracting
- G06T3/4007—Scaling of whole images or parts thereof, e.g. expanding or contracting based on interpolation, e.g. bilinear interpolation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T3/00—Geometric image transformations in the plane of the image
- G06T3/20—Linear translation of whole images or parts thereof, e.g. panning
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/10—Segmentation; Edge detection
- G06T7/11—Region-based segmentation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Image Processing (AREA)
Abstract
本发明涉及一种分块处理Wallis匀色方法及装置。通过对待处理的影像进行互不重叠分块并统计各影像块的灰度均值和标准偏差;根据角点的计算原则,计算各影像块的四个角点对应的灰度均值和标准偏差;根据每个像素到所在影像块边缘的距离和该像素所在影像块的四个角点的灰度均值和标准偏差,并利用双线性插值计算每个像素的灰度均值和标准偏差;根据每个像素的灰度均值和标准偏差,并结合参考影像的灰度均值和标准偏差,确定Wallis变换方程并根据该方程对待处理的影像进行匀色处理,实现了影像的匀色以便于大范围的影像镶嵌,解决了现有的Wallis变换在对大范围影像进行匀色时影像之间存在色彩和对比度都不一致的问题。
Description
技术领域
本发明涉及图像处理技术领域,特别是一种分块处理Wallis匀色方法及装置。
背景技术
在遥感领域中,大范围的镶嵌影像是影像分析和应用中的重要数据源。但是,受获取时刻光照、天气以及时相等条件的影响,待镶嵌影像之间存在着不同程度的色彩差异和对比度差异,使影像镶嵌更为复杂,也更为困难。因此,对区域范围内的多幅影像进行匀色,也称色彩一致性处理,消除影像之间色彩差异,使测区内影像在色调上保持一致,降低影像镶嵌的难度,具有重要的现实意义。
现有匀色方法可以分为两类,一类是非线性方法,如直方图匹配法和Gamma校正法;另一类是线性变换法,也是最受关注的方法,这类方法倾向于从影像的重叠区域选取不变像素作为样本参与统计,然后利用线性模型进行相对辐射校正。线性变换法中像素样本一般采用人工选取的方法,迭代加权的多元变化检测法,迭代慢特征分析法或者加权主成分分析法等进行筛选。上述的匀色方法虽然可以不同程度的解决色彩不一致问题,但是通常不能很好地消除对比度不一致现象。因此,一种基于均值和方差的特殊线性变换即Wallis变换,引入到色彩一致性处理中。Wallis变换是基于参考影像的处理方法,然而在进行大范围的影像镶嵌时,影像之间常有色彩不一致和对比度不一致同时存在的情况,因此该Wallis变换如何适用于大范围内影像之间的色彩一致性,需要进一步研究。
发明内容
本发明的目的是提供一种分块处理Wallis匀色方法及装置,用以解决现有的Wallis变换在对大范围影像进行匀色时影像之间存在色彩和对比度都不一致的问题。
为了实现影像的匀色以便于大范围的影像镶嵌,解决现有的Wallis变换在对大范围影像进行匀色时影像之间存在色彩和对比度都不一致的问题。本发明提供一种分块处理Wallis匀色方法,包括以下步骤:
1)对待处理的影像进行互不重叠分块并统计各影像块的灰度均值和标准偏差;
2)根据各影像块的灰度均值和标准偏差计算各影像块的四个角点对应的灰度均值和标准偏差;
3)根据每个像素到所在影像块边缘的距离和该像素所在影像块的四个角点的灰度均值和标准偏差,计算每个像素的灰度均值和标准偏差;
4)根据每个像素的灰度均值和标准偏差,并结合参考影像的灰度均值和标准偏差,对每个像素进行Wallis变换处理。
进一步地,为了角点赋值更加精确,步骤2)中角点若只属于一个影像块,则将该影像块的灰度均值和标准偏差赋给该角点;角点若为多个相邻影像块之间的公共角点,则将所属多个影像块的灰度均值和标准偏差的平均值赋给该角点。
进一步地,为了获得较好地匀色效果,步骤1)中对待处理的影像进行互不重叠分块,分块的个数为W×H,其中
W=r×w,H=r×h
式中,CV为影像的变异系数,CVRef为参考影像的变异系数,w和h分别为预设的行、列方向的参考分块数。
进一步地,为了实现不同影像间的色彩均衡,所述Wallis变换的公式为:
f(x,y)=[g(x,y)-m(x,y)]·(sf/s(x,y))+mf
式中,g(x,y)为待处理影像的灰度值,f(x,y)为处理后的结果影像的灰度值,mf为参考影像的灰度均值,sf为参考影像的标准偏差,m(x,y)和s(x,y)分别为像素的灰度均值和标准偏差。
进一步地,为了提高像素灰度均值和标准偏差的计算速度,步骤3)中每一个像素的灰度均值和标准偏差的计算为并行计算。
进一步地,为了实现对影像块均值和标准偏差的并行计算,步骤1)中各影像块的灰度均值和标准偏差为通过并行归约求和法求得。
为了更好的实现上述方法,本发明还提供一种分块处理Wallis匀色装置,包括存储器、处理器以及存储在存储器中并可在处理器上运行的计算机程序,所述处理器执行所述程序时实现以下步骤:
1)对待处理的影像进行互不重叠分块并统计各影像块的灰度均值和标准偏差;
2)根据各影像块的灰度均值和标准偏差计算各影像块的四个角点对应的灰度均值和标准偏差;
3)根据每个像素到所在影像块边缘的距离和该像素所在影像块的四个角点的灰度均值和标准偏差,计算每个像素的灰度均值和标准偏差;
4)根据每个像素的灰度均值和标准偏差,并结合参考影像的灰度均值和标准偏差,对每个像素进行Wallis变换处理。
进一步地,所述处理器包括至少两个线程块,一个影像块对应一个线程块,影像块中的像素与对应该影像块的线程块中的线程一一对应;步骤3)中一个像素的灰度均值和标准偏差的计算对应一个线程,所有线程并行运算。
进一步地,所述处理器在对各影像块的灰度均值和标准偏差求解时通过并行归约求和法进行运算。
进一步地,为了提高访问存储器的速度,所述存储器包括常量存储器,所述常量存储器用于存储四个角点的均值和标准偏差以及参考影像的灰度均值和标准偏差。
附图说明
图1是本发明的一种分块处理Wallis匀色方法中任意像素的计算原理图;
图2是本发明的GPU线程组织方式的原理图;
图3是本发明的归约求和法的原理图;
图4是本发明的并行归约求和的原理图;
图5是本发明的一种分块处理Wallis匀色方法的流程图。
具体实施方式
下面结合附图对本发明做进一步详细的说明。
本发明提供一种分块处理Wallis匀色方法,包括以下步骤:
1)对待处理的影像进行互不重叠分块并统计各影像块的灰度均值和标准偏差。
Wallis匀色法是一种特殊的线性变换法,通过将待处理影像的均值和标准偏差映射到参考影像的均值和标准偏差,实现不同影像间的色彩均衡。线性数学模型如下:
f(x,y)=g(x,y)r1+r0
式中,g(x,y)为待处理影像的灰度值,f(x,y)为处理后的结果影像的灰度值。r0为加性系数,r1为乘性系数,可以表示为:
r0=bmf+(1-b-r1)mg
式中,mg和mf分别为待处理影像和参考影像的灰度均值,sg和sf分别为待处理影像和参考影像的标准偏差。b为影像亮度系数,b∈[0,1]。c为影像方差扩展系数,c∈[0,1]。
Wallis变换的目的是将处理影像的灰度均值和标准偏差分别被强制到mf和sf,通常取b=1、c=1,此时线性数学模型变为:
f(x,y)=[g(x,y)-mg]·(sf/sg)+mf。
Wallis匀色法是根据影像整体的均值和标准偏差,利用同一个线性关系对影像的每一个像素进行处理。但是影像中的地物复杂多样,颜色信息也各不相同,影像整体的均值和标准偏差无法准确地反映局部地物色彩特征,采用同一个线性关系显然是不合理的。为了解决上述问题,本发明对待处理的影像进行互不重叠分块,分块的个数为W×H,其中
W=r×w,H=r×h
式中,CV为影像的变异系数,CVRef为参考影像的变异系数,w和h分别为预设的行、列方向的参考分块数。
采用分块策略时,匀色处理的质量受分块数目影响。若分块数目太多,即各影像块太小时,容易过度校正造成地物失真偏色并且计算量较大;若分块数目太少,即各影像块太大时,统计的均值方差不能准确反映地物分布,不能很好地消除影像之间的色彩差异。变异系数是标准偏差与均值之比,也称离散系数,可以描述影像内的物的丰富程度。当变异系数越大时,影像内的物种类越丰富,对应的分块数目应该越大即影像块越小,才能获得较好地匀色效果。
2)根据各影像块的灰度均值和标准偏差计算各影像块的四个角点对应的灰度均值和标准偏差。
其中,角点若只属于一个影像块,则将该影像块的灰度均值和标准偏差赋给该角点;角点若为多个相邻影像块之间的公共角点,则将所属多个影像块的灰度均值和标准偏差的平均值赋给该角点,最终计算得到各影像块的四个角点对应的灰度均值和标准偏差。
3)根据每个像素到所在影像块边缘的距离和该像素所在影像块的四个角点的灰度均值和标准偏差,计算每个像素的灰度均值和标准偏差。
如图1所示,对于影像块B(w,h)中的一点p(x,y)像素,其均值和标准偏差由角点Pw,h、Pw+1,h、Pw,h+1、Pw+1,h+1的均值和标准偏差以及该点到影像块边缘的距离Δx、Δy确定,利用双线性插值进行计算的公式如下:
式中,Δx为像素到m(w,h)对应的角点Pw,h的横向距离,Δy为像素到m(w,h)对应的角点Pw,h的纵向距离,m(x,y)和s(x,y)分别为像素的灰度均值和标准偏差,m(w,h)、m(w+1,h)、m(w,h+1)、m(w+1,h+1)分别为角点Pw,h、Pw+1,h、Pw,h+1、Pw+1,h+1的灰度均值,s(w,h)、s(w+1,h)、s(w,h+1)、s(w+1,h+1)分别为角点Pw,h、Pw+1,h、Pw,h+1、Pw+1,h+1的标准偏差,X和Y分别为该像素所在影像块B(w,h)的宽和高。
4)根据每个像素的灰度均值和标准偏差,并结合参考影像的灰度均值和标准偏差,对每个像素进行Wallis变换处理。
Wallis变换的方程为:
f(x,y)=[g(x,y)-m(x,y)]·(sf/s(x,y))+mf
式中,g(x,y)为待处理影像的灰度值,f(x,y)为处理后的结果影像的灰度值,mf为参考影像的灰度均值,sf为参考影像的标准偏差。
利用双线性插值计算每个像素的线性变换参数,可以保证相邻影像块之间的平滑性。此外,利用影像的各角点而非中心点参与计算,可以避免影像边缘的分块出现锯齿现象。
本发明还提供一种分块处理Wallis匀色装置,包括存储器、处理器以及存储在存储器中并可在处理器上运行的计算机程序,处理器执行程序时实现上述的匀色方法,本发明的分块处理Wallis匀色方法的主要计算任务包括三个步骤:计算各影像块的灰度均值和标准偏差、逐像素双线性插值和线性变换计算新的灰度值。分析可知:
(1)双线性插值和线性变换需要逐像素计算,是计算量最大的部分,各像素的计算是相互独立的,非常适合GPU并行处理,一个线程执行一个像素点的计算任务。
处理器包括至少两个线程块,一个影像块对应一个线程块,影像块中的像素与对应该影像块的线程块中的线程一一对应,双线性插值和线性变换为重复的密集计算任务,并行化比较简单,可以直接分配给各个线程同时计算,满负荷利用计算资源,可以大幅度缩短计算时间。线程是GPU的最小执行单元,具体执行时,按照“线程格网-线程块-线程”的层次结构进行组织,如图2所示。假设待处理影像大小为M×N,分配相同大小的线程格网,设置线程块大小为l×k,则线程块数目为一个影像块对应一个线程块,一个线程对应一个像素。各线程同时进行运算,并将结果按索引号赋给对应的像素。
(2)影像块灰度均值和标准偏差的本质均为累计求和计算,耦合度低,不能直接进行并行设计,因此处理器在对各影像块的灰度均值和标准偏差求解时通过并行归约求和法进行运算。
归约求和是一种缩减计算方法,基于对数步长交替两两求和,如图3所示,可以将求和的时间复杂度由O(N)(N为数据个数)降为O(log2N)。在每一个线程块内利用共享内存进行归约求和,可以达到并行加速的目的,并且交替策略可以避免存储片冲突并保持线程块的相邻线程处于活跃状态。
针对GPU的层次结构,本发明采用两遍归约求和策略。第一阶段内核执行n个并行归约,其中n是指线程块数,得到一个中间结果数组;第二个阶段通过调用一个线程块对这个中间数组进行归约,从而得到最终结果,如图4所示。具体步骤如下:
1、对输入数组落入每个线程中的数据进行求和。每个线程把它得到的累计值写入共享内存,交替因子为n*m,并在执行对数步长的归约前进行同步操作。
2、对共享内存中的值进行对数步长的归约操作。共享内存中后半部分的值被加到前半部分,即a[i]=a[i]+a[m/2],(0≤i<m/2),参与的线程依次减半。在此操作执行log2m次后,共享内存中第一个线程对应的值a[0]即为该线程块的和。共享内存的大小等于线程块的线程数量m,并且m必须是2的幂次。
3、线程块的和写入全局内存。
利用并行归约求和方法计算影像像素值之和以及像素个数,即可求出影像均值,同理可求标准偏差,本发明方法的流程如图5所示。
另外,为了达到尽可能高的并行加速比,本发明结合算法的自身特点,从配置划分、存储器带宽和指令吞吐量等方面进行优化。具体如下:
(1)合理组织线程
GPU中线程块独立被调度到流式多处理器中,来自同一个线程块的线程在同一个流式多处理器中执行。为了使流式多处理器的性能达到最优,流式多处理器中线程块的数量要小于8、线程数等于1536,线程块中的线程数要小于1024。分析可知,当线程块大小设置为256或512时,流式多处理器的性能最优,可以提高并行算法的运算速度。
(2)存储器优化
由于常量存储器的访问速度明显优于全局存储器,合理利用常量存储器代替全局存储器可以有效地减少内存流量和内存带宽。双线性插值中的四个角点的均值和标准偏差和线性变换中的参考均值和参考标准偏差数据量不大且每个线程都会读取,将这些参数分配为常量内存,可以大大加速访问数据速度和减少内存带宽,有效地提升程序的运行效率。
针对共享内存中的归约求和,为了减少不必要的线程同步,使用线程束同步优化。由于每个线程块中的线程束是按照锁步方式执行每条指令的,当线程块中的活动线程数低于硬件线程束的大小32时,无须再调用_syncthreads()内置函数。线程束同步可以避免每个Warp中出现分支导致效率低下,减小线程的闲置,提高并行度。
(3)指令优化
为了使指令吞吐量最大化,在满足精度的前提下,尽可能使用单精度float类型代替双精度double类型,使用硬件函数代替常规函数,用最少的指令完成相同的运算。
以上给出了本发明涉及的具体实施方式,但本发明不局限于所描述的实施方式。在本发明给出的思路下,采用对本领域技术人员而言容易想到的方式对上述实施例中的技术手段进行变换、替换、修改,并且起到的作用与本发明中的相应技术手段基本相同、实现的发明目的也基本相同,这样形成的技术方案是对上述实施例进行微调形成的,这种技术方案仍落入本发明的保护范围内。
Claims (10)
1.一种分块处理Wallis匀色方法,其特征在于,包括以下步骤:
1)对待处理的影像进行互不重叠分块并统计各影像块的灰度均值和标准偏差;
2)根据各影像块的灰度均值和标准偏差计算各影像块的四个角点对应的灰度均值和标准偏差;
3)根据每个像素到所在影像块边缘的距离和该像素所在影像块的四个角点的灰度均值和标准偏差,计算每个像素的灰度均值和标准偏差;
4)根据每个像素的灰度均值和标准偏差,并结合参考影像的灰度均值和标准偏差,对每个像素进行Wallis变换处理。
2.根据权利要求1所述的分块处理Wallis匀色方法,其特征在于,步骤2)中角点若只属于一个影像块,则将该影像块的灰度均值和标准偏差赋给该角点;角点若为多个相邻影像块之间的公共角点,则将所属多个影像块的灰度均值和标准偏差的平均值赋给该角点。
3.根据权利要求1或2所述的分块处理Wallis匀色方法,其特征在于,步骤1)中对待处理的影像进行互不重叠分块,分块的个数为W×H,其中
W=r×w,H=r×h
式中,CV为影像的变异系数,CVRef为参考影像的变异系数,w和h分别为预设的行、列方向的参考分块数。
4.根据权利要求3所述的分块处理Wallis匀色方法,其特征在于,所述Wallis变换的公式为:
f(x,y)=[g(x,y)-m(x,y)]·(sf/s(x,y))+mf式中,g(x,y)为待处理影像的灰度值,f(x,y)为处理后的结果影像的灰度值,mf为参考影像的灰度均值,sf为参考影像的标准偏差,m(x,y)和s(x,y)分别为像素的灰度均值和标准偏差。
5.根据权利要求1所述的分块处理Wallis匀色方法,其特征在于,步骤3)中每一个像素的灰度均值和标准偏差的计算为并行计算。
6.根据权利要求5所述的分块处理Wallis匀色方法,其特征在于,步骤1)中各影像块的灰度均值和标准偏差为通过并行归约求和法求得。
7.一种分块处理Wallis匀色装置,包括存储器、处理器以及存储在存储器中并可在处理器上运行的计算机程序,其特征在于,所述处理器执行所述程序时实现以下步骤:
1)对待处理的影像进行互不重叠分块并统计各影像块的灰度均值和标准偏差;
2)根据各影像块的灰度均值和标准偏差计算各影像块的四个角点对应的灰度均值和标准偏差;
3)根据每个像素到所在影像块边缘的距离和该像素所在影像块的四个角点的灰度均值和标准偏差,计算每个像素的灰度均值和标准偏差;
4)根据每个像素的灰度均值和标准偏差,并结合参考影像的灰度均值和标准偏差,对每个像素进行Wallis变换处理。
8.根据权利要求7所述的分块处理Wallis匀色装置,其特征在于,所述处理器包括至少两个线程块,一个影像块对应一个线程块,影像块中的像素与对应该影像块的线程块中的线程一一对应;步骤3)中一个像素的灰度均值和标准偏差的计算对应一个线程,所有线程并行运算。
9.根据权利要求8所述的分块处理Wallis匀色装置,其特征在于,所述处理器在对各影像块的灰度均值和标准偏差求解时通过并行归约求和法进行运算。
10.根据权利要求8或9所述的分块处理Wallis匀色装置,其特征在于,所述存储器包括常量存储器,所述常量存储器用于存储四个角点的均值和标准偏差以及参考影像的灰度均值和标准偏差。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811177734.1A CN109300083B (zh) | 2018-10-10 | 2018-10-10 | 一种分块处理Wallis匀色方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811177734.1A CN109300083B (zh) | 2018-10-10 | 2018-10-10 | 一种分块处理Wallis匀色方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109300083A true CN109300083A (zh) | 2019-02-01 |
CN109300083B CN109300083B (zh) | 2022-10-11 |
Family
ID=65162125
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811177734.1A Active CN109300083B (zh) | 2018-10-10 | 2018-10-10 | 一种分块处理Wallis匀色方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109300083B (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110031461A (zh) * | 2019-02-14 | 2019-07-19 | 江苏恒力化纤股份有限公司 | 一种涤纶长丝染色均匀度试验方法 |
CN111754590A (zh) * | 2020-05-14 | 2020-10-09 | 北京吉威空间信息股份有限公司 | 基于全球色彩特征库遥感影像自动匀色的方法 |
CN112669309A (zh) * | 2021-01-07 | 2021-04-16 | 航天宏图信息技术股份有限公司 | 遥感影像的局部匀色方法、装置、电子设备及存储介质 |
CN113421193A (zh) * | 2021-05-05 | 2021-09-21 | 桂林理工大学 | 均值-方差成本函数最小模型的多幅影像镶嵌辐射均衡方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105427372A (zh) * | 2015-06-11 | 2016-03-23 | 北京吉威时代软件股份有限公司 | 基于tin网的正射影像拼接色彩一致性处理技术 |
WO2017026051A1 (ja) * | 2015-08-11 | 2017-02-16 | オリンパス株式会社 | 画像処理装置、画像処理方法および画像処理プログラム |
CN106485664A (zh) * | 2015-08-27 | 2017-03-08 | 上海沃韦信息科技有限公司 | 一种基于小波变换和Wallis变换的卫星图像色彩平衡方法 |
-
2018
- 2018-10-10 CN CN201811177734.1A patent/CN109300083B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105427372A (zh) * | 2015-06-11 | 2016-03-23 | 北京吉威时代软件股份有限公司 | 基于tin网的正射影像拼接色彩一致性处理技术 |
WO2017026051A1 (ja) * | 2015-08-11 | 2017-02-16 | オリンパス株式会社 | 画像処理装置、画像処理方法および画像処理プログラム |
CN106485664A (zh) * | 2015-08-27 | 2017-03-08 | 上海沃韦信息科技有限公司 | 一种基于小波变换和Wallis变换的卫星图像色彩平衡方法 |
Non-Patent Citations (2)
Title |
---|
周丽雅等: "反差一致性保持的影像匀光算法", 《测绘科学技术学报》 * |
李凤英等: "一种镶嵌影像的分块直方图匹配匀光算法", 《遥感信息》 * |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110031461A (zh) * | 2019-02-14 | 2019-07-19 | 江苏恒力化纤股份有限公司 | 一种涤纶长丝染色均匀度试验方法 |
CN110031461B (zh) * | 2019-02-14 | 2022-03-18 | 江苏恒力化纤股份有限公司 | 一种涤纶长丝染色均匀度试验方法 |
CN111754590A (zh) * | 2020-05-14 | 2020-10-09 | 北京吉威空间信息股份有限公司 | 基于全球色彩特征库遥感影像自动匀色的方法 |
CN111754590B (zh) * | 2020-05-14 | 2024-04-02 | 北京吉威空间信息股份有限公司 | 基于全球色彩特征库遥感影像自动匀色的方法 |
CN112669309A (zh) * | 2021-01-07 | 2021-04-16 | 航天宏图信息技术股份有限公司 | 遥感影像的局部匀色方法、装置、电子设备及存储介质 |
CN112669309B (zh) * | 2021-01-07 | 2024-06-11 | 航天宏图信息技术股份有限公司 | 遥感影像的局部匀色方法、装置、电子设备及存储介质 |
CN113421193A (zh) * | 2021-05-05 | 2021-09-21 | 桂林理工大学 | 均值-方差成本函数最小模型的多幅影像镶嵌辐射均衡方法 |
Also Published As
Publication number | Publication date |
---|---|
CN109300083B (zh) | 2022-10-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109300083A (zh) | 一种分块处理Wallis匀色方法及装置 | |
Xu et al. | High-performance iterative electron tomography reconstruction with long-object compensation using graphics processing units (GPUs) | |
Zhang et al. | Image parallel processing based on GPU | |
CN107451955A (zh) | 一种k‑t算法重建天文图像中斑点图的并行化实现方法 | |
DE102021120605A1 (de) | Effiziente softmax-berechnung | |
CN107657599B (zh) | 基于混合粒度划分和动态负载分配的遥感图像融合系统并行实现方法 | |
Liu et al. | GPU-ClustalW: Using graphics hardware to accelerate multiple sequence alignment | |
CN109410136B (zh) | 基于最短传递路径的匀色方法及处理装置 | |
Spiliotis et al. | Parallel implementation of the image block representation using OpenMP | |
CN104123119B (zh) | 基于gpu的动态视觉测量特征点中心快速定位方法 | |
Fuentes-Alventosa et al. | GUD-Canny: a real-time GPU-based unsupervised and distributed Canny edge detector | |
US20090322781A1 (en) | Anti-aliasing techniques for image processing | |
CN110083446A (zh) | 一种gpu并行与零io模式下遥感影像实时处理方法及系统 | |
CN105869105B (zh) | 一种针对a+超分辨率技术的gpu加速方法 | |
Ke et al. | Efficient GPU computing framework of cloud filtering in remotely sensed image processing | |
Wang et al. | Egpuip: An embedded gpu accelerated library for image processing | |
Zhang et al. | Weighted mean shift object tracking implemented on gpu for embedded sustems | |
KR20190118023A (ko) | 이미지 자동 회귀 보간 장치 및 방법 | |
Shuhua et al. | GPU-based barrel distortion correction for acceleration | |
Gutierrez et al. | A fast level-set segmentation algorithm for image processing designed for parallel architectures | |
Colom et al. | Interactive VPL-based global illumination on the GPU using fuzzy clustering | |
Agulleiro et al. | Dynamic load scheduling on CPU-GPU for iterative tomographic reconstruction | |
Bethel et al. | Improving performance of structured-memory, data-intensive applications on multi-core platforms via a space-filling curve memory layout | |
Mahmud et al. | Image histogram equalization using parallel processing (GPU computation) in open CL to correct the contrast of images | |
Ke et al. | Optimized GPU implementation for dynamic programming in image data processing |
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 |