Artyom

והנה אתר שעושה את הכל: https://searchlight.semrock.com/ מה שלא ברור לי למה NBP עובד יותר טוב מ־Narrow Band UHC כי שניהם מכסים את הפיקים החשובים משמאל לימין H-Beta באדום, OIII בירוק ו־H Alpha בכחול.

סידרתי את הגרפים של 4 פילטרים בסקאלה אחת. זה מסתדבר מסביר המון הפילטרים NBP/Astronomik UHC וגם OIII של svbony נותנים רוחב פס מאוד דומה מסביב ל־O3. לעומתו SVBony/Optolong UHC מעבירים טווח הרבה יותר רחב למעשה חוץ מ UHC זולים יותר הפילטור ב־O3 מאוד דומה הטווח אדומים UHC זול ויקר מאוד דומים פה בא ההבדל עם NBP יש טווח סינון שונה באדום וגם NBP מעביר קצת גם בכחול. נראה לי צריך לעשות הפרדה מאוד ברורה בין UHC עם טווח רחב או צר. במילים אחרות UHC בטווח צר יכול להיות משמעותית יותר טוב. לא יודע אם המסקנה נכונה אבל זה פתאום עושה לי סדר למעשה ההבדל בין פילטרים של UHC יכול להיות בין תוצאה מצוינת לסבירה והבדל בין התייחסות (אם הוא חסר תועלת או מאוד שימושי)

אם כבר המגמה אחרת. דווקא כל הפילטרים זמינים גם ב־1.25 וגם זולים יותר. הבעיה שרוב הטלסקופים שזה נבדק עליהם הם 6" ומעלה ורק אני עם 4" 😉

זהו נראה ל־Astronomik UHC הוא באמת narrow band מצד שני עם 102 מ"מ אולי זה באמת יחסום יותר מידי אור. טוב הרבה נקודות למחשבה. או אולי פשוט לחכות לתצפית ולשאול פילטרים ולראות מה עובד אצלי ומה לא על כמה מטרות. כי התמונה המדהימה שראיתי הייתה עם O3 של סביבוני על Veil עם 8" של @שחר רייש. (מצד שני, הבדיקה בשטח זו תמיד הדרך הטובה ביותר ללא קשר לסוג של ציוד)

איזה פלטר UHC יש לך (איזו חברה דגם) כי יש מגוון רחב של פילטרים שונים שגם ההעברה לפי אורך הגל משתנה מאחד לשני מעניין לקחתי FLO ולקחתי משם מחירים וגרפים איפה שאפשר אז הנה 4 דגמים ה־UHC של Optolong (זהה לפי הגרף ל־UHC של SVBony) רוחב 48nm טווח מכ־460 עד כ־510 (מדידה מהגרף לפי FWHM) טווח "אדום" מ־630 מחיר 38£ ה־ES UHC רוחב 50 טווח מ־458 ל־508 לפי מפרט טווח אדום לא צוין אין גרף מחיר 55£ ה־Astronomik UHC-E רוחב כ־55nm (לפי FWHM מהגרף) טווח שונה מכ־470 ועד כ־530nm טווח אדום מ־650 מחיר 58£ ה־Astronomic UHC רוחב כ־26nm טווח מתחיל בכ־480 נגמר בכ־505 (לפי FWHM מהגרף) טווח אדום מ־640 המחיר 89£ אז האם אלה עם פס רחב שווים משהו או שווים בעיקר במפתחים קטנים?

המשך הדיון שהתפתח בקבוצת ווטאפ פרטית. מהם הפילטרים המומלצים לצפיה בערפיליות ובאילו תנאים. ההמלצות הנפוצות בחוכחת ההמונים הן UHC ו־OIII גם פילטר יחסית חדש DGM Optics NPB. UHC בד"כ מכיל שני טווחים מעבירים את אור עם הדגשה על מספר איזורים חשובים (מ־first light optics עבור Astronomik ומאתר של svbony עבור של UHC של סביבוני וזהה לחלוטין של Optolong) מה שאגב מראה שיש גם הבדל בין פילרים לדוגמה Astronomik חותך סביבות 470nm לעומת סביבוני ו־optolong סביבות 450nm פילטר OIII מעביר טווח מאוד קצר בסביבות 500nm עבור שני פסי פליטה של חמצן נלקח מ־FLO ואתר של סביבוני ההעברה של פילטר DGM Optics NPB הרבה יותר מורכבת (אגב לא מצאתי את הגרף באתר שלהם אלא רק פה https://okularum.eu/products/dgm-npb-filter-2 ) מחירים עלות פילטר UHC נעה בין 38 ל־89 GBP ב־FLO כ־$20 של סביבוני עבור 1.25 עלות פילטר OIII אפשר למצוא מכ־$50 של סביבוני וב־FLO בין 43GBP לכ־90 GBP לגודל 1.25 עלות NPB ב־2" של כ־170$ ובאופן תיאורתי 75$ עבור 1.25 אבל אני אישית לא מצאתי משהו במלאי ‎אז איזה פילטר מומלץ? לפי סקירה של @Ilan Shapira שפורסמה כאן ה־NPB גרם לכך שהפילטר OIII נשאר במגירה. גם לפי הדיון שהתפתח בקבוצה פילטר UHC בעצם חסר כל תועלת. לעומת זאת לפי הסקירה הזו: https://www.prairieastronomyclub.org/filter-performance-comparisons-for-some-common-nebulae/ הפילטר UHC הוא השימושי ביותר שעובד על הכי הרבה ערפיליות. גם OIII ועושה עבודה נהדרת אבל על קצת פחות מטרות. בסקירה הזו לא הייתה השוואה ל־NPB (כי הוא לא היה קיים בשוק) חשוב לציין שהסקירה הזו התייחסה לתנאים של שמיים חשוכים. בדיון זה נדבר על אסטרונומיה תצפיתית - לא צילום דעתכם? ניסיונכם? באיזה תנאים? על אילו מטרות? איזה פילטרי UHC טובים ואיזה לא?

לאחר שקיבלתי דיווחים שהאפליקציה רצה גם 1920x1080 על usb camera ועל שלי בדקתי את האפליקציה על טאבלט של אישתי... עובדת גם ברזולוציה מלאה וגם עושה stacking הלי בעיות ביצועים. כנראה הבעיה היא בחומרה של הטלפון שלי 😞

עדכון לגבי 4 מ"מ לאחר שימוש ב־60/400 היא צריכה יחס אורך מוקד הרבה יותר גבוה. כשאני עובד עם 80/660 ב־F/8.25 זה בסדר. אבל ב־F/6.5 התמונה מטושטשת מידי. בקיצור - רק לטלסקופים איטיים

צילום מסך אפליקציה מעודכת מצלמה לצילום תמונה בודדת. התווסף פרמטר gamma החלפתי טקסט באייקונים איפה שאפשר

שוחררה גרסה 0.0.2 https://github.com/artyom-beilis/android_live_stacker/releases/tag/0.0.2

טוב כמה הצלחות מהירות ומעודדות - הכל צולם עם AstroMaster 102/660 מוד ידני עם darks עם זמן חשיפה של 0.5s עם gamma 2.2 ו־gain=0, בלי עקיבה בורטל 8 רחובות, רכישה ב־800x600. מקוריים: M13 כ־120 פריימים - דקה של צילום M92 גם כ־120 פריימים טבעת - רק כ־30 פריימים (15 שניות) - ואז איבדתי אותה (צריך עקיבה 🙂 ) אותם רק אחרי קצת stretch ב־offline

כדי לא לשרוף כוכבים. אבל בגדול אני משתמש ב־gamma במקום - היא מגבירה סיגנל חלש אבל לא שורפת כוכבים בהירים. הפלט של מצלמה רק 8 ביט - אז עדיף להשתמש ב־gamma כדי לדחוס יותר טוב את הטווח הדינאמי של הסיגנל ב־8 ביט

טוב שתי תוצאות - קצת התקדמות לגונה - ~400 פריימים של 0.25s סה"כ כ־100 שניות, gamma 2.1 ו־gain 0 צבור M92 כ־40 פריימים של 0.5s סה"כ כ־20 שניות, gamma 2.1 ו־gain 0 מצלמה כמובן sv105 טלסקופ SVBony 60/400 עקיבה ידנית (השהייה והזזה מידי פעם) הכל עם darks נאסף עם live stacking רק בוצע stretch קל לאחר מכן. כיוון שהתגלתה רגישות יתר במעקב אחרי רגיסטרציה בין הפריימים, אני לא אשחרר גרסה היום. אנסה לתקן בימים הקרובים.

כמה עדכונים הוספתי שליטה ב־gamma במצלמה וגם הוספתי gamma correction ל־stacking - בצורה זו אני עושה gamma הפוכה כדי לעשות stacking בתחום לינארי ואז מחזיר gamma סטנדרטית. נראה שמשפר משמעותית ראות העפתי את הכותרת המיותרת כדי לעבוד במסך מלא. הוספתי מחיקה של darks ברגע שמשנים הגדרות (צריך באופן כללי לבנות presets ולשייך להם darks) הוספתי זריקת frameים במידה ואני לא מספיק לבצע עיבוד (בגדול לא אמורה להיות בעיה בזמני חשיפה ארוכים) אשתדל לשחרר גרסה חדשה היום. בנוסף הבנתי שאני צריך כבר מאריך ל־usb כדי לא לשים ספק מתח על החצובה.

עדכונים מצפית בורטל 5 מצפה משואה מה ראיתי ב־60 ראיתי לגונה לא רע בכלל ראיתי M13 אם כי פחות מובהק ממה שראיתי ב־102 מ"מ בקיצור טלסקופ נחמד לחלוטין - אבדוק במדבר ואז ארשום את הסיכום הסופי חצובה - בפעם הראשונה השתמשתי בה בחוץ. כל עוד אין רוח חזק, החצובה (עם איזון על בסיס חבל+גומיות) עובדת לא רע בכלל. אפשר לעבוד לצפות וכד'. כשהתחיל הרוח אז בהגדלה x66 התמונה מתחילה כבר לרצד.

גיליתי שברחובות ממש ב־5 דקות נסיעה ממני יש חברה שמתמחה בהדפסת 3d ונותנת מחירים מאוד נוחים ושירות מהיר. אז הדפסתי מתאם לחצובה לראש כדורי והנה התוצאה של עבודה עם ראש כדורי בגדול - עובד מצויין - אפשר להביא את הטלסקופ כמעט ל־90 מעלות. אבל התגלו שתי בעיות הזזת ראש כדורי במיוחד עם מעמיסים גם טלפון עם AstroHopper על הטלסקופ דורשת קצת יותר כוח משיטות איזון פשוטה עם הגומיה. משמעות היא שעל החצובה הגמישה הזו ה־backlash הופך למשמעותי - כמעלה - מה שמאוד מקשה על ניווט. אם דיאגונל יותר קל וללא טלפון ה-backlash כ־0.5 מעלה בציר הגובה. אם החצובה הייתה רצינית זאת לא הייתה בעיה, אבל זה לא מתאים חצובה כל־כך קלילה. בנוסף בגלל העדר מתח הופסת משקולת על הוו הופכת לקריטית מבחינת רעידות אבל ראש כדורי כציר אסטונומי מאוד נחמד הבעיה השניה שזה מסובב את הפיינד לצד ימין. עם זה היה RDF זה לא היה נורא אני מניח למרות שעיין העינית הייתה מפריעה. עם פיינדר אופטי זה חובה להתסכל בעין שמשאל שהיא לא דומיננטית - מאוד לא נוח. בדקתי עוד שתי שיטות נוספות: הזזת טלסקופ קדימה על מסילה: הבעיה פה שהטווח שזה מאזן מאוד קטן צריך הלזיז את הטלסקופ קדימה אחורה כדי לקבל יציבות כלשהי שימוש בעובדה שאפשר לסובב ראש ב־90 מעלות ולחבר את הטלסקופ והצד ולהוריד את מרכז הכובד לסיבוב: זה מוסיף המון לחץ צידי מגביר backlash - ניסיתי קצת, לא היה מוצלח במיוחד. שורה תחתונה: השיטה הטובה ביותר היא הפשוט והזולה ביותר חבל וגומיה! לחצובה רצינית אני הייתי ממש מתלבט בין ראש כדורי לבין איזון מכני עם גומי של pan handle כי פשוט ראש כדורי הרבה יותר כבד מעלה את משקל החצובה ב־1/3

הוראות שימוש שימו לב כרגע אפליקציה תומכת אך ורק מצלמות USB חיצוניות שעובדות בפרוטוקול UVC. היא לא תומכת עדיין במצלמה מובנית וגם לא במצלמות. מצלמות שתומכות ב־UVC כמו שידוע לי הן של svbony דגמים sv105 ו־sv205. נראה לי יש גם מצלמה של Bresser שגם עובדת עם UVC. כמו כן מלצמות ווב גנריות שהוסבו לצורכי אסטרונומיה קודם כל צריך להתקין אפליקציה בעזרת APK: נכנסים לקישור: https://github.com/artyom-beilis/android_live_stacker/releases/tag/0.0.1 ומורידים לטלפון. עושים התקנה. שימוש קודם מחברים מצלמה -בד"כ צריך מתאם מ־USB זכר ל־USB Type-C זכר (אם יש בטלפון/טאבט חיבור type-c) לדוגמה: מצלמת sv105/sv205 יכולה לדרוש ספק כוח חיצוני נוסף אם הטלפון/טאבלט לא ייתנו את ההספק הנדרש - פה סוללה לטעינת טלפון עושה את העבודה. קודם מחברים ספק ואז מחברים מצלמה לטלפון מפעילים את האפליקציה אחרי(!) חיבור המצלמה - היא תבקש אישור לגישה להתקן USB חיצוני בוחרים את הרזולוציה של העבודה - ההמלצה תתחילו מרזולוציה נמוכה כמו 640x480. אני הצלחתי להפעיל את sv105 יציב עד 800x600 אבל זה יכול להיות תלוי מכשיר שלכם. לפעמים לוקח זמן למצלמה "להתחמם" אם בפעם הראשונה לא מצליחים לראות וידאו, יוצאים נכנסים שוב. תפעול הכיוונון ניתן לכבות את Auto-Exposure ו־Auto Balance ע"י ביטול checkbox של AE - שזה auto-exposure ניתן לקבוע זמן חשיפה את מצב ה־white balane וה־gain ע"י 3 סליידרים כפתור capture ישמור תמונה נוכחית ה־Save Checkbox אומר לכם לשמור את התמונות שמהם עושים stacking לעיבוד ב־offline - שימו לב זה מעמיס ומייצר המון תמונות כפתור stack יתחיל תהליך live-stacking: בלוח הבקרה תראו מספרים הבאים שהם מראים משמאל לימין מספר תמונות שנרכשו ונשלחו ל־stacking, מספר תמונות שעברו stacking ומספר תמונות שהצליחו (כרגע עוד לא שמתי הגנה על קפיצות רגיסטרציה אז הם יהיו זהים) שימו לב אם עובדים ב־framerate גובה ו־frame גדול ואולי גם מפעילים שמירת כל תמונות ה־stack אז יכול להיות שתהליך של stacking יהיה איטי יותר של רכישת התמונות ותראו הבדל משמעותי בין מספר ראשון ושני כדי לעצור תהליך stacking לוחצים על Pause אחרי העצירה אפשר להמשיך את התהיך ע"י לחיצה על כפתור continue או לעצור סופית ולשמור את התוצאה ע"י לחיצה על Finish שימו לב שמתחת לכפתורים תוצג תמונה מוקטנת את הפריים שמצלמים בזמן ה־stacking - כך תוכלו לראות במידה והעקיבה שלכם לא טובה או עושים צילום ללא עקיבה שהאובייקט מתרחק משדה הראיה וצריך להשהות את הצילום ולתקן לחיצה על continue תמשיך את התהליך התמונות נשמרות ב־DCIM תחת ספרייה LiveStacker בהצלמה

טוב שוחררה גרסת א: https://github.com/artyom-beilis/android_live_stacker/releases/tag/0.0.1 למזוכיסטים שרוצים מנסות

עדכונים: תיקנתי חלק מבעיות יציבות הקשורות לפתיחת מצלמה מחדש (מסתבר שאם מסובבים את המסך נפתח activty חדש... וואו) אבל עדיין התמונה היציבה מתקבלת רק ב־800x600. אבל ברגע שניסית לצלם גיליתי שהאלגוריתם של stacking מקרטע ולא מבצע רגיסטרציה בין הפריימים כמו שצריך. צילמתי M92 ואחרי תיקון האלגוריתם הנה תוצאה על 35 פריימים של 0.5 שניות על טלסקופ של 102/660 - העיבוד בוצעה ב־Offline מסקנות: אני צריך לתקן פוקוס 🙂 ולנסות שוב בלילה - צריך להכין מסכה של בחטינוב. צריך לשמור נתוני ביניים ולא רק התמונה הסופית צריך להוסיך תמיכה ב-darks AstroMaster 102AZ, sv105, 35 frames, 0.5s, 100% gain, basic postprocessing (stretching), offline

עוד עדכונים: הצד הטוב: פתרתי בעיות ריצת ה־opencv שהיו לי. יש stacking באפליקציה ולפי בדיקות ראשוניות אין שום בעיות בזמני ריצה/חישובים - עשיתי live stacking גם בכ־30 FPS (אם כי על 800x600) - כך שכנראה עם זמני חשיפה אסטרונומיים של 100-500ms לא צפויה שום בעיה גם ב־full-hd (בהנחה שאצליח להפעיל אותו). הצד הפחות טוב יש בעיות יציבות במצלמה של svbony. אומנם מצלמת web כן עובדת ב־800x600, לעומתה sv105 שלי עובדת לסירוגין - כן, לא ולמה לא. כך שאתמול לא ממש הצלחתי לצלם אפילו שהיה לי ראות ברורה על m93. מה שכן ראיתי שמצלמת ה־web שלי באפליקציה אחרת עובדת אפילו ב־1920x1080 כך שיש עוד מקום לחפירות ושיפורים - אני אופטימי. אני צריך לעבוד על נושא יציבות החיבור וגודל הפריים.

העליתי את הקוד ל־github אפליקציה עצמה: https://github.com/artyom-beilis/android_live_stacker רישיון GPLv3 ספריית C ו־C++ שאמורה לעשות את העבודה - לתקשר עם מצלמת USB וגם לעשות stacking בהמשך: https://github.com/artyom-beilis/uvcctl רישיון MIT ספריית libuvc שהכנסתי לתוכה שינויים שיתמכו בשינוי גודל ה־buffer לטובת אנדרואיד: https://github.com/artyom-beilis/libuvc קוד האפליקציה כמובן בינתיים ברמת סקיצה ודורש המון refactoring לפי שאפשר יהיה להכריז עליו מבצעי

טוב עוד התקדמות. רצתי עם 800x600 על sv105 על טלסקופ 60/400. עשיתי stacking ב־offline עדיין לא מימשתי את ה־live stacking אבל הקוד כבר מוכן

האם זה עבור משקפת או עבור טלסקופ? כי ballhead לא ממש מתאים למשקפת. לטלסקופ קל כן, בהנחה שהוא מאוזן

על אובייקט כלשהו כדי שאפשר יהיה לראות תזוזה צריך סרטון פרוס לתמונות בודדות. יכול להיות כל פורמט jpeg/png וכד: זו השורה שטוענת: p = "%s%d.ppm" % (sys.argv[2], n) כאשר החלק הראשון מבציע על המסלול הבסיסי נגיד images/img והשני מספר - ו־n כמובן מספר בסדרה. אני עבדתי על תמונות 640x480 אבל כל דבר יעבוד - אולי צריך יהיה להגדיל חלון חיפוש בפונקציית crop פשוט (אני מניח שאתה עובד עם מצלמה אסטרונומית. אם פיקסל נגיד 3um אז ואורך המוקד 400 מ"מ אז גודל הפיסקל ברדיאן הוא 0.003/400 ובשניות הקשר חלק ב־pi ותכפל ב־180*60*60 כמובן אם אתה עובד עם פריים מוקטן פי 2 או 3 (כמוני 640 במקום 1920) אז תכפיל בהתאם וודאי, סה"כ תרגום מפריימים לשניות הוא חלוקה ב־30 - תסתכל בשורות בהם אני עושה PLOT אמרתי לא כתבתי סקרפיט גנרי אז תצטרך לעשות שיפצורים

הנה הקוד: קריאה python3 script.pt "title for graph" base_path_to_files כאשר הסקריפט מניח שהכל צולם ב־30 fps ב־10 שניות והתמונות ממוספרות מ־0 עד 299 (שורה 37). בשורה 64 יש סקלא שאומרת כמה שניות הקשת יש בכל פיקסל... אתה מוזמן לחפור ולשפר - מבחינתי הקוד ברשות הכלל (public domain) לא הכנתי אותו לשימוש כללי אז אם אתה רוצה תצטרך להתאים לצריכים ספציפיים שלך import numpy as np import math import sys from imageio import imread,imsave import matplotlib.pyplot as plt def norm_v(x,D): if x > D//2: return x-D return x def pcr(a,b): H,W = a.shape fft_a = np.fft.fft2(a) fft_b = np.fft.fft2(b) cb = np.ma.conjugate(fft_b) tmp = fft_a * cb tmp = tmp / np.absolute(tmp) res = np.fft.ifft2(tmp).real row,col = np.unravel_index(np.argmax(res),res.shape) return norm_v(row,H),norm_v(col,W) def _crop(img): C,H,W = img.shape S=480 dy = (H - S) // 2 dx = (W - S) // 2 return img[1,dy:dy+S,dx:dx+S] def main(): images=[] for n in range(300,0,-1): p = "%s%d.ppm" % (sys.argv[2], n) img = np.array(imread(p)).transpose(2, 0, 1).astype(np.float32) img = (img/255)**2.2 images.append(img) print("Loaded ",p) im0 = _crop(images[0]) x,y,r=[],[],[] for i in range(1,300): dy,dx=pcr(im0,_crop(images[i])) dr = math.sqrt(dx*dx+dy*dy) x.append(dx) y.append(dy) r.append(dr) print("Processing %d r=%f" % (i,dr)) #plt.plot(np.arange(0,299)/30.0,np.maximum(r,0.25)) plt.yscale('log') v = np.array([0] + r + [0]) avg=np.zeros((300-5,)) mxv=np.zeros((300-5,)) ctr = None for dr in range(5): if dr == 2: ctr = v[dr:dr+len(avg)] avg += v[dr:dr+len(avg)] mxv = np.maximum(v[dr:dr+len(avg)],mxv) avg/=5 scale = 4.64 plt.plot(np.arange(0,len(mxv))/30.0,np.maximum(1,mxv*scale)) plt.plot(np.arange(0,len(mxv))/30.0,np.maximum(1,ctr*scale)) minpoint = np.argmax(mxv*scale >= 10) / 30.0 maxpoint = np.argmax(mxv*scale >= 60*5) / 30.0 plt.plot([minpoint,minpoint],[1,np.max(r)*scale]) plt.plot([maxpoint,maxpoint],[1,np.max(r)*scale]) plt.xlabel('Time from end s') plt.ylabel('Distance from stable arcsec') plt.title("%s\n5' to 10'' = %3.1fs" % (sys.argv[1],maxpoint - minpoint)) plt.grid(True) plt.savefig(sys.argv[1].replace(' ','_').replace(',','').lower() + '.png') plt.show() if __name__ == "__main__": main()