Форумы GIS-Lab.info

Недавно поделились ссылкой на 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 м

Как так?

высота через 1 м;

а DEM где там?

Например, в грануле 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

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

Это рельеф, который использовался для создания RTC

Из 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.

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

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

Форумы GIS-Lab.info

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

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

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

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

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

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

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