Где взять карту?

[wowka1319]

Спасибо ericsson за развернутое объяснение. Короче эталонов нет, всё относительно. Однако возникают (чисто из любопытства) два вопроса:
1. Откуда взялась байка про гринвечский меридиан = нулевой меридиан? может существует какой-то старый ранее популярный и теперь неактуальный датум, где это истина.
2. Что помешало комиссии (или кто там) в WGS84 назначить за нулевой меридиан гривичский? Ведь это на мой взгляд вообще не имеет какого-то практического значения. Почему именно +5.31"?

Ну это так, из любопытства. А теперь по делу:
Я нашел нужный мне источник (откуда ArcMap у нас на работе качает тайлы). Тут (http://services.arcgisonline.com/ArcGIS ... /MapServer) его описание, по которому я слепил TMS драйвер:

<GDAL_WMS>
<Service name="TMS">
<ServerUrl>http://services.arcgisonline.com/ArcGIS ... /ServerUrl>>
</Service>
<DataWindow>
<UpperLeftX>-20037508,342787</UpperLeftX>
<UpperLeftY>20037508,342787</UpperLeftY>
<LowerRightX>20037508,342787</LowerRightX>
<LowerRightY>-20037508,342787</LowerRightY>
<TileLevel>12</TileLevel>
?<TileCountX>1</TileCountX>
?<TileCountY>1</TileCountY>
<YOrigin>top</YOrigin>
</DataWindow>
<Projection>EPSG:3857</Projection>
<BlockSizeX>256</BlockSizeX>
<BlockSizeY>256</BlockSizeY>
<BandsCount>3</BandsCount>
?<MaxConnections>10</MaxConnections>
<Cache />
</GDAL_WMS>

Там где стоят вопросы - это те строки, значение которых я не понял и поэтому просто скопировал из исходного xml (для другой карты).
Отсюда вопросы:
3. Что это за строки (особенно волнуют первые две, ибо полагаю что последняя указывает максимальное число попыток коннекта с сервером). Дословный перевод мне их понятен, но вот суть...

4. На каждый уровень детализации указан параметр "Resolution", указывающий, как я понял, сколько метров по длине и высоте охватывает пиксель тайла. Но возникает недоумевание: прямоугольная проекция сферического меркатора EPSG:3857 (он же EPSG:900913) - есть квадрат, где меридианы параллельны. А раз они параллельны (следовательно искажены), то чем ближе к полюсам тем разряженнее растояние между ними, а т.к. число пикселей неизменно, то ширина пикселя, двигаясь к полюсам, будет охватывать меньшее расстояние в метрах. Другими словами разрешение (оно же Resolution) не может быть постоянным... Тогда что же означает "Resolution"? Что со мной не так?

[ericsson]

Гринвичский меридиан, проходящий через меридианны круг Гринвичской обсерватории, и был международным нулевым очень долгое время, с 1884 года, когда в Вашингтоне прошла Международная меридианальная конференция.
А нынешний нулевой меридиан для WGS84 - это не Гринвичский, а IERS Reference meridian. И сдвиг этот произошел из-за американцев, которые построили первую спутниковую навигационную систему. Далее все случилось так: они взяли координаты своей опорной станции в американском датуме NAD27 (который является локальным, а не общемировым, аналогично советскому Пулково-1942) и при создании общемирового датума WGS просто приравняли эти самые координаты (то есть координаты WGS равны координатам NAD27 на месте референсной станции Университета Джона Хопкинса, где это все начинали разрабатывать). А потом, когда стали определять средствами этой системы координаты Гринвичского меридиана, обнаружился сдвиг. Далее его несколько раз уточняли (пока модель WGS совершенствовалась) но в итоге он оказался таким, как есть - чуть больше пяти угловых секунд, около ста метров.

TileCount - параметр, который используется для вычисления размера всего источника. Ставьте ноль спокойно...
А MaxConnections это не "MaxRetry". Он задает число одновременных подключений к серверу.

[ericsson]

Собственно, я сильно удивлен, что в NAD27 через океан от его зоны применения набежало всего сто метров...
Так что известная проблема с туристами, которые привязывают карты, сделанные в СК-42 или СК-95 в WGS84 без преобразования - не уникальна и имеет глубокие исторические корни

[trir]

offtop: а почему тогда в proj - нет IERS Reference meridian, а только greenwich

[ericsson]

Видимо, в PROJ4 есть как раз IRM, который назван ошибочно Greenwich, а вот "старого" Greenwich там как раз нет.

[ericsson]

А, забыл ответить на четвертый вопрос.
"Разрешение" в WMS/TMS (и не только) - это количество единиц проекции в пикселе, а не количество настоящих метров в нем. Так что если речь о Mercator, то все нормально - "разрешение" с севера на юг по горизонтали не меняется, меняется соотношение этого разрешения и реальных "метров" земной поверхности. На экваторе оно будет единицей, на полюсе стремиться к бесконечности. Так что считайте, что для Меркатора это разрешение в метрах на пиксель на экваторе.

[ericsson]

Тут еще с эталонами есть такая петрушка (привет соседней теме...): земная кора ползет, при том в разных направлениях, и все эти геодезические пункты на разных континентах чудесно разбредаются на заметное расстояние даже за год. Потому всякого рода ECEF, ITRS, ITRF - это штуки по большей части научного свойства, а не практического. Я такой сложной математикой уже не владею, но когда возился с RTKlib, пришлось припомнить кое-чего... Вот наших (и вообще любых) геодезистов в институтах начинать бы учить с RTKlib, а не с кнопконажимательного софта от Leica и Trimble.

[wowka1319]

Действительно, в SAS.Планета разрешение на экваторе совпадет с Resolution соответствующего уровня детализации (кстати в SAS эти уровни начинаются с 1, а не с 0. занятно...).
То есть смысл в локальных датумах как раз и заключается в том чтобы быть менее зависимыми от сдвигов земных плит?

ericsson, TileCount = 0 НЕ работает (результат на картинке). Этот параметр реально на что-то влияет. Я ставлю = 1. Так что же это?

[ericsson]

Смысл локальных датумов раньше был в том, что они лучше аппроксимировали поверхность Земли на территории какой-то страны. Когда еще не было ни единой гравитационной модели (геоида), ни общемировых датумов. Все это строилось на эллипсоидах, которые подгонялись под поверхность в отдельно взятой стране, по отдельно взятой геодезической триангуляционной сети. А про то, какова реальная скорость движения литосферных плит, тогда никто толком не знал, и не задумывались над этим - средств для регистрации этого не было. И считали все не на компьютере, потому когда местная система координат основана на каком-нибудь эллипсоиде Бесселя столетней давности - так и не страшно, потому что все ее эталонные точки собраны на ограниченной территории, и на движение плит и прочее какому-нибудь коммунальному хозяйству города или области - начхать.

TileCount - число тайлов и есть... если работает с единицей - ставьте единицу. Я не имею понятия, зачем этот параметр придуман - при обращении к WMS/TMS другими средствами он не используется, даже экстент (границы) знать не обязательно (ну выдаст сервер 404 на несуществующий тайл, и что?). Это что-то разработчики GDAL перемудрили и сами не знают, зачем (в документации написано что default value = 0, что указывает на то, что этим никто не пользуется).
Посмотрел свои старые минидрайверы, которыми у норвегов таскал их геологические карты - там у меня единица тоже стоит.

[wowka1319]

Как теперь перепроецировать полученный геотиф EPSG:3857 (сферический меркатор) в WGS84 - не могу разобрать в параметрах gdalwarp?
И еще: можно ли сделать такое же перепроецирование при таких же условиях, но только если исходная карта не в геотифе, а в JPG/PNG с известными bbox координатами? если да, то как?

[ericsson]

Первые же два параметра в списке http://www.gdal.org/gdalwarp.html вам в помощь.
Если исходная карта с известными координатами - ее сначала надо привязать.

[wowka1319]

Естественно я уже читал документацию, но не смог разобраться с параметрами (о чем уже писал).
Приведенный там пример мне не понятен:
gdalwarp -t_srs '+proj=utm +zone=11 +datum=WGS84' raw_spot.tif utm11.tif

а именно не понятен формат задания проекции.
Кто уже пользовался этой утилитой, подскажите напишите запрос gdalwarp, чтобы из geotiff в EPSG:385 получить geotiff в WGS84.

[Petruxin]

Если Вам просто надо перепроецировать, то не обязательно так глубоко копать. Скачайте QGIS тут например , установите и нажатием 2ух кнопок перепроецируйте что хотите.

[ericsson]

Самостоятельно разобраться, что такое srs definition не пробовали?
Пример:
http://spatialreference.org/ref/epsg/4326/
http://spatialreference.org/ref/epsg/4326/proj4/

[Дмитрий Барышников]

В -s_srs и -t_srs srs не обязательно писать строчку в формате PROJ.
Вот так тоже можно:
gdalwarp -t_srs 'EPSG:4326' raw_spot.tif utm11.tif