Re: [問題]計算10bits色溫

看板MATLAB作者YoursEver (銀英傳出webgame!?)時間8年前 (2016/08/10 12:55)推噓3(3推 0噓 1→)

留言4則, 3人參與討論串4/4 (看更多)

※ 引述《immortal1314 (Eden)》之銘言： : 恩好的 : 首先是這樣我對一個顯示面板進行量測 : 得到0~255gray level的 : 白光 W_Y W_x W_y Ex:某一階的 W_Y=400 W_x=0.3342 W_y=0.3153 : 紅光 R_x R_y Ex: R_x=0.6341 R_y=0.3412 : 綠光 G_x G_y Ex: G_x=0.3310 G_y=... : 藍光 B_x B_y Ex: ..... : 此為input : (x X y Y 大寫小寫代表的不同) : 然後從結果往回看 : 為了計算CCT 需要CIE色座標的x,y : 為了算CIE色座標的x,y 需要XYZ : 為了得到XYZ 需要算出R_Y G_Y B_Y : 於是乎目前流程是這樣的 : 輸入W的xyY 與 RGB的xy ,先算出RGB各自的Y , : 將RGB的xyY 轉換成各自對應的XYZ(XYZ才能做混色) : 混色後的 x=RX/(RX+GX+BX) : y=RY/(RY+GY+BY) : (x,y各有1021*1021*1021個組和)(省略了將0~255內插成0~1020的說明) : 再將x,y帶入CCT的公式得到色溫(K) : ******************************* : 簡單流程如下: : WxyY RGBxy-----------------> : R_Y G_Y B_Y----------------> : XYZ------------------------> : CIE座標上的xy--------------> : ----->CCT(color temperature) : ******************************* : 我要找出所有答案皆為某色溫(ex:5000K)的RGB各自的xyY : 也就是最一開始輸入的RGB_xyY 是哪些組合得到目標色溫的 : (目標色溫每一階0~255階找最接近那個這部分我還沒寫到code裡) : 遇到的問題: : 1.for 迴圈寫了3層雖然寫出了所有結果 but執行起來太慢 : 2.計算結果太大必須用實驗室的電腦才跑得起來光是xy就會吃掉14G的記憶體 : 急--> 3.即使計算完我也想不出怎麼往回找出是哪組R_x R_y R_Y得到我要的色溫 : G_x G_y G_Y : B_x B_y B_Y : 解決辦法: : 1.不用for 迴圈寫但又要能計算出所有結果 : 2.增加限制條件使它不用跑完所有的計算結果 : (1) 利用gamma2.2限制亮度 (R_Y+B_Y+G_Y,就是限制每一階的總亮度) : (2) 利用下一階不會超過上一階亮度的限制 (例如: G_Y 180階<G181階) : (3) 利用各階亮度的差不會太大限制搜索範圍 : 例如:不可能出現 R_Y的200階 + G_Y的10階 +B_Y : 相差太大的結果是不可能的 : 因為會導致下一階出現沒有值可用的情況 : ..................以上方法我都不會寫 : 所以目前是用最簡單的暴力解全部展開再取出所有目標與對應的RGBxyY : 卻發現不知道如何找出目標色溫對應的所有RGB_xyY : 非常感謝回復我的大神 QQ...我想很久還是百思不得解你的程式碼看起來好雜, 先把前面color transform的部分簡化一下, (1) 你的input data是 xyY representation, 你所需要的output是 XYZ representation, 所以, matlab下只需要靠 makecform 和 applycform 就可以完成, 假設你已經先排列好 W, R, G, B的 m-n-3 xyY 矩陣, 分別稱為: W_xyY, R_xyY, G_xyY, B_xyY; C = makecform( 'xyl2xyz' ); % convert xyY to XYZ; W_XYZ = applycform( W_xyY, C); R_XYZ = applycform( R_xyY, C); G_XYZ = applycform( G_xyY, C); B_XYZ = applycform( B_xyY, C); (2) 假設你已經知道你的目標色溫對應到某個表示法後的區間範圍以下是pseudo code: step-1: 建立一個mask來紀錄/標記目標色溫區域. 例: Indicator = zeros( size( IMG ) ); Indicator( intersect(find(IMG<upperbound), find(IMG>lowerbound)) = 1; step-2: 把目標區取出 Target_area = IMG .* Indicator; (3) 我不了解你說的: "利用各階亮度的差不會太大限制搜索範圍例如:不可能出現 R_Y的200階 + G_Y的10階 +B_Y" (4) 比較好的做法應該是把 xyz 換到 CCT approximation上, 然後直接對 CCT domain 切出你需要的值. 切法請直接套用(2). -- ※ 發信站: 批踢踢實業坊(ptt.cc), 來自: 140.109.22.216 ※ 文章網址: https://www.ptt.cc/bbs/MATLAB/M.1470804909.A.24A.html