Question

Hi,Dago,

I represent my question here.

I have a new idea about PPP, and I need to do the following transformation: Assume there is a station with position [ecefStaX, ecefStaY, ecefStaZ] in ECEF, and a satellite with position [ecefSatX, ecefSatY, ecefSatZ] in ECEF. The station's topocentric system is UEN, I want to rotate the satellite position vector around the U axes. I get a hint from another label "DgpsQuestion" and here is my codes:

// station position
Position staPos(ecefStaX, ecefStaY, ecefStaZ);   
// satellite position
Position satPos(ecefSatX, ecefSatY, ecefSatZ);  

// the vector between station position and satellite position
Triple ecefDXYZ(satPos[0]-staPos[0],satPos[1]-staPos[1],ecefDetZ=satPos[2]-staPos[2]);

// transform ecefDXYZ from ECEF to UEN
Triple uen = ecefDXYZ.R3(+staPos.longitude()).R2(-staPos.geodeticLatitude());

// rotate angle around U axes
uen.R3(angle);

// transform the vector from UEN to ECEF
Triple xyz = enu.R2(+staPos.geodeticLatitude()).R3(-staPos.longitude());

// get the desire satellite position
Position satPos2(staPos[0]+xyz[0],staPos[1]+xyz[1],staPos[2]+xyz[2]);

Is there any problem in my codes? Thank you.

Regards.

-- YanWei - 14 Jul 2009

Answer

If you answer a question - or have a question you asked answered by someone - please remember to edit the page and set the status to answered. The status is in a drop-down list below the edit box.

Hi Yan,

Again, I jump in ... This link FromGeodeticToTopocentricCoordinates might clear up some things for you. If so, please mark your question as solved/answered.

regards,

-- OctavianAndrei - 14 Jul 2009

Hi,Again,

The link is helpful. I make a mistake in my codes:

// rotate angle around U axes
uen.R1(angle);

Thank you.

wishes

-- YanWei - 14 Jul 2009