[SOLVED] Neighbourhood and Anisotropy Rotation

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

[SOLVED] Neighbourhood and Anisotropy Rotation

Postby mharley » Wed Jan 21, 2015 4:12 pm

Hi Didier
I did some 'experimentation' with the GUI.neigh() tool. The Anisotropy rotation using Angles appears to not be operative yet because after saving the neighbourhood object I see that it defaults to the unit rotation matrix when I display the neighbourhood as below.

Neighborhood characteristics
============================
Moving neighborhood option
Space dimension = 3
Minimum number of samples = 5
Maximum number of samples = 64
Number of angular sectors = 8
Maximum number of points per sector = 8
Maximum horizontal distance = 100.000000
Anisotropy : 1.500 1.000 0.800
Rotation :
[,1] [,2] [,3]
[1,] 1.000 0.000 0.000
[2,] 0.000 1.000 0.000
[3,] 0.000 0.000 1.000
>


In the tool the anisotropy rotation is around X, then around Y and lastly around Z, are these rotation angles ANTI-clockwise around the axes when looking along the axis towards the origin (i.e. do they correspond to the Mathematician convention within Isatis)?

I ran util.ang2mat process and recovered this matrix:
> rotmat
[,1] [,2] [,3]
[1,] 0.7500000 -0.2165064 0.6250000
[2,] 0.4330127 0.8750000 -0.2165064
[3,] -0.5000000 0.4330127 0.7500000

Which I manually entered into the GUI.neigh tool. But when I hit Save, it reverts to the unit rotation matrix.

Regards
MH
mharley
 
Posts: 3
Joined: Thu Nov 27, 2014 1:41 pm

Re: Neighbourhood and Anisotropy Rotation

Postby Didier Renard » Sat Feb 28, 2015 6:13 pm

This is now fixed. It will be published in Version 10.0.9
Didier Renard
 
Posts: 242
Joined: Thu Sep 20, 2012 4:22 pm


Return to Other Troubleshooting

Who is online

Users browsing this forum: No registered users and 1 guest

cron