Commit 513d872a authored by md0276b's avatar md0276b
Browse files

livrable T25

parent 8be2e9ef
---
title: "Space-time distribution of clupeiforms habitats in the French Atlantic and Mediterranean areas."
subtitle: "DEFIPEL project, task 2.5 deliverable"
author: Mathieu Doray, Ifremer / EMH
date: "`r Sys.Date()`"
abstract: "The space-time distribution of clupeiforms habitats in the French Atlantic (Bay of Biscay) and Mediterranean (Gulf of Lion) areas are characterised, based on data collected during routine integrated acoustic surveys. Clupeiforms habitats have been characterised in the Bay of Biscay and Gulf of Lion, based on data collected in springtime during the PELGAS survey (2000-2019), and in summertime during the PELMED survey (2003-2019)."
header-includes:
- \usepackage{float}
- \floatplacement{figure}{H}
output:
bookdown::pdf_document2: default
bookdown::html_document2: default
# pdf_document:
# extra_dependencies: flafter
# html_document:
# df_print: paged
---
```{r setup, include=FALSE}
knitr::opts_chunk$set(echo = FALSE)
#,fig.pos = "!H", out.extra = "")
library(knitr)
library(bookdown)
pref='/media/mathieu/IfremerMData/'
pref='G:/'
path.export.rasters.SPF=paste(pref,'PELGAS/Data/gridMaps/Poissons/gridMapBBspSizeDepth/thrQuantile98/xycor/res0.25/Rasters/',sep='')
rasterMosaics.list=list.files(path.export.rasters.SPF,pattern='rasterMosaic*')
rasterMosaics.list.ane=rasterMosaics.list[substr(rasterMosaics.list,14,21)=='ENGR-ENC']
rasterMosaics.list.ane.caption=gsub('rasterMosaic_','',rasterMosaics.list.ane)
rasterMosaics.list.ane.caption=gsub('.png','',rasterMosaics.list.ane.caption)
rasterMosaics.list.ane.caption=gsub('[(]','',rasterMosaics.list.ane.caption)
rasterMosaics.list.ane.caption=gsub('[]]','cm',rasterMosaics.list.ane.caption)
rasterMosaics.list.ane.caption=gsub('[,]','-',rasterMosaics.list.ane.caption)
rasterMosaics.list.ane.caption=gsub('[_]','-',rasterMosaics.list.ane.caption)
rasterMosaics.list.pil=rasterMosaics.list[substr(rasterMosaics.list,14,21)=='SARD-PIL']
rasterMosaics.list.pil.caption=gsub('rasterMosaic_','',rasterMosaics.list.pil)
rasterMosaics.list.pil.caption=gsub('.png','',rasterMosaics.list.pil.caption)
rasterMosaics.list.pil.caption=gsub('[(]','',rasterMosaics.list.pil.caption)
rasterMosaics.list.pil.caption=gsub('[]]','cm',rasterMosaics.list.pil.caption)
rasterMosaics.list.pil.caption=gsub('[,]','-',rasterMosaics.list.pil.caption)
rasterMosaics.list.pil.caption=gsub('[_]','-',rasterMosaics.list.pil.caption)
rasterMosaics.list.spr=rasterMosaics.list[substr(rasterMosaics.list,14,21)=='SPRA-SPR']
rasterMosaics.list.spr.caption=gsub('rasterMosaic_','',rasterMosaics.list.spr)
rasterMosaics.list.spr.caption=gsub('.png','',rasterMosaics.list.spr.caption)
rasterMosaics.list.spr.caption=gsub('[(]','',rasterMosaics.list.spr.caption)
rasterMosaics.list.spr.caption=gsub('[]]','cm',rasterMosaics.list.spr.caption)
rasterMosaics.list.spr.caption=gsub('[,]','-',rasterMosaics.list.spr.caption)
rasterMosaics.list.spr.caption=gsub('[_]','-',rasterMosaics.list.spr.caption)
#knitr::knit_hooks$set(plot = function(x, options) {
# hook_plot_tex(x, options)
#})
```
# Introduction
- Objectives:
- analyse springtime distribution of clupeiforms in the Bay of Biscay (BoB) in space and time, based on survey (PELGAS) data
- describe main spatial patterns and temporal trends in clupeiforms distribution in springtime from 2000 to 2019
# Material and methods
- Data
- 0.25° gridmaps of 2000-2019 clupeiforms biomass (anchovy, sardine, sprat), per 5cm length class and depth stratum from PELGAS survey
- Space-time MFA (\~ grouped PCA) on series of clupeiforms biomass gridmaps (FactoMineR package)
- Rows: gridmaps cells
- Columns: Clupeiforms biomass per 5cm length class and depth stratum
- Groups: years
- MFA results interpreted in:
- MFA space: individuals (cells) / variables / groups (time)
- Geographical space: MFA maps
# Results
## Deviance explained by MFA axis
```{r fig1, echo=FALSE, out.width="50%",fig.cap='Deviance explained'}
path.export.results.clup=paste(pref,'PELGAS/Resultats/gridMaps/MFA/Poissons/2000-2019/Clupeiforms/',sep='')
load(file=paste(path.export.results.clup,"MFAclupResults.RData",sep=''))
print(paste(dim(big.fssd.clup.MFAs3)[1],'map cells x',dim(big.fssd.clup.MFAs3)[2],'variables',
'grouped into',length(stnnc2.clup),'years'))
NMFA95=dimnames(res.fssd.MFA.clup$eig[round(res.fssd.MFA.clup$eig[,3])>=95,])[[1]]
print(paste('95% deviance explained by',substr(NMFA95[1],6,7),'MFA components'))
barplot(res.fssd.MFA.clup$eig[,3],main="Cumulative % of variance explained",
names.arg=1:nrow(res.fssd.MFA.clup$eig))
```
## MFA individuals (cells) plane
Map cells in MFA1:3 planes are presented in Figure \@ref(fig:fig2):
```{r fig2, echo=FALSE, out.width="70%",message=FALSE,fig.cap="Map cells in MFA1:3 plane"}
library(factoextra)
library(gridExtra)
p1 = fviz_mfa_ind(res.fssd.MFA.clup, col.ind = "contrib",title='')
p2 = fviz_mfa_ind(res.fssd.MFA.clup, col.ind = "contrib",axes = c(2, 3),title='')
grid.arrange(p1, p2, ncol = 2)
```
## Variables contributions to MFA1
```{r fig3, echo=FALSE, out.width="100%",fig.cap = "Variables with abs(correlation) with MFA1 higher or equal to 0.6. Years in top pannel, species as rows, size categories in bottom pannel."}
knitr::include_graphics(paste(path.export.results.clup,'MFA1clupVarCor.png',sep=''))
```
Variables well correlated with MFA1 (Figure \@ref(fig:fig3)):
- MFA1 positive correlation with small/medium anchovy and sardine since 2002 and sprat since 2014
- MFA1 negative correlation with large sardine near sea surface from 2000 to 2002
## Variables contributions to MFA2
```{r fig4, echo=FALSE, out.width="100%",fig.cap = "Variables with abs(correlation) with MFA2 higher or equal to 0.6. Years in top pannel, species as rows, size categories in bottom pannel."}
knitr::include_graphics(paste(path.export.results.clup,'MFA2clupVarCor.png',sep=''))
```
Variables well correlated with MFA2 (Figure \@ref(fig:fig4)):
- MFA2 positive correlation with large anchovy near seabed from 2003 to 2010, negative correlation with small sprat at beginning and end of series.
## Clupeiforms mean spatial patterns: MFA maps
```{r fig5, echo=FALSE, out.width="90%",fig.cap = "MFA1 and 2 loadings mpaps."}
knitr::include_graphics(paste(path.export.results.clup,'rasterLevelPlot_MFAcoordGeoSpace.png',sep=''))
```
Two main spatial gradients in MFA1&2 loadings (Figure \@ref(fig:fig5)):
- MFA1: Gironde mouth vs. NW / offshore areas decreasing gradient
- MFA2: positive MFA2 path South of Belle Ile
## Clupeiforms spatial patterns inter-annual variability
```{r fig6, echo=FALSE, out.width="90%",fig.cap = "Inter-annual inertia of MFA1 and 2 loadings."}
knitr::include_graphics(paste(path.export.results.clup,'rasterLevelPlot_1-2_MFAcoordInertia.png',sep=''))
```
Inter-annual variability (inertia) higher in Gironde mouth (MFA1) and in high MFA2 patch (Figure \@ref(fig:fig6)).
## Global time trends in clupeiforms spatial distribution
```{r fig7, echo=FALSE, out.width="90%",fig.cap = "Clupeiforms MFA1 and 2 time series averaged over whole Bay of Biscay, 2000-2019."}
knitr::include_graphics(paste(path.export.results.clup,'fishMFA12timeSeries.png',sep=''))
```
No significant trend in clupeiforms MFA1&2 time series averaged over the whole Bay of Biscay.(Figure \@ref(fig:fig7)).
## Local clupeiforms time trends
### MFA1 time series in characteristic areas
```{r fig8, echo=FALSE, out.width="90%",fig.cap = "Mean clupeiforms MFA1 and 2 time series in Gironde characteristic area, 2000-2019."}
knitr::include_graphics(paste(path.export.results.clup,'MFA1clupTimeSeriesCharAreas.png',sep=''))
```
Significant positive linear trend in clupeiforms MFA1 averaged over space in Gironde characteristic area for the 2000-2019 period (R2 = 24 %)
### Large sardine mean biomass time series in characteristic areas
```{r fig9, echo=FALSE, out.width="90%",fig.cap = "Large sardine mean biomass in offshore NW characteristic area and Bay of Biscay, 2000-2019."}
knitr::include_graphics(paste(path.export.results.clup,'BigPILbiomassTimeSeriesBoB&MFA1offWAreas.png',sep=''))
```
MFA1 has picked up and localised in offshore W area most of the decreasing trend in big sardine biomass
## Clupeiforms MFA annual anomalies map
```{r fig10, echo=FALSE, out.width="90%",fig.cap = "MFA1 annual anomalies map, 2000-2019."}
knitr::include_graphics(paste(path.export.results.clup,'rasterLevelPlot_MFA1pcoordGeoSpace.png',sep=''))
```
```{r fig11, echo=FALSE, out.width="90%",fig.cap = "MFA2 annual anomalies map, 2000-2019."}
knitr::include_graphics(paste(path.export.results.clup,'rasterLevelPlot_MFA2pcoordGeoSpace.png',sep=''))
```
## Conclusions
- Main Clupeiforms spatial pattern (21% deviance explained): MFA1 Gironde-South vs. NW offshore areas negative gradient.
- More small/medium anchovy and sardine (and sprat since 2014) near seabed in Gironde mouth area since 2003
- Less large sardine near sea surface in NW offshore area since 2002
- Secondary Clupeiforms spatial pattern (13% deviance explained): MFA2 South-offshore vs. Gironde-North onshore negative gradient.
- More small sprat near seabed in Gironde mouth and Northern coastal in from 2001 to 2003 and since 2007
- Clupeiforms time trends:
- No significant trend in clupeiforms MFA1&2 time series averaged over the whole Bay of Biscay
- Significant positive linear trend in clupeiforms MFA1 averaged over space in Gironde characteristic area for the 2000-2019 period: more small/medium Clupeiforms near seabed in Gironde mouth and Northern coastal.
- Overall decrease in Clupeiforms mean length in the BoB caused by 2 processes:
- Decrease of large sardine biomass near sea surface in NW offshore areas since 2002
- Increase of small/medium Clupeiforms biomass near seabed near Gironde mouth
# Annexes: gridmaps
```{r, echo=FALSE, out.width="90%",fig.cap = rasterMosaics.list.ane.caption}
knitr::include_graphics(paste(path.export.rasters.SPF,rasterMosaics.list.ane,sep=''))
#, fig.show = "hold",fig.subcap = c("1", "2", "3", "4")
#for (i in rasterMosaics.list.ane) {
# system(paste0("cd ",path.export.rasters.SPF," && ","cp"," '",i,"'",' file1.png'))
# #cat(paste0("![",i,"](",path.export.rasters.SPF, "file1.png)"),"\n")
# cat(paste0("![test](",path.export.rasters.SPF, "file1.png)"),"\n")
#}
```
```{r, echo=FALSE, out.width="90%",fig.cap = rasterMosaics.list.pil.caption}
knitr::include_graphics(paste(path.export.rasters.SPF,rasterMosaics.list.pil,sep=''))
```
```{r, echo=FALSE, out.width="90%",fig.cap = rasterMosaics.list.spr.caption}
knitr::include_graphics(paste(path.export.rasters.SPF,rasterMosaics.list.spr,sep=''))
```
---
title: "Distribution spatio-temporelle des habitats de clupéiformes dans les eaux françaises en Atlantique et Méditerranée."
subtitle: "Projet DEFIPEL, livrable de la tâche 2.5."
author: Mathieu Doray, Ifremer / EMH
date: "`r Sys.Date()`"
abstract: "The space-time distribution of clupeiforms habitats in the French Atlantic (Bay of Biscay) and Mediterranean (Gulf of Lion) areas are characterised, based on data collected during routine integrated acoustic surveys. Clupeiforms habitats have been characterised in the Bay of Biscay and Gulf of Lion, based on data collected in springtime during the PELGAS survey (2000-2019), and in summertime during the PELMED survey (2003-2019)."
header-includes:
- \usepackage{float}
- \floatplacement{figure}{H}
output:
bookdown::pdf_document2: default
bookdown::html_document2: default
# pdf_document:
# extra_dependencies: flafter
# html_document:
# df_print: paged
---
```{r setup, include=FALSE}
knitr::opts_chunk$set(echo = FALSE)
#,fig.pos = "!H", out.extra = "")
library(knitr)
library(bookdown)
pref='/media/mathieu/IfremerMData/'
pref='G:/'
path.export.rasters.SPF=paste(pref,'PELGAS/Data/gridMaps/Poissons/gridMapBBspSizeDepth/thrQuantile98/xycor/res0.25/Rasters/',sep='')
rasterMosaics.list=list.files(path.export.rasters.SPF,pattern='rasterMosaic*')
rasterMosaics.list.ane=rasterMosaics.list[substr(rasterMosaics.list,14,21)=='ENGR-ENC']
rasterMosaics.list.ane.caption=gsub('rasterMosaic_','',rasterMosaics.list.ane)
rasterMosaics.list.ane.caption=gsub('.png','',rasterMosaics.list.ane.caption)
rasterMosaics.list.ane.caption=gsub('[(]','',rasterMosaics.list.ane.caption)
rasterMosaics.list.ane.caption=gsub('[]]','cm',rasterMosaics.list.ane.caption)
rasterMosaics.list.ane.caption=gsub('[,]','-',rasterMosaics.list.ane.caption)
rasterMosaics.list.ane.caption=gsub('[_]','-',rasterMosaics.list.ane.caption)
rasterMosaics.list.pil=rasterMosaics.list[substr(rasterMosaics.list,14,21)=='SARD-PIL']
rasterMosaics.list.pil.caption=gsub('rasterMosaic_','',rasterMosaics.list.pil)
rasterMosaics.list.pil.caption=gsub('.png','',rasterMosaics.list.pil.caption)
rasterMosaics.list.pil.caption=gsub('[(]','',rasterMosaics.list.pil.caption)
rasterMosaics.list.pil.caption=gsub('[]]','cm',rasterMosaics.list.pil.caption)
rasterMosaics.list.pil.caption=gsub('[,]','-',rasterMosaics.list.pil.caption)
rasterMosaics.list.pil.caption=gsub('[_]','-',rasterMosaics.list.pil.caption)
rasterMosaics.list.spr=rasterMosaics.list[substr(rasterMosaics.list,14,21)=='SPRA-SPR']
rasterMosaics.list.spr.caption=gsub('rasterMosaic_','',rasterMosaics.list.spr)
rasterMosaics.list.spr.caption=gsub('.png','',rasterMosaics.list.spr.caption)
rasterMosaics.list.spr.caption=gsub('[(]','',rasterMosaics.list.spr.caption)
rasterMosaics.list.spr.caption=gsub('[]]','cm',rasterMosaics.list.spr.caption)
rasterMosaics.list.spr.caption=gsub('[,]','-',rasterMosaics.list.spr.caption)
rasterMosaics.list.spr.caption=gsub('[_]','-',rasterMosaics.list.spr.caption)
#knitr::knit_hooks$set(plot = function(x, options) {
# hook_plot_tex(x, options)
#})
```
# Introduction
- Objectives:
- analyse springtime distribution of clupeiforms in the Bay of Biscay (BoB) in space and time, based on survey (PELGAS) data
- describe main spatial patterns and temporal trends in clupeiforms distribution in springtime from 2000 to 2019
# Material and methods
- Data
- 0.25° gridmaps of 2000-2019 clupeiforms biomass (anchovy, sardine, sprat), per 5cm length class and depth stratum from PELGAS survey
- Space-time MFA (\~ grouped PCA) on series of clupeiforms biomass gridmaps (FactoMineR package)
- Rows: gridmaps cells
- Columns: Clupeiforms biomass per 5cm length class and depth stratum
- Groups: years
- MFA results interpreted in:
- MFA space: individuals (cells) / variables / groups (time)
- Geographical space: MFA maps
# Results
## Deviance explained by MFA axis
```{r fig1, echo=FALSE, out.width="50%",fig.cap='Deviance explained'}
path.export.results.clup=paste(pref,'PELGAS/Resultats/gridMaps/MFA/Poissons/2000-2019/Clupeiforms/',sep='')
load(file=paste(path.export.results.clup,"MFAclupResults.RData",sep=''))
print(paste(dim(big.fssd.clup.MFAs3)[1],'map cells x',dim(big.fssd.clup.MFAs3)[2],'variables',
'grouped into',length(stnnc2.clup),'years'))
NMFA95=dimnames(res.fssd.MFA.clup$eig[round(res.fssd.MFA.clup$eig[,3])>=95,])[[1]]
print(paste('95% deviance explained by',substr(NMFA95[1],6,7),'MFA components'))
barplot(res.fssd.MFA.clup$eig[,3],main="Cumulative % of variance explained",
names.arg=1:nrow(res.fssd.MFA.clup$eig))
```
## MFA individuals (cells) plane
Map cells in MFA1:3 planes are presented in Figure \@ref(fig:fig2):
```{r fig2, echo=FALSE, out.width="70%",message=FALSE,fig.cap="Map cells in MFA1:3 plane"}
library(factoextra)
library(gridExtra)
p1 = fviz_mfa_ind(res.fssd.MFA.clup, col.ind = "contrib",title='')
p2 = fviz_mfa_ind(res.fssd.MFA.clup, col.ind = "contrib",axes = c(2, 3),title='')
grid.arrange(p1, p2, ncol = 2)
```
## Variables contributions to MFA1
```{r fig3, echo=FALSE, out.width="100%",fig.cap = "Variables with abs(correlation) with MFA1 higher or equal to 0.6. Years in top pannel, species as rows, size categories in bottom pannel."}
knitr::include_graphics(paste(path.export.results.clup,'MFA1clupVarCor.png',sep=''))
```
Variables well correlated with MFA1 (Figure \@ref(fig:fig3)):
- MFA1 positive correlation with small/medium anchovy and sardine since 2002 and sprat since 2014
- MFA1 negative correlation with large sardine near sea surface from 2000 to 2002
## Variables contributions to MFA2
```{r fig4, echo=FALSE, out.width="100%",fig.cap = "Variables with abs(correlation) with MFA2 higher or equal to 0.6. Years in top pannel, species as rows, size categories in bottom pannel."}
knitr::include_graphics(paste(path.export.results.clup,'MFA2clupVarCor.png',sep=''))
```
Variables well correlated with MFA2 (Figure \@ref(fig:fig4)):
- MFA2 positive correlation with large anchovy near seabed from 2003 to 2010, negative correlation with small sprat at beginning and end of series.
## Clupeiforms mean spatial patterns: MFA maps
```{r fig5, echo=FALSE, out.width="90%",fig.cap = "MFA1 and 2 loadings mpaps."}
knitr::include_graphics(paste(path.export.results.clup,'rasterLevelPlot_MFAcoordGeoSpace.png',sep=''))
```
Two main spatial gradients in MFA1&2 loadings (Figure \@ref(fig:fig5)):
- MFA1: Gironde mouth vs. NW / offshore areas decreasing gradient
- MFA2: positive MFA2 path South of Belle Ile
## Clupeiforms spatial patterns inter-annual variability
```{r fig6, echo=FALSE, out.width="90%",fig.cap = "Inter-annual inertia of MFA1 and 2 loadings."}
knitr::include_graphics(paste(path.export.results.clup,'rasterLevelPlot_1-2_MFAcoordInertia.png',sep=''))
```
Inter-annual variability (inertia) higher in Gironde mouth (MFA1) and in high MFA2 patch (Figure \@ref(fig:fig6)).
## Global time trends in clupeiforms spatial distribution
```{r fig7, echo=FALSE, out.width="90%",fig.cap = "Clupeiforms MFA1 and 2 time series averaged over whole Bay of Biscay, 2000-2019."}
knitr::include_graphics(paste(path.export.results.clup,'fishMFA12timeSeries.png',sep=''))
```
No significant trend in clupeiforms MFA1&2 time series averaged over the whole Bay of Biscay.(Figure \@ref(fig:fig7)).
## Local clupeiforms time trends
### MFA1 time series in characteristic areas
```{r fig8, echo=FALSE, out.width="90%",fig.cap = "Mean clupeiforms MFA1 and 2 time series in Gironde characteristic area, 2000-2019."}
knitr::include_graphics(paste(path.export.results.clup,'MFA1clupTimeSeriesCharAreas.png',sep=''))
```
Significant positive linear trend in clupeiforms MFA1 averaged over space in Gironde characteristic area for the 2000-2019 period (R2 = 24 %)
### Large sardine mean biomass time series in characteristic areas
```{r fig9, echo=FALSE, out.width="90%",fig.cap = "Large sardine mean biomass in offshore NW characteristic area and Bay of Biscay, 2000-2019."}
knitr::include_graphics(paste(path.export.results.clup,'BigPILbiomassTimeSeriesBoB&MFA1offWAreas.png',sep=''))
```
MFA1 has picked up and localised in offshore W area most of the decreasing trend in big sardine biomass
## Clupeiforms MFA annual anomalies map
```{r fig10, echo=FALSE, out.width="90%",fig.cap = "MFA1 annual anomalies map, 2000-2019."}
knitr::include_graphics(paste(path.export.results.clup,'rasterLevelPlot_MFA1pcoordGeoSpace.png',sep=''))
```
```{r fig11, echo=FALSE, out.width="90%",fig.cap = "MFA2 annual anomalies map, 2000-2019."}
knitr::include_graphics(paste(path.export.results.clup,'rasterLevelPlot_MFA2pcoordGeoSpace.png',sep=''))
```
## Conclusions
- Main Clupeiforms spatial pattern (21% deviance explained): MFA1 Gironde-South vs. NW offshore areas negative gradient.
- More small/medium anchovy and sardine (and sprat since 2014) near seabed in Gironde mouth area since 2003
- Less large sardine near sea surface in NW offshore area since 2002
- Secondary Clupeiforms spatial pattern (13% deviance explained): MFA2 South-offshore vs. Gironde-North onshore negative gradient.
- More small sprat near seabed in Gironde mouth and Northern coastal in from 2001 to 2003 and since 2007
- Clupeiforms time trends:
- No significant trend in clupeiforms MFA1&2 time series averaged over the whole Bay of Biscay
- Significant positive linear trend in clupeiforms MFA1 averaged over space in Gironde characteristic area for the 2000-2019 period: more small/medium Clupeiforms near seabed in Gironde mouth and Northern coastal.
- Overall decrease in Clupeiforms mean length in the BoB caused by 2 processes:
- Decrease of large sardine biomass near sea surface in NW offshore areas since 2002
- Increase of small/medium Clupeiforms biomass near seabed near Gironde mouth
# Annexes: gridmaps
```{r, echo=FALSE, out.width="90%",fig.cap = rasterMosaics.list.ane.caption}
knitr::include_graphics(paste(path.export.rasters.SPF,rasterMosaics.list.ane,sep=''))
#, fig.show = "hold",fig.subcap = c("1", "2", "3", "4")
#for (i in rasterMosaics.list.ane) {
# system(paste0("cd ",path.export.rasters.SPF," && ","cp"," '",i,"'",' file1.png'))
# #cat(paste0("![",i,"](",path.export.rasters.SPF, "file1.png)"),"\n")
# cat(paste0("![test](",path.export.rasters.SPF, "file1.png)"),"\n")
#}
```
```{r, echo=FALSE, out.width="90%",fig.cap = rasterMosaics.list.pil.caption}
knitr::include_graphics(paste(path.export.rasters.SPF,rasterMosaics.list.pil,sep=''))
```
```{r, echo=FALSE, out.width="90%",fig.cap = rasterMosaics.list.spr.caption}
knitr::include_graphics(paste(path.export.rasters.SPF,rasterMosaics.list.spr,sep=''))
```
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment