[SOLVED] bad limits in projec.invert() induces error

If you have problems using RGeostats package or if you simply need some help on an unlisted subject.

[SOLVED] bad limits in projec.invert() induces error

Postby bez » Fri Mar 01, 2019 10:52 am

The problem appears when using SI.patches for instance (see below) but comes from projec.invert().

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.
bez
 
Posts: 21
Joined: Tue Oct 23, 2012 10:31 am

Re: bad limits in projec.invert() induces error

Postby Didier Renard » Fri Mar 01, 2019 2:45 pm

Has been corrected.
Will be available with version 11.2.9
Didier Renard
 
Posts: 337
Joined: Thu Sep 20, 2012 4:22 pm


Return to Other Troubleshooting

Who is online

Users browsing this forum: No registered users and 7 guests

cron