Modification of GMT_utils.jl and test_GMT.jl

Author: adbayonao

ext/GMT_utils.jl

@@ -158,7 +158,13 @@ function import_GeoTIFF(fname::String; fieldname = :layer1, negative = false, is
         if length(size(G.image)) == 3
             data = permutedims(G.image, [2, 1, 3])
         elseif length(size(G.image)) == 2
-            data[:, :, 1] = G.image'
+            if size(G.image)==(nx,ny)
+              data[:,:,1] = G.image
+            elseif size(G.image)==(ny,nx)
+              data[:,:,1] = G.image'
+            else
+              error("unknown size; ")
+            end
         end
 
     end

test/test_GMT.jl

@@ -11,4 +11,4 @@ test_fwd = import_GeoTIFF("test_files/length_fwd.tif", fieldname = :forward)
 @test  maximum(test_fwd.fields.forward) ≈ 33.17775km
 
 test2 = import_GeoTIFF("test_files/UTM2GTIF.TIF")
-@test   test2.fields.layer1[20, 20] == 105.0
+@test   test2.fields.layer1[20, 20] == 233.0