上个礼拜,和帮隔壁学院某研究蝙蝠的博士做了个图——根据全球哺乳动物的活动范围(多边形),在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的可能觉得很简单,但对别的邻域是一种很强有力的工具,甚至是新颖的。
最最后,胡适先生曾说,多研究些问题,少谈些“主义”。在现在这种“云”里“雾”里的浮躁年代,多解决些问题,多拓展些行业,少吹点牛皮,少扯点概念,或许更好。