[编程交流] 网格

; angtodms       by ymg                                                       ;
; Given an angle in radian and number of decimal after the seconds            ;
; Returns a string formatted in Degree Minutes and Seconds,                   ;
; with proper symblol for degree and leading 0 on the minutes and seconds     ;
;           0°00'00.0"    271°05'06.3" etc.                                   ;
(defun angtodms (a prec /)
    (setq a (angtos a 1 (+ prec 4)))
    (while (<  (vl-string-position (ascii "d") a) 3)           
        (setq a (strcat " " a))
    )
    (if (< (vl-string-position (ascii "'") a) 6)
 (setq a (strcat (substr a 1 4) "0" (substr a 5)))
    )
    (if (< (vl-string-position (ascii """) a) (+ 9 (if (= 0 prec) 0 (1+ prec))))
 (setq a (strcat (substr a 1 7) "0" (substr a ))
    )
    (vl-string-subst  "°" "d" a) 
)
(defun dtr (a) (* pi (/ a 180.0)))
;; lltodms                                                                    ;
;;                                                                            ;
;; Arguments:  ll, A list (lat lon) in decimal degree                         ;
;;           prec, Number of decimals on the second                           ;
;; Returns:    A string "S 24°34'08.2042" - E 52°22'11.5320""               ;
;;                                                                            ;
;; Example:  (lltodms '(-24.5689456 52.36987) 4)                              ;
;; Needs subroutine angtodms and dtr                                          ;
(defun lltodms (ll prec /)
   (setq lat (angtodms (abs (dtr (car ll))) prec)
         lon (angtodms (abs (dtr (cadr ll))) prec)
   )
   (if (minusp (car ll))
      (setq lat (strcat "S" lat))
      (setq lat (strcat "N" lat))
   )
   (if (minusp (car ll))
      (setq lon (strcat "W" lon))
      (setq lon (strcat "E" lon))
   )
   (strcat lat " - " lon)
)   

;;****************************************************************************;
;; lltodms          by ymg                                                    ;
;;                                                                            ;
;; Arguments: lat   Latitude in decimal degree                                ;
;;            lon   Longitude in decimal degree                               ;
;;           prec   Number of decimals on the second                          ;
;; Returns:      A string with lat and long in deg min sec                      ;
;;                                                                            ;
;; Example:  (lltodms -24.5689456 -52.369875 4)                               ;
;; Result :  " 24°34'08.2042"S -  52°22'11.5500"W"                          ;
;;****************************************************************************;
(defun lltodms (lat lon prec)
  
  ; Degree to radians routine                                               ;
  (defun dtr (a) (* pi (/ a 180.0)))
  ; angtodms       by ymg                                                    ;
  ; Given an angle in radian and number of decimal after the seconds         ;
  ; Returns a string formatted in Degree Minutes and Seconds,                ;
  ; with proper symblol for degree and leading 0 on the minutes and seconds  ;
  ;           0°00'00.0"    271°05'06.3" etc.                                ;
  (defun angtodms (a prec)
     (setq a (angtos a 1 (+ prec 4)))
     (while (<  (vl-string-position (ascii "d") a) 3)           
        (setq a (strcat " " a))
     )
     (if (< (vl-string-position (ascii "'") a) 6)
 (setq a (strcat (substr a 1 4) "0" (substr a 5)))
     )
     (if (< (vl-string-position (ascii """) a) (+ 9 (if (= 0 prec) 0 (1+ prec))))
 (setq a (strcat (substr a 1 7) "0" (substr a ))
     )
     (vl-string-subst  "°" "d" a) 
   )
   ;-------------------------------------------------------------------------;
   (strcat 
        (strcat (angtodms (abs (dtr lat)) prec)(if (minusp lat) "S" "N"))
         " - "
        (strcat (angtodms (abs (dtr lon)) prec)(if (minusp lon) "W" "E"))
   )    
)