Method_of_Gauss/File Reader.py at main · blackbird0831/Method_of_Gauss · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
import numpy as np
import jdcal as jd

# csv reader
def read_csv(name):
    file =  open(name, "r")
    data = file.read()

    data = data.split('\n')
    n = len(data)

    for i in range(len(data)):
        data[i] = data[i].split(' ')
        for j in range(len(data[i])):
            data[i][j] = data[i][j].split(":")
            for k in range(len(data[i][j])):
                data[i][j][k] = float(data[i][j][k])

    file.close()

    # t1 = jd.gcal2jd(data[0][0][0], data[0][1][0], data[0][2][0])
    # t2 = jd.gcal2jd(data[1][0][0], data[1][1][0], data[1][2][0])
    # t3 = jd.gcal2jd(data[2][0][0], data[2][1][0], data[2][2][0])

    # t1 = t1[0] + t1[1] + data[0][3][0]/24 + data[0][3][1]/24/60 + data[0][3][2]/24/3600
    # t2 = t2[0] + t2[1] + data[1][3][0]/24 + data[1][3][1]/24/60 + data[1][3][2]/24/3600
    # t3 = t3[0] + t3[1] + data[2][3][0]/24 + data[2][3][1]/24/60 + data[2][3][2]/24/3600
    t1 = data[0][0][0]
    ra1 = [data[0][1][0], data[0][2][0], data[0][3][0]]
    dec1 = [data[0][4][0], data[0][5][0], data[0][6][0]]
    R1 = [data[0][7][0], data[0][8][0], data[0][9][0]]

    t2 = data[1][0][0]
    ra2 = [data[1][1][0], data[1][2][0], data[1][3][0]]
    dec2 = [data[1][4][0], data[1][5][0], data[1][6][0]]
    R2 = [data[1][7][0], data[1][8][0], data[1][9][0]]

    t3 = data[2][0][0]
    ra3 = [data[2][1][0], data[2][2][0], data[2][3][0]]
    dec3 = [data[2][4][0], data[2][5][0], data[2][6][0]]
    R3 = [data[2][7][0], data[2][8][0], data[2][9][0]]
    return t1, ra1, dec1, R1, t2, ra2, dec2, R2, t3, ra3, dec3, R3

t1, ra1, dec1, R1, t2, ra2, dec2, R2, t3, ra3, dec3, R3 = read_csv("AsteroidData2.txt")
print(t1, ra1, dec1, R1, t2, ra2, dec2, R2, t3, ra3, dec3, R3)