Changeset 2894 in ntrip

Timestamp:

Jan 26, 2011, 4:55:35 PM (13 years ago)

Author:

mervart

Message:

 

Location:

trunk/BNC

Files:

• bncantex.cpp (modified) (5 diffs)
• bncantex.h (modified) (1 diff)
• bncmodel.cpp (modified) (2 diffs)
• bncmodel.h (modified) (1 diff)

trunk/BNC/bncantex.cpp

@@ -61 +61 @@
 }
 
-// 
+// Print 
 ////////////////////////////////////////////////////////////////////////////
 void bncAntex::print() const {
@@ -86 +86 @@
 }
 
-// 
+// Read ANTEX File
 ////////////////////////////////////////////////////////////////////////////
 t_irc bncAntex::readFile(const QString& fileName) {
@@ -182 +182 @@
           QTextStream inLine(&line, QIODevice::ReadOnly);
           inLine >> newFrqMap->neu[0] >> newFrqMap->neu[1] >> newFrqMap->neu[2];
+          newFrqMap->neu[0] *= 1e-3;
+          newFrqMap->neu[1] *= 1e-3;
+          newFrqMap->neu[2] *= 1e-3;
         }
         else if (line.indexOf("NOAZI") == 3) {
@@ -192 +195 @@
             inLine >> newFrqMap->pattern[ii];
           }
+          newFrqMap->pattern *= 1e-3;
         }
       }
@@ -199 +203 @@
   return success;
 }
+
+// Phase Center Offset (Receiver Antenna and GPS only)
+////////////////////////////////////////////////////////////////////////////
+double bncAntex::pco(const QString& antName, double eleSat) {
+
+  static const double f1 = t_CST::freq1;
+  static const double f2 = t_CST::freq2;
+  static const double c1 =   f1 * f1 / (f1 * f1 - f2 * f2);
+  static const double c2 = - f2 * f2 / (f1 * f1 - f2 * f2);
+
+  QMap<QString, t_antMap*>::const_iterator it = _maps.find(antName);
+  if (it != _maps.end()) {
+    t_antMap* map = it.value();
+    if (map->frqMapL1 && map->frqMapL2) {
+      double corr1 = -map->frqMapL1->neu[2] * sin(eleSat);
+      double corr2 = -map->frqMapL2->neu[2] * sin(eleSat);
+      return c1 * corr1 + c2 * corr2;
+    }
+  }
+
+  return 0.0;
+}

trunk/BNC/bncantex.h

@@ -37 +37 @@
   t_irc readFile(const QString& fileName);  
   void print() const;
+  double pco(const QString& antName, double eleSat);
 
  private:

trunk/BNC/bncmodel.cpp

@@ -267 +267 @@
       _antex = 0;
     }
+    _antennaName = settings.value("pppAntenna").toString();
   }
 
@@ -386 +387 @@
   }
 
-  return satData->rho + clk() + offset - satData->clk + tropDelay + wind;
+  double phaseCenter = 0.0;
+  if (_antex) {
+    phaseCenter = _antex->pco(_antennaName, satData->eleSat);
+  }
+
+  return satData->rho + phaseCenter + clk() 
+                      + offset - satData->clk + tropDelay + wind;
 }
 

trunk/BNC/bncmodel.h

@@ -155 +155 @@
   double                _sigClk0;
   bncAntex*             _antex;
+  QString               _antennaName;
 };