Skip to main content
756 questions
Newest Active Bountied Unanswered
Filter by
Sorted by
Tagged with
0 votes
0 answers
54 views

A scatter plot with an astropy wcs projection is limited to just the xlim/ylim of the axis, even if there are points outside that range: import numpy as np import matplotlib.pyplot as plt import ...
Matthew's user avatar
  • 2,310
0 votes
0 answers
96 views

Been investigating the 3I/Atlas (or C/2025 N1) rabbit hole. Of note is its ''anti-tail'' with some arguing it isn't a perspective effect. So, to test this I'm trying to use a ''syndyne-synchrone'' ...
John Manden's user avatar
  • 1
0 votes
0 answers
232 views

I’m trying to plot constellations with Python using Gaia star data, but the constellation lines don’t look like they do in Stellarium. I’ve tried plotting in both Alt/Az and RA/Dec, even used a ...
Ankit Sharma's user avatar
  • 1
4 votes
2 answers
110 views

Creating a SkyCoord object: c_end = SkyCoord( x=xf, y=yf, z=zf, v_x=vxf, v_y=vyf, v_z=vzf, frame='galactocentric' ) and attempting to put this into an Astropy QTable with ...
Gwald4's user avatar
  • 41
1 vote
2 answers
171 views

I am trying to create a slightly complex 2D plot with astronomical data, using Matplotlib, and Astropy. It involves a astronomical image, overlaid with a scatter plot, and Matplotlib patches. Using ...
aparajito's user avatar
  • 123
2 votes
1 answer
100 views

For example, in this code I plot a astronomy FITS image: import numpy as np from astropy.io import fits from astropy.visualization import simple_norm from astropy.wcs import WCS import matplotlib....
Firestar-Reimu's user avatar
  • 177
-1 votes
1 answer
68 views

I've been working on this for two weeks. I have a data frame with 3,600 columns and 200 rows, which correspond to a frequency sweep performed with a detector. The file's original format is .fit, so I'...
M1ctl4nt3cutl1's user avatar
  • 69
0 votes
1 answer
36 views

I keep getting KeyError: 'd_ρ' on my d_rho = ... line. I am trying to run the differentials so that I can put them in a table column (or single column per variable), to use for calculating velocity. ...
strueb's user avatar
  • 1
0 votes
1 answer
93 views

I am writing some code dealing with an aircraft locations with respect to an observer located at ground and I got a discrepancy between the distance calculated directly from geocentric Cartesian ...
raulgh's user avatar
  • 11
1 vote
2 answers
125 views

I am trying to do the transfomation from ITRS to topocentric AltAz frame using SkyCoord.transform_to() functionality: altazCoord = itrsCoord.transform_to(AltAz(obstime=time, location=obsGeoLocation) ...
MatejZigo's user avatar
  • 11
0 votes
1 answer
68 views

import astropy.units as u import numpy as np from astropy.io import fits from astropy.wcs import WCS import matplotlib.pyplot as plt from astropy.visualization import simple_norm from astropy....
Firestar-Reimu's user avatar
  • 177
2 votes
1 answer
142 views

I am trying to do a scatter plot of a few data points on a region of the sky and then add a circle to the plot. Note, that I don't want a plot of the full sky, just of the region around the data ...
NeStack's user avatar
  • 2,062
0 votes
1 answer
202 views

The thing is i'm working with two .fits images, but the pixscale of each one is different. Let's say Image 1 has a CDELT = +- 0.0013888888888889 [deg] and Image 2 has a CDELT = +- 0.00277778000000 [...
mr_starspotting's user avatar
  • 1
1 vote
1 answer
96 views

I'm using astropy, matplotlib, skyfield and astroquery to create sky charts with Python. With skyfield I do the stereographic projection of the sky. With astropy I use WCS, something like this: fig = ...
Unix's user avatar
  • 2,844
-1 votes
1 answer
347 views

My Python script is: import numpy as np import matplotlib.pyplot as plt import cv2 from fil_finder import FilFinder2D from astropy import units skeleton = cv2.imread("./zhang.png", 0) fil = ...
Abukeram's user avatar
  • 3
1 vote
1 answer
104 views

I use astropy to make calculations with units and I normally do these calculations in a jupyter notebook so the results are printed nicely. I would like to know how to control the number of digits ...
Shai Avr's user avatar
  • 1,440
2 votes
2 answers
499 views

I'm making a script in Python to get information for all objects from the NGC and IC catalogs. Actually, I already have this information from OpenNGC, however, coordinates don't have the same ...
Unix's user avatar
  • 2,844
2 votes
2 answers
854 views

I'm trying to generate star maps with the equatorial coordinates system (RAJ2000 and DEJ2000). However, I only get a grid system where meridians and parallels are in parallel, while parallels should ...
Unix's user avatar
  • 2,844
1 vote
1 answer
116 views

I have a list of object names without any spaces, but the Simbad.query_object of astroquery requires spaces. For example, with space works... see below: from astroquery.simbad import Simbad ...
HP AS's user avatar
  • 109
1 vote
1 answer
191 views

The sequence: fit_something = fitting.LevMarLSQFitter() fitted = fit_something(initial_guess_and_model,x,y,acc=accuracy,maxiter=maximumum_iterations) Will provide a fitted model, and in the ...
Claudio Pedrazzi's user avatar
  • 25
1 vote
1 answer
205 views

Background I want to convert the header information of a FITS data into a dictionary and save it as a binary file using Python. And of course, the binary data has to be loaded later. Problem The ...
jipparchus's user avatar
  • 11
1 vote
1 answer
1k views

I am facing a challenge that should be straightforward but has proven to be quite complex with the Astropy library. Specifically, I need to frequently compute the longitude and latitude from a ...
ysBach's user avatar
  • 361
0 votes
1 answer
318 views

I want to find the date in which an object will be at its maximum altitude in the sky for a set location. Given an object's right ascension and declination, an observers location and a range of dates, ...
Juan Martin Hunicken's user avatar
  • 1
0 votes
1 answer
323 views

I have used iraf in the past to do this, and I have no problem when the pixel scale does not vary (although I have to manually edit the task parameters for it to work correctly). Now I am working with ...
Carl's user avatar
  • 1
0 votes
1 answer
301 views

I have a list of MJD time. I need to convert it in list UTC but only the date (year/month/day) and no time, so a list that I can use to make a plot with time in x-axis. Thanks I tried with astropy and ...
Federica Giacchino's user avatar
  • 1
1 vote
1 answer
165 views

I'm running the latest version of Astropy in Python using a Mac. I'd like to change the directory where Astropy caches files. I'm running: astropy.config.set_temp_cache(path='/Users/josephprendergast/...
JoeP's user avatar
  • 11
0 votes
0 answers
130 views

I am trying to access the SDSS catalog through the FITS file I downloaded here, I am looking at the DR16Q_v4.fits file. Unfortunately I cannot seem to find where the actual spectrum data is located ...
Ducktective's user avatar
  • 1
2 votes
1 answer
886 views

I've followed the instructions from questions present on this site pertaining to this issue, but the code all seems off. For one thing, the altitude seems far too low, and it doesn't appear that the ...
TheForWarden's user avatar
  • 23
1 vote
1 answer
562 views

I want to transform the icrs coordinates (ar and dr) from epoch 1991.25 to epoch 2016.0, also in the icrs system. But I got this error. It's my first time using AstroPy. <SkyCoord (ICRS): (ra, dec) ...
B612's user avatar
  • 25
1 vote
2 answers
3k views

When I imported aplpy, I got this error: ImportError Traceback (most recent call last) Input In [12], in <cell line: 9>() 7 from astropy.coordinates import SkyCoord ...
laklak's user avatar
  • 11
0 votes
0 answers
288 views

I have the x,y,z positions of a gps satellite from its precise ephemerides and I am not sure how I can convert from these cartesian ECEF coordinates (WGS-84 G2139 IERS2010 frame specifically) to RA ...
woplda123's user avatar
  • 1
-1 votes
1 answer
93 views

There is no way to get the uncertainty from the Centroids (photutils.centroids) function. I need to be able to retrieve uncertainty due to using the tool to measure velocities requiring an amount of ...
Istwan's user avatar
  • 1
1 vote
1 answer
268 views

I will start this post by saying that although I can code a little, I'm not super good at it yet so apologies if the answer is really obvious and I do not see it. I have a series of .fits files that I ...
juno's user avatar
  • 11
0 votes
0 answers
122 views

I am making a 2d simulation of the solar system. For that you need a starting point for all planets. Right now I am storing that data in this form (ignore the color, radius, mass and Max_Trail_Length ...
Tmpecho's user avatar
  • 11
0 votes
0 answers
62 views

I have an astronomical image fits format with WCS (World Coordinate System) in its header. I want to do some morphological operations on the image and understandably I want to correct the WCS. Say I ...
niaei's user avatar
  • 2,659
0 votes
0 answers
123 views

I'm new to Python and just started learning this. I'm trying to figure out how to work with the Jupyter Notebook for searching for an object name in online databases. I wrote two versions of the code (...
Andy Ayr's user avatar
  • 134
2 votes
0 answers
836 views

I'm trying to import astropy.coordinates.representation in a Jupyter notebook and I get a circular import error: > ImportError: cannot import name 'SphericalRepresentation' from partially ...
Konstantinos Topaloglou's user avatar
  • 21
0 votes
1 answer
206 views

I'm trying to retrieve b and v magnitudes for a specific star from the Sloan Digital Sky Survey (SDSS) database, but I'm facing some challenges with my code. Here are my requirements: I need to input ...
yoongified's user avatar
  • 1
0 votes
0 answers
878 views

I'm trying to install poliastro on windows 10 through pip. I have pip updated to the latest version; however, every time I try to install poliastro, I keep getting the error message in the title, ...
Big Dave's user avatar
  • 1
0 votes
1 answer
543 views

I construct a table using Astropy as follows: data = astropy.table.Table(names=['Time', 'Rate (C/s)']). Then, I fill up the rows under the aforementioned columns using data.add_row. Finally, I intend ...
SKR's user avatar
  • 37
1 vote
1 answer
912 views

Can I use Astroquery/Astropy to get the radius and mass of solar system objects? Or another common library? I've been trying to use astroquery.horizons I've tried elements() and ephemeris() but no ...
DrDino's user avatar
  • 23
2 votes
0 answers
157 views

The little program below is almost a copy-paste for an example form Astropy website. The program works fine under PyCharm environment but when trying to make it executable using PyInstaller, I'm ...
Michel Michaud's user avatar
  • 47
0 votes
0 answers
224 views

I have been attempting to overplot a Quadrangle on a map that I have to replicate a slit with data on the sky, but I need to rotate the Quadrangle around some center point and have not been able to ...
Andrew's user avatar
  • 13
1 vote
2 answers
139 views

I want to generate a number of Gaussian sources in a fits array. My code is below. from astropy.io import fits from numpy import * xx=yy=1024 xc=xx/2.;yc=yy/2. A=10.0 gau=zeros([xx,yy]) for i in range(...
Viral Parekh's user avatar
  • 17
1 vote
1 answer
190 views

I'm trying to read multiple CSV files from nightly observations. Each CSV file has a time and magnitude column and I want to be able to stack all my nightly files into a single Timeseries object. It ...
Lipto's user avatar
  • 13
-3 votes
1 answer
279 views

I'm working on a CNN quasar detector that takes the optical spectrum as input, and classifies it as Quasar or non-quasar. I'm currently gathering data from this site:https://dr18.sdss.org/optical/...
Anirudh Menon's user avatar
  • 25
1 vote
0 answers
129 views

I have some astronomical light curves for which I'm trying to recover some periodicities. One of them (from TESS) has many thousand samples across ~23 days, with a ~1-day gap in the middle. The ...
John Wendeborn's user avatar
  • 93
2 votes
1 answer
780 views

I'm trying to plot regions using astropy, but at the same time plot points as scatter plot. But I'm not managing to understand what is happening with the coordinates or the matplotlib artists. First ...
Euler's user avatar
  • 151
0 votes
1 answer
424 views

I am working with astronomical data and need to work with custom units, specifically in Mega Jansky per steradian. I would like to create these custom units using the astropy.units library and also ...
user19416164's user avatar
  • 11
0 votes
0 answers
89 views

I am getting this error trying to read data from a FITS file. I can see that the file contains the correct data because DS9 viewer shows the image I would expect. The error I'm getting is below. Can ...
McMac Music's user avatar
  • 1

1
2 3 4 5
…
16