小编daz*_*sed的帖子

尽管几年前学习了python(2.7),但我才开始再次使用它.我正在使用python 3.4.4并尝试使用GDAL包将栅格作为数组读取.我按照这里概述的步骤操作: https ://sandbox.idre.ucla.edu/sandbox/tutorials/installing-gdal-for-windows, 但Python 3.4除外,并使用了GDAL二进制文件:http://www.gisinternals. COM/release.php

在测试基本功能时,我试图读取如图所示的tiff文件

import gdal as gdal
import numpy as np
import ogr
import osr
import os
import sys
e=('error has occurred')

# this allows GDAL to throw Python Exceptions
gdal.UseExceptions()

# open dataset
test = ('LE70130312004049EDC01_sr_adjacent_cloud_qa.tif')
print("file exists")


# getting metadata
gtif = gdal.Open(test)
print (gtif.GetMetadata())
print("metadata printed")

try:
    src_ds = gdal.Open(test)
    print("gdal.open success")
except (RuntimeError):
    print ('Unable to open INPUT.tif')
    print(e)

try:
    srcband = src_ds.GetRasterBand(1)
    print("get raster band 1 success")
except (RuntimeError): …