Update README.md and LICENSE.md

Author: giordano

LICENSE.md

@@ -1,6 +1,6 @@
 The Reproject.jl package is licensed under the MIT "Expat" License:
 
-> Copyright (c) 2019 Mosè Giordano
+> Copyright (c) 2019 Mosè Giordano, Rohit Kumar
 >
 > Permission is hereby granted, free of charge, to any person obtaining a copy
 > of this software and associated documentation files (the "Software"), to deal

Project.toml

@@ -1,6 +1,6 @@
 name = "Reproject"
 uuid = "d1dcc2e6-806e-11e9-2897-3f99785db2ae"
-authors = ["Mosè Giordano"]
+authors = ["Mosè Giordano", "Rohit Kumar"]
 version = "0.1.0"
 
 [deps]

README.md

@@ -2,3 +2,14 @@
 
 [![Build Status](https://travis-ci.com/giordano/Reproject.jl.svg?branch=master)](https://travis-ci.com/giordano/Reproject.jl)
 [![Coveralls](https://coveralls.io/repos/github/giordano/Reproject.jl/badge.svg?branch=master)](https://coveralls.io/github/giordano/Reproject.jl?branch=master)
+
+Implementation in [Julia](https://julialang.org/) of the
+[`reproject`](https://github.com/astropy/reproject) package by Thomas
+Robitaille, part of the Astropy project.
+
+License
+-------
+
+The `reproject` package is released under the terms of the BSD 3-Clause "New" or
+"Revised" License.  The `Reproject.jl` package received written permission to be
+released under the MIT "Expat" License.

src/core.jl

@@ -22,25 +22,25 @@ function reproject(input_data, output_projection; shape_out = nothing, order::In
     else
         array_in, wcs_out = parse_input_data(input_data, hdu_in)
     end
-    
+
     if output_projection isa FITS || output_projection isa String
         wcs_in, shape_out = parse_output_projection(output_projection, hdu_out)
     else
         wcs_in, shape_out = parse_output_projection(output_projection, shape_out)
     end
-    
+
     type_in = wcs_to_celestial_frame(wcs_in)
     type_out = wcs_to_celestial_frame(wcs_out)
-    
+
     if type_in == type_out && shape_out === nothing
         return array_in
     end
-    
+
     img_out = fill(NaN, shape_out)
-    
+
     itp = interpolator(array_in, order)
     shape_in = size(array_in)
-    
+
     for i in 1:shape_out[1]
         for j in 1:shape_out[2]
             pix_coord_in = [float(i), float(j)]
@@ -68,12 +68,12 @@ function reproject(input_data, output_projection; shape_out = nothing, order::In
 
             pix_coord_out = world_to_pix(wcs_out, [rad2deg(lon(coord_out)), rad2deg(lat(coord_out))])
 
-            if 1 <= pix_coord_out[1] <= shape_in[1] && 1 <= pix_coord_out[2] <= shape_in[2] 
+            if 1 <= pix_coord_out[1] <= shape_in[1] && 1 <= pix_coord_out[2] <= shape_in[2]
                 img_out[i,j] = itp(pix_coord_out[1], pix_coord_out[2])
             end
         end
     end
-    
+
     return img_out, (!isnan).(img_out)
 end
 
@@ -91,6 +91,6 @@ function interpolator(array_in::AbstractArray, order::Int)
     else
         itp = interpolate(array_in, BSpline(Quadratic(InPlace(OnCell()))))
     end
-    
+
     return itp
 end

src/parsers.jl

@@ -7,8 +7,8 @@ Parse input data and returns an Array and WCS object.
 
 # Arguments
 - `input_data`: image to reproject which can be name of a FITS file,
-                an ImageHDU or a FITS file. 
-- `hdu_in`: used to set HDU to use when more than one HDU is present. 
+                an ImageHDU or a FITS file.
+- `hdu_in`: used to set HDU to use when more than one HDU is present.
 """
 function parse_input_data(input_data::ImageHDU)
     return read(input_data), WCS.from_header(read_header(input_data, String))[1]
@@ -35,7 +35,7 @@ end
 Parse output projection and returns a WCS object and shape of output.
 
 # Arguments
-- `output_projection`: WCS information about the image to be reprojected which can be 
+- `output_projection`: WCS information about the image to be reprojected which can be
                        name of a FITS file, an ImageHDU or WCSTransform.
 - `shape_out`: shape of the output image.
 - `hdu_number`: specifies HDU number when file name is given as input.
@@ -65,13 +65,12 @@ end
 
 function parse_output_projection(output_projection::FITS, hdu_number)
     wcs_out = WCS.from_header(read_header(output_projection[hdu_number], String))[1]
-    
+
     if output_projection[hdu_number] isa ImageHDU
         shape_out = size(output_projection[hdu_number])
-    else    
+    else
         throw(ArgumentError("Given FITS file doesn't have ImageHDU"))
     end
-    
+
     return wcs_out, shape_out
 end
-

src/utils.jl

@@ -6,17 +6,17 @@ The reference frame supported in Julia are FK5, ICRS and Galactic.
 """
 function wcs_to_celestial_frame(wcs::WCSTransform)
     radesys = wcs.radesys
-    
+
     xcoord = wcs.ctype[1][1:4]
     ycoord = wcs.ctype[2][1:4]
-    
+
     if radesys == ""
         if xcoord == "GLON" && ycoord == "GLAT"
             radesys = "Gal"
         elseif xcoord == "TLON" && ycoord == "TLAT"
             radesys = "ITRS"
         end
     end
-    
+
     return radesys
 end

test/core.jl

@@ -7,13 +7,13 @@
         mkpath("data")
         download("https://astropy.stsci.edu/data/galactic_center/gc_msx_e.fits", "data/gc_msx_e.fits")
     end
-    
+
     imgin = FITS("data/gc_msx_e.fits")         # project this
     imgout = FITS("data/gc_2mass_k.fits")        # into this coordinate
-    
+
     hdu1 = astropy.io.fits.open("data/gc_2mass_k.fits")[1]
     hdu2 = astropy.io.fits.open("data/gc_msx_e.fits")[1]
-    
+
     function check_diff(arr1, arr2)
         diff = 1e-3
         count = 0
@@ -26,14 +26,14 @@
         end
         return  iszero(count)
     end
-    
+
     @test check_diff(reproject(imgin, imgout, order = 0)[1], rp.reproject_interp(hdu2, hdu1.header, order = 0)[1]) == true
     @test check_diff(reproject(imgout, imgin, order = 0)[1], rp.reproject_interp(hdu1, hdu2.header, order = 0)[1]) == true
     @test check_diff(reproject(imgin, imgout, order = 1)[1], rp.reproject_interp(hdu2, hdu1.header, order = 1)[1]) == true
     @test check_diff(reproject(imgin, imgout, order = 2)[1], rp.reproject_interp(hdu2, hdu1.header, order = 2)[1]) == true
     @test check_diff(reproject(imgin[1], imgout[1], shape_out = (1000,1000))[1],
             rp.reproject_interp(hdu2, astropy.wcs.WCS(hdu1.header), shape_out = (1000,1000))[1]) == true
-    
+
     wcs = WCSTransform(2; ctype = ["RA---AIR", "DEC--AIR"], radesys = "UNK")
     @test_throws ArgumentError reproject(imgin, wcs, shape_out = (100,100))
 

test/parsers.jl

@@ -14,30 +14,30 @@ using Reproject: parse_input_data, parse_output_projection
 
     indata = reshape(Float32[1:100;], 5, 20)
     write(f, indata; header=inhdr)
-    
+
     @testset "ImageHDU type" begin
         result = parse_input_data(f[1])
         @test result[1] isa Array
         @test result[2] isa WCSTransform
     end
-    
+
     @testset "Single HDU FITS file" begin
         result = parse_input_data(f, 1)
         @test result[1] isa Array
         @test result[2] isa WCSTransform
     end
     close(f)
-    
+
     @testset "String filename input" begin
         result = parse_input_data(fname, 1)
         @test result[1] isa Array
         @test result[2] isa WCSTransform
     end
-    
+
     f = FITS(fname, "w")
     write(f, indata; header=inhdr)
     write(f, indata; header=inhdr)
-    
+
     @testset "Multiple HDU FITS file" begin
         result = parse_input_data(f, 2)
         @test result[1] isa Array
@@ -47,7 +47,7 @@ using Reproject: parse_input_data, parse_output_projection
         result = parse_input_data(fname, 1)
         @test result[1] isa Array
         @test result[2] isa WCSTransform
-    end   
+    end
 end
 
 @testset "output parser" begin
@@ -65,33 +65,32 @@ end
 
     indata = reshape(Float32[1:100;], 5, 20)
     write(f, indata; header=inhdr)
-    
+
     @testset "ImageHDU type" begin
         result = parse_output_projection(f[1], (12,12))
         @test result[1] isa WCSTransform
         @test result[2] isa Tuple
         @test_throws DomainError parse_output_projection(f[1], ())
     end
     close(f)
-    
+
     @testset "String filename" begin
         result = parse_output_projection(fname, 1)
         @test result[1] isa WCSTransform
         @test result[2] isa Tuple
     end
-    
+
     wcs = WCSTransform(2;
                           cdelt = [-0.066667, 0.066667],
                           ctype = ["RA---AIR", "DEC--AIR"],
                           crpix = [-234.75, 8.3393],
                           crval = [0., -90],
                           pv    = [(2, 1, 45.0)])
-    
+
     @testset "WCSTransform input" begin
         result = parse_output_projection(wcs, (12,12))
         @test result[1] isa WCSTransform
         @test result[2] isa Tuple
         @test_throws DomainError parse_output_projection(wcs, ())
     end
 end
-

test/utils.jl

@@ -5,21 +5,21 @@ using Reproject: wcs_to_celestial_frame
                        )
     wcs2 = WCSTransform(2;
                        ctype = ["RA---AIR", "DEC--AIR"],
-                       equinox = 1888.67 
+                       equinox = 1888.67
                        )
     wcs3 = WCSTransform(2;
                        ctype = ["RA---AIR", "DEC--AIR"],
-                       equinox = 2000 
+                       equinox = 2000
                        )
     wcs4 = WCSTransform(2;
-                       ctype = ["GLON--", "GLAT--"], 
+                       ctype = ["GLON--", "GLAT--"],
                        )
     wcs5 = WCSTransform(2;
                        ctype = ["TLON", "TLAT"],
                       )
     wcs6 = WCSTransform(2;
                        ctype = ["RA---AIR", "DEC--AIR"],
-                       radesys = "UNK" 
+                       radesys = "UNK"
                       )
 
     @test wcs_to_celestial_frame(wcs1) == "ICRS"