MCP Stargazing

🚀 MCP - Stargazing

Calculate the altitude angle, rising time, and setting time of celestial bodies (the sun, moon, planets, stars, and deep - space objects) at any location on Earth, and support optional light pollution analysis.

🚀 Quick Start

MCP - Stargazing is a powerful tool that enables users to calculate various celestial parameters and conduct light pollution analysis. It supports multiple celestial body types and is timezone - aware, providing flexibility for users around the world.

✨ Features

• Altitude/Azimuth Calculation: Obtain the elevation and direction of celestial bodies.
• Rise/Set Time Calculation: Determine the time when a celestial body appears above the horizon and disappears.
• Light Pollution Analysis: Load and analyze light pollution maps (in GeoTIFF format).
• Supported Celestial Body Types:
  • Solar system objects (the sun, moon, planets)
  • Stars (e.g., "sirius")
  • Deep - space objects (e.g., "andromeda", "orion_nebula")
• Timezone - Aware: Support local time or UTC time.

📦 Installation

pip install astropy pytz numpy astroquery rasterio geopy

💻 Usage Examples

Basic Usage

Calculate Altitude/Azimuth

from src.celestial import celestial_pos
from astropy.coordinates import EarthLocation
import pytz
from datetime import datetime

# Observer location (New York)
location = EarthLocation(lat=40.7128, lon=-74.0060)

# Time (local time with timezone)
local_time = pytz.timezone("America/New_York").localize(datetime(2023, 10, 1, 12, 0))
altitude, azimuth = celestial_pos("sun", location, local_time)
print(f"Sun position: Altitude={altitude:.1f}°, Azimuth={azimuth:.1f}°")

Calculate Rise/Set Time

from src.celestial import celestial_rise_set

rise, set_ = celestial_rise_set("andromeda", location, local_time.date())
print(f"Andromeda: Rise time={rise.iso}, Set time={set_.iso}")

Load Light Pollution Map

from src.light_pollution import load_map

# Load a light pollution map in GeoTIFF format
vriis_data, bounds, crs, transform = load_map("path/to/map.tif")
print(f"Map bounds: {bounds}")

📚 Documentation

API Reference

celestial_pos(celestial_object, observer_location, time) (src/celestial.py)

• Input Parameters:
  • celestial_object: Name (e.g., "sun", "andromeda").
  • observer_location: EarthLocation object.
  • time: datetime (with timezone) or Astropy Time.
• Return Value: (Altitude in degrees, Azimuth in degrees)

celestial_rise_set() (src/celestial.py)

• Input Parameters:
  • celestial_object: Name.
  • observer_location: EarthLocation object.
  • date: Date object.
• Return Value: (Rise time, Set time)

Test Information

By default, all tests will run. To run the tests, execute:

pytest tests/

🔧 Technical Details

Future Work

• Support more celestial body types.
• Improve calculation accuracy.
• Add a graphical user interface.