Astropy: वेग के साथ फ्रेम और एक "इकाई" प्रतिनिधित्व कार्टेशियन में परिवर्तित करने में विफल रहता है

को निर्मित 23 दिस॰ 2017 · 3टिप्पणियाँ · स्रोत: astropy/astropy

#6944 में कुछ कोने-केस परीक्षणों ने उन फ़्रेमों के लिए कुछ समस्याग्रस्त व्यवहार का खुलासा किया जिनमें वेग है जहां वास्तविक प्रतिनिधित्व UnitSpherical । यह पता चला है कि यह सच है। नीचे दिए गए मामले इस मुद्दे को स्पष्ट करते हैं:

यह भी ध्यान दें कि यदि #6944 को मर्ज किया गया है, तो इस बग के लिए वर्कअराउंड के साथ कुछ स्थान हैं। यदि यह बग ठीक कर दिया गया है तो उन वर्कअराउंड को हटा दिया जाना चाहिए (उन्हें खोजने के लिए coordinates उप-पैकेज #7028 लिए खोजें)।

यह सही है:

>>> f = ICRS(1*u.deg, 2*u.deg)
>>> f.cartesian
<CartesianRepresentation (x, y, z) [dimensionless]
    ( 0.99923861,  0.01744177,  0.0348995)>

लेकिन ये दो मामले नहीं हैं:

>>> f = ICRS(1*u.deg, 2*u.deg,
          pm_dec=1*u.mas/u.yr, pm_ra_cosdec=2*u.mas/u.yr, radial_velocity=10*u.km/u.s)
>>> f.represent_as('cartesian', in_frame_units=True)
UnitConversionError: 'mas / (rad yr)' (frequency) and 'km / s' (speed) are not convertible

>>> f = ICRS(1*u.deg, 2*u.deg, 
                pm_dec=1*u.mas/u.yr, pm_ra_cosdec=2*u.mas/u.yr)
>>> g = f.transform_to(GCRS)
>>> g.cartesian
UnitConversionError: '1 / s' (frequency) and 'km / s' (speed) are not convertible

यह शायद एक उपयोगी सुराग है कि दूसरे असफल मामले में, f.cartesian काम करता है , लेकिन जीसीआरएस में परिवर्तन में कुछ ऐसा है जो इसे विफल करना शुरू कर देता है। यह संभव है कि ये सूक्ष्म रूप से अलग-अलग मुद्दे हों, लेकिन वे इतने करीब हैं कि उन्हें शायद एक साथ हल करना होगा।

@adrn या @mhvk , कोई विचार है कि यहां क्या हो रहा है? कोई यह तर्क दे सकता है कि कम से कम उपरोक्त परिवर्तन चरण कुछ हद तक खराब परिभाषित है, लेकिन किसी को भी कार्टेशियन प्रतिनिधित्व प्राप्त करने में सक्षम होना चाहिए, भले ही इकाइयां थोड़ी अजीब हों ...

Bug coordinates

स्रोत

eteq

सभी 3 टिप्पणियाँ

हाय @eteq ,

मैं इस मुद्दे पर काम करने की कोशिश कर रहा था। लेकिन उदाहरण के मामले में,


>>> f = ICRS(1*u.deg, 2*u.deg,
          pm_dec=1*u.mas/u.yr, pm_ra_cosdec=2*u.mas/u.yr, radial_velocity=10*u.km/u.s)
>>> f.represent_as('cartesian', in_frame_units=True)

उदाहरण में कोई distance नहीं है, बल्कि एक radial_velocity । मुझे समझ में नहीं आया कि इसका वास्तव में क्या मतलब है। इसके अलावा, अगर मैं फ्रेम में कुछ दूरी रखता हूं, तो यह पूरी तरह से काम करता है!

क्या आप कृपया मुझे इस मुद्दे को समझा सकते हैं, मैं भ्रमित हो रहा हूं।

shreyasbapat 1 अग॰ 2019

@shreyasbapat क्या आप कृपया यह देखने की कोशिश कर सकते हैं कि क्या https://github.com/astropy/astropy/pull/9064 इस समस्या को ठीक करता है? (संपादित करें: यह नहीं करता है)

astrojuanlu 2 अग॰ 2019

👎2

9086 ने मुझे इस मुद्दे की याद दिला दी। दो मामलों के लिए:

कोई दूरी नहीं लेकिन RV . के साथ

>>> f = ICRS(1*u.deg, 2*u.deg,
          pm_dec=1*u.mas/u.yr, pm_ra_cosdec=2*u.mas/u.yr, radial_velocity=10*u.km/u.s)
>>> f.represent_as('cartesian', in_frame_units=True)
UnitConversionError: 'mas / (rad yr)' (frequency) and 'km / s' (speed) are not convertible

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

आर.वी. समन्वय परिवर्तन द्वारा पेश किया गया

>>> f = ICRS(1*u.deg, 2*u.deg, 
                pm_dec=1*u.mas/u.yr, pm_ra_cosdec=2*u.mas/u.yr)
>>> g = f.transform_to(GCRS)
>>> repr(g)
# error
>>> g.data
<CartesianRepresentation (x, y, z) [dimensionless]
    (0.99923901, 0.01742676, 0.03489572)
 (has differentials w.r.t.: 's')>
g.data.differentials
<CartesianDifferential (d_x, d_y, d_z) in 1 / s
    (6.07736084e-13, -1.85772092e-11, -8.12625661e-12)>

यह कुछ ज्यादा ही सूक्ष्म है। उचित गति के बिना, g UnitSpherical , जो समझ में आता है, और यहां यह भी होना चाहिए, इसलिए परिवर्तन में कुछ गड़बड़ है। पूर्णता के लिए, इकाई क्षेत्र पर इसका प्रतिनिधित्व करना काम करता है:

g.represent_as('unitspherical')                                                               
Out[24]: 
<UnitSphericalRepresentation (lon, lat) in rad
    (0.01743826, 0.0349028)
 (has differentials w.r.t.: 's')>

In [25]: g.represent_as('unitspherical').differentials['s']                                            
Out[25]: 
<UnitSphericalDifferential (d_lon, d_lat) in rad / s
    (-1.85963079e-11, -8.13120751e-12)>

लेकिन फिर भी, g.cartesian को भी काम करना चाहिए - यहां, समस्या यह है कि फ्रेम उपयोगकर्ता के लिए तय करते हैं कि कौन सी इकाइयां समझदार हैं, इस बात पर जोर देते हुए कि कार्टेसियन वेग km/s (बेशक, यह वापस चला जाता है मेरी लंबे समय से चली आ रही शिकायत के लिए कि प्रतिनिधित्व में समन्वय इकाइयां मेरे लिए तय की गई हैं...): https://github.com/astropy/astropy/blob/c3dc7b38303b3615945aa36a11ebf7bd31d5cc0a/astropy/coordinates/baseframe.py#L1744 -L1753

इकाइयाँ स्वयं से प्राप्त होती हैं:

g.representation_info
{astropy.coordinates.representation.CartesianRepresentation: {'names': ['x',
   'y',
   'z'],
  'units': [None, None, None]},
 astropy.coordinates.representation.UnitSphericalRepresentation: {'names': ('ra',
   'dec'),
  'units': (Unit("deg"), Unit("deg"))},
 astropy.coordinates.representation.RadialRepresentation: {'names': ['distance'],
  'units': [None]},
 astropy.coordinates.representation.SphericalRepresentation: {'names': ('ra',
   'dec',
   'distance'),
  'units': (Unit("deg"), Unit("deg"), None)},
 astropy.coordinates.representation.PhysicsSphericalRepresentation: {'names': ['phi',
   'theta',
   'r'],
  'units': [Unit("deg"), Unit("deg"), None]},
 astropy.coordinates.representation.CylindricalRepresentation: {'names': ['rho',
   'phi',
   'z'],
  'units': [None, Unit("deg"), None]},
 astropy.coordinates.representation.CartesianDifferential: {'names': ('v_x',
   'v_y',
   'v_z'),
  'units': (Unit("km / s"), Unit("km / s"), Unit("km / s"))},
 astropy.coordinates.representation.UnitSphericalDifferential: {'names': ('pm_ra',
   'pm_dec'),
  'units': (Unit("mas / yr"), Unit("mas / yr"))},
 astropy.coordinates.representation.SphericalDifferential: {'names': ('pm_ra',
   'pm_dec',
   'radial_velocity'),
  'units': (Unit("mas / yr"), Unit("mas / yr"), Unit("km / s"))},
 astropy.coordinates.representation.UnitSphericalCosLatDifferential: {'names': ('pm_ra_cosdec',
   'pm_dec'),
  'units': (Unit("mas / yr"), Unit("mas / yr"))},
 astropy.coordinates.representation.SphericalCosLatDifferential: {'names': ('pm_ra_cosdec',
   'pm_dec',
   'radial_velocity'),
  'units': (Unit("mas / yr"), Unit("mas / yr"), Unit("km / s"))},
 astropy.coordinates.representation.RadialDifferential: {'names': ['d_distance'],
  'units': [None]},
 astropy.coordinates.representation.PhysicsSphericalDifferential: {'names': ['d_phi',
   'd_theta',
   'd_r'],
  'units': [None, None, None]},
 astropy.coordinates.representation.CylindricalDifferential: {'names': ['d_rho',
   'd_phi',
   'd_z'],
  'units': [None, None, None]}}

कोई यहां देखता है कि केवल CartesianDifferential ही गलत है, ऐसा इसलिए है क्योंकि इसमें पसंदीदा इकाइयां बिल्कुल सेट हैं।

जिसका शायद मतलब है कि #9086 में सुधार एक छोटे से संशोधन के साथ ठीक है: किसी को अनुरोधित इकाइयों में जाने का प्रयास करना चाहिए, लेकिन अगर यह विफल हो जाता है तो परवाह नहीं करनी चाहिए।

mhvk 2 मई 2020

क्या यह पृष्ठ उपयोगी था?

0 / 5 - 0 रेटिंग्स

Astropy: वेग के साथ फ्रेम और एक "इकाई" प्रतिनिधित्व कार्टेशियन में परिवर्तित करने में विफल रहता है

सभी 3 टिप्पणियाँ

9086 ने मुझे इस मुद्दे की याद दिला दी। दो मामलों के लिए:

कोई दूरी नहीं लेकिन RV . के साथ

आर.वी. समन्वय परिवर्तन द्वारा पेश किया गया

संबंधित मुद्दों