Debugging finding the bounds of Tilted Regions

GeoRegionsEcosystem · Jul 26, 2024 · e567197 · e567197
1 parent c39a4a8
commit e567197
Show file tree

Hide file tree

Showing 2 changed files with 3 additions and 5 deletions.
diff --git a/src/georegions/show.jl b/src/georegions/show.jl
@@ -20,6 +20,7 @@ function show(io::IO, geo::TiltRegion)
 		"    Position    (X,Y) : ",[geo.X,geo.Y], '\n',
 		"    Size      (ΔX,ΔY) : ",[geo.ΔX,geo.ΔY], '\n',
 		"    Tilt          (θ) : ",[geo.θ], '\n',
+		"    Bounds  (N,S,E,W) : ",[geo.N,geo.S,geo.E,geo.W], '\n',
         "        (is180,is360) : ",(geo.is180,geo.is360),"\n",
 	)
 end

diff --git a/src/georegions/tilt.jl b/src/georegions/tilt.jl
@@ -9,12 +9,9 @@ function getTiltBounds(
     θ  :: Real;
 )
 
-    W = X - ΔX * cosd(θ)
-    E = X + ΔX * cosd(θ)
-    S = Y - ΔY * sind(θ)
-    N = Y + ΔY * sind(θ)
+    lon,lat = getTiltShape(X,Y,ΔX,ΔY,θ)
 
-    return N,S,E,W
+    return maximum(lat), minimum(lat) , maximum(lon), minimum(lon)
 
 end