Форумы GIS-Lab.info

Приветствую всех на форуме!
Программу arcmap установил первый раз, но весь день безуспешно искал решение своей проблемы(
Имеется карта высот DTM разрядностью 32 бита. В Photoshop она не открывается.
Я решил открыть ее в arcmap и сохранить в растре с понижением битности до 24, чтобы photoshop ее понял.
При экспорте ставлю Tiff и 24 бита, файл сохраняется, но в 8 битах. ДА и на превью в arcmap сразу видно, что он 8бит, но в свойствах пишет, что это 32 бита.
Как мне решить проблему и сохранить в файл в высокой битности?
Заранее благодарен!

ngektor, можно воспользоваться утилитой gdal_translate из командной строки. Ключ "-ot" позволяет указать желаемую "битность" выходного растра.

antonv писал(а): ↑

16 фев 2021, 10:54

можно воспользоваться утилитой gdal_translate из командной строки. Ключ "-ot" позволяет указать желаемую "битность" выходного растра

Спасибо за ответ!
Я так полагаю, что там надо разбираться в программировании?))))
Нет ли у вас возможности перевести файл в нужную битность? Что мне это будет стоить?

Выложите карту высот DTM

lam писал(а): ↑

16 фев 2021, 14:20

Выложите карту высот DTM

Приветствую! вот ссылка на 32битный DTM, но весит он 1,8гб...
https://planetarymaps.usgs.gov/mosaic/m ... lend40.tif

ngektor, нет, в программировании можно не разбираться, достаточно набрать команду в командной строке. Примеры использования лежат здесь.

antonv писал(а): ↑

16 фев 2021, 17:30

ngektor, нет, в программировании можно не разбираться, достаточно набрать команду в командной строке. Примеры использования лежат здесь.

Приветствую. Все поставил, командой сделал перевод, файл переконвертировался
(C:\>gdal_translate -ot UInt16 C:\Users\kolyan\Desktop\Марс2\JEZ_hirise_soc_006.tif C:\Users\kolyan\Desktop\Марс2\JEZ_hirise_soc_00666.tif
Input file size is 21400, 21488
Warning 1: for band 1, nodata value has been clamped to 0, the original value being out of range.
0...10...20...30...40...50...60...70...80...90...100 - done.),
стал меньше на 30%, но все равно не читается, теперь нет карты ни в arcmap ни в ФШ. Причем в ФШ теперь пишет, что стало 16 бит ....

Не очень понятно, что вы делаете - вы просили 16 бит, вам их gdal и выдал.

Что за "карта высот DTM разрядностью 32 бита", откуда такая информация? Это картинка RGBс с альфа-каналом, или реально DTM? Если DTM, то какие там значения, int или float? Если картинка, то что вы с ней делать собрались, альфа-канал отрезать?

Для начала получите и покажите выхлоп команды gdal_info, поскольку качать 2Гб никто не будет.

gamm писал(а): ↑

17 фев 2021, 17:56

Не очень понятно, что вы делаете - вы просили 16 бит, вам их gdal и выдал.

Что за "карта высот DTM разрядностью 32 бита", откуда такая информация? Это картинка RGBс с альфа-каналом, или реально DTM? Если DTM, то какие там значения, int или float? Если картинка, то что вы с ней делать собрались, альфа-канал отрезать?

Для начала получите и покажите выхлоп команды gdal_info, поскольку качать 2Гб никто не будет.

Информация о 32 битах была в исходниках сайта и по информации из arcmap, там указано в свойствах файла.
https://drive.google.com/file/d/1IKoDDT ... sp=sharing

Данные, которые я получил от gralinfo:
C:\>gdalinfo C:\Users\kolyan\Desktop\Марс2\JEZ_hirise_soc_006.tif
Driver: GTiff/GeoTIFF
Files: C:\Users\kolyan\Desktop\╨Ь╨░╤А╤Б2\JEZ_hirise_soc_006.tif
C:\Users\kolyan\Desktop\╨Ь╨░╤А╤Б2\JEZ_hirise_soc_006.tif.ovr
Size is 21400, 21488
Coordinate System is:
PROJCRS["Equirectangular Mars 2000 Sphere IAU",
BASEGEOGCRS["D_Mars_2000_Sphere",
DATUM["Mars_2000_(Sphere)",
ELLIPSOID["Mars_2000_Sphere_IAU",3396190,0,
LENGTHUNIT["metre",1,
ID["EPSG",9001]]]],
PRIMEM["Reference_Meridian",0,
ANGLEUNIT["degree",0.0174532925199433,
ID["EPSG",9122]]]],
CONVERSION["Equidistant Cylindrical",
METHOD["Equidistant Cylindrical",
ID["EPSG",1028]],
PARAMETER["Latitude of 1st standard parallel",0,
ANGLEUNIT["degree",0.0174532925199433],
ID["EPSG",8823]],
PARAMETER["Longitude of natural origin",0,
ANGLEUNIT["degree",0.0174532925199433],
ID["EPSG",8802]],
PARAMETER["False easting",0,
LENGTHUNIT["metre",1],
ID["EPSG",8806]],
PARAMETER["False northing",0,
LENGTHUNIT["metre",1],
ID["EPSG",8807]]],
CS[Cartesian,2],
AXIS["easting",east,
ORDER[1],
LENGTHUNIT["metre",1,
ID["EPSG",9001]]],
AXIS["northing",north,
ORDER[2],
LENGTHUNIT["metre",1,
ID["EPSG",9001]]]]
Data axis to CRS axis mapping: 1,2
Origin = (4577366.000000000000000,1106618.000000000000000)
Pixel Size = (1.000000000000000,-1.000000000000000)
Metadata:
AREA_OR_POINT=Area
Image Structure Metadata:
INTERLEAVE=BAND
Corner Coordinates:
Upper Left ( 4577366.000, 1106618.000) ( 77d13'22.56"E, 18d40' 9.53"N)
Lower Left ( 4577366.000, 1085130.000) ( 77d13'22.56"E, 18d18'24.48"N)
Upper Right ( 4598766.000, 1106618.000) ( 77d35' 2.27"E, 18d40' 9.53"N)
Lower Right ( 4598766.000, 1085130.000) ( 77d35' 2.27"E, 18d18'24.48"N)
Center ( 4588066.000, 1095874.000) ( 77d24'12.41"E, 18d29'17.01"N)
Band 1 Block=21400x1 Type=Float32, ColorInterp=Gray
Min=-2665.371 Max=-1639.885
Minimum=-2665.371, Maximum=-1639.885, Mean=-2419.942, StdDev=215.685
NoData Value=-3.40282265508890445e+38
Overviews: 10700x10744, 5350x5372, 2675x2686, 1338x1343, 669x672, 335x336, 168x168
Metadata:
STATISTICS_COVARIANCES=46519.85400828344
STATISTICS_MAXIMUM=-1639.8854980469
STATISTICS_MEAN=-2419.9421734083
STATISTICS_MINIMUM=-2665.37109375
STATISTICS_SKIPFACTORX=1
STATISTICS_SKIPFACTORY=1
STATISTICS_STDDEV=215.68461699501

Мне нужен файл с плавной градацией серого для высокоточной фрезерной обработки материала. Если брать 8бит, то это всего 256 оттенков серого и будут видны ступеньки.

Заархивировал и положил на свой диск, теперь весит 380мб
https://drive.google.com/file/d/1T30psd ... sp=sharing

И вот как выглядит файл в arcmap:
https://drive.google.com/file/d/1vxs9-J ... sp=sharing

Для начала уберите портянку под кат

Что вы делаете, понять невозможно - зачем засовывать данные в фотошоп вместо чпу ???
В фотошопе разрядность или 8 бит, или 3 канала по 8 бит.
У вас там тип данных float, и все значения отрицательные. Если нужны целые, задайте выходной тип int32 и определите свое nodata. Если нужны положительные значения, отнимите минимум, в выхлопе он есть. Как это сделать, читайте в доках gdal.

gamm писал(а): ↑

17 фев 2021, 18:33

Что вы делаете, понять невозможно - зачем засовывать данные в фотошоп вместо чпу ???

мне нужна не вся карта, а определённая область.

gamm писал(а): ↑

17 фев 2021, 18:33

У вас там тип данных float, и все значения отрицательные. Если нужны целые, задайте выходной тип int32 и определите свое nodata. Если нужны положительные значения, отнимите минимум, в выхлопе он есть. Как это сделать, читайте в доках gdal.

Спасибо.

ngektor писал(а): ↑

17 фев 2021, 18:42

мне нужна не вся карта, а определённая область.

Выложите шейп файл полигона, в пределах которого необходимо обрезать DTM.
Или покажите на скриншоте.

Напишите свою почту.

lam писал(а): ↑

17 фев 2021, 22:49

Напишите свою почту

Моя почта 986666@inbox.ru

ngektor писал(а): ↑

17 фев 2021, 18:42

мне нужна не вся карта, а определённая область.

после фотошопа все равно больше 8 бит не получите.. вырезайте в Арке или прямо gdal-ом