Changeset 312

Timestamp:

02/02/2011 11:07:16 AM (12 years ago)

Author:

marcus

Message:

More station handling (defect).

File:

• SHX/trunk/src/SeismicHandler/modules/stations.py (modified) (4 diffs)

SHX/trunk/src/SeismicHandler/modules/stations.py

@@ -106 +106 @@
     def __getattr__(self, name):
         """
-        Special handling for start and end date.
-        """
-        if name.lower() not in ["start", "end"]:
+        Special handling for start and end date, combination of network and
+        station (netstation) and location, stream and component (lsc).
+        """
+        if name.lower() not in ["start", "end", "netstation", "lsc"]:
             return self.__getattribute__(name)
+
+        if name.lower() == "netstation":
+            return ".".join([self.__dict__["network"], self.__dict__["station"]])
+
+        if name.lower() == "lsc":
+            return ".".join([
+                self.__dict__["location"],
+                self.__dict__["stream"],
+                self.__dict__["component"]
+            ])
 
         # build UTCDateTime from o[n|ff]date[ms]
@@ -206 +217 @@
         Init / refresh station information.
         """
+        # read only information
+        data = []
+        for db in Settings.config.inventory.readonly:
+            data += self.__readDB(db)
+
         # data base access r/w
         self.dbreadwrite = Settings.config.inventory.database[0]
-        data = self.__readDB()
-
-        # read only information
-        data_read = []
-        for db in Settings.config.inventory.readonly:
-            data_read += self.__readDB(db)
-
-        # check for overlapping data from rw to readonly db
-        # first: identical datasets
-        double = set(data) & set(data_read)
-        for i in double:
-            if data[i] != data_read[i]:
-                self.logger.critital("Concurrent data found for channel: %s "
-                                     "WARNING: It's very likely that you use "
-                                     "outdated channel information!",
-                                     i
-                                    )
-            if data_read[i].fir and data_read[i].fir != data[i].fir:
-                self.logger.warning("Concurrent FIR data for channel: %s ",
-                                     i
-                                    )
-
-        print data, data_read
+        data += self.__readDB()
+        
+        for i in data:
+            if not i.netstation in self.stations:
+                self.stations[i.netstation] = {}
+
+            try:
+                self.stations[i.netstation][i.lsc].append(i)
+            except KeyError:
+                self.stations[i.netstation][i.lsc] = [i]
+
+#        print self.stations
 
     def __readDB(self, db=None):
@@ -250 +255 @@
 
         return self.dbsessions[db].query(ChannelMeta).order_by("ondate").all()
+
+    def __getitem__(self, channeldate):
+        """
+        Return channel meta data from channel code and time information.
+        """
+        channel, dt = channeldate
+        print "xxx", channel
+        # check for overlapping data from rw to readonly db
+        # first: identical datasets
+#        double = set(data) & set(data_read)
+#        print double
+#        for i in double:
+#            if data[i] != data_read[i]:
+#                self.logger.critital("Concurrent data found for channel: %s "
+#                                     "WARNING: It's very likely that you use "
+#                                     "outdated channel information!",
+#                                     i
+#                                    )
+#            if data_read[i].fir and data_read[i].fir != data[i].fir:
+#                self.logger.warning("Concurrent FIR data for channel: %s ",
+#                                     i
+#                                    )
+
 
     def add(self, station):
@@ -284 +312 @@
 if __name__ == "__main__":
     stations = Stations()
-
-    GRA1 = ChannelMeta(
-        network="GR",
-        station="GRA1",
-        location="",
-        stream="BH",
-        component="Z",
-        latitude=49.691888,
-        longitude=11.221720,
-        elevation=499.5,
-        depth=0,
-        gain=0.6e9,
-        zeros="(0.0,0.0) (0.0,0.0)",
-        poles="(-0.037004,0.037016) (-0.037004,-0.037016)",
-        start=UTCDateTime("2006-05-10T16:00:00.000100"),
-    )
-
-    GRA12 = ChannelMeta(
-        network="GR",
-        station="GRA1",
-        location="",
-        stream="BH",
-        component="Z",
-        latitude=49.691888,
-        longitude=11.221720,
-        elevation=499.5,
-        depth=0,
-        gain=0.6e9,
-        zeros="(0.0,0.0) (0.0,0.0)",
-        poles="(-0.037004,0.037016) (-0.037004,-0.037016)",
-        start=UTCDateTime("2006-05-10T16:00:00.000"),
-    )
-
-#    print GRA1 == GRA12, GRA1, GRA12
-
-#    stations.add(GRA1)
-
-    print GRA1.ondate, GRA1.start, GRA1.fir
-    
+    print stations[("GR.GRA1..BHZ", UTCDateTime())]