source: TI01-discovery/trunk/ingestAutomation/OAIBatch/SpaceTimeIngestFromMOLES.py @ 2263

Subversion URL: http://proj.badc.rl.ac.uk/svn/ndg/TI01-discovery/trunk/ingestAutomation/OAIBatch/SpaceTimeIngestFromMOLES.py@2263
Revision 2263, 8.4 KB checked in by selatham, 13 years ago (diff)

nocs config etc.

Line 
1#!/usr/bin/env python
2import cElementTree
3import elementtree.ElementTree as etree
4import molesReadWrite as MRW
5import sys
6import db_funcs
7import os
8
9#connect to db (in separate db functions module)
10connection = db_funcs.db_connect()
11
12def listify(item):
13        ''' listify checks if an item is a list, if it isn't it puts it inside a list and returns it. Always returns a list object.'''
14        if type(item) is list:
15            return item
16        else:
17            return [item]
18
19def id_exists(Mid):
20        sql = "select id from spatiotemp where id = '"+Mid+"';"
21        cursor = connection.cursor()
22        try:
23            cursor.execute(sql)
24        except:
25            print "Error: database error %s %s" %(sys.exc_type, sys.exc_value)
26        if len(cursor.fetchall()) <1:
27            return False
28        else:
29            return True
30
31
32def do_insert(Mid,west,south,east,north,startdate,enddate):
33        sql = "INSERT INTO spatiotemp (id, coordinates, startdate, enddate) VALUES ( '"+Mid+ "', sbox'(("+west+"d , "+south+"d), ("+east+"d , "+north+"d))', '"+startdate+"', '"+enddate+"');"
34        sql.replace("nulld","null")
35        print sql
36        cursor = connection.cursor()
37        try:
38            cursor.execute(sql)
39        except:
40            print "Error: database error %s %s" %(sys.exc_type, sys.exc_value)
41        connection.commit()
42
43def do_update(Mid,west,south,east,north,startdate,enddate):
44        sql = "UPDATE spatiotemp SET coordinates = sbox'(("+west+"d , "+south+"d), ("+east+"d , "+north+"d))', startdate='"+startdate+"', enddate= '"+enddate+"', update_time= now() WHERE id='"+Mid+"';"
45        sql.replace("nulld","null")
46        print sql
47        cursor = connection.cursor()
48        try:
49            cursor.execute(sql)
50        except:
51            print "Error: database error %s %s" %(sys.exc_type, sys.exc_value)
52        connection.commit()
53
54def main(indir):
55    if indir == "":
56        sys.exit("USAGE: argument 1 = full path of directory where MOLES records reside")
57    else:
58        print "INFO: moles records are in %s" %indir
59
60#this is a fix to the  ElementTree namespace problem that namespaces are usually represented as ns0, ns1, ns2 etc.
61    etree._namespace_map.update({'http://ndg.nerc.ac.uk/moles': 'moles', 'http://www.w3.org/1999/xlink':'xlink'})
62    numfilesproc = 0
63
64    filenames = os.listdir(indir)
65    for filename in filenames:
66        Mid = filename
67        print Mid
68        if filename.find('.xml') != -1:
69            full_filename = indir + "/" + filename
70            #print full_filename
71            no_bbox = False
72            no_dates = False
73            east = 'null'
74            west = 'null'
75            north = 'null'
76            south = 'null'
77            startdate="null"
78            enddate="null"
79            dgMeta=MRW.dgMetadata()
80            try:
81                dgMeta.fromXML(cElementTree.ElementTree(file=full_filename).getroot())
82            except:
83                print "WARNING: Cannot parse the XML moles document %s. Will not process" %full_filename
84                continue
85            try:
86                bbox_list=listify(dgMeta.dgMetadataRecord.dgDataEntity.dgDataSummary.dgDataCoverage.dgSpatialCoverage.BoundingBox)
87            except:
88                print "INFO: XML moles document %s does not contain a bounding box." %full_filename
89                no_bbox=True
90
91            try:
92                dates=dgMeta.dgMetadataRecord.dgDataEntity.dgDataSummary.dgDataCoverage.dgTemporalCoverage.DateRange
93                print dates
94            except:
95                print "INFO: XML moles document %s does not contain temporal info." %full_filename
96                no_dates=True
97            #if type(dates) =='NoneType':
98            #    no_dates=True
99
100            if no_bbox and no_dates:
101                print "INFO: XML moles document %s does not contain any spatiotemporal info." %full_filename
102                continue
103            if no_dates:
104                pass
105            else:
106                startdate=dates.DateRangeStart
107                enddate= dates.DateRangeEnd
108
109            if no_bbox:
110                pass
111            else:
112                # for bbox in bbox_list:
113                #parse the coordinates somewhat - only use the first bounding box.
114                #print bbox_list
115                bbox=bbox_list[0]
116                #print bbox
117                #west
118                try:
119                    west = bbox.LimitWest.strip()
120                except:
121                    print "ERROR:  Will not process File %s. Contains incorrect West bounding box limit." %full_filename
122                    continue
123                if west.endswith('E'):
124                    west=bbox.LimitWest.split('E')[0]
125                elif west.endswith('W'):
126                    if west.startswith('-'):
127                        west = bbox.LimitWest.split('W')[0]
128                    else:
129                        west = "-" +bbox.LimitWest.split('W')[0]
130                try:
131                    float(west)
132                except:
133                    print "ERROR:  Will not process File %s. Contains incorrect West bounding box limit." %full_filename
134                    continue
135                #print "West = %s" %west
136                #east
137                try:
138                    east = bbox.LimitEast.strip()
139                except:
140                    print "ERROR:  Will not process File %s. Contains incorrect East bounding box limit." %full_filename
141                    continue
142                if east.endswith('E'):
143                    east=bbox.LimitEast.split('E')[0]
144                elif east.endswith('W'):
145                    if east.startswith('-'):
146                        east = bbox.LimitEast.split('W')[0]
147                    else:
148                        east = "-" +bbox.LimitEast.split('W')[0]
149                try:
150                    float(east)
151                except:
152                    print "ERROR:  Will not process File %s. Contains incorrect East bounding box limit." %full_filename
153                    continue
154                #print "East = %s" %east
155                #north
156                try:
157                    north = bbox.LimitNorth.strip()
158                except:
159                    print "ERROR:  Will not process File %s. Contains incorrect North bounding box limit." %full_filename
160                    continue
161                if north.endswith('N'):
162                    north=bbox.LimitNorth.split('N')[0]
163                elif north.endswith('S'):
164                    if north.startswith('-'):
165                        north = bbox.LimitNorth.split('S')[0]
166                    else:
167                        north = "-" +bbox.LimitNorth.split('S')[0]
168                try:
169                    float(north)
170                except:
171                    print "ERROR: Will not process File %s. Contains incorrect North bounding box limit." %full_filename
172                    continue
173                #print "North = %s" %north
174                #south
175                try:
176                    south = bbox.LimitSouth.strip()
177                except:
178                    print "ERROR:  Will not process File %s. Contains incorrect South bounding box limit." %full_filename
179                    continue
180                if south.endswith('N'):
181                    south=bbox.LimitSouth.split('N')[0]
182                elif south.endswith('S'):
183                    if south.startswith('-'):
184                        south = bbox.LimitSouth.split('S')[0]
185                    else:
186                        south = "-" +bbox.LimitSouth.split('S')[0]
187                try:
188                    float(south)
189                except:
190                    print "ERROR: Will not process File %s. Contains incorrect North bounding box limit." %full_filename
191                    continue
192                #print "North = %s" %south
193
194            print "west= %s,south %s, east %s, north %s, startdate %s, enddate %s" %(west,south,east,north,startdate,enddate)
195            if id_exists( Mid ):
196                print "INFO: doc %s exists, updating\n" %Mid
197                do_update( Mid, west, south, east, north, startdate, enddate )
198            else:
199                print "INFO: doc %s does not exist, inserting new record\n" %Mid
200
201                do_insert( Mid, west, south, east, north, startdate, enddate )
202            numfilesproc += 1
203        else:
204            print "WARNING: File %s appears not to be XML. Will not be processed." %filename
205
206    print 'INFO: SpaceTimeIngestFromMOLES.py ran to end. files processed= %s' %(numfilesproc)
207
208if __name__=='__main__':
209    indir=sys.argv[1]
210    main(indir)
Note: See TracBrowser for help on using the repository browser.