الفلك

تحويل إحداثيات البكسل (في ملف FITS) إلى خط الاستواء

تحويل إحداثيات البكسل (في ملف FITS) إلى خط الاستواء


We are searching data for your request:

Forums and discussions:
Manuals and reference books:
Data from registers:
Wait the end of the search in all databases.
Upon completion, a link will appear to access the found materials.

لنفترض أن لدي ملف FITS مع مجال مرصود من النجوم ، كل واحد به إحداثيات x و y بالبكسل. افترض أن الصورة محاذاة بحيث يكون المحور x هو الارتفاع الأيمن والمحور y هو الانحراف.

أنا أبحث عن برنامج تعليمي لتحويل إحداثيات البكسل لكل نجم إلى خط استوائي ، حتى لو كان التحويل تقريبيًا.

لقد كنت أبحث منذ فترة ولكن لم أجد برنامجًا تعليميًا / كتابًا / موقعًا يعدد الخطوات اللازمة لإنتاج هذا التحول.


بافتراض أن كل بكسل (مربع) له نفس المقياس الزاوي (ليس معطى إذا كان مجال الرؤية كبيرًا) بقيمة $ theta $ درجة / بكسل.

ثم الانحراف (بالدرجات): $ delta simeq delta_0 + theta y $

الصعود الأيمن (بالدرجات): $ alpha simeq alpha_0 - theta x / cos delta $

حيث $ alpha_0 $ و $ delta_0 $ هما RA و Dec عند $ x = 0 $ ، $ y = 0 $. [علامة الطرح موجودة لأن الصعود الأيمن يزداد باتجاه يسار صورة السماء.]

يصبح هذا التقريب ضعيفًا مع زيادة مجال الرؤية.

تم وصف نهج أكثر تعقيدًا ، ولكنه دقيق ، في http://gtn.sonoma.edu/data_reduction/astrometry.php


إحداثيات¶

يصف هذا القسم إحداثيات السماء المستخدمة بواسطة الأدوات العلمية. يتم الرجوع إليه من وصف تنسيقات البيانات لشرح المعنى الدقيق للإحداثيات المخزنة.

ليس لدينا قسم منفصل لأنظمة الإحداثيات العالمية (WCS) وإحداثيات البكسل والإسقاطات التي تمت تغطيتها هنا أيضًا (انظر FITS WCS و WCSLIB للمراجع).

نناقش فقط إحداثيات السماء والصورة ثنائية الأبعاد هنا ، والإحداثيات الأخرى مثل على سبيل المثال لم يتم تناول الوقت أو محور الطاقة هنا.

تم اعتماد بعض الاتفاقيات من موقع astropy.coordinates ، وهو عبارة عن غلاف بايثون لمكتبة تنسيق ووقت IAU SOFA C ، وهو التنفيذ الرسمي لمعايير IAU الخاصة بعلم الفلك الأساسي. في بعض الحالات ، يتم تقديم أمثلة التعليمات البرمجية باستخدام astropy.coordinates للحصول على قيمة مرجعية يمكن استخدامها للتحقق من حزمة برامج معينة (في حالة عدم استنادها إلى astropy.coordinates).


تعليقات

الشيء المهم الآخر في تنسيق صورة FITS هو أنه يمكنه تسجيل مكان الصورة في السماء والإسقاط والحجم. يسمح ذلك بفتح الصورة في برامج مثل Aladin وبيانات الكتالوج المناسبة متراكبة مما يسمح بتحديد الكائنات في الصورة. أسهل طريقة للحصول على إحداثيات WCS في الصورة هي استخدام خدمة حل الألواح مثل تلك التي تقدمها Astrometry.net.
يتمثل الجانب السلبي لصور FITS في أن كيفية تسجيل الصور الملونة لا يتم تضمينها في المعيار ، ويجب أن تسجل البرامج المعلومات كمكعب مستوي 3 ولكن التفسير لكل مستوى لا يكون دائمًا هو نفسه. ملفات FITS كبيرة أيضًا ولكن يمكن ضغطها باستخدام برامج ضغط بدون فقد.

يجب أن تكون مسجلا للدخول لتكتب تعليق.

ريتشارد ، بعد أن اشترى للتو كاميرا ZWO CCD مبردة ، تمثل FITs تحديًا كبيرًا للمبتدئين. ما هو برنامج (Windows) الذي تقترح استخدامه "كمشاهد" أساسي؟ وأيها تستخدم لمعالجة الصور؟

يجب أن تكون مسجلا للدخول لتكتب تعليق.

إن "المصورين الفلكيين المبتدئين" لديهم "استبعاد كامل لتنسيق ملف .FITS" لأنهم لا يمتلكون كاميرات تصور في الواقع بتنسيق .FITS ، ولكن في أحسن الأحوال ، RAW. فلماذا المفاجأة بـ "الفصل الكامل"؟ إذا كان "المصورون الفلكيون المبتدئون" يمتلكون فقط الكاميرات التي تصور RAW ، فلماذا لا يرفضون FITS تمامًا حتى يحين الوقت الذي يخرجون فيه ويشترون كاميرا متخصصة قادرة على ملفات. لقد التقطت العديد من الصور الفلكية باستخدام كاميرا Canon Stock ووجدت أن هناك مشكلات أكبر من مجرد نوع الملف: المحاذاة القطبية الضعيفة ، والتظليل ، والبكسلات الميتة ، والبكسلات الساخنة ، والتتبع السيئ (وهذا هو غيظي مع نطاقات Celestron و Meade التي ليست كذلك مدهون بدرجات حرارة أقل من 32 * فهرنهايت! حتى لبنة طاقة التيار المتردد الضخمة من Celestron لم يتم تصنيفها للتعامل مع درجات الحرارة الأقل من درجة الحرارة! يبدو الأمر كما لو أن مصنعي التلسكوبات هؤلاء فوجئوا: "ماذا ، ستخرج عندما تكون درجة الحرارة أقل من درجة التجمد لعرضها أو التقاط صور ؟ نحن لم نفكر في ذلك ابدا ".

يجب أن تكون مسجلا للدخول لتكتب تعليق.

ريتشارد س. رايت الابن. مؤلف آخر

بالتأكيد ، RAW على ما يرام وليس هناك سبب للذهاب إلى FITS بالضرورة. صديقي ، أو معظم تجربتي هنا ، هي عندما ينتقلون إلى CCD أو CMOS المبرد ، وما إلى ذلك ، فإن DSLR لديها الكثير من الأجهزة الإلكترونية والذكاء ولديها بالفعل تنسيق ملف. لا تحفظ الكاميرات العلمية الصور من أي نوع في الواقع ، لكن البرنامج الموجود على جهاز الكمبيوتر الخاص بك يتحكم في الكاميرا ويحفظ البيانات. عادة كما يناسب. إن إجراء هذه القفزة (في تنسيقات الملفات) يكون أحيانًا حجر عثرة ، وآمل أن تساعد هذه المدونة في شرح الأساس المنطقي بشكل أفضل قليلاً لأي شخص يقوم بهذا الانتقال. قد يفكر المصورون الفلكيون الذين يبدأون في تصوير RAW في استخدام FITS أيضًا بسبب بيانات META الغنية ذات الصلة الفلكية التي يمكن تخزينها.

يجب أن تكون مسجلا للدخول لتكتب تعليق.

أتفق مع فرضية المؤلف بأن المبتدئين لا يقدرون أو يستخدمون FITS. يتخذ العديد من المبتدئين باستخدام DSLRs اختصار استخدام تنسيق jpg وتجاهل تنسيق الكاميرا الخام. من الأسهل عليهم مشاركة ملفات jpg الخاصة بهم على وسائل التواصل الاجتماعي. وإذا كانت هذه هي نهاية تطلعاتهم ، فلا بأس بذلك. ولكن إذا كانوا يريدون صورًا أفضل - معايرة ومحاذاة ومتكاملة ومعالجة - فإنهم بحاجة إلى التبديل إلى تنسيق أفضل من أجل معالجة أكثر تعقيدًا لصورتهم. عادة ما يكون للمصورين الذين يستخدمون برامج الصور الفلكية المخصصة خيار إنشاء ملفات FITS مباشرة من الكاميرا الخاصة بهم. تعتبر البيانات الإضافية الموجودة داخل الملف ذات قيمة لأدوات معالجة الصور. بالنسبة لي ، الاختبار الحمضي هو أغلفة المجلات: كم عدد مجلات علم الفلك التي تحتوي على أغلفة بها صور موجودة فقط بتنسيق jpg؟


المركز الوطني للملاحة الجوية وإدارة الفضاء

يحتوي الإصدار الحالي من FTOOLS على المهام التالية وتنفيذ نصوص Perl النصية ضمن الحزم الفرعية الثلاثة عشر. يمكن العثور على مساعدة إضافية لكل مهمة من خلال مرفق مساعدة IRAF أو مع مهمة Host FHELP بمجرد تحميل حزمة FTOOLS.

  • add_sis_rmf - [P] يُبلغ عن النسبة المئوية للتهم في ملف evtfile في كل شريحة
  • addascaspec - [P] دمج الملفات الطيفية ASCA
  • ascaarf - يولد ملف ASCA ARF لملف PHA للإدخال
  • ascaeffmap - يولد ملف صورة كاشف ASCA EFF لملف PHA للإدخال
  • ascaexpo - إنشاء خريطة تعرض لرصد ASCA.
  • ascalin - إنتاج منتجات بيانات الملف العلمي التي تمت معايرتها من ASCA
  • ascascreen - [P] ينظف ملفات ASCA العلمية الخام
  • ascaray - برنامج تتبع أشعة المرآة ASCA
  • ascatape - اقرأ ملفات FITS من شرائط توزيع بيانات ASCA المقدمة من NASA / GSFC / ADF.
  • الجانب - لحساب زوايا أويلر لأفضل نقطة عرض
  • الموقف - حساب المعلومات الموجزة لملف موقف ASCA
  • cleansis - الكشف عن قيم بكسل ASCA SIS الشاذة وإخراج ملف الأحداث التي تم تنظيفها بملحق ملحق "بكسل ساخن"
  • تصحيح RDD - قم بتطبيق تصحيح RDD على بيانات ASCA SIS FAINT أو BRIGHT أو BRIGHT2
  • dfe2mkf - أدخل معلومات إخراج FAINTDFE في ملف MKF الحالي
  • باهتة - تحويل بيانات الوضع الخافت لـ ASCA إلى تنسيق الوضع الساطع
  • faintdfe - حساب خطأ الإطار المظلم لبيانات وضع SIS FAINT
  • fasttime - يصحح أوقات بيانات وضع SIS FAST
  • fastwmap - إضافة WMAP إلى ملف طيفي لوضع FAST.
  • ؛ ghkcurve - استخراج منحنى الضوء من ملف التدبير المنزلي ASCA GIS.
  • ghkdump - عرض معلمات التدبير المنزلي GIS لملف GIS HK
  • gisrti - حساب القيمة الثابتة لوقت الصعود في ملف علمي لوضع PH GIS
  • gqaplot - أداة التحليل السريع لنظم المعلومات الجغرافية من ASCA
  • ldeadtime - إجراء تصحيح الوقت المحدد لمنحنيات ضوء FITS.
  • mkdtime - حساب كسور الوقت المحدد لنظام المعلومات الجغرافية
  • modalid - اطبع التكوين المشروط لمركبة الفضاء
  • sisclean - رسم رسم بياني للتردد لملف صورة SIS وتطبيق قطع حد الشدة
  • sisgbr - [P] تناسب QDP نسب التفرع المشتقة من بيانات FITS
  • sishist - التردد الرسومي لصورة تم إنشاؤها من ملف حدث SIS
  • sispi - يملأ عمود PI لملفات أحداث ASCA SIS.
  • sisrmg - إنشاء مصفوفة استجابة أداة SIS
  • sqaplot - أداة التحليل السريع ASCA SIS
  • timeconv - تحويل من قمر صناعي إلى مركز باري (استخدم بحذر في SUN)

مهام قاعدة بيانات المعايرة:

  • addarf - إضافة ملفين أو أكثر من ملفات الاستجابة المساعدة (ARFs)
  • addrmf - إضافة اثنين أو أكثر من RMFs معًا
  • ؛ brcaldb - تصفح مقتنيات بيانات Caldb
  • calcrpsf - حساب دالة انتشار النقطة الشعاعية
  • caldbflag - يقوم بتبديل حالة علم ملفات المعايرة في CIF
  • caldbinfo - للتحقق من تثبيت CALDB المحلي بشكل صحيح
  • chkcif - فحص الإدخالات في ملف فهرس المعايرة
  • cmppha - يحول TYPE II إلى ملف TYPE I OGIP pha
  • col2img - تحويل مجموعة بيانات استجابة ميزاء إلى صورة
  • crcif - ينشئ ملف فهرس معايرة فارغًا
  • dmprmf - يعرض ملف الاستجابة القياسي OGIP FITS
  • gcorpha - يعيد تعيين قنوات PHA من امتداد الطيف.
  • gcorrmf - يعيد تعيين مصفوفة إعادة توزيع الكاشف في مساحة القناة
  • genrsp - مولد استجابة طيفية عام
  • lstgood - قائمة & quotgood & quot الملفات الموجودة في ملف فهرس المعايرة
  • marfrmf - يضاعف ملف استجابة ثانوي بمصفوفة redist
  • ؛ mkcaldb - إنشاء قاعدة بيانات معايرة
  • ؛ mkcaldir - إنشاء بنية دليل Caldb
  • ؛ mkcalinit - إنشاء ملف caldbinit
  • mudcif - العديد من التحديثات لملف فهرس المعايرة (CIF)
  • quzcif - يستجوب كالدب لتحديد موقع مجموعة البيانات
  • rbnrmf - يضغط ملف RMF فعليًا
  • rbnrpsf - Rebins دالة انتشار النقطة الشعاعية
  • rpsf2eef - يحول مجموعة بيانات PSF (RPSF) الشعاعية إلى مجموعة بيانات دالة طاقة محاطة بدائرة شعاعية (REEF).
  • rpsfqdp - يحول ملف ملف التعريف الشعاعي بتنسيق OGIP إلى ASCII QDP
  • rsp2rmf - يحول مصفوفات استجابة تنسيق SF إلى تنسيق OGIP FITS
  • st2rpsf - يقرأ ملف STWFITS ويكتب ملف OGIP RPSF
  • stcal - يخزن ملفًا واحدًا أو أكثر من ملفات المعايرة في Caldb
  • stw2pha - يحول ملف stw PHA إلى تنسيق OGIP FITS
  • udcif - يضيف ملف معايرة إلى ملف CIF

مهام FITS الخاصة بأينشتاين:

  • sssarf - يقوم بإنشاء ملف استجابة مساعدة (ARF) لبيانات Einstein SSS PHA
  • sssbck - يبني ملفات PHA الخلفية لـ EINSTEIN SSS لمصدر معين
  • cdummyftool - & quotHello World & quot C نموذج / نموذج
  • fdummyftool - & quotHello World & quot نموذج / نموذج Fortran ftool
  • f2dhisto - قم بعمل رسم بياني ثنائي الأبعاد من عمودين في جدول
  • fadapt - قم بتحويل صورة مناسبة باستخدام مرشح قبعة علوية متكيف
  • farith - إجراء العمليات الحسابية على صورتين FITS
  • fboxcar - قم بلف صورة مناسبة باستخدام مرشح صندوق السيارة مستطيل الشكل
  • fcarith - قم بإجراء العمليات الحسابية على صورة FITS باستخدام ثابت
  • fgauss - لف صورة FITS بوظيفة جاوس بيضاوية الشكل
  • fim2lst - تحويل صورة ثنائية الأبعاد إلى قائمة بكسل (معكوس f2dhisto)
  • fimconcat - التقط صورتين وضعهما جنبًا إلى جنب كصورة واحدة
  • fimgbin - Rebin صورة مناسبة باستخدام حجم حاوية مستطيلة
  • fimgdmp - تفريغ محتويات صورة FITS في ملف ASCII
  • fimgmerge - دمج أي عدد من الصور في صورة معينة
  • fimgstat - حساب الإحصائيات (الحد الأقصى ، الحد الأدنى ، إلخ) لملف الصورة
  • fimgtrim - إعادة تعيين وحدات بكسل الصورة بقيم أعلى / أقل من قيم الحد
  • florentz - قم بتحويل صورة مناسبة باستخدام وظيفة Lorentzian الإهليلجية
  • fmaskfilt - تصفية قائمة الأحداث بناءً على صورة قناع الإدخال
  • fmrgmsk - دمج أقنعة مكانية أو أكثر
  • Fregcon - محول ملف منطقة SAOimage من أداة / كاشف إلى آخر
  • fsaoi - ترجمة ملف منطقة SAOImage إلى ملف إدخال لـ fselect
  • ؛ fvec2img - إنشاء صورة من عمود المتجهات

مهام Futils (للأغراض العامة):

  • faddcol - ينسخ الأعمدة من امتداد جدول FITS إلى آخر
  • fappend - إلحاق امتداد FITS بملف FITS آخر
  • fapropos - [P] يقوم بالبحث عن الكلمات الرئيسية لـ FTOOLS
  • fcalc - حساب قيم عمود باستخدام تعبير حسابي
  • fcatdiff - يقارن أعمدة ملف يناسب ويبلغ عن اختلافات الصفوف
  • fchecksum - تحقق من أو قم بتحديث الكلمات الأساسية للمجموع الاختباري في ملف FITS
  • fcolpar - إرجاع رقم العمود لاسم ملف الإدخال واسم العمود
  • fcreate - إنشاء جدول FITS من ملفات إدخال ASCII
  • fcurve - عمل رسم بياني لمنحنى ضوئي من عمود في جدول
  • fdelcol - حذف عمود محدد في جدول FITS
  • fdelhdu - حذف امتداد من ملف FITS
  • fdelrow - حذف الصفوف المحددة في جدول FITS
  • fdump - تفريغ محتويات جدول FITS في ملف ASCII
  • fextract - انسخ امتداد FITS من ملف إلى ملف جديد
  • ffilecat - ينسخ قيم الكلمات الأساسية من قائمة ملف FITS إلى جدول FITS
  • fhelp - [P (اختياري)] عرض معلومات تعليمات FTOOLS
  • fhisto - عمل مدرج تكراري لعمود في جدول
  • ؛ fimgcreate - إنشاء صورة مصفوفة أولية FITS من ملف قالب ASCII
  • findex - إنشاء ملف فهرس لعمود جدول FITS
  • finterp - استيفاء عمود من امتداد جدول FITS إلى آخر
  • fkeypar - انسخ كلمة رئيسية في رأس FITS إلى معلمة
  • fkeyprint - عرض الكلمات الرئيسية في رأس (رؤوس) FITS
  • fkeytab - نسخ كلمة رئيسية في رأس FITS إلى عنصر جدول FITS
  • flaunch - قاذفة FTOOLS رسومية
  • flcol - سرد معلومات العمود في ملحق جدول FITS
  • flist - سرد محتويات جدول FITS في ملف ASCII
  • flookup - تصفية ملف FITS من التحديد بناءً على جدول بحث
  • fmemsort - نوع ذاكرة سريع من جدول FITS
  • fmerge - دمج الصفوف من عدة جداول FITS في جدول FITS واحد
  • ؛ fmodcomm - تعديل الكلمات الأساسية للتعليق في رأس ملف FITS
  • ؛ fmodhead - تعديل الكلمات الأساسية للرأس في ملف FITS
  • fparhelp - [P] يوفر مساعدة لمعامل محدد لمفتاح قدم معين
  • fparkey - نسخ قيمة معلمة إلى كلمة رئيسية في رأس FITS
  • fparstamp - [P] انسخ محتويات ملف .par إلى رأس ملف FITS
  • fpartab - انسخ قيمة معلمة إلى عنصر جدول FITS
  • fplot - يرسم أعمدة من ملف FITS باستخدام QDP / PLT Plot Package
  • fscript - [P] أداة كتابة السيناريو العامة
  • fselect - إنشاء جدول جديد من الصفوف المحددة في الجدول
  • fsort - فرز جدول FITS في مكانه
  • fstatistic - حساب المتوسط ​​والانحراف المعياري والحد الأدنى والحد الأقصى للعمود
  • fstruct - سرد وصفًا لهيكل ملف FITS
  • ftabcopy - نسخ أعمدة محددة من جدول FITS إلى جدول جديد
  • ftabkey - انسخ عنصر جدول FITS إلى كلمة رئيسية في رأس FITS
  • ftabpar - نسخ عنصر جدول FITS إلى قيمة معلمة
  • fv - أداة عرض ومتصفح تفاعلي بملفات GUI FITS
  • fverify - تحقق من أن ملف FITS يتوافق مع معيار FITS
  • fversion - [P] تقرير رقم إصدار FTOOLS وتاريخ الإصدار
  • pconfig - [P] تكوين ملفات المعلمات إلى الإعدادات الافتراضية للمستخدم

المهام المحددة لمرصد أشعة جاما:

  • econvpha - يحول ملفات EGRET .spec إلى ملفات XSPEC .pha
  • econvrmf - يحول ملفات EGRET .resp إلى ملفات XSPEC .rmf.
  • fbbft2pha - يحول ملف BATSE BFITS إلى ملف PHA طيفي
  • fbdrm2rmf - يحول ملف BATSE DRM إلى تنسيق RMF
  • fesdb2rdf - يقرأ ملفات بيانات i / p QVP و EXP (اختياري) ويكتب ملف بيانات O / p FITS بتنسيق قياسي OGIP
  • fevpbtoa - يحول بيانات Comptel EVP المنسقة إلى ملف ASCII

المهام الخاصة بالفيزياء الفلكية عالية الطاقة:

  • addspec - يضيف أطياف PHA وملفات bkgd ويجمع بين مصفوفات الاستجابة
  • ascii2pha - ينشئ ملف OGIP FITS PHA من مجموعة بيانات ascii
  • bincurve - ناقلات الصناديق والمخرجات كمنحنى ضوئي
  • binspec - بن أطياف في أعمدة متجهة
  • chkarf - تحقق من صلاحية الكلمات الرئيسية في ملف ARF FITS
  • chkpha - تحقق من صلاحية الكلمات الأساسية في ملف PHA / PI FITS
  • chkrmf - تحقق من صلاحية الكلمات الأساسية في ملف RMF FITS
  • ecd2pha - يحول Einstein HEASARC CD-ROM PHA إلى ملف OGIP PHA
  • extrpsf - لاستخراج مجموعة بيانات PSF (RPSF) الشعاعية من ملف الحدث.
  • fadmap - إنتاج خرائط لبيانات SAS2 أو COSB
  • fixregion - تحويل ملف منطقة SAO إلى نظام إحداثيات آخر
  • grppha - يعالج ملف OGIP القياسي PHA FITS
  • hkexpand - قم بتوسيع ملف بيانات التدبير المنزلي بتنسيق مضغوط (HK)
  • hkscale - تحجيم ملف بيانات التدبير المنزلي الخاص بـ FITS إلى قيم مادية
  • hkunexpand - لضغط ملف بيانات التدبير المنزلي (HK) ذي التنسيق الموسع
  • mathpha - ينفذ عمليات حسابية على ملفات PHA
  • mekal - إنشاء ملف جدول XSPEC لنموذج انبعاث البلازما من ميكال
  • NH - إرجاع كثافة عمود الهيدروجين لـ Ra و Dec.
  • raysmith - إنشاء ملف جدول XSPEC لنموذج انبعاث البلازما من Raymond-Smith
  • rbnpha - لضغط مجموعة بيانات PHA
  • sf2pha - يحول ملف PHA بتنسيق SF إلى تنسيق OGIP FITS
  • sky2xy - التحويل من إحداثيات سماء إلى بكسل باستخدام كلمات WCS الأساسية
  • xy2sky - التحويل من إحداثيات البكسل إلى السماء باستخدام كلمات WCS الأساسية

مهام ROSAT الخاصة بمهمة:

  • abc - تطبيق تصحيحات Barycentric على بيانات حدث ROSAT.
  • ao - يقيّم التباين الهندسي لخلفية الأشعة السينية الشمسية المتناثرة
  • bct - إنتاج جدول تصحيح Barycenter من بيانات مدار ROSAT
  • calcbgdcor - حساب التصحيح لحساب الخلفية لطيف المصدر
  • castpart - تحسب خلفية جسيمات نموذج المصبوب في إحداثيات السماء
  • كشف - يكتشف المصادر من ملف الصورة
  • extpha - [P] يستخرج ملفات طيف المصدر من ملفات PSPC _src
  • fixrev0pha - يغير الكلمات الأساسية CDELT1 و CDELT2 في ملف ROSAT PHA
  • gevt2rdf - يحول مجموعة بيانات _events الألمانية إلى تنسيق RDF _bas
  • hiback - [P] ينشئ ملف _qsrc من ملف _src
  • hriaspcor - استخدم طريقة مورس لتصحيح أخطاء العرض
  • hriexpmap - يُنشئ خريطة تعرض لمجموعة بيانات رصد ROSAT HRI
  • hrifilt - ينشئ ملف ROSAT HRI makefilter
  • hrirpsf - يحسب النظرية HRI PSF
  • img2us - يحول من تنسيق ملف الصورة الألماني إلى تنسيق RFITS الأمريكي
  • lte - يقيّم التوافق متعدد الحدود الذي أجراه روتين RATEFIT
  • mkqsrc - يضيف حقولًا إضافية لفحص الجودة لملف روسات المصدر
  • mk_allevt - [P] يجمع بين جداول الأحداث القياسية والأمبير المرفوضة لملف RDF
  • pcarf - يولد ROSAT PSPC ARF لملف PHA للإدخال
  • pcbgdcor - ينفذ [P] حساب الخلفية. لـ ROSAT (يتم تشغيل pcparpha و pcarf و fixregion و calcbgdcor)
  • pcecor - التصحيح الزمني المعتمد على الطاقة لقناة PI في ملف الأحداث
  • pcexpmap - يُنشئ خريطة تعرض لمراقبة ROSAT PSPC
  • pcfilt - يُنشئ ملف ROSAT PSPC makefilter
  • pcparpha - يولد طيف PHA في الخلفية لـ ROSAT PSPC
  • pcpicor - [P] يصحح قناة PI في ملف أحداث PSPC
  • pcrpsf - حساب نظري PSF شعاعي لـ ROSAT PSPC
  • pcsasscor - يصحح قناة PI لخطأ الكسب المكاني SASS في ملف الأحداث
  • pctcor - يطبق التصحيح الزمني على بيانات PI في ملف أحداث PSPC
  • المعدل - لإنشاء منحنيات ضوئية مرصوفة على فترات 30 ثانية للنطاقات السبعة
  • Ratefit - يناسب منحنيات الضوء PSPC في سبعة نطاقات عريضة لارتفاع النبضة
  • rosbary - [P] إجراء تصحيح مركزي لـ ROSAT. (تشغيل abc & amp bct)
  • roscc2utc - يحول وقت المركبة الفضائية إلى UTC مع MJD ، yy: mm: dd: hr: min: sec
  • src2pha - إنشاء ملف PHA من امتدادات محددة في ملف RDF src
  • srcdetect - [P] الكشف عن المصدر الذي يقوم بتشغيل xselect ثم الكشف
  • الإمالة - يقيم الخلفية الاسمية للأشعة السينية الشمسية المتناثرة
  • انفجار - إزالة رشقات من الأحداث من قائمة الأحداث المطلوبة في الوقت المناسب
  • cktime - للتحقق من أن عمود الوقت مرتب بالوقت
  • deadtime - قم بتحديث عرض ملف PHA بحلول الموعد النهائي بناءً على GTIs
  • fltime - تصفية قائمة الأحداث ضمن فترات زمنية معينة (GTIs)
  • maketime - حساب الفواصل الزمنية (GTIs) من بيانات التدبير المنزلي (HK)
  • mgtime - دمج 2 أو أكثر من ملفات الفاصل الزمني (GTI)
  • sec2time - تحويل تعويض الوقت إلى الوقت المطلق
  • ؛ time2sec - تحويل الوقت المطلق إلى إزاحة الوقت

مهام VELA5B الخاصة بمهمة:

  • fvelalc - ينشئ منحنى ضوئي مناسب لمصدر Vela 5B معزول
  • fvelgallc - ينشئ منحنيات ضوئية متعددة لمصادر في حقول مشوشة
  • fvelmap - ينشئ مجموعة من خرائط مصدر FITS لمنطقة مستطيلة
  • velabox - يسرد ملفات البيانات الخام Vela 5B اللازمة لـ FVELALC
  • velamap - يسرد ملفات البيانات الخام Vela 5B اللازمة لـ FVELMAP

مهام Xronos (حزمة تحليل السلاسل الزمنية):

  • ascii2flc - [P] ينشئ ملف FITS لإدخال xronos من ملف ASCII
  • autocor - يحسب الارتباط التلقائي لسلسلة زمنية
  • crosscor - يحسب الارتباط المتبادل لسلسلتين زمنيتين
  • earth2sun - قم بتغيير الأوقات إلى مركز barycenter للنظام الشمسي
  • efold - ينشئ منحنى ضوئيًا مطويًا ، وصلابة ولون مؤامرات
  • efsearch - يكتشف أفضل فترة مربع كاي عن طريق الطي على مدى فترات
  • flc2ascii - [P] تفريغ إدخال أو إخراج xronos في ملف ASCII الذي يمكن قراءته بواسطة QDP
  • lcmath - يطرح أو يضيف اثنين من المنحنيات المضيئة.
  • lcstats - حساب المتغيرات الإحصائية في سلسلة زمنية
  • lcurve - ينشئ منحنيات خفيفة وصلابة ورسومات ملونة
  • listdata - قائمة ملف البيانات على الشاشة
  • powspec - يخلق كثافة طيف الطاقة
  • rbf2fits - تحويل & quot؛ تنسيق المخزن المؤقت لمعدل EXOSAT & quot في FITS.
  • المخطط الزمني - لحساب دالة الانحراف الزمني لسلسلة زمنية
  • xronwin - [P] ينشئ ملف نافذة مناسب لمهام xronos

المهام المحددة لمهمة XTE:

  • addhots - يضيف منحنى ضوضاء ضوضاء لقطة إلى منحنى ضوء الإدخال.
  • addine - يضيف منحنى جيبي إلى منحنى ضوء الإدخال.
  • asmappend - لفرز ملفات ASM وإلحاقها بجهاز MASTERFILE.
  • asmchannel - [P] يفصل ملف ألوان ASM إلى ملفات منفصلة لمنحنى الضوء لكل قناة طاقة ASM.
  • asmsrcupdate - [P] تشغيل ASMAPPEND على مجموعة من ملفات XTE ASM الرئيسية والملفات الواردة الخاصة بها.
  • bemerge - دمج ملفين طيفيين للأحداث المرصودة (.pha): النوع الأول أو الثاني.
  • chantrans - [P] ترجمة نطاق القناة النسبي إلى نطاق قناة مطلق من خلال قراءة ملفات XTE المنتجة .pha.
  • edsgcor - يحول ملف PCA / EDS Gain & amp Offset من ASCII إلى FITS.
  • fakelc - ينشئ منحنى ضوئيًا زائفًا بمتوسط ​​شدة ثابت.
  • fasebin - صندوق طور XTE (ينتج & quotfB & quot الملفات).
  • fbadd - إضافة عدة ملفات fB متناسبة معًا.
  • fbfsum - يسمح بحساب متوسط ​​ملفات fB في الطور عن طريق جمع الصفوف.
  • fbssum - يحسب متوسط ​​ملفات fB في الطاقة ويقوم بتحليل الطور.
  • fbsub - اطرح خطًا أساسيًا من مدرج تكراري fB.
  • fcollect - النصف الأول من XTEFILT يجمع بيانات التدبير المنزلي والمواقف والمدار من قاعدة بيانات XFF.
  • fxbary - تطبيق ملف تصحيح منحنى الضوء barycentric (.lc) FITS.
  • grosstimefilt - [P] تصفية جميع الأوقات التي تقع خارج النطاق المحدد وتعديل الكلمات الأساسية المناسبة في ملف XTE SA أو SE.
  • hxtarf - إنشاء ملف .arf لملاحظة HEXTE معينة.
  • hxtback - [P] يقرأ ملف HEXTE Science أو التدبير المنزلي FITS وينشئ ملفات FITS جديدة بناءً على موضع المجموعة.
  • hxtdeadlc - حساب وجدولة تصحيحات الوقت المحدد لملف .lc
  • hxtdeadpha - حساب وجدولة تصحيحات الوقت المحدد لملف .pha
  • make_se - [P] قم بإنشاء ملفات بيانات XTE / SE (حدث) من قائمة ملفات الإدخال التي قد تحتوي على ملفات بيانات Good_Xenon و / أو ملفات بيانات الوضع الشفاف.
  • mkgtif - ينشئ ملف امتداد GTI لبيانات HEXTE.
  • pca_bs - [P] قم بإجراء طرح تقريبي لخلفية المنحنى الخفيف لـ PCA.
  • pcabackest - ينشئ ويقدر طيف خلفية PCA من النموذج.
  • pcaclrsp - يقوم بإنشاء منحنى ضوء مصحح ميزاء من ملف PCA .lc.
  • pcagainset - إدراج قيم كسب وإزاحة PCA / EDS في ملف PCA .pha.
  • pcarmf - إنشاء ملف .rmf لسلسلة الأنود المحددة في وحدة PCU.
  • pcarsp - [P] تشغيل الأدوات اللازمة لإنشاء مصفوفات الاستجابة لبيانات طيفية XTE PCA.
  • perdgrm - يحسب إطارًا زمنيًا لمنحنيات الضوء التي تم أخذ عينات منها بشكل غير متساو.
  • rddescr - يوزع واصف القناة من ملف XTE .pha إلى ASCII.
  • recofmi - يعيد بناء ملف FMI (فهرس رئيسي) لقاعدة بيانات XTE.
  • runpcabackest - [P] يقوم بتشغيل PCABACKEST على سلسلة من ملفات البيانات StdMode2.
  • saexpd - يقرأ ناقلات المصفوفة العلمية ويخرج ملفات FITS العامة.
  • saextrct - ينشئ منحنى ضوئيًا و / أو طيفًا من متجهات مصفوفة علمية.
  • sebitmask - يُنشئ قناع بت لتصفية بيانات SE باستخدام FSELECT.
  • seexpd - يقرأ بيانات الأحداث العلمية ويخرج ملفات FITS العامة.
  • seextrct - يُنشئ منحنى ضوئيًا و / أو طيفًا من بيانات الأحداث العلمية.
  • sachip - [P] يتم تشغيل SAEXTRCT و FPLOT مرة واحدة من سطر الأوامر.
  • saplot - [P] يستعلم المستخدم دوريًا عن ملفات SA إضافية ثم يقوم بتشغيل SAEXTRCT و FPLOT ، لـ SOF.
  • sefilter - [P] يساعد المستخدم على إنشاء تعبير يتم إدخاله في SEBITMASK ويفحص النتيجة لتحديد ما إذا كان يمكن إدخال الملف مباشرةً في SEEXTRCT للتصفية ، أو ما إذا كان يجب إدخاله أولاً إلى FSELECT.
  • seplot - [P] يستعلم بشكل دوري المستخدم عن ملفات SE إضافية ثم يشغل SEEXTRCT و FPLOT لـ SOF.
  • حدد - [P] تشغيل SEBITMASK و FSELECT لإنشاء ملف FITS.
  • timetrans - [P] ترجمة النطاق الزمني النسبي إلى نطاق زمني مطلق من خلال قراءة ملفات XTE المنتجة .lc.
  • trans2fits - ينشئ ملف FITS للحدث العلمي من بيانات الوضع الشفاف.
  • xenon2fits - ينشئ ملف FITS للحدث العلمي من بيانات Xenon الجيدة.
  • xpcaarf - لإنشاء ملف .arf لإدخال XTE PCA .pha وملف .rmf.
  • xtecol - تطبيق تصحيح ميزاء HEXTE على منحنى ضوئي.
  • xtederive - يستخدم النصف الثاني من XTEFILT ملف الإخراج من fcollect لاشتقاق كميات مفيدة إضافية لملف XTEFilter.
  • xtefilt - [P] ينشئ ملف مرشح .xfl لملاحظة معينة.
  • xtetape - اقرأ أشرطة توزيع بيانات XTE المقدمة من NASA / GSFC / ADF.
  • xteversion - [P] اطبع إصدار وتاريخ حزمة XTE الفرعية.

مهام معالجة ملف المعلمة:

  • pconfig - [P] تكوين ملفات المعلمات إلى الإعدادات الافتراضية للمستخدم
  • pget - الحصول على قيمة معلمة واحدة أو أكثر في ملف معلمة
  • plist - سرد المعلمات من ملف أو أكثر من ملفات معلمات نمط IRAF
  • pquery - احصل على قيمة المعلمة في ملف معلمة بنمط IRAF
  • pset - عيّن معلمة واحدة أو أكثر في ملف معلمة بنمط IRAF
  • ؛ Punlearn - الحصول على نسخة من ملف معلمة النظام

واجهات المستخدم الرسومية لـ FTOOLS:

  • fv - عارض ملفات FITS تفاعلي
  • flaunch - قاذفة FTOOLS
  • xselect2 - الخلف الرسومي ، يدعم xselect حاليًا تحليل XTE فقط ولكن يتم إضافة مهام أخرى.
  • xdf - XTE Data Finder

[P] == نص Perl. يجب أن يكون لديك Perl مثبتًا لاستخدام هذه الأداة

إذا كانت FTOOLS مفيدة في بحثك ، فيرجى الرجوع إلى هذا الموقع (http://heasarc.gsfc.nasa.gov/ftools) واستخدام مرجع ASCL لـ HEASoft [ascl: 1408.004] أو مرجع ASCL لورقة FTOOLs الأصلية [ ascl: 9912.002]:

بلاكبيرن ، J.K. 1995 ، في ASP Conf. سر ، المجلد. 77 ، برامج وأنظمة تحليل البيانات الفلكية الرابع ، أد. ر.أ.شو ، إتش إي باين ، وجي جيه إي هايز (سان فرانسيسكو: ASP) ، 367.


كيف أقوم بقص صورة تناسبها مع الاحتفاظ بإحداثيات العالم للتخطيط في Astropy Python؟

كانت هذه القضية تزعجني لبعض الوقت. أحاول التعامل مع كمية كبيرة من البيانات التي تكون في شكل ملف يناسب (بترتيب 11000 × 9000 بكسل). ما أحتاج إلى فعله هو إنشاء مخطط تنسيق RA / Dec "مكبّر" (بشكل مثالي باستخدام astropy.wcs) للعديد من الكائنات في السماء مع ملامح من ملف يناسبها ، وتدرج رمادي (أو خريطة حرارية من ملف آخر).

مشكلتي هي أنه عندما أقوم بتقطيع البيانات من الصورة (إلى منطقة اهتمامي) أفقد الارتباط بإحداثيات السماء. هذا يعني أن الصورة المقطعة ليست في المكان الصحيح.

لقد قمت بتعديل مثال من المستندات الفلكية لتوفير معاناة بياناتي. (ملاحظة: أريد أن تغطي الخطوط العريضة مساحة أكبر من الصورة ، مهما كان الحل لذلك يجب أن يعمل على كلا البيانات)

هذا هو الكود الذي أواجه مشكلة فيه:

حاولت استخدام أسلاك WCS الخاصة بقطعتي المقطعة لإصلاح ذلك ، لكنني لست متأكدًا مما إذا كان بإمكاني تمريرها في أي مكان!


تعريفات رأس الملف المناسب

تنسيق الملف الأصلي لصور MaxIm DL هو FITS (نظام نقل الصور المرن) ، وهو معيار مستخدم على نطاق واسع في المجتمع الفلكي. يعد هذا اختيارًا ممتازًا لجميع أنواع الصور الفلكية نظرًا لأن التنسيق يدعم بيانات 16 بت والنقطة العائمة ، ويتضمن رأسًا مرنًا وقابل للتوسيع بدرجة كبيرة. يمكن للكلمات الأساسية غير القياسية القياسية والمستخدمة على نطاق واسع أن تنقل معلومات حول الصورة بين التطبيقات. بالإضافة إلى ذلك ، يتوفر كود المصدر "C" للمجال العام لقراءة وكتابة ملفات FITS على الويب.

أحد أوجه القصور في معيار FITS هو أنه يتم توقيع تنسيق الأرقام 16 بت ، مما ينتج عنه نطاق من -32768 إلى +32767. نطاق البيانات من العديد من الكاميرات هو 0 إلى 65535. للتغلب على هذه المشكلة ، يطرح MaxIm DL 32768 من كل بكسل قبل الحفظ. يتطلب معيار FITS استخدام الكلمات الأساسية BZERO و BSCALE عندما لا تكون القيمة في المصفوفة هي القيمة المادية ، حيث يقوم MaxIm DL بتعيين BZERO إلى 32768 و BSCALE إلى 1. يجب أن تضيف التطبيقات الأخرى التي تفسر هذه الكلمات الرئيسية بشكل صحيح 32768 تلقائيًا إلى بكسل الصورة ، مما يؤدي إلى عرض الصورة الصحيح. يمكن تجاوز هذا السلوك الافتراضي باستخدام الأمر File menu Settings.

طريقة ضغط FITS المستخدمة من قبل MaxIm DL هي ملكية خاصة وغير متوافقة مع حزم البرامج الأخرى. إذا كنت بحاجة إلى نقل الصور إلى حزمة أخرى ، فيجب عليك حفظ الصور على أنها غير مضغوطة. إذا كنت بحاجة إلى تحويل عدد كبير من الصور ، فاستخدم القائمة File أمر Batch Save and Convert.

MaxIm DL متوافق مع اقتراح SBIG لإضافات الكلمات الرئيسية FITS. هذا يعني أن العديد من الكلمات الرئيسية الاختيارية وغير القياسية تتم كتابتها بطريقة تتوافق مع عدد من حزم التصوير الفلكي الأخرى. راجع http://archive.sbig.com/pdffiles/SBFITSEXT_1r0.pdf لمزيد من المعلومات.

الكلمات الرئيسية الإلزامية لـ FITS هي كما يلي:

SIMPLE - دائمًا "T" ، للإشارة إلى رأس مناسب.

BITPIX - يشير إلى تنسيق الصفيف. تتضمن الخيارات 8 بت غير موقعة (8) ، 16 بت موقعة (16) ، 32 بت موقعة (32) ، 32 بت عائم IEEE (-32) ، و 64 بت عائم IEEE (-64). التنسيق القياسي هو 16 -64 يمكن قراءته بواسطة MaxIm DL ولكن لا تتم كتابته.

NAXIS - عدد المحاور في مصفوفة البيانات. يستخدم MaxIm DL 2 للصور أحادية اللون ، و 3 للصور الملونة.

NAXIS1 - يتوافق مع المحور X.

NAXIS2 - يتوافق مع المحور ص.

NAXIS3 - قيمة موجودة فقط للصور الملونة هي 3 دائمًا (مستويات اللون الأحمر والأخضر والأزرق موجودة بهذا الترتيب).

الكلمات الأساسية الاختيارية المحددة بواسطة معيار FITS والمستخدمة في MaxIm DL:

BSCALE - يجب ضرب هذه القيمة بقيم صفيف البيانات عند قراءة ملف FITS. يكتب MaxIm DL دائمًا قيمة 1 لهذه الكلمة الأساسية.

BZERO - يجب إضافة هذه القيمة إلى قيم مصفوفة البيانات عند قراءة ملف FITS. بالنسبة لملفات الأعداد الصحيحة ذات 16 بت ، يكتب MaxIm DL 32768 (ما لم يتم تجاوزه بواسطة مربع حوار الإعدادات).

DATE-OBS - تاريخ المراقبة بتنسيق ISO القياسي 8601 (المقاسات المتوافقة مع Y2K): CCYY-MM-DDThh: mm: ss.sss. يتم استخدام التوقيت العالمي في بداية التعريض الضوئي. ملاحظة: التنسيق البديل باستخدام DATE-OBS و TIME-OBS غير مكتوب ، لكن MaxIm DL سيفسره بشكل صحيح عند قراءته. الوقت مكتوب بدقة 10 مللي ثانية. السلوك الافتراضي هو الإبلاغ عن بدء وقت المراقبة ، ولكن يمكن لسائقي الكاميرا تغيير ذلك. & # 160 اعتبارًا من الإصدار 6.24 ، يضبط برنامج تشغيل DL Imaging الوقت لنقطة منتصف التعرض. & # 160

التاريخ - يشير إلى محفوظات معالجة الصورة. يمكن تكرار هذه الكلمة الرئيسية عدة مرات حسب الضرورة.

INSTRUME - معلومات الكاميرا. دخل المستخدم أو حصل عليه من برنامج تشغيل الكاميرا.

الكائن - اسم أو رقم فهرس الكائن الذي يتم تصويره ، إذا كان متاحًا من نافذة مراقبة المرصد أو حدده المستخدم في الإعدادات.

المراقب - المعلومات التي أدخلها المستخدم اسم المراقب.

TELESCOP - معلومات أدخلها المستخدم حول التلسكوب المستخدم.

الكلمات الأساسية للإضافة التي يمكن إضافتها أو قراءتها بواسطة MaxIm DL ، اعتمادًا على التكوين الحالي للجهاز والبرامج:

AIRMASS - طول المسار البصري النسبي عبر الغلاف الجوي.

ظروف مرصد AOCAMBT - ASCOM & # 160– درجة الحرارة المحيطة بالدرجات المئوية

AOCDEW - شروط مرصد أسكوم & # 160- نقطة الندى بالدرجات المئوية

ظروف مرصد AOCRAIN - ASCOM & # 160– معدل هطول الأمطار بالملليمتر / الساعة

ظروف مرصد AOCHUM - ASCOM - الرطوبة في المائة

ظروف مرصد AOCWIND - ASCOM & # 160– سرعة الرياح م / ث

شروط مرصد AOCWINDD - ASCOM & # 160- اتجاه الرياح بالدرجات (0..360)

ظروف مرصد AOCWINDG - ASCOM & # 160 - هبوب الرياح في متر / ثانية

ظروف مرصد AOCBAROM - ASCOM & # 160 - الضغط الجوي في hPa

AOCCLOUD - شروط مرصد ASCOM & # 160 - التغطية السحابية بالنسبة المئوية

AOCSKYBR – ASCOM Observatory Conditions – Sky brightness in Lux

AOCSKYQU – ASCOM Observatory Conditions – Sky quality (magnitudes per square arcsecond)

AOCSKYT – ASCOM Observatory Conditions – Sky temperature in degrees C

AOCFWHM – ASCOM Observatory Conditions – Seeing FWHM in arc seconds

APTDIA – diameter of the telescope in millimeters.

APTAREA – aperture area of the telescope in square millimeters. This value includes the effect of the central obstruction.

BAYERPAT – if present the image has a valid Bayer color pattern.

BOLTAMBT – Boltwood Cloud Sensor ambient temperature in degrees C.

BOLTCLOU – Boltwood Cloud Sensor cloud condition.

BOLTDAY – Boltwood Cloud Sensor daylight level.

BOLTDEW – Boltwood Cloud Sensor dewpoint in degrees C.

BOLTHUM – Boltwood Cloud Sensor humidity in percent.

BOLTRAIN – Boltwood Cloud Sensor rain condition.

BOLTSKYT – Boltwood Cloud Sensor sky minus ambient temperature in degrees C.

BOLTWIND – Boltwood Cloud Sensor wind speed in km/h.

CALSTAT – indicates calibration state of the image B indicates bias corrected, D indicates dark corrected, F indicates flat corrected.

CENTAZ – nominal Azimuth of center of image in degrees.

CENTALT – nominal Altitude of center of image in degress.

CBLACK – indicates the black point used when displaying the image (screen stretch).

CSTRETCH – initial display screen stretch mode.

CCD-TEMP – actual measured sensor temperature at the start of exposure in degrees C. Absent if temperature is not available.

COLORTYP – type of color sensor Bayer array or zero for monochrome.

CWHITE – indicates the white point used when displaying the image (screen stretch).

DATAMAX – pixel values above this level are considered saturated.

DAVRAD – Davis Instruments Weather Station solar radiation in W/m^2

DAVRAIN – Davis Instruments Weather Station accumulated rainfall in mm/day

DAVAMBT – Davis Instruments Weather Station ambient temperature in deg C

DAVDEW – Davis Instruments Weather Station dewpoint in deg C

DAVHUM – Davis Instruments Weather Station humidity in percent

DAVWIND – Davis Instruments Weather Station wind speed in km/h

DAVWINDD – Davis Instruments Weather Station wind direction in deg

DAVBAROM – Davis Instruments Weather Station barometric pressure in hPa

EXPTIME – duration of exposure in seconds.

DARKTIME – dark current integration time, if recorded. May be longer than exposure time.

EGAIN – electronic gain in photoelectrons per ADU.

FILTER – name of selected filter, if filter wheel is connected.

FLIPSTAT – status of pier flip for German Equatorial mounts.

FOCALLEN – focal length of the telescope in millimeters.

FOCUSPOS – Focuser position in steps, if focuser is connected.

FOCUSSSZ – Focuser step size in microns, if available.

FOCUSTEM – Focuser temperature readout in degrees C, if available.

IMAGETYP – type of image: Light Frame, Bias Frame, Dark Frame, Flat Frame, or Tricolor Image.

INPUTFMT – format of file from which image was read.

ISOSPEED – ISO camera setting, if camera uses ISO speeds.

JD or JD_GEO – records the geocentric Julian Day of the start of exposure.

JD-HELIO or JD_HELIO – records the Heliocentric Julian Date at the exposure midpoint.

MIDPOINT – UT of midpoint of exposure.

NOTES – user-entered information free-form notes.

OBJECT – name or designation of object being imaged.

OBJCTAZ – nominal azimuth of center of image

OBJCTDEC – Declination of object being imaged, string format DD MM SS, if available. Note: this is an approximate field center value only.

OBJCTHA – nominal hour angle of center of image

OBJCTRA – Right Ascension of object being imaged, string format HH MM SS, if available. Note: this is an approximate field center value only.

PEDESTAL – add this value to each pixel value to get a zero-based ADU. Calibration in MaxIm DL sets this to 100.

PIERSIDE – indicates side-of-pier status when connected to a German Equatorial mount.

READOUTM – records the selected Readout Mode (if any) for the camera.

ROTATANG – Rotator angle in degrees, if focal plane rotator is connected.

ROWORDER – Images taken by MaxIm DL are always TOP-DOWN.  

SBSTDVER – string indicating the version of the SBIG FITS extensions supported.

SET-TEMP – CCD temperature setpoint in degrees C. Absent if setpoint was not entered.

SITELAT – latitude of the imaging site in degrees, if available. Uses the same format as OBJECTDEC.

SITELONG – longitude of the imaging site in degrees, if available. Uses the same format as OBJECTDEC.

SNAPSHOT – number of images combined.

SWCREATE – string indicating the software used to create the file will be ”MaxIm DL Version x.xx”, where x.xx is the current version number.

SWMODIFY – string indicating the software that modified the file. May be multiple copies.

TILEXY – indicates tile position within a mosaic.

TRAKTIME – exposure time of the autoguider used during imaging.

XBAYROFF – X offset of Bayer array on imaging sensor.

YBAYROFF – Y offset of Bayer array on imaging sensor.

XBINNING – binning factor used on X axis

XORGSUBF – subframe origin on X axis

XPIXSZ – physical X dimension of the sensor's pixels in microns (present only if the information is provided by the camera driver). Includes binning.

YBINNING – binning factor used on Y axis

YORGSUBF – subframe origin on Y axis

YPIXSZ – physical Y dimension of the sensor's pixels in microns (present only if the information is provided by the camera driver). Includes binning.

In addition, when PinPoint Astrometry is used, World Coordinate System (WCS) information will be added to the FITS header.

Some of the above parameters can be viewed in the File Open dialog. All are visible in the FITS Header Window.


Thread: How to cut out a region from an SDSS FITS image file?

FITS files of Fields contain the science data in SDSS images, one per band (or filter).

Using PyFITS, you can select the image data for any region within a 1489 x 2048 pixel FITS image. Five times (one for each band/filter).

But how to select the (five) regions so they are each centered on the same (RA, Dec) point (within half a pixel)? .

I'm guessing you can work out the position, in pixel coordinates, of a particular (RA, Dec) point using the Header fields "RA" ("1st row - Right ascension of telescope boresigh"), "Dec" ("1st row - Declination of telescope boresight (d"), "SPA" ("1st row - Camera col position angle wrt north ("), etc . but how? . .

I just posted this in the Galaxy Zoo forum. I'll be interested to see how the responses differ .

The coordinate information in FITS files is encoded in sets of keywords describing the world coordinate system (WCS), which may take se several forms. For direct images, these typically specify a reference point in RA and dec, the pixel location of this reference point, the pixel scale (which may in general differ for each axis), and the orientation of the pixel grid on the sky. In SDSS files, they look like this:


CTYPE1 = 'RA---TAN'
CTYPE2 = 'DEC--TAN'
CUNIT1 = 'deg '
CUNIT2 = 'deg '
CRPIX1 = 1.02450000000000E+03 / Column Pixel Coordinate of Ref. بكسل
CRPIX2 = 7.44500000000000E+02 / Row Pixel Coordinate of Ref. بكسل
CRVAL1 = 2.47486643250000E+02 / RA at Reference Pixel
CRVAL2 = 2.43901742700000E+01 / DEC at Reference Pixel
CD1_1 = 6.89565533678154E-05 / RA degrees per column pixel
CD1_2 = 8.57010978890200E-05 / RA degrees per row pixel
CD2_1 = 8.56370751953130E-05 / DEC degrees per column pixel
CD2_2 = -6.9012795698925E-05 / DEC degrees per row pixel

As is common, the coordinates assume the so-called tangent-plane projection from the celestial sphere to the flat detector (which is fine for small fields of view unless the optics have very strong radial distortion). CDn_m specifies the change in celestial coordinate n due to a pixel change in image coordinate m CD1_2 and CD2_1 are 0 when the pixel grid is aligned with the coordinate axes.

The traditional aproach in astrometry is to define standard coordinates xi, eta:
Xi = cd1_1*(x-crpix1) + cd1_2*(y-crpix2)
Eta = cd2_1*(x-crpix1) + cd2_2*(y-crpix2)

and if there is no further distortion (or it has been corrected by resampling the image) transform those coordinates into angular ones:
cot δ sin (RA - CRVAL1) = (ξ) / (sin CRVAL2 + η cos CRVAL2) and cot δ cos (RA - CRVAL1) = (cot CRVAL2 - η sin CRVAL2) / (sin CRVAL2 + η cos CRVAL2)
(I see the Greek letters in those if they're lost, these are in the "Narrow-field astrometry" section here). I think I got the translation from RA, dec symbols to CRVALx properly. For very small coordinate differences, one can often simply pretend everything is a linear transformation from a notional rectangular RA/dec grid to pixel space.

The defining document for representing celestial coordinates in FITS is by Greisen and Calabretta most of it deals with assorted all-sky projections rather than narrow-field dorect imaging, where the transformation is rather simpler.

ETA: You may have seen that the various SDSS filter images on a given field have registration differences of several pixels, which is why this exercise needs doing to get matched subimages. For such small offsets, you can use the simplest Cartesian approximation to much better than a pixel accuracy. The important information is, in this case with the same pixel scale and orientation for each image, carried in the CRVAL and CRPIX sets of keywords.


An introduction to coordinate systems used in Astronomy

In geometry, a coordinate system is a system which uses one or more numbers, or إحداثيات, to uniquely determine the position of the points or other geometric elements on a manifold such as Euclidean space.

The following text briefly explains the coordinate systems being used in astronomy, some of which are listed below:

RA / DEC

RA (right ascension) and DEC (declination) are the longitudes and latitudes of the sky. RA corresponds to east / west direction, similar to longitude, while DEC measures north / south directions, like latitude.

WCS

World Coordinate System ( WCS ) is a set of transformations that map pixel locations in an image to their real-world units, such as their position on the sky sphere. These transformations can work both forward (from pixel to sky) and backward (from sky to pixel).

FITS WCS

The FITS “World Coordinate System” ( WCS ) standard defines some conventions and keywords to associate coordinates with each pixel of an image.

FITS

Flexible Image Transport System ( FITS ) is a digital file format useful for storage, transmission and processing of scientific and other images. It is the defacto standard used by many sky tessellation softwares — in this case, HEALPix.

HEALPix header files can contain the following three letters, each depicting the coordinate system being used:

  • ج:Celestial = ICRS = RA / DEC (equatorial) = FK5 J2000 (default)
  • جي:Galactic
  • ه:Ecliptic

International Celestial Reference System ( ICRS )

ICRS is the current standard celestial reference system adopted by the International Astronomical Union ( IAU ). انها origin is at the barycenter of the Solar System, with axes that are intended to be “fixed” with respect to space - this is referred to as International Celestial Reference Frame ( ICRF ). ICRS coordinates are approximately the same as equatorial coordinates.

Fifth Fundamental Catalogue ( FK5 )

FK5 is part of the “Catalogue of Fundamental Stars” which provides a series of six astrometric catalogues of high precision positional data for a small selection of stars to define a celestial reference frame. J2000 refers to the instant of 12 PM (midday) on 1st January, 2000. FK5 was published in 1991 and added 3,117 new stars.

Galactic coordinate system

The galactic coordinate system is a celestial coordinate system in spherical coordinates, with its origin at the Sun, the primary direction aligned with the approximate center of the Milky Way galaxy, and the fundamental plane parallel to an approximation of the galactic plane but offset to its north. GCS has its own Galactic longitude and Galactic latitude.

Ecliptic coordinate system

A celestial coordinate system commonly used for representing the positions and orbits of Solar System objects. The system’s origin can either be the center of the Sun or the center of the Earth, its primary direction is towards the vernal (northbound) equinox, and it follows a right-handed convention.


New Coordinate Transformation Tasks

Three new coordinate transformation tasks imcctran , skyctran , and wcsctran have been installed in the new IRAF V2.11 imcoords package. These tasks use the image header coordinate system to perform one or more of the following functions: 1) precess the image coordinate system ( imcctran ), 2) convert the image coordinate system from one celestial coordinate system to another ( imcctran ), 3) transform from pixel to world coordinates and vice versa ( skyctran and wcsctran ), or 4) locate objects detected in one image in another ( skyctran ). Skyctran can also be used to: 1) precess equatorial coordinates, or 2) convert coordinates from one celestial coordinate system to another. IRAF V2.11 supports equatorial, ecliptic, galactic, and supergalactic image celestial coordinate systems earlier versions support equatorial image coordinate systems only.

The following examples demonstrate the basic functionality of the new tasks. To reproduce these examples on a local system copy the image dev$wpix to a local directory and edit in the missing EQUINOX keyword as shown below.

Precess or Transform the Image Header Coordinate System

Imcctran converts the image coordinate system from one celestial coordinate system to another. The conversion shifts and rotates the image coordinate system but leaves the pixel coordinates of the reference point and the coordinate projection unchanged. Precession is treated as a special coordinate conversion.

The following example converts the coordinate system of wpix from equatorial B1950.0 to equatorial J2000.0 and then from equatorial J2000.0 to galactic. The initial image copy avoids overwriting the header of wpix which is used in later examples.

Transform from Pixel to World Coordinates andVice Versa

The tasks wcsctran and skyctran convert from pixel to world coordinates and vice versa using the image coordinate system. Wcsctran works on images of any dimension with any valid coordinate system. The task skyctran works only on two-dimensional images with valid celestial coordinate systems.

The first two examples show how to use wcsctran to convert from pixel coordinates to world (in this case equatorial) coordinates and back again. The format and units specifications ensure that the world coordinates are written and read in hours and degrees.

cl> wcsctran pix.coo eq.coo wpix logical world
>>> formats="%12.3H %12.2h"
cl> wcsctran eq.coo npix.coo wpix world
>>> logical units="h n" formats="%8.3f %8.3f"
cl> type eq.coo
cl> type npix.coo

cl> wcsctran STDIN STDOUT wpix logical
>>> world formats="%12.3H %12.2h"
. type in x and y pixel coordinates and hit return
. type <EOF> to quit
cl> wcsctran STDIN STDOUT wpix world
>>> logical units="h n" formats="%8.3f %8.3f"
. type in world coordinates in hours and degrees and hit return
. type <EOF> to quit

The next two examples show how to use skyctran to convert from pixel to galactic coordinates without modifying the image coordinate system which in this case is an equatorial coordinate system.

    cl> skyctran pix.coo gal.coo wpix galactic
    cl> type gal.coo

The final example shows how skyctran can be used to locate objects for which only galactic coordinates are available in an image with an equatorial coordinate system. The tvmark task is used to mark the detected objects on the image display.

Locate Objects Detected in One Image in Another

Skyctran can locate objects detected in one image in another image which has a different celestial coordinate system, e.g., equatorial and galactic as shown in the following example. Starfind is an imcoords task which automatically locates stellar objects in images.

    cl> imdelete wpix.tmp # if wpix.tmp already exists
    cl> imcopy wpix wpix.tmp
    cl> imcctran wpix.tmp galactic
    cl> display wpix.tmp 1 fi+
    cl> starfind wpix pix.wpix 1.25 100
    cl> skyctran pix.wpix pix.wpix.tmp wpix
    >>> wpix.tmp transform+
    cl> tvmark 1 pix.wpix.tmp col=204

Precess or Convert Celestial Coordinates

The skyctran task is a general celestial coordinate conversion tool as well as an image coordinate conversion tool as shown in the following two examples.


Transform pixel coordinates (in FITS file) to equatorial - Astronomy

I decided to create a library from the codes that I use for myself. I hope it will be useful for you too.

And here are our dependencies for this library:

astrolib | |----> astronomy.py | |----> catalog.py | |----> io.py | |----> visuals.py | |----> photometry.py For detailed information and help give help(module_name) command in the command line.

Here I created help documentation for the commands I use the frequently. You can work the same way for others. You can find detailed help in the code.

If you just want to get physical coordinates of sources

If you want to plot detected objects on the FITS file

threshnpixtnpixxminxmaxyminymaxxذx2y2xyerrx2erry2errxyأبthetacxxcyycxycfluxfluxcpeakpeakxcpeakycpeakxpeakypeakflag
float64int64int64int64int64int64int64float64float64float64float64float64float64float64float64float64float64float64float64float64float64float64float64float64float64int64int64int64int64int64
19.1716499329514482770798696721784.445069846708.5999344747.130377780347.14998870343-0.2123581408727.09243459217e-066.17207043114e-06-1.98509663699e-072.711598634722.63203310966-0.8084688782690.1403691917660.1399841904640.00833806581795723850.0724231.812517705.060546918033.35351567857097847090
19.1716499329406375491512318340501.969623282328.8562518266.851296200456.80527345789-0.2320989369691.135928646e-051.08330602851e-051.04181566067e-062.657352447512.56808257103-0.7359873056410.1461266130210.147114843130.00996751524508443818.09375444089.7513129.044921914977.50585945023295023280
19.1716499329288275254274438457265.640606005447.070895386.239602344276.07421432435-0.2796716715783.41360150936e-053.11016713517e-05-2.11652827665e-062.539399385452.4218313694-0.6416526436810.1605980545280.1649708002810.0147886537015185113.59375185279.531256102.232910167160.669921882664472664471
19.1716499329277256424442371389433.446658084379.954207576.353162226586.31800704473-0.2998433333794.85704037652e-054.50648124042e-05-1.72250253383e-062.576032400132.45666980743-0.7561203837390.1577552855010.1586330682040.0149736674502146899.203125147146.468754537.379882815280.504882814343804343800
19.1716499329198181424440899914432.153302841906.3693695745.903030420035.47633612435-0.1847532566650.0001646224050430.000146496989951-3.67874889536e-062.443748712542.32539439201-0.3568477332590.1695835739370.1827968508010.011442365124857640.042968857828.26953121876.255249022155.505371094329064329060
19.1716499329194176979994725740986.64280869732.1731027425.82860409715.77761345379-0.1236827776740.0001498288619250.000146736762267-1.92272602728e-062.435034275052.38260889053-0.6837545037270.1716456413270.1731605082750.0073489197529957640.687557787.36328121833.149780272154.96215829877329877320
19.1716499329191179421436494509428.481743905502.052570895.730493288785.83927472126-0.082716463680.0001858823948180.0001969069065816.81651164286e-062.425671100622.38450980186-1.076232552530.1745407283310.1712891757490.0049449265934549968.70312550177.3906251583.442504881840.379272464295024295020
19.1716499329186172748763528543755.619674466535.119803485.849523159595.86444076657-0.3025529664080.0001905489472940.000193111256106-1.3837525619e-052.481859445572.35676407814-0.7977221012120.1714114993810.1709754765030.017686616629448653.160156248840.031251514.218017581737.780151377565357565350
19.1716499329217195968983837859975.565085076850.9583006495.450414391988.61544221634-0.3695535180930.0002447681357190.000640536321842-2.73328688175e-052.942451238632.32547569275-1.456090331080.1840074509380.1164092123510.01578574813941191.07812541409.03906251270.405517581463.905395519768519768510
19.1716499329146133304317891904311.27460779897.8792616114.954607708995.00588141734-0.1893835352640.0002520843412260.000276073972364-1.24180158269e-052.274061441422.18840885162-0.8526742458340.2021246105430.2000543177130.015293640084632531.558593832735.93554691089.452880861271.765380863118983118980

Query and match detected objects with the GAIA DR1

مصدرRA_ICRSDE_ICRSe_RA_ICRSe_DE_ICRS__Gmag_pmRApmDEe_pmRAe_pmDEEpochPlx
degdegmasmasmagmas / yrmas / yrmas / yrmas / yryrmas
int64float64float64float32float32float32float64float64float32float32float32float32
6849534850623923456306.7733664604-23.51496447590.2620.22117.121--------2015.0--
6849534850623925120306.7712467689-23.51257927340.8700.74919.499--------2015.0--
6849534988062874624306.7217934709-23.51876625260.6900.58919.120--------2015.0--
6849534988062875264306.7263979238-23.51816992390.3250.26417.535--------2015.0--
6849535056782354560306.7183059223-23.51266948630.1960.13317.308--------2015.0--

If you just want to query the target object, reduce the radius like

مصدرRA_ICRSDE_ICRSe_RA_ICRSe_DE_ICRS__Gmag_pmRApmDEe_pmRAe_pmDEEpochPlx
degdegmasmasmagmas / yrmas / yrmas / yrmas / yryrmas
int64float64float64float32float32float32float64float64float32float32float32float32
6849629099386034304306.7717090460-23.47027877420.2160.1269.6038.0594.8160.9100.6192015.03.17

idxذradece_rae_decg_mean_magpmrapmdece_pmrae_pmdecepochplxfluxأبthetara_calcdec_calcra_diffdec_diff
float64float64float64float64float64float64float64float64float64float64float64float64float64float64float64float64float64float64float64float64float64float64
6.84962909939e+18784.445069846708.599934474306.771709046-23.47027877420.2160000056030.1260000020279.602999687198.0594.8160.9100000262260.6190000176432015.03.17000007629724231.81252.711598634722.63203310966-0.808468878269306.771497336-23.4704661458762.154325162674.53790621
6.84964298072e+18501.969623282328.856251826306.825119424-23.39661054310.1550000011920.096000000834510.1719999313nannannannan2015.0nan444089.752.657352447512.56808257103-0.735987305641306.824877954-23.396690687869.293691494288.518125998
6.84963092045e+18265.640606005447.07089538306.874457724-23.41640263490.349999994040.17100000381511.090999603319.918-7.0571.930999994281.282999992372015.03.65000009537185279.531252.539399385452.4218313694-0.641652643681306.874213526-23.4164391737879.11293615131.539743526
6.84963122969e+18433.446658084379.95420757306.839605839-23.40549171280.1959999948740.14000000059611.3380002975-1.037-1.8961.432999968530.9520000219352015.00.579999983311147146.468752.576032400132.45666980743-0.756120383739306.839362608-23.405562815875.629978918255.967955209
6.84962782808e+18432.153302841906.369369574306.845545923-23.50381569540.2319999933240.16300000250312.3310003281-6.099-11.5841.672000050541.09099996092015.01.1499999761657828.26953122.443748712542.32539439201-0.356847733259306.845308447-23.5039523172854.913985108491.838601999
6.84954694525e+18986.64280869732.173102742306.730808533-23.47662637960.1529999971390.11999999731812.3090000153nannannannan2015.0nan57787.36328122.435034275052.38260889053-0.683754503727306.73060855-23.4768569137719.935376696829.922643453
6.8496409535e+18755.619674466535.11980348306.775712653-23.43759406540.0970000028610.075000002980212.5030002594nannannannan2015.0nan48840.031252.481859445572.35676407814-0.797722101212306.775505562-23.4377531371745.527969957572.658003705
6.84953515986e+18975.565085076850.958300649306.734327716-23.49875473090.1049999967220.08200000226512.7259998322nannannannan2015.0nan41409.03906252.942451238632.32547569275-1.45609033108306.734126873-23.4989537639723.032314954716.51890772
6.84962820603e+18311.27460779897.879261611306.870049861-23.50105919440.1099999994040.072999998927112.875nannannannan2015.0nan32735.93554692.274061441422.18840885162-0.852674245834306.869815661-23.5011624192843.118385296371.609451281

Astrometry with IRAF's ccmap

First of all we need a FITS file that has been resolved with astrometry.net (For example: 108hecuba-001_R_affineremap.fits). Then our code will map the resources found here to the GAIA catalog and again perform astrometry with these coordinates via ccmap.

Ra/Dec or Long/Lat fit rmsRa/Dec or Long/Lat wcs rmsReference point (RA, DEC)Reference point (X, Y)X and Y scaleX and Y axis rotation
str16str16str16str16str28str18
0.0380.03820:27:17.965519.0310.673177.025
0.04850.0485-23:27:09.36627.6770.673176.983
(arcsec arcsec)(arcsec arcsec)(hours degrees)(pixels pixels)(arcsec/pixel arcsec/pixel)(degrees degrees)

You can see the results coordinates with

Ra/Dec or Long/Lat fit rmsRa/Dec or Long/Lat wcs rmsReference point (RA, DEC)Reference point (X, Y)X and Y scaleX and Y axis rotation
str16str16str16str16str28str18
0.0380.03820:27:17.965519.0310.673177.025
0.04850.0485-23:27:09.36627.6770.673176.983
(arcsec arcsec)(arcsec arcsec)(hours degrees)(pixels pixels)(arcsec/pixel arcsec/pixel)(degrees degrees)

But this solution does not include parallax and proper motion correction. To include these corrections

Ra/Dec or Long/Lat fit rmsRa/Dec or Long/Lat wcs rmsReference point (RA, DEC)Reference point (X, Y)X and Y scaleX and Y axis rotation
str16str16str16str16str28str18
0.03770.037720:27:17.965519.0310.673177.025
0.04860.0486-23:27:09.36627.6770.673176.983
(arcsec arcsec)(arcsec arcsec)(hours degrees)(pixels pixels)(arcsec/pixel arcsec/pixel)(degrees degrees)

Coordinates α and δ of a star in a fixed reference system change with time proportional to its proper motion μα and μδ. Let α0 and δ0 be its position at some time origin its values at time t are

Stellar Parallax Correction

Deltaα = π (x sin α − y cos α) / cos δ

Deltaδ = π [(x cos α + y sin α) sin δ − z cos δ]

where the coordinates x, y, z of the Earth are expressed in astronomical units and the corrections to the position are in arcseconds.

Ref: Kovalevsky, J., & Seidelmann, P. (2004). Fundamentals of Astrometry. كامبريدج: مطبعة جامعة كامبريدج. doi:10.1017/CBO9781139106832


شاهد الفيديو: تحويل الإحداثيات الجغرافية Latitude, Longitude إلى الإحداثيات الديكارتية Easting,Northing (يوليو 2022).


تعليقات:

  1. JoJotilar

    إكسسوارات المسرح تنجح

  2. Elkanah

    اعذروني على ما يجب أن أتدخل فيه ... وضع مشابه. دعوة المنتدى. اكتب هنا أو في PM.

  3. Darn

    الموضوع مثير للاهتمام ، وسوف أشارك في المناقشة. أعلم أنه يمكننا معًا الوصول إلى الإجابة الصحيحة.

  4. Melville

    نعم هذا كل الخيال



اكتب رسالة