Examples of gaia.cu1.mdb.cu3.localplanecoords.dm.LpcCentroid

Package gaia.cu1.mdb.cu3.localplanecoords.dm

Examples of gaia.cu1.mdb.cu3.localplanecoords.dm.LpcCentroid

gaia.cu1.tools.numeric.algebra.GVector2d

    synchronized(gaia) {
      scanAngle = gaia.getAttitude().getScanAngle(comrs, time);
    }
    
    StellarAstrometry astro = (refAstro != null)? refAstro: astrometry;
    GVector2d lsc = Conversion.fromCoMRStoLSC(comrs, astro.getAlpha(), astro.getDelta(), scanAngle);
    
    // Return field angles
    return lsc;
    
  }

View Full Code Here

    ObservationalDataSet obsData = new ObservationalDataSet();
    
    for(Transit tr: model.getFilteredTransits()) {
      
      // Obtain field angles
      GVector2d fieldAngles = model.getFieldAngles(tr);
      if(!useAc)
        fieldAngles = new GVector2d(fieldAngles.getX(), 0.0);
      
      // Add noise if boolean switch is true 
      GVector2d errors = new GVector2d ( sigTraAl, sigTraAc );
      if(addNoise) {
        GVector2d randomOffset = new GVector2d(sigTraAl * random.nextGaussian(), useAc? sigTraAc * random.nextGaussian() : 0);
        fieldAngles.add(randomOffset);
      }
      
      // Calculate scanning angle
      double scanAngle;

View Full Code Here

      boolean relativity = false;
      GaiaSimuTime gst = new GaiaSimuTime(new GaiaTime());
      gst.setJD(timeJD);
      GVector3d comrs = model.getStar().getAstrometry().getCoMRSPosition(gst, relativity);
      GVector3d comrsRef = model.getAstrometry().getCoMRSPosition(gst, relativity);
      GVector2d lpc = Conversion.fromCoMRStoLPC(comrs, comrsRef.getLongitude(), comrsRef.getLatitude());
      
      orbit.addTransit ( new TransitData ( lpc, new GVector2d(), true, gst, 0.0, false ) );
      
    }
    
    // Return the orbit
    return orbit;

View Full Code Here

  private void plotTrajectory ( Writer out ) throws IOException, SimuException {
    
    for ( TransitData tr: getTransits() ) {
      
      double timeJD = tr.getTime().getJD();
      GVector2d lpc = tr.getLPC();
      
      // Assemble output string
      String str = String.format("%20.13e %20.13e %20.13e\n",
          timeJD,
          lpc.getX(), lpc.getY());
      
      // Write to writer
      out.write ( str );
      
    }

View Full Code Here

    // Calculate local celestial coordinates via gnonomic projection with reference triad (GAIA-LL-061 eq. 2)
    GVector3d[] triad = getReferenceTriad(refAlpha0, refDelta0);
    double  a = triad[0].dot(comrs) / triad[2].dot(comrs) * GaiaParam.Nature.RADIAN_MILLIARCSECOND,
        d = triad[1].dot(comrs) / triad[2].dot(comrs) * GaiaParam.Nature.RADIAN_MILLIARCSECOND;
    
    return new GVector2d(a,d);
    
  }

View Full Code Here

    
    // Get reference triad
    GVector3d[] triad = getReferenceTriad(refAlpha0, refDelta0);
    
    // Perform inverse transformation to CoMRS coordinates using the triad
    GVector2d lpcRad = GVector2d.scale(GaiaParam.Nature.MILLIARCSECOND_RADIAN, lpc);
    
    GVector3d pos = new GVector3d(triad[2]);
    pos.add(GVector3d.scale(lpcRad.getX(), triad[0]));
    pos.add(GVector3d.scale(lpcRad.getY(), triad[1]));
    pos.scale(1/Math.sqrt(1+lpcRad.normSquared()));
    
    return pos;
    
  }

View Full Code Here

    
    // Transform to local plane coordinates
    double st = Math.sin(scanAngle);
    double ct = Math.cos(scanAngle);
    
    return new GVector2d( lsc.getX()*st - lsc.getY()*ct, lsc.getX()*ct + lsc.getY()*st );
    
  }

View Full Code Here

    double dz2 = Math.pow ( lscError.getY(), 2 );
    double st2 = Math.pow ( Math.sin(scanAngle), 2 );
    double ct2 = Math.pow ( Math.cos(scanAngle), 2 );
    
    // Calculate and return LPC error
    return new GVector2d ( Math.sqrt(dw2*st2 + dz2*ct2), Math.sqrt(dw2*ct2 + dz2*st2) );
    
  }

View Full Code Here

    
    // Transform to local plane coordinates
    double  st = Math.sin(scanAngle),
        ct = Math.cos(scanAngle);
    
    return new GVector2d( lpc.getX()*st + lpc.getY()*ct, - lpc.getX()*ct + lpc.getY()*st );
    
  }

View Full Code Here

          
          for ( int i = 0; i < obsData.getTransits().size(); i++ ) {
            
            if ( fitModel.useAcrossScan() ) {
              
              GVector2d bestFitLSC = new GVector2d(bestFit[2*i], bestFit[2*i+1]);
              GVector2d targetsLSC = new GVector2d(targets[2*i], targets[2*i+1]);
              GVector2d errorsLSC = new GVector2d(errors[2*i], errors[2*i+1]);
              double scanAngle = obsData.getTransits().get(i).getScanAngle();
              GVector2d bestFitLPC = Conversion.fromLSCtoLPC(bestFitLSC, scanAngle);
              GVector2d targetsLPC = Conversion.fromLSCtoLPC(targetsLSC, scanAngle);
              GVector2d errorsLPC = Conversion.fromLSCtoLPCError(errorsLSC, scanAngle);
              
              buffwriter.write(String.format("%20.12e %20.12e %20.12e %20.12e %20.12e %20.12e %20.12e %20.12e %20.12e %20.12e %20.12e %20.12e %20.12e\n",
                  times[i],
                  targets[2*i], targets[2*i+1],
                  errors[2*i], errors[2*i+1],
                  bestFit[2*i], bestFit[2*i+1],
                  targetsLPC.getX(), targetsLPC.getY(),
                  errorsLPC.getX(), errorsLPC.getY(),
                  bestFitLPC.getX(), bestFitLPC.getY()));
              
            } else {
              
              buffwriter.write(String.format("%20.12e %20.12e %20.12e %20.12e\n",

View Full Code Here

0 1 2 3

TOP

Related Classes of gaia.cu1.mdb.cu3.localplanecoords.dm.LpcCentroid

newExamples.archive.EpochAstrometryConverter

newExamples.archive.GogGbinTest

newExamples.archive.GogGbinTestOld

newExamples.archive.GSInvScanLaw

newExamples.archive.SimulateOnePlanetManyInclinations

newExamples.Conversion

newExamples.DataGenerator

newExamples.FittingRun

newExamples.Model

newExamples.ModelFunction

All source code are property of their respective owners. Java is a trademark of Sun Microsystems, Inc and owned by ORACLE Inc. Contact coftware#gmail.com.