diff --git a/test.py b/test.py index 0163ab0a6b4bd155ffc5acb15bd9caf163dec62c..e69de29bb2d1d6434b8b29ae775ad8c2e48c5391 100644 --- a/test.py +++ b/test.py @@ -1,34 +0,0 @@ -from astroquery.mast import Observations -from astropy.io import fits -from astropy.visualization import simple_norm -import matplotlib.pyplot as plt -from ccdproc import CCDData - -# 查询并下载观测数据 -obs_table = Observations.query_criteria(dataproduct_type=["image"], obs_collection="HST", instrument_name="ACS/WFC") -data_products_by_obs = Observations.get_product_list(obs_table[0]) -download_info = Observations.download_products(data_products_by_obs[:1], productType="SCIENCE") - -# 加载FITS文件 -fits_file_path = download_info['obs_id', 'local_path'][0]['local_path'] -hdu_list = fits.open(fits_file_path) -hdu = hdu_list[0] - -# 将HDU对象转换为CCDData对象,便于后续处理 -ccd_image = CCDData(hdu.data, unit="adu") - -# 使用ccdproc进行基本处理,比如减去平均值 -processed_image = ccd_image.subtract(ccd_image.mean()) - -# 显示处理前后的图像 -fig, (ax1, ax2) = plt.subplots(1, 2, figsize=(15, 7)) -norm = simple_norm(hdu.data, 'sqrt', percent=99) -ax1.imshow(hdu.data, norm=norm, origin='lower') -ax1.set_title('原始图像') -norm_processed = simple_norm(processed_image.data, 'sqrt', percent=99) -ax2.imshow(processed_image.data, norm=norm_processed, origin='lower') -ax2.set_title('处理后图像') -plt.show() - -# 关闭fits文件以释放资源 -hdu_list.close() \ No newline at end of file