8 | {/*
23 | );
24 | }
--------------------------------------------------------------------------------
/Orbit_Predictor-BackEnd/api/migrations/0003_remove_collision_conjunction_id_and_more.py:
--------------------------------------------------------------------------------
1 | # Generated by Django 5.1.3 on 2024-12-07 01:21
2 |
3 | import django.db.models.deletion
4 | from django.db import migrations, models
5 |
6 |
7 | class Migration(migrations.Migration):
8 |
9 | dependencies = [
10 | ('api', '0002_user'),
11 | ]
12 |
13 | operations = [
14 | migrations.RemoveField(
15 | model_name='collision',
16 | name='conjunction_id',
17 | ),
18 | migrations.RemoveField(
19 | model_name='probabilitycalc',
20 | name='conjuction_id',
21 | ),
22 | migrations.AddField(
23 | model_name='collision',
24 | name='cdm',
25 | field=models.ForeignKey(default=0, on_delete=django.db.models.deletion.CASCADE, related_name='collisions', to='api.cdm'),
26 | preserve_default=False,
27 | ),
28 | migrations.AddField(
29 | model_name='probabilitycalc',
30 | name='cdm',
31 | field=models.ForeignKey(default=0, on_delete=django.db.models.deletion.CASCADE, related_name='probability_calcs', to='api.cdm'),
32 | preserve_default=False,
33 | ),
34 | migrations.DeleteModel(
35 | name='Conjunction',
36 | ),
37 | ]
38 |
--------------------------------------------------------------------------------
/on-orbit-frontend/src/components/csa/csa.tsx:
--------------------------------------------------------------------------------
1 | import Link from "next/link"
2 |
3 | const Csa = () => {
4 | return(
5 |
6 |
7 |
Canadian Space Agency
8 |
9 |
10 |
Contact the
11 | CSA
12 |
13 |
14 |
15 |
About the
16 | CSA
17 |
18 |
19 |
20 |
Notices and copyright
21 |
22 |
23 |
24 |
25 | )
26 | }
27 |
28 | export default Csa
--------------------------------------------------------------------------------
/Orbit_Predictor-BackEnd/api/serializers/organization_serializer.py:
--------------------------------------------------------------------------------
1 | # organization_serializer.py
2 | from django.contrib.auth import get_user_model
3 | from rest_framework import serializers
4 | from ..models import Organization
5 | from .cdm_serializer import CDMSerializer # Adjust the import as needed
6 | from ..permissions import IsAdmin
7 |
8 | User = get_user_model()
9 |
10 | class OrganizationSerializer(serializers.ModelSerializer):
11 | users = serializers.PrimaryKeyRelatedField(
12 | many=True,
13 | queryset=User.objects.all()
14 | )
15 | cdms = CDMSerializer(many=True, read_only=True)
16 |
17 | def update(self, instance, validated_data):
18 | request = self.context.get("request", None)
19 | # Only admin (staff) users can override the alert_threshold or users
20 | if ('alert_threshold' in validated_data) or ('users' in validated_data):
21 | admin_permission = IsAdmin()
22 | if not admin_permission.has_permission(request, None):
23 | # If the request user fails the admin check, raise a ValidationError
24 | raise serializers.ValidationError({
25 | "detail": "Only admin users can modify alert_threshold or manage users."
26 | })
27 |
28 | return super().update(instance, validated_data)
29 |
30 | class Meta:
31 | model = Organization
32 | fields = ('id', 'name', 'alert_threshold', 'users', 'cdms')
33 |
--------------------------------------------------------------------------------
/on-orbit-frontend/package.json:
--------------------------------------------------------------------------------
1 | {
2 | "name": "on-orbit-frontend",
3 | "version": "0.1.0",
4 | "private": true,
5 | "scripts": {
6 | "dev": "next dev --turbopack",
7 | "build": "next build",
8 | "start": "next start",
9 | "lint": "next lint"
10 | },
11 | "dependencies": {
12 | "@emotion/react": "^11.13.3",
13 | "@emotion/styled": "^11.13.0",
14 | "@mui/material": "^6.1.6",
15 | "@supabase/supabase-js": "^2.48.1",
16 | "@types/d3": "^7.4.3",
17 | "@types/topojson-client": "^3.1.5",
18 | "cesium": "^1.126.0",
19 | "chart.js": "^4.4.8",
20 | "chartjs-adapter-date-fns": "^3.0.0",
21 | "class-variance-authority": "^0.7.1",
22 | "clsx": "^2.1.1",
23 | "d3": "^7.9.0",
24 | "d3-geo": "^3.1.1",
25 | "daisyui": "^4.12.20",
26 | "highcharts": "^12.1.2",
27 | "highcharts-react-official": "^3.2.1",
28 | "lucide-react": "^0.486.0",
29 | "next": "15.0.2",
30 | "react": "^18.3.1",
31 | "react-chartjs-2": "^5.3.0",
32 | "react-dom": "^18.3.1",
33 | "resium": "^1.19.0-beta.1",
34 | "satellite.js": "^5.0.0",
35 | "tailwind-merge": "^3.0.2",
36 | "tailwindcss-animate": "^1.0.7",
37 | "topojson-client": "^3.1.0"
38 | },
39 | "devDependencies": {
40 | "@types/node": "^20",
41 | "@types/react": "^18",
42 | "@types/react-dom": "^18",
43 | "eslint": "^8",
44 | "eslint-config-next": "15.0.2",
45 | "postcss": "^8",
46 | "tailwindcss": "^3.4.1",
47 | "typescript": "^5"
48 | }
49 | }
50 |
--------------------------------------------------------------------------------
/Orbit_Predictor-BackEnd/api/permissions.py:
--------------------------------------------------------------------------------
1 | # api/permissions.py
2 |
3 | from rest_framework import permissions
4 |
5 |
6 | class IsAdmin(permissions.BasePermission):
7 | """
8 | Allows access only to admin users.
9 | """
10 | def has_permission(self, request, view):
11 | return bool(request.user and request.user.role == 'admin')
12 |
13 |
14 | class IsCollisionAnalyst(permissions.BasePermission):
15 | """
16 | Allows access only to collision analysts.
17 | """
18 | def has_permission(self, request, view):
19 | return bool(request.user and request.user.role == 'collision_analyst')
20 |
21 |
22 | class IsUser(permissions.BasePermission):
23 | """
24 | Allows access only to regular users.
25 | """
26 | def has_permission(self, request, view):
27 | return bool(request.user and request.user.role == 'user')
28 |
29 |
30 | class CanViewCDM(permissions.BasePermission):
31 | """
32 | Custom permission to allow users to view CDMs based on their role and CDM privacy.
33 | - Admins and Collision Analysts can view all CDMs.
34 | - Regular Users can only view CDMs where privacy is True.
35 | """
36 | def has_object_permission(self, request, view, obj):
37 | if request.user.role in ['admin', 'collision_analyst']:
38 | return True
39 | return obj.privacy
40 |
41 | def has_permission(self, request, view):
42 | # Allow access only if the user is authenticated
43 | return bool(request.user and request.user.is_authenticated)
44 |
--------------------------------------------------------------------------------
/on-orbit-frontend/README.md:
--------------------------------------------------------------------------------
1 | This is a [Next.js](https://nextjs.org) project bootstrapped with [`create-next-app`](https://nextjs.org/docs/app/api-reference/cli/create-next-app).
2 |
3 | ## Getting Started
4 |
5 | First, run the development server:
6 |
7 | ```bash
8 | npm run dev
9 | # or
10 | yarn dev
11 | # or
12 | pnpm dev
13 | # or
14 | bun dev
15 | ```
16 |
17 | Open [http://localhost:3000](http://localhost:3000) with your browser to see the result.
18 |
19 | You can start editing the page by modifying `app/page.tsx`. The page auto-updates as you edit the file.
20 |
21 | This project uses [`next/font`](https://nextjs.org/docs/app/building-your-application/optimizing/fonts) to automatically optimize and load [Geist](https://vercel.com/font), a new font family for Vercel.
22 |
23 | ## Learn More
24 |
25 | To learn more about Next.js, take a look at the following resources:
26 |
27 | - [Next.js Documentation](https://nextjs.org/docs) - learn about Next.js features and API.
28 | - [Learn Next.js](https://nextjs.org/learn) - an interactive Next.js tutorial.
29 |
30 | You can check out [the Next.js GitHub repository](https://github.com/vercel/next.js) - your feedback and contributions are welcome!
31 |
32 | ## Deploy on Vercel
33 |
34 | The easiest way to deploy your Next.js app is to use the [Vercel Platform](https://vercel.com/new?utm_medium=default-template&filter=next.js&utm_source=create-next-app&utm_campaign=create-next-app-readme) from the creators of Next.js.
35 |
36 | Check out our [Next.js deployment documentation](https://nextjs.org/docs/app/building-your-application/deploying) for more details.
37 |
--------------------------------------------------------------------------------
/Orbit_Predictor-BackEnd/api/views/organization_views.py:
--------------------------------------------------------------------------------
1 | # organization_views.py
2 | from rest_framework import viewsets
3 | from rest_framework.decorators import action
4 | from rest_framework.permissions import IsAuthenticated
5 | from rest_framework.filters import SearchFilter
6 | from rest_framework.response import Response
7 |
8 | from ..models import Organization
9 | from ..serializers import OrganizationSerializer, UserSerializer, CDMSerializer
10 |
11 | class OrganizationViewSet(viewsets.ModelViewSet):
12 | """
13 | A viewset that provides the standard actions for the Organization model,
14 | plus custom endpoints to list users and CDMs.
15 | """
16 | queryset = Organization.objects.all()
17 | serializer_class = OrganizationSerializer
18 | permission_classes = [IsAuthenticated]
19 |
20 | filter_backends = [SearchFilter]
21 | search_fields = ['name']
22 |
23 | @action(detail=True, methods=['get'], url_path='users')
24 | def get_users(self, request, pk=None):
25 | """
26 | Returns all users belonging to this organization.
27 | """
28 | organization = self.get_object() # Retrieve the Organization instance
29 | serializer = UserSerializer(organization.users.all(), many=True)
30 | return Response(serializer.data)
31 |
32 | @action(detail=True, methods=['get'], url_path='cdms')
33 | def get_cdms(self, request, pk=None):
34 | """
35 | Returns all CDMs associated with this organization.
36 | """
37 | organization = self.get_object()
38 | serializer = CDMSerializer(organization.cdms.all(), many=True)
39 | return Response(serializer.data)
40 |
--------------------------------------------------------------------------------
/Orbit_Predictor-BackEnd/api/models/organization.py:
--------------------------------------------------------------------------------
1 | from django.db import models
2 | from django.conf import settings
3 | from django.db.models import Q
4 | from .cdm import CDM
5 |
6 | class Organization(models.Model):
7 | name = models.CharField(max_length=255, unique=True)
8 | alert_threshold = models.FloatField(default=1.0e-9)
9 | users = models.ManyToManyField(
10 | settings.AUTH_USER_MODEL,
11 | related_name='organizations',
12 | blank=True
13 | )
14 | cdms = models.ManyToManyField(
15 | CDM,
16 | related_name='organizations',
17 | blank=True
18 | )
19 |
20 | def save(self, *args, **kwargs):
21 | """
22 | When saving an Organization, automatically link any CDMs where this org's name
23 | matches sat1_operator_organization or sat2_operator_organization.
24 | """
25 | super().save(*args, **kwargs)
26 |
27 | # Find all CDMs where either sat1_operator_organization or sat2_operator_organization
28 | # matches this organization's name.
29 | matching_cdms = CDM.objects.filter(
30 | Q(sat1_operator_organization=self.name) | Q(sat2_operator_organization=self.name)
31 | )
32 | self.cdms.set(matching_cdms)
33 |
34 | def notify_if_collision_exceeds_threshold(self, collision):
35 | # TODO: Implement this method
36 | if collision.probability_of_collision >= self.alert_threshold:
37 | if collision.cdm in self.cdms.all():
38 | self.notify_users(collision)
39 |
40 | def notify_users(self, collision):
41 | # TODO: Implement this method
42 | for user in self.users.all():
43 | print(f"Notifying {user.username} about collision {collision}")
44 |
45 | def __str__(self):
46 | return self.name
47 |
--------------------------------------------------------------------------------
/on-orbit-frontend/src/components/about/about.tsx:
--------------------------------------------------------------------------------
1 | const About = () => {
2 | return (
3 |
4 |
5 |
About Us
6 |
7 |
8 | Welcome to the On-Orbit Collision Predictor. This project is a McGill Capstone project in collaboration with the Canadian Space Agency. Our goal is to provide accurate data on satellites in orbit and predict the probability of collisions to help ensure the safety and sustainability of space operations.
9 |
10 |
11 | Our team consists of dedicated students and professionals who are passionate about space and technology. We leverage advanced algorithms and data analysis techniques to deliver reliable predictions and insights.
12 |
13 |
14 | Thank you for visiting our site. If you have any questions or would like to learn more about our project, please feel free to contact us.
15 |
16 |
17 |
18 |
19 |
Team
20 |
21 |
Erik
22 |
Masa
23 |
Wasif
24 |
Yassine
25 |
26 |
27 |
28 | )
29 | }
30 |
31 | export default About
--------------------------------------------------------------------------------
/Orbit_Predictor-BackEnd/api/authentication.py:
--------------------------------------------------------------------------------
1 | # api/authentication.py
2 |
3 | import jwt
4 | from django.conf import settings
5 | from rest_framework import authentication, exceptions
6 | from .models import User
7 |
8 |
9 | class JWTAuthentication(authentication.BaseAuthentication):
10 | """
11 | Custom JWT Authentication class.
12 | """
13 |
14 | def authenticate(self, request):
15 | auth_header = authentication.get_authorization_header(request)
16 |
17 | if not auth_header:
18 | return None # No authentication credentials provided
19 |
20 | try:
21 | prefix, token = auth_header.decode('utf-8').split(' ')
22 | if prefix.lower() != 'bearer':
23 | return None # Invalid prefix
24 | except ValueError:
25 | raise exceptions.AuthenticationFailed('Invalid token header. No credentials provided.')
26 |
27 | return self.authenticate_credentials(token)
28 |
29 | def authenticate_credentials(self, token):
30 | try:
31 | payload = jwt.decode(token, settings.JWT_SECRET_KEY, algorithms=['HS256'])
32 | except jwt.ExpiredSignatureError:
33 | raise exceptions.AuthenticationFailed('Token has expired.')
34 | except jwt.InvalidTokenError:
35 | raise exceptions.AuthenticationFailed('Invalid token.')
36 |
37 | user_id = payload.get('user_id')
38 | role = payload.get('role')
39 |
40 | try:
41 | user = User.objects.get(id=user_id)
42 | except User.DoesNotExist:
43 | raise exceptions.AuthenticationFailed('User not found.')
44 |
45 | # Optionally, you can verify the role matches
46 | if user.role != role:
47 | raise exceptions.AuthenticationFailed('Invalid token payload.')
48 |
49 | return (user, token)
50 |
--------------------------------------------------------------------------------
/on-orbit-frontend/tailwind.config.ts:
--------------------------------------------------------------------------------
1 | import type { Config } from "tailwindcss";
2 |
3 | const config: Config = {
4 | darkMode: ["class"],
5 | content: [
6 | "./src/pages/**/*.{js,ts,jsx,tsx,mdx}",
7 | "./src/components/**/*.{js,ts,jsx,tsx,mdx}",
8 | "./src/app/**/*.{js,ts,jsx,tsx,mdx}",
9 | ],
10 | theme: {
11 | extend: {
12 | colors: {
13 | background: 'hsl(var(--background))',
14 | foreground: 'hsl(var(--foreground))',
15 | card: {
16 | DEFAULT: 'hsl(var(--card))',
17 | foreground: 'hsl(var(--card-foreground))'
18 | },
19 | popover: {
20 | DEFAULT: 'hsl(var(--popover))',
21 | foreground: 'hsl(var(--popover-foreground))'
22 | },
23 | primary: {
24 | DEFAULT: 'hsl(var(--primary))',
25 | foreground: 'hsl(var(--primary-foreground))'
26 | },
27 | secondary: {
28 | DEFAULT: 'hsl(var(--secondary))',
29 | foreground: 'hsl(var(--secondary-foreground))'
30 | },
31 | muted: {
32 | DEFAULT: 'hsl(var(--muted))',
33 | foreground: 'hsl(var(--muted-foreground))'
34 | },
35 | accent: {
36 | DEFAULT: 'hsl(var(--accent))',
37 | foreground: 'hsl(var(--accent-foreground))'
38 | },
39 | destructive: {
40 | DEFAULT: 'hsl(var(--destructive))',
41 | foreground: 'hsl(var(--destructive-foreground))'
42 | },
43 | border: 'hsl(var(--border))',
44 | input: 'hsl(var(--input))',
45 | ring: 'hsl(var(--ring))',
46 | chart: {
47 | '1': 'hsl(var(--chart-1))',
48 | '2': 'hsl(var(--chart-2))',
49 | '3': 'hsl(var(--chart-3))',
50 | '4': 'hsl(var(--chart-4))',
51 | '5': 'hsl(var(--chart-5))'
52 | }
53 | },
54 | borderRadius: {
55 | lg: 'var(--radius)',
56 | md: 'calc(var(--radius) - 2px)',
57 | sm: 'calc(var(--radius) - 4px)'
58 | }
59 | }
60 | },
61 | plugins: [require("daisyui"), require("tailwindcss-animate")],
62 | };
63 | export default config;
64 |
--------------------------------------------------------------------------------
/Orbit_Predictor-BackEnd/api/matlab/Utils/UnitTests/Product3x3_UnitTest.m:
--------------------------------------------------------------------------------
1 | classdef (SharedTestFixtures = ...
2 | {matlab.unittest.fixtures.PathFixture('..')}) ...
3 | Product3x3_UnitTest < matlab.unittest.TestCase
4 | % Product3x3_UnitTest - Unit test for Product3x3
5 | %
6 | % =========================================================================
7 | %
8 | % Copyright (c) 2023 United States Government as represented by the
9 | % Administrator of the National Aeronautics and Space Administration.
10 | % All Rights Reserved.
11 | %
12 | % =========================================================================
13 | %
14 | % Initial version: Jul 2023; Latest update: Jul 2023
15 | %
16 | % ----------------- BEGIN CODE -----------------
17 |
18 | methods (Test)
19 | % Tests 1000 random 3x3 matrices
20 | function testDiagonal (testCase)
21 | A = 2000 * rand(1000, 9) - 1000;
22 | B = 2000 * rand(1000, 9) - 1000;
23 |
24 | P = Product3x3(A, B);
25 |
26 | for i = 1:1000
27 | M = reshape(P(i, :), 3, 3)';
28 | T = reshape(A(i, :), 3, 3)' * reshape(B(i, :), 3, 3)';
29 |
30 | testCase.verifyEqual(M, T);
31 | end
32 | end
33 | end
34 |
35 | end
36 |
37 | % ----------------- END OF CODE -----------------
38 | %
39 | % Please record any changes to the software in the change history
40 | % shown below:
41 | %
42 | %---------------- CHANGE HISTORY ------------------
43 | % Developer | Date | Description
44 | %--------------------------------------------------
45 | % E. White | 07-03-2023 | Initial development
46 |
47 | % =========================================================================
48 | %
49 | % Copyright (c) 2023 United States Government as represented by the
50 | % Administrator of the National Aeronautics and Space Administration.
51 | % All Rights Reserved.
52 | %
53 | % =========================================================================
54 |
--------------------------------------------------------------------------------
/Orbit_Predictor-BackEnd/api/matlab/Utils/EquinoctialMatrices.m:
--------------------------------------------------------------------------------
1 | function [X,P,E,J,K,Q,QRemStat,QRaw,QRem,CRem] = EquinoctialMatrices(r,v,C,RemEqCov)
2 |
3 | % Calculate equinoctial state, covariance and Jacobian matrices
4 |
5 | % Cartestian state vector and covariance matrix
6 | X = [r; v]/1e3; % Pos/vel state vec. in km units
7 | P = C/1e6; % Pos/vel covariance in km units
8 |
9 | % Calculate the equinoctial state vector and associated matrices
10 | try
11 | % Calculate equinoctial elements from the pos-vel state
12 | [~,n,af,ag,chi,psi,lM,~] = ...
13 | convert_cartesian_to_equinoctial(X(1:3),X(4:6));
14 | lM = mod(lM,2*pi);
15 | % Equinoctial state vector, E
16 | E = [n,af,ag,chi,psi,lM]';
17 | % Jacobian matrix, J = dX/dE
18 | J = jacobian_equinoctial_to_cartesian(E,X);
19 | % Inverse of Jacobian matrix, K = dE/dX
20 | K = J \ eye(6,6);
21 | % Equinoctial state covariance, Q = K * P * K'
22 | Q = cov_make_symmetric(K * P * K');
23 | % Save raw covariance
24 | QRaw = Q;
25 | % Remediate eq. covariance, if required
26 | if RemEqCov
27 | % Calc Q remediation status and remediated Q matrix
28 | [~,~,~,~,QRemStat,~,~,QRem] = ...
29 | CovRemEigValClip(Q);
30 | if QRemStat
31 | Q = cov_make_symmetric(QRem);
32 | P = cov_make_symmetric(J * Q * J');
33 | CRem = 1e6*P;
34 | else
35 | CRem = C;
36 | end
37 | else
38 | % Calc Q remediation status only
39 | [~,~,~,~,QRemStat] = CovRemEigValClip(Q);
40 | QRem = Q;
41 | CRem = C;
42 | end
43 | catch
44 | warning('EquinoctialMatrices:CalculationFailure','Equinoctial state/covariance/Jacobian calculation failure');
45 | E = NaN(6,1);
46 | J = NaN(6,6);
47 | K = NaN(6,6);
48 | Q = NaN(6,6);
49 | QRemStat = NaN(1,1);
50 | QRaw = NaN(6,6);
51 | QRem = NaN(6,6);
52 | if RemEqCov
53 | CRem = NaN(6,6);
54 | else
55 | CRem = C;
56 | end
57 | end
58 |
59 | return
60 | end
61 |
--------------------------------------------------------------------------------
/Orbit_Predictor-BackEnd/api/urls.py:
--------------------------------------------------------------------------------
1 | from django.urls import include, path
2 | from rest_framework.routers import DefaultRouter
3 | from .views import (
4 | CollisionListCreateView, CollisionDetailView, UserViewSet,
5 | ProbabilityCalcListCreateView, ProbabilityCalcDetailView,
6 | CDMSerializerListCreateView, CDMCalcDetailView, RegisterView, LoginView, CDMViewSet, RefreshTokenView, CDMCreateView, OrganizationViewSet,
7 | CollisionTradespaceView, CollisionLinearTradespaceView, CurrentUserView, CDMPrivacyToggleView, UserNotificationToggleView
8 | )
9 |
10 | router = DefaultRouter()
11 | router.register(r'cdms', CDMViewSet, basename='cdm')
12 | router.register(r'organizations', OrganizationViewSet, basename='organization')
13 | router.register(r'users', UserViewSet, basename='user')
14 |
15 | urlpatterns = [
16 | path('collisions/', CollisionListCreateView.as_view(), name='collision-list-create'),
17 | path('collisions/
/', CollisionDetailView.as_view(), name='collision-detail'),
18 | path('probabilities/', ProbabilityCalcListCreateView.as_view(), name='probability-list-create'),
19 | path('probabilities//', ProbabilityCalcDetailView.as_view(), name='probability-detail'),
20 | # path('cdms/', CDMSerializerListCreateView.as_view(), name='cdm-list-create'),
21 | path('cdms//', CDMCalcDetailView.as_view(), name='cdm-detail'),
22 | path('cdms/create/', CDMCreateView.as_view(), name='cdm-create'),
23 | path('cdms//privacy/', CDMPrivacyToggleView.as_view(), name='cdm-privacy-toggle'),
24 | path('tradespace/', CollisionTradespaceView.as_view(), name='collision-tradespace'),
25 | path('tradespace/linear/', CollisionLinearTradespaceView.as_view(), name='collision-linear-tradespace'),
26 | path('register/', RegisterView.as_view(), name='register'),
27 | path('login/', LoginView.as_view(), name='login'),
28 | path('refresh/', RefreshTokenView.as_view(), name='refresh_token'),
29 | path('users/current_user/', CurrentUserView.as_view(), name='current_user'),
30 | path('users/notifications/', UserNotificationToggleView.as_view(), name='user-notification-toggle'),
31 | path('', include(router.urls)),
32 | ]
33 |
--------------------------------------------------------------------------------
/Orbit_Predictor-BackEnd/api/matlab/Utils/PeakOverlapMD2.m:
--------------------------------------------------------------------------------
1 | function [MD2,Xu,Ps,Asdet,Asinv,POPconverged] = PeakOverlapMD2(t,t10,Eb10,Qb10,t20,Eb20,Qb20,HBR,EMD2,params)
2 |
3 | % Calculate the effective Maha. distance using the peak overlap point.
4 | %
5 | % Actually, calculate MD^2 + log(det(As)), to account for slight variations
6 | % in the determinant of the As matrix, which shows up in the denominator.
7 |
8 | % Find the peak overlap position (POP) at the nominal TCA
9 |
10 | % Calc pos/vel mean states and associated Jacobians at specified time
11 | [Jb1,xb1] = jacobian_E0_to_Xt(t,Eb10);
12 | [Jb2,xb2] = jacobian_E0_to_Xt(t,Eb20);
13 |
14 | % Calc peak overlap position
15 | % params.verbose = true;
16 | [POPconverged,~,~,~,POP] = PeakOverlapPos(t, ...
17 | xb1,Jb1,t10,Eb10,Qb10, ...
18 | xb2,Jb2,t20,Eb20,Qb20, ...
19 | HBR,params);
20 |
21 | if POPconverged > 0
22 |
23 | % Relative POP-corrected distance
24 | Xu = POP.xu2-POP.xu1; ru = Xu(1:3);
25 |
26 | % Construct the POP-corrected joint pos/vel covariance
27 | if params.XCprocessing
28 | % DCP 6x1 sensitivity vectors for cartesian pos/vel
29 | % state at the current eph time - see Hall (2021) eq 57
30 | GCp = POP.Js1*params.GEp;
31 | GCs = POP.Js2*params.GEs;
32 | Ps = POP.Ps1 + POP.Ps2 - params.sigpXsigs * (GCs*GCp'+GCp*GCs');
33 | else
34 | Ps = POP.Ps1 + POP.Ps2;
35 | end
36 |
37 | % Calculate the inverse of A, remediating NPDs with eigenvalue clipping
38 | As = Ps(1:3,1:3);
39 | % Lclip = (params.Fclip*HBR)^2;
40 | Lclip = 0; % Required for one Obj 48901 DCP case
41 | % Slower version
42 | % [~,~,~,~,~,Asdet,Asinv] = CovRemEigValClip(As,Lclip);
43 | % Faster version
44 | [Veig,Leig] = eig(As); Leig = diag(Leig); Leig(Leig < Lclip) = Lclip;
45 | Asdet = Leig(1) * Leig(2) * Leig(3);
46 | Asinv = Veig * diag(1./Leig) * Veig';
47 |
48 | % Calculate effective Maha distance
49 | MD2 = ru' * Asinv * ru;
50 |
51 | % Calculate the "effective" or modified Maha distance,
52 | % by adding ln(|As|), or ln(|As|) - 2*ln(|vu])
53 | if EMD2 == 1
54 | MD2 = MD2 + log(Asdet);
55 | elseif EMD2 == 2
56 | MD2 = MD2 + log(Asdet) - 2*log(norm(Xu(4:6)));
57 | end
58 |
59 | else
60 |
61 | % Return null values for no POP convergence
62 | MD2 = NaN; Xu = []; Ps = []; Asdet = []; Asinv = [];
63 |
64 | end
65 |
66 | return
67 | end
--------------------------------------------------------------------------------
/Orbit_Predictor-BackEnd/api/matlab/Utils/Product3x3.m:
--------------------------------------------------------------------------------
1 | function [Out] = Product3x3(a,b)
2 | % Product3x3 - vectorized 3x3 matrix multiplication routine
3 | %
4 | % Syntax: [Out] = Product3x3(a,b);
5 | %
6 | % =========================================================================
7 | %
8 | % Copyright (c) 2023 United States Government as represented by the
9 | % Administrator of the National Aeronautics and Space Administration.
10 | % All Rights Reserved.
11 | %
12 | % =========================================================================
13 | %
14 | % Input:
15 | %
16 | % a - Left matrix [nx9]
17 | % b - Right matrix [nx9]
18 | %
19 | % =========================================================================
20 | %
21 | % Output:
22 | %
23 | % Out - Matrix product [nx9]
24 | %
25 | % =========================================================================
26 | %
27 | % Notes:
28 | %
29 | % Inputs and outputs are formatted as [col1, col2, col3]:
30 | % [(1,1) (2,1) (3,1) (1,2) (2,2) (3,2) (1,3) (2,3) (3,3)]
31 | %
32 | % =========================================================================
33 | %
34 | % Initial version: May 2022; Latest update: Jul 2023
35 | %
36 | % ----------------- BEGIN CODE -----------------
37 |
38 | Out = [a(:,1).*b(:,1)+a(:,2).*b(:,4)+a(:,3).*b(:,7) ...
39 | a(:,1).*b(:,2)+a(:,2).*b(:,5)+a(:,3).*b(:,8) ...
40 | a(:,1).*b(:,3)+a(:,2).*b(:,6)+a(:,3).*b(:,9) ...
41 | a(:,4).*b(:,1)+a(:,5).*b(:,4)+a(:,6).*b(:,7) ...
42 | a(:,4).*b(:,2)+a(:,5).*b(:,5)+a(:,6).*b(:,8) ...
43 | a(:,4).*b(:,3)+a(:,5).*b(:,6)+a(:,6).*b(:,9) ...
44 | a(:,7).*b(:,1)+a(:,8).*b(:,4)+a(:,9).*b(:,7) ...
45 | a(:,7).*b(:,2)+a(:,8).*b(:,5)+a(:,9).*b(:,8) ...
46 | a(:,7).*b(:,3)+a(:,8).*b(:,6)+a(:,9).*b(:,9)];
47 | end
48 |
49 | % ----------------- END OF CODE -----------------
50 | %
51 | % Please record any changes to the software in the change history
52 | % shown below:
53 | %
54 | % ----------------- CHANGE HISTORY ------------------
55 | % Developer | Date | Description
56 | % ---------------------------------------------------
57 | % L. Baars | 05-19-2022 | Initial development
58 | % E. White | 07-12-2023 | Added compliant documentation
59 |
60 | % =========================================================================
61 | %
62 | % Copyright (c) 2023 United States Government as represented by the
63 | % Administrator of the National Aeronautics and Space Administration.
64 | % All Rights Reserved.
65 | %
66 | % =========================================================================
67 |
--------------------------------------------------------------------------------
/Orbit_Predictor-BackEnd/api/models/collision.py:
--------------------------------------------------------------------------------
1 | from pathlib import Path
2 | from django.db import models
3 | from .cdm import CDM
4 | import matlab.engine
5 |
6 | class Collision(models.Model):
7 | cdm = models.ForeignKey(CDM, on_delete=models.CASCADE, related_name='collisions')
8 | probability_of_collision = models.FloatField()
9 | sat1_object_designator = models.CharField(max_length=50)
10 | sat2_object_designator = models.CharField(max_length=50)
11 |
12 | @classmethod
13 | def create_from_cdm(cls, cdm):
14 | if not cdm:
15 | raise ValueError("A valid CDM object must be provided.")
16 | matlabPathFile = Path(__file__).resolve().parent.parent / 'matlab'
17 |
18 | eng = matlab.engine.start_matlab()
19 | eng.addpath(str(matlabPathFile))
20 | r1 = matlab.double([cdm.sat1_x, cdm.sat1_y, cdm.sat1_z])
21 | v1 = matlab.double([cdm.sat1_x_dot, cdm.sat1_y_dot, cdm.sat1_z_dot])
22 | cov1 = matlab.double([
23 | [cdm.sat1_cov_rr, cdm.sat1_cov_rt, cdm.sat1_cov_rn],
24 | [cdm.sat1_cov_tr, cdm.sat1_cov_tt, cdm.sat1_cov_tn],
25 | [cdm.sat1_cov_nr, cdm.sat1_cov_nt, cdm.sat1_cov_nn]
26 | ])
27 |
28 | r2 = matlab.double([cdm.sat2_x, cdm.sat2_y, cdm.sat2_z])
29 | v2 = matlab.double([cdm.sat2_x_dot, cdm.sat2_y_dot, cdm.sat2_z_dot])
30 | cov2 = matlab.double([
31 | [cdm.sat2_cov_rr, cdm.sat2_cov_rt, cdm.sat2_cov_rn],
32 | [cdm.sat2_cov_tr, cdm.sat2_cov_tt, cdm.sat2_cov_tn],
33 | [cdm.sat2_cov_nr, cdm.sat2_cov_nt, cdm.sat2_cov_nn]
34 | ])
35 | HBR = cdm.hard_body_radius
36 | RelTol = 1e-08
37 | HBRType = 'circle'
38 | probability_of_collision_matlab = eng.Pc2D_Foster(r1, v1, cov1, r2, v2, cov2, HBR, RelTol, HBRType, nargout=1)
39 | probability_of_collision = float(probability_of_collision_matlab)
40 |
41 | if probability_of_collision > 1.0:
42 | probability_of_collision = 1.0
43 |
44 | return cls.objects.create(
45 | cdm=cdm,
46 | probability_of_collision=probability_of_collision,
47 | sat1_object_designator=cdm.sat1_object_designator,
48 | sat2_object_designator=cdm.sat2_object_designator,
49 | )
50 |
51 | def save(self, *args, **kwargs):
52 | # Ensure satellite IDs are copied from the related CDM
53 | if self.cdm:
54 | self.sat1_object_designator = self.cdm.sat1_object_designator
55 | self.sat2_object_designator = self.cdm.sat2_object_designator
56 | super().save(*args, **kwargs)
57 |
58 | def __str__(self):
59 | return f"Collision {self.id} between {self.sat1_object_designator} and {self.sat2_object_designator}"
60 |
--------------------------------------------------------------------------------
/Orbit_Predictor-BackEnd/api/migrations/0004_remove_user_password_hash_cdm_privacy_user_groups_and_more.py:
--------------------------------------------------------------------------------
1 | # Generated by Django 5.1.3 on 2024-12-07 02:05
2 |
3 | from django.db import migrations, models
4 |
5 |
6 | class Migration(migrations.Migration):
7 |
8 | dependencies = [
9 | ('api', '0003_remove_collision_conjunction_id_and_more'),
10 | ('auth', '0012_alter_user_first_name_max_length'),
11 | ]
12 |
13 | operations = [
14 | migrations.RemoveField(
15 | model_name='user',
16 | name='password_hash',
17 | ),
18 | migrations.AddField(
19 | model_name='cdm',
20 | name='privacy',
21 | field=models.BooleanField(default=True),
22 | ),
23 | migrations.AddField(
24 | model_name='user',
25 | name='groups',
26 | field=models.ManyToManyField(blank=True, help_text='The groups this user belongs to. A user will get all permissions granted to each of their groups.', related_name='user_set', related_query_name='user', to='auth.group', verbose_name='groups'),
27 | ),
28 | migrations.AddField(
29 | model_name='user',
30 | name='is_active',
31 | field=models.BooleanField(default=True),
32 | ),
33 | migrations.AddField(
34 | model_name='user',
35 | name='is_staff',
36 | field=models.BooleanField(default=False),
37 | ),
38 | migrations.AddField(
39 | model_name='user',
40 | name='is_superuser',
41 | field=models.BooleanField(default=False, help_text='Designates that this user has all permissions without explicitly assigning them.', verbose_name='superuser status'),
42 | ),
43 | migrations.AddField(
44 | model_name='user',
45 | name='last_login',
46 | field=models.DateTimeField(blank=True, null=True, verbose_name='last login'),
47 | ),
48 | migrations.AddField(
49 | model_name='user',
50 | name='password',
51 | field=models.CharField(default='password', max_length=128, verbose_name='password'),
52 | preserve_default=False,
53 | ),
54 | migrations.AddField(
55 | model_name='user',
56 | name='role',
57 | field=models.CharField(choices=[('admin', 'Administrator'), ('collision_analyst', 'Collision Analyst'), ('user', 'User')], default='user', max_length=50),
58 | ),
59 | migrations.AddField(
60 | model_name='user',
61 | name='user_permissions',
62 | field=models.ManyToManyField(blank=True, help_text='Specific permissions for this user.', related_name='user_set', related_query_name='user', to='auth.permission', verbose_name='user permissions'),
63 | ),
64 | ]
65 |
--------------------------------------------------------------------------------
/Orbit_Predictor-BackEnd/api/matlab/Utils/GetAlfanoTestCases.m:
--------------------------------------------------------------------------------
1 | function conj = GetAlfanoTestCases(params)
2 | %==========================================================================
3 | %
4 | % Use the Alfano (2009) test case Excel files to a conjunction list in
5 | % the form output by function OCMDB_get_conjunctions.
6 | %
7 | % REFERENCE:
8 | %
9 | % S.Alfano, "Satellite Conjunction Monte Carlo Analysis" AAS 09-233 (2009).
10 | %
11 | %==========================================================================
12 |
13 | % Initializations and defaults
14 |
15 | if (nargin == 0); params = []; end
16 |
17 | if ~isfield(params,'verbose'); params.verbose = []; end
18 | if isempty(params.verbose); params.verbose = 0; end
19 |
20 | if ~isfield(params,'case_list'); params.case_list = []; end
21 | if isempty(params.case_list); params.case_list = (1:12); end
22 |
23 | if ~isfield(params,'data_path'); params.data_path = []; end
24 | if isempty(params.data_path); params.data_path = ...
25 | 'Alfano_2009_Test_Cases'; end
26 |
27 | % Number of A09 conjunctions to retrieve
28 |
29 | Nconj = numel(params.case_list);
30 |
31 | % Load the database. This creates a structure called "DB".
32 |
33 | if params.verbose
34 | disp(' ');
35 | disp(['Loading ' num2str(Nconj) ' Alfano (2009) conjunction test cases from path:']);
36 | disp([' ' params.data_path]);
37 | end
38 |
39 | % Allocate the conjunction arrays
40 | conj.case = NaN(1,Nconj); % Test case number
41 | conj.X1 = NaN(6,Nconj); % ECI state
42 | conj.C1 = NaN(6,6,Nconj); % ECI covariance
43 | conj.X2 = NaN(6,Nconj); % ECI state
44 | conj.C2 = NaN(6,6,Nconj); % ECI covariance
45 | conj.HBR = NaN(1,Nconj); % Combined hard-body radii
46 |
47 | % Loop through conjunctions
48 |
49 | for nc = 1:Nconj
50 |
51 | % Get this A09 test case
52 | conj.case(nc) = params.case_list(nc);
53 | if params.verbose
54 | disp([' Loading case number ' num2str(conj.case(nc)) ' from Excel files']);
55 | end
56 | tc = get_alfano_test_case(conj.case(nc),params.data_path);
57 |
58 | % Combined HBR
59 | conj.HBR(nc) = tc.HBR;
60 |
61 | % Primary object ECI state
62 | X1 = [tc.R1o' tc.V1o'];
63 |
64 | % Primary object ECI covariance
65 | C1ECI = tc.P1o;
66 | % C1ECI = cov_make_symmetric(C1ECI);
67 |
68 | % Secondary object ECI state
69 | X2 = [tc.R2o' tc.V2o'];
70 |
71 | % Secondary object ECI covariance
72 | C2ECI = tc.P2o;
73 | % C2ECI = cov_make_symmetric(C2ECI);
74 |
75 | % Define the outputs in MKS units
76 | conj.X1(:,nc) = X1;
77 | conj.X2(:,nc) = X2;
78 | conj.C1(:,:,nc) = C1ECI;
79 | conj.C2(:,:,nc) = C2ECI;
80 |
81 | end
82 |
83 | return;
84 | end
85 |
--------------------------------------------------------------------------------
/Orbit_Predictor-BackEnd/api/models/user.py:
--------------------------------------------------------------------------------
1 | # api/models/user.py
2 | import bcrypt
3 | from django.db import models
4 | import uuid
5 | from django.contrib.auth.models import AbstractBaseUser, PermissionsMixin, BaseUserManager
6 |
7 |
8 | class UserManager(BaseUserManager):
9 | def create_user(self, email, password=None, role='user', **extra_fields):
10 | """
11 | Creates and saves a User with the given email, password, and role.
12 | """
13 | if not email:
14 | raise ValueError('The Email field must be set')
15 | email = self.normalize_email(email)
16 | user = self.model(email=email, role=role, **extra_fields)
17 | if password:
18 | user.set_password(password)
19 | else:
20 | raise ValueError('The Password field must be set')
21 | user.save(using=self._db)
22 | return user
23 |
24 | def create_superuser(self, email, password=None, role='admin', **extra_fields):
25 | """
26 | Creates and saves a superuser with the given email, password, and role.
27 | """
28 | extra_fields.setdefault('is_staff', True)
29 | extra_fields.setdefault('is_superuser', True)
30 |
31 | if role != 'admin':
32 | raise ValueError('Superuser must have role=admin.')
33 | if not extra_fields.get('is_staff'):
34 | raise ValueError('Superuser must have is_staff=True.')
35 | if not extra_fields.get('is_superuser'):
36 | raise ValueError('Superuser must have is_superuser=True.')
37 |
38 | return self.create_user(email, password, role, **extra_fields)
39 |
40 |
41 | class User(AbstractBaseUser, PermissionsMixin):
42 | ROLE_CHOICES = (
43 | ('admin', 'Administrator'),
44 | ('collision_analyst', 'Collision Analyst'),
45 | ('user', 'User'),
46 | )
47 |
48 | id = models.UUIDField(primary_key=True, default=uuid.uuid4, editable=False)
49 | email = models.EmailField(unique=True)
50 | role = models.CharField(max_length=50, choices=ROLE_CHOICES, default='user')
51 | is_active = models.BooleanField(default=True)
52 | is_staff = models.BooleanField(default=False) # Required for admin access
53 | created_at = models.DateTimeField(auto_now_add=True)
54 |
55 | notifications = models.BooleanField(default=True)
56 |
57 | interested_cdms = models.ManyToManyField(
58 | 'CDM',
59 | related_name='interested_users',
60 | blank=True
61 | )
62 |
63 | objects = UserManager()
64 |
65 | USERNAME_FIELD = 'email'
66 | REQUIRED_FIELDS = []
67 |
68 | def __str__(self):
69 | return self.email
70 |
71 | # Override set_password and check_password to use bcrypt
72 | def set_password(self, raw_password):
73 | salt = bcrypt.gensalt()
74 | hashed = bcrypt.hashpw(raw_password.encode('utf-8'), salt)
75 | self.password = hashed.decode('utf-8')
76 |
77 | def check_password(self, raw_password):
78 | return bcrypt.checkpw(raw_password.encode('utf-8'), self.password.encode('utf-8'))
79 |
--------------------------------------------------------------------------------
/Orbit_Predictor-BackEnd/api/matlab/Utils/UnitTests/CheckAndResizeCov_UnitTest.m:
--------------------------------------------------------------------------------
1 | classdef (SharedTestFixtures = { ...
2 | matlab.unittest.fixtures.PathFixture('..')}) ...
3 | CheckAndResizeCov_UnitTest < matlab.unittest.TestCase
4 |
5 | properties (TestParameter)
6 | numR = {1; 1};
7 | cov = {[1 2 3; 4 5 6; 7 8 9];
8 | [9 8 7; 6 5 4; 3 2 1]};
9 | expCov = {[1 2 3 4 5 6 7 8 9];
10 | [9 8 7 6 5 4 3 2 1]};
11 | end
12 |
13 | methods (Test)
14 | function test_nx9_valid_input(testCase)
15 | % Tests if an nx9 with the right number of rows is successfully
16 | % accepted
17 | numR = 2;
18 | cov = [1 2 3 4 5 6 7 8 9;
19 | 2 3 4 5 6 7 8 9 10];
20 | expCov = cov;
21 |
22 | actCov = CheckAndResizeCov(numR, cov);
23 |
24 | testCase.verifyEqual(actCov, expCov);
25 | end
26 |
27 | function test_nx9_num_rows_mismatch(testCase)
28 | % Tests if the function returns an error if there is a mismatch
29 | % between the "n" in the nx9 and number of rows
30 | numR = 1;
31 | cov = [1 2 3 4 5 6 7 8 9;
32 | 2 3 4 5 6 7 8 9 10];
33 |
34 | verifyError(testCase,@() CheckAndResizeCov(numR, cov),'CheckAndResizeCov:RowCountMismatch2D');
35 | end
36 |
37 | function test_3x3(testCase)
38 | % Tests if the function properly reformats a 3x3 into a 1x9
39 | numR = {1,2};
40 | cov = {[1 2 3;
41 | 4 5 6;
42 | 7 8 9], ...
43 | [9 8 7;
44 | 6 5 4;
45 | 3 2 1]};
46 | expCov = {[1 2 3 4 5 6 7 8 9], ...
47 | [9 8 7 6 5 4 3 2 1;
48 | 9 8 7 6 5 4 3 2 1]};
49 |
50 | for i = 1:length(numR)
51 | actCov = CheckAndResizeCov(numR{i}, cov{i});
52 |
53 | testCase.verifyEqual(actCov, expCov{i});
54 | end
55 | end
56 |
57 | function test_6x6(testCase)
58 | % Test if the function properly reformats a 6x6 into a 1x9
59 | numR = {1,2};
60 | cov = {[1 2 3 4 5 6;
61 | 2 7 8 9 0 1;
62 | 3 8 2 3 4 5;
63 | 4 9 3 6 7 8;
64 | 5 0 4 7 9 0;
65 | 6 1 5 8 0 1], ...
66 | [9 8 7 6 5 4;
67 | 8 3 2 1 0 9;
68 | 7 2 8 7 6 5;
69 | 6 1 7 4 3 2;
70 | 5 0 6 3 1 0;
71 | 4 9 5 2 0 9]};
72 | expCov = {[1 2 3 2 7 8 3 8 2], ...
73 | [9 8 7 8 3 2 7 2 8;
74 | 9 8 7 8 3 2 7 2 8]};
75 |
76 | for i = 1:length(numR)
77 | actCov = CheckAndResizeCov(numR{i}, cov{i});
78 |
79 | testCase.verifyEqual(actCov, expCov{i});
80 | end
81 | end
82 | end
83 | end
--------------------------------------------------------------------------------
/Orbit_Predictor-BackEnd/api/migrations/0006_cdm_sat1_catalog_name_cdm_sat1_covariance_method_and_more.py:
--------------------------------------------------------------------------------
1 | # Generated by Django 5.1.3 on 2025-02-03 02:45
2 |
3 | from django.db import migrations, models
4 |
5 |
6 | class Migration(migrations.Migration):
7 |
8 | dependencies = [
9 | ('api', '0005_merge_20250203_0240'),
10 | ]
11 |
12 | operations = [
13 | migrations.AddField(
14 | model_name='cdm',
15 | name='sat1_catalog_name',
16 | field=models.CharField(blank=True, max_length=100, null=True),
17 | ),
18 | migrations.AddField(
19 | model_name='cdm',
20 | name='sat1_covariance_method',
21 | field=models.CharField(blank=True, max_length=100, null=True),
22 | ),
23 | migrations.AddField(
24 | model_name='cdm',
25 | name='sat1_international_designator',
26 | field=models.CharField(blank=True, max_length=100, null=True),
27 | ),
28 | migrations.AddField(
29 | model_name='cdm',
30 | name='sat1_object_name',
31 | field=models.CharField(blank=True, max_length=100, null=True),
32 | ),
33 | migrations.AddField(
34 | model_name='cdm',
35 | name='sat1_object_type',
36 | field=models.CharField(blank=True, max_length=100, null=True),
37 | ),
38 | migrations.AddField(
39 | model_name='cdm',
40 | name='sat1_operator_organization',
41 | field=models.CharField(blank=True, max_length=100, null=True),
42 | ),
43 | migrations.AddField(
44 | model_name='cdm',
45 | name='sat1_reference_frame',
46 | field=models.CharField(blank=True, max_length=100, null=True),
47 | ),
48 | migrations.AddField(
49 | model_name='cdm',
50 | name='sat2_catalog_name',
51 | field=models.CharField(blank=True, max_length=100, null=True),
52 | ),
53 | migrations.AddField(
54 | model_name='cdm',
55 | name='sat2_covariance_method',
56 | field=models.CharField(blank=True, max_length=100, null=True),
57 | ),
58 | migrations.AddField(
59 | model_name='cdm',
60 | name='sat2_international_designator',
61 | field=models.CharField(blank=True, max_length=100, null=True),
62 | ),
63 | migrations.AddField(
64 | model_name='cdm',
65 | name='sat2_object_name',
66 | field=models.CharField(blank=True, max_length=100, null=True),
67 | ),
68 | migrations.AddField(
69 | model_name='cdm',
70 | name='sat2_object_type',
71 | field=models.CharField(blank=True, max_length=100, null=True),
72 | ),
73 | migrations.AddField(
74 | model_name='cdm',
75 | name='sat2_operator_organization',
76 | field=models.CharField(blank=True, max_length=100, null=True),
77 | ),
78 | migrations.AddField(
79 | model_name='cdm',
80 | name='sat2_reference_frame',
81 | field=models.CharField(blank=True, max_length=100, null=True),
82 | ),
83 | ]
84 |
--------------------------------------------------------------------------------
/Orbit_Predictor-BackEnd/api/matlab/Utils/CheckAndResizeCov.m:
--------------------------------------------------------------------------------
1 | function [cov] = CheckAndResizeCov(numR, cov)
2 | % CheckAndResizeCov - Resizes the covariance passed in into an n x 9
3 | % matrix to be used in vectorized covariance
4 | % processing.
5 | %
6 | % Description:
7 | %
8 | % Reformats the input covariance into an n x 9 matrix representing the
9 | % 3x3 position covariance in the following format:
10 | % cov = [C(1,1) C(1,2) C(1,3) C(2,1) C(2,2) C(2,3) C(3,1) C(3,2) C(3,3)]
11 | % Where "n" is the number of rows passed in.
12 | %
13 | % Depending on the input format, this function will respond in several
14 | % different ways to reformat the cov matrix:
15 | % - If the input is an nx9, then this function verifies that n matches
16 | % numR.
17 | % - If the input is a 3x3, then this function will reformat the matrix
18 | % into a 1x9 and then repeat this vector numR times to create n rows.
19 | % - If the input is a 6x6, then this function will take the upper left
20 | % 3x3 component and reformat it into a 1x9. Then it will repeat this
21 | % vector numR times to create n rows.
22 | % - If the input is a 3x3xn, then this function will verify that n
23 | % matches numR and will reformat each 3x3 into a corresponding 1x9.
24 | % - If the input is a 6x6xn, then this function will verify that n
25 | % matches and will take the upper left 3x3 component of each 6x6 into a
26 | % corresponding 1x9.
27 | % - Any other input is considered an error.
28 |
29 | covSize = size(cov);
30 | % If a 2D covariance matrix was passed in
31 | if size(covSize,2) == 2
32 | if covSize(2) ~= 9 && (covSize(1) ~= 3 || covSize(2) ~= 3) && ...
33 | (covSize(1) ~= 6 || covSize(2) ~= 6)
34 | error('CheckAndResizeCov:Invalid2DCov','2D Covariance matrix must have 9 columns or be a 3x3 or 6x6 matrix!');
35 | end
36 | % Resize down to a 3x3 if a 6x6 was passed in
37 | if covSize(1) == 6 && covSize(2) == 6
38 | cov = cov(1:3,1:3);
39 | covSize = size(cov);
40 | end
41 | if covSize(1) == 1
42 | cov = repmat(cov,numR,1);
43 | elseif covSize(1) == 3 && covSize(2) == 3
44 | cov = reshape(permute(cov,[2 1]),1,9);
45 | cov = repmat(cov,numR,1);
46 | elseif covSize(1) ~= numR
47 | error('CheckAndResizeCov:RowCountMismatch2D','2D Covariance cannot be resized to match r matrix');
48 | end
49 | % If a 3D covariance matrix was passed in
50 | elseif size(covSize,2) == 3
51 | % The 3D matrix should have the dimension 3x3xnumR or 6x6xnumR
52 | if (covSize(1) ~= 3 || covSize(2) ~= 3 || covSize(3) ~= numR) && ...
53 | (covSize(1) ~= 6 || covSize(2) ~= 6 || covSize(3) ~= numR)
54 | error('CheckAndResizeCov:Invalid3DCov','3D covariance matrix must be of size 3x3xnumR or 6x6xnumR');
55 | end
56 | % Resize down to a 3x3xnumR if a 6x6xnumR was passed in
57 | if covSize(1) == 6 && covSize(2) == 6
58 | cov = cov(1:3,1:3,:);
59 | end
60 | cov = reshape(permute(cov,[3 2 1]),numR,9);
61 | else
62 | error('CheckAndResizeCov:InvalidCov','Improperly sized covariance was detected');
63 | end
64 | end
--------------------------------------------------------------------------------
/Orbit_Predictor-BackEnd/api/matlab/Utils/jacobian_E0_to_Xt.m:
--------------------------------------------------------------------------------
1 | function [JT,XT] = jacobian_E0_to_Xt(T,E0,fr,mu,errflag)
2 | % =========================================================================
3 | %
4 | % Use epoch equinoctial orbital element state, Ecent0, to calculate the
5 | % cartesian states, XcentT = (rvec',vvec')', for a set of offset times from
6 | % epoch, T = t-t0. Here Ecent0 represents the "center-of-expansion"
7 | % equinoctial elements used for the first-order Taylor series expansion of
8 | % the motion:
9 | % X(t) = Xcent(t) + Jcent(t) * (E0 - Ecent0)
10 | % or
11 | % E0 = Ecent0 + [Jcent(t)]^-1 * (X(t)-Xcent(t))
12 | % with
13 | % J(t) = [dX(t)/dE0] at E0 = Ecent0 (i.e., X(t) = Xcent(t))
14 | %
15 | % =========================================================================
16 | %
17 | % INPUT:
18 | %
19 | % T = Time offsets from initial (s). [NTx1] or [1xNT]
20 | % E0 = Equinoctial elements at initial time t0. [6x1]
21 | % [n,af,ag,chi,psi,lam]' at t = t0.
22 | % See Vallado and Alfano (2015) for details.
23 | % fr = Equinoctial element retrograde factor (optional, default = +1)
24 | % mu = Gravitational constant (optional).
25 | % errflag = Error flag (optional, default = 2)
26 | % 0 => No error or warning issued for F nonconvergence (not recommended)
27 | % 1 => Warning issued for F nonconvergence (not recommended)
28 | % 2 => Error issued for F nonconvergence (recommended)
29 | %
30 | % =========================================================================
31 | %
32 | % OUTPUT:
33 | %
34 | % JT = EpochEquinoctial-to-EphemerisCartesian transformations [6x6xNT]
35 | % XT = Center-of-expansion cartesian state ephemeris (km & km/s) [6xNT]
36 | %
37 | % =========================================================================
38 | %
39 | % REFERENCE:
40 | %
41 | % Vallado and Alfano (2015), AAS 15-537
42 | % Broucke and Cefola (1972), Celestial Mechanics, Vol.5, pp.303-310
43 | %
44 | % =========================================================================
45 |
46 | % Defaults and intializations
47 |
48 | Nargin = nargin; % Nargout = nargout;
49 |
50 | na = 3;
51 | if Nargin < na || isempty(fr)
52 | % Default to prograde equinoctial elements
53 | fr = 1;
54 | end
55 |
56 | na = na+1;
57 | if Nargin < na || isempty(mu)
58 | % Earth gravitational constant (EGM-96) [km^3/s^2]
59 | mu = 3.986004418e5;
60 | end
61 |
62 | na = na+1;
63 | if Nargin < na || isempty(errflag)
64 | errflag = 2;
65 | end
66 |
67 | % Number of ephemeris times
68 | NT = numel(T);
69 |
70 | % Calculate X(t)
71 | [rT,vT] = convert_equinoctial_to_cartesian( ...
72 | E0(1),E0(2),E0(3),E0(4),E0(5),E0(6),T, ...
73 | fr,mu,[],[],errflag);
74 | XT = [rT; vT];
75 |
76 | % Initialize dE(t)/dE(t0) STM
77 | phi = eye(6,6);
78 |
79 | % Initialize E(t)
80 | ET = E0;
81 |
82 | % Initialize output Jacobian array
83 | JT = NaN(6,6,NT);
84 |
85 | % Loop over times and calculate output
86 |
87 | for nT=1:NT
88 | % Calculate the equinoctial mean longitude at this time
89 | ET(6) = E0(6)+T(nT)*E0(1);
90 | % Define off-diagonal dE(t)/dE(t0) STM element
91 | phi(6,1) = T(nT);
92 | % Calculate dX(t)/dE(t) STM
93 | J = jacobian_equinoctial_to_cartesian(ET,XT(:,nT),fr,mu);
94 | % Calculate the dX(t)/dE0 STM
95 | JT(:,:,nT) = J*phi;
96 | end
97 |
98 | return
99 | end
--------------------------------------------------------------------------------
/Orbit_Predictor-BackEnd/api/sample_data/oct5_data/cdm0.json:
--------------------------------------------------------------------------------
1 | [
2 | {
3 | "CCSDS_CDM_VERS": "1.0",
4 | "CREATION_DATE": "2024-10-04T03:34:25.108",
5 | "ORIGINATOR": "CSpoc",
6 | "MESSAGE_ID": "7407_conj0_1338",
7 | "TCA": "2024-10-05T17:10:36.475",
8 | "MISS_DISTANCE": "14923",
9 | "SAT1_OBJECT": "OBJECT1",
10 | "SAT1_OBJECT_DESIGNATOR": "12345",
11 | "SAT1_CATALOG_NAME": "SATCAT",
12 | "SAT1_INTERNATIONAL_DESIGNATOR": "2019-030E",
13 | "SAT1_EPHEMERIS_NAME": "NONE",
14 | "SAT1_COVARIANCE_METHOD": "CALCULATED",
15 | "SAT1_MANEUVERABLE": "NO",
16 | "SAT1_REFERENCE_FRAME": "ITRF",
17 | "SAT1_X": "2569.097065",
18 | "SAT1_Y": "2245.654904",
19 | "SAT1_Z": "6281.497978",
20 | "SAT1_X_DOT": "4.418769571",
21 | "SAT1_Y_DOT": "4.833547743",
22 | "SAT1_Z_DOT": "-3.526774282",
23 | "SAT1_CR_R": "16.466809999999988",
24 | "SAT1_CT_R": "-124109.12391596804",
25 | "SAT1_CT_T": "18938.469999999998",
26 | "SAT1_CN_R": "15.272398047558655",
27 | "SAT1_CN_T": "2452.061015626061",
28 | "SAT1_CN_N": "3.1734509999999996",
29 | "SAT1_CRDOT_R": "0.02171203999999999",
30 | "SAT1_CRDOT_T": "-411.05878366331143",
31 | "SAT1_CRDOT_N": "-0.002652613771943246",
32 | "SAT1_CRDOT_RDOT": "0.02171203999999999",
33 | "SAT1_CTDOT_R": "-0.0002704973468097328",
34 | "SAT1_CTDOT_T": "0.06805182225726791",
35 | "SAT1_CTDOT_N": "-1.681653824600242e-05",
36 | "SAT1_CTDOT_RDOT": "-7.643144107555492e-08",
37 | "SAT1_CTDOT_TDOT": "1.6743309999999926e-05",
38 | "SAT1_CNDOT_R": "2.5077447026675033e-05",
39 | "SAT1_CNDOT_T": "-0.0002619771444092925",
40 | "SAT1_CNDOT_N": "1.5047272707261375e-05",
41 | "SAT1_CNDOT_RDOT": "4.0999922806762933e-10",
42 | "SAT1_CNDOT_TDOT": "-2.6774131967338362e-11",
43 | "SAT1_CNDOT_NDOT": "8.919619999999988e-06",
44 | "SAT2_OBJECT": "OBJECT1",
45 | "SAT2_OBJECT_DESIGNATOR": "67891",
46 | "SAT2_CATALOG_NAME": "SATCAT",
47 | "SAT2_INTERNATIONAL_DESIGNATOR": "2021-030E",
48 | "SAT2_EPHEMERIS_NAME": "NONE",
49 | "SAT2_COVARIANCE_METHOD": "CALCULATED",
50 | "SAT2_MANEUVERABLE": "YES",
51 | "SAT2_REFERENCE_FRAME": "ITRF",
52 | "SAT2_X": "4684.39413",
53 | "SAT2_Y": "-1485.2901920000004",
54 | "SAT2_Z": "26229.795955999998",
55 | "SAT2_X_DOT": "16.037539142",
56 | "SAT2_Y_DOT": "12646.667095485998",
57 | "SAT2_Z_DOT": "5518.846451436",
58 | "SAT2_CR_R": "71140.79",
59 | "SAT2_CT_R": "-173447055517246.4",
60 | "SAT2_CT_T": "2955756000.0",
61 | "SAT2_CN_R": "73262356.0198991",
62 | "SAT2_CN_T": "-16541557533.518105",
63 | "SAT2_CN_N": "2172.732",
64 | "SAT2_CRDOT_R": "3395.637",
65 | "SAT2_CRDOT_T": "-10036658308451.043",
66 | "SAT2_CRDOT_N": "26304.7422813006",
67 | "SAT2_CRDOT_RDOT": "3395.637",
68 | "SAT2_CTDOT_R": "-1420.343747333301",
69 | "SAT2_CTDOT_T": "17713388.417774897",
70 | "SAT2_CTDOT_N": "-48.19052769752484",
71 | "SAT2_CTDOT_RDOT": "-20.40639119554754",
72 | "SAT2_CTDOT_TDOT": "0.02723944",
73 | "SAT2_CNDOT_R": "2654.7336674472876",
74 | "SAT2_CNDOT_T": "86797934.59951398",
75 | "SAT2_CNDOT_N": "234.35235697465833",
76 | "SAT2_CNDOT_RDOT": "-99.86165235094477",
77 | "SAT2_CNDOT_TDOT": "-0.0027195986896554718",
78 | "SAT2_CNDOT_NDOT": "0.14935299999999999"
79 | }
80 | ]
--------------------------------------------------------------------------------
/Orbit_Predictor-BackEnd/api/sample_data/oct5_data/cdm3.json:
--------------------------------------------------------------------------------
1 | [
2 | {
3 | "CCSDS_CDM_VERS": "1.0",
4 | "CREATION_DATE": "2024-10-05T03:15:53.042",
5 | "ORIGINATOR": "CSpoc",
6 | "MESSAGE_ID": "1671_conj3_9268",
7 | "TCA": "2024-10-05T17:10:36.475",
8 | "MISS_DISTANCE": "14579",
9 | "SAT1_OBJECT": "OBJECT1",
10 | "SAT1_OBJECT_DESIGNATOR": "12345",
11 | "SAT1_CATALOG_NAME": "SATCAT",
12 | "SAT1_INTERNATIONAL_DESIGNATOR": "2019-030E",
13 | "SAT1_EPHEMERIS_NAME": "NONE",
14 | "SAT1_COVARIANCE_METHOD": "CALCULATED",
15 | "SAT1_MANEUVERABLE": "NO",
16 | "SAT1_REFERENCE_FRAME": "ITRF",
17 | "SAT1_X": "2569.097065",
18 | "SAT1_Y": "2245.654904",
19 | "SAT1_Z": "6281.497978",
20 | "SAT1_X_DOT": "4.418769571",
21 | "SAT1_Y_DOT": "4.833547743",
22 | "SAT1_Z_DOT": "-3.526774282",
23 | "SAT1_CR_R": "9.460428",
24 | "SAT1_CT_R": "100.95118385368234",
25 | "SAT1_CT_T": "494.3005",
26 | "SAT1_CN_R": "13.374609096870207",
27 | "SAT1_CN_T": "259.87230439211817",
28 | "SAT1_CN_N": "3.337079",
29 | "SAT1_CRDOT_R": "0.0004868924999999991",
30 | "SAT1_CRDOT_T": "-0.23848413183125663",
31 | "SAT1_CRDOT_N": "-0.00024371647114500403",
32 | "SAT1_CRDOT_RDOT": "0.0004868924999999991",
33 | "SAT1_CTDOT_R": "-0.00010277619533014733",
34 | "SAT1_CTDOT_T": "-0.00031933176551537717",
35 | "SAT1_CTDOT_N": "-1.550394012729719e-05",
36 | "SAT1_CTDOT_RDOT": "3.824140680442102e-10",
37 | "SAT1_CTDOT_TDOT": "1.0876239999999834e-05",
38 | "SAT1_CNDOT_R": "1.837409542882393e-05",
39 | "SAT1_CNDOT_T": "0.00016486648335839135",
40 | "SAT1_CNDOT_N": "1.5561972107337534e-05",
41 | "SAT1_CNDOT_RDOT": "-2.0723211585924206e-10",
42 | "SAT1_CNDOT_TDOT": "-2.1257833532321998e-11",
43 | "SAT1_CNDOT_NDOT": "8.516542999999602e-06",
44 | "SAT2_OBJECT": "OBJECT1",
45 | "SAT2_OBJECT_DESIGNATOR": "67891",
46 | "SAT2_CATALOG_NAME": "SATCAT",
47 | "SAT2_INTERNATIONAL_DESIGNATOR": "2021-030E",
48 | "SAT2_EPHEMERIS_NAME": "NONE",
49 | "SAT2_COVARIANCE_METHOD": "CALCULATED",
50 | "SAT2_MANEUVERABLE": "YES",
51 | "SAT2_REFERENCE_FRAME": "ITRF",
52 | "SAT2_X": "4665.49413",
53 | "SAT2_Y": "-1347.5901919999997",
54 | "SAT2_Z": "25913.695956",
55 | "SAT2_X_DOT": "15.837539142",
56 | "SAT2_Y_DOT": "12646.667095485998",
57 | "SAT2_Z_DOT": "5518.846451436",
58 | "SAT2_CR_R": "68518.15000000001",
59 | "SAT2_CT_R": "-163998711741309.6",
60 | "SAT2_CT_T": "2927633000.0",
61 | "SAT2_CN_R": "71959560.65875487",
62 | "SAT2_CN_T": "-23899082735.30432",
63 | "SAT2_CN_N": "2171.335000000001",
64 | "SAT2_CRDOT_R": "3363.2929999999983",
65 | "SAT2_CRDOT_T": "-9846467682226.9",
66 | "SAT2_CRDOT_N": "34722.3962751339",
67 | "SAT2_CRDOT_RDOT": "3363.2929999999983",
68 | "SAT2_CTDOT_R": "-1346.8315032614892",
69 | "SAT2_CTDOT_T": "15849903.560941886",
70 | "SAT2_CTDOT_N": "-47.94457892416641",
71 | "SAT2_CTDOT_RDOT": "-18.264878606315833",
72 | "SAT2_CTDOT_TDOT": "0.026997380000000005",
73 | "SAT2_CNDOT_R": "2657.4141866193004",
74 | "SAT2_CNDOT_T": "85617046.08932368",
75 | "SAT2_CNDOT_N": "234.1003829657977",
76 | "SAT2_CNDOT_RDOT": "-98.47141484241098",
77 | "SAT2_CNDOT_TDOT": "-0.0027234595338783167",
78 | "SAT2_CNDOT_NDOT": "0.14930469999999996"
79 | }
80 | ]
--------------------------------------------------------------------------------
/Orbit_Predictor-BackEnd/api/sample_data/oct5_data/cdm2.json:
--------------------------------------------------------------------------------
1 | [
2 | {
3 | "CCSDS_CDM_VERS": "1.0",
4 | "CREATION_DATE": "2024-10-04T18:19:27.022",
5 | "ORIGINATOR": "CSpoc",
6 | "MESSAGE_ID": "2352_conj2_5250",
7 | "TCA": "2024-10-05T17:10:36.475",
8 | "MISS_DISTANCE": "14475",
9 | "SAT1_OBJECT": "OBJECT1",
10 | "SAT1_OBJECT_DESIGNATOR": "12345",
11 | "SAT1_CATALOG_NAME": "SATCAT",
12 | "SAT1_INTERNATIONAL_DESIGNATOR": "2019-030E",
13 | "SAT1_EPHEMERIS_NAME": "NONE",
14 | "SAT1_COVARIANCE_METHOD": "CALCULATED",
15 | "SAT1_MANEUVERABLE": "NO",
16 | "SAT1_REFERENCE_FRAME": "ITRF",
17 | "SAT1_X": "2569.097065",
18 | "SAT1_Y": "2245.654904",
19 | "SAT1_Z": "6281.497978",
20 | "SAT1_X_DOT": "4.418769571",
21 | "SAT1_Y_DOT": "4.833547743",
22 | "SAT1_Z_DOT": "-3.526774282",
23 | "SAT1_CR_R": "11.30329",
24 | "SAT1_CT_R": "-1736.0459177032214",
25 | "SAT1_CT_T": "1406.0959999999998",
26 | "SAT1_CN_R": "11.450941880675588",
27 | "SAT1_CN_T": "204.01998943137414",
28 | "SAT1_CN_N": "3.319463999999999",
29 | "SAT1_CRDOT_R": "0.0015412149999999955",
30 | "SAT1_CRDOT_T": "-2.1589377491633517",
31 | "SAT1_CRDOT_N": "-0.0001754131058154928",
32 | "SAT1_CRDOT_RDOT": "0.0015412149999999955",
33 | "SAT1_CTDOT_R": "-0.00014407361589835086",
34 | "SAT1_CTDOT_T": "0.0006101929798184829",
35 | "SAT1_CTDOT_N": "-1.3096352890776559e-05",
36 | "SAT1_CTDOT_RDOT": "-5.944793822298559e-10",
37 | "SAT1_CTDOT_TDOT": "1.2788699999999509e-05",
38 | "SAT1_CNDOT_R": "1.2313530455547843e-05",
39 | "SAT1_CNDOT_T": "0.00024296227126581737",
40 | "SAT1_CNDOT_N": "1.6925774251709703e-05",
41 | "SAT1_CNDOT_RDOT": "-3.218738348956514e-10",
42 | "SAT1_CNDOT_TDOT": "-1.4154215132856788e-11",
43 | "SAT1_CNDOT_NDOT": "8.802304999999936e-06",
44 | "SAT2_OBJECT": "OBJECT1",
45 | "SAT2_OBJECT_DESIGNATOR": "67891",
46 | "SAT2_CATALOG_NAME": "SATCAT",
47 | "SAT2_INTERNATIONAL_DESIGNATOR": "2021-030E",
48 | "SAT2_EPHEMERIS_NAME": "NONE",
49 | "SAT2_COVARIANCE_METHOD": "CALCULATED",
50 | "SAT2_MANEUVERABLE": "YES",
51 | "SAT2_REFERENCE_FRAME": "ITRF",
52 | "SAT2_X": "4663.5941299999995",
53 | "SAT2_Y": "-1304.8901919999998",
54 | "SAT2_Z": "25819.095956",
55 | "SAT2_X_DOT": "15.837539142",
56 | "SAT2_Y_DOT": "12646.667095485998",
57 | "SAT2_Z_DOT": "5518.846451436",
58 | "SAT2_CR_R": "68469.38",
59 | "SAT2_CT_R": "-163366625861765.34",
60 | "SAT2_CT_T": "2918959000.0000005",
61 | "SAT2_CN_R": "71915714.42231454",
62 | "SAT2_CN_T": "-23710013898.2132",
63 | "SAT2_CN_N": "2170.8850000000007",
64 | "SAT2_CRDOT_R": "3353.289999999999",
65 | "SAT2_CRDOT_T": "-9788098835180.406",
66 | "SAT2_CRDOT_N": "34493.37382121202",
67 | "SAT2_CRDOT_RDOT": "3353.289999999999",
68 | "SAT2_CTDOT_R": "-1348.4149760644132",
69 | "SAT2_CTDOT_T": "15903045.190279152",
70 | "SAT2_CTDOT_N": "-47.94694154153545",
71 | "SAT2_CTDOT_RDOT": "-18.32563855207348",
72 | "SAT2_CTDOT_TDOT": "0.027011589999999967",
73 | "SAT2_CNDOT_R": "2659.732660844721",
74 | "SAT2_CNDOT_T": "85254697.53509147",
75 | "SAT2_CNDOT_N": "234.04542544087747",
76 | "SAT2_CNDOT_RDOT": "-98.04438819670375",
77 | "SAT2_CNDOT_TDOT": "-0.0027233793844917732",
78 | "SAT2_CNDOT_NDOT": "0.1492901"
79 | }
80 | ]
--------------------------------------------------------------------------------
/Orbit_Predictor-BackEnd/api/sample_data/oct5_data/cdm1.json:
--------------------------------------------------------------------------------
1 | [
2 | {
3 | "CCSDS_CDM_VERS": "1.0",
4 | "CREATION_DATE": "2024-10-04T12:11:49.036",
5 | "ORIGINATOR": "CSpoc",
6 | "MESSAGE_ID": "5372_conj1_5222",
7 | "TCA": "2024-10-05T17:10:36.475",
8 | "MISS_DISTANCE": "14544",
9 | "SAT1_OBJECT": "OBJECT1",
10 | "SAT1_OBJECT_DESIGNATOR": "12345",
11 | "SAT1_CATALOG_NAME": "SATCAT",
12 | "SAT1_INTERNATIONAL_DESIGNATOR": "2019-030E",
13 | "SAT1_EPHEMERIS_NAME": "NONE",
14 | "SAT1_COVARIANCE_METHOD": "CALCULATED",
15 | "SAT1_MANEUVERABLE": "NO",
16 | "SAT1_REFERENCE_FRAME": "ITRF",
17 | "SAT1_X": "2569.097065",
18 | "SAT1_Y": "2245.654904",
19 | "SAT1_Z": "6281.497978",
20 | "SAT1_X_DOT": "4.418769571",
21 | "SAT1_Y_DOT": "4.833547743",
22 | "SAT1_Z_DOT": "-3.526774282",
23 | "SAT1_CR_R": "12.438180000000003",
24 | "SAT1_CT_R": "-2859.9593342082176",
25 | "SAT1_CT_T": "3143.8580000000006",
26 | "SAT1_CN_R": "11.9965192167563",
27 | "SAT1_CN_T": "617.3295806639882",
28 | "SAT1_CN_N": "3.2434519999999987",
29 | "SAT1_CRDOT_R": "0.003560762999999993",
30 | "SAT1_CRDOT_T": "-11.17429351063103",
31 | "SAT1_CRDOT_N": "-0.0006064604505376215",
32 | "SAT1_CRDOT_RDOT": "0.003560762999999993",
33 | "SAT1_CTDOT_R": "-0.00017415206257147947",
34 | "SAT1_CTDOT_T": "-0.0013124245331266138",
35 | "SAT1_CTDOT_N": "-1.3809345175536252e-05",
36 | "SAT1_CTDOT_RDOT": "1.7064174107178696e-09",
37 | "SAT1_CTDOT_TDOT": "1.4082519999999975e-05",
38 | "SAT1_CNDOT_R": "1.9231443904931584e-05",
39 | "SAT1_CNDOT_T": "4.198231150531984e-05",
40 | "SAT1_CNDOT_N": "1.5660229650232342e-05",
41 | "SAT1_CNDOT_RDOT": "-1.582394681859056e-10",
42 | "SAT1_CNDOT_TDOT": "-2.1936230727864984e-11",
43 | "SAT1_CNDOT_NDOT": "8.604338999999885e-06",
44 | "SAT2_OBJECT": "OBJECT1",
45 | "SAT2_OBJECT_DESIGNATOR": "67891",
46 | "SAT2_CATALOG_NAME": "SATCAT",
47 | "SAT2_INTERNATIONAL_DESIGNATOR": "2021-030E",
48 | "SAT2_EPHEMERIS_NAME": "NONE",
49 | "SAT2_COVARIANCE_METHOD": "CALCULATED",
50 | "SAT2_MANEUVERABLE": "YES",
51 | "SAT2_REFERENCE_FRAME": "ITRF",
52 | "SAT2_X": "4663.89413",
53 | "SAT2_Y": "-1329.8901919999998",
54 | "SAT2_Z": "25882.795955999998",
55 | "SAT2_X_DOT": "15.837539142",
56 | "SAT2_Y_DOT": "12646.667095485998",
57 | "SAT2_Z_DOT": "5518.846451436",
58 | "SAT2_CR_R": "68744.55000000002",
59 | "SAT2_CT_R": "-164516051071393.47",
60 | "SAT2_CT_T": "2924869999.999999",
61 | "SAT2_CN_R": "72054166.0937517",
62 | "SAT2_CN_T": "-22723563949.867527",
63 | "SAT2_CN_N": "2171.175",
64 | "SAT2_CRDOT_R": "3360.1049999999996",
65 | "SAT2_CRDOT_T": "-9827853287582.252",
66 | "SAT2_CRDOT_N": "33367.9076961158",
67 | "SAT2_CRDOT_RDOT": "3360.1049999999996",
68 | "SAT2_CTDOT_R": "-1353.7353720994918",
69 | "SAT2_CTDOT_T": "16012071.648502622",
70 | "SAT2_CTDOT_N": "-47.961060064859275",
71 | "SAT2_CTDOT_RDOT": "-18.451090163742993",
72 | "SAT2_CTDOT_TDOT": "0.027021619999999986",
73 | "SAT2_CNDOT_R": "2658.8393484284798",
74 | "SAT2_CNDOT_T": "85496953.755255",
75 | "SAT2_CNDOT_N": "234.09259233918314",
76 | "SAT2_CNDOT_RDOT": "-98.32962134287756",
77 | "SAT2_CNDOT_TDOT": "-0.0027230776774019596",
78 | "SAT2_CNDOT_NDOT": "0.1492"
79 | }
80 | ]
--------------------------------------------------------------------------------
/Orbit_Predictor-BackEnd/api/sample_data/oct5_data/cdm4.json:
--------------------------------------------------------------------------------
1 | [
2 | {
3 | "CCSDS_CDM_VERS": "1.0",
4 | "CREATION_DATE": "2024-10-05T10:59:22.033",
5 | "ORIGINATOR": "CSpoc",
6 | "MESSAGE_ID": "6530_conj4_1024",
7 | "TCA": "2024-10-05T17:10:36.475",
8 | "MISS_DISTANCE": "14510",
9 | "SAT1_OBJECT": "OBJECT1",
10 | "SAT1_OBJECT_DESIGNATOR": "12345",
11 | "SAT1_CATALOG_NAME": "SATCAT",
12 | "SAT1_INTERNATIONAL_DESIGNATOR": "2019-030E",
13 | "SAT1_EPHEMERIS_NAME": "NONE",
14 | "SAT1_COVARIANCE_METHOD": "CALCULATED",
15 | "SAT1_MANEUVERABLE": "NO",
16 | "SAT1_REFERENCE_FRAME": "ITRF",
17 | "SAT1_X": "2569.097065",
18 | "SAT1_Y": "2245.654904",
19 | "SAT1_Z": "6281.497978",
20 | "SAT1_X_DOT": "4.418769571",
21 | "SAT1_Y_DOT": "4.833547743",
22 | "SAT1_Z_DOT": "-3.526774282",
23 | "SAT1_CR_R": "11.506860000000005",
24 | "SAT1_CT_R": "1248.036661134601",
25 | "SAT1_CT_T": "259.5581000000001",
26 | "SAT1_CN_R": "17.916002509223198",
27 | "SAT1_CN_T": "245.7711485383819",
28 | "SAT1_CN_N": "3.8349209999999996",
29 | "SAT1_CRDOT_R": "0.00022726939999999785",
30 | "SAT1_CRDOT_T": "-0.05799527808110657",
31 | "SAT1_CRDOT_N": "-0.00021179419934301688",
32 | "SAT1_CRDOT_RDOT": "0.00022726939999999785",
33 | "SAT1_CTDOT_R": "-0.00015491341398556023",
34 | "SAT1_CTDOT_T": "-0.0015060774553316179",
35 | "SAT1_CTDOT_N": "-2.0952096822822975e-05",
36 | "SAT1_CTDOT_RDOT": "1.4689758301709094e-09",
37 | "SAT1_CTDOT_TDOT": "1.3466079999999983e-05",
38 | "SAT1_CNDOT_R": "2.6661413068295898e-05",
39 | "SAT1_CNDOT_T": "0.00031937485126386103",
40 | "SAT1_CNDOT_N": "2.1609511370945364e-05",
41 | "SAT1_CNDOT_RDOT": "-3.382421549510018e-10",
42 | "SAT1_CNDOT_TDOT": "-3.128348087699916e-11",
43 | "SAT1_CNDOT_NDOT": "1.0368329999999764e-05",
44 | "SAT2_OBJECT": "OBJECT1",
45 | "SAT2_OBJECT_DESIGNATOR": "67891",
46 | "SAT2_CATALOG_NAME": "SATCAT",
47 | "SAT2_INTERNATIONAL_DESIGNATOR": "2021-030E",
48 | "SAT2_EPHEMERIS_NAME": "NONE",
49 | "SAT2_COVARIANCE_METHOD": "CALCULATED",
50 | "SAT2_MANEUVERABLE": "YES",
51 | "SAT2_REFERENCE_FRAME": "ITRF",
52 | "SAT2_X": "4659.49413",
53 | "SAT2_Y": "-1319.7901920000004",
54 | "SAT2_Z": "25850.995956",
55 | "SAT2_X_DOT": "15.837539142",
56 | "SAT2_Y_DOT": "12646.667095485998",
57 | "SAT2_Z_DOT": "5518.846451436",
58 | "SAT2_CR_R": "67978.79000000001",
59 | "SAT2_CT_R": "-162130280857593.2",
60 | "SAT2_CT_T": "2922881000.0",
61 | "SAT2_CN_R": "71696592.84559105",
62 | "SAT2_CN_T": "-25896048808.354958",
63 | "SAT2_CN_N": "2171.0660000000007",
64 | "SAT2_CRDOT_R": "3357.8220000000006",
65 | "SAT2_CRDOT_T": "-9814496912200.408",
66 | "SAT2_CRDOT_N": "37010.43898900366",
67 | "SAT2_CRDOT_RDOT": "3357.8220000000006",
68 | "SAT2_CTDOT_R": "-1332.471339849396",
69 | "SAT2_CTDOT_T": "15492739.465907317",
70 | "SAT2_CTDOT_N": "-47.89920848575927",
71 | "SAT2_CTDOT_RDOT": "-17.854404929205973",
72 | "SAT2_CTDOT_TDOT": "0.026952349999999986",
73 | "SAT2_CNDOT_R": "2657.53461602947",
74 | "SAT2_CNDOT_T": "85412367.9110327",
75 | "SAT2_CNDOT_N": "234.04549101615999",
76 | "SAT2_CNDOT_RDOT": "-98.23003979101016",
77 | "SAT2_CNDOT_TDOT": "-0.0027243267415416354",
78 | "SAT2_CNDOT_NDOT": "0.14929579999999992"
79 | }
80 | ]
--------------------------------------------------------------------------------
/on-orbit-frontend/src/app/login/page.tsx:
--------------------------------------------------------------------------------
1 | "use client"
2 |
3 | import { useState } from 'react';
4 | import { useRouter } from 'next/navigation';
5 | import Footer from '@/components/footer/page';
6 | import Csa from '@/components/csa/csa';
7 |
8 | export default function Login() {
9 | const [email, setEmail] = useState('');
10 | const [password, setPassword] = useState('');
11 | const [error, setError] = useState(null);
12 | const router = useRouter();
13 |
14 | const handleSubmit = async (e: React.FormEvent) => {
15 | e.preventDefault();
16 | setError(null);
17 |
18 | try {
19 | const response = await fetch('http://localhost:8000/api/login/', {
20 | method: 'POST',
21 | headers: {
22 | 'Content-Type': 'application/json',
23 | },
24 | body: JSON.stringify({ email, password }),
25 | });
26 |
27 | if (!response.ok) {
28 | throw new Error(`Login failed. Please check your credentials.`);
29 | }
30 |
31 | const data = await response.json();
32 |
33 | localStorage.setItem('token', data.access_token);
34 |
35 | router.push('/dashboard');
36 | } catch (err) {
37 | if (err instanceof Error) {
38 | setError(err.message);
39 | } else {
40 | setError('An unknown error occurred.');
41 | }
42 | }
43 | };
44 |
45 | return (
46 |
47 |
48 |
49 |
Login
50 |
51 | {error && (
52 |
53 | {error}
54 |
55 | )}
56 |
57 |
93 |
94 | {/*
*/}
95 |
96 |
;
92 | if (error) return Error: {error}
;
93 |
94 | return (
95 |
96 |
97 |
109 |
110 |
Tradespace Trajectory
111 |
112 |
113 |
114 |
115 |
116 |
117 |
96 |
97 |
98 | Dashboards /
99 | Overview
100 |
101 |
102 |
103 | {loading ? (
104 |
105 |
106 |
Loading user info...
107 |
108 |
109 |
110 | ) : error ? (
111 |
112 |
113 |
Error
114 |
{error}
115 |
121 |
122 |
123 | ) : user ? (
124 |
125 |
126 |
127 | User Profile
128 |
129 |
130 |
131 | Email: {user.email}
132 |
133 |
134 | Role: {user.role}
135 |
136 |
137 | Created At:{" "}
138 | {new Date(user.created_at).toLocaleString()}
139 |
140 |
141 |
Enable notifications:
142 |
155 |
156 |
157 |
158 |
159 | ) : (
160 |
161 |
No user data available.
162 |
163 | )}
164 |
165 |
37 | {/* Top navbar for other pages */}
38 | {!isDashboardOrCesium && (
39 |
69 | )}
70 |
71 | {/* Sidebar navbar for dashboard + cesium-view */}
72 | {isDashboardOrCesium && (
73 |
74 |
125 |
126 | )}
127 |
62 |
63 |
Sign Up
64 | {error &&
{error}
}
65 |
145 |
146 |
147 |