Skip to content

Commit

Permalink
Update README.md
Browse files Browse the repository at this point in the history
  • Loading branch information
dncpax committed Feb 11, 2021
1 parent 13b626c commit d94029c
Showing 1 changed file with 7 additions and 7 deletions.
14 changes: 7 additions & 7 deletions README.md
Original file line number Diff line number Diff line change
Expand Up @@ -26,15 +26,15 @@ Exemplo: tiles=29SNB, 29SNC, 29SPB, 29SPC (para o alqueva)\

## gdal_processa_pasta3.py
Processa bandas numa pasta sentinel2 em jp2 16 bits para tif 8bits etrs89 usando um script de comandos gdal externo (gdal_sentinel2_rgbi_reproj_msk4.sh).\
Parameters: pasta com ficheiros a processar e onde se guardam os resultados.\
Parameters: pasta com ficheiros a processar e onde se guardam os resultados.

## gdal_sentinel2_rgbi_reproj_msk4.sh
Script para converter as 4 bandas de jp2 16bits para tif 8bits, e projectar para etrs89, e construir 2 mosaicos: rgb e irg.\
Os mosaicos são TILE_rgb_mask.vrt e TILE_irg_mask.vrt.\
Ou seja, este script é chamado pelo gdal_processa_pasta para cada tile sentinel, e produz para cada TILE: 4 jpeg 8bit, 1 vrt rgb, 1 vrt irg, 1 máscara para cada vrt. Calcula pirâmides e estatísticas para os 2 vrt.
Sintaxe: indicar nome da tile, eg: T29SPB_20190813T112121"\
Têm de existir 4 ficheiros jp2 para converter, eg: b02, b03, b04, b08"\
Cria máscara automaticamente usando o comando nearblack do gdal.\
Cria máscara automaticamente usando o comando nearblack do gdal.

## gdal_stretch2.py
Cria vrt's com constrast stretch, a partir de originais em formato vrt.\
Expand All @@ -46,7 +46,7 @@ Calcula para cada banda o min e max a clipar com a fórmula do desvio padrão: M
Usamos 2.8stddev porque se ajusta bem ao sentinel e é simples de calcular.\
Recalcula as stats considerando a mascara, sendo diferente das calculadas pelo gdalinfo.\
Aplica o stretch usando o SCALE do formato VRT, sem alterar as imagens originais.\
E cria 2 vrt rgb+alpha e irg+alpha: TILE_*_viz.vrt\
E cria 2 vrt rgb+alpha e irg+alpha: TILE_*_viz.vrt

## gdal_criartileindex.py
Processa tileindexes para termos um tileindex time-aware para o mapserver, usando um script de comandos ogr (gdal_tileindexes.sh).\
Expand All @@ -55,13 +55,13 @@ Também cria 2 tileindexes na pasta dos jp2 (1 para rgb, 1 para irg). Estes são
Antes de importar apaga registos que já existam para esta pasta/data para evitar duplicados.\
Ou seja, para cada pasta de data temos 2 vrt (rgb e irg), e um shapefile com a listagem destas imagens.
Depois temos um shapefile global onde estão todas as imagens de todas as datas. Este é usado pelo MapServer.\
Parameters: pasta com ficheiros a processar (data, eg 20190804)\
Parameters: pasta com ficheiros a processar (data, eg 20190804)

## gdal_tileindexes.sh
Cria tileindexes para os vrt numa data usando o comando gdaltindex do gdal.\
E actualiza tileindexes globais, usando o comando ogrinfo e sql para copiar os registos.\
No fim, exporta o shapefile global para json para criar uma lista de datas que possa ser usada no viewer html.\
Sintaxe: par1=pasta com ficheiros eg 20190803.\
Sintaxe: par1=pasta com ficheiros eg 20190803.

## sentinelpt.map
Mapfile que publica todos os vrt RGB e IRG em WMS e com suporte temporal (parametro TIME).
Expand All @@ -71,5 +71,5 @@ Inicialmente tentava aplicar um stretch às imagens, mas esse stretch é agora a
# Estrutura em disco
Os scripts assumem uma estrutura em disco onde se colocam os ficheiros.\
dados - pasta mãe com os scripts, e onde ficam os tileindex globais\
|-- 20200501 - pasta com ficheiros para o dia 20200501\
onde são colocados os ficheiros descarregados pelo script e todos os resultados
        |-- 20200501 - pasta criada pelos scripts com ficheiros para o dia 20200501\
        |-- aaaammdd - idem

0 comments on commit d94029c

Please sign in to comment.