Commit e7d813dc authored by Jeff Piollé's avatar Jeff Piollé

fixed subset extraction for pointcollection

parent cee883f3
......@@ -37,6 +37,7 @@ class PointCollection(AbstractFeature):
latitudes=latitudes,
depths=depths,
times=times)
print "WRITE", self._fields.keys()
def extract_subset(self, boundaries=None, slices=None, indices=None):
"""Extract a subset from the trajectory. The created subset is a new
......@@ -155,6 +156,7 @@ class PointCollection(AbstractFeature):
for geof in self._geolocation_fields:
if self._geolocation_fields[geof]:
output.write_field(self._geolocation_fields[geof])
for dataf in self._fields.keys():
output.write_field(self._fields[dataf])
output.close()
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment