Archive for the ‘Map’ Category

全球哺乳动物分布图

上个礼拜,和帮隔壁学院某研究蝙蝠的博士做了个图——根据全球哺乳动物的活动范围(多边形),在1×1经纬网上,获取每个网格上活动(栖息?)的哺乳动物种类数目。

原始数据是这样的:

我处理(Dissolve)了下,有5398种哺乳动物,和维基上说的差不多。我关注的是活动范围,其中,活动范围最小的是Melomys rubicola(一种老鼠?)范围小到几乎可以忽略不计啊,在(143,-9)附近;最大的是Orcinus orca(虎鲸),基本上有水的地方都是它的地盘,这就叫四海为家吧。

先生成经纬网(ArcPy):

然后按Intersect进行SpatialJoin(本来想还是用ArcPy直接遍历每个几何进行overlaps、contains和within比较,无奈数据太大,要先进行几何的envelope比较,这个不涉及空间操作,但是写起来多稍微麻烦了点呢,而且为什么SearchCursor获取的几何不能保存在list中?不解)。

之后统计网格数目,再将统计表Join到网格数据中,完事。

SpatialJoin那步效率太差,反正我晚上回去开始跑,第二天早上来看是跑完的。后来看了下,最后有190多万个多边形。。。想想也是,光光经纬网格就有360×180了,平均下来每个网格上才3种哺乳动物吧。

最后的图:

没受过专业的制图训练,见谅了,虽然我们专业前面是有Cartography这种说法的,但现在还有人(老师)知道有制图这回事吗?

Scene里面拉伸下:

最后,上上个礼拜在3sNews上分享了下Geospatial Revolution的视频,还不算是无人问津。我觉得这种制图,就是EP4中介绍的地理空间技术在科研方面的应用之一。学GIS的可能觉得很简单,但对别的邻域是一种很强有力的工具,甚至是新颖的。

最最后,胡适先生曾说,多研究些问题,少谈些“主义”。在现在这种“云”里“雾”里的浮躁年代,多解决些问题,多拓展些行业,少吹点牛皮,少扯点概念,或许更好。

南京,PX离我们有多远?

首先,要明白PX是什么?

http://zh.wikipedia.org/wiki/%E5%AF%B9%E4%BA%8C%E7%94%B2%E8%8B%AF

其次,了解一下什么地方的什么人对PX做过什么?

http://zh.wikipedia.org/wiki/%E5%8E%A6%E9%97%A8PX%E9%A1%B9%E7%9B%AE

然后,南京还有PX啊?

http://news.sina.com.cn/s/2006-05-18/07578959048s.shtml

再然后,南京市政府怎么看PX?

http://www.njbbs.gov.cn/simple/?t25324.html

继续然后,同学们怎么看?

http://bbs.nju.edu.cn/bbstcon?board=V_Suggestions&file=M.1307188576.A&start=0

等等,我列的肯定不全,希望有人能找更多的资料。

例如,你可以Google一下:金陵化工、PX、栖霞、对二甲苯,等等关键字。

不过如果你有干爹,那就不用这么折腾了。

疑问:PX危险吗?

这我不知道,如果你知道,那你找出依据来告诉我吧。

 

但是,PX离我们有多远?

这个我知道:


View Larger Map

推荐看大图

你知道PX离我们有多近了吗?

Earthquakes with 1,000 or More Deaths since 1900

USGS上的有各种灾难的数据,包括地震的,例如这个:Earthquakes with 1,000 or More Deaths since 1900,1900年以来,死亡1000人以上的地震数据。当然,还有许多其他的,比如:所有8级以上的地震、死亡5000人以上的地震、按地区分类的以及各种专题图等等。在这么多数据中,我选择了Earthquakes with 1,000 or More Deaths since 1900这个,结合ArcGIS 10中的时态数据功能来做地理数据的时态显示。

我的思路是按数据处理、地理数据符号化、时态显示三部分来做:

  • 数据处理:

将网页的数据复制到文本中,进行一些必要的手动处理,使数据更为标准,将每两行作为一条记录进行读取:

以行为单位,需要提取的数据是:日期、位置、经纬度、死亡人数、震级。

使用ArcPy生成地理数据,创建要素类,并添加日期、位置、经纬度、死亡人数、震级等字段:

  • 地理数据符号化:

以震级和死亡人数为字段,进行渲染。颜色由绿至红表示死亡人数从小到大;圆的半径越大,表示震级越高。显然,两者不成正比的:

  • 时态显示:

为图层启用时间功能,打开时间滑块窗口,就能进行时态数据显示了:

 

视频: