Examples of gov.fnal.eag.healpix.PixTools

gov.fnal.eag.healpix.PixTools
contains methods translated from HEALPix Fortran90 with increased map resolution in comparison to original Fortran code. @author N Kuropatkin

    System.out.println(" test query_disc is done");


  }


  public void testQuery_Triangle() throws Exception {
    PixTools pt = new PixTools();
    long nside = 4;
    int nest = 0;
    long ipix = 0;
    int[] pixel1 = { 57, 58, 59, 60, 61, 62, 74, 75, 76, 77, 78, 90, 91,
        92, 93, 107, 108, 109, 123, 124, 140 };
    int[] pixel2 = { 88, 89, 90, 91, 92, 105, 106, 107, 108, 121, 122, 123,
        138, 139, 154 };
    int[] pixel3 = { 49, 64, 65, 80, 81, 95, 96, 97, 111, 112, 113, 126,
        127, 128, 129, 142, 143, 144, 145 };
    int[] pixel4 = { 36, 52, 53, 67, 68, 69, 83, 84, 85, 86, 98, 99, 100,
        101, 102, 114, 115, 116, 117, 118, 119, 129, 130, 131, 132,
        133, 134, 135 };
    int[] pixel5 = { 58, 127, 56, 126, 125, 50, 123, 124, 119, 48, 122,
        121, 118, 117, 74, 175, 120, 115, 116, 191, 72, 174, 173, 114,
        113, 190, 189, 66 };
    int[] pixel6 = { 110, 123, 124, 125, 140, 141, 156 };
    int[] pixel7 = { 53, 68, 69 };
    long pix1 = 62;
    long pix2 = 57;
    long pix3 = 140;
    System.out.println("Start test Query Triangle !!!!!!!!!!!!!!!!!!!!");
    Vector3d v11 = pt.pix2vect_ring(nside, pix1);


    Vector3d v22 = pt.pix2vect_ring(nside, pix2);


    Vector3d v33 = pt.pix2vect_ring(nside, pix3);


    //    System.out.println("nside="+nside+" triangle pixels "+pix1+" "+pix2+"
    // "+pix3);
    int inclusive = 0;


    ArrayList pixlist;
    pixlist = pt.query_triangle(nside, v11, v22, v33, nest, inclusive);


    int nlist = pixlist.size();


    for (int i = 0; i < nlist; i++) {
      ipix = ((Long) pixlist.get(i)).longValue();


      assertEquals("pixel = " + ipix, pixel1[i], ipix, 1e-10);
    }
    pix1 = 92;
    pix2 = 88;
    pix3 = 154;
    v11 = pt.pix2vect_ring(nside, pix1);
    v22 = pt.pix2vect_ring(nside, pix2);
    v33 = pt.pix2vect_ring(nside, pix3);


    inclusive = 0;
    ArrayList pixlist1;
    pixlist1 = pt.query_triangle(nside, v11, v22, v33, nest, inclusive);


    nlist = pixlist1.size();


    for (int i = 0; i < nlist; i++) {
      ipix = ((Long) pixlist1.get(i)).longValue();
      assertEquals("pixel = " + ipix, pixel2[i], ipix, 1e-10);
    }
    pix1 = 49;
    pix2 = 142;
    pix3 = 145;
    v11 = pt.pix2vect_ring(nside, pix1);
    v22 = pt.pix2vect_ring(nside, pix2);
    v33 = pt.pix2vect_ring(nside, pix3);


    inclusive = 0;
    ArrayList pixlist2;
    pixlist2 = pt.query_triangle(nside, v11, v22, v33, nest, inclusive);


    nlist = pixlist2.size();


    for (int i = 0; i < nlist; i++) {
      ipix = ((Long) pixlist2.get(i)).longValue();
      assertEquals("pixel = " + ipix, pixel3[i], ipix, 1e-10);
    }
    pix1 = 36;
    pix2 = 129;
    pix3 = 135;
    v11 = pt.pix2vect_ring(nside, pix1);
    v22 = pt.pix2vect_ring(nside, pix2);
    v33 = pt.pix2vect_ring(nside, pix3);


    inclusive = 0;


    pixlist2 = pt.query_triangle(nside, v11, v22, v33, nest, inclusive);


    nlist = pixlist2.size();


    for (int i = 0; i < nlist; i++) {
      ipix = ((Long) pixlist2.get(i)).longValue();
      assertEquals("pixel = " + ipix, pixel4[i], ipix, 1e-10);
    }
    pix1 = 36;
    pix2 = 129;
    pix3 = 135;
    nest = 1;
    v11 = pt.pix2vect_ring(nside, pix1);
    v22 = pt.pix2vect_ring(nside, pix2);
    v33 = pt.pix2vect_ring(nside, pix3);
    inclusive = 0;


    pixlist2 = pt.query_triangle(nside, v11, v22, v33, nest, inclusive);


    nlist = pixlist2.size();


    for (int i = 0; i < nlist; i++) {
      ipix = ((Long) pixlist2.get(i)).longValue();
//            System.out.println("ipix="+ipix);
      assertEquals("pixel = " + ipix, pixel5[i], ipix, 1e-10);
    }
    pix1 = 123;
    pix2 = 156;
    pix3 = 110;
    nest = 0;
    v11 = pt.pix2vect_ring(nside, pix1);
    v22 = pt.pix2vect_ring(nside, pix2);
    v33 = pt.pix2vect_ring(nside, pix3);
    inclusive = 0;


    pixlist2 = pt.query_triangle(nside, v11, v22, v33, nest, inclusive);


    nlist = pixlist2.size();


    for (int i = 0; i < nlist; i++) {
      ipix = ((Long) pixlist2.get(i)).longValue();
      assertEquals("pixel = " + ipix, pixel6[i], ipix, 1e-10);
      //      System.out.println("i="+i+" pixel#="+ipix);
    }
    pix1 = 69;
    pix2 = 53;
    pix3 = 68;
    nest = 0;
    v11 = pt.pix2vect_ring(nside, pix1);
    v22 = pt.pix2vect_ring(nside, pix2);
    v33 = pt.pix2vect_ring(nside, pix3);
    inclusive = 0;


    pixlist2 = pt.query_triangle(nside, v11, v22, v33, nest, inclusive);


    nlist = pixlist2.size();


    for (int i = 0; i < nlist; i++) {
      ipix = ((Long) pixlist2.get(i)).longValue();

View Full Code Here

    System.out.println(" test query_triangle is done");


  }


  public void testQuery_Polygon() throws Exception {
    PixTools pt = new PixTools();
    long nside = 4;
    int nest = 0;
    long ipix = 0;
    int inclusive = 0;
    int[] result = { 51, 52, 53, 66, 67, 68, 69, 82, 83, 84, 85, 86, 98,
        99, 100, 101, 115, 116, 117 };
    int[] result1 = { 55, 70, 71, 87 };
    int[] result2 = { 137, 152, 153, 168 };
    int[] result3 = { 27, 43, 44, 58, 59, 60, 74, 75, 76, 77, 89, 90, 91,
        92, 93, 105, 106, 107, 108, 109, 110, 121, 122, 123, 124, 125,
        138, 139, 140, 141, 154, 156 };


    Vector3d[] vertices = new Vector3d[6];
    System.out.println("Start test query_polygon !!!!!!!!!!!!!!!!!!!!!!");
    ArrayList vlist = new ArrayList();
    Vector3d v = pt.pix2vect_ring(nside, 53);
    vlist.add((Object) v);
    v = pt.pix2vect_ring(nside, 51);
    vlist.add((Object) v);
    v = pt.pix2vect_ring(nside, 82);
    vlist.add((Object) v);
    v = pt.pix2vect_ring(nside, 115);
    vlist.add((Object) v);
    v = pt.pix2vect_ring(nside, 117);
    vlist.add((Object) v);
    v = pt.pix2vect_ring(nside, 86);
    vlist.add((Object) v);


    ArrayList pixlist;
    pixlist = pt.query_polygon(nside, vlist, nest, inclusive);
    //    System.out.println(" List size="+pixlist.size());
    int nlist = pixlist.size();
    //    System.out.println(" Pixel list:");
    for (int i = 0; i < nlist; i++) {
      ipix = ((Long) pixlist.get(i)).longValue();
      assertEquals("pixel = " + ipix, result[i], ipix, 1e-10);
      //      System.out.println("i="+i+" pixel # "+ipix);
    }


    /* Yet another test */


    ArrayList vlist1 = new ArrayList();
    v = pt.pix2vect_ring(nside, 71);
    vlist1.add((Object) v);
    v = pt.pix2vect_ring(nside, 55);
    vlist1.add((Object) v);
    v = pt.pix2vect_ring(nside, 70);
    vlist1.add((Object) v);
    v = pt.pix2vect_ring(nside, 87);
    vlist1.add((Object) v);
    pixlist = pt.query_polygon(nside, vlist1, nest, inclusive);


    nlist = pixlist.size();


    for (int i = 0; i < nlist; i++) {
      ipix = ((Long) pixlist.get(i)).longValue();
      //      System.out.println("i="+i+" pixel # "+ipix);
      assertEquals("pixel = " + ipix, result1[i], ipix, 1e-10);
    }


    /* Yet another test */
    ArrayList vlist2 = new ArrayList();
    v = pt.pix2vect_ring(nside, 153);
    vlist2.add((Object) v);
    v = pt.pix2vect_ring(nside, 137);
    vlist2.add((Object) v);
    v = pt.pix2vect_ring(nside, 152);
    vlist2.add((Object) v);
    v = pt.pix2vect_ring(nside, 168);
    vlist2.add((Object) v);
    pixlist = pt.query_polygon(nside, vlist2, nest, inclusive);


    nlist = pixlist.size();


    for (int i = 0; i < nlist; i++) {
      ipix = ((Long) pixlist.get(i)).longValue();
      assertEquals("pixel = " + ipix, result2[i], ipix, 1e-10);
      //      System.out.println("i="+i+" pixel # "+ipix);
    }
    /* Yet another test */


    ArrayList vlist3 = new ArrayList();
    v = pt.pix2vect_ring(nside, 110);
    vlist3.add((Object) v);
    v = pt.pix2vect_ring(nside, 27);
    vlist3.add((Object) v);
    v = pt.pix2vect_ring(nside, 105);
    vlist3.add((Object) v);
    v = pt.pix2vect_ring(nside, 154);
    vlist3.add((Object) v);
    v = pt.pix2vect_ring(nside, 123);
    vlist3.add((Object) v);
    v = pt.pix2vect_ring(nside, 156);
    vlist3.add((Object) v);
    pixlist = pt.query_polygon(nside, vlist3, nest, inclusive);


    nlist = pixlist.size();


    for (int i = 0; i < nlist; i++) {
      ipix = ((Long) pixlist.get(i)).longValue();

View Full Code Here

  }
    public void testMaxResolution() {
        System.out.println(" Start test MaxRes !!!!!!!!!!!!!!!!!!!!!");


        long nside = 1048576;
        PixTools pt = new PixTools();
        double res = pt.PixRes(nside);
        System.out.println("Minimum size of the pixel side is "+res+" arcsec.");
        assertEquals("res = " + res, 0.2, res, 1e-1);
        System.out.println(" End of MaxRes test _______________________");
    }

View Full Code Here

        assertEquals("res = " + res, 0.2, res, 1e-1);
        System.out.println(" End of MaxRes test _______________________");
    }
    public void testQueryDiskRes() throws Exception {
        System.out.println(" Start test DiskRes !!!!!!!!!!!!!!!!!!!!!");
        PixTools pt = new PixTools();


        int nest = 0;
        long ipix = 0;
 
        int inclusive = 0;
        double theta= Math.PI;
        double phi = Math.PI;
        double radius = Math.toRadians(0.2/3600.); //  One arcse
        long nside = pt.GetNSide(radius);
        System.out.println(" calculated nside="+nside);
        long cpix = pt.ang2pix_ring(nside,theta,phi);
        Vector3d vc = pt.pix2vect_ring(nside, cpix);
        ArrayList pixlist;
        pixlist = pt.query_disc(nside, vc, radius, nest, inclusive);


        int nlist = pixlist.size();
        for (int i = 0; i < nlist; i++) {
            ipix = ((Long) pixlist.get(i)).longValue();
            Vector3d v = pt.pix2vect_ring(nside,ipix);
            double dist = pt.AngDist(v,vc);
            assertTrue(dist<=2.*radius);
        }
        cpix = pt.ang2pix_nest(nside,theta,phi);
        Vector3d vc1 = pt.pix2vect_nest(nside, cpix);
        long cpixnest= pt.vect2pix_nest(nside,vc1);
        ArrayList pixlist1;
        nest = 1;
        radius *=4;
        pixlist1 = pt.query_disc(nside, vc1, radius, nest, inclusive);
        int nlist1 = pixlist1.size();
        for (int i = 0; i < nlist1; i++) {
            ipix = ((Long) pixlist1.get(i)).longValue();
            Vector3d v = pt.pix2vect_nest(nside,ipix);
            double dist = pt.AngDist(v,vc1);
            assertTrue(dist<=2.*radius);
        }
        System.out.println(" test query disc Hi Res is done -------------------"); 
    }

View Full Code Here

    }
    public void testGetNside() {
        System.out.println(" Start test GetNside !!!!!!!!!!!!!!!!!!!!!");


        double pixsize = 0.3;
        PixTools pt = new PixTools();
        long nside = pt.GetNSide(pixsize);
        System.out.println("Requared nside is "+nside);
        assertEquals("nside = " + nside, 1048576, nside, 1e-1);
        System.out.println(" End of GetNSide test _______________________");
    }

View Full Code Here

0 1 2

TOP

Related Classes of gov.fnal.eag.healpix.PixTools

gov.fnal.eag.healpix.test.PixToolsTest

java.util.Vector

java.util.ArrayList

javax.vecmath.Vector3d

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.