ऑक्टेव में FITS (लचीली छवि परिवहन प्रणाली) प्रारूप

ऑक्टेव के साथ अपने परिचित को जारी रखते हुए, मैंने FITS फाइलों के साथ प्राथमिक संचालन करने की कोशिश करने का फैसला किया। जैसा कि यह निकला, "डिफ़ॉल्ट" ऑक्टेव में फिट्स के साथ काम करने के लिए कोई उपकरण नहीं हैं। लेकिन Google ने मेरी मदद की। इसलिए, Google ने मुझे sourceforge पर ऑक्टेव पृष्ठ पर लाया। मैंने जिस FITS पैकेज को डाउनलोड किया था, उसे अनपैक किया और संकलन करने का प्रयास किया। लेकिन वहाँ यह था! जब चल रहा हो ।/configure, कुछ नहीं हुआ। Src डायरेक्टरी में देखने के बाद मुझे कुछ भी दिलचस्प नहीं लगा, केवल एक प्राथमिक मेकफिल। इसकी सामग्री ने मुझे बताया कि मुझे ऑक्टेव पैकेज संकलक - मोक्टफाइल को स्थापित करने की भी आवश्यकता थी। मैनड्राइव में, यह ऑक्टेव-डेवेल पैकेज में है। आवश्यक पैकेज स्थापित करने के बाद, मैंने मेक रन किया और तीन फाइलें प्राप्त कीं:

  read_fits_image.oct save_fits_image.oct save_fits_image_multi_ext.oct 

उन्हें कहाँ रखा जाए, मुझे पता नहीं था, लेकिन

  \ _oct का पता लगाएं 

मुझे निर्देशिका में लाया

  /usr/lib/octave/3.4.2/oct/i586-mandriva-linux-gnu/ 

जहां (रूट की ओर से, बेशक) फाइलों की नकल की गई थी। ऑक्टेव को शुरू करने के बाद, मैंने जो पहली चीज़ चेक की थी कि क्या ये कमांड काम करती हैं। वे काम करते हैं। फ़ाइल फ़ाइल नाम को पढ़ने के लिए। चर चर में, आपको करने की आवश्यकता है:

var = read_fits_image("filename.fits"); 

मुख्य बात यह है कि कमांड के अंत में अर्धविराम को भूलना नहीं है, इसलिए नंबरों की धाराएं प्राप्त करने के लिए नहीं! अर्धविराम के बारे में, मुझे "याद" था जब मैं इसे लगाना भूल गया था। मुझे सुखद आश्चर्य हुआ कि, मैटलैब के विपरीत, जो उपयोगकर्ता को डेटा की अंतहीन स्ट्रीम से भर देता है, ऑक्टेव आउटपुट को पृष्ठों में तोड़ता है (जाहिर है, स्क्रीन पर प्रदर्शित करने के लिए कम उपयोग किया जाता है)। फ़ाइल के हेडर को पढ़ने के लिए, निष्पादित करें:

 [var, head] = read_fits_image("filename.fits"); 

व्हॉस कमांड से पता चलता है कि सिर 80 अक्षरों की लंबाई वाली स्ट्रिंग स्ट्रिंग है (FITS प्रारूप के प्रतिबंधों को याद रखें), जिसमें फ़ाइल का पूरा हेडर केवल पंक्ति-वार नाम / मान / टिप्पणी में विभाजित किए बिना लिखा गया है (यदि आप ऑक्टेव में हेडर के साथ काम करने की योजना बनाते हैं) , आपको यह सब स्वयं करना होगा, या आवश्यक कार्यों को FITS पैकेज में जोड़ना होगा, क्योंकि इसके लिए धनराशि है) मैंने imshow का उपयोग करके फ़ाइल को स्क्रीन पर प्रदर्शित करने का प्रयास किया है। अरे नहीं! मैंने देखा ... बिल्कुल कुछ नहीं। लेकिन imagesc ने मुझे एक चित्र दिखाया और यहां तक ​​कि कलर मैप को भी कॉलोर्मैप कमांड्स में बदल दिया। यदि आप चित्र को लघुगणक में प्रदर्शित करते हैं:

 imagesc(log(var)); 

प्रदर्शन अधिक सुखद होगा। मेरे फिटव्यू को याद करते हुए , मुझे अफसोस हुआ कि " स्क्रीबर " के रूप में gnuplot का उपयोग करना बहुत ही नर्वस-व्रैकिंग है: ऑक्टेव के लिए आवश्यक डेटा को gnnnlot को भेजने के लिए प्रतीक्षा करने में काफी समय लगता है, और तब तक gnuplot उन्हें स्क्रीन पर प्रदर्शित करता है। और अंत में, फॉर्म ऑब्जेक्ट की एक निर्देशिका में फ़ाइलों को समेटने के लिए एक साधारण सा प्रोग्राम। जो पहले नंबर से शुरू होता है और अंतिम नंबर के साथ समाप्त होता है।

 function ret = sum_fits(first, last) ret = [1 2 3]; frst = 1; i = 0; for num = [ first : last ] i++; name = sprintf("object_%04d.fit", num); %    II = read_fits_image(name); %   med = median(median(II)); %      printf("%d:\tread file %s, median = %d\n", i, name, med); fflush(1); %    ! II -= med; %  "" (bias'  ) if(frst == 1) %  frst = 0; ret = II; else ret += II; end end end 

और क्षेत्र के रोटेशन के केंद्र के निर्देशांक निर्धारित करने के लिए प्राप्त छवियों के लिए एक उदाहरण:

 II = sum_fits(1,147); a = II; a(find(a < 1)) = 1; imagesc(log(a)) 

हमें निम्न चित्र मिलते हैं:
हां, मैं लगभग भूल गया था: ऑक्टेव में भी मैट-जैसे इंटरफेस हैं, उदाहरण के लिए, QtOctave: सामान्य तौर पर, इंटरफ़ेस डिज़ाइन बहुत सभ्य है, हालांकि, यह बहुत कष्टप्रद है कि कमांड लाइन को अलग से निकाला जाता है और आपको कुछ लिखना शुरू करने के लिए इसमें अपना माउस प्रहार करने की आवश्यकता है। इसलिए, ऑक्टेव CLI बहुत बेहतर और अधिक सुविधाजनक है।