Python'da yoğunluk haritası nasıl çizilir?

Question

Bir 2D haritadaki düzenli aralıklı noktaların x, y değerlerini içeren bir .txt dosyası var, 3. koordinat o noktadaki yoğunluk.

4.882812500000000E-004 4.882812500000000E-004 0.9072267 
1.464843750000000E-003 4.882812500000000E-004 1.405174 
2.441406250000000E-003 4.882812500000000E-004 24.32851 
3.417968750000000E-003 4.882812500000000E-004 101.4136 
4.394531250000000E-003 4.882812500000000E-004 199.1388 
5.371093750000000E-003 4.882812500000000E-004 1278.898 
6.347656250000000E-003 4.882812500000000E-004 1636.955 
7.324218750000000E-003 4.882812500000000E-004 1504.590 
8.300781250000000E-003 4.882812500000000E-004 814.6337 
9.277343750000000E-003 4.882812500000000E-004 273.8610 

aşağıdaki komutlarla gnuplot bu yoğunluk haritası, çizmek Ne zaman:

set palette rgbformulae 34,35,0 
set size square 
set pm3d map 
splot "dens_map.map" u 1:2:(log10($3+10.)) title "Density map"` 

bana bu güzel görüntüyü verir:

Şimdi istiyorum matplotlib ile aynı sonuç.

Answer 1

Renk haritasının seçilmesi ve renk ekseninin logaritmik bir şekilde normalleştirilmesi dahil olmak üzere daha eksiksiz bir cevabı hedefliyorum.

import matplotlib.pyplot as plt 
import matplotlib.cm as cm 
from matplotlib.colors import LogNorm 
import numpy as np 
x, y, z = np.loadtxt('data.txt', unpack=True) 
N = int(len(z)**.5) 
z = z.reshape(N, N) 
plt.imshow(z+10, extent=(np.amin(x), np.amax(x), np.amin(y), np.amax(y)), 
     cmap=cm.hot, norm=LogNorm()) 
plt.colorbar() 
plt.show() 

Verilerinizin basit bir yeniden şekillendirme ile 2d dizisine dönüştürülebildiğini varsayalım. Eğer durum bu değilse, bu formdaki verileri almak için biraz daha fazla çalışmanız gerekir. Eğer veri bir ızgara üzerinde (eğer göründüğü gibi) yatıyorsa, imhow ve değil pcolormesh kullanarak daha verimli. Aşağıdaki resimde yukarıdaki kod parçacığı sonuçları, sen öyle istiyorsun oldukça yakın gelir:

Answer 2

@HYRY'den gelen yorum iyidir, ancak tam bir minimum çalışma yanıtı (resimle!) Daha iyidir.

verileri numunenin benziyorsa plt.pcolormesh

import pylab as plt 
import numpy as np 

# Sample data 
side = np.linspace(-2,2,15) 
X,Y = np.meshgrid(side,side) 
Z = np.exp(-((X-1)**2+Y**2)) 

# Plot the density map using nearest-neighbor interpolation 
plt.pcolormesh(X,Y,Z) 
plt.show() 

kullanarak, numpy komutunu numpy.genfromtext kullanarak yükleyebilirsiniz.