From b6c7c400eef6ef799af630835375e48a179a6622 Mon Sep 17 00:00:00 2001
From: "SK. SALAHUDDIN" <127140921+sksalahuddin2828@users.noreply.github.com>
Date: Tue, 11 Jul 2023 17:22:22 +0600
Subject: [PATCH] Create 3D_Earth_Rotation_and_Satellite_Locations.py

---
 3D_Earth_Rotation_and_Satellite_Locations.py | 47 ++++++++++++++++++++
 1 file changed, 47 insertions(+)
 create mode 100644 3D_Earth_Rotation_and_Satellite_Locations.py

diff --git a/3D_Earth_Rotation_and_Satellite_Locations.py b/3D_Earth_Rotation_and_Satellite_Locations.py
new file mode 100644
index 0000000..9420a18
--- /dev/null
+++ b/3D_Earth_Rotation_and_Satellite_Locations.py
@@ -0,0 +1,47 @@
+import matplotlib.pyplot as plt
+from mpl_toolkits.basemap import Basemap
+from matplotlib.animation import FuncAnimation
+import warnings
+
+# Suppress warning messages
+warnings.filterwarnings("ignore")
+
+# Set up the figure and axes
+fig = plt.figure()
+ax = fig.add_subplot(111)
+
+# Create a Basemap instance
+m = Basemap(projection='ortho', lat_0=0, lon_0=0, resolution='l')
+
+# Satellite orbit data
+satellites = [
+    {"lat": 20, "lon": 40},
+    {"lat": -30, "lon": 120},
+    {"lat": 40, "lon": -60},
+    {"lat": -50, "lon": 180},
+    {"lat": 10, "lon": -150},
+    {"lat": 60, "lon": 30},
+    {"lat": -70, "lon": -90},
+    {"lat": 0, "lon": 60},
+    {"lat": 30, "lon": -120},
+    {"lat": -20, "lon": 150},
+]
+
+def update(frame):
+    ax.cla()  # Clear the axes
+    m = Basemap(projection='ortho', lat_0=0, lon_0=frame, resolution='l')  # Update lon_0 for rotation
+    m.drawcoastlines(color='gray')
+    m.bluemarble()  # Display the Blue Marble image
+
+    # Plotting satellite positions
+    for satellite in satellites:
+        lat = satellite["lat"]
+        lon = satellite["lon"]
+        x, y = m(lon, lat)
+        m.plot(x, y, 'ro', markersize=5, latlon=False)
+
+    ax.set_title('3D Earth Rotation and Satellite Locations')
+
+# Set up the animation with a faster rotation speed (interval=20)
+ani = FuncAnimation(fig, update, frames=range(0, 360, 2), interval=20)
+plt.show()