Precision snow removal analysis
A method for precision snow removal analysis provides a melt rate in inches per day for roadways in a geometric region is provided. Elevation data is used to calculate upward looking hemispherical viewsheds at each roadway coordinate to determine the total solar radiation. The snowmelt capacity and slope are determined at each coordinate and the data is displayed on a geographical map of the roadway. The processing speed of the computer is improved since the snowmelt capacity is only determined for a roadway subset of the region coordinates.
1. A method for precision snow removal analysis comprising the steps of:
a) acquiring Digital Elevation Model (DEM) coordinate data of a geographical area in an electronic file format;
b) acquiring a roadway surface subset of Geographic Information System (GIS) coordinate data in the geographical area in an electronic file format;
c) calculating with a computer the upward looking hemispherical viewsheds for each coordinate in the DEM that is geographically within the roadway surface subset of the GIS coordinate data;
d) calculating with the computer the total solar radiation for each coordinate in the DEM that coincides geographically within the roadway surface subset of the GIS coordinate data;
e) calculating with the computer the snowmelt capacity in inches per day, using the total solar radiation, for each coordinate in the DEM that resides within the roadway surface subset of the GIS coordinate data; and
f) displaying the snowmelt capacity in inches per day for each coordinate in the DEM that resides within the roadway surface subset of the GIS coordinate data on a colorized map with a computer display device.
2. The method of claim 1 wherein the acquiring step a) utilizes a Light Detection and Ranging (LiDAR) method to acquire the DEM coordinate data.
3. The method of claim 1 wherein the calculating step c) comprises a three step method of calculating a viewshed, a sunmap, and a skymap for each coordinate in the DEM that resides within the roadway surface subset of the GIS coordinate data.
4. The method of claim 3 wherein each skymap is defined by 16 zenith divisions and 16 azimuth divisions.
5. The method of claim 1 wherein the calculating step d) includes calculating a direct solar radiation and a diffuse solar radiation.
6. The method of claim 1 wherein the displaying step f) includes assigning a color to each coordinate in the DEM that resides within the roadway surface subset of the GIS coordinate data to visually indicate the snowmelt capacity in inches per day.
7. The method of claim 6 wherein red and yellow colors indicate higher calculated snow melt rates and green and blue colors indicate a lower or zero calculated snow melt rates.
8. The method of claim 1 wherein the calculating step e) also includes calculating a slope for each coordinate in the DEM that resides within the roadway surface subset of the GIS coordinate data and using the slope and snowmelt calculations for classifying the roadway vulnerability.
9. A method for creating a snowmelt map for roadways comprising the steps of:
a) acquiring a first coordinate data set having a horizontal accuracy of 15 cm to 20 cm and an elevation accuracy of 0.3 meter to 1.0 meter;
b) acquiring a second coordinate data set of a roadway surface that resides geographically within the first coordinate data set;
c) calculating with a computer the upward looking hemispherical viewsheds for each roadway coordinate in the second coordinate data set using the corresponding elevation data in the first coordinate data set;
d) calculating with the computer the total solar radiation for each roadway coordinate in the second coordinate data;
e) calculating with the computer the snowmelt capacity in inches per day for each roadway coordinate in the second coordinate data; and
f) displaying a visual representation of the snowmelt capacity in inches per day for each roadway coordinate in the second coordinate data with a computer display device.
10. The method of claim 9 wherein the acquiring step a) utilizes a Light Detection and Ranging (LiDAR) method to acquire the first coordinate data set.
11. The method of claim 9 wherein the calculating upward looking hemispherical viewsheds step c) comprises a three step method of calculating a viewshed, a sunmap, and a skymap for each coordinate in the second coordinate data set.
12. The method of claim 11 wherein each skymap is defined by 16 zenith divisions and 16 azimuth divisions.
13. The method of claim 9 wherein the calculating step d) includes calculating a direct solar radiation and a diffuse solar radiation.
14. The method of claim 9 wherein the displaying step f) includes assigning a color to each coordinate in the second coordinate data set to visually indicate the snowmelt capacity in inches per day.
15. The method of claim 14 wherein red and yellow colors indicate higher calculated snow melt rates and green and blue colors indicate a lower or zero calculated snow melt rates.
16. The method of claim 9 wherein the calculating step d) also includes calculating a slope for each roadway coordinate in the second coordinate data and using the slope and snowmelt calculations for classifying the roadway vulnerability.
17. A method for creating a snowmelt dataset for roadways where a snowmelt capacity in inches per day is calculated for only a subset of a Digital Elevation Model (DEM) coordinate dataset and where the subset contains only the coordinates that are associated with a roadway in 1-meter resolution.