Questions on the astronomical applications of Mathematica.

learn more… | top users | synonyms

39
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++, ...
23
votes
3answers
3k views

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

Given the following world images: ...
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 ...
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 ...
19
votes
1answer
263 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 ...
15
votes
2answers
223 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
279 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 ...
13
votes
3answers
1k 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
2answers
377 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
2answers
342 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
2answers
369 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 ...
11
votes
1answer
260 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, ...
9
votes
1answer
181 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 ...
9
votes
1answer
259 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
401 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
343 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 ...
8
votes
2answers
421 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 ...
7
votes
2answers
247 views

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

What astronomical event do AstronomicalData's "NextRiseTime", "NextSetTime", ...
7
votes
2answers
537 views

Only final result from NDSolve

Finally I started to play with differential equations in Mathematica. And I have faced the problem, which seems too me so basic that I'm afraid this question is going to be closed soon. However, ...
7
votes
2answers
579 views
7
votes
0answers
92 views

MoonPosition and TimeZone support

With this question, I noticed some some weird results of MoonPosition. Consider the following code: ...
6
votes
2answers
485 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 ...
6
votes
1answer
205 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. ...
6
votes
0answers
111 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, ...
5
votes
2answers
262 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
223 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
132 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 ...
4
votes
1answer
42 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. ...
4
votes
1answer
466 views

Combining Gravity Turn and Orbit Models

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

Obsolete AstronomicalData - how to convert to PlanetData

As of Mathematica version 10, AstronomicalData is obsolete, deprecated in favor of PlanetData, ...
4
votes
0answers
106 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 ...
2
votes
2answers
83 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 ...
2
votes
1answer
165 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 ...
2
votes
2answers
859 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 ...
2
votes
0answers
42 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: ...
1
vote
1answer
204 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 ...
1
vote
2answers
105 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 ...
1
vote
2answers
127 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 ...
1
vote
1answer
70 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 ...
1
vote
1answer
124 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 ...
1
vote
0answers
81 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 ...
1
vote
0answers
571 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 ...
0
votes
1answer
104 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 ...
0
votes
0answers
46 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 : ...
0
votes
0answers
71 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 ...
0
votes
0answers
115 views

Correctly set up FindMinimum

I am trying to setup the following minimization to align two astronomical images. ...
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 ...