check lengths match

r-spatial · Dec 5, 2024 · 3a55e12 · 3a55e12
1 parent 3e09b51
commit 3a55e12
Showing 1 changed file with 2 additions and 0 deletions.
diff --git a/src/geodetic.cpp b/src/geodetic.cpp
@@ -50,6 +50,8 @@ Rcpp::NumericVector CPL_geodetic_azimuth(Rcpp::List sfc, double semi_major, doub
     Rcpp::List sfc2(sfc2_);
     if (sfc2.size() < 1)
       stop("bearing needs at least 2 points"); // #nocov
+    if (sfc.size() != sfc2.size())
+      stop("length of x and y must match");
     std::vector<LWGEOM *> lw2 = lwgeom_from_sfc(sfc2);
     for (int i = 0; i < ret.size(); i++) {
       ret[i] = lwgeom_azumith_spheroid((LWPOINT*) lw[i], (LWPOINT*) lw2[i], &s);