如何用pyephem计算折射校正后的RA/DEC坐标？
使用pyephem计算折射校正后的RA/DEC坐标需要以下步骤：
- 导入pyephem库：在Python代码中导入pyephem库，以便使用其中的函数和类。
import ephem- 创建一个观测者对象：使用ephem.Observer()函数创建一个观测者对象，该对象表示观测者的位置和观测时间。
observer = ephem.Observer()- 设置观测者的位置和观测时间：通过设置观测者对象的属性来指定观测者的位置和观测时间。例如，可以使用
observer.lat和observer.lon属性设置观测者的纬度和经度，使用observer.date属性设置观测时间。
observer.lat = '纬度'
observer.lon = '经度'
observer.date = '观测时间'- 创建一个天体对象：使用ephem.FixedBody()函数创建一个天体对象，该对象表示要观测的天体。
body = ephem.FixedBody()- 设置天体的坐标：通过设置天体对象的属性来指定天体的坐标。例如，可以使用
body._ra和body._dec属性设置天体的赤经和赤纬。
body._ra = '赤经'
body._dec = '赤纬'- 计算折射校正后的坐标：使用观测者对象的
body.compute(observer)方法计算折射校正后的天体坐标。
body.compute(observer)- 获取折射校正后的坐标：使用天体对象的
body.a_ra和body.a_dec属性获取折射校正后的赤经和赤纬坐标。
corrected_ra = body.a_ra
corrected_dec = body.a_dec以上是使用pyephem计算折射校正后的RA/DEC坐标的基本步骤。pyephem是一个强大的天文计算库，可以用于计算天体的位置和运动，包括折射校正等。在实际应用中，可以根据具体需求进行更多的参数设置和计算操作。
腾讯云相关产品和产品介绍链接地址：
- 腾讯云服务器（CVM）：提供弹性计算能力，满足各类业务需求。详情请参考腾讯云服务器
- 腾讯云数据库（TencentDB）：提供高性能、可扩展的数据库服务。详情请参考腾讯云数据库
- 腾讯云人工智能（AI）：提供丰富的人工智能服务和解决方案，包括图像识别、语音识别、自然语言处理等。详情请参考腾讯云人工智能
- 腾讯云物联网（IoT）：提供全面的物联网解决方案，帮助连接和管理物联网设备。详情请参考腾讯云物联网
- 腾讯云存储（COS）：提供安全、可靠的对象存储服务，适用于各类数据存储需求。详情请参考腾讯云存储
- 腾讯云区块链（BCS）：提供高性能、可扩展的区块链服务，支持企业级应用场景。详情请参考腾讯云区块链
- 腾讯云元宇宙（Metaverse）：提供虚拟现实（VR）和增强现实（AR）技术支持，构建沉浸式体验。详情请参考腾讯云元宇宙
相关·内容
我尝试将对象的J2000 RA/DEC坐标转换为“观察位置”，即包含折射效果的地心RA/DEC坐标。按照文档()，我这样做了：
from math import pi
import ephem
from datetime import datetime
ra = 20.370473492 / 12. * pi
dec = 40.256674958 / 180. * pi
tt = datetime(2016, 07, 27, 23, 30, 0)
lowell = ephem.Observer()
lowell.lon = '-111:32.1'
lowell
浏览 2提问于2016-08-02得票数 1
回答已采纳
给定以下信息:观测者位于、经度、和纬度，由GPS给出。以方位和天顶观测天空中由GPS给出的时间上的物体。
我如何才能得到银河坐标的点在天空上？这能用PyEphem来完成吗？
我觉得这让我差点就到了：
from datetime import datetime
import ephem
observer = ephem.Observer()
observer.pressure = 0
observer.lon = str(4.95569830927)
observer.lat = str(52.3545603701)
observer.elevation = 56.8426188165
obse
浏览 1提问于2014-08-28得票数 1
回答已采纳
2回答
编辑:刚刚发现，我需要把地球上某个位置的纬度、经度和海拔转换成J2000坐标，与ra/dec或月球无关。对此我很抱歉。你的回答确实给了我很多深刻的见解。请参阅下面编辑的问题。
问:如何将纬度、经度和高程转换为J2000坐标(XYZ)。在爱姆有转换吗？我查了一下医生，但找不到我需要的东西(或者因为我在这个领域缺乏知识而忽略了一些东西)。谢谢
*****************旧(漠视) ******************
我有月球位置在右侧扬升(RA)和倾斜(Dec)，我想把它们转换成X，Y，Z坐标。有内置的PyEphem函数吗？还有，它背后的数学是什么？谢谢。
编辑:我使用的是J2000坐标
浏览 4提问于2012-08-14得票数 2
回答已采纳
我使用PyEphem已经有一段时间了，最近在我的软件上出现了一个问题，在这个问题中，尽管局部地平线上的高度是恒定的，但计算出的赤经并不是恒定的。我可以想象，只要局部地平线上的高度也是恒定的，那么赤纬将是恒定的，只有右赤经变化？ 我正在使用radec_of()函数和一个从经纬度坐标创建的观察者类，elevation和preassure设置为0。 self.QTH.lat = str(lat)
self.QTH.lon = str(lon)
self.QTH.pressure = 0
self.QTH.elevation = 0 和我的函数，从本地/水平坐标返回赤道坐标。 # Returns e
浏览 35提问于2021-06-29得票数 0
我使用pyEphem来计算卫星的RA/Decs，我被上面计算和描述的不同值搞糊涂了
这段代码
sat=ephem.readtle("SATNAME ", \
"1 38356U 12030A 14148.90924578 .00000000 00000-0 10000-3 0 5678",\
"2 38356 0.0481 47.9760 0002933 358.9451 332.7970 1.00270012 3866")
gatech = ephem.Observ
浏览 0提问于2014-05-31得票数 1
最近，MPCORB数据库改变了它的标准时期。然后pyephem没能计算出小行星的位置。
以2015-11-23.55UT上的小行星2001 BP25为例。使用新旧标准时期，我计算了它的位置如下：
import ephem
import datetime
tobs = ephem.Date( datetime.datetime(2015, 11, 23, 13, 12, 13) ) # observation time
obj = ephem.readdb( "2001 BP25 0,e,3.47304,134.76688,294.08950,2.4271529,0.26065039
浏览 0提问于2015-11-26得票数 1
这看起来是一项非常基本的任务，但我无法解决它。
PyEphem文档：
介绍如何使用.alt和.az属性中的仰角和方位角执行相反的转换，即从Body和Observer对象转换为外观拓扑中心位置。
但是，我应该如何从仰角和方位角开始，并获得RA和Dec？
例如，这里有一组坐标，我想在赤道坐标系中得到RA和Dec：
az = 3.30084818 #rad
el = 0.94610742 #rad
lat = 34.64 #deg
lon = -103.7 #deg
alt = 35800.26 #m
ut = 2455822.20000367 #julian date
谢谢!
浏览 4提问于2012-01-23得票数 4
回答已采纳
我正试着从alt/azi到RA/Dec，在一个固定的地点，在天空中的一个点，尝试pyEphem。我尝试了几种不同的方法，我得到了一些正确的答案，在几度以内，但我期待更好，我不知道问题出在哪里。
我一直在用Canopus作为测试用例(我不是专门追求星星，所以我不能使用内置的目录)。所以就我而言，我知道
stn = ephem.Observer()
# yalgoo station, wa
stn.long = '116.6806'
stn.lat = '-28.3403'
stn.elevation = 328.0
stn.pressure = 0 # no re
浏览 2提问于2015-02-12得票数 3
我希望获得各种XYZ坐标中的对象坐标，而不是它们出现在天空中的位置。我感兴趣的是：
欧洲经委会："“(不是与地球自转)
ECEF："“(非惯性，也称为”地球中心自转“)
和某种太阳系或黄道坐标(太阳或重心中心)。
我已经安装了PyEphem，运行了几个示例，并开始理解什么是，但到目前为止，我发现的示例似乎集中在天空中的事物是如何出现的。
要么链接到一个有用的“如何”，要么一个简单的例子将是有帮助的。现在，我不知道如何处理标准选项：
In [55]: mars.
mars.a_dec mars.g_dec m
浏览 0提问于2015-12-30得票数 3
回答已采纳
我一直在用PyEphem来处理行星相对于地球上某一特定位置的运动。然而，我有时注意到，我从PyEphem得到的结果显示为不连续的衰落，而右扬升则是连续的。RA和Dec坐标每30分钟在这些地块上进行采集。
我认为恒星物体会以连续的方式运动，但是当倾斜是负的时，它看起来是不连续的。
知道为什么会这样吗？我可以张贴我的脚本，如果这也有帮助。
代码：
from ephem import *
import datetime
import re
def cSTI(number):
degrees = int(number[0])
minutes = int(number[1]
浏览 0提问于2018-09-11得票数 1
回答已采纳
2回答
我只想确保我理解PyEphem对angles做了什么。
假设我有一个观测者在30°S，60°W，今天中午看着天空。所以我这么做了：
obs = ephem.Observer()
obs.lat, obs.lon = '-30', '-60'
obs.date = '2012/6/22 12:00:00'
现在我想找出当观察者观察某个方位和高度时，他看到的是天球的哪一部分(在RA和赤纬中)。我可以这样做：
obs.radec_of(az,alt)
这就是我感到困惑的地方。据我所知，PyEphem适用于所有角度的弧度，所以我应该将az和alt作为弧度
浏览 1提问于2012-06-23得票数 4
回答已采纳
我在az/alt中生成视界上方的随机位置(az=0-360，alt=0-90)，并用radec_to()将其计算为RA和DEC。为了检查结果，我重新转换了它们。
所以我不明白的是，为什么我得到一半左右的坐标回到地平线下？
import ephem
from datetime import datetime
import random
home = ephem.Observer()
home.lon = '-70.4' # +E
home.lat = '-24.62' # +N
home.elevation = 1189 # meters
home.da
浏览 6提问于2021-06-09得票数 1
回答已采纳
2回答
在尝试进行坐标转换时，我对PyEphem有一个问题。在这个小例子中，我试着计算(l，b) = (52°,68.5°)的赤道坐标。
我认为它应该是这样工作的：
import ephem
galactic = ephem.Galactic(52.0/360.0,68.5/360.0)
equatorial = ephem.Equatorial(galactic, epoch=ephem.J2000)
print('%.13f %.13f' % (equatorial.ra*360.0, equatorial.dec*360.0))
输出是1640.9226879684597 -1
浏览 1提问于2015-03-12得票数 1
回答已采纳
我想计算五个拉格朗日点在天空中的位置(L1 - L5)，当提供协调世界时的时间戳和地球上的纬度/经度/高度观测站位置时，计算为RA/DEC或AZ/EL。我看过pyephem，pyorbital和skyfield，但据我所知，它们不支持这一点。如果有人能给我指点用python计算这个的参考/代码，那将不胜感激。
附注:如果不是pyephem/skyfield/pyorbital的一部分，我认为这将是一个有趣的添加功能？
浏览 4提问于2018-06-20得票数 1
有没有一种方法可以使用PyEphem来计算目标列表的上升/设置值？我的目标有一个坐标的astropy.table(在RA和Dec中)，我想使用PyEphem (或其他软件包)为观察者参考系中的每个对象生成另外两列，具有上升和设置时间。
我知道PyEphem有一个预先安装的物体目录，比如明亮的恒星和行星，但我想做的是利用它的location.next_setting(ephem.Object())和location.next_rising(ephem.Object())函数来计算我的目标列表何时可以观察到，而不仅仅是明亮的恒星和行星。
我知道有PyEphem的替代品，到目前为止，我已经尝试安装A
浏览 1提问于2016-03-18得票数 4
我无法解决为什么我得到一个不同的方位/高度之间的PyEphem和天域之间的彗星计算，但正确的上升和倾斜匹配。
代码示例：
#!/usr/bin/env python3
# External verification:
# https://in-the-sky.org/data/object.php?id=0088P
# https://theskylive.com/88p-info
# https://heavens-above.com/comet.aspx?cid=88P
# https://www.calsky.com/cs.cgi/Comets/3?
import datetim
浏览 7提问于2020-07-22得票数 0
回答已采纳
以下是Pyephem文档中提供的示例：
iss = ephem.readtle("ISS (ZARYA)","1 25544U 98067A 03097.78853147 .00021906 00000-0 28403-3 0 8652","2 25544 51.6361 13.7980 0004256 35.6671 59.2566 15.58778559250029")
gatech = ephem.Observer()
gatech.lon, gatech.lat = '-84.39733
浏览 0提问于2013-03-01得票数 1
回答已采纳
1回答
我有一组给定的数据，其中第一列和第二列分别是天空坐标(ra，dec)，第三列和第四列是笛卡尔坐标系中的坐标(x，y)。
我需要使用坐标x和y制作一个二维插值曲面，并使用Ra和Dec制作另一个曲面。问题是存在遮罩区域，如上图所示。我可以通过绘制丢失的数据来说明它们(在目录中有非NaN值)。这就是我到目前为止尝试过的，但没有给出正确的答案：
from scipy.interpolate import griddata
import numpy as np
import matplotlib.pyplot as plt
data = np.loadtxt('test.asc
浏览 0提问于2014-09-02得票数 3
1回答
在特定的日期和时间，如何计算星体在某个观测站(经度、纬度)的位置(ra，dec)？我需要在计算中包括的所有元素(恒星的适当运动，大气压力和温度)的全部坐标缩减.)
我试过了，但我不确定我能不能完成。
import ephem
polaris = ephem.readdb("Polaris,f|M|F7,2:31:48.704,89:15:50.72,2.02,1")
polaris.compute('2016/3/1 16:22:56')
print polaris.a_ra
print polaris.a_dec
我也尝试了天体计划，我认为我更接近于解
浏览 5提问于2016-04-12得票数 2
回答已采纳
目前，我正在制作一个脚本，在给定一组天体坐标的情况下，它将在接下来的几天告诉您，该点将在特定的望远镜中可见。标准很简单，在水平坐标系中，物体的高度必须在30到65度之间(这里的变量"crit“代表这个高度，但以弧度表示)。因此，我为望远镜设置了一组名为"Ant“的参数，然后使用Pyephem：
#imported ephem as ep
obj= ep.FixedBody()
obj._ra= E.ra
obj._dec= E.dec
obj._epoch = E.epoch
Ant.date = ep.now()
for d in range(days):
for i
浏览 3提问于2016-07-12得票数 1
回答已采纳
我正在使用PyEphem高精度地计算月球的未来位置，并对量化不确定性并将其作为脚本的一部分返回感兴趣。PyEphem中是否有返回身体位置不确定性的功能，或者我是否需要自己计算它们？
我正在考虑月球在J2000，RA和DEC的位置，并试图获得从当前时间到未来1个月的月球位置。
谢谢!
浏览 1提问于2018-06-19得票数 0
1回答
我有一个望远镜的文本文件，有三列，RA，DEC和TIME。RA和DEC是天空坐标，而时间是在该坐标下观察的时间。虽然RA和DEC在许多条目中重复，但每个条目中的时间不是，即，一个可以指向(RA，DEC)=(120,30)，首先2小时，然后再4小时，最后再1小时。其思想是使用python编写另一个文本文件，其中每对坐标只出现一次，时间是该坐标中所有时间的总和。当然，有些坐标只观察一次，所以它们没有重复。 我做了一些尝试，目前我设法编写了相同的输入文件，但每个坐标的第一个输入总是好的结果。问题是，如果一对坐标重复5次，文本文件将写入该坐标5次，即第一个坐标中的时间和，其余的我不知道，但它们不应该
浏览 18提问于2019-06-06得票数 0
回答已采纳
我正在重构一些使用PyEphem来使用Skyfield的旧代码，并且我得到了一个人的GHA/dec结果的细微差别。
def sf(year):
from skyfield.api import N, W, wgs84
from skyfield.api import load
from skyfield.units import Angle
ts = load.timescale()
eph = load('de421.bsp')
greenwich = wgs84.latlon(58.47722 * N, 0.0 * W)
浏览 21提问于2022-04-09得票数 0
回答已采纳
我在地球上有一个定向天线，我想知道它在ra，dec坐标中指向哪里。我想使用新的skyfield库，因为pyephem已经不再开发了。
很容易计算出我在地球上的topos位置：
planets = load('de421.bsp')
earth = planets['earth']
paris = earth + Topos('48.839059 N', '2.310147 E')
但是我不知道如何从这一点来表示az，alt对。我见过from_altaz方法，但我不能让它工作。
如果我从topos中尝试：
antenna = par
浏览 1提问于2019-02-22得票数 2
我正在尝试设置pyephem来做一个正确的上升和下降，然后用它来计算物体的高度和方位。然而，我在分配这些特征时遇到了问题。下面是我要做的：
obj = ephem.FixedBody()
print (x) #x was assigned to 5.575
obj._ra = x
print (obj._ra)
此命令的输出为：
5.575
21:17:41.75
据我所知，我应该以小时为单位输入RA值，以度为单位输入Dec值。我就是这么做的。但是由于某些原因，我将5.575赋给了变量obj._ra，而这个变量并没有保存正确的值。我的偏差值也有同样的问题。有没有人能解释一下我哪里错了？
浏览 2提问于2016-05-09得票数 1
1回答
我有，我想在我的MCMC代码中使用它。关键是实现的速度，以避免减慢我的马尔可夫链蒙特卡罗抽样。The problem：在目录中，我在第一和第二列中有两个名为ra和dec的参数，它们是天空坐标
data=np.loadtxt('Final.Cluster.Shear.NegligibleShotNoise.Redshift.cat')
ra=data[:,0]
dec=data[:,1]
然后，在七列和八列中，X和Y位置，即网格坐标，它们是网格空间中的点。
Xpos=data[:,6]
Ypos=data[:,7]
在我编写的函数中，需要调用一个Xcenter和Ycenter位置
浏览 3提问于2014-08-28得票数 1
回答已采纳
1回答
我正在寻找一种方法来定位在我的FITS图像上的像素坐标，它对应于物体在度上的ra和dec位置，在过采样后。如果我不是过度抽样的话，这会很简单，但我必须这么做。给出一个未经修改的图像，我可以：
from astropy.wcs import WCS
ra, dec = (43.603, 31.029)
w = WCS('myimage.fits')
x, y = w.all_world2pix(ra, dec, 1) #this gives me the pixel coordinates of the object at (ra, dec) position
然而，当我过采样
浏览 0提问于2018-08-16得票数 3
2回答
识别表中的RA和Dec重复项
、、、、
我正在尝试创建一个由特定区域中的对象组成的表。我正在收集irac，2 2mass和wise数据。我已经将所有数据合并到一个表中，现在正在尝试消除基于RA和Dec坐标的重复数据。我的最终目标是拥有一个没有任何副本的完整表。因此，如果对象在irac数据中，则没有与2mass或wise数据等效的对象。我刚开始使用蟒蛇。如果能得到复制件的行#那就太好了。
import pandas as pd
df = pd.read_csv('filename.csv')
duplicate = []
for num in df['ra' and 'dec']:
浏览 0提问于2019-07-25得票数 2
回答已采纳
3回答
我在天空中有一套一万点的收盘点。它们是在空中使用RA (右提升)和DEC (斜度)绘制的。在绘制图时，它们呈圆形。
我想要做的是，将圆圈分割成8个等份的，然后一次删除每个部分，然后使用其余的部分进行一些计算。
为了做到这一点，我想到了这个例子，即用弧线对它们进行切片。
我知道，弧的方程式是由：
S = r * theta
哪里
r --> radius
theta --> angle (in our case 45 degrees)
我想做这样的事：
slice1 = []
for a,b in zip(ra,dec):
if a>some value an
浏览 3提问于2015-07-14得票数 6
回答已采纳
1回答
我正在使用galpy package (http://galpy.readthedocs.io/en/latest/orbit.html)执行orbit集成。我想通过改变orbit的势模型和初始条件来执行许多轨道积分(~1000000 integrations)。对于每个轨道集成，我希望访问可观测坐标系(ra，dec，dist，pm_ra，pm_dec，vrad)中的Orbit，基于太阳的distance和velocity，以(degrees, kpc, mas yr-1, km/sec)为物理单位。
getOrbit()函数允许即时访问Orbit，但仅在R,phi,z,Vr,Vt,vz坐标系
浏览 3提问于2018-05-14得票数 0
我用AstroPy SkyCoord来完成从赤道坐标到银河系坐标的转换，有上百万的数据，它非常慢。任何人都有加快速度的想法，否则就需要永远运行整个数据集。守则如下：
from astropy import units as u
from astropy.coordinates import SkyCoord
import numpy as np
ra1 = np.loadtxt('data.txt',usecols=(0,))
dec1 = np.loadtxt('data.txt',usecols=(1,))
size = len(ra1)
for i in
浏览 2提问于2016-03-22得票数 2
回答已采纳
1回答
我正在努力寻找最好的方法，从坐标列表中删除不需要的目标。我的坐标(Ra，Dec)是用astropy.coordinates.SkyCoord形成的，但我有大量无法观测的目标，它们的倾斜度太低，所以我要做的是对我的列表进行排序，删除所有倾斜度低于-10度的目标(比如我的望远镜在北半球)。
这是我的代码行，它生成列表，它称为radecs，用于简化，并从天体球面坐标中获取Ra & Dec。
radecs = astropy.coordinates.SkyCoord(ra=phi*u.rad, dec=(0.5*np.pi - theta)*u.rad)
这是一个例子，说明我的列表是如何在Py
浏览 9提问于2016-03-14得票数 1
回答已采纳
我想要的是:我试着读取一个包含6,000个坐标(ra，and dec)的列表，对于每一个坐标，它们的周围都是78个点。我正在申请一个角度(ang)，然后试图找到新的RA和DEC。z_sq=(x2 + y2)似乎有问题，因为我得到了不支持的操作数类型(S)，用于**或pow()：‘’和'int‘。
import matplotlib.pyplot as plt
import numpy as np
import pylab as py
coords=np.genfromtxt('HETDEX_reg.txt',dtype=None,usecols=(0,1,2),nam
浏览 1提问于2015-08-07得票数 0
这是这个问题的后续问题：pyephem problem with alt/az to ra/dec and back 我想在观测者上方随机分布恒星。我在最后一个问题中做到了这一点。 我现在得到的问题是，星星似乎向中心聚集。 我怎么才能摆脱这个呢？ import matplotlib.pyplot as plt
import numpy as np
import math
import ephem
from datetime import datetime
import random
home = ephem.Observer()
home.lon = '-70.4'
浏览 137提问于2021-06-14得票数 1
回答已采纳
2回答
优化二维阵列
、、、、
我用Python2.7编写代码。我有两个坐标元组的二维数组。
array1 = [[[00_RA,00_DEC] [01_RA,01_DEC] ... [0N_RA,0N_DEC]]
[[10_RA,10_DEC] [11_RA,11_DEC] ... [1N_RA,1N_DEC]]
...
[[M0_RA,M0_DEC] [M1_RA,M1_DEC] ... [MN_RA,MN_DEC]]]
array2 = [[[00_ra,00_dec] [01_ra,01_dec] ... [0n_ra,0n_dec]]
浏览 3提问于2015-06-12得票数 2
回答已采纳
我正在寻找PyEphem python模块，以便为占星术研究找到准确的月亮位置。
import ephem
'Setting Delhi as location with 77E12 longitude, 28:36 lattitude
Delhi = ephem.Observer()
Delhi.lon = ephem.degrees('77:12')
Delhi.lon
Delhi.lat = ephem.degrees('28:36')
Delhi.date = '2015/3/22
浏览 1提问于2015-07-19得票数 0
2回答
我正绕着GroupBy转，但我仍然需要一些帮助。假设我有一个DataFrame，列为Group，给出对象组号、一些参数R和球面坐标RA和Dec。下面是一个模拟DataFrame：
df = pd.DataFrame({
'R' : (-21.0,-21.5,-22.1,-23.7,-23.8,-20.4,-21.8,-19.3,-22.5,-24.7,-19.9),
'RA': (154.362789,154.409301,154.419191,154.474165,154.424842,162.568516,8.355454,8.346
浏览 7提问于2017-09-19得票数 3
回答已采纳
我试图使用矩阵乘法将一些矩阵的ra和dec坐标转换为银河系坐标，但我得到了一个错误。代码如下：
import pandas as pd
import numpy as np
from astropy.io import fits
import astropy
import matplotlib.pyplot as plt
import math
fits_filename="Gaia_DR2/gaiadr2_20pc.fits"
df=pd.DataFrame()
with fits.open(fits_f
浏览 14提问于2019-08-02得票数 0
回答已采纳
2回答
我正在使用numpy和astropy用Python编写代码。在代码中，我希望创建与我的数据集相似的随机numpy数组。之后，我想把这些球面坐标的随机阵列转换成笛卡尔坐标。不幸的是，我不断地得到一个值错误，我完全不明白为什么会发生这种情况--我试图做一些虚拟测试，比如它们的形状相同，所有数组的值都是合理的，类型相同，等等，但是我被卡住了。这是我的代码：
from astropy.coordinates import SkyCoord
from astropy import units as u
import numpy as np
R = 445 + np.random.randn(26
浏览 1提问于2018-06-26得票数 0
回答已采纳
1回答
我想绘制一些二维地图，显示天文源在天空中的位置。
即：
我有一些数据ASCII表与赤道坐标的来源(RA -右提升，DEC -倾斜).举个例子：
这里RA以小时(h)、分钟(m)和秒(S)为单位；DEC -以度(°)、分钟(')和秒(")为单位。
现在，我正在用TOPCAT绘制地图，看起来如下：
除了TOPCAT的工作不能脚本化之外，一切都很好。但这是至关重要的，因为我有大量的数据。
我决定使用matplotlib (带有astropy)，但发现在HH:MM:SS和°：‘：中设置轴是一个很难的问题，如上面的图中所示。
“不讲道理的教程”说：
轴对象ax知道需
浏览 2提问于2017-03-30得票数 0
我试图用projection aitoff来绘制两个坐标数组(大约480万个元素)，比如：
fig = plt.figure(figsize=(12,10))
ax = fig.add_subplot(111, projection="aitoff")
ax.grid(alpha=0.3)
ax.hist2d(ra, dec, bins = 100)
但我得到了以下错误：
TypeError: Changing axes limits of a geographic projection is not supported. Please consider using Car
浏览 4提问于2022-03-23得票数 0
回答已采纳
基本目标:给定一组天空坐标(ra，dec)，这个Python函数应该返回以这些坐标为中心的天空图像。
详细目标:该函数将使用天文查询和simbad执行坐标查询。理想情况下，它将能够访问以(ra，dec)为中心的给定区域内的对象数据库。因此，还将输入一些大小参数(FoV，或radius)。但是，当通过浏览器对simbad进行查询时，有一个选项可以绘制由AladinLite生成的以输入坐标为中心的星区，并且函数应该返回这个图像(希望是一个适合的图像)。
我尝试过的是:我在这里所做的一切都发现了天体查询和simbad的使用。但是，我在文档中找不到任何关于如何通过天文查询获取来自AladinLite的
浏览 3提问于2020-06-11得票数 1
回答已采纳
我需要计算给定日期/时间的视方位角和仰角以及亚月经度。我得到的az/el角通常与其他来源(MoonCalc.org、Horizons等)一致。但是，对于次月经/次月经，还没有很好的比较来源。更重要的是，我怀疑我使用dec/ra值得到的经度/经度，因为ra在很长一段时间内几乎没有变化。
这是我所做的基本决定：
roc.date='2018/1/1 01:00:00'
moon=ephem.Moon(roc)
print('rocMoonTest: %s UTC-4, lat/lon = %0.4f [+N], %0.4f [+E]' %
(roc.date
浏览 0提问于2018-09-15得票数 1
我试图将数据库中每颗恒星的Ra (坐标)值与同一列表中其他恒星的所有其他Ra值进行比较，并仅提取1度内的恒星。我很好地上传了数据，并进行了转换。实际上，我的问题是将所有值与所有其他值进行比较。我为循环使用了一个布尔值。
for i,r in enumerate(lyrae):
q=r['ra']-lyrae['ra']
m=(q<1) & (q>0)
c=q[m]
print(r,c)
for i,r in enumerate(lyrae):
v=r['dec']-lyr
浏览 0提问于2017-11-07得票数 0
1回答
我正在使用Astropy.coordinates将两个天文星表与RA，DEC坐标进行匹配。
我可以通过遵循astropy documentation并执行以下操作来查找最近邻居的列表：
from astropy.coordinates import SkyCoord
from astropy import units as u
cat1 = SkyCoord(ra=ra1*u.degree, dec=dec1*u.degree)
cat2 = SkyCoord(ra=ra2*u.degree, dec=dec2*u.degree)
idx, d2d, d3d = cat1.match_
浏览 10提问于2018-09-04得票数 0
回答已采纳
1回答
有效地计算数千个坐标对之间的距离
、、、、
我在python中打开了一个目录，它有大约7万行不同对象的数据(ra、dec坐标和对象名称)。我还有另外一个列表，其中有大约15,000个感兴趣的对象，它们也出现在前面提到的目录中。对于这15,000个对象中的每一个，我想看看在70000个大列表中是否有任何其他对象在对象的10弧秒内具有ra，dec坐标。如果这是真的，我只想标记这个对象并转到下一个对象。然而，这一过程需要很长时间，因为目前感兴趣的对象之间的距离(15 000次中的70 000次)是计算出来的。这需要几天时间！我怎样才能更有效地完成同样的任务？下面是我当前的代码，其中all_objects是所有感兴趣的15000个对象名称的列表
浏览 1提问于2019-02-14得票数 2
我试图使用Simbad类从astroquery获取所有帧的坐标，就像 (基本数据部分)显示的那样
我有以下代码：
from astroquery.simbad import Simbad
def get():
Simbad.reset_votable_fields()
Simbad.remove_votable_fields('coordinates')
Simbad.add_votable_fields('ra(:;A;ICRS;J2000)', 'dec(:;D;ICRS;2000)')
Simbad.ad
浏览 5提问于2016-12-28得票数 4
回答已采纳
1回答
备注：我在这个网站上看到了所有类似的问题，但我仍然无法做到这一点。
我为天文坐标创建了一个类，并试图重载<<操作符，以便在类中存储用户输入。这两个坐标(正确的上升和倾斜)是双星的矢量。
当我调用参数化构造函数时，我得到以下两个错误：“调用类类型的对象而没有适当的运算符()或将函数转换为指针-函数类型”和“不匹配对‘(坐标)的调用”(std::vector&，std::vector&)'“。似乎我必须通过参考传递向量，但我没有设法。我该怎么做呢？
以下是代码：
#include <iostream>
#include <vector>
浏览 5提问于2022-05-11得票数 1
2回答
我在中有一组点，我需要将它们投影到一个平面上，即天象或。
astropy显然能够执行这种类型的投影，。
问题是，我不知道如何将这些模块应用于我的数据，而且文档也相当稀少。
对于本例，假设可以通过以下方法生成赤道ra, dec坐标(以度数表示)：
ra = np.random.uniform(130., 135., 1000)
dec = np.random.uniform(-55., -57., 1000)
浏览 5提问于2016-02-22得票数 4
回答已采纳
1回答
绘制非标准轴
、、、、
我正在用RA和Dec坐标绘制一幅图像。图像与(RA，Dec) = (0,0)点重叠。对于y轴(斜度)，这是很好的，但对于x轴，以度工作，我发现很难正确地绘制这个图。
我希望从(RA，Dec) = (358，-2)开始，从4平方度到(2，2)。有没有办法让python/matplotlib理解x轴上的360度=0度？
我就是这样开始的：
lowerleft_ra = 358.
lowerleft_dec = -2.
lx = lowerleft_ra
ly = lowerleft_dec
uy = lowerleft_dec + 4.
ux = lowerleft_ra + 4.
plt.im
浏览 1提问于2016-01-17得票数 1
回答已采纳
扫码
添加站长 进交流群
领取专属 10元无门槛券
手把手带您无忧上云
相关资讯
热门标签
云服务器
ICP备案
对象存储
实时音视频
即时通信 IM活动推荐
腾讯云开发者
