Python 中二维 Numpy 数组的等值线图

Question

my aim is to get a contour plot in Python for an (100,100) array imported and created with Fortran.

I imported the array from Fortran in the following way:

x=np.linspace(0.02,10,100), y=np.linspace(0.47,4,100)

f = (np.fromfile(('/path/result.dat'
), dtype=np.float64).reshape((len(x), len(y)), order="F"))

So the result is dependent from x and y and gives a value for every combination of x and y. How can I create a corresponding contour plot? So far what I tried was:

X, Y= np.meshgrid(x, y)


plt.contourf(X, Y, f, colors='black')
plt.show()

But the resulting contour plot shows values that dont make sense. I also tried imshow() but it did not work. If you could help me, I will be very grateful!

Answer 1

The arrangement of X , Y , and f plays a role here. Without looking at how the result.dat was generated, though, it is difficult to answer this question. Intuition tells me, the values of f(x,y) may not match with the meshgrid.

The improper values might be arising because the values of X and Y don't correspond to the values of f . Try order = "C" or order = "A" . Also, your x and y should really be defined before reshaping the data.

x=np.linspace(0.02,10,100)
y=np.linspace(0.47,4,100)
f = np.fromfile(('/path/result.dat'), dtype=np.float64).reshape((len(x), len(y)), order="<>")

Maybe try reordering X and Y if this doesn't work.