https://github.com/scivision/satcam

Identify what times a satellite should cross a camera pixel and triangulate satellites
https://github.com/scivision/satcam

geoscience matlab

Last synced: 14 days ago
JSON representation

Identify what times a satellite should cross a camera pixel and triangulate satellites

• Host: GitHub
• URL: https://github.com/scivision/satcam
• Owner: scivision
• License: apache-2.0
• Created: 2015-05-05T18:21:29.000Z (over 9 years ago)
• Default Branch: main
• Last Pushed: 2022-04-11T01:44:38.000Z (over 2 years ago)
• Last Synced: 2024-05-02T03:09:32.263Z (7 months ago)
• Topics: geoscience, matlab
• Language: MATLAB
• Homepage:
• Size: 159 KB
• Stars: 6
• Watchers: 3
• Forks: 1
• Open Issues: 0
• Metadata Files:
  • Readme: README.md
  • License: LICENSE

Awesome Lists containing this project

README

        
# Satellite Camera Time/Space synchronization

Using (approximately!) known satellite positions to help verify absolute camera timing when satellite(s) are in field-of-view (FOV).

* RunIridium91Apr11: plot where satellites exist vs. time on x,y pixel map

* RunSatCrossTime: plot when satellite crossed expected pixel (for fine absolute timing verification)

* PlotIridiumTLE: Compare high accuracy Iridium position data with TLE and Optical data

## Examples

I would be happy to discuss your specific needs.

### PlotIridiumTLE.py

You will need to request the high accuracy Iridium data for your data

from the JHU APL AMPERE team. 

You will also need the [TLE data available from AGI](https://www.agi.com/resources/satdb/satdbpc.aspx). 

This file is sort of purpose-specific for now, but can be modified to better fit your needs.

## Algorithms

### RunSatCrossTime

1.  Using TLE and SGP4 propagator, get az/el for satellite from an observer vs. time

2.  Load video frames from the selected times.

3.  Pick a pixel to plot w.r.t. time. Variable "satpix" is an Nx2 matrix of x,y pixel coordinates, and N is the number of time steps.

4.  You will get a 1-D plot of intensity vs. time for that pixel. Did the satellite pixel-crossing correspond to the maximum intensity? If not, there is a timing error that you can quantify by how far in time the maximum intensity was from the TLE/SGP4 prediction.

This code was a quick one-off, if you are actually interested in this let's talk and fix the codebase.