Frana M.Javello e Sentinel 1

Qualche settimana fa un conoscente mi ha mandato um video di una corona di frana sul monte Javello legata all'evento piovoso del 6 novembre

Ho provato a vedere se riuscivo a vedere l'estensione della frana dalle immagini di Sentinel 2. Sono riuscito ad inviduare un aspetto interessante ovvero che sul versante era presente una zona disboscata che e' stata tagliata nel 2023 ed i corrispondenza di questa area 

Immagine dell'area disboscata

Prima del taglio forestale

dopo il taglio forestale

Come si puo' vedere dopo l'evento l'impluvio risulta essere interessato da trasporto solido (e' la striscia bianca che scende verso valla a circa meta' dell'immagine leggermente a sinistra)

Nonostante cio' in ottico della frana nemmeno l'ombra. Ho provato, cosa che non avevo mai provato, ad utilizzare Sentinel 1.

Per fare change detection con immagini radar si devono cercare immagini appartenenti alla stessa orbita. il primo tentativo e' fallito perche' il sensore illuminava il versante sbagliato del M.Ferrato e quindi non si mostravano modifiche al versante

Ho preso quindi la coppia 

S1A_IW_GRDH_1SDV_20231031T171457_20231031T171522_051012_062696_0499.SAFE
S1A_IW_GRDH_1SDV_20231112T171456_20231112T171521_051187_062C9B_8EBB.SAFE

Ciascuna immagine e' stata trattata con in SNAP con

Radar/Apply Orbit
Radar/Radiometric/S1 Remove Thermal Noise
Radar/Radiometric/Calibrate (attivando sia Sigma0 che Beta0)
Radar/Speckle Filtering/Single Product Speckle Filtering
Radar/Geometric/Terrain Correction/Range Doppler Terrain Correction
Raster/Subset (per inquadrare solo l'area di interesse)
Raster/Data Conversion/Linear to/fromDb
Radar/Coregistration (per passare da due prodotti ad un prodotto solo)

Per il change detection sono utilizzati diversi algoritmi (sia basati su dati SLC che GRD) come . Ho provato il metodo piu' semplice che corrisponde al fare il rapporto (non la differenza) tra le bande Beta0  prima e dopo l'evento

rapporto beta0

scontornando manualmente le aree

Il metodo non e' particolarmente diagnostico. L'algoritmo indicava una frana anche a monte di Schignano ma persone del luogo mi hanno detto che non si sono verificati movimenti franosi. Diciamo che puo' essere una indicazione ma non altro

in effetti sembra che Sentinel 1 riesca a verificare delle zone di potenziale frana dove Sentinel 2 non riusciva a vedere. Spero di riuscire a trovare il tempo per un sopralluogo di persona 

1) Strategies for landslide detection using open-access synthetic aperture radar backscatter change in Google Earth Engine Alexander L. Handwerger , Shannan Y. Jones , Pukar Amatya  Hannah R. Kerner  , Dalia B. Kirschbaum , and Mong-Han Huang https://doi.org/10.5194/nhess-2021-283

2) Exploring event landslide mapping using Sentinel-1 SAR backscatter products Michele Santangelo, Mauro Cardinali, Francesco Bucci, Federica Fiorucci, Alessandro Cesare Mondini  Geomorphology Volume 397, 15 January 2022, 108021

Un sacco di anni fa

 20 anni fa ero Flip....150.217.73.242

Maratona Firenze 2023

 Anche quest'annno mi sono intrufolato nelle foto della maratona

DS0152 e geofono

Volevo provare l'oscilloscopio DS0152 e visto che ho in prestito un geofono ho guardato il segnale in uscita. Francamente pensavo che il sensore facesse uscire poche decine di millivolt ma come si vede dalla foto si passa da qualche decina di Volts fino a qualche Volt (ho letto anche fino a 9V ma con sollecitazioni meccaniche che non sono ragionevoli per un geofono)

 Anche la lettura del periodo del geofono (15 Hz) risulta coerente con le caratteristiche di fabbrica

MilkV camera

per mandare in auto esecuzione lo script di streaming si edita il file /mnt/system/auto.sh aggiungendo

/mnt/data/install/CviIspTool.sh 64M

Non riesco a vedere lo stream con VLC su Linux ma usando mpv non ci sono problemi

mpv --no-cache --untimed --no-demuxer-thread --video-sync=audio --vd-lavc-threads=1 rtsp://192.168.42.1:8554/stream0

Piana Firenze-Prato-Pistoia Sentinel 1 6-7 novembre 2023

 

Dettagli 

R 3.6 docker personalizzato

 Un esempio di Dockerfile complesso per R 3.6

FROM rstudio/r-base:3.6-bookworm

ARG DEBIAN_FRONTEND=noninteractive

RUN apt-get -y update

RUN apt-get -y install gdal-bin libgdal-dev libpng-dev libudunits2-dev libfontconfig1-dev libmagick++-dev openjdk-17-jdk apache2 mc nano locales && apt-get clean

RUN service apache2 restart

#ENTRYPOINT ["/usr/sbin/apache2", "-k", "start"]

ENV APACHE_RUN_USER www-data

ENV APACHE_RUN_GROUP www-data

ENV APACHE_LOG_DIR /var/log/apache2

EXPOSE 80

CMD ["apachectl", "-D", "FOREGROUND"]

CMD apachectl -D FOREGROUND

# installa client Oracle

RUN wget https://download.oracle.com/otn_software/java/sqldeveloper/sqlcl-23.2.0.178.1027.zip && unzip sqlcl*.zip

ENV PATH=sqlcl/bin:$PATH

#imposta la localizzazione italiana

RUN sed -i -e 's/# it_IT.UTF-8 UTF-8/it_IT.UTF-8 UTF-8/' /etc/locale.gen && locale-gen

ENV LANG it_IT.UTF-8

ENV LANGUAGE it_IT:it

ENV LC_ALL it_IT.UTF-8

# imposta la timezone a Europe/Rome

ENV TZ="Europe/Rome"

# scarica il sorgente della libreria non presente su CRAN

RUN curl -O "https://cran.r-project.org/src/contrib/Archive/randomForest/randomForest_4.6-10.tar.gz"

#setup R configs

RUN echo "r <- getOption('repos'); r['CRAN'] <- 'http://cran.us.r-project.org'; options(repos = r);" > ~/.Rprofile

RUN Rscript -e "install.packages('zoo')"

RUN Rscript -e "install.packages('leaflet')"

RUN Rscript -e "install.packages('sf')"

RUN Rscript -e "install.packages('leafpop')"

RUN Rscript -e "install.packages('forecast')"

RUN Rscript -e "install.packages('stringr')"

RUN Rscript -e "install.packages('lubridate')"

RUN Rscript -e "install.packages('randomForest_4.6-10.tar.gz', repos = NULL, type='source')"

RUN Rscript -e "install.packages('rpart')"

RUN Rscript -e "install.packages('e1071')"

RUN Rscript -e "install.packages('kknn')"

RUN Rscript -e "install.packages('jpeg')"

RUN Rscript -e "install.packages('data.table')"

RUN Rscript -e "install.packages('tidyr')"

RUN Rscript -e "install.packages('magick')"

RUN Rscript -e "install.packages('leaflet.extras')"

RUN Rscript -e "install.packages('leafem')"

RUN Rscript -e "install.packages('raster')"

# copia i file

ADD go.sh /root/go.sh

ADD omnia_apuane_spool.sql /root/omnia_apuane_spool.sql

ADD omnia_apuane.R /root/omnia_apuane.R

ADD omnia_arno.R /root/omnia_arno.R

ADD omnia_laguna.R /root/omnia_laguna.R

#crea folders in html

RUN mkdir /var/www/html/mappe

RUN mkdir /var/www/html/mappe/omnia_apuane

RUN mkdir /var/www/html/Acque

RUN mkdir /var/www/html/Acque/LAGUNA

RUN mkdir /var/www/html/Acque/ARNO

RUN mkdir /var/www/html/Acque/ARNO/bollettini

# per il ROracle

ADD instantclient-basic-linux.x64-21.11.0.0.0dbru.zip instantclient-basic-linux.x64-21.11.0.0.0dbru.zip

ADD instantclient-sdk-linux.x64-21.11.0.0.0dbru.zip instantclient-sdk-linux.x64-21.11.0.0.0dbru.zip

# ATTENZIONE il comando ADD decompatta in automatico

# i file tar.gz. Per questo motivo la riga successiva usa COPY

COPY ROracle_1.3-2_R_x86_64-linux-gnu.tar.gz /root/ROracle_1.3-2_R_x86_64-linux-gnu.tar.gz

RUN unzip instantclient-basic-linux.x64-21.11.0.0.0dbru.zip -d /opt/oracle/

RUN unzip instantclient-sdk-linux.x64-21.11.0.0.0dbru.zip -d /opt/oracle/

RUN apt-get install -y libaio1

ENV PATH="${PATH}:/opt/oracle/instantclient_21_11"

ENV LD_LIBRARY_PATH="$LD_LIBRARY_PATH:/opt/oracle/instantclient_21_11"

RUN ldconfig

# compila ROracle

RUN Rscript -e 'install.packages("/root/ROracle_1.3-2_R_x86_64-linux-gnu.tar.gz",repos=NULL,type="source")'

RUN Rscript -e "install.packages('RCurl')"

RUN Rscript -e "install.packages('stringr')"

#cancella i file temporanei

RUN rm randomForest_4.6-10.tar.gz

RUN rm sqlcl*.zip

RUN rm /root/*.tar.gz

RUN rm /root/*.zip