Аэрофото под большим углом к вертикали

[ericsson]

Набрел тут вот на этот снимок
http://www.wwii-photos-maps.com/aerialm ... SG-128.jpg
По крайней мере, на настоящий момент, это - единственное немецкое аэрофото московской части Лосиного острова, доступное в сети. Проблема в том, что сделано оно под углом градусов 60 к вертикали.
Хочется привязать. Точки - не проблема вообще, лесные кварталы с тех пор не менялись, застройка Преображенского, Борогодского и окрестностей Ярославского шоссе хоть и эволюционировала, но многие улицы остались.
Попробовал набить точек в GlobalMapper. При методе трансформации triangulation, конечно, корёжится все за пределами многоугольника, описанного по крайним точкам, а все что на максимальном удалении от камеры (там аж до Садового достает) - все в кашу. С полиномиальной трансформацией все несколько лучше, однако даже внутри плотного множества точек совпадение с современной съемкой страдает сильно.
Вопроса, на самом деле, два:
- есть ли не слишком глубокоспециализированный софт, который бы позволял подобрать параметры съемки и подогнать снимок, манипулируя ими?
- если нет, то каковы были бы ваши действия по привязке этого (я собираюсь отрезать южную часть снимка где-то по линии Измайловского шоссе, Б. Семеновской и дальше к центру, чтобы работать с частями, снятыми под допустимыми углами)?

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

Поиграться параметрами камеры можно в GoogleEath:
- подгружаете снимок как картинку
- играете параметрами камеры, пока снимок с местностью не совпадет
- при совпадении - получите искомые углы и параметры камеры
Единственное - играться можно долго и следует учитывать, что снимок обрезан немного справа и снизу.

[gamm]

Набрел тут вот на этот снимок
Попробовал набить точек в GlobalMapper. При методе трансформации triangulation, конечно, корёжится все за пределами многоугольника, описанного по крайним точкам, а все что на максимальном удалении от камеры (там аж до Садового достает) - все в кашу. С полиномиальной трансформацией все несколько лучше, однако даже внутри плотного множества точек совпадение с современной съемкой страдает сильно.
Вопроса, на самом деле, два:
- есть ли не слишком глубокоспециализированный софт, который бы позволял подобрать параметры съемки и подогнать снимок, манипулируя ими?
- если нет, то каковы были бы ваши действия по привязке этого (я собираюсь отрезать южную часть снимка где-то по линии Измайловского шоссе, Б. Семеновской и дальше к центру, чтобы работать с частями, снятыми под допустимыми углами)?

не тем вы занимаетесь. Снимок скорее всего сделан камерой с центральной проекцией, местность достаточно плоская (гор не видно), поэтому имеем незатейливое проективное преобразование (дробно-линейное), что сводится к решению системы линейных уравнений методом наименьших квадратов. Не заморачиваясь никакими моделями камеры, тем более что "точек море". А вот после проективной коррекции начинайте шаманить. Если точек очень много (сотни), то можно сразу начать с дробно-полиномиального преобразования (это которое RPC).

Если пришлете десяток-другой точек (снимок<->плоские координаты, лучше "конвертом" по всему полю, и обязательно по краям должны быть), то я попробую найти програмульку и посчитать параметры проективного преобразования.Если точек много, то и RPC можно попробовать. Кстати, это все должно быть в GRASS.

Когда я был маленьким, это вообще делали наклоняя снимок под нужным углом на специальном устройстве с крутящимися ручками

[Александр Мурый]

Кстати, это все должно быть в GRASS.

Вот и я думал, как такое вытворить в GRASS, читая мануал к i.ortho.photo. Но что-то про указанные выше преобразования так и не нашёл.. Плюс там надо задавать параметры камеры и т.п.

[gamm]

Кстати, это все должно быть в GRASS.

Вот и я думал, как такое вытворить в GRASS, читая мануал к i.ortho.photo. Но что-то про указанные выше преобразования так и не нашёл.. Плюс там надо задавать параметры камеры и т.п.

значит, остается дробно-линейное или дробно-полиномиальное, кстати для последнего нужны еще и высоты, будет коррекция за высоту.

[ericsson]

Ну собственно..
PixelX,PixelY,X,Y,PointName,Error
это в пулково/гаусс-крюгер

Код: Выделить всё

292.37777461,582.16715950,7417304.470,6194711.439,"Point 1",243.63
1415.15507095,1563.72987772,7419231.484,6190132.251,"Point 2",183.52
2204.07873731,2293.74731090,7421272.693,6185248.901,"Point 3",227.13
714.39800215,2308.27121351,7414383.404,6184232.930,"Point 4",228.01
1530.94168045,520.78240867,7420585.977,6194475.000,"Point 5",93.71
2291.52728850,359.22797518,7422387.082,6194741.196,"Point 8",370.64
1498.48834764,2398.21940037,7418032.131,6183798.943,"Point 7",178.01
2222.15888236,1160.76980052,7421976.584,6192044.975,"Point 9",211.79
361.39194427,1543.16941439,7415310.583,6190255.482,"Point 10",143.30
223.95556390,1949.41555400,7413307.427,6187381.142,"Point 11",217.88
588.91181551,1316.22823672,7416800.461,6191513.335,"Point 12",166.16
497.97974013,309.37622877,7418206.033,6195425.727,"Point 13",274.84

это в пулково/географическая

Код: Выделить всё

292.37777461,582.16715950,37.6790907709,55.8675310028,"Point 1",242.73
1415.15507095,1563.72987772,37.7112237921,55.8267402025,"Point 2",183.63
2204.07873731,2293.74731090,37.7451966219,55.7832278340,"Point 3",226.44
714.39800215,2308.27121351,37.6357433038,55.7729359876,"Point 4",227.40
1530.94168045,520.78240867,37.7315612856,55.8659591904,"Point 5",93.41
2291.52728850,359.22797518,37.7602456057,55.8686424839,"Point 8",369.52
1498.48834764,2398.21940037,37.6939968728,55.7696705021,"Point 7",177.77
2222.15888236,1160.76980052,37.7544650138,55.8443661484,"Point 9",211.31
361.39194427,1543.16941439,37.6486409824,55.8271754940,"Point 10",143.32
223.95556390,1949.41555400,37.6176029544,55.8010120218,"Point 11",217.51
588.91181551,1316.22823672,37.6720223447,55.8387285912,"Point 12",165.96
497.97974013,309.37622877,37.6932719954,55.8740982922,"Point 13",274.31

еще раз - ссылка на картинку
http://www.wwii-photos-maps.com/aerialm ... SG-128.jpg

[gamm]

Ну собственно..
PixelX,PixelY,X,Y,PointName,Error

ушло в личку ...