Mercurial > hg > ChinaGisRestApi
view RestDbInterface.py @ 76:e0a79d926902
Path for kml included
author | fknauft |
---|---|
date | Tue, 19 Oct 2010 12:56:34 +0200 |
parents | 2f4c427dec44 |
children | 9c66c0ab395c |
line wrap: on
line source
''' Created on 19.5.2010 @author: casties ''' from OFS.Folder import Folder from Products.PageTemplates.PageTemplateFile import PageTemplateFile from Products.ZSQLExtend import ZSQLExtend import logging import time from zope.interface import implements from zope.publisher.interfaces import IPublishTraverse from ZPublisher.BaseRequest import DefaultPublishTraverse #from zope.publisher.interfaces import NotFound #from zope.app import zapi #from zope.component import queryMultiAdapter import Shared.DC.ZRDB.DA from Products.ZSQLMethods.SQL import SQLConnectionIDs class RestDbInterface(Folder): """Object for RESTful database queries path schema: /db/{schema}/{table}/ omitting table gives a list of schemas omitting table and schema gives a list of schemas """ implements(IPublishTraverse) meta_type="RESTdb" manage_options=Folder.manage_options+( {'label':'Config','action':'manage_editRestDbInterfaceForm'}, ) # management templates manage_editRestDbInterfaceForm=PageTemplateFile('zpt/editRestDbInterface',globals()) # data templates XML_index = PageTemplateFile('zpt/XML_index', globals()) XML_schema = PageTemplateFile('zpt/XML_schema', globals()) XML_schema_table = PageTemplateFile('zpt/XML_schema_table', globals()) HTML_index = PageTemplateFile('zpt/HTML_index', globals()) HTML_schema = PageTemplateFile('zpt/HTML_schema', globals()) HTML_schema_table = PageTemplateFile('zpt/HTML_schema_table', globals()) def __init__(self, id, title, connection_id=None): """init""" self.id = id self.title = title # database connection id self.connection_id = connection_id # create template folder self.manage_addFolder('template') # create data folder #self.manage_addFolder('daten') def getCursor(self): """returns fresh DB cursor""" conn = getattr(self,"_v_database_connection",None) if conn is None: # create a new connection object try: if self.connection_id is None: # try to take the first existing ID connids = SQLConnectionIDs(self) if len(connids) > 0: connection_id = connids[0][0] self.connection_id = connection_id logging.debug("connection_id: %s"%repr(connection_id)) da = getattr(self, self.connection_id) da.connect('') # we copy the DAs database connection conn = da._v_database_connection #conn._register() # register with the Zope transaction system self._v_database_connection = conn except Exception, e: raise IOError("No database connection! (%s)"%str(e)) cursor = conn.getcursor() return cursor def executeSQL(self, query, *args): """execute query with args on database and return all results. result format: {"fields":fields, "rows":data}""" cur = self.getCursor() cur.execute(query, args) # description of returned fields fields = cur.description # get all data in an array data = cur.fetchall() cur.close() return {"fields":fields, "rows":data} def publishTraverse(self,request,name): """change the traversal""" # get stored path path = request.get('restdb_path', []) logging.debug("publishtraverse: name=%s restdb_path=%s"%(name,path)) if name == "index_html": # end of traversal return self.index_html #TODO: should we check more? else: # traverse if len(path) == 0: # first segment if name == 'db': # virtual path -- continue traversing path = [name] request['restdb_path'] = path else: # try real path tr = DefaultPublishTraverse(self, request) ob = tr.publishTraverse(request, name) return ob else: path.append(name) # continue traversing return self def index_html(self,REQUEST,RESPONSE): """index method""" # ReST path was stored in request path = REQUEST.get('restdb_path',[]) # type and format are real parameter format = REQUEST.get('format','HTML').upper() type = REQUEST.get('type',None) # id and doc are used for GoogleMaps content id = REQUEST.get('id',[]) doc = REQUEST.get('doc',None) logging.debug("index_html path=%s format=%s type=%s"%(path,format,type)) RESPONSE.setHeader('Access-Control-Allow-Origin', '*') if type is not None: # non-empty type -- look for template pt = getattr(self.template, "%s_%s"%(format,type), None) if pt is not None: return pt(format=format,type=type,path=path) if len(path) == 1: # list of schemas return self.showListOfSchemas(format=format) elif len(path) == 2: # list of tables return self.showListOfTables(format=format,schema=path[1]) elif len(path) == 3: # GIS if format=="GIS": return self.showGoogleMap(schema=path[1],table=path[2],id=id,doc=doc) if format=="KML_URL": return self.getKmlUrl(schema=path[1],table=path[2],id=id,doc=doc) # table return self.showTable(format=format,schema=path[1],table=path[2]) # don't know what to do return str(REQUEST) def showTable(self,format='XML',schema='public',table=None): """returns PageTemplate with tables""" logging.debug("showtable") pt = getattr(self.template, '%s_schema_table'%format, None) if pt is None: return "ERROR!! template %s_schema_table not found"%format data = self.getTable(schema,table) return pt(data=data,tablename=table) def getTable(self,schema='public',table=None,username='guest'): """return table data""" logging.debug("gettable") data = self.executeSQL("select * from %s.\"%s\" "%(schema,table)) return data def showListOfTables(self,format='XML',schema='public'): """returns PageTemplate with list of tables""" logging.debug("showlistoftables") pt = getattr(self.template, '%s_schema'%format, None) if pt is None: return "ERROR!! template %s_schema not found"%format data = self.getListOfTables(schema) return pt(data=data,schema=schema) def getListOfTables(self,schema='public',username='guest'): """return list of tables""" logging.debug("getlistoftables") # get list of fields and types of db table qstr="""select c.relname FROM pg_catalog.pg_class c LEFT JOIN pg_catalog.pg_namespace n ON n.oid = c.relnamespace WHERE c.relkind IN ('r','') AND n.nspname NOT IN ('pg_catalog', 'pg_toast') AND pg_catalog.pg_table_is_visible(c.oid)""" #qstr="select attname, format_type(pg_attribute.atttypid, pg_attribute.atttypmod) from pg_attribute, pg_class where attrelid = pg_class.oid and pg_attribute.attnum > 0" data=self.executeSQL(qstr) return data def showListOfSchemas(self,format='XML'): """returns PageTemplate with list of schemas""" logging.debug("showlistofschemas") pt = getattr(self.template, '%s_index'%format, None) if pt is None: return "ERROR!! template %s_index not found"%format data = self.getListOfSchemas() return pt(data=data) def getListOfSchemas(self,username='guest'): """return list of schemas""" logging.debug("getlistofschemas") # TODO: really look up schemas data={'fields': (('schemas',),), 'rows': [('public','chgis','mpdl'),]} return data # Methods for GoogleMaps creation def showGoogleMap(self,schema='public',table='mpdl',id=[],doc=None): logging.debug("showGoogleMap") data = self.getDataForGoogleMap(schema,table,id,doc) kmlFileName=self.getKMLname(data=data,table=table) initializeStringForGoogleMaps="""onload=\"initialize(\'http://chinagis.mpiwg-berlin.mpg.de/chinagis/REST/daten/"""+kmlFileName+"""\')\""""#+str(data) initializeStringForGoogleMaps=initializeStringForGoogleMaps.replace("None","0") googleMap_page=str(self.htmlHead())+str(self.getGoogleMapString(initializeStringForGoogleMaps)) return googleMap_page def kml_header(self): return """<?xml version=\"1.0\" encoding=\"UTF-8\"?> <kml xmlns=\"http://www.opengis.net/kml/2.2\" xmlns:gx=\"http://www.google.com/kml/ext/2.2\" xmlns:kml=\"http://www.opengis.net/kml/2.2\" xmlns:atom=\"http://www.w3.org/2005/Atom\">" <Document> <open>1</open> <Style id="sh_dot_red"> <IconStyle> <scale>1.2</scale> <Icon> <href>http://chinagis.mpiwg-berlin.mpg.de/chinagis/images/dot_red.png</href> </Icon> </IconStyle> <ListStyle> </ListStyle> </Style> <Style id="sn_dot_red"> <IconStyle> <Icon> <href>http://chinagis.mpiwg-berlin.mpg.de/chinagis/images/dot_red.png</href> </Icon> </IconStyle> <ListStyle> </ListStyle> </Style> <StyleMap id="msn_dot_red"> <Pair> <key>normal</key> <styleUrl>#sn_dot_red</styleUrl> </Pair> <Pair> <key>highlight</key> <styleUrl>#sh_dot_red</styleUrl> </Pair> </StyleMap> """ def htmlHead(self): return """ <!DOCTYPE html> <html> <head> <meta http-equiv="content-type" content="text/html; charset=utf-8"> <meta http-Equiv="Cache-Control" Content="no-cache"> <meta http-Equiv="Pragma" Content="no-cache"> <meta http-Equiv="Expires" Content="0"> <link rel="stylesheet" type="text/css" href="../../../../mpiwgstyles.css"> <style type="text/css"> html { height: 100% } body { height: 100%; margin: 0px; padding: 0px } #map_canvas { height: 100% } </style> <style type="text/css"> td {height:10; font-size:12;} th {height:10; font-size:12;} b {height:10; font-size:12;} h4 {height:10; font-size:15;} </style> <style type="text/css"> .tooltip { background-color:#ffffff; font-weight:bold; border:2px #006699 solid; } </style> <script type="text/javascript" src="http://maps.google.com/maps/api/js?sensor=false"> </script> <script type="text/javascript"> var map = null; function initialize(data) { kmlURL=initialize.arguments[0]; var myMapOpt = { zoom: 5, center: new google.maps.LatLng(29,116), mapTypeId: google.maps.MapTypeId.TERRAIN }; var map = new google.maps.Map(document.getElementById("map_canvas"), myMapOpt); var myLayerOpt = { map: map, preserveViewport: false }; var newLayer=new google.maps.KmlLayer(kmlURL, myLayerOpt); } </script> </head> """ def getKmlUrl(self,schema='public',table='mpdl',id=[],doc=None): logging.debug("getKmlUrl") data = self.getDataForGoogleMap(schema,table,id,doc) kml=self.getKMLname(data=data,table=table) return """http://chinagis.mpiwg-berlin.mpg.de/chinagis/REST/daten/"""+kml def getDataForGoogleMap(self,schema='public',table='mpdl',id=[],doc=None): logging.debug("getDataForGoogleMap") qstr="SELECT * FROM "+schema+".\""+table+"\"" data={} try: if id!=[]: qstr=qstr+" WHERE " for id_item in id.split(","): if schema=='mpdl' or table.find('mpdl')>-1: qstr=qstr+" mpdl_xmlsource_id = '"+id_item+ "' OR" else: qstr=qstr+" cast(id as text) LIKE '"+id_item+ "' OR" qstr=str(qstr).rsplit(" ",1)[0] #to remove last " OR " data=self.ZSQLSimpleSearch(qstr) return data except: return qstr def getAttributeList(self, table=""): attributeList={} if table != "": sql_for_attributes="SELECT column_name FROM information_schema.columns WHERE table_name = '"+ str(table) + "' AND column_name NOT LIKE 'the_geom'" # print sql_for_attributes results=self.ZSQLSimpleSearch(sql_for_attributes) # print results for result in results: if attributeList=={}: attributeList=[getattr(result,"column_name")] else: attributeList.append(getattr(result,"column_name")) return attributeList def getKMLname(self,data=[],table=""): logging.debug("getKMLname") kml4Marker=self.kml_header() initializeStringForGoogleMaps="" gisID="" xCoord=0 yCoord=0 SQL="" unknownPointCounter=0 kmlFileName="marker"+str(time.time())+".kml" attributeList=self.getAttributeList(table) # Mapping a set of points from table-based SQL-query: if data!=None: try: SQL="SELECT \"attribute with gis_id\" FROM public.metadata WHERE tablename LIKE %r" % (table) gisIDattribute=self.ZSQLSimpleSearch(SQL) except: return "table not registered within metadata" for dataset in data: try: xCoord=getattr(dataset,'longitude') yCoord=getattr(dataset,'latitude') if str(xCoord)=='None': xCoord=0 if str(yCoord)=='None': yCoord=0 except: try: xCoord=getattr(dataset,'x_coord') yCoord=getattr(dataset,'y_coord') if str(xCoord)=='None': xCoord=0 if str(yCoord)=='None': yCoord=0 except: continue2getPoint4GISid="" try: if float(xCoord)==0 and float(yCoord)==0: gisID=getattr(dataset,getattr(gisIDattribute[0],'attribute with gis_id')) coords=self.getPoint4GISid(gisID) if coords!=None: xCoord=coords[0] yCoord=coords[1] except: kml4Marker=kml4Marker + "No Coordinates for:"+gisID +", "+ str(xCoord)+", "+ str(xCoord)+"\n" markerTitel="" kml4Marker+="\n<Placemark>\n" kml4Marker+="<description>\n<![CDATA[<b>" markerTitel="<name>%s</name>\n" % (dataset[2]) i=0 for value in dataset: # if values != (None, None): # if str(values).find('name')>-1: # markerTitel="<name>%s</name>\n" % (values[1]) # continue # elif str(values).find('place')>-1: # markerTitel="<name>%s</name>\n" % (values[1]) # continue # kml4Marker=kml4Marker+": "+str(value) value_string=str(value).replace("'","__Apostroph__") value_string=str(value_string).replace('"','__DoubleApostroph__') value_string=str(value_string).replace(';','__$$__') value_string=str(value_string).replace('&','&') if str(value_string).find('http')>-1: value_string='<A HREF=%r target=_blank>Link to Document</A>' % (value_string) kml4Marker+=attributeList[i]+": "+value_string+"</a><br>\n" i+=1 kml4Marker+="]]>\n</description>\n" kml4Marker+="<styleUrl>#msn_dot_red</styleUrl>\n" kml4Marker+="<Point>" try: if float(xCoord)!=0 and float(yCoord)!=0: kml4Marker+="<coordinates>%s,%s,0</coordinates>" % (xCoord,yCoord) else: unknownPointCounter+=1 kml4Marker+="<coordinates>%s,%s,0</coordinates>" % (125,25+int(unknownPointCounter)*0.4) markerTitel="<name>Un-referenced Location!</name>\n" except: unknownPointCounter+=1 kml4Marker+="<coordinates>%s,%s,0</coordinates>" % (125,25+int(unknownPointCounter)*0.4) markerTitel="<name>Un-referenced Location!</name>\n" kml4Marker+="</Point>\n" kml4Marker+= markerTitel kml4Marker+="</Placemark>\n" kml4Marker=kml4Marker+"</Document>\n</kml>" # kml4Marker=str(kml4Marker).replace('&','$$') # kml4Marker=str(kml4Marker).replace(';','__$$__') # kml4Marker=str(kml4Marker).replace('#','__SHARP__') isLoadReady='false' while isLoadReady=='false': isLoadReady=self.RESTwrite2File(self.daten,kmlFileName,kml4Marker) return kmlFileName def getGoogleMapString(self,kml): logging.debug("getGoogleMapString") printed= '<body %s> '%kml +"""\n <div id="map_canvas" style="width: 98%; height: 95%"> </div> \n </body> \n </html>""" return printed def getPoint4GISid(self,gis_id): j=0 coords=(0,0) if gis_id != None: gis_id_orig=gis_id while (True): j=j+1 if (j>100): # FJK: just to prevent endless loops break if (gis_id.isdigit()): # FJK: regular exit from while-loop break else: gis_id=gis_id.strip('abcdefghijklmnopqrstuvwxyz_') # FJK: to strip all letters gis_id=gis_id.strip() # FJK: to strip all whitespaces resultpoint = [0,0] results = None #try: if int(gis_id)>0: SQL="SELECT x_coord,y_coord FROM chgis.chgis_coords WHERE gis_id LIKE cast("+ str(gis_id) +" as text);" results=self.ZSQLSimpleSearch(SQL) #print results if results != None: for result in results: resultpoint=[getattr(result,str('x_coord')),getattr(result,str('y_coord'))] if resultpoint !=[0,0]: return resultpoint else: coords=self.getCoordsFromREST_gisID(gis_id_orig) if len(coords)>0: try: SQL="INSERT INTO chgis.chgis_coords (gis_id,x_coord,y_coord) VALUES (%s,%s,%s); ANALYZE chgis.chgis_coords;"%(gis_id,coords[1],coords[0]) returnstring=self.ZSQLSimpleSearch(SQL) return coords except: SQL="INSERT INTO chgis.chgis_coords (gis_id,x_coord,y_coord) VALUES (%s,%s,%s); ANALYZE chgis.chgis_coords;"%(gis_id,coords[0][1],coords[0][0]) returnstring=self.ZSQLSimpleSearch(SQL) return coords[0] #except: # return "no coords found" else: return coords[0] def getCoordsFromREST_gisID(self,gis_id): coordlist=[] h=0 while (h<5 and coordlist==[]): h+=1 urlresponse=self.urlFunctions.zUrlopenParseString(self.urlFunctions.zUrlopenRead("http://chgis.hmdc.harvard.edu/xml/id/"+gis_id)) baseDocElement=self.urlFunctions.zUrlopenDocumentElement(urlresponse) childnodes=self.urlFunctions.zUrlopenChildNodes(baseDocElement) itemnodes=self.urlFunctions.zUrlopenGetElementsByTagName(baseDocElement,'item') itemspatialnodes=[] for i in range(0,self.urlFunctions.zUrlopenLength(itemnodes)): itemnode=self.urlFunctions.zUrlopenGetItem(itemnodes,i) itemspatialnodes=self.urlFunctions.zUrlopenGetElementsByTagName(itemnode,'spatial') if itemspatialnodes!=[]: for j in range(0,self.urlFunctions.zUrlopenLength(itemspatialnodes)): coord=[] itemspatialnode= self.urlFunctions.zUrlopenGetItem(itemspatialnodes,j) itemspatiallatnodes=self.urlFunctions.zUrlopenGetElementsByTagName(itemspatialnode,'degrees_latitude') for k in range(0,self.urlFunctions.zUrlopenLength(itemspatiallatnodes)): itemspatiallatnode= self.urlFunctions.zUrlopenGetItem(itemspatiallatnodes,k) coord.append(self.urlFunctions.zUrlopenGetTextData(itemspatiallatnode)) itemspatiallngnodes=self.urlFunctions.zUrlopenGetElementsByTagName(itemspatialnode,'degrees_longitude') for k in range(0,self.urlFunctions.zUrlopenLength(itemspatiallngnodes)): itemspatiallngnode= self.urlFunctions.zUrlopenGetItem(itemspatiallngnodes,k) coord.append(self.urlFunctions.zUrlopenGetTextData(itemspatiallngnode)) coordlist.append(coord) gis_id= "_"+gis_id return coordlist # End for GoogleMaps creation def RESTwrite2File(self,datadir, name,text): # try: fileid=name if fileid in datadir.objectIds(): datadir.manage_delObjects(fileid) newfile=open(name,'w') newfile.write(text) newfile.close() file4Read=open(name,'r') fileInZope=datadir.manage_addFile(id=fileid,file=file4Read) return "Write successful" # except: # return "Could not write" def manage_editRestDbInterface(self, title=None, connection_id=None, REQUEST=None): """Change the object""" if title is not None: self.title = title if connection_id is not None: self.connection_id = connection_id #checkPermission=getSecurityManager().checkPermission REQUEST.RESPONSE.redirect('manage_main') manage_addRestDbInterfaceForm=PageTemplateFile('zpt/addRestDbInterface',globals()) def manage_addRestDbInterface(self, id, title='', label='', description='', createPublic=0, createUserF=0, REQUEST=None): """Add a new object with id *id*.""" ob=RestDbInterface(str(id),title) self._setObject(id, ob) #checkPermission=getSecurityManager().checkPermission REQUEST.RESPONSE.redirect('manage_main')