Changeset 6878 in ntrip


Ignore:
Timestamp:
Jun 3, 2015, 5:29:19 PM (5 years ago)
Author:
stuerze
Message:

SSR VTEC: consideration of the fact that max degree and order have to be included

Location:
trunk/BNC/src
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • trunk/BNC/src/RTCM3/RTCM3coDecoder.cpp

    r6857 r6878  
    463463      t_vTecLayer layer;
    464464      layer._height = ionoLayer.Height;
    465       layer._C.ReSize(ionoLayer.Degree, ionoLayer.Order);
    466       layer._S.ReSize(ionoLayer.Degree, ionoLayer.Order);
    467       for (unsigned iDeg = 0; iDeg < ionoLayer.Degree; iDeg++) {
    468         for (unsigned iOrd = 0; iOrd < ionoLayer.Order; iOrd++) {
     465      layer._C.ReSize(ionoLayer.Degree+1, ionoLayer.Order+1);
     466      layer._S.ReSize(ionoLayer.Degree+1, ionoLayer.Order+1);
     467      for (unsigned iDeg = 0; iDeg <= ionoLayer.Degree; iDeg++) {
     468        for (unsigned iOrd = 0; iOrd <= ionoLayer.Order; iOrd++) {
    469469          layer._C[iDeg][iOrd] = ionoLayer.Cosinus[iDeg][iOrd];
    470470          layer._S[iDeg][iOrd] = ionoLayer.Sinus[iDeg][iOrd];
  • trunk/BNC/src/satObs.cpp

    r6859 r6878  
    299299    const t_vTecLayer& layer = vTec._layers[ii];
    300300    *out << setw(2)  << ii+1 << ' '
    301          << setw(2)  << layer._C.Nrows() << ' '   
    302          << setw(2)  << layer._C.Ncols() << ' '   
     301         << setw(2)  << layer._C.Nrows()-1 << ' '
     302         << setw(2)  << layer._C.Ncols()-1 << ' '
    303303         << setw(10) << setprecision(1) << layer._height << endl 
    304304         << setw(10) << setprecision(4) << layer._C
     
    334334    in >> dummy >> maxDeg >> maxOrd >> layer._height;
    335335
    336     layer._C.ReSize(maxDeg, maxOrd);
    337     layer._S.ReSize(maxDeg, maxOrd);
    338 
    339     for (int iDeg = 0; iDeg < maxDeg; iDeg++) {
    340       for (int iOrd = 0; iOrd < maxOrd; iOrd++) {
     336    layer._C.ReSize(maxDeg+1, maxOrd+1);
     337    layer._S.ReSize(maxDeg+1, maxOrd+1);
     338
     339    for (int iDeg = 0; iDeg <= maxDeg; iDeg++) {
     340      for (int iOrd = 0; iOrd <= maxOrd; iOrd++) {
    341341        inStream >> layer._C[iDeg][iOrd];
    342342      }
    343343    }
    344     for (int iDeg = 0; iDeg < maxDeg; iDeg++) {
    345       for (int iOrd = 0; iOrd < maxOrd; iOrd++) {
     344    for (int iDeg = 0; iDeg <= maxDeg; iDeg++) {
     345      for (int iOrd = 0; iOrd <= maxOrd; iOrd++) {
    346346        inStream >> layer._S[iDeg][iOrd];
    347347      }
Note: See TracChangeset for help on using the changeset viewer.