- Timestamp:
- Jul 21, 2023, 7:25:42 AM (16 months ago)
- Location:
- trunk/BNC/src
- Files:
-
- 3 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/BNC/src/bncantex.cpp
r9675 r10130 118 118 } 119 119 120 // Print 121 //////////////////////////////////////////////////////////////////////////// 122 QString bncAntex::snxCodeSinexString(const std::string& antName) { 123 124 if (antName.find("NULLANTENNA") != string::npos) { 125 return QString(" ----------"); 126 } 127 128 QString antNameQ = antName.c_str(); 129 if (_maps.find(antNameQ) == _maps.end()) { 130 return QString(" ----------"); 131 } 132 else { 133 return QString(" %1").arg(_maps[antNameQ]->snxCode, 10, QLatin1Char(' ')); 134 } 135 } 136 120 137 // Read ANTEX File 121 138 //////////////////////////////////////////////////////////////////////////// … … 182 199 inLine >> newAntMap->zen1 >> newAntMap->zen2 >> newAntMap->dZen; 183 200 } 201 else if (line.indexOf("SINEX CODE") == 60) { 202 QTextStream inLine(&line, QIODevice::ReadOnly); 203 inLine >> newAntMap->snxCode; 204 } 184 205 185 206 // Start of Frequency -
trunk/BNC/src/bncantex.h
r9587 r10130 41 41 void print() const; 42 42 QString pcoSinexString(const std::string& antName, t_frequency::type frqType); 43 QString snxCodeSinexString(const std::string& antName); 43 44 double satCorr(const QString& prn, t_frequency::type frqType, 44 45 double eleSat, double azSat, bool& found) const; … … 63 64 public: 64 65 t_antMap() { 65 zen1 = 0.0; 66 zen2 = 0.0; 67 dZen = 0.0; 66 zen1 = 0.0; 67 zen2 = 0.0; 68 dZen = 0.0; 69 snxCode = ""; 68 70 } 69 71 ~t_antMap() { … … 78 80 double zen2; 79 81 double dZen; 82 QString snxCode; 80 83 QMap<t_frequency::type, t_frqMap*> frqMap; 81 84 bncTime validFrom; -
trunk/BNC/src/bncsinextro.cpp
r10127 r10130 32 32 33 33 _opt = opt; 34 _sampl = 34 _sampl = sampl; 35 35 36 36 _antex = 0; … … 54 54 int daysec = int(fmod(GPSWeeks, 86400.0)); 55 55 int dayOfYear = datTim.date().dayOfYear(); 56 QString yy = datTim.toString("yy");57 QString creationTime = QString("%1:%2:%3").arg(yy )56 QString yyyy = datTim.toString("yyyy"); 57 QString creationTime = QString("%1:%2:%3").arg(yyyy) 58 58 .arg(dayOfYear, 3, 10, QLatin1Char('0')) 59 59 .arg(daysec , 5, 10, QLatin1Char('0')); … … 75 75 int nominalStartSec = daysec - (int(fmod(double(daysec), double(intr)))); 76 76 int nominalEndSec = nominalStartSec + intr - _sampl; 77 QString endTime = QString("%1:%2:%3").arg(yy )77 QString endTime = QString("%1:%2:%3").arg(yyyy) 78 78 .arg(dayOfYear , 3, 10, QLatin1Char('0')) 79 79 .arg(nominalEndSec , 5, 10, QLatin1Char('0')); … … 89 89 corr = settings.value("PPP/corrFile").toString(); 90 90 } 91 92 91 QString signalPriorities = QString::fromStdString(_opt->_signalPriorities); 93 92 if (!signalPriorities.size()) { … … 117 116 << creationTime.toStdString() << " " << ac.toStdString() << " " 118 117 << startTime.toStdString() << " " << endTime.toStdString() << " P " 119 << epo.toStdString() << " 0 " << " T "<< endl;118 << _opt->_roverName.substr(0,9) << endl; 120 119 _out << "*-------------------------------------------------------------------------------" << endl; 121 120 _out << "+FILE/REFERENCE" << endl; … … 124 123 _out << " OUTPUT " << "Total Troposphere Zenith Path Delay Product" << endl; 125 124 _out << " SOFTWARE " << BNCPGMNAME << endl; 126 _out << " INPUT " << "OBS: " << _opt->_roverName.substr(0,9) << ", SSR: " 127 _out << " VERSON NUMBER " << QString("%1").arg(solId, 10, QLatin1Char('')).toStdString() << endl;125 _out << " INPUT " << "OBS: " << _opt->_roverName.substr(0,9) << ", SSR: " << corr.toStdString() << endl; 126 _out << " VERSON NUMBER " << QString("%1").arg(solId, 3, QLatin1Char('0')).toStdString() << endl; 128 127 _out << "-FILE/REFERENCE" << endl; 129 128 _out << "*-------------------------------------------------------------------------------" << endl; 130 129 131 130 QString systems; 131 if (settings.value("PPP/lcGPS" ).toString() != "no") {systems += "G, ";} 132 if (settings.value("PPP/lcGLONASS").toString() != "no") {systems += "R, ";} 133 if (settings.value("PPP/lcGalileo").toString() != "no") {systems += "E, ";} 134 if (settings.value("PPP/lcBDS" ).toString() != "no") {systems += "C";} 132 135 QString blqFileName = QString::fromStdString(_opt->_blqFileName); 133 136 QString blqFileBaseName; … … 141 144 } 142 145 } 146 bool eleWeighting = false; 147 if (settings.value("PPP/eleWgtCode").toInt() > 0 || 148 settings.value("PPP/eleWgtCode").toInt() > 0 ) { 149 eleWeighting = true; 150 } 151 143 152 _out << "+TROP/DESCRIPTION" << endl; 144 153 _out << "*KEYWORD______________________ VALUE(S)______________" << endl; 145 _out << " SAMPLING INTERVAL " 146 << setw(4) << _sampl << endl; 147 _out << " SAMPLING TROP " 148 << setw(4) << _sampl << endl; 149 _out << " ELEVATION CUTOFF ANGLE " 150 << setw(4) << int(_opt->_minEle * 180.0/M_PI) << endl; 154 _out << " TROPO SAMPLING INTERVAL " << setw(22) << 1 << endl; 155 _out << " DATA SAMPLING TROP " << setw(22) << _sampl << endl; 156 _out << " ELEVATION CUTOFF ANGLE " << setw(22) << int(_opt->_minEle * 180.0/M_PI) << endl; 157 _out << " OBSERVATION WEIGHTING " << ((eleWeighting) ? "SINEL" : "") << endl; 158 _out << " GNSS SYSTEMS " << systems.toStdString() << endl; 159 _out << " TIME SYSTEM " << "G" << endl; 160 _out << " TROPO MODELING METHOD " << "KALMAN FILTER" << endl; 151 161 if (! blqFileName.isEmpty()) { 152 162 _out << " OCEAN TIDE LOADING MODEL " << blqFileBaseName.toStdString() + blqFileExtension.toStdString() << endl; 153 163 } 154 164 _out << " TROP MAPPING FUNCTION " << "Saastamoinen" << endl; 155 _out << " SOLUTION_FIELDS_1 " << "TROTOT STDEV" << endl; 165 _out << " TROPO PARAMETER NAMES " << "TROTOT STDEV" << endl; 166 _out << " TROPO PARAMETER UNITS " << "1e+03 1e+03" << endl; 156 167 _out << "-TROP/DESCRIPTION"<< endl; 157 168 _out << "*-------------------------------------------------------------------------------" << endl; … … 185 196 186 197 _out << "+TROP/COORDINATES" << endl; 187 _out << "*S ITE PT SOLN T STA_X_______ STA_Y_______ STA_Z_______ SYSTEM REMARK" << endl;198 _out << "*STATION__ PT SOLN T __DATA_START__ __DATA_END____ __STA_X_____ __STA_Y_____ __STA_Z_____ SYSTEM REMRK" << endl; 188 199 _out << " " << _opt->_roverName.substr(0,9) << " A " 189 << solId.toStdString() << " P" 200 << QString("%1").arg(solId, 4, QLatin1Char(' ')).toStdString() << " P " 201 << startTime.toStdString() << " " << endTime.toStdString() 190 202 << QString(" %1").arg(_opt->_xyzAprRover(1), 12, 'f', 3, QLatin1Char(' ')).toStdString() 191 203 << QString(" %1").arg(_opt->_xyzAprRover(2), 12, 'f', 3, QLatin1Char(' ')).toStdString() 192 204 << QString(" %1").arg(_opt->_xyzAprRover(3), 12, 'f', 3, QLatin1Char(' ')).toStdString() 193 << " IGS20 " << ac.toStdString() << endl; 205 << QString(" %1").arg("IGS20", 6, QLatin1Char(' ')).toStdString() 206 << QString(" %1").arg(ac, 5, QLatin1Char(' ')).toStdString() << endl; 194 207 _out << "-TROP/COORDINATES"<< endl; 195 208 _out << "*-------------------------------------------------------------------------------" << endl; … … 197 210 198 211 _out << "+SITE/ECCENTRICITY" << endl; 199 _out << "* UP______ NORTH___ EAST____" << endl;200 _out << "*S ITE PT SOLN T DATA_START__ DATA_END____ AXE ARP->BENCHMARK(M)_________" << endl;212 _out << "* UP______ NORTH___ EAST____" << endl; 213 _out << "*STATION__ PT SOLN T __DATA_START__ __DATA_END____ AXE MARKER->ARP(m)__________" << endl; 201 214 _out << " " << _opt->_roverName.substr(0,9) << " A " 202 << solId.toStdString() << " P " 203 << startTime.toStdString() << " " << endTime.toStdString() << " UNE" 204 << QString("%1").arg(_opt->_neuEccRover(3), 9, 'f', 4, QLatin1Char(' ')).toStdString() 205 << QString("%1").arg(_opt->_neuEccRover(1), 9, 'f', 4, QLatin1Char(' ')).toStdString() 206 << QString("%1").arg(_opt->_neuEccRover(2), 9, 'f', 4, QLatin1Char(' ')).toStdString() << endl; 215 << QString("%1").arg(solId, 4, QLatin1Char(' ')).toStdString() << " P " 216 << startTime.toStdString() << " " << endTime.toStdString() 217 << QString(" %1").arg("UNE", 3, QLatin1Char(' ')).toStdString() 218 << QString("%1").arg(_opt->_neuEccRover(3), 8, 'f', 4, QLatin1Char(' ')).toStdString() 219 << QString("%1").arg(_opt->_neuEccRover(1), 8, 'f', 4, QLatin1Char(' ')).toStdString() 220 << QString("%1").arg(_opt->_neuEccRover(2), 8, 'f', 4, QLatin1Char(' ')).toStdString() << endl; 207 221 _out << "-SITE/ECCENTRICITY" << endl; 208 222 _out << "*-------------------------------------------------------------------------------" << endl; … … 211 225 _out << "+SITE/RECEIVER" << endl; 212 226 _out << "*STATION__ PT SOLN T __DATA_START__ __DATA_END____ DESCRIPTION_________ S/N_________________ FIRMW______" << endl; 213 _out << " " << _opt->_roverName.substr(0,9) << " A " << solId.toStdString() << " P " 227 _out << " " << _opt->_roverName.substr(0,9) << " A " 228 << QString("%1").arg(solId, 4, QLatin1Char(' ')).toStdString() << " P " 214 229 << startTime.toStdString() << " " << endTime.toStdString() 215 230 << QString(" %1").arg(_opt->_recNameRover.c_str(), 20,QLatin1Char(' ')).toStdString() 216 << " ----- " << " -----------" << endl;231 << " --------------------" << " -----------" << endl; 217 232 _out << "-SITE/RECEIVER" << endl; 218 233 _out << "*-------------------------------------------------------------------------------" << endl; 219 234 } 220 221 _out << "+SITE/ANTENNA" << endl;222 _out << "*STATION__ PT SOLN T __DATA_START__ __DATA_END____ DESCRIPTION_________ S/N_________________ PCV_MODEL_" << endl;223 _out << " " << _opt->_roverName.substr(0,9) << " A " << solId.toStdString() << " P "224 << startTime.toStdString() << " " << endTime.toStdString() << " "225 << _opt->_antNameRover << " -----" << endl;226 _out << "-SITE/ANTENNA" << endl;227 _out << "*-------------------------------------------------------------------------------" << endl;228 229 235 if (!_opt->_antexFileName.empty()) { 230 236 _antex = new bncAntex(_opt->_antexFileName.c_str()); 231 if (_opt->_LCsGPS.size()) { 232 _out << "+SITE/GPS_PHASE_CENTER" << endl; 233 _out << "* UP____ NORTH_ EAST__ UP____ NORTH_ EAST__" << endl; 234 _out << "*DESCRIPTION_________ S/N__ L1->ARP(m)__________ L2->ARP(m)__________ AZ_EL____" << endl; 235 _out << QString(" %1").arg(_opt->_antNameRover.c_str(), 20,QLatin1Char(' ')).toStdString() 236 << " -----"; 237 for (int i = 0; i < frqStrList.size(); ++i) { 238 if (frqStrList.at(i).contains('G')) {//cout << frqStrList.at(i).toStdString() << endl; 239 _out << _antex->pcoSinexString(_opt->_antNameRover, t_frequency::toInt(frqStrList.at(i).toStdString())).toStdString(); 240 } 241 } 242 _out << " ---------" << endl; 243 _out << "-SITE/GPS_PHASE_CENTER" << endl; 244 _out << "*-------------------------------------------------------------------------------" << endl; 245 } 246 if (_opt->_LCsGLONASS.size()) { 247 _out << "+SITE/GLONASS_PHASE_CENTER" << endl; 248 _out << "* UP____ NORTH_ EAST__ UP____ NORTH_ EAST__" << endl; 249 _out << "*DESCRIPTION_________ S/N__ L1->ARP(m)__________ L2->ARP(m)__________ AZ_EL____" << endl; 250 _out << QString(" %1").arg(_opt->_antNameRover.c_str(), 20,QLatin1Char(' ')).toStdString() 251 << " -----"; 252 for (int i = 0; i < frqStrList.size(); ++i) { 253 if (frqStrList.at(i).contains('R')) {//cout << frqStrList.at(i).toStdString() << endl; 254 _out << _antex->pcoSinexString(_opt->_antNameRover, t_frequency::toInt(frqStrList.at(i).toStdString())).toStdString(); 255 } 256 } 257 _out << " ---------" << endl; 258 _out << "-SITE/GLONASS_PHASE_CENTER" << endl; 259 _out << "*-------------------------------------------------------------------------------" << endl; 260 } 261 if (_opt->_LCsGalileo.size()) { 262 _out << "+SITE/GALILEO_PHASE_CENTER" << endl; 263 _out << "* UP____ NORTH_ EAST__ UP____ NORTH_ EAST__" << endl; 264 _out << "*DESCRIPTION_________ S/N__ L1->ARP(m)__________ L2->ARP(m)__________ AZ_EL____" << endl; 265 _out << QString(" %1").arg(_opt->_antNameRover.c_str(), 20,QLatin1Char(' ')).toStdString() 266 << " -----"; 267 for (int i = 0; i < frqStrList.size(); ++i) { 268 if (frqStrList.at(i).contains('E')) {//cout << frqStrList.at(i).toStdString() << endl; 269 _out << _antex->pcoSinexString(_opt->_antNameRover, t_frequency::toInt(frqStrList.at(i).toStdString())).toStdString(); 270 } 271 } 272 _out << " ---------" << endl; 273 _out << "-SITE/GALILEO_PHASE_CENTER" << endl; 274 _out << "*-------------------------------------------------------------------------------" << endl; 275 } 276 if (_opt->_LCsBDS.size()) { 277 _out << "+SITE/BEIDOU_PHASE_CENTER" << endl; 278 _out << "* UP____ NORTH_ EAST__ UP____ NORTH_ EAST__" << endl; 279 _out << "*DESCRIPTION_________ S/N__ L1->ARP(m)__________ L2->ARP(m)__________ AZ_EL____" << endl; 280 _out << QString(" %1").arg(_opt->_antNameRover.c_str(), 20,QLatin1Char(' ')).toStdString() 281 << " -----"; 282 for (int i = 0; i < frqStrList.size(); ++i) { 283 if (frqStrList.at(i).contains('C')) {//cout << frqStrList.at(i).toStdString() << endl; 284 _out << _antex->pcoSinexString(_opt->_antNameRover, t_frequency::toInt(frqStrList.at(i).toStdString())).toStdString(); 285 } 286 } 287 _out << " ---------" << endl; 288 _out << "-SITE/BEIDOU_PHASE_CENTER" << endl; 289 _out << "*-------------------------------------------------------------------------------" << endl; 290 } 237 _out << "+SITE/ANTENNA" << endl; 238 _out << "*STATION__ PT SOLN T __DATA_START__ __DATA_END____ DESCRIPTION_________ S/N_________________ PCV_MODEL_" << endl; 239 _out << " " << _opt->_roverName.substr(0,9) << " A " 240 << QString("%1").arg(solId, 4, QLatin1Char(' ')).toStdString() << " P " 241 << startTime.toStdString() << " " << endTime.toStdString() 242 << QString(" %1").arg(_opt->_antNameRover.c_str(), 20,QLatin1Char(' ')).toStdString() 243 << " --------------------" << _antex->snxCodeSinexString(_opt->_antNameRover).toStdString() << endl; 244 _out << "-SITE/ANTENNA" << endl; 245 _out << "*-------------------------------------------------------------------------------" << endl; 291 246 delete _antex; 292 247 _antex = 0; 293 } 294 248 } 295 249 _out << "+TROP/SOLUTION" << endl; 296 250 _out << "*STATION__ ____EPOCH_____ TROTOT STDDEV " << endl;
Note:
See TracChangeset
for help on using the changeset viewer.