Skip to content

Commit

Permalink
update data
Browse files Browse the repository at this point in the history
  • Loading branch information
gisma committed Sep 21, 2023
1 parent 20bf1d0 commit deae2a6
Show file tree
Hide file tree
Showing 3 changed files with 5 additions and 4 deletions.
Binary file added mc_session/data/all_GW1000A-WIFIDB2F.xlsx
Binary file not shown.
Binary file added mc_session/data/all_GW1000A-WIFIFC29.xlsx
Binary file not shown.
9 changes: 5 additions & 4 deletions mc_session/mc3.qmd
Original file line number Diff line number Diff line change
Expand Up @@ -322,8 +322,8 @@ fnDTM = paste0(wd,"copernicus_DEM.tif")
# Weather Data adapt if you download a new file
# https://www.ecowitt.net/home/index?id=20166)
# https://www.ecowitt.net/home/index?id=149300
fn_dataFC29 = paste0(wd,"all_GW1000A-WIFIFC29(202308270000-202308292319).xlsx")
fn_dataDB2F =paste0(wd,"all_GW1000A-WIFIDB2F(202308270000-202308292219).xlsx")
fn_dataFC29 = paste0(wd,"all_GW1000A-WIFIFC29.xlsx")
fn_dataDB2F =paste0(wd,"all_GW1000A-WIFIDB2F.xlsx")
# station data as derived by the field group
fn_pos_data= paste0(wd,"stations_prelim.shp")
Expand Down Expand Up @@ -454,11 +454,12 @@ k <- krige(A20230829220000 ~ altitude, m, st_as_stars(DTM),
plot(k)
crs(voronoi)=crs
v=sf::st_as_sf(voronoi)
# map it
mapview(raster(interpNN) ,col=rainbow(25)) +
mapview( raster(interpIDW),col=rainbow(25)) +
mapview(k,col=rainbow(25))+
sf::st_as_sf(voronoi)
mapview(k,col=rainbow(25))+ v
```

### Kriging with the relation to the DGM
Expand Down

0 comments on commit deae2a6

Please sign in to comment.