Python gdal 读取栅格数据
2015-04-03 11:03
706 查看
程序仅供参考。
import sys import gdal from gdalconst import * from osr import SpatialReference fn = r'E:\work\Ex1016\xian_mask1.tif' print fn ds = gdal.Open(fn,GA_ReadOnly) if ds is None: print 'cannot open ',fn sys.exit(1) print 'size is ',ds.RasterXSize,'x',ds.RasterYSize,'x',ds.RasterCount print 'projection is:',ds.GetProjection() geotransform = ds.GetGeoTransform() if not geotransform is None: print 'origin(x,y) is:',geotransform[0],',',geotransform[3] print 'Pixel size is:',geotransform[1],',',geotransform[5] band = ds.GetRasterBand(1) print 'Band Type=',gdal.GetDataTypeName(band.DataType) data=band.ReadAsArray() print 'data is ',data min = band.GetMinimum() max = band.GetMaximum() print 'Nodatavalue is:',band.GetNoDataValue() if min is None or max is None: (min,max)=band.ComputeRasterMinMax(1) print 'Min=%.3f,Max=%.3f' % (min,max) if band.GetOverviewCount() > 0: print 'Band has', band.GetOverviewCount(), 'overviews.' if not band.GetRasterColorTable() is None: print 'Band has a color table with ',band.GetRasterColorTable().GetCount(),\ 'entries'
相关文章推荐
- Python3+gdal 读取tiff格式数据
- GDAL获取栅格数据各个像素对应的经纬度(Python版)
- [GDAL]2.读取栅格和矢量数据
- Gdal浅析(一)读取栅格数据到Dataset
- GDAL读取S-57海图数据中文属性值乱码问题解决(续)
- GDAL\OGR读取数据示例 C#版本
- 关于GDAL读取MODIS卫星数据
- 用 Python 的输入输出功能读取和写入数据
- GDAL读取S-57海图数据中文属性值乱码问题解决
- GDAL支持的栅格数据格式
- 关于GDAL读取MODIS卫星数据
- 使用gdal读取图像数据,然后用构建gdi+位图显示
- 探索 Python,第 8 部分: 用 Python 的输入输出功能读取和写入数据
- 利用GDAL读取图片,将数据以BIP的方式存在内存中。
- Python 通过thrift接口连接Hbase读取存储数据
- GDAL\OGR读取数据示例 C#版本
- 用 Python 的输入输出功能读取和写入数据
- Delphi中使用python脚本读取Excel数据
- 用 Python 的输入输出功能读取和写入数据
- 使用python读取csv文件,并将数据更新至mysql