Data collection
Drone survey over the study area; orthomosaic and DSM generated in Agisoft Metashape and Spatix.
A drone survey of an industrial and residential block in Kamrup Metro District, Assam, turned into a rooftop-level solar suitability assessment: solar radiation, slope, aspect and shadow are standardised, weighted and overlaid to rank every roof for photovoltaic installation.
Weights were assigned from expert judgement and published literature on GIS-based rooftop PV suitability analysis. Solar radiation carries the most weight because it acts directly on generation; slope, shading and aspect follow.
Hillshade was computed on the DSM at five times of day to capture how surrounding buildings and roof structures throw shade across the rooftops. Each slice was reclassified so that sunlit roof area scores higher.
| Time of day | Sun azimuth | Sun altitude |
|---|---|---|
| 08:00 | 92.7° | 40.2° |
| 10:00 | 115.0° | 66.4° |
| 12:00 | 214.4° | 76.3° |
| 14:00 | 259.5° | 52.7° |
| 16:00 | 274.2° | 25.9° |
The weighted overlay index is classified from the most to the least suitable roof. Higher scores mean more radiation, less shadow, a flatter surface and a better orientation.
Note: the source document describes the scheme both as these five named classes (Step 9) and as a four-class dissolve (1–4) when rooftop areas are summarised per class; the map atlas below follows the five named classes. Suitability values on the underlying rasters are shown on a 1–4 scale in the overlay section and a 1–5 scale in the reclassification section.
Nine steps chain the drone capture, terrain derivatives and overlay into one decision-ready surface.
Drone survey over the study area; orthomosaic and DSM generated in Agisoft Metashape and Spatix.
DSM and orthomosaic imported to QGIS / ArcGIS Pro; building rooftop polygons digitised from the orthomosaic as analysis units.
Slope and aspect derived from the DSM — roof inclination and orientation.
Area Solar Radiation tool (ArcGIS Pro) on the DSM, default atmospheric parameters, annual totals.
Hillshade from the DSM at five times of day; reclassified so sunlit roof area scores higher.
Polygon area computed with the Field Calculator in QGIS.
Mean solar radiation sampled per rooftop polygon from the radiation raster.
Each criterion standardised into suitability classes; combined by weighted overlay.
Suitability index classed into Suitable / Moderately Suitable / Less Suitable / Unsuitable.
Every criterion layer and the final overlay, as produced in the analysis. Each map is shown with its legend — what every colour on the map stands for — alongside the legend as exported from QGIS in the source document. The suitability raster was finally polygonised and dissolved by class so each rooftop carries a measurable area.
The source legend lists North twice; the last row carries the darkest red.
Reclassified hillshade: sunlit roof scores 1, shaded roof scores 0.
Counts how many of the five shadow slices each roof cell is sunlit in.
Note the four classes here, against the five named classes in the classification section.
The source document carries no legend for this vector layer. Read it with Fig 12, where yellow is the Suitable class: yellow shows the higher-suitability roof area and dark the lower, dissolved per rooftop polygon.
For the sample building used as the worked check, the polygonised suitability classes were measured with the QGIS Field Calculator. The four class areas sum to the building's classified roof area — the same summary applied across the whole study area to estimate installable rooftop space.
| Class | Area |
|---|---|
| Suitable | 147.28 m² |
| Moderately Suitable | 68.81 m² |
| Less Suitable | 28.51 m² |
| Unsuitable | 62.64 m² |
| Total sample roof | 307.24 m² |
The overlay separates rooftops into distinct suitability bands rather than a single blanket score, so the large flat industrial roofs and the smaller residential roofs read differently on the final map.
The suitability map was validated by visual comparison against the high-resolution orthomosaic. Randomly selected roofs were checked to confirm the assigned class matched the roof's radiation, slope, aspect and available area. Consistency between the GIS classes and the observed roof conditions supported the reliability of the assessment.
Because the raster was polygonised and dissolved by class, the areas can be summed per class. That turns the map into a quantified estimate of the rooftop area available for PV installation across the study area — the step that makes the output usable for planning.
The framework is repeatable and low-cost: a drone flight plus open GIS tooling is enough to rank roofs and prioritise where panels should go first.
“The results identified rooftops with varying levels of suitability, enabling the prioritization of buildings with high solar potential, favourable roof characteristics and adequate installation space.”