링크해 주신 글은 저도 읽어 보니, David Fanning 박사님 글 답지 않게 전개가 복잡하네요.
원래 이 글의 내용은 위경도 자료로 된 데이터를 25km x 25km 의 북극 Polar Stereographic 자료로 변환하는 과정을 보여주기 위한 것입니다. 원하시는 작업과 반대인 것 같습니다.
그렇지만, 맨 뒤쪽에 Faster Gridding 챕터를 보면, 본래 목적이 위경도 영상을 Polar Stereographic으로 변환하기 위한 작업임에도 불구하고, Polar Stereographic 의 x, y 좌표를 위경도로 바꾸는 과정이 나옵니다. 스마트하게 문제를 풀기위해서 좌표를 역으로 접근하다 보니 그렇게 하는 것인데, 이 부분이 도움이 되실 것 같습니다.
David Fanning 박사님의 Coyote 라이브러리가 필요하고, cgScaleVector 함수는 현재 Scale_Vector 함수로 대체되었습니다. 이 부분은 수정하셔야 할 것 같습니다.