28 #include <CCfits/CCfits>
40 <<
"has started is not allowed";
49 <<
"has started is not allowed";
58 <<
"has started is not allowed";
71 fits = std::make_shared<CCfits::FITS>(filename, CCfits::RWmode::Write);
75 auto& info = *table.getColumnInfo();
78 for (
size_t column_index = 0; column_index < info.size(); ++column_index) {
79 column_name_list.
push_back(info.getDescription(column_index).name);
80 column_unit_list.
push_back(info.getDescription(column_index).unit);
85 CCfits::HduType hdu_type = (
m_format ==
Format::BINARY) ? CCfits::HduType::BinaryTbl : CCfits::HduType::AsciiTbl;
87 auto extension_map = fits->extension();
88 auto extension_i = extension_map.find(
m_hdu_name);
89 bool new_hdu = (extension_i == extension_map.end() || extension_i->second->version() != 1);
91 CCfits::Table* table_hdu;
93 table_hdu =
dynamic_cast<CCfits::Table*
>(extension_i->second);
95 table_hdu = fits->addTable(
m_hdu_name, 0, column_name_list, column_format_list, column_unit_list, hdu_type);
98 for (
size_t column_index = 0; column_index < info.size(); ++column_index) {
99 auto& desc = info.getDescription(column_index).description;
100 table_hdu->addKey(
"TDESC" +
std::to_string(column_index + 1), desc,
"");
102 auto shape_str =
getTDIM(table, column_index);
103 if (!shape_str.empty()) {
104 table_hdu->addKey(CCfits::Column::TDIM() +
std::to_string(column_index + 1), shape_str,
"");
109 table_hdu->writeComment(c);
123 fits = std::make_shared<CCfits::FITS>(
m_filename, CCfits::RWmode::Write);
127 auto& info = *table.getColumnInfo();
128 for (
size_t column_index = 0; column_index < info.size(); ++column_index) {