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

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

Cope with no dates OR no coords.

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("'null'",None)
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("'null'",None)
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 = 0
74            west = 0
75            north = 0
76            south = 0
77            startdate="January 1, 01 BC"
78            enddate="January 1, 01 BC"
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            try:
91                dates=dgMeta.dgMetadataRecord.dgDataEntity.dgDataSummary.dgDataCoverage.dgTemporalCoverage.DateRange
92                print dates
93            except:
94                print "INFO: XML moles document %s does not contain temporal info." %full_filename
95                no_dates=True
96
97            if no_bbox and no_dates:
98                print "INFO: XML moles document %s does not contain any spatiotemporal info." %full_filename
99                continue
100            if no_dates:
101                pass
102            else:
103                startdate=dates.DateRangeStart
104                enddate= dates.DateRangeEnd
105
106            if no_bbox:
107                pass
108            else:
109                # for bbox in bbox_list:
110                #parse the coordinates somewhat - only use the first bounding box.
111                print bbox_list
112                bbox=bbox_list[0]
113                print bbox
114                #west
115                west = bbox.LimitWest.strip()
116                if west.endswith('E'):
117                    west=bbox.LimitWest.split('E')[0]
118                elif west.endswith('W'):
119                    if west.startswith('-'):
120                        west = bbox.LimitWest.split('W')[0]
121                    else:
122                        west = "-" +bbox.LimitWest.split('W')[0]
123                try:
124                    float(west)
125                except:
126                    print "ERROR:  Will not process File %s. Contains incorrect West bounding box limit." %full_filename
127                    continue
128                #print "West = %s" %west
129                #east
130                east = bbox.LimitEast.strip()
131                if east.endswith('E'):
132                    east=bbox.LimitEast.split('E')[0]
133                elif east.endswith('W'):
134                    if east.startswith('-'):
135                        east = bbox.LimitEast.split('W')[0]
136                    else:
137                        east = "-" +bbox.LimitEast.split('W')[0]
138                try:
139                    float(east)
140                except:
141                    print "ERROR:  Will not process File %s. Contains incorrect East bounding box limit." %full_filename
142                    continue
143                #print "East = %s" %east
144                #north
145                north = bbox.LimitNorth.strip()
146                if north.endswith('N'):
147                    north=bbox.LimitNorth.split('N')[0]
148                elif north.endswith('S'):
149                    if north.startswith('-'):
150                        north = bbox.LimitNorth.split('S')[0]
151                    else:
152                        north = "-" +bbox.LimitNorth.split('S')[0]
153                try:
154                    float(north)
155                except:
156                    print "ERROR: Will not process File %s. Contains incorrect North bounding box limit." %full_filename
157                    continue
158                #print "North = %s" %north
159                #south
160                south = bbox.LimitSouth.strip()
161                if south.endswith('N'):
162                    south=bbox.LimitSouth.split('N')[0]
163                elif south.endswith('S'):
164                    if south.startswith('-'):
165                        south = bbox.LimitSouth.split('S')[0]
166                    else:
167                        south = "-" +bbox.LimitSouth.split('S')[0]
168                try:
169                    float(south)
170                except:
171                    print "ERROR: Will not process File %s. Contains incorrect North bounding box limit." %full_filename
172                    continue
173                #print "North = %s" %south
174
175            if id_exists( Mid ):
176                print "INFO: doc %s exists, updating\n" %Mid
177                do_update( Mid, west, south, east, north, startdate, enddate )
178            else:
179                print "INFO: doc %s does not exist, inserting new record\n" %Mid
180                do_insert( Mid, west, south, east, north, startdate, enddate )
181            numfilesproc += 1
182        else:
183            print "WARNING: File %s appears not to be XML. Will not be processed." %filename
184
185    print 'INFO: SpaceTimeIngestFromMOLES.py ran to end. files processed= %s' %(numfilesproc)
186
187if __name__=='__main__':
188    indir=sys.argv[1]
189    main(indir)
Note: See TracBrowser for help on using the repository browser.