Actually it is pretty easy:
Here is the log of pdbset
[ccp4@roo job_55]$ pdbset xyzin part.pdb
...
Logical name: XYZIN File name: part.pdb
PDB file is being opened on unit 1 for INPUT.
MATRICES DERIVED FROM CRYST1 CARD IN COORDINATE FILE
RF RO
0.016 0.009 -0.000 -0.000 61.922 -30.961 0.000 -0.000
-0.000 0.019 -0.000 0.000 0.000 53.626 0.000 0.000
0.000 -0.000 0.004 0.000 0.000 0.000 248.752 -0.000
-0.000 0.000 -0.000 1.000 -0.000 0.000 -0.000 1.000
rota euler 10 20 30
Data line--- rota euler 10 20 30
end
Data line--- end
Logical name: XYZOUT File name: XYZOUT
PDB file is being opened on unit 2 for OUTPUT.
Coordinates will be transformed as follows:
( 0.714610 -0.613092 0.336824 ) ( x ) ( 0.000 )
( 0.633718 0.771281 0.059391 ) ( y ) + ( 0.000 )
( -0.296198 0.171010 0.939693 ) ( z ) ( 0.000 )
....
Of course you still have to worry about the orthogonalisation code used for the SELFROT search.