0
我无法使用起始点和停止点的坐标切片图像。到目前为止,我有以下代码在两个aribtray点上切片图像
hdulist = fits.open(filename)
hdr= hdulist[0].header
import numpy as np
import scipy
import matplotlib.pyplot as plt
from matplotlib.pyplot import figure, show
from astropy.io import fits
from scipy import interpolate
data=hdulist[0].data
#Make a line with "num" points
D, B = input('Enter the coordinates of the starting point:').split(',')
E, C = input("Enter the coordinates of the stopping point: ").split(',')
x0= float(D)
x1= float(E)
y0= float(B)
y1= float(C)
x = np.arange(data.shape[1])
y = np.arange(data.shape[0])
#length = int((np.hypot(x1-x0, y1-y0))) (can be used instead of num_points)
num_points = 1000
xvalues = np.linspace(x0, x1, num_points)
yvalues = np.linspace(y0, y1, num_points)
f = scipy.interpolate.interp2d(x, y, data) #default is linear
# Extract the values along the line
profile = f(xvalues, yvalues) #this gives me a 2D array, I think it needs to be 1D
#c = profile.flatten()
print(profile.shape)
'轮廓'不是线性的,而是立方体。有没有办法让我的轮廓线性化,这样我就可以在起点和终点之间的点上切分图像?我只需要制作'轮廓'1D而不是2D。
我要绘制这样的:
import numpy as np
from numpy import random
from matplotlib.pyplot import figure, show
vels = np.linspace(0, 530, len(profile))
fig = figure()
frame = fig.add_subplot(1,1,1)
frame.plot(vels, profile)
frame.set_ylabel('y-axis')
frame.set_xlabel('x-axis')
frame.grid(True)
show()
print(vels.shape)
print(profile.shape)
print(len(profile))
我的代码不工作,因为我得到的阴谋没有显示线路的片,但立方体的一个切片。
非常感谢,我的代码现在按照我想要的方式工作! – Thomas