From 53936f1af7f695aa21e1ff2e05c31a3f7ae56ddf Mon Sep 17 00:00:00 2001 From: Anshul Singhvi Date: Mon, 7 Oct 2024 20:43:39 -0700 Subject: [PATCH] fix geodesic again --- ext/GeometryOpsProjExt/segmentize.jl | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/ext/GeometryOpsProjExt/segmentize.jl b/ext/GeometryOpsProjExt/segmentize.jl index 40935ec44..e933f17a2 100644 --- a/ext/GeometryOpsProjExt/segmentize.jl +++ b/ext/GeometryOpsProjExt/segmentize.jl @@ -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)) @@ -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)