转:opencv 坐标系变换 logPolar()和linearPolar()

转载:https://blog.csdn.net/zhubaohua_bupt/article/details/83932171

logPolar是把数据从笛卡尔转到对数极坐标系,linearPolar把数据从笛卡尔坐标系转到极坐标坐标系。

说到笛卡尔坐标系到极坐标系的转换,大家并不陌生,公式如下:


image.png
image.png

笛卡尔转到对数极坐标系的转换多了一层对数,可表示如下

image.png

这里log以自然数e为底。

接下来,首先分析一下logPolar和linearPolar输入输出参数的含义以及用法,然后给出一些实验结果。

1 logPolar()

image

这里把imgproc.hpp里的函数接口粘了过来,从接口说明里我们可以看到,

输入输出图像类型尺寸一致,输入参数还包括center():转换中心,M: scale 参数和flags(变换后图像插值方法)

logPolar()是这样做的转换:


I = (dx,dy) = (x - center.x,y - center.y) \\
\rho = M \cdot log_e(\texttt{magnitude} (I)) ,\\
\phi = Ky \cdot \texttt{angle} (I)_{0..360 deg} \\

M = src.cols / log_e(maxRadius) \\
Ky = src.rows / 360 \\

上述说明中,除了基本的变换公式外,M和Ky这两个参数需要注意一下,这两个参数决定着变换后半径[图片上传失败...(image-ca21af-1604301036485)]

和角度[图片上传失败...(image-c1d07e-1604301036485)]

的尺度,也就是说,在变换后的极坐标图像上,x、y方向上单位像素与半径[图片上传失败...(image-bb889a-1604301036485)]

和角度[图片上传失败...(image-33dad6-1604301036485)]

的关系。

  • M决定输出图像x轴([图片上传失败...(image-cbe3ef-1604301036484)]

    )的尺度,Ky决定y轴([图片上传失败...(image-dc4d21-1604301036484)]

    )的尺度。

  • 注意到,函数没有让我们手动设置Ky的值,Ky =src.rows / 360已经写死,这样的话,输出极坐标图像y轴从y=0到yy=row-1的范围表示角度从0到360度的变化。

参数M没有写死,我们可以通过接口改变,那M = src.cols / log_e(maxRadius), 代表什么意思呢, maxRadius是什么呢?

`回过头来看笛卡尔向极坐标的变换,现在把输入图像看作是一个二维笛卡尔坐标系,如下图我们可以自己设定center,假设center为图像中心,那么maxRadius=[图片上传失败...(image-85e64c-1604301036484)]

,那么我们就可以看出,此时乘M的作用就是在输出图像的x轴上,將log(`[图片上传失败...(image-9479bd-1604301036485)]

)从x=0拉伸到x=col-1。

![image](https://upload-images.jianshu.io/upload_images/4688503-dce4d2219e4dbbc4?imageMogr2/auto-orient/strip%7CimageView2/2/w/1240)

image

                  笛卡尔坐标系图输入图像                                         对数极坐标系下图像

2 linearPolar()

image

linearPolar和logPolar类似,比logPolar简单些,其变换如下图。

image
image
      笛卡尔坐标系图输入图像                                                        极坐标系下图像

3实验

下面来看几个比较特殊的实验,为了方便,用python实现:

def polar(img):
  h, w = img.shape[:2]
  maxRadius = math.hypot(w/2,h/2)
  m = w / math.log(maxRadius)
  log_polar = cv2.logPolar(img, (w/2, h/2), m, cv2.WARP_FILL_OUTLIERS + cv2.INTER_LINEAR) 
  linear_polar = cv2.linearPolar(img, (w/2, h/2), maxRadius, cv2.WARP_FILL_OUTLIERS + cv2.INTER_LINEAR) 
  cv2.imshow("log_polar",log_polar)
  cv2.imshow("linear_polar",linear_polar)
  cv2.waitKey(1000)

第一组: logPolar()

image
image
image
                     src                                   m = `src.cols / log_e(maxRadius)` /2        m = `src.cols /log_e(maxRadius)`

第二组:

image
          ![image](https://upload-images.jianshu.io/upload_images/4688503-ffa8c87293bfa9ce?imageMogr2/auto-orient/strip%7CimageView2/2/w/1240)
image
                                                                                        logPolar()                                                      linearPolar()

         src                                                  m = `src.cols /log_e(maxRadius)`                                   `maxRadius`    

第三组:

image
    ![image](https://upload-images.jianshu.io/upload_images/4688503-54b63a5058210a71?imageMogr2/auto-orient/strip%7CimageView2/2/w/1240)

         ![image](https://upload-images.jianshu.io/upload_images/4688503-a8d0689816fe805f?imageMogr2/auto-orient/strip%7CimageView2/2/w/1240)

                                                                            logPolar()                                                      linearPolar()

                src                                    m = `src.cols /log_e(maxRadius)`                              `maxRadius`
  • 可以想一下在这几个例子中,变换后最右侧的尖儿为什有的能连延伸到最右侧,有的不能。
  • 从笛卡尔坐标系来看,这三组的src貌似形状很完美,但经过转换后可以看出后两组的瑕疵了,如果后两组也是完美对称的,那么在极坐标系下,第二组纵向应该是笔直的,第三组横行应该是平行的。
©著作权归作者所有,转载或内容合作请联系作者
  • 序言:七十年代末,一起剥皮案震惊了整个滨河市,随后出现的几起案子,更是在滨河造成了极大的恐慌,老刑警刘岩,带你破解...
    沈念sama阅读 225,764评论 6 524
  • 序言:滨河连续发生了三起死亡事件,死亡现场离奇诡异,居然都是意外死亡,警方通过查阅死者的电脑和手机,发现死者居然都...
    沈念sama阅读 96,848评论 3 409
  • 文/潘晓璐 我一进店门,熙熙楼的掌柜王于贵愁眉苦脸地迎上来,“玉大人,你说我怎么就摊上这事。” “怎么了?”我有些...
    开封第一讲书人阅读 173,181评论 0 370
  • 文/不坏的土叔 我叫张陵,是天一观的道长。 经常有香客问我,道长,这世上最难降的妖魔是什么? 我笑而不...
    开封第一讲书人阅读 61,430评论 1 303
  • 正文 为了忘掉前任,我火速办了婚礼,结果婚礼上,老公的妹妹穿的比我还像新娘。我一直安慰自己,他们只是感情好,可当我...
    茶点故事阅读 70,451评论 6 403
  • 文/花漫 我一把揭开白布。 她就那样静静地躺着,像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上,一...
    开封第一讲书人阅读 53,879评论 1 316
  • 那天,我揣着相机与录音,去河边找鬼。 笑死,一个胖子当着我的面吹牛,可吹牛的内容都是我干的。 我是一名探鬼主播,决...
    沈念sama阅读 42,163评论 3 432
  • 文/苍兰香墨 我猛地睁开眼,长吁一口气:“原来是场噩梦啊……” “哼!你这毒妇竟也来了?” 一声冷哼从身侧响起,我...
    开封第一讲书人阅读 41,189评论 0 281
  • 序言:老挝万荣一对情侣失踪,失踪者是张志新(化名)和其女友刘颖,没想到半个月后,有当地人在树林里发现了一具尸体,经...
    沈念sama阅读 47,727评论 1 328
  • 正文 独居荒郊野岭守林人离奇死亡,尸身上长有42处带血的脓包…… 初始之章·张勋 以下内容为张勋视角 年9月15日...
    茶点故事阅读 39,741评论 3 350
  • 正文 我和宋清朗相恋三年,在试婚纱的时候发现自己被绿了。 大学时的朋友给我发了我未婚夫和他白月光在一起吃饭的照片。...
    茶点故事阅读 41,852评论 1 358
  • 序言:一个原本活蹦乱跳的男人离奇死亡,死状恐怖,灵堂内的尸体忽然破棺而出,到底是诈尸还是另有隐情,我是刑警宁泽,带...
    沈念sama阅读 37,444评论 5 352
  • 正文 年R本政府宣布,位于F岛的核电站,受9级特大地震影响,放射性物质发生泄漏。R本人自食恶果不足惜,却给世界环境...
    茶点故事阅读 43,162评论 3 341
  • 文/蒙蒙 一、第九天 我趴在偏房一处隐蔽的房顶上张望。 院中可真热闹,春花似锦、人声如沸。这庄子的主人今日做“春日...
    开封第一讲书人阅读 33,569评论 0 25
  • 文/苍兰香墨 我抬头看了看天上的太阳。三九已至,却和暖如春,着一层夹袄步出监牢的瞬间,已是汗流浃背。 一阵脚步声响...
    开封第一讲书人阅读 34,735评论 1 278
  • 我被黑心中介骗来泰国打工, 没想到刚下飞机就差点儿被人妖公主榨干…… 1. 我叫王不留,地道东北人。 一个月前我还...
    沈念sama阅读 50,443评论 3 383
  • 正文 我出身青楼,却偏偏与公主长得像,于是被迫代替她去往敌国和亲。 传闻我的和亲对象是个残疾皇子,可洞房花烛夜当晚...
    茶点故事阅读 46,931评论 2 368