The drainage of the water seems to be an importante factor of flood in the lowest area. In my prediction I tried to add columns with range of Altitudes. I subset the altitude [40m to 2xxx m] to 5 or 6 levels and included for learning machine. This feature has improved to score but not enougth. Could any one share how can we take the slope in account? Thanks
d2 = data.frame(lon=df$X,lat = df$Y,elevation = df$elevation)
rast = raster::rasterFromXYZ(d2)
raster::crs(rast) = c("+proj=longlat +datum=WGS84 +no_defs +ellps=WGS84 +towgs84=0,0,0")
slope= raster::terrain(rast,opt = 'slope', unit = 'degrees',neighbors = 8)
slope_values = raster::extract(slope,insert your data coordinates in spatialpoints format)
Hope this helps
Many Thanks! It is realy helpful.
How did you extract the latitude and longitude of each rectangle
Columns X(Longitude) and Y(Latitude) are the Coordinates for each rectangle.
Thanks, Holar you are showing a good example on the platform, which will help the beginners and community a lot.
Surely we have more discussion in the end month. Enjoy your Top Position
is there any way to do it in python ?
Check rasterio library
Hello, Can you explain to us how did you subset the altitude to 5 or 6 levels ? . Thank you so much
you need to plot the region using different color depending on Altitude intervals (ranges). Changes the interval following a logic of water movement. Then try to cluster/group surfaces with high risk of Flood with the same color. Finally, you can obtain different color for Lac, Rivers, Mountain, plain, .... It is not perfect but maybe useful. ..
Thank you so much
Help guys!....has anyone managed to calculate slope in Python!?
Compute as arctan. As Holar mentioned before (check code) we have 8 neighbours for 1 square (but squares on edges are exceptions).
Due to wikipedia (https://en.wikipedia.org/wiki/Grade_(slope)) and/or common sense: 1) evaluare elevation difference between neighbours (8 values). 2) evaluate distance between neighbours (8 values too).
I used geopy for task #2. Task #1 evaluated in clear python (I mean numpy and pandas ofc).
train['slope_A'] = (np.arctan(train['elevgrad_A']/(train.apply(lambda row: geodesic(point_1, point_2).m, axis=1))))*180/np.pi
train['slope_A'] - 1 of 8 slopes (pandas)
np.arctan - slope evaluation in rads (np = numpy)
train['elevgrad_A'] - 1 of 8 difference between neighbours
train.apply(lambda row: geodesic(point_1, point_2).m, axis=1 - evaluation for distance between neighbours (geopy)
*180/np.pi - translation to degrees (not sure that this is necessary)
Thanks very much! @Randolf......I will give it a try!