code: def geom2localortho(geom):
      
      cx, cy = geom.Centroid().GetPoint_2D()
      lon, lat, z = cT_helper(cx, cy, 0, geom.GetSpatialReference(), wgs_srs)
      local_srs = localortho(lon,lat)
      local_geom = geom_dup(geom)
      geom_transform(local_geom, local_srs)
      return local_geom
docs: Convert existing geom to local orthographic projection
  
      Useful for local cartesian distance/area calculations