[問題] cv2.circle與記憶體位置

[znmkhxrw (QQ)]

想請問一下，為什麼以下程式碼會錯誤：

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://www.ptt.cc/bbs/Python/M.1525977030.A.530.html
※ 編輯: znmkhxrw (219.68.160.241), 05/11/2018 02:31:26
推 clang: numpy 預設是回傳一個view 而不是一份拷貝，然後cv2底層
→ clang: 會去檢查記憶體有沒有對齊，所以在沒有copy()的時候他的
→ clang: stride是反的，會過不了檢查
a=np.array([1,2])
b=a[::-1]
print(id(a))
print(id(b))

不一樣耶@@
c大是這個意思嗎??

推 niki25672729: a =[300, 400]
→ niki25672729: b = a[::-1]
→ niki25672729: b[0] is a[1]    =>    True
→ niki25672729: 不知道跟這個有沒有關係？
我還不太熟 "is" 跟 "==" 的關係

不過就你舉的例子來說 以我以前的sense a[1] == 400, b=[400,300] 所以 b[0]==400

相等好像很合理??

→ laputaflutin: is比id, == 只是比較值是否相同
SOGA 謝謝

推 bibo9901: a[::-1]其實是 a.__getitem__(slice(None,None,-1))
→ bibo9901: 具體怎麼實現是看a的底層實作. numpy回傳的是不同的view
→ bibo9901: ,也就是表面上是ndarray, 其實是一個C ptr + 一個stride
→ bibo9901: 真正的資料並沒有改變, 只是取用的方式不同而已
→ bibo9901: 所以如果cv要求資料必須是連續的或對齊的, 當然就炸了
→ bibo9901: 而copy會複製出一份新的資料, 而且儲存方式是連續的
→ bibo9901: 可以做個簡單實驗; a=np.random.rand(3); b=a[::-1]
→ bibo9901: b[0] = 7; 然後你會看到a的內容也變了
→ bibo9901: 不要拿內建list的行為做類比, 他們只是語法一樣而已
→ bibo9901: 請仔細讀numpy或opencv的doc
→ bibo9901: 可以用a.strides和a.__array_interface__['data'][0]
→ bibo9901: 來取得stride和真實資料的記憶體位置.
好，謝謝reference，這樣就說得通了，就是原理不是我想的那樣
※ 編輯: znmkhxrw (210.242.52.37), 05/11/2018 13:39:11