Skip to content

Commit

Permalink
fix geodesic again
Browse files Browse the repository at this point in the history
  • Loading branch information
asinghvi17 committed Oct 8, 2024
1 parent a41c257 commit 53936f1
Showing 1 changed file with 2 additions and 2 deletions.
4 changes: 2 additions & 2 deletions ext/GeometryOpsProjExt/segmentize.jl
Original file line number Diff line number Diff line change
@@ -1,6 +1,6 @@
# This holds the `segmentize` geodesic functionality.

import GeometryOps: GeodesicSegments, _fill_linear_kernel!
import GeometryOps: GeodesicSegments, _segmentize, _fill_linear_kernel!
import Proj

function GeometryOps.GeodesicSegments(; max_distance, equatorial_radius::Real=6378137, flattening::Real=1/298.257223563, geodesic::Proj.geod_geodesic = Proj.geod_geodesic(equatorial_radius, flattening))
Expand All @@ -10,7 +10,7 @@ end
# This is the same method as in `transformations/segmentize.jl`,
# but it constructs a Proj geodesic line every time.
# Maybe this should be better...
function _segmentize(method::Geodesic, geom, T::Union{GI.LineStringTrait, GI.LinearRingTrait}; max_distance)
function _segmentize(method::Geodesic, geom, ::Union{GI.LineStringTrait, GI.LinearRingTrait}; max_distance)
proj_geodesic = Proj.geod_geodesic(method.equatorial_radius, method.flattening)
first_coord = GI.getpoint(geom, 1)
x1, y1 = GI.x(first_coord), GI.y(first_coord)
Expand Down

0 comments on commit 53936f1

Please sign in to comment.