A toy example is :
- Code: Select all
toto <- db.create(x1=runif(100,-1,1),x2=runif(100,-1,1),z1=runif(100))
projec.define(db=toto)
SI.patches(toto)
The computation works.
But it crashes when inverting coordinates.
It works without projection.
The problem is that the input data gets longitudes around 0.
But projec.invert only accepts longitude between 1 and 359, which is restrictive.
My suggestion is that projec.invert() should consider longitude from 0 to 360 and latitude from -90 and +90.