Questions on the astronomical applications of Mathematica.

learn more… | top users | synonyms

6
votes
1answer
112 views

Finding planetary conjunctions with Mathematica (project-level)

Summary: What's an efficient way of finding planetary conjuctions with Mathematica using the SPICE kernels, and without using AstronomicalData (which I believe is less efficient). NASA publishes ...
2
votes
0answers
119 views

Creating an orrery

I am interested in creating an orrery, ideally with SystemModeler. I have been experimenting with the MultiBody library since that allows you to do animations. The most relevant example seems to be ...
4
votes
1answer
267 views

Using ParametricNDSolve and NMinimize to find accurate interplanetary trajectory arrival periapse radius

I have an interplanetary trajectory simulation that calculates the required velocity (using a Lambert solver) and escape angle for a spacecraft to travel from Earth to Mars. This model does not use a ...
3
votes
1answer
107 views

Computing planet conjunctions with 2D circular orbits still hard?

All methods I've seen for computing planetary/stellar conjuctions (when two planets or a planet and a star appear close together in the sky because their angular separation, as measured from Earth, ...
10
votes
3answers
623 views

Only final result from NDSolve

Finally, I started to play with differential equations in Mathematica. And I have faced the problem, which seems to me so basic that I'm afraid this question is going to be closed soon. However, ...
9
votes
1answer
162 views

GeoProjection for astronomical data - wrong ticks

Here is an excellent answer to what I wanted to do, and it works like a charm. Unfortunately, when I want to change it a bit, it fails. For setting the stage, here is the datafile, and this is my ...
17
votes
2answers
459 views

Using map projections with astronomical data

I noticed that all important "Geoprojections" are available in projections for a spherical reference models: GeoProjectionData function; 1 - I am trying using the sinusoidal projection for ...
25
votes
3answers
3k views
1
vote
0answers
38 views

PlanetData “VelocityAroundSun” giving inconsistent answers

My goal is to create a ListPlot showing the Earth's orbital speed over the course of a year, with one datapoint per day. I'm seeing problems long before I even get to the point of plotting the data! ...
6
votes
1answer
118 views

ParametricNDSolve[ ] and NMinimize[ ] in Combination

I have the following code that is used to simulate a spacecraft orbiting Mars, and uses ParametricNDSolve to vary the spacecraft's initial conditions using a parameter value ...
5
votes
1answer
552 views

p-iteration algorithm to solve Lambert's problem for interplanetary trajectories

I've been reading up on using the p-iteration method to solve Lambert's problem for choosing the correct interplanetary trajectory between two planets given the ...
16
votes
1answer
383 views

Right Ascension and Declination in astronomical functions

In the process of answering this question, I was forced to confront the various astronomical coordinate systems used by Mathematica. Background In astronomy, positions of celestial objects (stars, ...
4
votes
0answers
69 views

Is there a problem with some of the curated Star Data?

I'm not really sure if I'm doing something wrong, or if there's a problem with some of the Mathematica curated data on stars. My problem is this: I'm trying to find the location (latitude and ...
19
votes
2answers
4k views

Creating a simulation of our Solar System

As the title says, my objective is to create a simulation of the motion of the planets in our Solar System using Mathematica. All the theoretical background regarding the equations of motion of ...
2
votes
2answers
99 views

Finding Minimium Periapse Distance Between Spacecraft and Planet

I have the following simple model that simulates a spacecraft orbiting Mars. What I'm hoping to do is be able to find the periapse radius of the spacecraft as I did below, except instead of using ...
4
votes
1answer
577 views

Combining Gravity Turn and Orbit Models

I have a mathematical model for the motion of an orbiting spacecraft about Earth: ...
2
votes
0answers
44 views

Inaccurate Heliocentric coordinates of PlanetaryMoonData

I have noticed that the heliocentric coordinates obtained with PlanetaryMoonData are very inaccurate. Fore example, consider the Galilean moons of Jupiter: ...
2
votes
1answer
80 views

How to find µ (Gravitational constant mass product) with PlanetData?

I would like to use PlanetData to get µ, the product of the Gravitational constant G and a ...
0
votes
0answers
48 views

how to make one plot for each line of a table of 127 datas (127 plots) , or one plot with 127 curves onto it?

Good evening, I encounter a difficult problem. I have 127 lines of data, with for each line, 9 columns. I imported from a txt file : ...
8
votes
2answers
431 views

Generate a list of symbols for the planets

I would like to generate a list of the symbols for the planets to include in a nice table with other data. How can I programmatically generate (eg. using a pure function) the symbols from a list of ...
19
votes
1answer
278 views

Segmentation of sunspots

I recently went for a walk with my camera and snapped a few pictures of the sun. The pictures didn't turn out too well, but this was to be expected without a telescope. I thought it might be a nice ...
7
votes
0answers
97 views

MoonPosition and TimeZone support

With this question, I noticed some some weird results of MoonPosition. Consider the following code: ...
4
votes
1answer
47 views

Accessing Coordinates of Astronomical Data, e.g. Declination of all Galaxies

I am trying to access the coordinate information of all galaxies, stars, etc, available to the new astronomical data incorparated to Mathematica. However, the suggested (documented) approach fails. ...
5
votes
1answer
135 views

Finding minimum fly-by radius between Mars and spacecraft from interpolating function

I've written an interplanetary trajectory solver/plotter that plots the path taken by a spacecraft on an Earth-Mars mission, but have run into a little trouble when the spacecraft actually reaches ...
9
votes
1answer
237 views

Geocentric, Earth-fixed, coordinates of Sun, Moon (and others objects eventually)

I know I can get "HelioCoordinates" of the Earth and the Moon with ...
7
votes
0answers
115 views

How to get star coordinates for a given time

In the documentation of StarData, there is an example that shows how to get the altitude of Sirius at a given time (in the horizontal coordinate system). However, ...
0
votes
0answers
86 views

Plotting stellar spectrum from FITS file

I would like to make a plot of wavelength v. relative frequency from the stellar spectrum data in a FITS file (in order to fit a curve to one of the spectral lines). I used ...
15
votes
2answers
246 views

Align stars in multiple images

I hope that you find this question interesting and that this will help others with similar problems as this configuration is not so particular. My problem I have multiple images of a particular ...
40
votes
3answers
2k views

Mollweide maps in Mathematica

Context In my field of research, many people use the following package: healpix (for Hierarchical Equal Area isoLatitude Pixelization) which has been ported to a few different languages (F90, C,C++, ...
4
votes
0answers
116 views

HEALPix grid on a sphere

Are there any Mathematica tools for using the HEALPix ( Hierarchical Equal Area isoLatitude Pixelation) grid on a sphere? It is commonly used in astronomy, but I have not found any Mathematica ...
1
vote
1answer
159 views

Astronomy: Transform coordinates from horizon to equatorial

I'm using Mathematica 10. Giving a longitude & latitude on Earth, giving a local time, how can I transform Horizon Coordinte (az , alt) to Equatorial Coordinate (RA , DEC). Use function so I can ...
14
votes
1answer
285 views

Is there a way to get AbsoluteTime to correctly report absolute time differences?

Absolute time supposedly reports the "total number of seconds" between two dates, but it does not. For example because of the leap second applied at 2012-06-30T23:59:60Z ...
0
votes
1answer
105 views

AstronomicalData not including units?

I like using mathematica to do my physics homework, and AstronomicalData is useful for that, but it doesn't seem to work properly with units. [the question is, at ...
1
vote
1answer
219 views

Combining solar system planets orbits and positions in Graphics3D and Show

I am building an animated model of the solar system, where I would like to display the OrbitPaths of the planets along with their ...
11
votes
2answers
352 views

How to get the distance of Jupiter from the sun?

I am aware of the command AstronomicalData["Jupiter", "Distance"] from which I can get the distance between the earth and Jupiter. Is there any way to use ...
13
votes
3answers
2k views

How do I compute sidereal time and JD?

There was once a package that provided a wide range of functions for computing astronomical information. This package now appears to be out of date; but as near as I can tell, many of the useful ...
5
votes
1answer
129 views

Obsolete AstronomicalData - how to convert to PlanetData

As of Mathematica version 10, AstronomicalData is obsolete, deprecated in favor of PlanetData, ...
1
vote
2answers
111 views

Problem evaluating a complicated integral

Apologies for stating my problem poorly in the first instance, thank you for the help in narrowing down the issue. Problem For those interested, it's for the exact unbiased inverse of the ...
6
votes
2answers
512 views

Combining three R, G, B matrices into a single matrix of 3-tuples?

Working with Mathematica 8 on OSX here. I have three $1024 \times 1024$ matrices representing the counts from consecutive exposures onto a CCD with $R$, $G$, and $B$ filters on it, and am looking to ...
2
votes
2answers
941 views

Circular Orbit Velocity in Vector Components

I'm trying to find out how to calculate the initial velocity needed to keep a spacecraft in a circular orbit around the Earth given its initial position and an initial angle theta in a Cartesian ...
6
votes
1answer
213 views

How can I get this Graphics3D to rotate around a particular sphere?

I made this Graphics3D of the Sun with the plotted orbits of a few astronomical objects. ...
0
votes
0answers
118 views

Correctly set up FindMinimum

I am trying to setup the following minimization to align two astronomical images. ...
2
votes
1answer
168 views

Using DeleteCases with multiple AstronomicalData properties

I am attempting to use AstronomicalData and DeleteCases together to remove invalid cases of two properties while listing the ...
1
vote
2answers
129 views

Computing lists using AstronomicalData

Edit: Initially in my question I asked how I could search the database of AstronomicalData["Exoplanets"] and return 2 filtered values that only contained valued ...
8
votes
2answers
345 views

Impact Crater Simulation

I have been trying to model crater formation on a given planetary surface - $500\,\textrm{km}^2$. The locations of impacts are random, however, if an impact is within $30\,\textrm{km}$ of another, the ...
0
votes
0answers
240 views

Earth to Mars Interplanetary Travel Escape Velocity

I'm trying to create a script that simulates interplanetary travel from Earth to Mars using a Hohmann transfer, but I'm having a few difficulties getting the correct escape velocity. I'm using an x-y ...
5
votes
1answer
227 views

AstronomicalData and Planetary Heliocentric (x,y,z) Velocity Components

I'm trying to find data on the heliocentric velocities of planets in our solar system and was directed to Mathematica's AstronomicalData in a previous topic. ...
1
vote
0answers
582 views

Solar System Orbital Parameters [duplicate]

I recently made an n-body simulation and thought it would be interesting to try and model it on the interior of our solar system (Sun to Mars), but I cannot find initial conditions for such an ...
19
votes
1answer
1k views

The Orbit and Perigee of the Flamsteed comet

Historical context This year we have the 330-th anniversary of the Battle of Vienna - one of the great formative events of European history, it took place on September 12, 1683. Kara Mustafa, Grand ...
13
votes
2answers
389 views

How is AstronomicalData interpreting its date arguments?

It is easy to confirm that the AstronomicalData is not interpreting its arguments as UTC (or "GMT", which is, in any case, ambiguous) dates by verifying that the ...