Questions on the astronomical applications of Mathematica.

learn more… | top users | synonyms

43
votes
3answers
3k 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++, ...
30
votes
2answers
2k 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 ...
25
votes
3answers
4k views

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

Given the following world images: ...
22
votes
2answers
5k 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 ...
19
votes
1answer
584 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, ...
19
votes
1answer
256 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 ...
19
votes
1answer
342 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 ...
17
votes
2answers
666 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 ...
16
votes
3answers
870 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, ...
15
votes
2answers
349 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 ...
15
votes
1answer
304 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 ...
14
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 ...
13
votes
1answer
2k 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
13
votes
2answers
405 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 ...
11
votes
5answers
2k 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. ...
11
votes
2answers
395 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
1answer
344 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: Can it be done in ...
11
votes
1answer
212 views

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

AstronomicalData objects are consistently in SI units: ...
11
votes
1answer
464 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 ...
10
votes
2answers
430 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 ...
9
votes
1answer
221 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 ...
9
votes
1answer
212 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 ...
9
votes
1answer
286 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 ...
9
votes
3answers
430 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
2answers
449 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
649 views
7
votes
2answers
269 views

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

What astronomical event do AstronomicalData's "NextRiseTime", "NextSetTime", ...
7
votes
2answers
218 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 ...
7
votes
2answers
580 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
193 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 ...
7
votes
1answer
282 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. ...
7
votes
1answer
158 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
0answers
136 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 ...
7
votes
0answers
111 views

MoonPosition and TimeZone support

With this question, I noticed some some weird results of MoonPosition. Consider the following code: ...
6
votes
2answers
263 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
0answers
83 views

Heliocentric coordinates in PlanetData

I tried to use planets' HeliocentricLongitude and I wonder if anyone has any insight on what coordinate system they use. I took vernal equinox times from http://ns1763.ca/equinox/eqindex.html and ...
5
votes
2answers
278 views

Why is Neptune missing from AstronomicalData?

An old notebook I've recently started working with has a strange error: Prior to version 8, a figure I generated using Drop[AstronomicalData["Planet"],-1] to ...
5
votes
1answer
291 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?
5
votes
2answers
195 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: ...
5
votes
1answer
257 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. ...
5
votes
1answer
163 views

Obsolete AstronomicalData - how to convert to PlanetData

As of Mathematica version 10, AstronomicalData is obsolete, deprecated in favor of PlanetData, ...
5
votes
1answer
682 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 ...
5
votes
1answer
149 views

Confused by result of MoonPosition when running on Linux [closed]

I've just switched from Windows to Linux and am running Mathematica 10.1 on Ubuntu 14.04. MoonPosition[] Gives me: ...
5
votes
1answer
139 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: ...
5
votes
1answer
146 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 ...
5
votes
0answers
92 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 ...
4
votes
2answers
135 views

How to find the Andromeda Galaxy

I would like to find the Andromeda Galaxy ($M31$) with the unaided eye. One way is to use the following image and graphic but I can’t find it just by looking in a dark sky. It is hard to detect ...
4
votes
1answer
332 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 ...
4
votes
1answer
1k views

Combining Gravity Turn and Orbit Models

I have a mathematical model for the motion of an orbiting spacecraft about Earth: ...
4
votes
1answer
115 views

Impact Crater Sim. with added criteria

I have been trying with out success to edit this impact crater simulation script into producing a variable which calculates crater density. In addition I have been trying to figure out how to ...