Моделирование лесных пожаров на основе ГИС-данных

[seidju]

Здравствуйте! Я пишу диссертацию на тему моделирования лесных пожаров. В настоящее время исследую возможность моделирования на реальных данных, полученных на основе электронных карт. Американские коллеги поделились интересной программой, которая на основе DEM-карт создает трехмерную модель ландшафта, на которую впоследствии наносятся горючие лесные материалы и моделируется пожар, причем карты должны быть высокого разрешения.
Я с ГИС никогда не работал, ибо занимаюсь высокопроизводительными вычислениями, так вот американцы мне скинули пример карты на основе SRTM по моим координатам, однако разрешение это карты не очень выско, в результате не удается корректно создать модель. Я пытался самостоятельно скачать SRTM и открыть его в программе, однако, как мне сказали с картой предварительно нужно что-то сделать в ArcGIS, а уже потом добавлять в программу. Вот что мне написали:
Nonetheless, I had already processed the SRTM data so you will find that on Box at (I have been assuming the Zunits of this dataset are meters, if that is not the case it needs further processing to convert z units to meters (Use ArcGIS Spatial Analyst/Times tool to perform conversion. If you wish to resample the data, and I am not saying that is the correct thing to do, you can do that using ARcGIS using the resample tool. You can resample to a smaller resolution and try that.
Пробовал использовать утилиту Spatial Analyst, но как то ничего не вышло.
Подскажите пожалуйста, какие предварительные действия необходимо совершить с картой?
А так же как привести данные с ASTER к виду, как прикладываемая карта SRTM?
Буду признателен любой помощи! Спасибо!
P.S. Прикладываю саму карту http://dfiles.ru/files/fpw2riotu

[nadiopt]

ASTER GDEM 2 ни к какому виду приводить не надо, она такая же, как SRTM, только разрешение выше. Ресэмпл, как вам правильно сказали, тут не годится.
http://gis-lab.info/qa/aster-gdem.html

[seidju]

Хорошо, ASTER такая же как и SRTM, но как привести эти данные к виду, что во вложении? Выбрать систему координат как я понимаю необходимо, а так же единицу измерения (метры), как это сделать?

[nadiopt]

мне не посмотреть, что во вложении:( а вы скачали Гдем? она по идее уже в системе координат, как мне помнится, WGS UTM (ну или просто WGS)

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

Американские коллеги поделились интересной программой, которая на основе DEM-карт создает трехмерную модель ландшафта, на которую впоследствии наносятся горючие лесные материалы и моделируется пожар, причем карты должны быть высокого разрешения.

А что за программа, интересно? И какие конкретно у неё требования к моделям рельефа ("высокого разрешения" — это слишком общо).

Скачал приложенный файл "srtm_55_01_UTM.tif" с параметрами:

Спойлер

Driver: GTiff/GeoTIFF
Files: srtm_55_01_UTM.tif
Size is 4505, 7894
Coordinate System is:
PROJCS["WGS_1984_UTM_Zone_46N",
GEOGCS["WGS 84",
DATUM["WGS_1984",
SPHEROID["WGS 84",6378137,298.257223563,
AUTHORITY["EPSG","7030"]],
AUTHORITY["EPSG","6326"]],
PRIMEM["Greenwich",0],
UNIT["degree",0.0174532925199433],
AUTHORITY["EPSG","4326"]],
PROJECTION["Transverse_Mercator"],
PARAMETER["latitude_of_origin",0],
PARAMETER["central_meridian",93],
PARAMETER["scale_factor",0.9996],
PARAMETER["false_easting",500000],
PARAMETER["false_northing",0],
UNIT["metre",1,
AUTHORITY["EPSG","9001"]],
AUTHORITY["EPSG","32646"]]
Origin = (308124.367862528888509,6655205.483634565025568)
Pixel Size = (71.000000000000000,-71.000000000000000)
Metadata:
AREA_OR_POINT=Area
Image Structure Metadata:
COMPRESSION=LZW
INTERLEAVE=BAND
Corner Coordinates:
Upper Left ( 308124.368, 6655205.484) ( 89d33'36.27"E, 59d59'21.35"N)
Lower Left ( 308124.368, 6094731.484) ( 90d 0'10.06"E, 54d57'45.08"N)
Upper Right ( 627979.368, 6655205.484) ( 95d17'44.67"E, 60d 0'50.85"N)
Lower Right ( 627979.368, 6094731.484) ( 94d59'59.84"E, 54d58'58.87"N)
Center ( 468051.868, 6374968.484) ( 92d28' 0.04"E, 57d30'58.12"N)
Band 1 Block=128x128 Type=Int16, ColorInterp=Gray
NoData Value=-32768

Как я понял, нужно просто перепроецировать в какой-либо ГИС (не обязательно в ArcGIS) нужные куски ASTER GDEM в UTM 46-ю зону и всё. А высоты там уже исходно в метрах.

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

Farsite небось какой-нибудь.

[seidju]

Программа называется WFDS, FARSITE для наших условий использовать не получится ибо нет у нас топливных карт, плюс модель WFDS все таки более приближена к реальности. В WFDS задается входной файл сцены, что там будет - зависит от пользователя, собственно мне дали утилиту которая создает этот самый входной файл на основе карт рельефа.
http://www2.bfrl.nist.gov/userpages/wmell/public.html - в самом низу

[seidju]

Вот интересующий меня участок на основе ASTER http://dfiles.ru/files/pjb6jco05

[seidju]

Собственно я бы и хотел узнать как это сделать, более менее работаю с ArcGIS, но обилие инструментов меня приводит в ужас, а google не сильно помогает ввиду отсутсвия знаний в этой области.

Американские коллеги поделились интересной программой, которая на основе DEM-карт создает трехмерную модель ландшафта, на которую впоследствии наносятся горючие лесные материалы и моделируется пожар, причем карты должны быть высокого разрешения.

А что за программа, интересно? И какие конкретно у неё требования к моделям рельефа ("высокого разрешения" — это слишком общо).

Скачал приложенный файл "srtm_55_01_UTM.tif" с параметрами:

Спойлер

Driver: GTiff/GeoTIFF
Files: srtm_55_01_UTM.tif
Size is 4505, 7894
Coordinate System is:
PROJCS["WGS_1984_UTM_Zone_46N",
GEOGCS["WGS 84",
DATUM["WGS_1984",
SPHEROID["WGS 84",6378137,298.257223563,
AUTHORITY["EPSG","7030"]],
AUTHORITY["EPSG","6326"]],
PRIMEM["Greenwich",0],
UNIT["degree",0.0174532925199433],
AUTHORITY["EPSG","4326"]],
PROJECTION["Transverse_Mercator"],
PARAMETER["latitude_of_origin",0],
PARAMETER["central_meridian",93],
PARAMETER["scale_factor",0.9996],
PARAMETER["false_easting",500000],
PARAMETER["false_northing",0],
UNIT["metre",1,
AUTHORITY["EPSG","9001"]],
AUTHORITY["EPSG","32646"]]
Origin = (308124.367862528888509,6655205.483634565025568)
Pixel Size = (71.000000000000000,-71.000000000000000)
Metadata:
AREA_OR_POINT=Area
Image Structure Metadata:
COMPRESSION=LZW
INTERLEAVE=BAND
Corner Coordinates:
Upper Left ( 308124.368, 6655205.484) ( 89d33'36.27"E, 59d59'21.35"N)
Lower Left ( 308124.368, 6094731.484) ( 90d 0'10.06"E, 54d57'45.08"N)
Upper Right ( 627979.368, 6655205.484) ( 95d17'44.67"E, 60d 0'50.85"N)
Lower Right ( 627979.368, 6094731.484) ( 94d59'59.84"E, 54d58'58.87"N)
Center ( 468051.868, 6374968.484) ( 92d28' 0.04"E, 57d30'58.12"N)
Band 1 Block=128x128 Type=Int16, ColorInterp=Gray
NoData Value=-32768

Как я понял, нужно просто перепроецировать в какой-либо ГИС (не обязательно в ArcGIS) нужные куски ASTER GDEM в UTM 46-ю зону и всё. А высоты там уже исходно в метрах.

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

Используйте инструмент "Project Raster".