NumPy없이 OGR 포인트에서 GDAL 래스터의 픽셀 값을 얻습니까?

나는 조경을 가로 지르는 풍부한 야생 수분 조절제의 계산 모델을 연구하고 있습니다. 모델 자체가 완성되었으므로 이제 사후 처리 단계로 어려움을 겪고 있습니다.

GDAL pollinator 공급 래스터가 다음과 같이 보입니다 (색이 밝을수록 pollinator가 픽셀을 더 많이 방문합니다).

그리고 풍경의 샘플 위치를 나타내는 점의 OGR shapefile이 있습니다.

이 포인트 아래의 픽셀에 대한 분석을 시도하고 있지만 그렇게하려면 포인트 아래의 픽셀 값을 추출 할 수 있어야합니다.

파이썬을 통해 OGR 및 GDAL 만 사용하여 포인트 아래에서 픽셀 값을 추출 할 수 있습니까? ReadAsArray()출력 래스터가 매우 크고 (메모리에 너무 커서) 래스터 전체를 메모리로 읽는 것을 피하고 싶습니다.

이 게시물 은 비슷하지만 명령 줄 호출이 필요하다는 것을 알았 습니다 .

gdal.Dataset 또는 gdal.Band ReadRaster 메소드를 사용할 수 있습니다 . 참고 항목 GDAL 및 OGR API 자습서 아래의 예를. ReadRaster는 numpy를 사용 / 요구하지 않으며 반환 값은 원시 이진 데이터이며 표준 python struct 모듈을 사용하여 압축을 풀어야 합니다.

예를 들면 :

from osgeo import gdal,ogr
import struct

src_filename = '/tmp/test.tif'
shp_filename = '/tmp/test.shp'

src_ds=gdal.Open(src_filename) 
gt=src_ds.GetGeoTransform()
rb=src_ds.GetRasterBand(1)

ds=ogr.Open(shp_filename)
lyr=ds.GetLayer()
for feat in lyr:
    geom = feat.GetGeometryRef()
    mx,my=geom.GetX(), geom.GetY()  #coord in map units

    #Convert from map to pixel coordinates.
    #Only works for geotransforms with no rotation.
    px = int((mx - gt[0]) / gt[1]) #x pixel
    py = int((my - gt[3]) / gt[5]) #y pixel

    structval=rb.ReadRaster(px,py,1,1,buf_type=gdal.GDT_UInt16) #Assumes 16 bit int aka 'short'
    intval = struct.unpack('h' , structval) #use the 'short' format code (2 bytes) not int (4 bytes)

    print intval[0] #intval is a tuple, length=1 as we only asked for 1 pixel value

또는 사용하지 않기로 한 이유는을 사용 numpy하여 전체 배열을 읽지 않기 때문 ReadAsArray()이므로 아래는 numpy전체 래스터 를 사용 하고 읽지 않는 예입니다 .

from osgeo import gdal,ogr
import struct

src_filename = '/tmp/test.tif'
shp_filename = '/tmp/test.shp'

src_ds=gdal.Open(src_filename) 
gt=src_ds.GetGeoTransform()
rb=src_ds.GetRasterBand(1)

ds=ogr.Open(shp_filename)
lyr=ds.GetLayer()
for feat in lyr:
    geom = feat.GetGeometryRef()
    mx,my=geom.GetX(), geom.GetY()  #coord in map units

    #Convert from map to pixel coordinates.
    #Only works for geotransforms with no rotation.
    px = int((mx - gt[0]) / gt[1]) #x pixel
    py = int((my - gt[3]) / gt[5]) #y pixel

    intval=rb.ReadAsArray(px,py,1,1)
    print intval[0] #intval is a numpy array, length=1 as we only asked for 1 pixel value