這篇文章將為大家詳細(xì)講解有關(guān)怎么利用python實(shí)現(xiàn)AR,小編覺得挺實(shí)用的,因此分享給大家做個(gè)參考,希望大家閱讀完這篇文章后可以有所收獲。
為來賓等地區(qū)用戶提供了全套網(wǎng)頁設(shè)計(jì)制作服務(wù),及來賓網(wǎng)站建設(shè)行業(yè)解決方案。主營業(yè)務(wù)為成都網(wǎng)站建設(shè)、網(wǎng)站設(shè)計(jì)、來賓網(wǎng)站設(shè)計(jì),以傳統(tǒng)方式定制建設(shè)網(wǎng)站,并提供域名空間備案等一條龍服務(wù),秉承以專業(yè)、用心的態(tài)度為用戶提供真誠的服務(wù)。我們深信只要達(dá)到每一位用戶的要求,就會(huì)得到認(rèn)可,從而選擇與我們長期合作。這樣,我們也可以走得更遠(yuǎn)!大致步驟如下:
識(shí)別參考平面
在這一步中,我們所需要做的事就是提取參考圖像和目標(biāo)圖像的sift特征,然后使用RANSAC算法穩(wěn)健地估計(jì)單應(yīng)性矩陣。
代碼如下:
#計(jì)算特征 sift.process_image('D:輸入圖片/book_frontal.JPG', 'im0.sift') l0, d0 = sift.read_features_from_file('im0.sift') sift.process_image('D:輸入圖片/book_per.JPG', 'im1.sift') l1, d1 = sift.read_features_from_file('im1.sift') #匹配特征,并計(jì)算單應(yīng)性矩陣 matches = sift.match_twosided(d0, d1) ndx = matches.nonzero()[0] fp = homography.make_homog(l0[ndx, :2].T) ndx2 = [int(matches[i]) for i in ndx] tp = homography.make_homog(l1[ndx2, :2].T) model = homography.RansacModel() H, inliers = homography.H_from_ransac(fp, tp, model)
由上面代碼可得到單應(yīng)性矩陣,它能夠?qū)⒁环鶊D像中標(biāo)記物的點(diǎn)映射到另一幅圖像中的對(duì)應(yīng)點(diǎn)。還需要建立X-Y(Z=0)三維坐標(biāo)系,標(biāo)記物在Z=0平面上,原點(diǎn)在標(biāo)記物的某個(gè)位置上。
從單應(yīng)性推導(dǎo)出從參考面坐標(biāo)系到目標(biāo)圖像坐標(biāo)系的轉(zhuǎn)換
在進(jìn)行坐標(biāo)轉(zhuǎn)換之前,為了檢驗(yàn)單應(yīng)性矩陣結(jié)果的正確性,需要將一些三維物體放置在目標(biāo)圖像上,本實(shí)驗(yàn)使用了一個(gè)立方體。產(chǎn)生立方體的代碼如下:
def cube_points(c, wid): p = [] p.append([c[0]-wid, c[1]-wid, c[2]-wid]) p.append([c[0]-wid, c[1]+wid, c[2]-wid]) p.append([c[0]+wid, c[1]+wid, c[2]-wid]) p.append([c[0]+wid, c[1]-wid, c[2]-wid]) p.append([c[0]-wid, c[1]-wid, c[2]-wid]) p.append([c[0]-wid, c[1]-wid, c[2]+wid]) p.append([c[0]-wid, c[1]+wid, c[2]+wid]) p.append([c[0]+wid, c[1]+wid, c[2]+wid]) p.append([c[0]+wid, c[1]-wid, c[2]+wid]) p.append([c[0]-wid, c[1]-wid, c[2]+wid]) p.append([c[0]-wid, c[1]-wid, c[2]+wid]) p.append([c[0]-wid, c[1]+wid, c[2]+wid]) p.append([c[0]-wid, c[1]+wid, c[2]-wid]) p.append([c[0]+wid, c[1]+wid, c[2]-wid]) p.append([c[0]+wid, c[1]+wid, c[2]+wid]) p.append([c[0]+wid, c[1]-wid, c[2]+wid]) p.append([c[0]+wid, c[1]-wid, c[2]-wid] return array(p).T
先計(jì)算出照相機(jī)的標(biāo)定矩陣,就可以得出兩個(gè)視圖間的相對(duì)變換
代碼如下:
#計(jì)算照相機(jī)標(biāo)定矩陣,使用圖像的分辨率為747*1000 K = my_calibration((747, 1000)) #位于邊長為0.2,Z=0平面上的三維點(diǎn) box = cube_points([0, 0, 0.1], 0.1) #投影第一幅圖像上底部的正方形 cam1 = camera.Camera(hstack((K, dot(K, array([[0], [0], [-1]]))))) #底部正方形上的點(diǎn) box_cam1 = cam1.project(homography.make_homog(box[:, :5])) #使用H將點(diǎn)變換到第二幅圖像上 box_trans = homography.normalize(dot(H,box_cam1)) #從cam1和H中計(jì)算第二個(gè)照相機(jī)矩陣 cam2 = camera.Camera(dot(H, cam1.P)) A = dot(linalg.inv(K), cam2.P[:, :3]) A = array([A[:, 0], A[:, 1], cross(A[:, 0], A[:, 1])]).T cam2.P[:, :3] = dot(K, A) #使用第二個(gè)照相機(jī)矩陣投影 box_cam2 = cam2.project(homography.make_homog(box))
在圖像(像素空間)中投影我們的3D模型并繪制它。
#底部正方形的二維投影 figure() imshow(im0) plot(box_cam1[0, :], box_cam1[1, :], linewidth=3) title('2D projection of bottom square') axis('off') #使用H對(duì)二維投影進(jìn)行變換 figure() imshow(im1) plot(box_trans[0, :], box_trans[1, :], linewidth=3) title('2D projection transfered with H') axis('off') #三維立方體 figure() imshow(im1) plot(box_cam2[0, :], box_cam2[1, :], linewidth=3) title('3D points projected in second image') axis('off')
實(shí)驗(yàn)結(jié)果如下
什么是AR
AR全稱Augmented Reality,意為增強(qiáng)現(xiàn)實(shí)技術(shù)。
它是一種將真實(shí)世界信息和虛擬世界信息“無縫”集成的新技術(shù),是把原本在現(xiàn)實(shí)世界的一定時(shí)間空間范圍內(nèi)很難體驗(yàn)到的實(shí)體信息(視覺信息,聲音,味道,觸覺等),通過電腦等科學(xué)技術(shù),模擬仿真后再疊加,將虛擬的信息應(yīng)用到真實(shí)世界,被人類感官所感知,從而達(dá)到超越現(xiàn)實(shí)的感官體驗(yàn)。真實(shí)的環(huán)境和虛擬的物體實(shí)時(shí)地疊加到了同一個(gè)畫面或空間同時(shí)存在。
增強(qiáng)現(xiàn)實(shí)技術(shù),不僅展現(xiàn)了真實(shí)世界的信息,而且將虛擬的信息同時(shí)顯示出來,兩種信息相互補(bǔ)充、疊加。在視覺化的增強(qiáng)現(xiàn)實(shí)中,用戶利用頭盔顯示器,把真實(shí)世界與電腦圖形多重合成在一起,便可以看到真實(shí)的世界圍繞著它。
增強(qiáng)現(xiàn)實(shí)技術(shù)包含了多媒體、三維建模、實(shí)時(shí)視頻顯示及控制、多傳感器融合、實(shí)時(shí)跟蹤及注冊(cè)、場景融合等新技術(shù)與新手段。增強(qiáng)現(xiàn)實(shí)提供了在一般情況下,不同于人類可以感知的信息。
如何利用python實(shí)現(xiàn)AR
步驟和實(shí)現(xiàn)姿態(tài)估計(jì)無太大差別。由上述內(nèi)容計(jì)算出照相機(jī)的位置和姿態(tài),使用這些信息來放置計(jì)算機(jī)圖像學(xué)模型。這里我們放置了一個(gè)紅色的小茶壺。
在運(yùn)行代碼之前,我們需要先安裝PyGame和PyOpenGL,下載鏈接(https://www.lfd.uci.edu/~gohlke/pythonlibs/)。
核心代碼如下:
def set_projection_from_camera(K): #從照相機(jī)標(biāo)定矩陣中獲得視圖 glMatrixMode(GL_PROJECTION) glLoadIdentity() fx = K[0,0] fy = K[1,1] fovy = 2*math.atan(0.5*height/fy)*180/math.pi aspect = (width*fy)/(height*fx) near = 0.1 far = 100.0 gluPerspective(fovy,aspect,near,far) glViewport(0,0,width,height) def set_modelview_from_camera(Rt): #從照相機(jī)姿態(tài)中獲取模擬視圖矩陣 glMatrixMode(GL_MODELVIEW) glLoadIdentity() Rx = np.array([[1,0,0],[0,0,-1],[0,1,0]]) R = Rt[:,:3] U,S,V = np.linalg.svd(R) R = np.dot(U,V) R[0,:] = -R[0,:] t = Rt[:,3] M = np.eye(4) M[:3,:3] = np.dot(R,Rx) M[:3,3] = t M = M.T m = M.flatten() glLoadMatrixf(m) def draw_background(imname): #使用四邊形繪制背景圖像 bg_image = pygame.image.load(imname).convert() bg_data = pygame.image.tostring(bg_image,"RGBX",1) glMatrixMode(GL_MODELVIEW) glLoadIdentity() glClear(GL_COLOR_BUFFER_BIT | GL_DEPTH_BUFFER_BIT) glEnable(GL_TEXTURE_2D) glBindTexture(GL_TEXTURE_2D,glGenTextures(1)) glTexImage2D(GL_TEXTURE_2D,0,GL_RGBA,width,height,0,GL_RGBA,GL_UNSIGNED_BYTE,bg_data) glTexParameterf(GL_TEXTURE_2D,GL_TEXTURE_MAG_FILTER,GL_NEAREST) glTexParameterf(GL_TEXTURE_2D,GL_TEXTURE_MIN_FILTER,GL_NEAREST) glBegin(GL_QUADS) glTexCoord2f(0.0,0.0); glVertex3f(-1.0,-1.0,-1.0) glTexCoord2f(1.0,0.0); glVertex3f( 1.0,-1.0,-1.0) glTexCoord2f(1.0,1.0); glVertex3f( 1.0, 1.0,-1.0) glTexCoord2f(0.0,1.0); glVertex3f(-1.0, 1.0,-1.0) glEnd() glDeleteTextures(1) def draw_teapot(size): #在原點(diǎn)處繪制紅色茶壺 glEnable(GL_LIGHTING) glEnable(GL_LIGHT0) glEnable(GL_DEPTH_TEST) glClear(GL_DEPTH_BUFFER_BIT) glMaterialfv(GL_FRONT,GL_AMBIENT,[0,0,0,0]) glMaterialfv(GL_FRONT,GL_DIFFUSE,[0.5,0.0,0.0,0.0]) glMaterialfv(GL_FRONT,GL_SPECULAR,[0.7,0.6,0.6,0.0]) glMaterialf(GL_FRONT,GL_SHININESS,0.25*128.0) glutSolidTeapot(size) width,height = 1000,747 def setup(): pygame.init() pygame.display.set_mode((width,height),OPENGL | DOUBLEBUF) pygame.display.set_caption("OpenGL AR demo")
實(shí)驗(yàn)結(jié)果如圖:
運(yùn)行代碼遇到的錯(cuò)誤及解決方案
錯(cuò)誤如下:An error ocurred while starting the kernelfreeglut ERROR: Function called without first calling ‘glutInit'.
原因:經(jīng)大神指點(diǎn)得知這個(gè)錯(cuò)誤是freeglut和glut共存的緣故,它們倆定義了相同的方法,這個(gè)是動(dòng)態(tài)鏈接庫的重疊問題,我的在ana\Lib\site-packages\OpenGL\DLLS文件夾里面。
你需要?jiǎng)h除freeglut.vc15.dll這個(gè)文件。我這是已經(jīng)刪除完的樣子。
關(guān)于“怎么利用python實(shí)現(xiàn)AR”這篇文章就分享到這里了,希望以上內(nèi)容可以對(duì)大家有一定的幫助,使各位可以學(xué)到更多知識(shí),如果覺得文章不錯(cuò),請(qǐng)把它分享出去讓更多的人看到。
另外有需要云服務(wù)器可以了解下創(chuàng)新互聯(lián)scvps.cn,海內(nèi)外云服務(wù)器15元起步,三天無理由+7*72小時(shí)售后在線,公司持有idc許可證,提供“云服務(wù)器、裸金屬服務(wù)器、高防服務(wù)器、香港服務(wù)器、美國服務(wù)器、虛擬主機(jī)、免備案服務(wù)器”等云主機(jī)租用服務(wù)以及企業(yè)上云的綜合解決方案,具有“安全穩(wěn)定、簡單易用、服務(wù)可用性高、性價(jià)比高”等特點(diǎn)與優(yōu)勢(shì),專為企業(yè)上云打造定制,能夠滿足用戶豐富、多元化的應(yīng)用場景需求。
分享文章:怎么利用python實(shí)現(xiàn)AR-創(chuàng)新互聯(lián)
網(wǎng)址分享:http://aaarwkj.com/article24/ppcje.html
成都網(wǎng)站建設(shè)公司_創(chuàng)新互聯(lián),為您提供網(wǎng)站策劃、軟件開發(fā)、用戶體驗(yàn)、外貿(mào)建站、微信公眾號(hào)、網(wǎng)站排名
聲明:本網(wǎng)站發(fā)布的內(nèi)容(圖片、視頻和文字)以用戶投稿、用戶轉(zhuǎn)載內(nèi)容為主,如果涉及侵權(quán)請(qǐng)盡快告知,我們將會(huì)在第一時(shí)間刪除。文章觀點(diǎn)不代表本網(wǎng)站立場,如需處理請(qǐng)聯(lián)系客服。電話:028-86922220;郵箱:631063699@qq.com。內(nèi)容未經(jīng)允許不得轉(zhuǎn)載,或轉(zhuǎn)載時(shí)需注明來源: 創(chuàng)新互聯(lián)
猜你還喜歡下面的內(nèi)容