作者znmkhxrw (QQ)
看板Python
標題[問題] cv2.circle與記憶體位置
時間Fri May 11 02:30:28 2018
想請問一下，為什麼以下程式碼會錯誤：
import cv2
img_rgb = cv2.imread("image.jpg")[:,:,::-1]
cv2.circle(img_rgb, (616,44),4,[255, 0, 0], thickness=-1)
TypeError: Layout of the output array img is incompatible with cv::Mat
(step[ndims-1] != elemsize or step[1] != elemsize*nchannels)
但是！我只要
加入.copy()就對了 即
....[:,:,::-1].copy()
來龍去脈如下
==========================================================
首先舉個例子釐清一件事情：
a = [1,2]
b = a[::-1]
則 b 就是 [2,1]，而且記憶體位置不同！ id(a) != id(b)
邏輯就是把 b 指向某個
不同於a的記憶體位置，值為a = [1,2][::-1]
因此跟 b = a[::-1].copy() 應該是一樣的
再來，回到原始問題
img_rgb = cv2.imread("image.jpg")[:,:,::-1]
^^^^^^^ ^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
將這變數 值為cv2.imread("image.jpg")[:,:,::-1]
指向 只是cv2.imread載入進來是bgr順序，所以用[:,:,::-1]變成rgb順序而已
cv2.circle(img_rgb, (616,44),4,[255, 0, 0], thickness=-1)
^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
在img_rgb上某個位置加入某個大小的圓點
但就是
錯的！
不過
在[:,:,::-1]後面加個.copy()就對了
---------------------------------------------------------
debug三個多小時終於找到問題在這....
不過完全違背我對 iterator[index] 的認知
照理說
iterator[index] 就是造一份新的記憶體位置
但是在cv2.circle到底發生什麼事！？
很想知道原因QQ 謝謝解惑！
--
※ 發信站: 批踢踢實業坊(ptt.cc), 來自: 219.68.160.241
※ 文章網址: https://webptt.com/m.aspx?n=bbs/Python/M.1525977030.A.530.html
※ 編輯: znmkhxrw (219.68.160.241), 05/11/2018 02:31:26
1F：推 clang: numpy 預設是回傳一個view 而不是一份拷貝，然後cv2底層 05/11 04:00
2F：→ clang: 會去檢查記憶體有沒有對齊，所以在沒有copy()的時候他的 05/11 04:00
3F：→ clang: stride是反的，會過不了檢查 05/11 04:01
a=np.array([1,2])
b=a[::-1]
print(id(a))
print(id(b))
不一樣耶@@
c大是這個意思嗎??
4F：推 niki25672729: a =[300, 400] 05/11 11:00
5F：→ niki25672729: b = a[::-1] 05/11 11:01
6F：→ niki25672729: b[0] is a[1] => True 05/11 11:01
7F：→ niki25672729: 不知道跟這個有沒有關係？ 05/11 11:01
我還不太熟 "is" 跟 "==" 的關係
不過就你舉的例子來說 以我以前的sense a[1] == 400, b=[400,300] 所以 b[0]==400
相等好像很合理??
8F：→ laputaflutin: is比id, == 只是比較值是否相同 05/11 11:54
SOGA 謝謝
9F：推 bibo9901: a[::-1]其實是 a.__getitem__(slice(None,None,-1)) 05/11 12:54
10F：→ bibo9901: 具體怎麼實現是看a的底層實作. numpy回傳的是不同的view 05/11 12:56
11F：→ bibo9901: ,也就是表面上是ndarray, 其實是一個C ptr + 一個stride 05/11 12:57
12F：→ bibo9901: 真正的資料並沒有改變, 只是取用的方式不同而已 05/11 12:58
13F：→ bibo9901: 所以如果cv要求資料必須是連續的或對齊的, 當然就炸了 05/11 12:59
14F：→ bibo9901: 而copy會複製出一份新的資料, 而且儲存方式是連續的 05/11 13:00
15F：→ bibo9901: 可以做個簡單實驗; a=np.random.rand(3); b=a[::-1] 05/11 13:02
16F：→ bibo9901: b[0] = 7; 然後你會看到a的內容也變了 05/11 13:03
17F：→ bibo9901: 不要拿內建list的行為做類比, 他們只是語法一樣而已 05/11 13:04
18F：→ bibo9901: 請仔細讀numpy或opencv的doc 05/11 13:04
19F：→ bibo9901: 可以用a.strides和a.__array_interface__['data'][0] 05/11 13:12
20F：→ bibo9901: 來取得stride和真實資料的記憶體位置. 05/11 13:14
好，謝謝reference，這樣就說得通了，就是原理不是我想的那樣
※ 編輯: znmkhxrw (210.242.52.37), 05/11/2018 13:39:11