牛虻的天空

清明节二三事

gadfly621 — Thu, 07 Apr 2011 14:49:51 +0800

清明好像“三节”合一，不宜出行，就在附近转了转，二三事记录一下。

1）三个小朋友，5 、6岁的样子在大人的引领下在公园游玩，一个瘦点的小男孩对另外两个说，你知道世界上有几只毒蛇吗？很诡异的问题，我想了想真不知道，“五种，不，十种”他想了想，其他人估计都很疑惑这小孩不简单啊，他说：有泰山蛇、黄山蛇、香山蛇…..

2）一个10多岁得胖乎乎的孩子买了个香肠，自己没吃，对妈妈说：“妈妈你吃一口”。心想这孩子真孝顺。“一小口啊，别吃多了”他补充道。呵呵……

3)周末，在中体锻炼，来的基本上都是附近科学院小区里的人，更衣间两个大人对话，“你小孩去美国还是加拿大啊”“TA想去美国，可是加拿大好移民”另一个3、4十岁的人说，“出去让TA学会做人，国内有自然科学嘛？！除了以前航天科工、原子弹有点真事，其他全是骗子……”，想想自己也在他们口中行骗的道路上渐行渐远……

MODIS数据处理的两个软件

gadfly621 — Thu, 31 Mar 2011 01:22:47 +0800

1 官方软件 MRT (MODIS Reprojection Tool)

基本的介绍和下载地址：https://lpdaac.usgs.gov/lpdaac/tools/modis_reprojection_tool。
功能：可以实现基本投影变换及重采样处理。
界面：

安装：（1）界面运行需要安装JRE (Java Runtime Environment)，在win下要重启,Linux 下不用
批处理：win/linux下均需要shell环境，相对较复杂

2 ITT Modis数据处理软件包The MODIS Conversion Toolkit (MCTK)

其基本的介绍和下载地址：http://www.ittvis.com/language/en-US/UserCommunity/CodeLibrary/id/1193.aspx。
功能：可以实现基本投影变换及重采样处理。
界面：

安装：（1）需要ENVI环境，实际是ENVI的二次开发
批处理：IDL环境即可，比较方便

BJ54投影坐标转换为经纬度地理坐标

gadfly621 — Tue, 15 Mar 2011 23:26:18 +0800

同学发过来两列数据，告诉是BJ54的坐标系，让转化为经纬度坐标，如下表所示，很无语也很无奈，搞搞吧。

表1 BJ54投影坐标

控制拐点	X坐标	Y坐标
1	4924000	20428200
2	4924000	20429400
3	4922800	20428200
4	4922800	20429400
1	4878997	433164
2	4878997	435075
3	4875884	435075
4	4875884	433164
5	4878997	433164

BJ54在测绘中应用很广泛，尤其是我国早期的测绘成果大多以此为基础，但是在遥感影像中通常使用WGS-84坐标系统，所以平时用的并不多，没办法硬着头皮试一下吧。首先发现一个奇怪的问题，这些坐标尽管是在一个地区的，但是给的Y坐标却位数并不一样，第一组前面多一个“20”，很纳闷，后来想起学过BJ54坐标是有时是有分带号的，查了一下果然如此。

表2 BJ54坐标系统参数

Geographic Coordinate System:	Projection coordinate system
Name: GCS_Beijing_1954 Alias: Abbreviation: Remarks: Angular Unit: Degree (0.017453292519943299) Prime Meridian: Greenwich (0.000000000000000000) Datum: D_Beijing_1954 Spheroid: Krasovsky_1940 Semimajor Axis: 6378245.000000000000000000 Semiminor Axis: 6356863.018773047300000000 Inverse Flattening: 298.300000000000010000	Projection: Gauss_Kruger Parameters: False_Easting: 500000.000000 False_Northing: 0.000000 Central_Meridian: 117.000000 Scale_Factor: 1.000000 Latitude_Of_Origin: 0.000000 Linear Unit: Meter (1.000000)

转化的过程中，使用两种软件进行转化，包括proj软件法跟ARCMAP软件转化，可是转化过程中又发现一个问题，貌似同学X Y坐标输入反了，后来仔细查了一下原来BJ54系统定义Y坐标是横轴，X坐标为纵轴，这个问题狠狠郁闷了一下。

1. PROJ方法：

PROJ是一款很好的坐标转化工具，非常适合学习与使用，应该也是开源的。通过使用命令就可以实现转化(斜体为输入，黑体为输出结果)，其参数可以参考表2：

C:\>proj.exe -I +proj=tmerc +x_0=500000 +y_0=0 +lon_0=117 +lat_0=0 +units=m +k=1 +ellps=krass

428200 4924000

116d5’52.872″E 44d26’50.038″N

429400 4924000

116d6’47.135″E 44d26’50.463″N

2. ARCMAP方法：

将数据参考前述文章，转化为点矢量文件，将输出的矢量文件在ARCMAP打开后将DATA FRAME PROPERTIES ->COORDINATE SYSTEM –>PREDEFINED->Geographic Coordinate System:设置为BJ54大地经纬度坐标的形式进行显示。然后在矢量文件属性表中，增加LON，LAT两列，分别选择右键单击calculate geometry，就可以将经纬度坐标添加到两列中，然后EXPORT输出即可。

参考文献：

move blog again!

gadfly621 — Mon, 07 Mar 2011 21:18:10 +0800

move my blog from gaoshuai.info to wordpress.com. Welcome everyone!

越南印象

gadfly621 — Mon, 07 Mar 2011 20:48:33 +0800

北京时间2010年10月31日从北京出发到底河内的内排机场，跟首都机场的设施和服务相比，河内的“首都”机场就显得寒酸了很多，而且效率很低，出境排队等了半个小时，取行李大于半个小时，好在这个季节晚上9点多的天气还不错，不冷不热，不然真有些受不了。由于飞机晚点，晚饭时间被大大推迟，好在当地旅行团安排的比较周到，晚上吃的还算是不错，尤其是有名的越南米粉，很受用。

从机场到市区的酒店道路并不是很好，跟北京相比差很多，尤其是越南到现在为止还没有高速公路。而且在路上摩托车与其他机动车不分，而且摩托超多，行驶十分缓慢。不过好在道路两旁的建筑很特别，跟中国的很不同，临街的门面很小，但是很瘦长，而且都3 4 层高的小楼，有点像几个集装箱磊在一起的那种感觉。据说是因为在河内地皮很贵，仅次于纽约跟东京。

上午参加这次会议的主办地越南的国家会议中心(NCC)，修的很气派，前几天东盟+3的会议就是在这里召开的，除了体会了越南的效率确实有些低之外，还有就是coffee time的美食还算是不错。下午去了越南有名的购物地，有点类似于五道口批发市场的地方，东西必须仔细挑，有次有好，优劣不一，而且越南盾很不值钱，搞的比较糊涂。

会议组织的welcome party算是比较有意思的一个项目，但是跟去年在北京相比十分没的比，无论是从美食还是节目都是如此，不过中国跟中国台北的节目还是很值得一看，呵呵，很有意思。

走马观花了越南的著名景点，胡志明主席的故居及现居、文庙、还剑湖，当然也少不了购物以及其他，在一个类似中国县城的地方，参观了充满中国改革开放以前的一些建筑，历史，就仿佛还在国内一般，要不是沿街越南风味的门头，以及越南语还真以为是在中国南方的某个城市。

行程最后一天去了下龙湾，很美丽的一个地方

十一回家

gadfly621 — Mon, 07 Mar 2011 20:44:22 +0800

本来报名参加关于农业问题的博客大赛的，结果因为开题的事情+其他事情耽误，在此瞎想一篇，慰藉一下自己心灵吧。

从来没有十一回家的习惯，因为担心人多，可是今年应景回去了一次，人多自不必多言，被人推到车上，足足站了8个半小时到家。回去无怪乎看望一下父母他们，尽管电话里说会不会都行，有事就尽量别回，可是我知道他们是期望我回去的，这个世上恐怕再也找不到比父母更能宽容自己的人了。没有带钱，也没有带任何东西，孤零零一个人站着回来，在公交车上靠着太阳照的暖暖的玻璃窗，美美的睡了一觉，回家的感觉真好。

家里的农活很多，刚刚过八月十五，正是收获的季节，收玉米种小麦，都要在这一两个星期结束，好在我回家之前，玉米已经收回来了，家里堆得满满的，好像是6亩地的吧，老爸说今年收成不好，都让

Hello world!

gadfly621 — Fri, 08 Oct 2010 08:56:17 +0800

Welcome to WordPress.com. This is your first post. Edit or delete it and start blogging!

摘抄一条评论

gadfly621 — Mon, 27 Sep 2010 22:02:59 +0800

一直以来不怎么喜欢科学网，至于原因在twitter上说过，好话不说二遍。不过今天看了一篇头版报道《NASA发布全球空气质量地图中国空气全球最差》稍微仔细看了一下nasa报道跟刊物文章作者原文，感觉还是有相当出入，哗众取宠之嫌。因为不是研究大气，所以也不好评论。看了一下跟帖，感觉这条还比较好（应该是搞遥感的），摘抄如下：

因为搞一点遥感（但不是大气）所以还是比较感兴趣，所以就将NASA的报道跟related links 里的作者的全文看了一下。
首先，觉得科学网的编辑明显在学《自然》，将PM2.5的分布图当成了大气污染的程度分布图，在NASA的报道里没有拔高的这样的高度，作者为了文章发表在引言里夸张几句也就罢了。
再说原文：
第一作者是大气的，其他作者明显是偏重环境健康领域多一些，对于遥感可能只是拿过来就用，事实也是如此。文章得到的结果主要基于MODIS的产品和他们自己的模型。前者主要是modis+misr,用的是MODIS AOD产品和MODIS BRDF产品(这也是NASA为什么高兴的原因了)，还是比较有争议的，其产品的精度到底有多高，反正我见到的质疑的文章就很多，不仅在中国，美国自己人也质疑，产品尺度那么大，验证点那么少，不足为怪，不信可以问lix。作者的计算机模型没接触过，但是其验证也主要是在北美，北美人当然在家门口研究，中国点很少，主要是beijing,shanghai,guangzhou，这几个点就给出了整个中国的统计偏差，很不靠谱。
一个不靠谱的产品乘以一个不靠谱的模型就得到了那副图（这是作者原文里唯一的数学公式，不信自己看）。

支持电子农务

gadfly621 — Tue, 07 Sep 2010 09:19:08 +0800

unix-center的蒋老大发出号召，要为农村做点事情，想来自己一直在农村长大，一直到大城市与农村的家之间奔波，还是很喜欢农村的环境的，捐电脑暂时不行，不过可以做一些其他的事情，例如这次博客大赛，时间短，奖品丰厚，呵呵。

欢迎光临主页：http://www.dznw.net/home/?1021

改变科研文化？

gadfly621 — Sun, 05 Sep 2010 11:25:56 +0800

刚刚拜读了施一公、饶毅最近在《SCIENCE》上的檄文 (中文)《China’s Research Culture》。列举了中国科研的一些现状，包括巨型科研经费的管理分配，科学家研究员科研态度各方面以及科研圈中人的态度来谈中国科研的现状，陈述了事实，也讲了些道理，但是老实说对于改变这种现状并不抱有大的期望，原因很简单因为我们的科研文化其实反映的就是我们的文化，往大里说是中华民族文化的一部分，不可分割的整体，整体左倾，而你要想部分右倾除非你把它割下来。老实说对于新中国又没有正确的科研文化一直以来很疑问，因为我们知道的中国的牛的有成就的科学家都是民国培养的，共和国更喜欢工程师而不是科学家，所以既希望改变这种几十年有点不太正确的传统，又希望改变这种错综复杂的交织的现状是相当不现实的，但是说出来总是有意义的，尤其是从有成就的科学家口中说出来尤其如此，起码让大家感到一种希望。