在用pcolor進(jìn)行顯示圖片時(shí)，要調(diào)用colormap，caxis，shading進(jìn)行設(shè)置，才有可能還原圖片本來(lái)的色彩

pcolor的用法是：

pcolor(N)

其中N是一個(gè)矩陣，矩陣的行列數(shù)表示圖片有多少個(gè)色塊。忘了說(shuō)了，pcolor不是用像素來(lái)表示圖片大小，而是用色塊來(lái)表示。如果N是一個(gè)a*b的矩陣，表示這個(gè)圖片長(zhǎng)寬各位(a-1)*(b-1)個(gè)色塊。

a = [1 2 3 45 6 7 89 10 11 12];pcolor(a); shading faceted;

colormap在我的上一篇博客中有介紹http://blog.csdn.net/whoispo/article/details/49913513

這里要說(shuō)明的是：N矩陣的值無(wú)論有多大多小，都可以等比例變換到colormap中（1到colormap的行數(shù)）。fill中也有類(lèi)似的功能，imshow沒(méi)有這樣的功能，也就是imshow的矩陣的值必須在colormap中。

caxis就表示N中顏色的范圍，用[min,max]表示，對(duì)于N中的所有小于等于min的值，都映射到colormap的第一行，N中所有大于等于max的值，都映射到colormap的最后一行。一般我們用默認(rèn)的值，也就是min是N中最小的值，max是N中最大的值。

最后說(shuō)一下shading，shading有三種模式：

shading flat
shading faceted
shading interp

flat和faceted只是色塊和色塊直接顯示不顯示黑線的問(wèn)題，沒(méi)有太大區(qū)別。色塊的顏色用的是N中較小索引的顏色，也就是第(i,j)個(gè)色塊，用的是N(i,j)的顏色值。色塊內(nèi)部顏色一致，所以N中最后一行和最后一列都沒(méi)有用到。

interp是一個(gè)插值方法。色塊內(nèi)部顏色不是一致的，第(i,j)個(gè)色塊內(nèi)部顏色，用的是第N(i,j)， N(i,j+1), N(i+1,j), N(i+1,j+1)的點(diǎn)的顏色進(jìn)行插值出來(lái)的。因此會(huì)用到N中的所有數(shù)。

好了理論部分說(shuō)完了，就說(shuō)一下使用時(shí)會(huì)出現(xiàn)的問(wèn)題吧。

問(wèn)題1：
問(wèn)題描述：pcolor顯示圖片一團(tuán)黑。

[X, map] = imread('corn.tif'); pcolor(X); colormap(map);

問(wèn)題原因：這是因?yàn)閙atlab默認(rèn)采用shading faceted。圖片中的黑色，實(shí)際上是很多黑線。你放大看一看就明白了。

問(wèn)題解決方法：
把黑線去掉就行了，有兩種方法

solution1:

[X, map] = imread('corn.tif'); h = pcolor(X); colormap(map); set(h, 'edgecolor', 'none');

solution2:

[X, map] = imread('corn.tif'); h = pcolor(X); colormap(map); shading flat;

問(wèn)題2：圖片顯示還是不理想，不是原圖。
問(wèn)題原因：[X, map] = imread(‘corn.tif’) 返回的map不是pcolor要用到的colormap，而是imshow用的colormap。我上面說(shuō)了，imshow是不對(duì)矩陣做等比例變換的，imshow的矩陣中0對(duì)應(yīng)map中的第一行，1對(duì)應(yīng)第二行，……

你可以這樣試一試

[X, map] = imread('corn.tif'); imshow(X, map);

這樣就很完滿(mǎn)地顯示了圖片。但是這里我想用pcolor顯示圖片，不用imshow。這個(gè)問(wèn)題的根本原因其實(shí)就是上面返回的map矩陣，后面許多都是沒(méi)用的行，也就是都是0的行。這個(gè)問(wèn)題對(duì)于imshow不會(huì)出現(xiàn)問(wèn)題，但是對(duì)于pcolor就會(huì)出問(wèn)題。

因此保留map中前面的非零行作為pcolor的colormap，就沒(méi)有問(wèn)題了。

[X, map] = imread('corn.tif'); map = map(1:128,:); pcolor(X); colormap(map); shading flat;

看起來(lái)長(zhǎng)寬有些不對(duì)啊，用axis equal可以了

[X, map] = imread('corn.tif'); map = map(1:128,:); pcolor(X); colormap(map); shading flat; axis equal;

其實(shí)上面的這些我認(rèn)為都不是什么太大的問(wèn)題，自己嘗試一下就明白了。下面是我認(rèn)為比較重要的問(wèn)題。

我在試上面的例子的時(shí)候，用的是shading interp。因?yàn)槲抑坝玫腸olormap都是系統(tǒng)自帶的winter,summer之類(lèi)的colormap，這些colormap變化是平滑的，所以認(rèn)為用shading interp會(huì)使得圖片顏色更加平滑，放大不會(huì)出現(xiàn)鋸齒。結(jié)果：

[X, map] = imread('corn.tif'); map = map(1:128,:); pcolor(X); colormap(map); shading interp; axis equal;

我放大之后才想明白這是什么問(wèn)題：

這個(gè)原因就是colormap不是平滑的，所以插值出的顏色也是不平滑的。上面是一個(gè)色塊的4個(gè)角點(diǎn)，matlab會(huì)在中間進(jìn)行插值，插值的顏色會(huì)索引colormap，但是你看一下這個(gè)colormap:

[X, map] = imread('corn.tif'); map = map(1:128,:); pcolor(X); colormap(map); shading interp; axis equal;x = [0, 0, 1, 1]; y = [0, 1, 1, 0]; figure; fill(x, y, [0, 0, 2, 2]); colormap(map);

我想到了一個(gè)方法也許能解決這個(gè)問(wèn)題。把這個(gè)colormap進(jìn)行排序，把它變得平滑一些。因?yàn)閏olormap變化了，所以重新修改圖片矩陣的數(shù)值。這樣也許就能用shading interp了。

我現(xiàn)在還沒(méi)有試這個(gè)方法，也許以后有時(shí)間再試。

其實(shí)如果不出問(wèn)題的話，shading interp真是個(gè)好東西。用這個(gè)方法顯示出的圖片無(wú)論放大多少都不會(huì)出現(xiàn)鋸齒模糊，就像矢量圖片一樣。

總結(jié)

以上是生活随笔為你收集整理的关于matlab中pcolor显示图片时的shading设置问题的全部?jī)?nèi)容，希望文章能夠幫你解決所遇到的問(wèn)題。

如果覺(jué)得生活随笔網(wǎng)站內(nèi)容還不錯(cuò)，歡迎將生活随笔推薦給好友。