Привязка ECW

Вопросы по нескольким пакетам сразу, или вопросы, которые непонятно к какой ГИС отнести
Ответить
krnark
Новоприбывший
Сообщения: 3
Зарегистрирован: 14 апр 2008, 23:43
Репутация: 0

Привязка ECW

Сообщение krnark » 15 апр 2008, 00:36

Вопрос наверное чайниковский, я в этой области и есть чайник, хотя что такое к примеру проекция и датум представление имею.

У меня такая задача, сделать из снимков GoogleEarth один файл (уровня города и окрестностей) с привязкой. Решил сделать в формате ECW, там координатная иформация в самом файле хранится все ГИС его понимают, да бесплатные просмотрщики есть. В качестве исходного материала у меня есть куски сделанные из кеша GoogleEarth с привязкой в формате map OziExplorer (там 9 точек прописаны). В качестве кодировщика ECW я использовал ER Mapper 7.1 Т.к. формат map Ozi напрямую ER Mapper не понимает я превел его в GeoTIFF программой Global Mapper 9. В ней указывал сделать к tiff еще отдельно файлы twf, там данные такого вида

0.59716445006055929
0.00000000000000000
0.00000000000000000
-0.59547252617757873
3833470.14047288710000000
6332835.74489106980000000

как я понял это размер пиксела в метрах и координаты верхнего левого угла (кстати скажите в чем это координаты, я не понял)
Потом из этих tiff ER Mapper атоматически собрал мозаику в один ECW с заголовочным файлом вида

DatasetHeader Begin
Version = "6.0"
Name = "poltava2.ers"
LastUpdated = Sun Apr 13 21:34:05 GMT 2008
DataFile = "poltava2.ecw"
DataSetType = Translated
DataType = Raster
ByteOrder = LSBFirst
CoordinateSpace Begin
Datum = "WGS84"
Projection = "MRAFRICA"
CoordinateType = EN
Rotation = 0:0:0.0
CoordinateSpace End
RasterInfo Begin
CellType = Unsigned8BitInteger
CellInfo Begin
Xdimension = 0.59716430835621026
Ydimension = 0.59547252617757473
CellInfo End
NrOfLines = 44295
NrOfCellsPerLine = 45056
RegistrationCoord Begin
Eastings = 3833469.83633191650000000
Northings = 6358015.26831095950000000
RegistrationCoord End
NrOfBands = 3
RasterInfo End
DatasetHeader End

В результате координаты в ER Mapper вычисляются правильно (совпадение с GoogleEarth до сотой секунды) а расстояния не правильные, причем Global Mapper по этому файлу все считает правильно. Стал разбиратся, меряю (в ER Mapper) объект длина которого известна 1км, получается в точках что-то около 2700 в метрах 1600, что и понятно 2700*0.59 Ладно, померял при помощи Global Mapper (который правильно все меряет) файл по горизонтали вертикали и всчитал
Xdimension=0,38543146306818181
Ydimension= 0.38602551077999774
И поставил координаты верхнего левого угла в долг/шир вот так
RegistrationCoord Begin
Latitude = 49:40:39.00
Longitude = 34:26:11.93
RegistrationCoord End

В результате меряет теперь ER Mapper правильно, а координаты вычисляет неверно (в 2 минуты расхождение). Я так понимаю что он координаты вычсляет на основе координат верхнего левого угла и размера пиксела. Возвращаю старые размеры пиксела координаты становятся правильными, размеры соответственно нет. Скажите что я делаю не правильно.

Аватара пользователя
Максим Дубинин
MindingMyOwnBusiness
Сообщения: 9129
Зарегистрирован: 06 окт 2003, 20:20
Репутация: 748
Ваше звание: NextGIS
Откуда: Москва
Контактная информация:

Сообщение Максим Дубинин » 15 апр 2008, 03:46

как все запутано
стал разбиратся, меряю (в ER Mapper) объект длина которого известна 1км, получается в точках что-то около 2700
Реальная длина и длина измеренная в Меркаторе на сфере (исходной СК данных гугл) две большие разницы
viewtopic.php?t=1774
пристегивайтесь, турбулентность прямо по курсу

krnark
Новоприбывший
Сообщения: 3
Зарегистрирован: 14 апр 2008, 23:43
Репутация: 0

Сообщение krnark » 15 апр 2008, 06:47

Это все хорошо, однако у меня перед глазами этот самый Global Mapper который по этому файлу меряет расстояния совершенно правильно, стоит у него тот же mercator. Опять же на квадрате в 17км даже при самых диких нестыковках по проекциям, никак не можно по расстояниям врать в полтора раза. Тут у меня где-то явно грубая ошибка но не пойму в чем.

И подскажите что это за система координат, которая
Eastings = 3833469.83633191650000000
Northings = 6358015.26831095950000000

krnark
Новоприбывший
Сообщения: 3
Зарегистрирован: 14 апр 2008, 23:43
Репутация: 0

Сообщение krnark » 15 апр 2008, 22:38

Проблема решается преобразованием в UTM, но как-то мне это не очень нравится, появляются дополнительные погрешности. Опять же у этого Global Mapper без всяких манипуляций полная тождественность с GoogleEarth. Какой-то он сильно умный... или наооборот.

Аватара пользователя
Максим Дубинин
MindingMyOwnBusiness
Сообщения: 9129
Зарегистрирован: 06 окт 2003, 20:20
Репутация: 748
Ваше звание: NextGIS
Откуда: Москва
Контактная информация:

Сообщение Максим Дубинин » 16 апр 2008, 04:35

Проблема решается преобразованием в UTM
естественно решается, просто считать начинаете в проекции более корректно передающей длины.
Только это не проблема, одни проекции для одних целей, другие для других

а Global Mapper ваш видимо действительно думает, что самый умный и производит вычисления например по расстояниям большого круга, зная ск данных (в данном случае меркатор) легко получить географические координаты. Если это так, то это на мой взгляд проблема, так как я не хочу, чтобы программа принимала за меня решение как считать, наверняка отключаться должно в настройках.
пристегивайтесь, турбулентность прямо по курсу

Ответить

Вернуться в «Общий - ПО»

Кто сейчас на конференции

Сейчас этот форум просматривают: нет зарегистрированных пользователей и 3 гостя