Привязка ECW
-
- Новоприбывший
- Сообщения: 3
- Зарегистрирован: 14 апр 2008, 23:43
- Репутация: 0
Привязка ECW
Вопрос наверное чайниковский, я в этой области и есть чайник, хотя что такое к примеру проекция и датум представление имею.
У меня такая задача, сделать из снимков 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 минуты расхождение). Я так понимаю что он координаты вычсляет на основе координат верхнего левого угла и размера пиксела. Возвращаю старые размеры пиксела координаты становятся правильными, размеры соответственно нет. Скажите что я делаю не правильно.
У меня такая задача, сделать из снимков 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
- Откуда: Москва
- Контактная информация:
как все запутано
viewtopic.php?t=1774
Реальная длина и длина измеренная в Меркаторе на сфере (исходной СК данных гугл) две большие разницыстал разбиратся, меряю (в ER Mapper) объект длина которого известна 1км, получается в точках что-то около 2700
viewtopic.php?t=1774
пристегивайтесь, турбулентность прямо по курсу
-
- Новоприбывший
- Сообщения: 3
- Зарегистрирован: 14 апр 2008, 23:43
- Репутация: 0
Это все хорошо, однако у меня перед глазами этот самый Global Mapper который по этому файлу меряет расстояния совершенно правильно, стоит у него тот же mercator. Опять же на квадрате в 17км даже при самых диких нестыковках по проекциям, никак не можно по расстояниям врать в полтора раза. Тут у меня где-то явно грубая ошибка но не пойму в чем.
И подскажите что это за система координат, которая
Eastings = 3833469.83633191650000000
Northings = 6358015.26831095950000000
И подскажите что это за система координат, которая
Eastings = 3833469.83633191650000000
Northings = 6358015.26831095950000000
-
- Новоприбывший
- Сообщения: 3
- Зарегистрирован: 14 апр 2008, 23:43
- Репутация: 0
- Максим Дубинин
- MindingMyOwnBusiness
- Сообщения: 9129
- Зарегистрирован: 06 окт 2003, 20:20
- Репутация: 748
- Ваше звание: NextGIS
- Откуда: Москва
- Контактная информация:
естественно решается, просто считать начинаете в проекции более корректно передающей длины.Проблема решается преобразованием в UTM
Только это не проблема, одни проекции для одних целей, другие для других
а Global Mapper ваш видимо действительно думает, что самый умный и производит вычисления например по расстояниям большого круга, зная ск данных (в данном случае меркатор) легко получить географические координаты. Если это так, то это на мой взгляд проблема, так как я не хочу, чтобы программа принимала за меня решение как считать, наверняка отключаться должно в настройках.
пристегивайтесь, турбулентность прямо по курсу
Кто сейчас на конференции
Сейчас этот форум просматривают: нет зарегистрированных пользователей и 3 гостя