Rotation Angle Estimation from local system to geocentric system and reference point determination software(V1.0)

%readme1
pillar coordinate data (based on local surveying and GNSS rings) are in the TRANS.m and SolveOAetalWeightSS25m

%how to use
1.run TRANS.m 
If there is a error, you can skip it. This is a inintial step.
    
2. dont execute clear command in the commond window

3. run SolveOAetalWeightSS25.m again and again until the correction values of 
OA phi and lamda become as follows:
>>
SolveOAetalWeightSS25
OA phi and lamda:
0.00000000000000008
0.00000000000000018
0.00000000000000008
128 25 23.5789579030
31 6 51.4484144612
121 13 10.1769533389
OA = 2.2414079618 ;
phi = 0.5430468263;
lamda = 2.1156792809;
sOA= 0.0001503301;
sphi= 0.0004199848;
slamda= 0.0001994914;
>>

4. the rotation martix is in the variable "RT", plese check the result.

%readme2

mean code: DelBatch.m :  Batch to exclude the target point outliers 

input csv file: Target point data in local system

LL0.mat: temp files

compensate or compenstate : used for reweighting method

Readxyzazel.m: read target point data and the AZ EL of the telescope

ipd2023.m: get partials

getFL.m: calc C used in O-C