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