Как отобразить цифровую модель рельефа (DEM) (.raw) с использованием Python?

Question

Я хочу отобразить файл DEM (.raw) с помощью Python, но в результате может быть что-то не так.

Ниже мой код:

img1 = open('DEM.raw', 'rb')
rows = 4096
cols = 4096

f1 = np.fromfile(img1, dtype = np.uint8, count = rows * cols)
image1 = f1.reshape((rows, cols)) #notice row, column format
img1.close()

image1 = cv2.resize(image1, (image1.shape[1]//4, image1.shape[0]//4))
cv2.imshow('', image1)
cv2.waitKey(0)
cv2.destroyAllWindows()

И я получил этот результат: показать результат

Оригинальный файл DEM находится здесь: DEM.raw

Answer 1

В вашем коде нет ничего плохого, это то, что есть в вашем файле. Вы можете преобразовать его в JPEG или PNG с помощью ImageMagick в командной строке следующим образом:

magick -size 4096x4096 -depth 8 GRAY:DEM.raw result.jpg

И вы получите почти то же самое:

Проблема в другом.

Взяв подсказку Фреда (@fmw42) и поэкспериментируя, я имею в виду "экспериментируя осторожно и научно", я могу получить более вероятный результат, если обработать ваше изображение как 4096x2048 пикселей с 16 битами на дюйм и MSB с первым порядком байтов:

magick -size 4096x2048 -depth 16 -endian MSB gray:DEM.raw  -normalize result.jpg