How to make a virtual layer from ogr geometry?

Question

I know how to do this:

nodes = driver.Open(node_nav,0)
node_lyr = nodes.GetLayer()  
node_lyr.SetSpatialFilter(None)

But how do can use the layer methods on a geometry I created from scratch?

multip=ogr.Geometry(ogr.wkbMultiPoint)
for n in nodes:
    p=ogr.Geometry(ogr.wkbPoint)
    x,y = n.node
    x = float(x)/dim
    y = float(y)/dim

    p.AddPoint(x,y)
    multip.AddGeometry(p)

This is the traceback for this attempt:

File "/usr/lib/python2.7/dist-packages/osgeo/ogr.py", line 2927, in <lambda>
 __getattr__ = lambda self, name: _swig_getattr(self, Geometry, name)
File "/usr/lib/python2.7/dist-packages/osgeo/ogr.py", line 55, in _swig_getattr
raise AttributeError(name)
AttributeError: SetSpatialFilter

I don't want to save the geometry in a shape file (like it is described here for example).

Answer 1

You could write the geometry to an in-memory dataset using the '/vsimem/path_file.shp' virtual file system.

e.g.

drv = ogr.GetDriverByName( 'ESRI Shapefile' )
ds = drv.CreateDataSource( "/vsimem/blahblah.shp" )