Собственная Wikimapia

[Shaggy]

Да, снимки становятся бесплатными на NASA и USGS, но провести данную работу по моей теме ни кто не сможет, т.к. всем лень, а если ещё и надо тратить деньги, так это вообще ни в какие ворота ) лучше в ресторане прожрать.
Так, что хотелось бы отметить своё авторство и свои трудозатраты.

ImageMagick видно мощная утилита, надо разбираться. Если есть готовые наработки - прошу поделиться.
ВЗ на каждом тайле наверно перебор. Думается надо сделать один ВЗ на3 тайла (потом число оптимизировать) потом например повторение через три тайла по горизонтали и по вертикали.

Вот если бы такой скриптик где надо менять только циферки на сколько тайлов надпись растягивать и какие делать пропуски...

[paleogis]

Да, снимки становятся бесплатными на NASA и USGS, но провести данную работу по моей теме ни кто не сможет, т.к. всем лень, а если ещё и надо тратить деньги, так это вообще ни в какие ворота ) лучше в ресторане прожрать.

Никто Вам ничего не должен. Ваша тема, Ваши проблемы. Мне тоже с USGS приходилось покупал съемку за деньги, и заниматься копирастией, даже в голову не приходило, тем более я как и Вы не правообладатель съемки. На Вашем месте я бы вообще с копирастией не заморачивался. Будьте проще, и люди к Вам потянутся.

[ericsson]

"Выполнить данную работу" - это скачать снимки с серверов NASA? И при чем тут деньги - снимки уже бесплатно там доступны, если вы их купили.
Получается, не "никто не может", а "никому не надо". Что должно быть прозрачным намеком на то, что свои watermark на эти снимки лепить - бессмысленно, они никому не нужны. Скачать готовый тайловый источник ведь не проще, чем исходный снимок.

[Shaggy]

Хорошо, хорошо, заклевали сразу )))

Я вообще мало, что понимаю в координатах и проекциях стремлюсь быстрее быстрее перейти прикладному использованию ГИС в связи с этим наверно делаю много ошибок. Они теперь стали проявляться при работе с более большими площадями чем я делал до этого.

Привязываю снимок с USGS с помощью QGIS по координатам строений взятых с карты Google с помощью SAS.Планета. В начале привязки растра выбираю систему координат WGS84 EPSG 4326, хотя SAS.Планета работает в меркатор, но в QGIS этих меркаторов хоть пруд пруди, какой на самом деле надо выбирать?

Если координаты углов снимка с USGS они какой системе координат\проекции?

В QGIS целевую систему координат выбираю тоже WGS84 EPSG 4326.

Режу на тайлы ГлобалМеппером там ставлю систему координат WGS84 проекцию Меркатор. Ощущение, что делаю не правильно всё время было, но для маленьких частей снимка всё хорошо работало, теперь при привязки большого снимка явные большие ошибки.

Как правильно сделать?

[ericsson]

Раз вы понимаете, в чем проблема (в том, что вы делаете какие-то шаги, не понимая до конца их смысл) - значит она уже частично решена.

Если координаты углов снимка с USGS они какой системе координат\проекции?

Могут быть в какой угодно. Для начала, вопросы:
- вы уверены, что у этих снимков уже нет привязки?
- если ее действительно нет, то что это конкретно за данные?
(Я могу предположить, что это сканы исторических снимков Corona, в таком случае они не приведены к какой-либо проекции.)

В QGIS целевую систему координат выбираю тоже WGS84 EPSG 4326.

Если мое предположение выше верно, то в этом случае их можно привязывать в какой угодно системе координат. Для выбора есть два критерия:
- система должны быть близка к той, которая используется в снимках (а если там никакой нет, то более-менее соответствовать ей по геометрии, что точно не так в случае выбора EPSG:4326, растягивающей изображение вдоль параллелей), чтобы потери информации были минимальны;
- система может быть выбрана той, в которую вы потом собираетесь экспортировать данные, если это не противоречит предыдущему пункту.
Теперь еще вопрос: в какой системе координат вы собираетесь показывать эти данные? Скорее всего, в "Меркатор на сфере", как у Гугла (EPSG:900913 она же EPSG:3857). Если так и есть, то и привязывайте сразу в ней. В ней же и будете экспортировать в TMS.

[Shaggy]

Привязываю снимки Declassified Data раздел Declass 2.
Снимок в формате tiff, но кадр разделён на 2 части с перехлёстом. QGIS говорит, что tiff не содержит геоданные.

На сколько я понял реальная точка на поверхности земли в разных проекциях имеет разные координаты.

Каждый кадр (состоящий из 2 tiff файлов) комплектуется файлом с методанными следующего вида
Data Set Attribute Attribute Value
Entity ID DZB00401800040H009001
Acquisition Date 1965/05/30
Mission 4018
Frame 9
Image Type Black and White
Camera Type KH-7 High Resolution Surveillance
Camera Resolution 2 to 4 feet
Film Type 9 inch by variable
Generation 3
Segment Count 1
Polarity Negative
Operations Number 00040
Down Load Available Yes
Center Latitude 46°36'57.60"N
Center Longitude 74°32'27.60"E
NW Corner Lat 46°49'55.20"N
NW Corner Long 74°27'00.00"E
NE Corner Lat 46°48'00.00"N
NE Corner Long 74°43'55.20"E
SE Corner Lat 46°24'00.00"N
SE Corner Long 74°37'55.20"E
SW Corner Lat 46°25'55.20"N
SW Corner Long 74°21'00.00"E
Center Latitude dec 46.616
Center Longitude dec 74.541
NW Corner Lat dec 46.832
NW Corner Long dec 74.45
NE Corner Lat dec 46.8
NE Corner Long dec 74.732
SE Corner Lat dec 46.4
SE Corner Long dec 74.632

О том, в какой проекции эти координаты не написано, возможно, если в earthexplorer отображается место расположение выбранных кадров на карте гугл то это координаты, судя по вашему предыдущему сообщению меркатор на сфере, хорошо, я попробую.

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

Поэтому я взял один файл и в том месте, где я хорошо знаю строения наметил 5 точек и привязал описанным выше способом. Всё отлично получилось, но для маленькой площади 5х3 км. При привязке по этим же 5ти точкам всего файла я решил, что как то не точно привязалось ближе к краю файла.

Тогда я к своим точкам добавил 2 точки из файла с методанными - результат плачевный, даже мои точки сместились на сотни метров от объектов. Поэтому я понял, что без понятия проекции я не могу дальше двигаться. В связи с этим вопросы:

1. Т.к. я для привязки в QGIS беру координаты строений с карты гугл при помощи SAS.Планета то СК WGS84 т.к. всё американское? Проекция меркатор на сфере?

2. Выходная СК тоже должна быть WGS84 меркатор на сфере? Для использования в SAS.Планета и Геомиксере?

3. Какую проекцию ставить в ГлобалМеппер (использую его только для нарезки тайлов)? В нём есть просто проекция меркатор, ни чего про сферу нет. И он явно изменяет растровый файл (видно не вооружённым взглядом).

4. Зачем каждая программа меняет проекцию? Нельзя ли все пересчёты растра сделать при привязке, а потом его просто резать на тайлы.

Старался изложить всё наиболее подробно. Если, что то упустил - скажите, допишу.

[ericsson]

Ну то есть я угадал.
Значится так:
- те координаты, которые обозначены в метаданных снимка (то что вы там процитировали) - они, конечно, в WGS84, но они очень приблизительные, их использовать для привязки вообще нельзя, я встречал данные, которые по этим координатам были смещены на много километров, они просто "плюс-минус валенок", чтобы в каталоге поиск работал;

- идея использовать точки откуда-нибудь - она верная, но учтите, что и Гугл - не идеал.

- точек на эти снимки вам понадобится не пять, а многие десятки, потому что они не все сняты строго в надир, местность не идеально ровная и так далее - причин этому дофига.

- воровать снимки у Гугла для получения координат вам не обязательно - достаточно в старом интерфейсе google maps найти в левом сайдбаре ссылочку Maps Labs внизу, а далее включить LatLng Marker или LatLng Tooltip (это позволит видеть координаты точки в WGS84 либо в тултипе, либо в виде маркера, который ставится через контекстное меню, только делать это нужно на максимальном приближении, а то они округляются).

- барьер в 2Гб - маловато. Я не пользуюсь QGIS для привязки, использую GlobalMapper 64bit, там такой проблемы у меня нет.

- Геомиксером не пользуюсь, пусть знатоки ответят, но там еще может быть не Меркатор на сфере, а Меркатор на эллипсоиде WGS84. Умеет ли Геомиксер эту проекцию показывать - не знаю.

- В GlobalMapper отдельно задается проекция и датум. Проекция - Меркатор, датум - Google. Можно инициализировать через ввод EPSG:3857.

- Кэш SAS вроде бы в этой же проекции.

- Привязывать и трансформировать растр в новую проекцию лучше один раз в одной программе. Разберетесь как это делать в QGIS - делайте там (я не пользовался, представляю себе процесс только по документации). Разберетесь с GlobalMapper - делайте в нем (в нем есть полезная фича - даже если у вас указана целевая проекция Меркатора, вы можете в диалоге привязки вводить широту и долготу для точек в географических координатах - он это умеет распознавать).
И там, и там это возможно. Три программы для всего этого не нужны, нужна всего одна.

Вопрос 4 не понял и не вижу в нем смысла.

[Shaggy]

ГеоМиксер работает в EPSG 3395

Создал с помощью TilingTools тайлы, эта прога работает в десятки раз быстрее ГлобалМеппера, почему такая разница?
Кстати TilingTools не понимает ни русские пути, ни русские имена файлов и когда сам лежит и в пути до него есть русские буквы, так же не работает!

Опробовал всё на маленьком кусочке. Но на будущее уже есть вопрос по работе с файлами более 2х гигов, ГлобалМеппер их нормально воспринимает но вот после привязки говорит, что сохранять в геотифф будет более более 3х часов это нормально? QGIS сохраняет привязанный файл размером 1,5 гига буквально считанные минуты.

[ericsson]

Ну, это давно известный факт, что GM на операциях с файлами - тормоз. Разработчики никак не выберут приличную быструю библиотеку для этого дела.

[Shaggy]

Решил посмотреть Mapserver ни где не могу найти скомпилированную версию под linux, не поможете?

[rhot]

Решил посмотреть Mapserver ни где не могу найти скомпилированную версию под linux, не поможете?

а сами чё?

[Shaggy]

а сами не умеем, есть возможность поставить на работающий сервер просто...

[gimran]

Покажите инструкцию админу линуксовому.

[ Сообщение с мобильного устройства ]

[Denis Rykov]

MapServer есть в репах, какой у вас дистр?

[ Сообщение с мобильного устройства ]

[Shaggy]

У меня платный хостинг
http://hosting.nic.ru/tariff201.shtml

не знаю , что там крутится, соответственно нет пакетов