Novo (e melhor) MDE gratuito disponível - FABDEM
Вставка
- Опубліковано 3 січ 2025
- Preprint do artigo (DEMIX): arxiv.org/abs/...
CÓDIGO GEE:
var fabdem = ee.ImageCollection("projects/sat-io/open-datasets/FABDEM"),
countries = ee.FeatureCollection("FAO/GAUL/2015/level0"),
imageVisParam = {"opacity":1,"bands":["b1"],"max":1550.551025390625,"gamma":1},
geometry =
/* color: #98ff00 */
/* displayProperties: [
{
"type": "rectangle"
}
] */
ee.Geometry.Polygon(
[[[-46.110284265802356, -23.071890035886575],
[-46.110284265802356, -23.38864002817647],
[-45.621392664239856, -23.38864002817647],
[-45.621392664239856, -23.071890035886575]]], null, false);
print('FABDEM Collection size :',fabdem.size())
//Explanation on setting default Projection here / 1494038930643042309
var elev = fabdem.mosaic().setDefaultProjection('EPSG:3857',null,30)
//you can also use this incase you don't want to specify CRS
//var elev = fabdem.mosaic().setDefaultProjection(glo30.first().projection())
// Add the elevation to the map. Play with the visualization tools
// to get a better visualization.
Map.addLayer(elev, {}, 'elev', false);
// Use the terrain algorithms to compute a hillshade with 8-bit values.
var shade = ee.Terrain.hillshade(elev);
Map.addLayer(shade, {}, 'hillshade', false);
// Create an "ocean" variable to be used for cartographic purposes
var ocean = elev.lte(0);
Map.addLayer(ocean.mask(ocean), {palette:'000022'}, 'ocean', false);
// Create a custom elevation palette from hex strings.
var elevationPalette = ['006600', '002200', 'fff700', 'ab7634', 'c4d0ff', 'ffffff'];
// Use these visualization parameters, customized by location.
var visParams = {min: 1, max: 3000, palette: elevationPalette};
// Create a mosaic of the ocean and the elevation data
var visualized = ee.ImageCollection([
// Mask the elevation to get only land
elev.mask(ocean.not()).visualize(visParams),
// Use the ocean mask directly to display ocean.
ocean.mask(ocean).visualize({palette:'000022'})
]).mosaic();
// Note that the visualization image doesn't require visualization parameters.
Map.addLayer(visualized.clip(countries), {}, 'elev palette');
// Export the image, specifying scale and region.
// alterar apenas o description
Export.image.toDrive({
image: elev,
description: 'fabdem_sao_jose_campos',
scale: 30,
region: geometry
});
Boa tarde professora Ligia Batista! Tudo bem? Atualmente o MDT FABDEM continua com a primeira posição de melhor MDT?
Olá! Segundo artigo publicado, fica entre o Copernicus DEM e o FABDEM, dependendo da aplicação. O Fabdem na verdade é o copernicus dem com correções. Mais detalhes no vídeo ua-cam.com/video/-vrZzolC_Xs/v-deo.html e no artigo de um grupo internacional dedicado a este tipo de estudo: Bielski, C.; López-Vázquez, C.; Grohmann, C.H.; Guth. P.L.; Hawker, L.; Gesch, D.; Trevisani, S.; Herrera-Cruz, V.; Riazanoff, S.; Corseaux, A.; Reuter, H.; Strobl, P., 2024. Novel approach for ranking DEMs: Copernicus DEM improves one arc second open global topography. IEEE Transactions on Geoscience & Remote Sensing. doi.org/10.1109/TGRS.2024.336...
SENSACIONAL. Muito bacana o vídeo.
Muito obrigada pela audiência!
O artigo do professor Carlos Grohmann, ele fez um vídeo ano passado no canal dele sobre o artigo
Sim. Eu assisti...
Muito bacana conhecer esse DEM. No meu bacharelado de Geografia nos anos 2010 eu trabalhei muito com modelos de elevação, principalmente o SRTM, ASTER e mais pro final o TopoDATA do INPE que era um refino do SRTM de 90 pra 30 m de resolução também.
Estou bem defasado de lá pra cá! hahahaha agora só drone!
Enfim, fiquei curioso, você sabe dizer se já foi feito algum comparativo entre este FABDEM e o TopoDATA (além da vantagem de ser um dado mais atual no caso do FABDEM)???
Obrigado!!!
O artigo que compara vários modelos globais está em ieeexplore.ieee.org/document/10440392
mas o topodata não está entre os modelos comparados, só o srtm mesmo...
Obrigado professora pelo vídeo e pelo script do GEE. Consegui usar e baixar o modelo para minha área de interesse, mas o arquivo saiu deslocado espacialmente, não está geocodificado para a minha área de estudo tem algumas centenas de metros de diferença. Alguém mais teve esse problema?
Olá Gustavo. Me lembro de outra pessoa falando da mesma situação, mas quando faço o teste fica tudo ok. O GEE exporta como sistema padrão de referência o EPSG 3857. A camada está configurada assim?
Muito bom!
Obrigada!
Muito obrigado professora! excelente vídeo.
no código que você disponibilizou tive que ajustar essa parte, pois o polígono ficou fixo nessas coordenadas ( mesmo após eu excluir e inserir um novo ), eu não sei muito de programação então apenas ajustei as cordeadas pro meu polígono de interesse e deu certo.
geometry =
/* color: #98ff00 */
/* displayProperties: [
{
"type": "rectangle"
}
] */
ee.Geometry.Polygon(
[[[-46.110284265802356, -23.071890035886575],
[-46.110284265802356, -23.38864002817647],
[-45.621392664239856, -23.38864002817647],
[-45.621392664239856, -23.071890035886575]]], null, false);
print('FABDEM Collection size :',fabdem.size())
Olá, obrigada pelo feedback. Que bom que conseguiu ajustar!