Python: Contour features in a polygon shapefiles with a straightforward solution

Question

Dear Members list,

First of all, I apologize to post this question modified and improved from a previous post. Recently i am working on shapefiles polygon in order to compute basic contour features:

1. Area,
2. Perimeter,
3. Area convex hull,
4. Perimeter convex hull,
5. major axis length = gives the length of major axis,
6. minor axis length = gives the length of minor axis,

where major and minor axis length are computed following the Figure:

Using osgeo.gdal, ogr and shapely is possible to load and calculate all indeces but not the major and minor axis length. Reading online solution can be using

1. scikit-image = Measure region properties
2. OpenCV

I am looking a straightforward solution in order to make my code easy and elegant. Some blogs suggest to make an ellipse approximation to my polygon in order to retrieve the major and minor axis length. Is it the best solution?

Any references would be quite helpful. Thanks in advance

---

import osgeo.gdal, ogr
from shapely.geometry import Polygon

shp = osgeo.ogr.Open('../examples/mypoly.shp')
layer = shp.GetLayer()
feature = layer.GetFeature(0)
geometry = feature.GetGeometryRef()
# get area
Area = geometry.GetArea()
pts = geometry.GetGeometryRef(0)
points = []
for p in range(pts.GetPointCount()):
   points.append((pts.GetX(p), pts.GetY(p)))
polygon = Polygon(points)
# get Perimeter
Perimeter = polygon.length
# convex Hull
ConvexHull = polygon.convex_hull
# get Perimeter convex Hull
PerimeterConvexHull = ConvexHull.length
# get Area convex Hull
AreaConvexHull = ConvexHull.area

these are the coordinate vertices of my polygon

polygon = Polygon([(560023.4495758876400000 6362057.3904932579000000),(560023.4495758876400000 6362060.3904932579000000),(560024.4495758876400000 6362063.3904932579000000),(560026.9495758876400000 6362068.3904932579000000),(560028.4495758876400000 6362069.8904932579000000),(560034.9495758876400000 6362071.8904932579000000),(560036.4495758876400000 6362071.8904932579000000),(560037.4495758876400000 6362070.3904932579000000),(560037.4495758876400000 6362064.8904932579000000),(560036.4495758876400000 6362063.3904932579000000),(560034.9495758876400000 6362061.3904932579000000),(560026.9495758876400000 6362057.8904932579000000),(560025.4495758876400000 6362057.3904932579000000),(560023.4495758876400000 6362057.3904932579000000)])

in order to test my code from here:

polygon = Polygon(points)
# get Perimeter
Perimeter = polygon.length
# convex Hull
ConvexHull = polygon.convex_hull
# get Perimeter convex Hull
PerimeterConvexHull = ConvexHull.length
# get Area convex Hull
AreaConvexHull = ConvexHull.area

Answer 1

Hi，i have some work on shapefiles polygon likes your work. And i use the minimum_rotated_rectangle in shapely to define the majoraxis_length and minoraxis_length. and my code is it:

polygon = Polygon(points)  
ConvexHull = polygon.convex_hull  #  convex Hull
#minimum_ rotated_ Rectangle is the minimum boundary rectangle, exterior is the external point coordinates of the minimum boundary matrix, and coords is the coordinates of the points
p1 = Point(ConvexHull.minimum_rotated_rectangle.exterior.coords[0])
p2 = Point(ConvexHull.minimum_rotated_rectangle.exterior.coords[1])
p3 = Point(ConvexHull.minimum_rotated_rectangle.exterior.coords[2])
majoraxis_length  =  p1.distance(p2)
minoraxis_length = p2.distance(p3)