I have sets of lat/lon points for vortices in the ocean, and what I want to do is draw curves around them (they are roughly elliptical), and determine their eccentricity and orientation as they move through the ocean in time. I'm using python, and I'm on a mac.

I am having a really dreadful time navigating the sea (pun intended) of special libraries that are intended for this purpose.

Can anyone provide me with a simple list of instructions to do this?

I found some tools that point to Fiona and Shapely as easy-to-use tools, but I cannot for the life of me figure out how to get this array of lat-lon points into a .shp file which Fiona can open.

1) How do I create a .shp file out of these data points? If there is another filetype I should use, please let me know how. This is as far as I've gotten:

import Shapefile
w = shapefile.Writer(shapeType = 8) #8 is multipoint

From here I'm totally stuck, and can't find documentation that helps me past this. I don't need any additional features, literally just the points.

I have tried the following as well: w = shapefile.Writer()

for i in range(vortex.shape[0]):
    w.point(vortex[i, 0], vortex[i, 1])

for i in range(0, 327):
    print w.shapes()[i].points


Vortex is a numpy array with dimensions (327, 2), with each row containing a lat and a lon value (duh).

This prints out a list of the correct points, so I'm on the right track. I found this snippet of code, which I think should plot my points:

import fiona
import shapely.geometry as geometry
input_shapefile = 'sample_vortex.shp'
shapefile = fiona.open(input_shapefile)
points = [geometry.shape(point['geometry'])
      for point in shapefile]

import pylab as pl
x = [p.coords.xy[0] for p in points]
y = [p.coords.xy[1] for p in points]
pl.plot(x,y,'o', color='#f16824')

But it actually does not, and yields a blank plot, which tell me, I think, that my .shp file isn't correct.

Any ideas?

2) Do you know how to draw these ellipses and get the eccentricity (and how much the major axis deviates from horizontal, while we're at it)?

closed as too broad by PolyGeo Dec 8 '17 at 20:35

Please edit the question to limit it to a specific problem with enough detail to identify an adequate answer. Avoid asking multiple distinct questions at once. See the How to Ask page for help clarifying this question. If this question can be reworded to fit the rules in the help center, please edit the question.

  • Where is your lng/lat list of coordinates coming from CSV file or hard coded in python? – artwork21 Dec 8 '17 at 17:42
  • It's a numpy array. I added details about my efforts above. – BenL Dec 8 '17 at 17:47
  • Check out this q/a, gis.stackexchange.com/questions/52705/… – artwork21 Dec 8 '17 at 17:54
  • This isn't helpful actually; I have a bunch of numbers; I do not have data in any geospatial format. All of these answers deal directly with data already in that format, or require the use of osgeo, which for whatever reason doesn't work on my machine (despite following multiple sets of instructions to install it, I get the following error every time I try to import it dlopen(/Users/ben/anaconda/lib/python2.7/site-packages/osgeo/_gdal.so, 2): Library not loaded: @rpath/libjpeg.8.dylib Referenced from: /Users/ben/anaconda/lib/libgdal.20.dylib Reason: image not found – BenL Dec 8 '17 at 17:58
  • As per the Tour there should be only one question asked per question. – PolyGeo Dec 8 '17 at 20:35

Fiona and Shapely may be red herrings in your case. You can plot the points in your vortex array without using either of them. You could use SciPy's ConvexHull to get the envelope of the points. Shapely can also compute a convex hull but requires more data gymnastics: you'd need to make a MultiPoint object of your points and call its convex_hull method.

I'm at a loss about fitting an ellipse to the convex hull, but this algorithm (in Python) looks like it could be a start: http://shortrecipes.blogspot.com/2008/12/python-ellipse-fitting.html.

Not the answer you're looking for? Browse other questions tagged or ask your own question.