График Smooth Contour полностью заполнен

У меня есть данные со значениями (X, Y, Z). Я попытался построить график плотности со значениями Z для интенсивности. Однако сюжет, который я получаю, не является гладким и имеет многогранник, то есть не полностью заполнен.

Ниже приведен код с данными.

но я хочу получить ровный и полностью заполненный участок


import numpy as np
from scipy.interpolate import griddata
import matplotlib.pyplot as plt
import xlrd
location = "~/Desktop/Data.xlsx"
data = xlrd.open_workbook(location)
sheet = data.sheet_by_index(0)
sample=2000

x=np.array(sheet.col_values(0))[0:sample]
y=np.array(sheet.col_values(1))[0:sample]
z=np.hamming(9000)[0:sample]
print z
def plot_contour(x,y,z,resolution = 500,contour_method='cubic'):
    resolution = str(resolution)+'j'
    X,Y = np.mgrid[min(x):max(x):complex(resolution),   min(y):max(y):complex(resolution)]
    points = [[a,b] for a,b in zip(x,y)]
    Z = griddata(points, z, (X, Y), method=contour_method)
    return X,Y,Z

X,Y,Z = plot_contour(x,y,z,resolution = 500,contour_method='linear')

plt.style.context("seaborn-deep")

plt.contourf(X,Y,Z)
plt.colorbar()
plt.show()

Это результат:

введите здесь описание изображения

Это то, чего я хочу добиться с помощью контурного графика:

введите здесь описание изображения


python matplotlib contourf
person gfdsal    schedule 03.10.2019    source источник

Ответы (1)


arrow_upward
1
arrow_downward

plt.contourf() здесь не главная проблема, он просто работает с имеющимися у него данными. Проблема в линейной интерполяции в scipy.interpolate.griddata().

Я рекомендую не использовать griddata, а вместо этого использовать один из следующих методов:

Все эти методы будут заполнять сетку. Если вы построите результат с помощью plt.imshow(), вы получите тип графика, который вы показываете в своем вопросе, который не является графиком plt.contourf().

Вот демо-блокнот, показывающий все эти подходы ( включая griddata).

person kwinkunks    schedule 04.10.2019
comment
да, я использовал scipy.interpolate.Rbf(), но, видимо, график был пикселизирован (независимо от разрешения. Проблема в том, что данных действительно недостаточно, поэтому я смог последовать вашему совету и изучить интерполяцию GP, чтобы получить больше данных , Получил то, что искал. Ура, что вернулся. - person gfdsal; 05.10.2019