What's the math behind this?

trandoductin · **Posted:** Mon Apr 18, 2016 10:49 pm **(#1)**

look at this picture

I am trying to write a script that asks user for N sections they want to have and then divide
a square (or rectangle) image into that many NxN squares (or rectangles).

I thought if i loop from 1 to N for Y and then inside that loop i loop from 1 to N for X
selecting the square or rectangle..
using math
I can perspective it into the image shown above...
but i am stuck at the math...
so basically I got cos(180 degrees) down to cos(0 degrees) ..i put degrees here but i actually convert it to radians before calling cos...
But I am still stumpled at the math to get the mapping of these squares/rectangles into the shape shown above.
i am assuming if i have the math down, a rectangle image will turn to an oval shape and a square will turn into a circle...
anyone can offer any help on the math behind this?

Odinbc · **Posted:** Tue Apr 19, 2016 1:20 pm **(#2)**

Wolfram MathWorld is worth a look.

ofnuts · **Joined:** Oct 25, 2010 **Posts:** 4812

trandoductin wrote:

look at this picture
[ Image ]

I am trying to write a script that asks user for N sections they want to have and then divide
a square (or rectangle) image into that many NxN squares (or rectangles).

I thought if i loop from 1 to N for Y and then inside that loop i loop from 1 to N for X
selecting the square or rectangle..
using math
I can perspective it into the image shown above...
but i am stuck at the math...
so basically I got cos(180 degrees) down to cos(0 degrees) ..i put degrees here but i actually convert it to radians before calling cos...
But I am still stumpled at the math to get the mapping of these squares/rectangles into the shape shown above.
i am assuming if i have the math down, a rectangle image will turn to an oval shape and a square will turn into a circle...
anyone can offer any help on the math behind this?

Assuming center at 0,0:

Code: Select all

a=atan2(y,x) # not the same as atan(y/x)!
X=x*cos(a)
Y=y*sin(a)

paynekj · **Posted:** Tue Apr 19, 2016 3:24 pm **(#4)**

I can create pretty much the same thing using the plug-in-lens-distortion run a few times on a grid:

Attachment:

Untitled.png [ 168.07 KiB | Viewed 2677 times ]

Code: Select all

(define (script-fu-kp_lens aimg adraw)
  (let* (
          (space 32)
          (thick 2)
        )
   (gimp-ellipse-select aimg 0 0 (car (gimp-image-width aimg)) (car (gimp-image-height aimg)) CHANNEL-OP-REPLACE TRUE FALSE 0.0)
   (plug-in-grid RUN-NONINTERACTIVE aimg adraw thick space 0 '(0 0 0) 255 thick space 0 '(0 0 0) 255 0 2 6 '(0 0 0) 255)
   (plug-in-lens-distortion RUN-NONINTERACTIVE aimg adraw 0.0 0.0 100.0 100.0 0.0 0.0)
   (plug-in-lens-distortion RUN-NONINTERACTIVE aimg adraw 0.0 0.0 100.0 100.0 0.0 0.0)
   (plug-in-lens-distortion RUN-NONINTERACTIVE aimg adraw 0.0 0.0 100.0 100.0 0.0 0.0)
   (plug-in-lens-distortion RUN-NONINTERACTIVE aimg adraw 0.0 0.0 100.0 100.0 0.0 0.0)
   (plug-in-lens-distortion RUN-NONINTERACTIVE aimg adraw 0.0 0.0 100.0 100.0 0.0 0.0)
   (plug-in-lens-distortion RUN-NONINTERACTIVE aimg adraw 0.0 0.0 100.0 100.0 0.0 0.0)

; flush the image
   (gimp-displays-flush)
  )
)

Kevin

trandoductin · **Posted:** Tue Apr 19, 2016 3:57 pm **(#5)**

ofnuts wrote:

trandoductin wrote:

look at this picture
[ Image ]

I am trying to write a script that asks user for N sections they want to have and then divide
a square (or rectangle) image into that many NxN squares (or rectangles).

I thought if i loop from 1 to N for Y and then inside that loop i loop from 1 to N for X
selecting the square or rectangle..
using math
I can perspective it into the image shown above...
but i am stuck at the math...
so basically I got cos(180 degrees) down to cos(0 degrees) ..i put degrees here but i actually convert it to radians before calling cos...
But I am still stumpled at the math to get the mapping of these squares/rectangles into the shape shown above.
i am assuming if i have the math down, a rectangle image will turn to an oval shape and a square will turn into a circle...
anyone can offer any help on the math behind this?

Assuming center at 0,0:

Code: Select all

a=atan2(y,x) # not the same as atan(y/x)!
X=x*cos(a)
Y=y*sin(a)

I tried this formula and found an atan2 defined somewhere else on the web...and got this

Code: Select all

; script: SphereFold.scm
; author: Tin Tran
; date: April 18, 2016
; experimental
(define (atan2 y x)
    (if (> x 0)
      (atan (/ y x))
      (if (< x 0)
        (if (< y 0)
          (- (atan (/ y x)) *pi*)
          (+ (atan (/ y x)) *pi*) )
        (cond ;; x is zero
          ((> y 0) (/ *pi* 2))
          ((< y 0) (- (/ *pi* 2)))
          (else ;; x==y==0 is typically undefined but we return 0 instead
            0 ) ) ) ) )
(define (script-fu-sphere-fold image layer 
            sections         
         )
   
      (let* 
         (
         (copied-layer 0)
         (width (car (gimp-image-width image)))
         (height (car (gimp-image-height image)))
         (y-section 0)
         (x-section 0)
         (section-angle (/ 180 sections))
         (section-width (/ width sections))
         (section-height (/ height sections))
         (x-radius (/ width 2))
         (y-radius (/ height 2))
         (x-angle-start 0)
         (y-angle-start 0)
         (x-angle-end 0)
         (y-angle-end 0)
         (x-select-start 0)
         (y-select-start 0)
         (x-select-end 0)
         (y-select-end 0)
         (transformed-item 0)
         (Opposite 0)
         (Adjacent 0)
         (top-left-angle 0)
         (top-right-angle 0)
         (bottom-left-angle 0)
         (bottom-right-angle 0)
         )
         ;(gimp-image-undo-disable image); DN = NO UNDO
         
         (gimp-image-undo-group-start image)                   ;undo-group in one step
         ;(gimp-image-resize image (* width 2) (* height 2) (/ width 2) (/ height 2))
         ;(gimp-layer-resize-to-image-size layer)
         (set! y-angle-start 180)
         (set! y-select-start 0)
         (while (< y-section sections)
            
            (set! y-angle-end (- y-angle-start section-angle))
            (set! y-select-end (+ y-select-start section-height))
            
             (set! x-section 0)
            (set! x-angle-start 180)
            (set! x-select-start 0)
            (while (< x-section sections)
               (set! x-angle-end (- x-angle-start section-angle))
               (set! x-select-end (+ x-select-start section-width))
               
               (gimp-image-select-rectangle image CHANNEL-OP-REPLACE x-select-start y-select-start section-width section-height)
               
               ; top-left angle
               (set! Adjacent (- x-select-start x-radius)) 
               (set! Opposite (- y-select-start y-radius))
               (set! top-left-angle (atan2 Opposite Adjacent))
               ; top-right angle
               (set! Adjacent (- x-select-end x-radius)) 
               (set! Opposite (- y-select-start y-radius))
               (set! top-right-angle (atan2 Opposite Adjacent))
               
               ; bottom-left-angle
               (set! Adjacent (- x-select-start x-radius)) 
               (set! Opposite (- y-select-end y-radius))
               (set! bottom-left-angle (atan2 Opposite Adjacent))
               ; bottom-right-angle
               (set! Adjacent (- x-select-end x-radius)) 
               (set! Opposite (- y-select-end y-radius))
               (set! bottom-right-angle (atan2 Opposite Adjacent))
               
               (set! transformed-item (car (gimp-item-transform-perspective layer
                                               (- x-radius (*(cos top-left-angle) (abs (- x-radius x-select-start)) ))
                                       (- y-radius (*(sin top-left-angle) (abs (- y-radius y-select-start)) ))
                                       
                                       (- x-radius (*(cos top-right-angle) (abs (- x-radius x-select-end)) ))
                                       (- y-radius (*(sin top-right-angle) (abs (- y-radius y-select-start)) ))
                                       
                                       (- x-radius (*(cos bottom-left-angle) (abs (- x-radius x-select-start)) ))
                                       (- y-radius (*(sin bottom-left-angle) (abs (- y-radius y-select-end)) ))
                                       
                                       (- x-radius (*(cos bottom-right-angle) (abs (- x-radius x-select-end)) ))
                                       (- y-radius (*(sin bottom-right-angle) (abs (- y-radius y-select-end)) ))
                                       
                                    )
                                 )
               )
               
               (gimp-floating-sel-to-layer transformed-item)
               
               (set! x-angle-start (- x-angle-start section-angle))
               (set! x-select-start (+ x-select-start section-width))
               (set! x-section (+ x-section 1))
            )
            (set! y-angle-start (- y-angle-start section-angle))
            (set! y-select-start (+ y-select-start section-height))
            (set! y-section (+ y-section 1))
         )
         
         ;(gimp-image-undo-enable image) ;DN = NO UNDO
         (gimp-image-undo-group-end image)                     ;undo group in one step
         (gimp-displays-flush)
       )
   
   
    
) ;end of define
(script-fu-register
  "script-fu-sphere-fold"         ;function name
  "<Image>/Script-Fu/Create New/Sphere Fold ..."    ;menu register
  "Creates Map Fold effect"       ;description
  "Tin Tran"                          ;author name
  "copyright info and description"         ;copyright info or description
  "2016"                          ;date
  "RGB*, GRAY*"                        ;mode
  SF-IMAGE      "Image" 0                   
  SF-DRAWABLE   "Layer" 0
  SF-ADJUSTMENT "Fold sections" '(3 2 30 1 5 0 0)
)

any idea why it looks like it's the inverted sphere and divided into 4 quarters..instead of what i want?

I thought maybe there's something wrong with my scheme script and thought i'd try a test on codeskulptor with the same formula and I am not getting a quarter of a circle like i was hoping codeskulptor page that uses formula

ofnuts · **Joined:** Oct 25, 2010 **Posts:** 4812

trandoductin wrote:

any idea why it looks like it's the inverted sphere and divided into 4 quarters..instead of what i want?

I thought maybe there's something wrong with my scheme script and thought i'd try a test on codeskulptor with the same formula and I am not getting a quarter of a circle like i was hoping codeskulptor page that uses formula

Absolutely no idea. To figure that out I would add print statements (or gimp messages) at a few places in the code to see what the values are and how often some code blocks are executed.

PhotoMaster · **Posted:** Sat Jun 04, 2016 2:24 pm **(#7)**

I would look at Filters>Map>Map Object
Select Map to Sphere.

Check out viewtopic.php?f=10&t=430 or
viewtopic.php?f=10&t=429 for some ideas.
These are two of my tutorials here on GC that I think will help you.