在本教程中,此方法有效:
var a = 1;
var b = 2;
function () {
return !(a > b);
}();
但是,一旦我插入自己的文件,它就会停止工作。我正在使用google collab。这段代码是正确的-一直有效
from matplotlib import pyplot as plt
from astropy.io import fits
from astropy.wcs import WCS
from astropy.utils.data import get_pkg_data_filename
filename = get_pkg_data_filename('tutorials/FITS-images/HorseHead.fits')
hdu = fits.open(filename)[0]
wcs = WCS(hdu.header)
fig = plt.figure()
fig.add_subplot(111, projection=wcs)
plt.imshow(hdu.data, origin='lower', cmap=plt.cm.viridis)
然后我在教程的初始代码中使用它
from google.colab import drive
drive.mount('/content/gdrive')
root_path = '/content/gdrive/My Drive/research/HFI_SkyMap_353-psb-field-IQU_2048_R3.00_full.fits'
我尝试转换为浮点数-不适用于我。然后我尝试逐列转换数据(只有3个)
filename = root_path
hdu = fits.open(filename)[1]
没用。
错误消息:
a = np.array(hdu.data['I_STOKES'], dtype = float)[1:]
b = np.array(hdu.data['Q_STOKES'], dtype = float)[1:]
c = np.array(hdu.data['U_STOKES'], dtype = float)[1:]