注册 登录  
 加关注
   显示下一条  |  关闭
温馨提示!由于新浪微博认证机制调整,您的新浪微博帐号绑定已过期,请重新绑定!立即重新绑定新浪微博》  |  关闭

面朝大海 春暖花开

一点心得,转载本博客文章,请注明转帖,谢谢!

 
 
 

日志

 
 
关于我

中国科学院博士 主要从事遥感机理、定量反演、数据处理以及GIS应用研究。ArcGIS、Envi 、ERDAS、Ecognition软件、IDL语言、6S、SAIL

网易考拉推荐

【转载】IDL代码共享-Sinusodial投影坐标转为经纬度坐标  

2014-10-24 19:57:56|  分类: idl |  标签: |举报 |字号 订阅

  下载LOFTER 我的照片书  |
;+
; :Name:
;    SNSOID_TO_LATLON
;
; :Description:
;    Convert coordinate from sinusodial to latitude/longitude.
;
; :Returns:
;
; :Params:
;    X  - Abscissa of the input point.
;    Y  - Ordinate of the input point.
;    LAT - A named variable that will contain the returned latitude in degree.
;    LON - A named variable that will contain the returned longitude in degree.
;    
; :Keywords:
;    R  - Radius of curvature in meridian. The default value is 6371007.181000 meters.
;
;    
; :Examples:
;
; :Author: Dabin Ji
; :Email: dabinj@gmail.com
; :Date: 2011-10-5
;-
PRO SNSOID_TO_LATLON, X, Y, LAT, LON, R = R
    COMPILE_OPT IDL2

    pi = 3.1415926D
    
    X = DOUBLE(X)
    Y = DOUBLE(Y)
    
    IF ~KEYWORD_SET(R) THEN BEGIN
        R = 6371007.181D    ;meters
    ENDIF ELSE BEGIN
        R = DOUBLE(R)
    ENDELSE

    lat = y / R;
    lon = x / (R * COS(lat))
    
    ;Convert unit of lat_lon from radian to degree.
    lat = lat * 180.0 / pi
    lon = lon * 180.0 / pi
END
  评论这张
 
阅读(217)| 评论(0)
推荐 转载

历史上的今天

评论

<#--最新日志,群博日志--> <#--推荐日志--> <#--引用记录--> <#--博主推荐--> <#--随机阅读--> <#--首页推荐--> <#--历史上的今天--> <#--被推荐日志--> <#--上一篇,下一篇--> <#-- 热度 --> <#-- 网易新闻广告 --> <#--右边模块结构--> <#--评论模块结构--> <#--引用模块结构--> <#--博主发起的投票-->
 
 
 
 
 
 
 
 
 
 
 
 
 
 

页脚

网易公司版权所有 ©1997-2017