Масштабирование от 16 до 8 бит в диапазоне с использованием GDAL

Я хочу масштабировать 16-битное изображение без знака в 8-битное изображение с помощью GDAL (OSGeo4W). Однако я хочу ограничить значения пикселей, чтобы они оставались в пределах 8-битных изображений. Я хочу убедиться, что визуальная пропорция в 16-битном изображении представлена в 8-битном изображении (0-255). Использование gdal_translate для преобразования из 16-битного в 8-битное сокращает значения пикселей и не создает RGB-изображение. Я не уверен, какую функцию и настройку вызывать из gdal_translate, или есть ли лучший вариант.

GDALINFO

C:\>gdalinfo C:\Projects\Stormwater\ForPCI\images\1537TORO_6cm_04761_RGBI.tif
Driver: GTiff/GeoTIFF
Files: C:\Projects\Stormwater\ForPCI\images\1537TORO_6cm_04761_RGBI.tif
Size is 17310, 11310
Coordinate System is `'
Metadata:
  TIFFTAG_DATETIME=2015:05:03 12:27:13
  TIFFTAG_IMAGEDESCRIPTION=UltraCam-Lvl03
--------------
CAM_ID: UC-SXp-1-50215465 [3]
IMG_NO: 4761
RECORD_GUID: 8fc3d7a9-d2e5-40e5-babf-939f1f803dcf
IMG_GUID: C850D0A4-39CA-4BC8-ACD7-81A363D78C6E
FILE_GUID: 86C4F8E4-A2FD-43FB-B3B1-4C4ED8DF0E42
LICENSE_ID: 820342819
SOFTWARE: UltraCam Aerial Radiometry Core 12.6.1408.2501
PIXEL_SIZE_WIDTH: 6 [micron]
PIXEL_SIZE_HEIGHT: 6 [micron]
APERTURE: F_8
EXPOSURE_TIME: 0.002000000000 [s]
HIGH_ISO_MODE_CAPTURE: off
HIGH_ISO_MODE_PROCESSING: off

----- Inner Orientation -----
PRINCIPAL_DISTANCE: 100.500000000000 [mm]
PRINCIPAL_POINT_X: -0.000000000000 [mm]
PRINCIPAL_POINT_Y: 0.180000000000 [mm]
SENSOR_AREA_WIDTH: 103.860000000000 [mm]
SENSOR_AREA_HEIGHT: 67.860000000000 [mm]
-----------------------------

----- Exposure Annotation Data -----
MID_EXPOSURE_CORRECTION: 0.000762000000 [s]
FMS_SENSOR_CODE: UCXp
FMS_CAMERA_PORT: 1
FMS_EXPOSURE_NUMBER: 4761
FMS_PROJECT: 1537TORO
FMS_AREA: 1537TORO_6cm_6030
FMS_LINE_NUMBER: 47
FMS_SEGMENT_NUMBER: 1
FMS_WAY_POINT_NUMBER: 0
GPS_DATE: 150412
GPS_TIME: 162517
GPS_LATITUDE: N43.745319 [degree]
GPS_LONGITUDE: W079.358953 [degree]
GPS_ALTITUDE: 1159 [m]
GPS_POSITION_SOLUTION: GPS
GPS_TRACK_OVER_GROUND: 73 [degree]
GPS_ABOVE_GROUND_LEVEL: 1008.7 [m]
GPS_GROUND_SPEED: 80.4 [mps]
------------------------------------

----- Level-3 Parameters -----
IMG_TYPE: High resolution multi channel RGBI
ROTATION: 0 [degree]
------------------------------------
TIFFTAG_RESOLUTIONUNIT=1 (unitless)
TIFFTAG_SOFTWARE=UltraCam Aerial Radiometry Core 12.6.1408.2501
Image Structure Metadata:
INTERLEAVE=PIXEL
Corner Coordinates:
Upper Left  (    0.0,    0.0)
Lower Left  (    0.0,11310.0)
Upper Right (17310.0,    0.0)
Lower Right (17310.0,11310.0)
Center      ( 8655.0, 5655.0)
Band 1 Block=17310x1 Type=UInt16, ColorInterp=Red
Band 2 Block=17310x1 Type=UInt16, ColorInterp=Green
Band 3 Block=17310x1 Type=UInt16, ColorInterp=Blue
Band 4 Block=17310x1 Type=UInt16, ColorInterp=Undefined

gdal imagery gdal-translate

— Райан Гарнетт
источник

Не могли бы вы предоставить некоторые подробности о типе изображений, которые вы используете?

— Аарон

Например, прикрепив отчет gdalinfo к исходному файлу. Напишите также точную команду gdal_translate, которую вы использовали.

— user30184 14.11.16

gdalinfo добавлен в вопрос как ОБНОВЛЕНИЕ. Образы - стереоизображения 3-го уровня. Была использована команда gdal_translate: gdal_translate -ot Byte input.tif output.tif

— Райан Гарнетт

Возможно, реальные данные 16-битных изображений используют только узкую часть полного диапазона. Проверка гистограммы, например, с помощью QGIS может выявить это. Вы также можете просто попробовать с параметром -scale gdal.org/gdal_translate.html . Если вы хотите получить RGB-вывод из 4-полосного оригинала, выберите три из них с параметрами -b, например-b 1 -b 2 -b 3

— user30184

Если вы не хотите, чтобы значения выше 255 были сокращены, вам нужно уменьшить их. Для этого gdal_translate предоставляет опцию -scale:

Из руководства :

-scale [src_min src_max [dst_min dst_max]]: изменить масштаб значений входных пикселей от диапазона src_min до src_max до диапазона от dst_min до dst_max. Если опущен, выходной диапазон составляет от 0 до 255. Если опущен, входной диапазон автоматически вычисляется из исходных данных.

Так что вы можете использовать следующее:

gdal_translate -of GTiff -ot Byte -scale 0 65535 0 255 src_dataset dst_dataset

Почему 255 и 65535?

У вас есть 2 ^ nbits значений. Поскольку минимальное значение обычно равно 0, максимальное значение равно 2^nbits-1.

8 бит -> 2 ^ 8 = 256 значений -> 0-255
16 бит -> 2 ^ 16 = 65536 значений -> 0-65535

— pLumo
источник