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

看板Python作者 (QQ)時間6年前 (2018/05/11 02:30), 6年前編輯推噓3(3017)
留言20則, 4人參與, 6年前最新討論串1/2 (看更多)
想請問一下,為什麼以下程式碼會錯誤: 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

05/11 04:00, 6年前 , 1F
numpy 預設是回傳一個view 而不是一份拷貝,然後cv2底層
05/11 04:00, 1F

05/11 04:00, 6年前 , 2F
會去檢查記憶體有沒有對齊,所以在沒有copy()的時候他的
05/11 04:00, 2F

05/11 04:01, 6年前 , 3F
stride是反的,會過不了檢查
05/11 04:01, 3F
a=np.array([1,2]) b=a[::-1] print(id(a)) print(id(b)) 不一樣耶@@ c大是這個意思嗎??

05/11 11:00, 6年前 , 4F
a =[300, 400]
05/11 11:00, 4F

05/11 11:01, 6年前 , 5F
b = a[::-1]
05/11 11:01, 5F

05/11 11:01, 6年前 , 6F
b[0] is a[1] => True
05/11 11:01, 6F

05/11 11:01, 6年前 , 7F
不知道跟這個有沒有關係?
05/11 11:01, 7F
我還不太熟 "is" 跟 "==" 的關係 不過就你舉的例子來說 以我以前的sense a[1] == 400, b=[400,300] 所以 b[0]==400 相等好像很合理??

05/11 11:54, 6年前 , 8F
is比id, == 只是比較值是否相同
05/11 11:54, 8F
SOGA 謝謝

05/11 12:54, 6年前 , 9F
a[::-1]其實是 a.__getitem__(slice(None,None,-1))
05/11 12:54, 9F

05/11 12:56, 6年前 , 10F
具體怎麼實現是看a的底層實作. numpy回傳的是不同的view
05/11 12:56, 10F

05/11 12:57, 6年前 , 11F
,也就是表面上是ndarray, 其實是一個C ptr + 一個stride
05/11 12:57, 11F

05/11 12:58, 6年前 , 12F
真正的資料並沒有改變, 只是取用的方式不同而已
05/11 12:58, 12F

05/11 12:59, 6年前 , 13F
所以如果cv要求資料必須是連續的或對齊的, 當然就炸了
05/11 12:59, 13F

05/11 13:00, 6年前 , 14F
而copy會複製出一份新的資料, 而且儲存方式是連續的
05/11 13:00, 14F

05/11 13:02, 6年前 , 15F
可以做個簡單實驗; a=np.random.rand(3); b=a[::-1]
05/11 13:02, 15F

05/11 13:03, 6年前 , 16F
b[0] = 7; 然後你會看到a的內容也變了
05/11 13:03, 16F

05/11 13:04, 6年前 , 17F
不要拿內建list的行為做類比, 他們只是語法一樣而已
05/11 13:04, 17F

05/11 13:04, 6年前 , 18F
請仔細讀numpy或opencv的doc
05/11 13:04, 18F

05/11 13:12, 6年前 , 19F
可以用a.strides和a.__array_interface__['data'][0]
05/11 13:12, 19F

05/11 13:14, 6年前 , 20F
來取得stride和真實資料的記憶體位置.
05/11 13:14, 20F
好,謝謝reference,這樣就說得通了,就是原理不是我想的那樣 ※ 編輯: znmkhxrw (210.242.52.37), 05/11/2018 13:39:11
文章代碼(AID): #1Qz8_6Km (Python)
文章代碼(AID): #1Qz8_6Km (Python)