Saturday, October 1, 2016

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

No comments:

Post a Comment