Questions tagged [astronomy]

Questions on the astronomical applications of Mathematica.

Filter by
Sorted by
Tagged with
49
votes
4answers
4k 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++, ...
36
votes
2answers
3k 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 ...
29
votes
2answers
9k 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 ...
28
votes
1answer
1k views

Interstellar Image Effects

I'd like to know how to create this type of gravitational lensing effect in an image (here are some star fields from the hubble) "I wrote down the equations and I tested them in Mathematica by ...
26
votes
4answers
7k views

Composition: how to make a day and night world map?

Given the following world images: ...
21
votes
1answer
1k 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, ...
20
votes
1answer
517 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 ...
19
votes
3answers
2k 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, I'...
19
votes
2answers
1k 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 ...
17
votes
2answers
5k views

What is the closest star to Earth other than the sun?

How to find the Sun's nearest neighbor star using Mathematica? I tried but it didn't work
16
votes
1answer
404 views

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

AbsoluteTime 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:...
15
votes
3answers
3k 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 ...
15
votes
2answers
991 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 ...
14
votes
1answer
703 views

Speed up NDSolve compared to Python (calls to LSODA)

I migrated a numerical model code from Python to Mathematica and am surprised how much faster the Python version runs. Profiling of the Python version tells me that it is about 100 times faster (120 ...
14
votes
1answer
670 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 ...
13
votes
5answers
4k views

Beautiful dance of the Earth and Venus around the sun

I have seen very nice dynamic graphics created by users of Ma. I would like to get help whether the following graphics can be created using Ma. I was not able to create the dynamic picture here. ...
13
votes
2answers
562 views

Processing images from NASA's Juno Jupiter Orbiter

Question Can you help me improve what I’ve done to create an enhanced-color picture from an image provided by NASA’s JunoCam? [My work is shown in the penultimate image of Jupiter in this post.] I’...
13
votes
1answer
1k 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 ...
13
votes
2answers
444 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 ...
12
votes
1answer
697 views

Can I create a radar chart of sunrise/sunset times for a year

Can I create a radar chart (or just a simple curve) of sunrise/sunset times, for a given location thougout a year? It would be something like this, but with a smooth curve: (source: hassanhameed.com)...
11
votes
2answers
555 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 ...
11
votes
3answers
494 views

How to speed up plotting with PlotPoints?

I'm plotting an ellipse 1/(1 - e Cos[x]) with e=0.9 using PolarPlot. Taking $x\in[0,2\pi)$, ...
11
votes
1answer
241 views

How can I get curated data to be expressed consistently in SI units?

AstronomicalData objects are consistently in SI units: ...
11
votes
1answer
394 views

Calculate moonrise

Mathematica has Sunrise[], but no Moonrise[]. But it has MoonPosition[], so I thought it ...
10
votes
2answers
530 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 ...
10
votes
1answer
367 views

Plotting how sunrise times change over the year

I created a somewhat lengthy Python code to show how the time of sunrise varies over a year, resulting in the following figure (daylight savings were disregarded): The dataset can be found here. The ...
10
votes
1answer
618 views

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

I made this Graphics3D of the Sun with the plotted orbits of a few astronomical objects. ...
10
votes
1answer
396 views

Finding planetary conjunctions with Mathematica (project-level) [closed]

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 ...
10
votes
1answer
343 views

3D orbits of moons around their respective planets

I have seen a lot of solar system models (both animated and not), but I am trying to replicate the same thing with moons. Either just one planet and its lunar system or the our whole solar system with ...
9
votes
1answer
333 views

How do I determine astronomical transit times?

Mathematica, through the "legacy" Scientific Astronomer package, used to have the ability to easily determine accurate transit times of astronomical objects simply, with ...
8
votes
2answers
651 views

Extracting and displaying data from GW150914 hdf5 files

Has anyone looked at the LIGO data using Mathematica or know how to do this? I'm not familiar with the hdf5 format but I've been able to import the data, consisting primarily of 131,072 floating ...
8
votes
2answers
519 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 ...
8
votes
2answers
746 views

Problems getting clean sunrise, sunset, civil, nautical and astronomical twilight data

If in a notebook you enter: ...
8
votes
1answer
407 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 ...
8
votes
3answers
505 views

Why is AstronomicalData so slow?

I'm having trouble getting some pretty simple functions using AstronomicalData to perform at anywhere near the speed I need them to. For example, I have many ...
8
votes
0answers
86 views

Extracting, Formatting, and Exporting Mathematica AstronomicalData - Is There an Easier Way?

If any of you have tried to use Mathematica AstronomicalData and its newest replacements like StarData, you may already be familiar with the large data issues, and the almost impossible download times....
7
votes
3answers
1k views

RK4 Gravity Simulator

I have the following RK4 solver which splits the two 2nd order ODEs, used to calculate x and y positions under the influence of a gravitating body where $$x''(t)=\frac{G m x(t)}{(x(t)^2+y(t)^2)^{3/2}}$...
7
votes
2answers
320 views

What do AstronomicalData's “SetTime” and “RiseTime” values refer to?

What astronomical event do AstronomicalData's "NextRiseTime", "NextSetTime", ...
7
votes
2answers
973 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 ...
7
votes
1answer
186 views

Mass distribution in NDSolve

I'm trying to model the solar system with NDSolve, basically a numerical n-body solution. I model the initial positions of bodies from a data set thusly: ...
7
votes
1answer
275 views

How to get star coordinates for a given time

Fixed in 10.3.1 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 ...
7
votes
1answer
122 views

StarData: proper motion of a star in HelioCoordinates

I am not an astronomer, but am trying to calculate the trajectory of an interstellar space probe, for example to Proxima Centauri. The HelioCoordinates seem appropriate, as typically used within the ...
7
votes
0answers
162 views

MoonPosition and TimeZone support

Bug introduced in 10.0 and fixed in 10.2 MoonPosition is new in 10.0 With this question, I noticed some some weird results of ...
6
votes
1answer
395 views

GeoGraphics - mapping the Moon

The link here shows a nice way to depict the map of the Moon in Mathematica. But, is there a way to center the Mollweide projection onto the far side of the Moon?
6
votes
2answers
550 views

Solving the path of Earth around Sun

This maybe isn't a universally helpful question. Maybe a little of a code dump. But here goes. I'm trying to solve the path Earth moves around Sun from Earths mass, Suns mass, Earths initial velocity ...
6
votes
1answer
131 views

Does SolarEclipse only cover a limited date range?

If I use the example from the documentation for SolarEclipse and simply replace the years there with older ones, say from around Shakespeare's time, ...
6
votes
2answers
599 views

Series solution of the Lane-Emden equation

I have the following problem: I'd like to show how the Lane-Emden equation would look like if someone would solve it as a series. Here's the code: ...
6
votes
1answer
150 views

TimeZone Anomaly

I ran into the following situation when I was trying to determine the time of sunrise (using the built-in Sunrise function) for a variable location on a variable ...
6
votes
1answer
217 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 ...
6
votes
1answer
376 views

Can NDSolve address spherical gravitational collapse?

As far as I know spherical gravitational collapse -of central importance to theoretical physics- is thought to be out of the scope of Wolfram Mathematica. However there are cases, as in this article,...