schain-jc Commit - r1199:13d983fa985c

1

import os

2

import sys

3

import glob

4

import fnmatch

5

import datetime

6

import time

7

import re

8

import h5py

9

import numpy

10

11

import pylab as plb

12

from scipy.optimize import curve_fit

13

from scipy import asarray as ar, exp

14

15

SPEED_OF_LIGHT = 299792458

16

SPEED_OF_LIGHT = 3e8

17

18

try:

19

from gevent import sleep

20

except:

21

from time import sleep

22

23

from .utils import folder_in_range

24

25

from schainpy.model.data.jrodata import Spectra

26

from schainpy.model.proc.jroproc_base import ProcessingUnit, Operation, MPDecorator

27

from schainpy.utils import log

28

from schainpy.model.io.jroIO_base import JRODataReader

29

30

def pol2cart(rho, phi):

31

x = rho * numpy.cos(phi)

32

y = rho * numpy.sin(phi)

33

return(x, y)

34

35

FILE_STRUCTURE = numpy.dtype([ # HEADER 48bytes

36

('FileMgcNumber', '<u4'), # 0x23020100

37

('nFDTdataRecors', '<u4'),

38

('OffsetStartHeader', '<u4'),

39

('RadarUnitId', '<u4'),

40

('SiteName', 'S32'), # Null terminated

41

])

42

43

44

class FileHeaderBLTR():

45

46

def __init__(self, fo):

47

48

self.fo = fo

49

self.size = 48

50

self.read()

51

52

def read(self):

53

54

header = numpy.fromfile(self.fo, FILE_STRUCTURE, 1)

55

self.FileMgcNumber = hex(header['FileMgcNumber'][0])

56

self.nFDTdataRecors = int(header['nFDTdataRecors'][0])

57

self.RadarUnitId = int(header['RadarUnitId'][0])

58

self.OffsetStartHeader = int(header['OffsetStartHeader'][0])

59

self.SiteName = header['SiteName'][0]

60

61

def write(self, fp):

62

63

headerTuple = (self.FileMgcNumber,

64

self.nFDTdataRecors,

65

self.RadarUnitId,

66

self.SiteName,

67

self.size)

68

69

header = numpy.array(headerTuple, FILE_STRUCTURE)

70

header.tofile(fp)

71

''' ndarray.tofile(fid, sep, format) Write array to a file as text or binary (default).

72

73

fid : file or str

74

An open file object, or a string containing a filename.

75

76

sep : str

77

Separator between array items for text output. If "" (empty), a binary file is written,

78

equivalent to file.write(a.tobytes()).

79

80

format : str

81

Format string for text file output. Each entry in the array is formatted to text by

82

first converting it to the closest Python type, and then using "format" % item.

83

84

'''

85

86

return 1

87

88

89

RECORD_STRUCTURE = numpy.dtype([ # RECORD HEADER 180+20N bytes

90

('RecMgcNumber', '<u4'), # 0x23030001

91

('RecCounter', '<u4'), # Record counter(0,1, ...)

92

# Offset to start of next record form start of this record

93

('Off2StartNxtRec', '<u4'),

94

# Offset to start of data from start of this record

95

('Off2StartData', '<u4'),

96

# Epoch time stamp of start of acquisition (seconds)

97

('nUtime', '<i4'),

98

# Millisecond component of time stamp (0,...,999)

99

('nMilisec', '<u4'),

100

# Experiment tag name (null terminated)

101

('ExpTagName', 'S32'),

102

# Experiment comment (null terminated)

103

('ExpComment', 'S32'),

104

# Site latitude (from GPS) in degrees (positive implies North)

105

('SiteLatDegrees', '<f4'),

106

# Site longitude (from GPS) in degrees (positive implies East)

107

('SiteLongDegrees', '<f4'),

108

# RTC GPS engine status (0=SEEK, 1=LOCK, 2=NOT FITTED, 3=UNAVAILABLE)

109

('RTCgpsStatus', '<u4'),

110

('TransmitFrec', '<u4'), # Transmit frequency (Hz)

111

('ReceiveFrec', '<u4'), # Receive frequency

112

# First local oscillator frequency (Hz)

113

('FirstOsciFrec', '<u4'),

114

# (0="O", 1="E", 2="linear 1", 3="linear2")

115

('Polarisation', '<u4'),

116

# Receiver filter settings (0,1,2,3)

117

('ReceiverFiltSett', '<u4'),

118

# Number of modes in use (1 or 2)

119

('nModesInUse', '<u4'),

120

# Dual Mode index number for these data (0 or 1)

121

('DualModeIndex', '<u4'),

122

# Dual Mode range correction for these data (m)

123

('DualModeRange', '<u4'),

124

# Number of digital channels acquired (2*N)

125

('nDigChannels', '<u4'),

126

# Sampling resolution (meters)

127

('SampResolution', '<u4'),

128

# Number of range gates sampled

129

('nHeights', '<u4'),

130

# Start range of sampling (meters)

131

('StartRangeSamp', '<u4'),

132

('PRFhz', '<u4'), # PRF (Hz)

133

('nCohInt', '<u4'), # Integrations

134

# Number of data points transformed

135

('nProfiles', '<u4'),

136

# Number of receive beams stored in file (1 or N)

137

('nChannels', '<u4'),

138

('nIncohInt', '<u4'), # Number of spectral averages

139

# FFT windowing index (0 = no window)

140

('FFTwindowingInd', '<u4'),

141

# Beam steer angle (azimuth) in degrees (clockwise from true North)

142

('BeamAngleAzim', '<f4'),

143

# Beam steer angle (zenith) in degrees (0=> vertical)

144

('BeamAngleZen', '<f4'),

145