Некоторые вопросы по использованию данных SRTM

[AN1440]

1. Если следовать рекомендациям статьи "Импорт данных SRTM с помощью Arcview Spatial Analyst", получаем изображение в проекции "географическая" (градусы). В виде АркВью она выглядит сплющенной и употребить ее в таком виде трудно.
Можно ли перепроектировать?
Попытался вручную перевести данные SRTM-файла в текстовый формат, а затем загрузить его в АркВью и получить точечную тему (по пути можно и пересчитать из WGS84 в Красовского, пересчитать отрицательные значения ну и все такое). Получилось нормально - квадрат 1°х1°, равномерно заполненный 1442401 точками, атрибуты в таблице - широта, долгота и высота. Теперь АркВью может его перепроектировать во что угодно. Только вот SpatialAnalist почему-то наотрез не желает такую тему Interpolate to Grid (пункт в меню недоступен).
-----------------------------------
2. Как-то странно расчитывается освещенность поверхности. На общей картине выхватываются отдельные вершины и обрабатывается их освещенная сторона. Все остальное - и теневая сторона, и горизонтальные поверхности - имеет Value=0. Причем, могут стоять рядом 2 с виду одинаковые вершины, одна из них освещена, а другая - Value=0.
Вот если делать картинку по методике из статьи "Визуализация рельефа в Arcview GIS" (исходя из темы горизонталей), рельеф действительно получается "лестницей", но освещенность расчитывается вполне естественно. Можно поставить точку, проваленную относительно плоскости, и получить изображение конусообразной воронки на теневой картине.
-----------------------------------
3. В статье упоминалось, что 1201х1201 имеет лишний ряд.
Надо полагать, и лишнюю колонку? Ряд, надо полагать, верхний, а колонка - правая?

[Максим Дубинин]

1. Перепроектировать растр можно, не переводя в точки, ибо это - нонсенс и потеря информации из-за последующей интерполяции, в:
-Arcinfo
-ERDAS
-Arcview (нужно спец. расширение - есть на Arcscripts)
-ArcGIS
Для большинства из вышеперечисленных есть описание как именно это делать в упомянутом ПО. Поищите в вопросах/ответах, если не найдете, скажите.

2. Рассчитывается где и на каких данных? Много раз делал hillshade для SRTM - все корректно.

3. Не то чтобы лишний, а дублирующий. Ряд - нижний, колонка - правая. Надо обновит статью, спасибо, что заметили.

[AN1440]

3. Не то чтобы лишний, а дублирующий. Ряд - нижний, колонка - правая. Надо обновит статью, спасибо, что заметили.

- но ведь НИЖНИЙ левый угол - привязочный???

-Arcview (нужно спец. расширение - есть на Arcscripts)

- там нашлось "ProjectRaster" - но, как я понял, оно предназначено для АркВью-9.1, а у меня - 3.30
(там в комментариях:
#Requirement
# Attach this script to a tool in ArcToolbox.
# Developed in ArcGIS ArcView licence Version 9.1
# You May have to change the Toolbox path from D:/Program Files/ArcGIS/ArcToolbox/Toolboxes/Data Management
# Tools.tbx
# to C:/Program Files/ArcGIS/ArcToolbox/Toolboxes/Data anagement Tools.tbx or vice-versa
)

[Максим Дубинин]

но ведь НИЖНИЙ левый угол - привязочный???

Кто сказал? Нет, регистрация изображения осуществляется по верхнему левому (http://gis-lab.info/qa/tfw.html).

а у меня - 3.30

Тоже неплохо.
http://arcscripts.esri.com/details.asp?dbid=11368

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

[AN1440]

Кто сказал? Нет, регистрация изображения осуществляется по верхнему левому (http://gis-lab.info/qa/tfw.html).

-Не, ну я имел в виду, что в названии квадрата фигурирует нижний левый угол... То есть, "у них там" может по нему и привязка? Что в АркВью используется верхний - это известно... Вот я и подумал - а не второй ли сверху ряд надо использовать?...
=============
Вот это расширение - работает. Спасибо! Я понял, что оно было где-то там же, но мимо меня проскочило...
--------------
Странное дело! На перепроектированном гриде тени рисуются просто великолепно! А на исходном - получалась натуральная порнография...

[Максим Дубинин]

в названии квадрата фигурирует нижний левый угол

Да, но как говорится "на заборе тоже написано", то есть название файла может быть любым.

"у них там" может по нему и привязка

Нет, привязка, или просто георегистрация в данном случае осуществляется Arcview, который делает это по этому углу и никакому другому. Так что все ок.

На перепроектированном гриде тени рисуются просто великолепно!

А, стало ясно. На неспроектированном растре надо подбирать отличный от того, что по умолчанию коэффициент вертикального масштабирования (Z), так как в нем единицы X,Y (градусы) не соответствуют Z (метры). После перепроектирования связь восстанавливается. Когда вы делали свой hillshade, то по умолчанию наверняка использовали коэффициент = 1, что так же нонсенс, так как 1градус явно не одно и тоже, что и 1 метр. В случае градусов надо было использовать что-то типа 0,002 и тогда все было бы хорошо.

[AN1440]

... так как в нем единицы X,Y (градусы) не соответствуют Z (метры)...

Мдааа... Получается, малейший подъемчик превращался в вертикальную стену, верхушки гор торчали в космос, далеко за орбиту Шаттла, с которого их снимали... Не хотел бы я быть на том Шаттле...
Просто этот самый коэффициент в процессе обработки нигде не встретился в явном виде - вот и не пришло в голову...

[Максим Дубинин]

Просто этот самый коэффициент в процессе обработки нигде не встретился в явном виде - вот и не пришло в голову...

Да, по умолчанию hillshade вызываемый через интерфейс в Spatial Analyst не дает задавать Z. Надо действовать сторонними готовыми скриптами или пользоваться Avenue.

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

aGrid.HillShade (anAzimuth, anAltitude, zFactor)

Получилась конструктивная беседа, надо бы подновить статьи.

[AN1440]

Приходится поднять тему снова:

Перепроектирование грида с помощью расширения "REPROJECT GRIDS" происходит из одного вида АркВъю в другой, согласно заданным в них проекциям. В исходном виде при выборе проекции "географическая" получается "No projection". Это соответствует исходному состоянию СРТМ-данных? Они, как известно, в "градусах на эллипсоиде ВГС-84".
"No projection" - это оно и есть?

[Максим Дубинин]

Не совсем, "no projection" это т.н. Nil, в зависимости от того как понимает это расширение Nil, оно может назначать ему все что угодно, но в принципе логично было бы назначить Geographic/WGS84, так что скорее всего это одно и тоже.