ALOS PALSAR (12,5 m DEM) Global Radar Imagery, 2006-2011

Получение, новые источники, открытость
Ответить
dime1
Гуру
Сообщения: 938
Зарегистрирован: 29 май 2011, 19:41
Репутация: 93
Откуда: Киев

ALOS PALSAR (12,5 m DEM) Global Radar Imagery, 2006-2011

Сообщение dime1 » 10 ноя 2017, 15:09

Недавно поделились ссылкой на https://www.asf.alaska.edu/sar-data/palsar/ , а скачать геоданные можно с https://vertex.daac.asf.alaska.edu/# (интерфейс кривой, потому приблизительно обвести область интереса, а потом уже можно перемещать, видоизменять многоугольник).

Спутник - Advanced Land Observing Satellite-1 (ALOS), also known as DAICHI;
Съёмочная аппаратура - Phased Array type L-band Synthetic Aperture Radar (PALSAR);
Охват - 87.8 degrees north latitude and 75.9 degrees south latitude;
Дата съёмки - 2006-2008;
Разрешение - 12,5 м, высота дана дискретностью в 1 м;
Уровни обработки - L1.0 (Unprocessed/raw), L1.5 (Geo-referenced amplitude image), RTC (Radiometrically and terrain corrected geo-coded GeoTIFF image);
Артефакты в виде пропусков и дыр - не было замечено.

Скачал геоданные HiRes RTC и сравнивал их с SRTMGL1 (1 arc sec) и у меня следующие вопросы:
1) SAR-L должен же просвечивать кроны деревьев? Но в конечном продукте чётко прослеживается, что это не так и лесные массивы различимы.
2) Где-то читал, что это просто upscaling геоданных SRTM. Или там какая-то хитрая методика data fusion, где за базис берутся геоданные SRTM, но проводится отдельная съёмка для их уточнения и улучшения?
3) У меня при сопоставлении с SRTMGL1 разница абсолютных высот по профилям всегда +20 м :| Как так?
Последний раз редактировалось dime1 10 ноя 2017, 16:10, всего редактировалось 1 раз.

alexandr cherepanov
Гуру
Сообщения: 530
Зарегистрирован: 30 ноя 2006, 13:31
Статьи: 3
Проекты: 1
Репутация: 111
Откуда: Moscow

Re: ALOS PALSAR (12,5 m DEM) Global Radar Imagery, 2006-2011

Сообщение alexandr cherepanov » 10 ноя 2017, 16:00

высота через 1 м;
а DEM где там?

Александр Мурый
Гуру
Сообщения: 5150
Зарегистрирован: 26 сен 2009, 16:26
Статьи: 3
Проекты: 5/1
Репутация: 776
Ваше звание: званий не имею
Откуда: Москва

Re: ALOS PALSAR (12,5 m DEM) Global Radar Imagery, 2006-2011

Сообщение Александр Мурый » 10 ноя 2017, 16:37

Например, в грануле ALPSRP073190860 вот растр с рельефом:

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

$ gdalinfo AP_07319_FBD_F0860_RT1.dem.tif
Driver: GTiff/GeoTIFF
Files: AP_07319_FBD_F0860_RT1.dem.tif
Size is 6444, 5640
Coordinate System is:
PROJCS["WGS 84 / UTM zone 37N",
    GEOGCS["WGS 84",
        DATUM["WGS_1984",
            SPHEROID["WGS 84",6378137,298.257223563,
                AUTHORITY["EPSG","7030"]],
            AUTHORITY["EPSG","6326"]],
        PRIMEM["Greenwich",0,
            AUTHORITY["EPSG","8901"]],
        UNIT["degree",0.0174532925199433,
            AUTHORITY["EPSG","9122"]],
        AUTHORITY["EPSG","4326"]],
    PROJECTION["Transverse_Mercator"],
    PARAMETER["latitude_of_origin",0],
    PARAMETER["central_meridian",39],
    PARAMETER["scale_factor",0.9996],
    PARAMETER["false_easting",500000],
    PARAMETER["false_northing",0],
    UNIT["metre",1,
        AUTHORITY["EPSG","9001"]],
    AXIS["Easting",EAST],
    AXIS["Northing",NORTH],
    AUTHORITY["EPSG","32637"]]
Origin = (554616.062500000000000,4859517.750000000000000)
Pixel Size = (12.500000000000000,-12.500000000000000)
Metadata:
  AREA_OR_POINT=Point
  TIFFTAG_DATETIME=2015:08:21 03:49:20
  TIFFTAG_SOFTWARE=Created with GAMMA Software www.gamma-rs.ch
Image Structure Metadata:
  INTERLEAVE=BAND
Corner Coordinates:
Upper Left  (  554616.062, 4859517.750) ( 39d40'47.73"E, 43d53'12.25"N)
Lower Left  (  554616.062, 4789017.750) ( 39d40'22.16"E, 43d15' 7.15"N)
Upper Right (  635166.062, 4859517.750) ( 40d40'56.67"E, 43d52'34.95"N)
Lower Right (  635166.062, 4789017.750) ( 40d39'53.42"E, 43d14'30.65"N)
Center      (  594891.062, 4824267.750) ( 40d10'30.05"E, 43d33'55.20"N)
Band 1 Block=6444x1 Type=Int16, ColorInterp=Gray
  NoData Value=0
Редактор материалов, модератор форума

dime1
Гуру
Сообщения: 938
Зарегистрирован: 29 май 2011, 19:41
Репутация: 93
Откуда: Киев

Re: ALOS PALSAR (12,5 m DEM) Global Radar Imagery, 2006-2011

Сообщение dime1 » 10 ноя 2017, 18:45

Да, там в архиве если по размерам файлов посмотреть, то HH и VV -поляризация отдельно, а ".dem.tif" - это и есть RTC-файл ЦМВ.

alexandr cherepanov
Гуру
Сообщения: 530
Зарегистрирован: 30 ноя 2006, 13:31
Статьи: 3
Проекты: 1
Репутация: 111
Откуда: Moscow

Re: ALOS PALSAR (12,5 m DEM) Global Radar Imagery, 2006-2011

Сообщение alexandr cherepanov » 10 ноя 2017, 18:48

Это рельеф, который использовался для создания RTC :D
Из product guide https://media.asf.alaska.edu/uploads/RT ... e_v1.2.pdf
.dem.tif - digital elevation model used for terrain correction in GeoTIFF format
60 N to 57 S - SRTM GL1 EGM96 Up samplingproduct 12.5 m
The DEMs were pre-processed by ASF to a consistent raster format (GeoTIFF) from the original source formats: height (*.hgt), ESRI ArcGrid (*.adf), etc. Many of the NASA-provided DEMs were provided as orthometric heights with EGM96 vertical datum. These were converted by ASF to ellipsoid heights using the ASF MapReady tool named geoid_adjust.
The SRTM GL1 data showed a number of artifacts that needed to be corrected. In order to detect the anomalies, the slope to the neighboring pixel was calculated in both image direction. Any absolute local slope greater than 400 m was marked as anomaly. The spikes and holes forming the anomalies in the DEM were removed by adaptive triangulation using the Quick Terrain Modeler software.

dime1
Гуру
Сообщения: 938
Зарегистрирован: 29 май 2011, 19:41
Репутация: 93
Откуда: Киев

Re: ALOS PALSAR (12,5 m DEM) Global Radar Imagery, 2006-2011

Сообщение dime1 » 13 ноя 2017, 13:01

В итоге, какая может быть польза этих продуктов и варианты их использования?

Если .dem.tif - это апсэмплинг-версия до 12,5 м из SRTMGL1 с обработкой склонов, аномалий, еще и приведённая к высоте эллипсоида с помощью "geoid_adjust" ASF MapReady tool. И откорректированный продукт поляризации HH.tif, а также .inc_map.tif маска layover/shadow.
Получается, будет интересна только тем, кто работает с активными ДДЗ. А рельеф этот, если привести обратно к EGM96, то не факт, что качество ЦМВ будет лучше, чем исходный SRTMGL1.
И чего они этим SAR-L PALSAR не уточнили высоту рельефа под лесными массивами...

Ответить

Вернуться в «Данные»

Кто сейчас на конференции

Сейчас этот форум просматривают: нет зарегистрированных пользователей и 2 гостя