Top Banner
Simulating Complex Satellites and a Space-Based Surveillance Sensor Simulation Cody R. Singletary and Francis K. Chun Department of Physics, U.S. Air Force Academy ABSTRACT Maintaining space situational awareness requires the ability to track earth-bound satellites in order to know and predict their position, movement, size, and physical features. However, there are many satellites in orbit that are simply too small or too far away to resolve by conventional optical imaging. We can use photometric techniques to gather information about the body in question, but the problem comes in how we interpret the light curve data. Light curves are created by measuring the intensity of reflected sunlight off of the object as it passes overhead. The intensity is dependent on a variety of factors to include the size, shape, orientation, and material composition of the satellite. When we attempt to solve the inverse problem for light curves, we are attempting to extract information about these different factors. Forward modeling of photometric light curves provides a way to generate a large amount of data under controlled conditions for working the inverse problem and is an effective way to test Non- Resolved Space Object Identification (NRSOI) techniques. Currently, there are few implementations of such modeling programs, one of which only allows simple geometric shapes with the option of antennas. We present our modification to that existing code to create complex models plus our new code to calculate shadowing on the complex object. Then we show the results from the new model and a comparison to the original tool. The next generation of space surveillance sensors will be on satellites. Space based sensors avoid many of the problems of ground based sensors, such as, waiting for lighting conditions to match satellite passes and a night sky. These sensors are restricted only by sun exclusion angle and line of sight around the Earth. This allows for more effective techniques and a much longer time on target. We present an addition to the existing code to consider a sensor in orbit. The code is generalized to provide flexibility in testing different orbital parameters and provides pass prediction with predictive forward modeling for any orbit. 1. INTRODUCTION Space situational awareness is a mission of growing importance for the Air Force and the United States. As more and more countries develop their space capabilities, it will become vital for the U.S. to maintain awareness of all objects orbiting the earth, especially as it pertains to U.S. security and space capability. Being able to track satellites in order to maintain a catalog, will need to evolve into characterizing those satellites. High-resolution imagery of satellites will certainly provide a means of characterization, but as satellites become smaller or as the range to satellites increases, it will become harder and harder to obtain the requisite high-resolution imagery (Fig. 1). Thus other means for characterizing satellites will need to be developed. Research into non-resolvable space object identification techniques is beginning, especially with the adaptation of standard astronomical observational techniques to satellites. Analysis of a satellite’s time-varying photometric signature can provide some information on its form, fit, and function. A satellite’s photometric light curve is the intensity of reflected sunlight off of the object as it passes overhead. The intensity is dependent on a variety of factors to include the size, shape, orientation, and material composition of the satellite. When we attempt to invert a satellite’s light curve, we are attempting to extract information about these different factors. Forward modeling of photometric light curves provides a way to generate a large amount of data under controlled conditions for working the inverse problem and is an effective way to test Non-Resolved Space Object Identification (NRSOI) techniques. Currently there are two ways of generating this data. One is to use ray-tracing on a complex model of the satellite, a task which requires enormous computing power and time. The second way, which we are using here is to create simpler models of satellite bodies and use albedo-area calculations with the Cook-Torrance BRDF model. This is a much faster, albeit less accurate method.
12

Simulating Complex Satellites and a Space-Based ... · PDF fileSimulating Complex Satellites and a Space-Based Surveillance Sensor Simulation Cody R. Singletary and Francis K. Chun

Mar 09, 2018

Download

Documents

trinhkhue
Welcome message from author
This document is posted to help you gain knowledge. Please leave a comment to let me know what you think about it! Share it to your friends and learn new things together.
Transcript
Page 1: Simulating Complex Satellites and a Space-Based ... · PDF fileSimulating Complex Satellites and a Space-Based Surveillance Sensor Simulation Cody R. Singletary and Francis K. Chun

Simulating Complex Satellites and a Space-Based Surveillance Sensor Simulation

Cody R. Singletary and Francis K. Chun Department of Physics, U.S. Air Force Academy

ABSTRACT

Maintaining space situational awareness requires the ability to track earth-bound satellites in order to know and predict their position, movement, size, and physical features. However, there are many satellites in orbit that are simply too small or too far away to resolve by conventional optical imaging. We can use photometric techniques to gather information about the body in question, but the problem comes in how we interpret the light curve data. Light curves are created by measuring the intensity of reflected sunlight off of the object as it passes overhead. The intensity is dependent on a variety of factors to include the size, shape, orientation, and material composition of the satellite. When we attempt to solve the inverse problem for light curves, we are attempting to extract information about these different factors. Forward modeling of photometric light curves provides a way to generate a large amount of data under controlled conditions for working the inverse problem and is an effective way to test Non-Resolved Space Object Identification (NRSOI) techniques. Currently, there are few implementations of such modeling programs, one of which only allows simple geometric shapes with the option of antennas. We present our modification to that existing code to create complex models plus our new code to calculate shadowing on the complex object. Then we show the results from the new model and a comparison to the original tool. The next generation of space surveillance sensors will be on satellites. Space based sensors avoid many of the problems of ground based sensors, such as, waiting for lighting conditions to match satellite passes and a night sky. These sensors are restricted only by sun exclusion angle and line of sight around the Earth. This allows for more effective techniques and a much longer time on target. We present an addition to the existing code to consider a sensor in orbit. The code is generalized to provide flexibility in testing different orbital parameters and provides pass prediction with predictive forward modeling for any orbit.

1. INTRODUCTION

Space situational awareness is a mission of growing importance for the Air Force and the United States. As more and more countries develop their space capabilities, it will become vital for the U.S. to maintain awareness of all objects orbiting the earth, especially as it pertains to U.S. security and space capability. Being able to track satellites in order to maintain a catalog, will need to evolve into characterizing those satellites. High-resolution imagery of satellites will certainly provide a means of characterization, but as satellites become smaller or as the range to satellites increases, it will become harder and harder to obtain the requisite high-resolution imagery (Fig. 1). Thus other means for characterizing satellites will need to be developed. Research into non-resolvable space object identification techniques is beginning, especially with the adaptation of standard astronomical observational techniques to satellites. Analysis of a satellite’s time-varying photometric signature can provide some information on its form, fit, and function.

A satellite’s photometric light curve is the intensity of reflected sunlight off of the object as it passes overhead. The intensity is dependent on a variety of factors to include the size, shape, orientation, and material composition of the satellite. When we attempt to invert a satellite’s light curve, we are attempting to extract information about these different factors.

Forward modeling of photometric light curves provides a way to generate a large amount of data under controlled conditions for working the inverse problem and is an effective way to test Non-Resolved Space Object Identification (NRSOI) techniques. Currently there are two ways of generating this data. One is to use ray-tracing on a complex model of the satellite, a task which requires enormous computing power and time. The second way, which we are using here is to create simpler models of satellite bodies and use albedo-area calculations with the Cook-Torrance BRDF model. This is a much faster, albeit less accurate method.

Page 2: Simulating Complex Satellites and a Space-Based ... · PDF fileSimulating Complex Satellites and a Space-Based Surveillance Sensor Simulation Cody R. Singletary and Francis K. Chun

Report Documentation Page Form ApprovedOMB No. 0704-0188

Public reporting burden for the collection of information is estimated to average 1 hour per response, including the time for reviewing instructions, searching existing data sources, gathering andmaintaining the data needed, and completing and reviewing the collection of information. Send comments regarding this burden estimate or any other aspect of this collection of information,including suggestions for reducing this burden, to Washington Headquarters Services, Directorate for Information Operations and Reports, 1215 Jefferson Davis Highway, Suite 1204, ArlingtonVA 22202-4302. Respondents should be aware that notwithstanding any other provision of law, no person shall be subject to a penalty for failing to comply with a collection of information if itdoes not display a currently valid OMB control number.

1. REPORT DATE SEP 2009 2. REPORT TYPE

3. DATES COVERED 00-00-2009 to 00-00-2009

4. TITLE AND SUBTITLE Simulating Complex Satellites and a Space-Based Surveillance Sensor Simulation

5a. CONTRACT NUMBER

5b. GRANT NUMBER

5c. PROGRAM ELEMENT NUMBER

6. AUTHOR(S) 5d. PROJECT NUMBER

5e. TASK NUMBER

5f. WORK UNIT NUMBER

7. PERFORMING ORGANIZATION NAME(S) AND ADDRESS(ES) U.S. Air Force Academy,Department of Physics,USAF Academy,CO,80840

8. PERFORMING ORGANIZATIONREPORT NUMBER

9. SPONSORING/MONITORING AGENCY NAME(S) AND ADDRESS(ES) 10. SPONSOR/MONITOR’S ACRONYM(S)

11. SPONSOR/MONITOR’S REPORT NUMBER(S)

12. DISTRIBUTION/AVAILABILITY STATEMENT Approved for public release; distribution unlimited

13. SUPPLEMENTARY NOTES 2009 Advanced Maui Optical and Space Surveillance Technologies Conference, 1-4 Sep, Maui, HI.

14. ABSTRACT Maintaining space situational awareness requires the ability to track earth-bound satellites in order toknow and predict their position, movement, size, and physical features. However, there are many satellitesin orbit that are simply too small or too far away to resolve by conventional optical imaging. We can usephotometric techniques to gather information about the body in question, but the problem comes in howwe interpret the light curve data. Light curves are created by measuring the intensity of reflected sunlightoff of the object as it passes overhead. The intensity is dependent on a variety of factors to include the size,shape, orientation, and material composition of the satellite. When we attempt to solve the inverse problemfor light curves, we are attempting to extract information about these different factors. Forward modelingof photometric light curves provides a way to generate a large amount of data under controlled conditionsfor working the inverse problem and is an effective way to test Non- Resolved Space Object Identification(NRSOI) techniques. Currently, there are few implementations of such modeling programs, one of whichonly allows simple geometric shapes with the option of antennas. We present our modification to thatexisting code to create complex models plus our new code to calculate shadowing on the complex object.Then we show the results from the new model and a comparison to the original tool. The next generation ofspace surveillance sensors will be on satellites. Space based sensors avoid many of the problems of groundbased sensors, such as, waiting for lighting conditions to match satellite passes and a night sky. Thesesensors are restricted only by sun exclusion angle and line of sight around the Earth. This allows for moreeffective techniques and a much longer time on target. We present an addition to the existing code toconsider a sensor in orbit. The code is generalized to provide flexibility in testing different orbitalparameters and provides pass prediction with predictive forward modeling for any orbit.

15. SUBJECT TERMS

Page 3: Simulating Complex Satellites and a Space-Based ... · PDF fileSimulating Complex Satellites and a Space-Based Surveillance Sensor Simulation Cody R. Singletary and Francis K. Chun

16. SECURITY CLASSIFICATION OF: 17. LIMITATION OF ABSTRACT Same as

Report (SAR)

18. NUMBEROF PAGES

10

19a. NAME OFRESPONSIBLE PERSON

a. REPORT unclassified

b. ABSTRACT unclassified

c. THIS PAGE unclassified

Standard Form 298 (Rev. 8-98) Prescribed by ANSI Std Z39-18

Page 4: Simulating Complex Satellites and a Space-Based ... · PDF fileSimulating Complex Satellites and a Space-Based Surveillance Sensor Simulation Cody R. Singletary and Francis K. Chun

Fig.

The initial Directoratedetermine observationobject idenviewed frocubesats [2sided cylindifficult to

Our first tarectangularthe code foviable mod

In additiondifferent asspace-basedata for a s

2.1. Compu

To model tindividual distance froorientationshape requthe overall

When we bfor the algofrom the recause shadused in bodcurrently c

The verticeThese vertito arrange

. 1. Degradatio

photometric me of the Air Fosatellite pass gn location. Th

ntification probom closely spac2]. The initial nders, and othe extract inform

ask was to crear solar panels) or shadowing edel, and determ

n to the complespects of obsered sensors planspace based sen

utational Algo

the intensity offacets of the siom the satellite

n of the facet wuires accounting facet intensity

began adaptingorithm. We mueflectance areadowing are lit. dy centered cocalculating the

es of the shadoices define thethe vertices in

on of high-reso

modeling tool wrce Research L

geometry and ahis tool has beeblems such as aced observers [tool included vr regular geom

mation? In answ

ate a technique and make mor

effects. Once tmine whether w

ex body modelirving with a spnned for the futnsor as it obser

2. M

rithm

f a satellite’s phimulated objece to the sensor,

with respect to tg for the shado

y by subtracting

g the code for cust find the fac

a on the shadowThese facets aordinates to printensity.

owing facet are perimeter of thorder that they

olution images

we started withLaboratory (AFallows a user ton used in the p

analyzing when[1]. A version very simple con

metric shapes. Wwering this que

to combine there complex struhis was done, w

we can extract m

ing, we felt it wpace-based sensture. We againrved other sate

MODELING

hotometric lighct for each time, the solar phasthe sun and theowing from oneg out the shado

complex shapecets which cauwed facet. Ourare previously sroject the lit “sh

e now in the plahe shadow, buy appear on the

with increasin

h was developeFRL). The simo define satellitpast by AFRL tn single facets of this tool hanvex shapes suWe ask, can westion, we had

e simple shapeuctures. Once we compared tmore informati

would be necessor, which is esn modified the ellites. This wo

OF COMPLE

ht curve, the cue step in the runse angle, the ree observer. Coe part of the saowing area from

s, we found these shadowing,

r algorithm begspecified by thhadowing face

ane of the facetut are typically e perimeter. Th

ng distance and

d in MATLABmulation uses twte size, shape, mto investigate vbecome visibles also been useuch as cubes, ree learn more wto modify the

es (e.g. rectanguwe had these cthe new resultsion. We discus

ssary to have a specially poignoriginal photomork is discussed

EX SHAPES

urrent tool sumn. The reflecteeflective propermbining simpl

atellite onto anom the facet’s a

ere were a fewcalculate the a

gins by findinghe user for the tets” into the pla

t for which weunordered. Thhe shadowed f

d/or decreasing

B® by the Direwo line elemenmaterial, and avarious non-rese or illuminateed to investigatectangular para

with complex bsoftware.

ular parallelepcomplex bodies to the old andss this work in

tool with whicnant considerinmetric modelind in Section 3.

ms up the refleced light of courrties of the facele shapes into aother part. We

albedo-area pro

w specific, majoarea of the shad whether any ftime being. Thanes of the face

e’re calculatinghe Qhull algorifacet is split int

g satellite size.

ected Energy nts (TLEs) to attitude as well solvable space ed by the sun aste glints from allelepipeds, Nodies, or is it t

iped with s, we had to ad

d to see if this iSection 2.

ch we can test ng the potentialng tool to gene

cted light from rse depends onet, and the a more complee attempt to redoduct.

or portions needow, and subtrfacets which cahe sun unit vecet for which we

g the intensity.ithm is then appto a grid of man

as

s

N-too

dapt is a

l erate

n the

x duce

ded ract it an tor is e are

plied ny

Page 5: Simulating Complex Satellites and a Space-Based ... · PDF fileSimulating Complex Satellites and a Space-Based Surveillance Sensor Simulation Cody R. Singletary and Francis K. Chun

small squaintersection

Out calculadimensionsbody modepanels of 1seen in Fig

2.2. Analy

We presensimulated. occurring aaxial spinnimplementeach time sresults for

In all our pdifferent licorrect behnon-shadowfew possib

The first pa27842 (841is a dawn psatellite waorientationof AEOS arise in periand the expnon-zero re

ares and we findns is summed a

F

ations were coms, a diffuse albel which had th1m×3m×0.05mgs. 3-5 at differ

sis and Results

nt three differen These passes

at orientations ner and a nadir-ted, followed bstep to determithe shadowing

passes aside froght curve from

havior qualitatiwing lightcurv

ble sources. Th

ass (Fig. 3, Eng1 km altitude, 9pass observed fas given a spinn with its spin aand the observeiodic shadowinperimental sateesiduals thus in

d the intersectiand multiplied

Fig. 2. Schema

mpleted with tbedo of 0.3, andhe control mod

m dimensions anrent orientation

s

nt passes of ourwere simulateof 9°, 45°, and-ram pointing s

by a pass with tine magnitude g version comp

om one (not shm that of the noively. These reves. We have nhese will be dis

gagement #1) w90.2º inclinatiofrom AEOS be

nning attitude aaxis toward theer begins to seeng. Interestinglellite are identindicating a diff

ion of the shadby the area of

atic flow for ca

two models, a cd no specular r

del dimensions nd with the samns.

r complex bodyd using parame

d 90° to the termstabilized bodythe shadowing residuals. The

pared to the non

own here), theon-shadowing cesiduals reach not yet determiscussed later.

we present is fon, RAAN of 8etween 1537-15along the axial e observer. Ase the lit and shly however, foical with zero rfference betwee

dow polygon wf each square in

alculating the in

control model teflection. Theas a main body

me albedo. Th

y which are jueters for the AEminator. For thy. We simulatecode impleme

e residuals woun-shadowing v

e shadowing, orcode. While wa maximum ofined where the

for the complex84.54º, and ecc553 UT which direction. As the satellite madowed side or most of the p

residual. Only en the light cur

with the grid of n the grid, givin

ntensity for a g

that was a parae second modely with the addi

he experimenta

st a sample of EOS telescopehe satellite attied a pass witho

ented, then comuld show if ourversion.

r complex-bodwe believe not qf 0.12 stellar msupposed erro

x body in an orcentricity of 0.0

crosses the terseen in Fig. 3,

moves through if the satellite, a

pass, the light cduring the last

rves.

the facet. Finang us the shado

 

given facet.

allelepiped withl was an experiition of two facl, complex bod

the total numbe on Mt. Haleakitude, we ran cout the shadow

mpared the magr code is produ

dy code does shquantitatively c

magnitudes of dor occurs, but c

rbit associated 0010). On Marrminator at a 9the satellite be

its orbit, it passaccounting forcurves of the cot third of the pa

ally, the numbeowed area.

h 1m×1m×3mimental, complce-perpendiculdy model can b

ber of passes wkala with passeases for both a

wing code gnitudes reportucing different

how a slightly correct, shows difference froman attribute it t

with SSN Objrch 28, 2009, t° angle. The

egins with an ses nearly overr the seeminglyontrolled satellass does one se

er of

lex lar be

e es an

ted for

the m the

to a

ect there

rhead y late lite ee

Page 6: Simulating Complex Satellites and a Space-Based ... · PDF fileSimulating Complex Satellites and a Space-Based Surveillance Sensor Simulation Cody R. Singletary and Francis K. Chun

Fig. 3. Engagement #1. The top panel shows the simulated results of the photometric tool. The image of the globe includes the terminator, a sensor site indicated by the red plus sign, and the satellite pass in yellow moving from north to south. The complex body image is the orientation of the satellite as seen from the sensor site. The plot of intensity (I) versus time shows the overall intensity in white and individual facet intensities in blue. Finally, the bottom panel is a plot of residuals or difference between the complex body with shadowing and the complex body without shadowing.

The second pass (Fig. 4, Engagement #2) we present is for the complex body in an orbit associated with SSN Object 28773 (560 km altitude, 31.41º inclination, RAAN of 335.71º, and eccentricity of 0.0007). On February 5, 2007, there is a dusk pass observed from AEOS between 0430-0500 UT which crosses the terminator at a 90° angle. This satellite was also given an axial spinner attitude. Initially, the satellite begins the pass oriented with its spin axis toward the observer. The residuals in the first third of the light curve have a higher frequency than the actual rotation rate of the satellite, however the last two-thirds of the residual plot shows shadowing occurring at the rotation frequency of the satellite. Further analysis is required to understand the difference in the shadowing frequency.

Page 7: Simulating Complex Satellites and a Space-Based ... · PDF fileSimulating Complex Satellites and a Space-Based Surveillance Sensor Simulation Cody R. Singletary and Francis K. Chun

 

Fig. 4. Engagement #2. Same format as Fig. 3.

The final pass (Fig. 5, Engagement #3) we present is for the complex body in an orbit associated with SSN Object 27844 (1171 km altitude, 42.2º inclination, RAAN of 83.55º, and eccentricity of 0.0011). On March 13, 2009, there was a dawn pass observed at AEOS between 1545-1558 UT which crosses the terminator at a 45° angle. Contrary to the previous cases, this satellite was in a nadir-ram stabilized attitude with its body x-axis pointing to the ram and its body z-axis pointing to nadir. For this engagement, the satellite begins the pass oriented with its ram axis nearly toward the observer. As the satellite moves through its orbit, the observer begins to see an increasing and systematic shadowing effect as indicated in the residual plot. One does not see any periodic shadowing due to satellite spin since the satellite is not spinning. Additionally, as the satellite passes the observer and across the terminator, the “solar panel” appendages appear to reduce the overall intensity from the side facets of the main body as one would expect. 

Page 8: Simulating Complex Satellites and a Space-Based ... · PDF fileSimulating Complex Satellites and a Space-Based Surveillance Sensor Simulation Cody R. Singletary and Francis K. Chun

 

 

Fig. 5. Engagement #3. Same format as Fig. 3.

We may draw three basic conclusions from this data shown above. First, the shadowing code appears to account for one facet shadowing another facet. Second, the shadowing effect appears consistent with the satellite attitude as seen by the marked differences in the residuals between the axial-spinner and nadir-ram pointing attitudes. Finally, we were able to “measure” the shadowing effect by comparing a complex body with shadowing to one without shadowing. However, when observing an actual satellite in orbit and measuring its time-varying photometric light curve, one must devise a way to determine whether the object is indeed a complex shape exhibiting shadowing, or a simple object with no shadowing.

 

3. SPACE-BASED SENSOR MODELING

3.1 Pass Prediction and Visualization

The photometric tool created by AFRL and modified as described above to accommodate complex shapes has one limitation. Currently, it can only model satellites as observed from fixed ground sites. Theoretically, this tool can model a sensor site at a fixed longitude and latitude, but at an altitude representative of a satellite orbit. Of course, in doing so, one does not even come close to modeling the effects of the sensor movement on the target’s photometric light curve. With future space-based sensors planned, we decided to modify the AFRL tool to model a satellite’s photometric light curve as measured by another orbiting satellite. Most of the tool’s functions and routines were still applicable, except that now we had to modify the code to accept orbit parameters for the sensor satellite and predict when the sensor satellite will be able to observe the target satellite under various constraints such as a sun-exclusion angle and no measurements against the hard earth background. The first step in the modification was to rewrite the pass prediction code. The current pass prediction algorithm takes two-line elements (TLEs) of the observing satellite and propagates them using SGP4. It then finds the TLEs of the target satellite and propagates it using the observing satellite’s coordinates as the “ground observer” coordinates. Once both orbits were propagated, we had to calculate all ephemeris points where the Earth and the sun

Page 9: Simulating Complex Satellites and a Space-Based ... · PDF fileSimulating Complex Satellites and a Space-Based Surveillance Sensor Simulation Cody R. Singletary and Francis K. Chun

excluded thoriginal co Once the avisualizatioin Fig. 6. IEarth’s surwhile the tthan the tarwhile the psatellite orbreference. It isn’t unti

Fig. 6at a lothe satvector

The secondprojection the same oblue when when viewThe sensorbeginning the hard ea

he view of the ode, so no modi

accepted passeson of the geomIt is interactiverface for sky loarget satellite orget satellite. Tportion that is vbit. The globOne can see th

il the end of th

. Globe visualower altitude tatellite target is r is green and th

d visualization along with the

orbital data as inviewing is not

wing is not possr-to-target vectof the orbits, th

arth backgroun

target. The epification of thi

s were computemetry. The firste with the user ocation referencorbit is the red-That portion ofviewable is in ybe visualizationhat for most of

he orbit does th

lization of the arget satellite. yellow when

he Sun unit ve

method plots te vector pointinn Fig. 6 with tht possible and dsible and red wtor is green andhe satellite sen

nd.

phemeris wheres aspect was ne

ed, we needed t is a three-dimable to rotate tce. As seen in-yellow trace. f the target sateyellow. The sun also shows a f the target satehe hard-earth ba

viewing possib The orbital paviewing is posctors are orang

the orbital pathng from the senhe Earth locatedark blue when

when possible. d the Sun’s locnsor cannot view

e the target waecessary.

a way to visuamensional viewthe view to any

n Fig. 6, the satFor this case,

ellite orbit thatun unit vector igreen line from

ellite orbit, the ackground exc

bility (in RA-Dath of the satelssible and red wge arrows.

hs on a right asnsor to target aned in the centern possible, whiThe beginningation and excluw the target sa

as actually lit w

alize them. Thw from far awayy angle and theellite sensor’s the satellite se

t cannot be viewis indicated bym the satellite ssatellite sensorlusion restrict

Dec projectionllite sensor is dwhen it is not p

scension, declinnd the locationr. The orbital pile the orbital pg the orbits are usion area is ye

atellite due to th

was previously

here are two basy, the “globe” ve celestial spheorbit is indicatnsor was in a hwed by the sate

y the orange arrsensor to the tar can observe tthe viewing.

 

n) for a satellitedark blue. Thepossible. The

nation (RA-Den of sun excluspath of the satepath of the targ

indicated by thellow. One canhe restriction o

calculated by t

sic methods fovisualization shere projected toted by the bluehigher altitude ellite sensor is rows from the target satellite fothe target satell

e sensor lookine orbital path osensor-to-targ

ec) Mollweide ion. Fig. 7 sho

ellite sensor is get satellite is bhe green squarn see that at th

of viewing agai

the

or hown

o the e trace

orbit in red target

for lite.

ng of

get

ows light

black res. he inst

Page 10: Simulating Complex Satellites and a Space-Based ... · PDF fileSimulating Complex Satellites and a Space-Based Surveillance Sensor Simulation Cody R. Singletary and Francis K. Chun

Fig. 7centerwhen when and th

We found percent of approximawanted to s

Once the bmodeling. bidirectionminimal chorbit, follocoordinate

3.2 Photom

Once the mproducing have a chacurves thatshow a disobserver bto the anglasymptotespasses for object with0.15, for a

All runs ofapproximaplots, one oshowed theangles. Thalso showepoints belo

7. Visualizatior. The orbital possible, whilpossible. The

he Sun’s locatio

the pass predicruntime depen

ately six hours. simulate differ

basic pass predi All that was r

nal reflectance hange to the BRowed by the tars.

metric Data fro

modifications tocorrect resultsracteristic diffet depend only otinct diffusive eing 90° to thee between a fas of 90° and 18each of three s

h albedo of 0.3total of 6 case

f the six cases aately 300 differof a cube with e defined lowehis feature is med intensities bow the lower bo

on of the viewpath of the sae the orbital pgreen squares

on and exclusio

ction to be the nding on other s

We added a frent satellite sh

iction functionequired was todistribution funRDF calculatiorget satellite’s o

m a Space-Bas

o the original s. Earlier resea

ferences [1]. Spon the shape ofcutoff angle at

e sun. For a tetace and an edge80° correspondshapes: cube, te, and separatels.

appeared to prorent passes for specular compr bound curves

more evident inrighter than thound curve wh

wing possibilityatellite sensor path of the targ

indicate the beon area is yello

most computatselected optionfeature to save hapes for the sa

nality was compo calculate the vnction (BRDF)on code. Howeorbit, using the

sed Sensor

software was carch has shownpecifically, it af the object andt 90° solar phatrahedron, the ce [3]. Cylinderding to the end etrahedron, andly with an analy

oduce phase aneach case, whi

ponent and the s which appearn the cube than e lower bound

hich we are cur

y in RA-Dec Mis light blue w

get satellite is eginning of theow.

tionally intensins. For examplthe pass inform

ame passes, thu

pleted, we movvector from the) calculation asever to do so, we ECI coordina

omplete, we ran that the phaseappears that evd not on its sizese angle correscutoff solar phrs have two lowcaps and cylin

d cylinder. Eacytical BRDF w

ngle-magnitudeich gave a veryother a purely

r to asymptotethe cylinder dcurve as expec

rrently investig

Mollweide prowhen viewing black when vie orbits. The s

ive part of the le, a pass predimation to a fileus saving time

ved on to modie sensor to the s the observatiowe had to first ates of the satel

an test cases to e angle-magnituery object has e or material psponding to the

hase angle is 54wer bound or mndrical sides, rech of these wer

with diffuse alb

e distributions y wide distribut

diffuse cylindto dimmer mag

data. The cube cted. There ar

gating.

ojection with tis not possible

iewing is not psensor-to-targe

code, typicallyiction for six m

e so that one caand obtaining

ification of thetarget and pas

on vector. Doipropagate the llite sensor as t

verify that theude distributiona distinct loweroperties. For e face of the cu4.7° (or minimum intenespectively. Wre also run as a

bedo of 0.3 and

as expected. Wtion of phase a

der are shown ingnitude at the cwhich had a sp

re however som

 

the Earth in the and dark blupossible and ret vector is gree

y taking 90-95 months took an then load it iresults faster.

e photometric s that into the ing so would rsatellite sensorthe observer

e modificationsn for various ser bound brightinstance, cube

ube visible to a) correspond

nsity curves witWe ran six mont

a purely diffused specular albed

We plotted angles. Two ofn Fig. 8. All correct phase pecular compo

me unexpected

he ue ed en

if one

equire r’s

s were hapes tness es an ding th ths of e do of

f these

onent data

Page 11: Simulating Complex Satellites and a Space-Based ... · PDF fileSimulating Complex Satellites and a Space-Based Surveillance Sensor Simulation Cody R. Singletary and Francis K. Chun

Fig. 8The pcompo

In our testiobservationbased sensan orbitingin a stabilizdetermine higher conand attitudsatellites th

Modificatito be worklight curvebodies, finFor the spaon-target fometrics to

The authorand their oPerformanDetachmenResearch PMr. Joe Ko

1. Hall D.

. Phase angle-hase angle is onent on the cu

ing, we noticedn pass. This isor would. Add

g sensor sees mzed nadir-ram satellite shapes

nfidence one cades on each pashan a ground b

ons to an existking properly. es. For the shadding a method ace-based sensofor different cladefine these as

rs would like tooverall commence Computing nt 15, Air ForcProgram under oesters of the S

et. al., “Separa

-magnitude disthe x axis wh

ube compared t

d the orbital sens due to the satditionally, beca

many more targpointing attituds, the more pha

an have in the iss, an orbiting sbased sensor.

ting photometriThere is still hdowing modelto generally ca

or model, futurasses of orbits, spects.

o thank Dr. Donts, suggestionModernization

ce Research Labwhich much o

Sensors Directo

ating Attitude

stributions for ahile the magnitto the purely di

nsor has the inellite sensor mause of the widet satellite attitde. Since it apase angle coverdentification. sensor could po

4. CO

ic tool to allowhowever some a, future work wall the shadowre work includclassifying av

5. ACKNO

oyle Hall and Mns, and recommn Program Offiboratory for th

of this researchorate, Air Forc

6. R

and Shape Effe

a cube (left plotude is the y aiffuse cylinder

nherent propertymoving much fade geometry chtudes than a grppears that pharage one can oAdditionally, botentially allow

ONCLUSION

w for complex sanalysis requirwould include c

wing function, ades combining iailable coverag

OWLEDGEM

Mr. Paul Kervinmendations. Wfice, the Maui Hheir support to th was conductee Research Lab

REFERENCES

fects for Non-re

ot) and cylindeaxis. Notice thr.

y of a wider coaster and farthehanges inherentround-based sese angle-magn

obtain to discerbecause of the

w for quicker id

NS

satellite shapesred to fully undcreating more

and creating a git with the shadge for different

MENTS

n for the use ofe also thank th

High Performanthe U.S. Air Fod. Finally, thisboratory.

S

esolved Object

er (right plot) frhe structure fro

overage of phaser during a passt in the orbits pnsor, especiall

nitude distributirn the lower bo

wider coveragdentification o

s and a space-bderstand the simand different tygeneralized comdowing versiont classes of orb

f their photomehe Department nce Computingorce Academy s research was

ts,” The 2007 A

from 300 passeom the specula

se angles for evs than a groundpassing each oty when the tarions can be use

ound curve, thege of phase angf shape for unk

based sensor apmulated photomypes of compomposite body Gn, researching tbits, and develo

etric modeling of Defense Higg Center and Cadet Summefunded in part

AMOS Technic

es. ar

very d-ther, get is ed to

e gles known

ppear metric osite GUI. time-oping

tool gh

er t by

cal

Page 12: Simulating Complex Satellites and a Space-Based ... · PDF fileSimulating Complex Satellites and a Space-Based Surveillance Sensor Simulation Cody R. Singletary and Francis K. Chun

Conference Proceedings, Kihei, HI, 2007.

2. Hall D., “Optical CubeSat Discrimination,” The 2008 AMOS Technical Conference Proceedings, Kihei, HI, 2008. 3. "Tetrahedron." Wikipedia, The Free Encyclopedia. 26 Aug 2009, 21:44 UTC. 26 Aug 2009 <http://en.wikipedia.org/w/index.php?title=Tetrahedron&oldid=310245038>.