zsmagyari/ClimTest

1library(shiny)2library(leaflet)3library(raster)4library(ncdf4)5library(ncdf4.helpers)6library(bslib)7 8r_colors <- rgb(t(col2rgb(colors()) / 255))9names(r_colors) <- colors()10fl <- nc_open("tg_ens_mean_0.1deg_reg_v29.0e.nc")11rast <- raster("tg_ens_mean_0.1deg_reg_v29.0e.nc",band=1,varname="tg")12proj4string(rast)=CRS("+init=EPSG:4326")13pal <- colorNumeric(c("#0C2C84", "#41B6C4", "#FFFFCC"), values(rast), na.color = "transparent")14varNames <- nc.get.variable.list(fl)15 16ui <- page_sidebar(17        theme = bs_theme(bootswatch = "minty"),18        sidebar = sidebar(19        selectInput("vars", "Variables", varNames),20        numericInput("band","Day",1,min=1,max=365),  21        p(),22        p(),23        p(),24        p()25      ),26     card(27         full_screen=TRUE,28         card_header="Map",29         leafletOutput("mymap", height=1000)30       )31    )32 33server <- function(input, output, session) {34 35  observeEvent(input$band, {36    rast <- raster("tg_ens_mean_0.1deg_reg_v29.0e.nc",band=input$band,varname="tg")37    proj4string(rast)=CRS("+init=EPSG:4326")38    leafletProxy("mymap") %>% addRasterImage(rast, colors=pal, opacity = 0.5)39  })40  41  output$mymap <- renderLeaflet({42    leaflet() %>%43      addTiles() %>%44      addRasterImage(rast,colors = pal, opacity = 0.8)45  })46}47 48shinyApp(ui, server)