2011-03-21 45 views
2

我正在處理衛星圖像。我有一個程序(Python)的由雲的頂部的溫度,以對比衛星圖像,該程序是如下:如何使用python matplotlib估計衛星圖像中的像素數

from pylab import * 
#import pylab imread,imshow 


sat=imread('1101092045G13I04.tif') 
imshow(sat) 
map=imread('map.tif') 
map=mean(map,2,)/3 
#contour(map,cmap=cm.gray) 

imshow((sat+map)**2,cmap=cm.gray,origin=1) 

frio=where(logical_and((418.-sat)-273.15>=-110,(418.-sat)-273.15<=-31),(418.-sat)-273.15,0) 
quente=where(logical_and((660.-sat)/2-273.15>-31,(660.-sat)/2-273.15<=40),(660.-sat)/2-273.15,0) 
imshow=((frio+quente)) 

temperatura=[-80,-70,-60,-50,-40,-30] 

#cores -> 'r','b','k','c','g','m','y' 

amar='#ffff00' 
verm='#ff3333' 
verd='#00ff00' 
lara='#ff9900' 
aaaa='#ff00ff' 

contourf((frio+quente),temperatura,transparent='true',colors=[aaaa,verm,lara,amar,verd]) 
colorbar(cmap=cm.hot,shrink=0.6,orientation='horizontal',ticks=[-80, -70, -60, -50, -40, -30]) 


show() 
#savefig('testeII.png') 

現在我需要估計的像素的數量爲溫度的每個比例。我該怎麼做?請幫幫我。

+3

順便說一句,'映射'在一個不幸的變量名稱中,因爲它會影響內置的Python函數'map'。 – sastanin 2011-03-21 12:38:09

+0

「每個溫度的等級」都是不明確的。你已經列出了6種溫度和5種顏色,但是你們似乎也根據冷熱來分組。 – Eric 2011-03-21 14:07:43

回答

6

現在我需要估計每個溫度範圍的像素數。

我想histogram函數是你所需要的。