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

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


    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

TOP

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

Copyright © 2018 www.massapicom. All rights reserved.
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.