Saturday, October 1, 2016

Matlab+GDAL: 空间数据重采样

Introduction

GDAL包含10种以上的空间重采样方法。我们就利用这10种以上的方法进行空间重采样过程。

Procedure

变量med是重采样方法,方便自由切换。重采样参数来自标准文件。如Fig. 1所示,从上至下依次是输出栅格、输入栅格、标准文件。
Fig. 1

References

[1] gdalwarp.

Matlab+GDAL: 返回栅格数据属性

Introduction

利用GDAL取得栅格数据的属性,这些属性返回至Matlab以作后用。
这部分代码是以Matlab调用Python函数形式完成,当然,栅格属性还有很多种,从参考文献[1]可以仔细了解,需要哪种属性就在原有的.py代码上稍加修改就可以。

Procedure

以返回栅格数据的空间投影信息为例,返回信息如下:
1
PROJCS["unnamed",GEOGCS["WGS 84",DATUM["WGS_1984",SPHEROID["WGS 84",6378137,298.257223563,AUTHORITY["EPSG","7030"]],AUTHORITY["EPSG","6326"]],PRIMEM["Greenwich",0],UNIT["degree",0.0174532925199433],AUTHORITY["EPSG","4326"]],PROJECTION["Albers_Conic_Equal_Area"],PARAMETER["standard_parallel_1",25],PARAMETER["standard_parallel_2",47],PARAMETER["latitude_of_center",0],PARAMETER["longitude_of_center",105],PARAMETER["false_easting",0],PARAMETER["false_northing",0],UNIT["metre",1,AUTHORITY["EPSG","9001"]]]

References

Monday, September 26, 2016

CDR: Matlab制图

Introduction

CDR在制作标准数学曲线方面并不容易,它强在编辑功能强大,另一方面,Matlab制作标准的数学曲线却是非常容易、快速,但M本身编辑功能比较繁琐,那么M制作的曲线可以直接导入CDR编辑吗?强强联合,答案是可以的。

Procedure

首先,在M中制作好所需要的曲线,这里以正弦曲线为例。制作的图形用export_fig函数以eps文件输出,参考下面代码,输出图片名称是:sinewave.eps。
1
export_fig sinewave -eps
注意:图片的制作环节可能设置背景透明的效果更佳。
第二步,在CDR中导入sinewave.eps,注意导入过程选择如Fig. 1所示。
Fig. 1
导入之后,图形出现在CDR的编辑区,如Fig. 2所示该图形是透明的,之后就可以按照自己的要求尽情编辑了。
Fig. 2

Saturday, September 24, 2016

Tool: 人民网地方领导留言板

Introduction

人民网地方领导留言板已经建立超过10年,这里或许也是一条维护权益的途径。
在这里比较一下北京和天津的留言回复情况,截图在2016年9月25日。
先来看看首都的情况(见Fig. 1),书记、市长的回复都是0条,最高回复数量出现在西城区,是个2。我想普通人看到这样的场面,一般就不再抱有希望了。
Fig. 1
再来看看有地方文件加持的天津版(见Fig. 2),滨海新区留言回复比例96.06%,看起来还可以。
Fig. 2

References

[2] 关于人民网等网上留言办理工作的规定.