谁知道ENVI里面图像增强中(Enhance)高斯增强是什么原理,文档或者链接都行啊

谁知道ENVI里面图像增强中(Enhance)高斯增强是什么原理,文档或者链接都行啊,第1张

高斯滤波(高斯平滑)是图像处理,计算机视觉里面最常见的操作。平时,我们都是用matlab或者opencv的函数调用:imfilter或者cvSmooth,并不关心底层的实现。然而当开发者要实做高斯滤波的时候,往往就会很迷惘,往往会被以下几个问题困扰:

给定sigma,即标准偏差,怎么确定离散化后滤波器的窗口大小?

给定窗口大小,怎么计算高斯核的sigma,即标准方差?

怎么实现可分离滤波器?

我在google上搜了一下,还真没几个人把实现的细节讲清楚了。这里,我尝试结合三份源码,做个小小的总结。

三份源码分别是:

opencv中的cvfiltercpp

autopano-sift-c中的GaussianConvolutionc

GIMP中的blur-gaussc和unsharp-maskc

在图像处理中,高斯滤波一般有两种实现方式,一是用离散化窗口滑窗卷积,另一种通过傅里叶变换。最常见的就是第一种滑窗实现,只有当离散化的窗口非常大,用滑窗计算量非常大(即使用可分离滤波器的实现)的情况下,可能会考虑基于傅里叶变化的实现方法。这里我们只讨论第一种方法。

二维高斯函数的形式是这样的:

f(x,y) = A e^{- \left(\frac{(x-x_o)^2}{2\sigma_x^2} + \frac{(y-y_o)^2}{2\sigma_y^2} \right)}

有着如下的形状,形状很激凸,怎么会叫高斯平滑呢,分明是高斯激凸嘛:

基本上,离散化的主旨就是保留高斯函数中心能量最集中的中间部分,忽略四周能量很小的平坦区域。这只是个很感性的描述,具体实现起来,就会出现千奇百怪的版本。下面结合三份源码,看看现实世界里的高斯平滑到底长的什么样子。

首先是第一个问题:给定sigma,怎么计算窗口大小?

直接上opencv的源码,在cvFilter函数中:

param1 = cvRound(sigma1(depth == CV_8U 3 : 4)2 + 1)|1;

opencv认为半径为3sigma的窗口就是高斯函数能量最集中的区域。(为什么在图像深度不是8U的时候,使用4sigma半径的窗口就不得而知了,有高人指点一下?)

autopan0-sift-c是图像拼接软件hugin里面的sift实现,在实现DoG的时候需要做不同尺度的高斯平滑,实现如下,在GaussianConvolution_new1函数中:

dim = 1 + 2 ((int) (30 sigma));

可见autopano也是实现的3sigma半径的窗口。

在GIMP里,实现比较奇特,在blur_gaussc的make_rle_curve函数里面,

const gdouble sigma2 = 2 sigma sigma;

const gdouble l = sqrt (-sigma2 log (10 / 2550));

int n = ceil (l) 2;

if ((n % 2) == 0)

n += 1;

我是没看懂那个 log (10 / 2550)是干嘛的。。。惭愧。效果来看,这个实现的窗口半径是约等于22sigma。

然后是第二个问题:给定窗口大小,怎么计算sigma?

opencv的实现,在cvFiltercpp的init_gaussian_kernel函数中:

sigmaX = sigma > 0 sigma : (n/2 – 1)03 + 08;

再次不可解。。。乘以03还可以接受,加上08是为嘛啊?

autopano没有实现这个特性。

GIMP的实现:

/ we want to generate a matrix that goes out a certain radius

from the center, so we have to go out ceil(rad-05) pixels,

inlcuding the center pixel Of course, that’s only in one direction,

so we have to go the same amount in the other direction, but not count

the center pixel again So we double the previous result and subtract

one

The radius parameter that is passed to this function is used as

the standard deviation, and the radius of effect is the

standard deviation 2 It’s a little confusing

/

radius = fabs (radius) + 10;

std_dev = radius;

radius = std_dev 2;

/ go out ‘radius’ in each direction /

matrix_length = 2 ceil (radius – 05) + 1;

注释讲的很清楚了,基本上就是认为sigma应该等于窗口半径的一半。

看完这三分源码,结论就是,关于sigma和半径,你爱怎么算就怎么算吧,差不多就行。。。(额。。费了半天劲,就得到这么一句废话啊)。

第三个问题是可分离滤波器:

由于高斯函数可以写成可分离的形式,因此可以采用可分离滤波器实现来加速。所谓的可分离滤波器,就是可以把多维的卷积化成多个一维卷积。具体到二维的高斯滤波,就是指先对行做一维卷积,再对列做一维卷积。这样就可以将计算复杂度从O(MMNN)降到O(2MMN),M,N分别是图像和滤波器的窗口大小。问题是实现时候怎么计算一维的卷积核呢?

其实很简单,按照前面计算出来的窗口大小,计算所有离散点上一维高斯函数的权值,最后别忘了将权值之和归一化到1

有码有真相,来自opencv:

for( i = 0; i <= n/2; i++ )

{

double t = fixed_kernel (double)fixed_kernel[i] : exp(scale2Xii);

if( type == CV_32FC1 )

{

cf[(n/2+i)step] = (float)t;

sum += cf[(n/2+i)step]2;

}

else

{

cd[(n/2+i)step] = t;

sum += cd[(n/2+i)step]2;

}

}

sum = 1/sum;

for( i = 0; i <= n/2; i++ )

{

if( type == CV_32FC1 )

cf[(n/2+i)step] = cf[(n/2-i)step] = (float)(cf[(n/2+i)step]sum);

else

cd[(n/2+i)step] = cd[(n/2-i)step] = cd[(n/2+i)step]sum;

}

一.添加北京54坐标

通过修改ENVI相应的系统文档,可以将北京54坐标系添加进ENVI软件系统,具体方法和步骤如下:

1、找到ENVI安装目录下相应的坐标系存储文档,如“C:\RSI\IDL63\products\envi43\map_proj\datumtxt”,并打开。

2、将以下描述语句添加在文档最后一行:“Beijing-54, Krassovsky, -12, -113, -41”,其中“Beijing-54”是新坐标的名称; “Krassovsky”是它所使用的椭球体名称;“-12, -113, -41”是新坐标系与WGS-84坐标系的坐标差。

3、将修改后的文档保存,关闭。

4、重新启动ENVI。新的坐标系添加成功。

二、自定义高斯-克吕格(Kauss Kruger)投影

具体方法和步骤如下:

1.选择 Map Customize Map Projection,或在任何地图投影选择对话框中,点击“New”按钮 。

2、当出现 Customized Map Projection Definition 对话框时,假设为定义高斯投影的13带,13带的中央经度为75度,具体设置参数如下(上图所示):

Projection name :GK Zone 13(Beijing-54)

Projection Type:Transverse Mercator

Projection Datum:Beijing-54

False easting:500000

False northing:0

Latitude of projection origin: 0

Longitude of central meridian: 75

Scale factor:1

3.选择 Projection Add New Projection将投影添加到 ENVI 所用的投影列表中。

4.存储新的或更改过的投影信息,选择 File Save Projections。选择 File Cancel 退出该功能。

三、地图投影转换

使用 Convert Map Projection 选项可以将经过地理坐标定位的文件转换成其他投影,该转换通过将文件纠正到其他的投影来完成。

1、 选择 Map Convert Map Projection。选择所需输入文件,点击“OK”。将出现Convert Map Projection Parameters 对话框。三、地图投影转换

使用 Convert Map Projection 选项可以将经过地理坐标定位的文件转换成其他投影,该转换通过将文件纠正到其他的投影来完成。

1、 选择 Map Convert Map Projection。选择所需输入文件,点击“OK”。将出现Convert Map Projection Parameters 对话框。

2、点击Chang Proj…切换按纽,在“ Projection Select”列表中,通过点击需要转换的投影类型GK Zone 13(Beijing-54)。OK

选择File->Open,选择_MTLtxt文件打开。

ENVI自动显示RGB显示真彩色图像,打开DataManager对话框,可以看到ENVI自动读取元数据信息,包括中心波长信息、波段名称等。

欢迎分享,转载请注明来源:品搜搜测评网

原文地址:https://pinsoso.cn/meirong/1786594.html

(0)
打赏 微信扫一扫微信扫一扫 支付宝扫一扫支付宝扫一扫
上一篇 2023-10-18
下一篇2023-10-18

随机推荐

  • 关于雪花秀精华露的使用方法,雪花秀润致焕活肌底精华露怎么用呢?

    本文主要讲述关于“关于雪花秀精华露的使用方法,雪花秀润致焕活肌底精华露怎么用呢?”的内容。雪花秀润燥精华是雪花秀明星产品中销量第一的产品,可谓是当家花旦级别的产品,对于这样一款效果不错的产品我们应该掌握正确的使用方法达到最为理想的护理效果。

    2024-04-15
    64100
  • 娇兰和雅诗兰黛哪个好?雅诗兰黛哪些产品好用?

    雅诗兰黛大家应该并不陌生,他家的气垫BB,粉底液,口红,眼霜和精华都挺出名的,价格也适中,很多人都能接受。但是,娇兰是贵妇级的,有着悠久的历史,以香水起家,护肤品和彩妆都很出名,其中一个黑兰面霜八千多块,并不是很多人都能接受的。1、娇兰和雅

    2024-04-15
    57800
  • 化妆品公司有哪几家

    国际知名化妆品品牌前十名有:兰蔻、雅诗兰黛、迪奥、资生堂、倩碧、海蓝之谜、SK2、香奈儿、赫莲娜、伊丽莎白雅顿。一、兰蔻兰蔻1935年诞生于法国,是由Armand Petitjean(阿曼达·珀蒂让)创办的品牌。作为全球知名的高端化妆品品牌

    2024-04-15
    55300
  • 妮维雅洗面奶男士怎么样

    妮维雅这个品牌,一直专业的致力于男性护肤产品的开发和研究,其中也根据男性的生活特点和肤质的特点,来研究出了很多适合的产品,而关于妮维雅男士洗面奶,也是可以帮助去除多余油脂,非常适合男性肌肤使用泡沫,也非常的丰富。在使用之后,会感觉我们的皮肤

    2024-04-15
    51300
  • 精华露和保湿露的区别

    精华液是护肤品中有效含量最高的,同时也是效果最好的,一般精华液都比较贵,欧美的动辄都是几千块,国产的也要好几百。但是这些都真的好用吗?其实很多精华液并没有那么好,相反有些不太知名的精华液倒是效果非常的好,后来才知道,好用是有道理的,比如研发

    2024-04-15
    37600
  • 精华露和精华液一样吗?

    精华水、精华液、精华露、是一样的吗,还是有区别的? 精华水就是比较粘稠比较营养的化妆水。 精华液就是用在化妆水之后那一步,精华露要比精华液稍微稀一点。精华素和精华液一样吗?什么区别 精华液的分子比较小,质地较稀,易吸收,用了之后

    2024-04-15
    49000
  • 艾灸盒怎么用?

    艾灸盒怎么用?艾灸盒有很多种,我们平时最常见的要属随身灸艾灸盒,这种艾灸盒有铜制的也有铁制的,一般大多数选择用铁质的,下面给大家详细介绍一下铁制的随身灸,它的具体用法:首先把艾灸盒打开,把艾条剪成一小段,插到里面的支架上面,然后点燃艾灸条,

    2024-04-15
    40000

发表评论

登录后才能评论
保存