[gdal-dev] dfMaxAngleStepSizeDegrees in ApproximateArcAngles

Joaquim Manuel Freire Luís jluis at ualg.pt
Thu Jun 16 09:27:55 PDT 2022


Even, I wrapped OGR_G_ApproximateArcAngles() in Julia and I think there is something wrong with the  dfMaxAngleStepSizeDegrees argument
(here shortcut to ‘inc’)

This is fine (just repeating the default 0.0 value). Note, I’m dropping the z (=0) argin but that’s only the julia wrapper interface.

D = GMT.Gdal.arcellipse(0, 0, 1., 1., 0., -90., inc=0.)
BoundingBox: [6.123233995736766e-17, 1.0, 0.0, 1.0]
24×2 GMTdataset{Float64, 2}
Row │ x            y
     │ Float64      Float64
─────┼────────────────────────
   1 │ 6.12323e-17  1.0
   2 │ 0.0682424    0.997669
   3 │ 0.136167     0.990686
   4 │ 0.203456     0.979084
   5 │ 0.269797     0.962917
  ⋮  │      ⋮           ⋮
  20 │ 0.962917     0.269797
  21 │ 0.979084     0.203456
  22 │ 0.990686     0.136167
  23 │ 0.997669     0.0682424
  24 │ 1.0          0.0

But if I use for example 10 (or any other number) the dfEndAngle is no longer respected

D = GMT.Gdal.arcellipse(0, 0., 1., 1., 0., -90., inc=10.)
BoundingBox: [6.123233995736766e-17, 0.9993908270190958, -0.17364817766693033, 1.0]
26×2 GMTdataset{Float64, 2}
Row │ x            y
     │ Float64      Float64
─────┼─────────────────────────
   1 │ 6.12323e-17   1.0
   2 │ 0.0697565     0.997564
   3 │ 0.139173      0.990268
   4 │ 0.207912      0.978148
   5 │ 0.275637      0.961262
   6 │ 0.34202       0.939693
   7 │ 0.406737      0.913545
   8 │ 0.469472      0.882948
  ⋮  │      ⋮           ⋮
  19 │ 0.951057      0.309017
  20 │ 0.970296      0.241922
  21 │ 0.984808      0.173648
  22 │ 0.994522      0.104528
  23 │ 0.999391      0.0348995
  24 │ 0.999391     -0.0348995
  25 │ 0.994522     -0.104528
  26 │ 0.984808     -0.173648


Joaquim
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/gdal-dev/attachments/20220616/f43d111a/attachment.htm>


More information about the gdal-dev mailing list