cjaca / hyperspectral-image-analysis Goto Github PK

    raster = rasterio.open(img_path, 'r')
    arr = raster.read()
    ref_count = 1
    ref_cord = [[round(arr.shape[1]), round(arr.shape[2])]]
    ……
    # ref_vect is reference pixel
    ref_vect = [[0 for x in range(bands)] for y in range(ref_count)]
    for i in range(ref_count):
        for j in range(bands):
            ref_vect[i][j] = arr[j, int(ref_cord[i][0])-1, int(ref_cord[i][1])-1]