Страница 1 из 2
Ошибка перепроецирования растра
Добавлено: 02 ноя 2012, 08:52
Игорь Лебедь
При перепроецировании растра выходит следующая ошибка. Исходный растр 1.rsw, выходной - 2.tif, в чем может быть дело? Версия QGIS - NextGIS 1.8.0.

- rsw1.png (67.9 КБ) 14218 просмотров
Может лучше стоит использовать какой-то специальный модуль для перепроецирования?
Re: Ошибка перепроецирования растра
Добавлено: 02 ноя 2012, 11:11
Дмитрий Барышников
Скопируйте пожалуйста строку команды, которая у вас на скриншоте закрывает диалог в эту тему.
Re: Ошибка перепроецирования растра
Добавлено: 02 ноя 2012, 11:57
Игорь Лебедь
gdalwarp -s_srs EPSG:4284 -t_srs "+proj=tmerc +lat_0=0 +lon_0=41.123456 +k=1 +x_0=123456 +y_0=-12345.678 +ellps=krass +towgs84=23.92,-141.27,-80.9,-0,0.35,0.82,-0.12 +units=m +no_defs" -r near -multi -of GTiff D:/Work/49_Kalitva/Melovaya_edit2/Melovaya_3_GK-42-7.rsw D:/Work/49_Kalitva/Melovaya_edit2/Melovaya_3.tif
Re: Ошибка перепроецирования растра
Добавлено: 02 ноя 2012, 14:10
Дмитрий Барышников
Попробуйте в эту строчку добавить -et 0.0 - там есть кнопка отредактировать.
Просто похожая ошибка всплывала
http://gis-lab.info/qa/orbview3-ortho-gdal.html (ближе к концу статьи).
Re: Ошибка перепроецирования растра
Добавлено: 02 ноя 2012, 14:44
Игорь Лебедь
Спасибо, к сожалению, не помогло - то же сообщение. Я только не уверен, что именно туда добавил -et 0.0 - перед "-r near -multi". Может, неправильно, подскажите, пожалуйста!
Re: Ошибка перепроецирования растра
Добавлено: 02 ноя 2012, 15:11
Дмитрий Барышников
Место вставки не важно. Раз ошибка та же, значит параметр не влияет. Видимо сами точки проблемные. Попробуйте уменьшит их количество. Можно вообще попробовать с тремя точками и аффинным преобразованием. Может с растром проблемы.
Re: Ошибка перепроецирования растра
Добавлено: 02 ноя 2012, 16:16
Максим Дубинин
тут скорее надо просить исходник для примера куда-нибудь выложить
Re: Ошибка перепроецирования растра
Добавлено: 02 ноя 2012, 16:48
Дмитрий Барышников
Чуть меньше 3000 GCP точек несколько пугают меня.
Re: Ошибка перепроецирования растра
Добавлено: 04 ноя 2012, 18:30
Игорь Лебедь
Максим Дубинин писал(а):тут скорее надо просить исходник для примера куда-нибудь выложить
Пожалуйста! Я правда и не знаю, что можно с ним сделать.
http://depositfiles.com/files/yf9rhl2l3
Re: Ошибка перепроецирования растра
Добавлено: 04 ноя 2012, 19:22
gimran
Уберите галочку с метода интерполяции. Если не ошибаюсь в этом ошибка.
Re: Ошибка перепроецирования растра
Добавлено: 04 ноя 2012, 19:37
Игорь Лебедь
gimran писал(а):Уберите галочку с метода интерполяции. Если не ошибаюсь в этом ошибка.

- gislab2.png (49.12 КБ) 14106 просмотров
Leide, nichts.
Re: Ошибка перепроецирования растра
Добавлено: 04 ноя 2012, 20:40
Максим Дубинин
У вас исходный растр уже в Transverse Mercator.
Код: Выделить всё
c:\temp>gdalinfo Melovaya_3_GK-42-7.rsw
Driver: RMF/Raster Matrix Format
Files: Melovaya_3_GK-42-7.rsw
Size is 11421, 2624
Coordinate System is:
PROJCS["unnamed",
GEOGCS["Pulkovo 1942",
DATUM["Pulkovo_1942",
SPHEROID["Krassowsky 1940",6378245,298.2999999999985,
AUTHORITY["EPSG","7024"]],
AUTHORITY["EPSG","6284"]],
PRIMEM["Greenwich",0],
UNIT["degree",0.0174532925199433],
AUTHORITY["EPSG","4284"]],
PROJECTION["Transverse_Mercator"],
PARAMETER["latitude_of_origin",0],
PARAMETER["central_meridian",39],
PARAMETER["scale_factor",1],
PARAMETER["false_easting",7500000],
PARAMETER["false_northing",0],
UNIT["metre",1,
AUTHORITY["EPSG","9001"]]]
Origin = (7574502.710000000000000,5482427.790000000000000)
Pixel Size = (0.674542682926829,-0.674542682926829)
Metadata:
AREA_OR_POINT=Area
Image Structure Metadata:
INTERLEAVE=PIXEL
Corner Coordinates:
Upper Left ( 7574502.710, 5482427.790) ( 40d 1'40.34"E, 49d28' 9.57"N)
Lower Left ( 7574502.710, 5480657.790) ( 40d 1'39.14"E, 49d27'12.29"N)
Upper Right ( 7582206.662, 5482427.790) ( 40d 8' 2.90"E, 49d28' 6.00"N)
Lower Right ( 7582206.662, 5480657.790) ( 40d 8' 1.58"E, 49d27' 8.72"N)
Center ( 7578354.686, 5481542.790) ( 40d 4'50.99"E, 49d27'39.19"N)
Band 1 Block=11421x1 Type=Byte, ColorInterp=Red
Band 2 Block=11421x1 Type=Byte, ColorInterp=Green
Band 3 Block=11421x1 Type=Byte, ColorInterp=Blue
Вы же указываете ему:
-s_srs EPSG:4284
Что никак не соответствует действительности.
Вот вам и ответ на ваш вопрос, уберите -s_srs, зачем он, если у вас СК прописана уже в растре и тем более, это не 4284.
Re: Ошибка перепроецирования растра
Добавлено: 04 ноя 2012, 21:39
Игорь Лебедь
Максим Дубинин писал(а):У вас исходный растр уже в Transverse Mercator.
Код: Выделить всё
c:\temp>gdalinfo Melovaya_3_GK-42-7.rsw
Driver: RMF/Raster Matrix Format
Files: Melovaya_3_GK-42-7.rsw
Size is 11421, 2624
Coordinate System is:
PROJCS["unnamed",
GEOGCS["Pulkovo 1942",
DATUM["Pulkovo_1942",
SPHEROID["Krassowsky 1940",6378245,298.2999999999985,
AUTHORITY["EPSG","7024"]],
AUTHORITY["EPSG","6284"]],
PRIMEM["Greenwich",0],
UNIT["degree",0.0174532925199433],
AUTHORITY["EPSG","4284"]],
PROJECTION["Transverse_Mercator"],
PARAMETER["latitude_of_origin",0],
PARAMETER["central_meridian",39],
PARAMETER["scale_factor",1],
PARAMETER["false_easting",7500000],
PARAMETER["false_northing",0],
UNIT["metre",1,
AUTHORITY["EPSG","9001"]]]
Origin = (7574502.710000000000000,5482427.790000000000000)
Pixel Size = (0.674542682926829,-0.674542682926829)
Metadata:
AREA_OR_POINT=Area
Image Structure Metadata:
INTERLEAVE=PIXEL
Corner Coordinates:
Upper Left ( 7574502.710, 5482427.790) ( 40d 1'40.34"E, 49d28' 9.57"N)
Lower Left ( 7574502.710, 5480657.790) ( 40d 1'39.14"E, 49d27'12.29"N)
Upper Right ( 7582206.662, 5482427.790) ( 40d 8' 2.90"E, 49d28' 6.00"N)
Lower Right ( 7582206.662, 5480657.790) ( 40d 8' 1.58"E, 49d27' 8.72"N)
Center ( 7578354.686, 5481542.790) ( 40d 4'50.99"E, 49d27'39.19"N)
Band 1 Block=11421x1 Type=Byte, ColorInterp=Red
Band 2 Block=11421x1 Type=Byte, ColorInterp=Green
Band 3 Block=11421x1 Type=Byte, ColorInterp=Blue
Вы же указываете ему:
-s_srs EPSG:4284
Что никак не соответствует действительности.
Вот вам и ответ на ваш вопрос, уберите -s_srs, зачем он, если у вас СК прописана уже в растре и тем более, это не 4284.
Спасибо огромное! Кажется, разобрался. 4284 - это система проекта по умолчанию, прописалась и сюда, получается нормально либо если указать 28407 - действительную систему исходного растра, либо если снять галку "исходная проекция". Расшифруйте пожалуйста, что значит "-s_srs"?
Re: Ошибка перепроецирования растра
Добавлено: 04 ноя 2012, 21:50
Максим Дубинин
да, либо убрать исходную совсем, либо указать, но правильную
-a_srs, -s_srs, -t_srs
Re: Ошибка перепроецирования растра
Добавлено: 04 ноя 2012, 21:55
Игорь Лебедь
Ясно, спасибо. То есть перепроецирование в QGIS основано на ogr2ogr, да?