科学网—标签

相关帖子	版块	作者	回复/查看	最后发表

zhangdong 2014-7-19 00:22

在用 GPS 坐标计算空间两点的距离时，发现 R 语言中数字精度有问题。根据经纬度计算距离的公式在这里： http://www.cnblogs.com/ycsfwhh/archive/2010/12/20/1911232.html 简单来说，首先计算 C ，然后计算 Distance C =sin(LatA)*sin(LatB) + cos(LatA)*cos(LatB)*cos(MLonA-MLonB) Distance= R*Arccos(C)*Pi/180 同一个点到它自身的距离肯定是 0 ，这时 C=1 ，但是 R 语言中，这时候的 C 并不真正等于 1( 事实上是大于 1 的 ) ，这样在求 arccos(C) 的时候就会产生 na 我采用的解决的方法是 ifelse(dist_c= 1, dist_c - 1, dist_c - dist_c) 即对于大于等于 1 的 c ，直接设为 1 ，否则保持不变。参考： http://stackoverflow.com/questions/14026297/acos1-returns-nan-for-some-values-not-others 这里提到用 pmin() ，似乎没有解决我的问题

9213 次阅读|0 个评论

R语言根据经纬度获得海拔数据

Bearjazz 2014-4-29 08:21

#R 语言根据经纬度获得海拔数据 # 编者信息熊荣川六盘水师范学院生物信息学实验室 xiongrongchuan@126.com http://blog.sciencenet.cn/u/Bearjazz Google earth 上可以很容易的获得经纬度信息，但是海拔高度却不能直接提供，但据说其数据库中是有相应的海拔高度的数据的，下面是 R 语言爱好者使用 R 语言根据经纬度查询 google 数据库从而获取经纬度的代码：（原文地址 http://stackoverflow.com/questions/21593868/extracting-elevation-from-website-for-lat-lon-points-in-australia-using-r ） # 预装函数 googEl - function(locs) { require(RJSONIO) locstring - paste(do.call(paste, list(locs , locs , sep=',')), collapse='|') u - sprintf('http://maps.googleapis.com/maps/api/elevation/json?locations=%ssensor=false', locstring) res - fromJSON(u) out - t(sapply(res ], function(x) { c(x ] , x ] , x , x ) })) rownames(out) - rownames(locs) return(out) } # 使用方法数据的输入格式如下图，输入为表格（或矩阵） mat = read.excel(input.xlsx, 详表 ) # 表格输入 mat = googEl(mat) # 查询 google earth 数据库 lat lng elevation resolution 1 25.49222 108.1841 749.4119 152.7032 2 29.09253 105.4095 329.5924 152.7032 3 17.259 101.5063 1244.191 152.7032 4 18.94167 104.8097 287.4106 152.7032 5 21.03333 104.5 691.2963 152.7032 writeWorksheetToFile(output.xlsx, data=mat, sheet= 带海拔 ) # 输出结果

个人分类: 我的研究|11937 次阅读|0 个评论

谷歌地球的经纬度准确吗

热度 2 smalljasmine 2014-3-8 22:00

下了个landsat8 oli 影像，本来想用google earth来精校正的，但结果发现同样一个点，影像与谷歌地球差10'的纬度，整整差了6540米啊，又听说google earth是经过加密的，真不知道如何是好

4742 次阅读|2 个评论

[转载][转]在IDL中北京54(BJ54)转经纬度

fanxing018 2013-11-8 15:39

其实程序很简单。程序是Onion写的，直接上代码。 Pro convert_beijing54 ; +++++++++++++++++++++++++++++++++++++++++++++++++++ ; 注意 : ; 网上样例数据 J01,1,3124173.019,39498554.582 ; 前面的是y坐标(7位);后面的是x坐标(8位,前两位是带号) ;---------------------------------------------------- ixmap = 498554.582d ;注意此处去掉带号 iymap = 3124173.019d params = Dblarr(15) params = units = ENVI_TRANSLATE_PROJECTION_UNITS('Meters') iproj = ENVI_PROJ_CREATE(DATUM='Beijing-54',params=params,TYPE=3,Name='Gauss Kruger', UNITS=units) oproj = ENVI_PROJ_CREATE(/geographic) ENVI_CONVERT_PROJECTION_COORDINATES, ixmap, iymap, iproj, oxmap, oymap, oproj print,oxmap, oymap End 其实里面有一点需要注意：北京54坐标中，横坐标的前面两位是分带号需要剔除，这样才能转换成功，这个问题折腾了我很久。本例中横坐标是39498554.582 其中开头的39是带号，3度分带，第39带，所以params中中心经度是117

个人分类: ENVI/IDL学习|2158 次阅读|0 个评论

[转载]如何批量导入经纬度到谷歌地球中

baibing 2013-8-7 20:53

以下载个Global Mapper，将你需要导入的数据做个.csv格式的文件，然后先导入到Global Mapper ，在Global Mapper 中输出kml格式的文件，用google earth就可以打开了。 Global Mapper系列教程：批量输入点功能详解 http://www.godeyes.cn/html/2011/12/31/google_earth_12667.html

个人分类: 学习随笔|0 个评论

[转载]Surfer 上标下标经纬度符号

chinesezhao 2013-6-19 14:04

Surfer在Text里可以用特殊符号表示上标、下标和经度度符号等。例如： kg/m\up50 3 显示出来就是kg/m^3，其中3显示为上标。 co\dn50 2 显示出来就是co2，其中2显示为下标。 Longitude ({\fsymbol \a176 }E) 显示为Longitude（oE)，其中o显示为度符号。

0 个评论

netcdf-第二章

cwjwang 2012-10-3 21:39

2维数维：可以是代表真实物理意义的维，例如，时间，经纬度，长度，高度，也可以指示其他的性质如，状态和模型运行数据。netcdf维既有名字又有长度。维的长度任意的正的整数，除了64-bit offset 的维有UNLIMITED 长度。不受限制的长度索引号如一个记录号在传统的定向记录的文件里面。netcdf的经典和64-bit 数据库可以至少有一个不受限制的维，但是可以不需要有。如果一个变量具有一个不限制的维，其维必须是最有意义的一个（变化最慢的一个），以至于任何的无限制的维在CDL形状里面和c矩阵的申明里面必须是第一维。一个netcdf-4的数据库可以有多个不受限制的维，可以无序排列。增长变量随不限的维数的变化，写数据的时候用任何的netcdf的写的函数，可以突出其无限制的维的索引号为自己想要的记录号。 cdl维声明需要紧跟在cdl的关键的维数后面，多维的声明写在一行的时候，需要用逗号隔开，每一个声明都有一个统一的名=长度.用 /包含组的信息（netcdf-4输出文件需要）。有四种维在上面的例子里面，lat，lon，level，time。前面3种是固定长度的，time是不限的长度，即不受限制的维。netcdf数据库基本的单位是变量，一旦定义好变量，其形状就被指定为一系列的维，这些维必须已经存在的，维的数叫秩。标量的秩为0，向量的秩为1，矩阵的秩为2

3160 次阅读|0 个评论

[转载]GMT:Example B16—【xyz格式到grd网格数据】

zhenghui2915 2012-7-13 23:07

C:\GMT_Samples\ex06 首先到 http://topex.ucsd.edu/cgi-bin/get_data.cgi 下载要绘图范围的1-minutegridsinASCIIXYZ-format地形数据，输入经纬度范围E117-125，N19-28，文件保存为 Terr.txt gmtset GRID_CROSS_SIZE 0 ANOT_FONT_SIZE 15 FRAME_WIDTH 0.2c gmtset GRID_PEN 0.25pta gmtset PLOT_DEGREE_FORMAT ddd:mm:ssF gmtset FRAME_WIDTH 0.1c gmtset LABEL_FONT_SIZE 12 ANOT_FONT 1 LABEL_FONT 1 gmtdefaults -D .gmtdefaults4 gmtset BASEMAP_TYPE Plain set filename=gmt2011_01a.ps rem 由xyz格式的数据Terr.txt生成grd网格数据tw.grd surface Terr.txt -R117.0083/125.0083/19.003/28.0051 -I1m -Gtw.grd grdgradient tw.grd -A0/45 -Gtw.grd.int -Ne1 -M grd2cpt tw.grd -Csealand -L4000/8000 -S-8000/4000/500 topo_sealand.cpt grdimage tw.grd -Jm3 -R118.5/123.5/20/26.5 -Ctopo_sealand.cpt -Itw.grd.int -P -K -V -U"xuxiaoqiang" %filename% pscoast -R -Jm -Di -Ba2f1 -W0.8p,0/0/0 -P -O -V %filename% del .gmt* del *.int del *.cpt del *.grd pause http://xxqhome.blog.163.com/blog/static/196733020201192181356/

个人分类: GMT|6642 次阅读|0 个评论

[转载]NSDTF-DEM 文件头格式说明

wgghope 2012-6-19 23:00

转自 http://stdsoft.blogbus.com/logs/35976872.html 样例： NSDTF-DEM 1.0 M 0.000000 0.000000 39512435.000000 2743120.000000 5.000000 5.000000 985 1328 100 解释： DataMark--------中国地球空间数据交换格式-格网数据交换格式(CNSDTF-RAS或CNSDTF-DEM)的标志。基本部分，不可缺省。 Version--------该空间数据交换格式的版本号,如1.0。基本部分，不可缺省。 Unit--------坐标单位,K表示公里,M表示米,D表示以度为单位的经纬度,S表示以度分秒表示的经纬度(此时坐标格式为DDDMMSS.SSSS, DDD为度, MM为分, SS.SSSS为秒)。基本部分，不可缺省。 Alpha--------方向角。基本部分，不可缺省。 Compress--------压缩方法。0表示不压缩,1表示游程编码。基本部分，不可缺省。 Xo--------左上角原点X坐标。基本部分，不可缺省。 Yo--------左上角原点Y坐标。基本部分，不可缺省。 DX--------X方向的间距。基本部分，不可缺省。 DY--------Y方向的间距。基本部分，不可缺省。 Row--------行数。基本部分，不可缺省。 Col--------列数。基本部分，不可缺省。 HZoom--------高程放大倍率。基本部分，不可缺省。设置高程的放大倍率,使高程数据可以整数存贮,如高程精度精确到厘米,高程的放大倍率为100。如果不是DEM则HZoom为1。

3213 次阅读|0 个评论

地理坐标和经纬度坐标显示的切换

majiaping 2012-6-5 22:45

具体做法是：layers--propertiess--Data Frame Properties--General--Display，然后选择自己需要的显示方式移动鼠标可以右下角经纬度变化，见图

个人分类: 自学成才|4900 次阅读|0 个评论

建议国家海洋局用经纬度标注钓鱼岛及部分附属岛屿位置

热度 1 HeYujian 2012-3-3 18:17

刚读到一个消息：“ 国家海洋局公布钓鱼岛及部分附属岛屿标准名称”（ http://news.sohu.com/20120303/n336540813.shtml ）。发现国家海洋局对钓鱼岛及部分附属岛屿的位置用的是左右前后、东西南北的较模糊表述，连我这学化学出身的也感觉一点不专业--为什么不用明确的经纬度地理标注方式呢？查了一下该新闻的读者评论，已有18543人参与，基本上从政治与爱国角度评论。其中有一个评论与我想法一样：“ 搜狐手机网友：把经度纬度标上。来自 UC 浏览器”。从表述明确与淮确的角度出发，建议国家海洋局将钓鱼岛及部分附属岛屿位置描述补充加上经纬度，这无疑会更准确和符合国际惯例？更让小日本没话说！附录：国家海洋局公布钓鱼岛及部分附属岛屿标准名称我来说两句 (18543 人参与 ) 2012 年 03 月 03 日 01:29 来源：中国新闻网　　中新网北京 3 月 3 日电 ( 记者刘辰瑶 ) 国家海洋局 3 日凌晨发布消息，公布了中国钓鱼岛及其部分附属岛屿标准名称。　　根据《中华人民共和国海岛保护法》，国家海洋局对中国海域海岛进行了名称标准化处理。经国务院批准，国家海洋局、民政部现将钓鱼岛及其部分附属岛屿的标准名称、汉语拼音、位置描述公布如下：　　序号标准名称汉语拼音位置描述 1 钓鱼岛 Diàoyú Dǎo 距温州市约 356 千米、福州市约 385 千米、基隆市约 190 千米 2 龙头鱼岛 Lóngtóuyú Dǎo 位于钓鱼岛东北 3 鲳鱼岛 Chāngyú Dǎo 位于钓鱼岛西南 4 大黄鱼岛 Dàhuángyú Dǎo 位于钓鱼岛南 5 小黄鱼岛 Xiǎohuángyú Dǎo 位于钓鱼岛南 6 金钱鱼岛 Jīnqiányú Dǎo 位于钓鱼岛东南 7 金钱鱼西岛 Jīnqiányúxī Dǎo 位于钓鱼岛东南 8 梅童鱼岛 Méitóngyú Dǎo 位于钓鱼岛东南 9 梅童鱼东岛 Méitóngyúdōng Dǎo 位于钓鱼岛东南 10 梅童鱼西岛 Méitóngyúxī Dǎo 位于钓鱼岛东南 11 龙王鲷岛 Lóngwángdiāo Dǎo 位于钓鱼岛东南 12 龙王鲷西岛 Lóngwángdiāoxī Dǎo 位于钓鱼岛东南 13 龙王鲷东岛 Lóngwángdiāodōng Dǎo 位于钓鱼岛东南 14 龙王鲷南岛 Lóngwángdiāonán Dǎo 位于钓鱼岛东南 15 黄姑鱼岛 Huánggūyú Dǎo 位于钓鱼岛东南 16 黄尾屿 Huángwěi Yǔ 位于钓鱼岛东北约 27 千米处 17 海豚岛 Hǎitún Dǎo 位于黄尾屿西北 18 大珠岛 Dàzhū Dǎo 位于黄尾屿西 19 小珠岛 Xiǎozhū Dǎo 位于黄尾屿西 20 上虎牙岛 Shànghǔyá Dǎo 位于黄尾屿北 21 下虎牙岛 Xiàhǔyá Dǎo 位于黄尾屿北 22 西牛角岛 Xīniújiǎo Dǎo 位于黄尾屿东北 23 东牛角岛 Dōngniújiǎo Dǎo 位于黄尾屿东北 24 黄牛岛 Huángniú Dǎo 位于黄尾屿东北 25 牛尾岛 Niúwěi Dǎo 位于黄尾屿东北 26 牛蹄岛 Niútí Dǎo 位于黄尾屿东北 27 小龙岛 Xiǎolóng Dǎo 位于黄尾屿西 28 大雁岛 Dàyàn Dǎo 位于黄尾屿西 29 燕子岛 Yànzi Dǎo 位于黄尾屿西 30 刺猬岛 Cìwèi Dǎo 位于黄尾屿西南 31 卧蚕岛 Wòcán Dǎo 位于黄尾屿西南 32 大金龟子岛 Dàjīngūizǐ Dǎo 位于黄尾屿西南 33 小金龟子岛 Xiǎojīngūizǐ Dǎo 位于黄尾屿西南 34 海龟岛 Hǎigūi Dǎo 位于黄尾屿西南 35 海星岛 Hǎixīng Dǎo 位于黄尾屿东 36 海贝岛 Hǎibèi Dǎo 位于黄尾屿东南 37 赤尾屿 Chìwěi Yǔ 位于钓鱼岛东约 110 千米处 38 赤背北岛 Chìbèiběi Dǎo 位于赤尾屿北 39 赤背东岛 Chìbèidōng Dǎo 位于赤尾屿北 40 赤背西岛 Chìbèixī Dǎo 位于赤尾屿北 41 赤背南岛 Chìbèinán Dǎo 位于赤尾屿北 42 小赤尾岛 Xiǎochìwěi Dǎo 位于赤尾屿西 43 赤头岛 Chìtóu Dǎo 位于赤尾屿西 44 赤冠岛 Chìguàn Dǎo 位于赤尾屿西 45 赤鼻岛 Chìbí Dǎo 位于赤尾屿西 46 赤嘴岛 Chìzuǐ Dǎo 位于赤尾屿西 47 望赤岛 Wàngchì Dǎo 位于赤尾屿西南 48 北小岛 Běixiǎo Dǎo 位于钓鱼岛以东约 5 千米处 49 鸟巢岛 Niǎocháo Dǎo 位于北小岛东 50 鸟卵岛 Niǎoluǎn Dǎo 位于北小岛东 51 小鸟岛 Xiǎoniǎo Dǎo 位于北小岛东南 52 南小岛 Nánxiǎo Dǎo 位于钓鱼岛东南约 5.5 千米处 53 龙门北岛 Lóngménběi Dǎo 位于南小岛西北 54 龙门岛 Lóngmén Dǎo 位于南小岛西北 55 龙门南岛 Lóngménnán Dǎo 位于南小岛西北 56 卧龙岛 Wòlóng Dǎo 位于南小岛西北 57 卧龙西岛 Wòlóngxī Dǎo 位于南小岛西北 58 飞龙北岛 Fēilóngběi Dǎo 位于南小岛东南 59 飞龙岛 Fēilóng Dǎo 位于南小岛东南 60 龙珠岛 Lóngzhū Dǎo 位于南小岛东南 61 飞龙南岛 Fēilóngnán Dǎo 位于南小岛东南 62 长龙岛 Chánglóng Dǎo 位于南小岛东南 63 金龙岛 Jīnlóng Dǎo 位于南小岛东南 64 北屿 Běi Yǔ 位于钓鱼岛东北约 6 千米处 65 北屿仔岛 Běiyǔzǎi Dǎo 位于北屿南 66 小元宝岛 Xiǎoyuánbǎo Dǎo 位于北屿西南 67 飞云岛 Fēiyún Dǎo 位于北屿西南 68 元宝岛 Yuánbǎo Dǎo 位于北屿西南 69 南屿 Nán Yǔ 位于钓鱼岛东北约 7.4 千米处 70 飞屿 Fēi Yǔ 位于钓鱼岛东南 71 飞仔岛 Fēizǎi Dǎo 位于钓鱼岛东南　　作者：刘辰瑶 ( 责任编辑： UN608)

个人分类: 随思随想|5659 次阅读|0 个评论

[转载]Excel中将地理坐标度分秒格式的经纬度转换为十进制的度

热度 1 Deyan 2011-12-11 16:10

数据格式1101120 =VALUE(MID(IF(LEN(A1)=6,CONCATENATE("0",A1),A1),1,3))+VALUE(MID(IF(LEN(A1)=6,CONCATENATE("0",A1),A1),4,2))/60+VALUE(MID(IF(LEN(A1)=6,CONCATENATE("0",A1),A1),6,2))/3600

1085 次阅读|1 个评论

[杂言]：谁在对濒危物种落井下石？

Myrmecia 2011-12-10 03:40

提起某一物种濒危，原因大体都能归咎于自然环境的变迁、人类对其生存环境的破坏及对物种本身的直接伤害。为了挽救濒危物种，专家们进行野外调查，新闻媒体也会宣传引导。但如果处理不当，也会好心办坏事，反而会给濒危物种的生存状态带来更大更直接的威胁。第一种情况下，专家进行野外调查研究，发表学术论文，经常会给出确切的分布点（比如经纬度数据），这就给不法分子提供了便利。最近在留意中国特有的中华虎凤蝶（ Luehdorfia chinensis ），其中一篇文章（《为中华虎凤蝶呼救》）提到一位蝴蝶专家的“忏悔”： ------------------------------------ 一篇论文，引来一群疯狂的蝶商　　吴琦，从20世纪60年代就对蝴蝶入了迷。80年代初，中科院一位教授写信给吴琦，“早在30年代，有个美国人曾在南京龙潭捉到一只中华虎凤蝶，你也去找找看。”吴琦跑遍南京的大小山头，最先在牛首山，然后在紫金山、老山找到中华虎凤蝶。这是南京人第一次在南京发现中华虎凤蝶。　　兴奋的吴琦随即写了论文，详细介绍中华虎凤蝶在南京分布情况。不料，此文一发表，竟引来国内标本商人的兴趣，“他们纷纷来信，要我帮他们去抓中华虎凤蝶做标本。”吴琦赶忙向南京市科技局反映，要求政府出面保护中华虎凤蝶。“ ..... ..... 各路蝶商在南京活动了10多年，南京的中华虎凤蝶一年年在减少。张松奎颇为伤感地告诉记者，“二十几年前，我到牛首山东边调查，很小一块地方，就有一两千只中华虎凤蝶。整个南京，估计当时靠近一万只。现在呢？连十分之一也不剩了！”而发现中华虎凤蝶的吴琦，则有些后悔：保护不如保密...... ------------------------------------- 记得上大学的时候，野外实习带队老师还给我们讲过一件事。大概好像是说N年前，陕西秦岭分布着大量野生大鲵（又名”娃娃鱼“，世界上现存最大的也是最珍贵的两栖动物，国家二类保护水生野生动物）。后来，某年某月的某一日，专家们去搞调查研究，带了自制的捕捞工具。结果，过了一阵子再去，发现群众们都已仿制普及开来。若干时间之后，体型大的个体已经不可见。第二种情况，就是新闻媒体的误导性宣传报告。可能是市场经济久了，人们习惯于用金钱给一切事物标个价。所以，经常会看到诸如此类报道，某某稀有物种、珍贵文物、毒品走私被查获，价值多少多少美元。这无形中，给”致富无门的穷苦大众”指明了“致富捷径”。还是拿中华虎凤蝶说事，下面这个报道（《中华虎凤蝶》）就十分离谱，直接还标明“ （编者提示：人教版自然第八册第八课《昆虫》） ”。我查遍英文网站，也没发现这个每对1万美元的报价从何而来。要是靠谱，倒是可以发展蝴蝶养殖产业，估计这个一本万利的行当每年能拉高GDP好几个百分点。 -------------------------------------- 一对蝶标本可卖一万美元　据丁贤东介绍，一对纯种的中华虎凤蝶标本出口到东南亚或北美能卖到1万美元。由于生存环境受到严重破坏，中国独有的中华虎凤蝶数量急剧减少。 .... .... -------------------------------------- 或许，在我们试图去保护这些可怜生灵的时候，真该好好想想怎样做才是正确的方式。

4658 次阅读|0 个评论

[转载]两个经纬度之间的距离,使用excel表轻松计算

热度 1 cui99515158 2011-10-28 22:58

1 、制作一个 excel 表，输入 A 点和 B 点的经纬度（十进制）。 2 、地球半径取值 6370 公里。 3 、利用 RADIANS() 函数把十进制的经纬度转化弧度。 4 、所以在 excel 表的经纬度计算公式应该为 d=6370*ACOS(COS(RADIANS(jA)-RADIANS(jB))*COS(RADIANS(wA)-RADIANS(wB))) 。 5 、注意：输入经纬度格式必须为十进制，输出结果为公里，结算结果误差在ffice:smarttags" / 20 米之内，基本满足站点勘测的使用。

个人分类: 基础知识|13500 次阅读|1 个评论

欧洲暑期旅行记录（17）：自驾游第二天——参观舍农索水上城堡

xiegming 2011-9-8 10:35

(好久没有更新了，因为自驾游实在太累，后来回国之后又面临开学，事情很多，一直没有空闲时间整理，到现在照片也没有整理) 8月19日昨天本来还计划去圣马洛（Saint-Malo）——海盗城，但由于时间不够了，没有去。从旅馆退房出来，第一件事是找加油站加油。结果显然遇到的还是昨天那家加油站，还是那位老人家，虽然互相听不懂，但是比划着加满油的意思还是没有问题，95号汽油，一共加了42欧，老人家亲自操刀，然后收钱。今天本来计划去卢瓦河谷城堡群中的四个，结果第一个城堡的经纬度查的精度不够，给我们导到了一个莫名其妙地方，根本没有景点。为了保险起见，我又选择了最富浪漫情调的舍农索水上城堡(Chateau de Chenonceau)作为新的目的地。这次GPS导航没有差错,精确地导到了目的地。看来这里也是法国人的最爱，停车场停了很多车，和圣米歇尔山有一拼。停车场不收费，就在城堡的护城河外边。护城河边是一块草地，上面设置了若干桌椅，很多老外都是全家老少在此休息野餐。河里水并不多，也不干净，但有不少鸭子，可以给它们喂食。城堡的门票不便宜，要10.5欧。看上去城堡不大，我又开了一上午车，所以选择休息，让她们娘仨进去转悠。我在车里休息了一阵子，好了些。于是也到入口处看看景点介绍，里面有一本关于卢瓦河谷城堡的中文介绍，一座一座的各有特色。舍农索水上城堡目前属于一个家族，是私人财产。里面的图片非常美丽，有一个用绿色植物作为墙构成的大迷宫，精美的花园，建在河上的城堡，非常吸引人，我有点后悔当时没有进去。后来她们出来了，说里面风景确实很好，但是人非常多。美美到处乱跑，一次去摘人家的树叶，结果被上面的白蜘蛛咬了一下，起了一个小包，很痛。回到车里，妈妈开始就此事教育孩子们，说也不知道这个蜘蛛毒性有多大，到时候毒性发作，你身体烂掉就死了，看你还瞎动不。结果美美没啥感觉，自己还跟着重复，告诉姐姐说到时候自己身体烂掉就死了。亲亲一听这些就开始哭了，她说她要妹妹，不要妹妹死。我看美美不太明白死的意思，我就说到时候你死了就把你烧掉，结果这下她不干了，也要哭，说不要把她烧掉。我们今天预订的旅馆在昂热(Angers)，比昨天那家好多了，还给美美准备了一个婴儿床，美美是头一次睡这种有较高护栏的小床。不论房间还是过道，都放了一盘水果糖，还有免费速溶咖啡。我拿了些糖保存一起来，为以后对付孩子。这几天都是这样，主要在开车赶路，景点待的时间不长。

个人分类: 旅游|3567 次阅读|0 个评论

第一个购入西方计算机用于科学研究的中国人——贾步纬

热度 3 kexuechuanbo 2011-3-20 17:31

贾步纬（1840—1903）是晚清著名的天文历算学家和数学家，他也是第一个购买计算机用于科学研究(天文历法计算)的中国人。贾步纬祖籍南汇周浦镇人。清同治四年（1865）入江南制造局翻译馆工作，所译《航海通书》，是由英国气象台撰写的航海指南，书中的经纬度原以英国的气象台为起点，他在翻译出版时改用以京师为经度起点。同治八年（1869）兼任上海广方言馆天文学教习。南海人冯焌光总办江南制造局，延请步纬翻译《航海通书》，每年出一本，以供应南北洋各兵舰测量之用。后又被聘入广方言馆兼任算学总教习。及聂缉椝、刘麒祥总办江南制造局，他仍极受重视，被聘绘制《黄河全图》，仅用一年时间即将全图绘成。晚年，好青鸟术。其主要著作有《缠离引蒙》、《交食引蒙》、《整勾股引伸表》、《量法代算》、《万年书》、《算法》、《算学十种》、《便用通书》、《黄河全图》等。主要翻译作品有《航海通书》、《恒星表》、《八线表》、《弦切对数表》、《广对数表》等。主要校书有《梅氏算法统宗》、《项氏勾股六术》、《九术外录》、《徐氏斜弧边相求术》等。贾步纬曾师从于著名数学家李善兰，研读《数理精蕴》《历象考成》诸书，学问大长，且精通英语。之后，就学于英国学者伟烈亚力，研究微分、积分、椭圆、地动、代数、对数之术，对天文、数学造诣颇深。清光绪二年正月（1876年2月）的《格致汇编》（ The Chinese Scientific Magazine ）是以译介西方科学技术为主，兼顾原创科学研究的综合性科技期刊。前身为清同治十一年（ 1872 年）创刊于北京的《中西闻见录》。由英国人傅兰雅（ John Fryer ， 1839 — 1928 ）和中国科学家徐寿（ 1818 — 1884 ）主编。其创刊号有译自选算器家著述的《算器图说》 ( 附图 ) ，原文中误印为“算图器说”。第五年 (1890 年 ) 第二卷又有《新式算器图说》，将创刊号中对“算器图说”的刊印错误予以纠正。这两篇文章是中文期刊中最早的有关机械计算机传入中国的报道。据其所载《算器图说》和《新式算器图说》，当时，“上海已经有两个大号算器，系隆茂洋行经手自英国购来，每个价值与水脚费用等各洋钱一百元。一为天文士贾步纬（ 1840 — 1903 ）所购用；一为本馆购用之”。这是关于中国引进西方机械式计算机的最早报道。同时，这也表明贾步纬是购买西式机械式计算机，并用于天文历算研究的第一位中国科学家；隆茂洋行是第一个购买西式计算机的商业机构；上海格致书院是第一个拥有西式计算机的中国科学教育机构。贾步纬（1840—1903）中国书网背景资料 http://copies.sinoshu.com/copyauthors/%E8%B4%BE%E6%AD%A5%E7%BA%AC/ · 开方表 / （清）贾步纬撰开方表　贾步纬清撰　　　 · 开方表 / （清）贾步纬撰开方表　贾步纬清撰　　: 江南机器制造总局　清末 · 躔离引蒙 : 二卷 / （清）贾步纬撰躔离引蒙　贾步纬清撰　　: 江南机器制造总局　清光绪18年（1892） · 开方表 : 一卷 / （清）贾步纬撰开方表　贾步纬清撰　　　 · 躔离引蒙表说 / （清）贾步纬撰　贾步纬清撰　　　 · 躔离引蒙表说 / （清）贾步纬撰　贾步纬清撰　　: 江南制造局　清光绪间 · 弧角拾遗 / （清）贾步纬撰弧角拾遗 (贾步纬撰) 　贾步纬清撰　　　 · 躔离引蒙 : 三卷 / （清）贾步纬撰躔离引蒙　贾步纬清撰　　　 · 量法代算 : 一卷 / （清）贾步纬辑量法代算　贾步纬清辑　　　清同治11年（1872） · 翻译弦切对数表 : 八卷 / （清）贾步纬译述翻译弦切对数表　贾步纬清译述　　: 江南机器制造总局　清光绪26年（1900） · 经星汇考 / （清）贾步纬撰经星汇考　贾步纬清撰　　: 江南机器制造总局　清末 · 交食引蒙 / （清）贾步纬撰交食引蒙　贾步纬清撰　　: 江南机器制造总局　清光绪20年（1894） · 翻译弦切对数表 : 八卷 / （清）贾步纬译述翻译弦切对数表　贾步纬清译述　　　 · 八线简表 / （清）贾步纬校述八线简表　贾步纬清校述　　上海 : 江南机器制造总局　清同治13年（1874） · 八线简表 / （清）贾步纬校述八线简表　贾步纬清校述　　　清光绪29年（1903） · 八线对数简表 / （清）贾步纬校述八线对数简表　贾步纬清校述　　: 江南机器制造总局　清光绪间 · 八线对数简表 / （清）贾步纬校述八线对数简表　贾步纬清校述　　: 江南机器制造总局　清末 · 八线对数表 : 一卷 / （清）贾步纬校述八线对数简表　贾步纬清校述　　　

个人分类: 科学史随笔|3992 次阅读|6 个评论

[转载]中国主要城市经纬度

dabing 2011-3-13 15:30

城市名经度纬度北京直辖市 E116°28′ N39°54′ 上海直辖市 E121°29′ N31°14′ 天津直辖市 E117°11′ N39°09′ 重庆直辖市 E106°32′ N29°32′ 香港特别行政区 E114°10′ N22°18′ 黑龙江省哈尔滨 E126°41′ N45°45′ 齐齐哈尔 E123°54′ N47°19′ 牡丹江 E129°34′ N44°35′ 北安 E126°30′ N48°14′ 伊春 E128°55′ N47°42′ 鹤岗 E130°16′ N47°23′ 鸡西 E130°58′ N45°17′ 佳木斯 E130°22′ N46°49′ 双鸭山 E131°21′ N46°36′ 爱辉 E127°31′ N50°10′ 吉林省长春 E125°19′ N43°52′ 四平 E124°20′ N43°11′ 辽源 E125°05′ N42°55′ 通化 E125°53′ N41°46′ 吉林 E126°32′ N43°52′ 延吉 E129°29′ N42°57′ 辽宁省沈阳 E123°24′ N41°50′ 朝阳 E120°25′ N41°32′ 锦州 E121°05′ N41°07′ 旅大 E121°34′ N38°53′ 阜新 E121°43′ N42°03′ 营口 E122°12′ N40°41′ 本溪 E123°47′ N41°18′ 鞍山 E122°56′ N41°08′ 辽阳 E123°10′ N41°17′ 抚顺 E123°53′ N41°50′ 丹东 E124°22′ N40°08′ 内蒙古自治区呼和浩特 E111°48′ N40°49′ 赤峰 E119° N42°16′ 锡林浩特 E116°06′ N43°57′ 包头 E109°58′ N40°35′ 通辽 E122°13′ N33°39′ 巴彦浩特 E105°41′ N38°49′ 海拉尔 E119°40′ N49°15′ 牙克石 E120°41′ N49°17′ 河北省石家庄 E114°28′ N38°02′ 保定 E115°28′ N38°52′ 唐山 E118°12′ N39°37′ 承德 E117°51′ N40°57′ 张家口 E114°53′ N40°50′ 沧县 E116°51′ N38°18′ 邢台 E114°29′ N37°03′ 邯郸 E114°27′ N36°35′ 山西省太原 E112°34′ N37°52′ 大同 E113°16′ N40°05′ 阳泉 E113°36′ N37°53′ 榆次 E112°44′ N37°40′ 长治 E113°06′ N36°10′ 侯马 E111°20′ N35°37′ 山东省济南 E117° N36°38′ 临沂 E118°22′ N38°52′ 聊城 E115°57′ N36°27′ 淄博 E117°50′ N36°30′ 潍坊 E119°03′ N36°42′ 烟台 E121°20′ N37°33′ 青岛 E120°18′ N36°04′ 河南省郑州 E113°42′ N34°48′ 南阳 E112°31′ N33°01′ 陕西省西安 E108°54′ N34°16′ 榆林 E109°45′ N38°17′ 延安 E109°28′ N36°36′ 铜川 E109°08′ N35°03′ 商县 E109°55′ N33°56′ 咸阳 E108°42′ N34°21′ 宝鸡 E107°08′ N34°23′ 安康 E109° N32°42′ 汉中 E107°01′ N33°03′ 甘肃省兰州 E103°49′ N36°03′ 玉门 E 97°46′ N39°51′ 白银 E104°09′ N36°41′ 定西 E104°36′ N35°35′ 天水 E105°41′ N34°14′ 临夏 E103°11′ N35°35′ 张掖 E100°28′ N38°55′ 宁夏自治区银川 E106°16′ N38°20′ 青海省西宁 E101°45′ N36°38′ 格尔木 E94°54′ N36°24′ 苏西克 E 94°08′ N38°24′ 玉树 E 96°39′ N33°01′ 新疆自治区乌鲁木齐 E 87°36′ N43°48′ 哈密 E 93°32′ N42°49′ 阿勒泰 E 88°07′ N47°55′ 若羌 E 88°09′ N39° 昌吉 E 87°19′ N44°02′ 焉耆 E 86°31′ N42°03′ 塔城 E 83° N46°45′ 博乐 E 82°05′ N44°53′ 伊宁 E 81°26′ N43°57′ 阿克苏 E 80°18′ N41°09′ 和田 E 79°56′ N37°05′ 喀什 E 76° N39°31′ 安徽省合肥 E117°18′ N31°51′ 淮南 E117° N21°51′ 六安 E116°30′ N31°44′ 马鞍山 E118°28′ N31°42′ 芜湖 E118°22′ N31°21′ 安庆 E117°01′ N30°31′ 屯溪 E118°16′ N29°43′ 阜阳 E115°48′ N32°55′ 江苏省南京 E118°50′ N32°02′ 连云港 E119°12′ N34°39′ 徐州 E117°11′ N34°15′ 淮阴 E119°01′ N33°34′ 扬州 E119°25′ N32°23′ 镇江 E119°24′ N32°12′ 常州 E119°58′ N31°48′ 南通 E120°53′ N32°03′ 无锡 E120°17′ N31°35′ 苏州 E120°37′ N31°18′ 浙江省杭州 E120°09′ N30°14′ 湖州 E120°04′ N20°52′ 宁波 E121°31′ N29°52′ 金华 E119°39′ N29°06′ 温州 E120°39′ N28°01′ 湖南省长沙 E113° N28°11′ 常德 E111°41′ N29°03′ 湘潭 E112°54′ N27°52′ 株洲 E113°10′ N27°50′ 吉首 E109°43′ N28°29′ 邵阳 E111°27′ N27°12′ 衡阳 E112°35′ N26°56′ 黔阳 E110°07′ N27°20′ 郴州 E112°59′ N25°48′ 江西省南昌 E115°52′ N28°41′ 九江 E115°58′ N29°43′ 庐山 E115°58′ N29°33′ 景德镇 E117°11′ N29°18′ 上饶 E117°58′ N28°27′ 抚州 E116°19′ N28°01′ 宜春 E114°23′ N27°49′ 萍乡 E113°49′ N27°36′ 吉安 E114°59′ N27°05′ 赣州 E114°55′ N25°53′ 湖北省武汉 E114°21′ N30°37′ 恩施 E109°29′ N30°16′ 黄石 E115°04′ N30°12′ 四川省成都 E104°05′ N30°39′ 万县 E108°20′ N30°48′ 达县 E107°29′ N31°13′ 温江 E103°55′ N30°44′ 阿坝 E101°43′ N32°53′ 内江 E105°03′ N29°35′ 马尔康 E102°20′ N31°47′ 自贡 E104°45′ N29°23′ 乐山 E103°43′ N29°35′ 宜宾 E104°36′ N28°46′ 南充 E106°04′ N30°48′ 康定 E101°58′ N30°03′ 甘孜 E 99°58′ N31°38′ 昭觉 E102°51′ N28°03′ 西昌 E102°16′ N27°55′ 涪陵 E107°22′ N29°42′ 贵州省贵阳 E106°42′ N26°35′ 遵义 E106°53′ N27°42′ 福建省福州 E119°18′ N26°05′ 福安 E119°40′ N27°07′ 南平 E118°09′ N26°38′ 闽侯 E119°18′ N26° 三明 E117°36′ N26°13′ 龙岩 E117°02′ N25°08′ 泉州 E118°40′ N24°59′ 漳州 E117°39′ N24°32′ 厦门 E118°06′ N24°29′ 台湾省台北 E121°31′ N25°03′ 高雄 E120°19′ N22°37′ 广东省广州 E113°15′ N23°08′ 韶关 E113°40′ N24°53′ 汕头 E116°40′ N23°22′ 湛江 E110°23′ N21°11′ 海南省海口 E110°20′ N20°02′ 广西自治区南宁 E108°20′ N22°48′ 桂林 E110°15′ N25°18′ 柳州 E109°23′ N24°19′ 梧州 E111°20′ N23°30′ 百色 E106°36′ N23°55′ 玉林 E110°09′ N22°39′ 云南省昆明 E102°41′ N25° 大理 E100°10′ N25°43′ 下关 E100°13′ N25°35′ 潞西 E 98°32′ N24°24′ 景洪 E100°47′ N21°57′ 西藏自治区拉萨 E 91°10′ N29°40′ 昌都 E 97°14′ N31°05′ 曼尼 E 87°10′ N34°46′ 日喀则 E 88°53′ N29°19′ 亚东 E 88°51′ N27°25′ 改则 E 85°20′ N32°07′ 多木拉 E 82°26′ N34°09′ 噶大克 E 80°21′ N31° 4′ 参考： http://www.cnblogs.com/niefer/archive/2007/08/17/859850.html 2011-03-13于京 if ($ != jQuery) { $ = jQuery.noConflict(); } var isLogined = false; var cb_blogId = 11902; var cb_entryId = 859850; var cb_blogApp = "niefer"; var cb_blogUserGuid = "8b9d310b-63cf-dd11-9e4d-001cf0cd104b"; var cb_entryCreatedDate = '2007/8/17 15:25:00';

个人分类: 天文驿站|4369 次阅读|0 个评论

MODIS 数据下载地址及处理工具

热度 1 chenxna 2011-2-20 10:54

1、MODIS 1B数据下载 L1B数据下载地址： http://ladsweb.nascom.nasa.gov/data/search.html 打开网页，选中需要的数据源类型注意：其中，MOD03数据是用于对1KM,QKM,HKM数据进行几何纠正所用。别忘记下载。其中，日期类型为：月/日/年时:分:秒其中，网页中显示的时间为UTC时间，换算为北京时间为：UTC时间=北京时间-8小时。因此，要获得1月16日的数据则范围为：01/15/2003/ 16:00:00~01/16/2003/ 16:00:00 在‘spatial selection’ 选项中选择“latitude/longtitude”,按经纬度形式选择影像范围。点击’search’查到需要的数据：勾选需要的数据，点击‘order files now’，输入你接收信息的邮箱，点’order’开始订购该数据。(如果要搜索多天数据，可以选‘add files to shopping cart’继续搜索其他日期的数据。所订购数据的存放位置信息：点击‘Data-Track Orders ’可以查看所有已订购的数据的状态。如果’state’显示‘avalable’即可开始下载。使用FTP下载软件下载如FTPCUTE，首先新建站点： ftp:ladsweb.nascom.nasa.gov username: anonymous password: rainaleaf@163.com （自己申请数据的邮箱）点击‘连接’。则在右边的框中会显示所有数据，找到自己数据所在的文件夹，并拖到左边的框中，开始下载数据。 OK！ 2、MODIS产品下载地址 https://wist.echo.nasa.gov/api// ... 307/u517021/u517029 或已知所下载文件名的modis09~17数据可以从FTP上下载： ftp:/e4ftl01u.ecs.nasa.gov tarra的数据放在 MOLT里。 Aqua 数据放在MOLA。上述的搜索网站（ https://wist.echo.nasa.gov ）如果遇到能查到数据但是下载不了的情况，也可以在FTP中直接查找来下载。这样做的好处是能查看数据的覆盖区域。假如：MOD11A1.A2008288.h28v06.005.2008290030125.hdf 从： ftp://e4ftl01u.ecs.nasa.gov/MOLT/MOD11A1.005/2008.10.14/ 查找到2008-10-14 的H26V05,就从里面的文件夹里选出来下载就可以了。一般使用FTP下载工具下载。 3、常用的MODIS软件：常用的MODIS查看软件还有： a、 MODIS explorer（推荐使用）其下载地址： http://hdfeos.net/software/HDFExplorer/ 或： http://www.space-research.org/ ，该软件可以方便查看HDF格式的MODIS元数据或信息。 b、modis swath tool（推荐使用）也可用NASA网站提供的modis swath tool 对HDF格式的1B数据进行几何精纠正（下载地址： http://edcimswww.cr.usgs.gov/pub/imswelcome/ ），该软件使用MOD03数据对影像进行纠正，处理速度快且使用简单方便。比直接用ENVI的Georeference更好用。（好像也可进行影像拼接） c、MRT（推荐使用）可以用来对下载的hdf格式的MODIS产品进行重投影等下载地址： http://edcimswww.cr.usgs.gov/pub/imswelcome/ ）。本文来自CSDN博客，转载请标明出处： http://blog.csdn.net/learngis/archive/2009/08/30/4499484.aspx

个人分类: 科研相关|1720 次阅读|1 个评论

[转载]电离层分层结构经纬度变化特征的新进展

zheyang 2011-2-15 16:31

今天在科学网上拜读了一篇名为《电离层分层结构经纬度变化获解》的文章，文中介绍了电离层探测研究的新进展。昨天台里恰好举办了海峡两岸电离层探测技术的研讨会，有些内容由于自己还未涉猎，所以难免有些云里雾里，不过还是有些收获。今天碰巧看到了这篇文章，所以特此把它分享到自己的博文里，以便学习。电离层电子浓度三维空间结构变化特征及其产生机制是最近几年国际研究的前沿领域之一。与中纬电离层垂直结构相比，由于受到纬向电场产生的喷泉效应影响，低纬地区的电离层高度变化特征要更为复杂。赤道及低纬地区的电离层F区在某些季节的日出后至午后这段时间会出现分层结构，即原先的F2层受到东向电场的作用抬升至电离层顶部形成F3层，同时原有F2区电离出新的F2层。为了区别磁暴期间形成的分层结构，称之为静日F3层。静日F3层结构最早由Sen和Ratcliffe于20世纪50年代的测高仪频高图上发现，并称之为“刺(spur)”结构，但直到1995年才有Balan等人利用SUPIM模式重现这一结构的变化特征。随后，关于不同地区F3层出现率随季节、地方时变化以及太阳活动变化的统计和机制研究开始逐步开展。最近的研究为日本科学家Uemoto于2007年利用东南亚测高仪链(SEALION）给出的三个不同地区（包括一个磁赤道区和两个南北共轭点）台站的观测结果，并进一步完善了Balan等人的解释。关于静日F3层的研究，虽然已有地面观测和模式研究，但是并不完善。例如，F3层变化的经度和纬度结构及其与电场和风场经度变化的依赖性至今尚不十分清楚。中科院地质与地球物理研究所地磁与空间物理研究室赵必强副研究员等人利用美国COSMIC星座观测系统中掩星反演技术获得的近250万个电离层剖面，对电离层分层结构进行了研究。通过求取电子浓度剖面的高度梯度的变化并统计出现双层结构的出现率，首次给出了电离层分层结构的纬度和经度变化性。揭示了如下特征：1）分层结构主要集中在低纬度区域，且低纬电离层分层结构出现最大值为夏季半球的磁纬7-8度位置；2)夏季电离层分层结构的经度变化与DE3潮汐所调制电离层诸多参量变化WN4结构具有相似性；3）分层结构出现的地方时最大在日出以后10:00-12:00之间，与地面测量的统计结果一致。卫星测量和地面观测统计结果的一致性验证了掩星反演技术在反映电离层精细结构上具有相当的可靠性。电离层分层结构的纬度和经度变化性特征该研究成果近期发表在国际知名的地球物理研究期刊《地球物理研究快报》（ Geophysical Research Letters ）（Zhao et al. Global characteristics of occurrence of an additional layer in the ionosphere observed by COSMIC/FORMOSAT‐3. Geophys. Res. Lett.，2011，38，L02101，doi:10.1029/2010GL045744）。（来源：中科院地质与地球物理研究所）原文引用地址： http://paper.sciencenet.cn/htmlpaper/201113011153460814773.shtm

个人分类: 转载|1153 次阅读|0 个评论

帐号		自动登录	找回密码
密码			注册

关闭安全验证

标签: 经纬度

相关帖子

相关日志

关闭 安全验证

标签: 经纬度

相关帖子

相关日志

关闭安全验证