{"id":10419,"date":"2023-05-18T10:30:00","date_gmt":"2023-05-18T08:30:00","guid":{"rendered":"https:\/\/www.sigterritoires.fr\/?p=10419"},"modified":"2023-05-16T13:25:04","modified_gmt":"2023-05-16T11:25:04","slug":"qgis-plugin-green-view-index-comparison-with-aerial-images","status":"publish","type":"post","link":"https:\/\/www.sigterritoires.fr\/index.php\/en\/qgis-plugin-green-view-index-comparison-with-aerial-images\/","title":{"rendered":"QGis plugin Green View Index: comparison with aerial images"},"content":{"rendered":"\n

The QGIS Green View Index plugin is a powerful tool to measure the vegetation cover of a given area. It is interesting to see what it brings compared to the classical processing based on vertical aerial photos.<\/p>\n\n\n\n

We publish a series of three articles:<\/p>\n\n\n\n

1-Concepts of the Green View index and installation of the plugin<\/a><\/p>\n\n\n\n

2-Tutorial of the Green View Index plugin.<\/a><\/p>\n\n\n\n

3-Comparison of the results of the plugin and the aerial photo processing (this article)<\/p>\n\n\n\n\n\n\n\n

We are going to compare the green index calculated on two aerial photos of the city center of Brest with the processing of the Green View Index plugin of QGis.<\/p>\n\n\n\n

For this we downloaded two aerial images from the Brest M\u00e9tropole website.<\/p>\n\n\n\n

Each of them covers an area of about 1 km\u00b2.<\/p>\n\n\n\n

\"\"<\/a><\/figure>\n\n\n\n

Generating Sample Points<\/h1>\n\n\n\n

For the calculation of the Green View Index, it is necessary to define an area of interest (AOI). We will therefore create a new layer of polygon type corresponding to each aerial image.<\/p>\n\n\n\n

\"\"<\/a><\/figure>\n\n\n\n

Use the menu Layer-> Create a layer -> New Shapefile<\/a> layer<\/em><\/strong><\/p>\n\n\n\n

\"\"<\/a><\/figure>\n\n\n\n

For the second photo we repeat the operation.<\/p>\n\n\n\n

\"\"<\/a><\/figure>\n\n\n\n

So we have our two AOIs. We can run the first script of the Green View Index plugin, Generate Sample Points, to create a layer of points randomly distributed on each of our AOIs.<\/p>\n\n\n\n

\"\"<\/a><\/figure>\n\n\n\n

1-Using AOI only<\/h2>\n\n\n\n

We use the default settings, i.e. creating 100 points distributed over the AOI. We did not use a road network.<\/p>\n\n\n\n

\"\"<\/a><\/figure>\n\n\n\n

We run the script on each of the AOIs and we obtain two layers of type point: pointsSouth and pointsNorth.<\/p>\n\n\n\n

\"\"<\/a><\/figure>\n\n\n\n

We can then execute the second script of the plugin,Download Google Street View Images.<\/p>\n\n\n\n

\"\"<\/a><\/figure>\n\n\n\n

We have kept the default settings here too: the six directions around the viewpoint and the three view tilts.<\/p>\n\n\n\n

\"\"<\/a><\/figure>\n\n\n\n

In the Log tab you can see the number and the points for which there was no image available.<\/p>\n\n\n\n

\"\"<\/a><\/figure>\n\n\n\n

For the North zone, 5 points have no images. For the South zone 19 points do not give any result.<\/p>\n\n\n\n

While browsing the image directories, we noticed a rather peculiar phenomenon. Indeed, some viewpoints located on buildings return images\u2026 of the interior of the buildings!<\/p>\n\n\n\n

\"\"<\/a><\/figure>\n\n\n\n

Another \u00ab\u00a0anomaly\u00a0\u00bb that we noticed is the presence of shots in the middle of the night.<\/p>\n\n\n\n

\"\"<\/a><\/figure>\n\n\n\n

It is obvious that these images must be excluded from the calculation. In summary we obtained<\/p>\n\n\n\n

NORTH Zone: Out of 100 viewpoints, 5 points without images, 4 points inside buildings and 1 night viewpoint, that is 10% of non-usable points.<\/p>\n\n\n\n

South Zone: Out of 100 viewpoints, 19 points without images and 2 points inside buildings, that is 21% of non-usable points.<\/p>\n\n\n\n

2- Using AOI and a road network.<\/h2>\n\n\n\n

We therefore downloaded the Brest road network to test the joint use of the AOIs and the network.<\/p>\n\n\n\n

\"\"<\/a><\/figure>\n\n\n\n

The following image shows the position of the 100 points resulting from the use of AOIs and the road network (\u22c7) as well as the position of the 100 points resulting from the use of only the AOI (\u2206).<\/p>\n\n\n\n

\"\"<\/a><\/figure>\n\n\n\n

The same checks of the images downloaded from Google Street View were carried out and obtained:<\/p>\n\n\n\n

NORTH Zone: Out of 100 view points, 3 points without images, no points inside buildings and no night view points, that is 3% of points not usable.<\/p>\n\n\n\n

South Zone: Out of 100 viewpoints, no points without images and 1 point inside buildings, that is to say 1% of non-usable points.<\/p>\n\n\n\n

It seems clear that, in dense urban areas, the use of the road network is highly necessary. The difference in results will be seen later.<\/p>\n\n\n\n

Green View Index calculation<\/h1>\n\n\n\n

The Green View Index is calculated for each area and each type of point generation<\/p>\n\n\n\n

\"\"<\/a><\/figure>\n\n\n\n

The results of the average GVI obtained are as follows:<\/p>\n\n\n\n

North Zone<\/p>\n\n\n\n

All AOI points\u2026\u2026\u2026\u2026\u2026\u2026\u2026\u2026\u2026\u2026\u2026\u2026 0.061<\/p>\n\n\n\n

AOI points without anomalies\u2026\u2026.. 0.064<\/p>\n\n\n\n

All AOI and roadway points\u2026\u2026……. 0.051<\/p>\n\n\n\n

South Zone<\/p>\n\n\n\n

All AOI points\u2026\u2026\u2026\u2026\u2026\u2026\u2026\u2026\u2026\u2026\u2026\u2026 0.034<\/p>\n\n\n\n

AOI points without anomalies\u2026\u2026. 0.034<\/p>\n\n\n\n

All AOI and road points\u2026\u2026………….. 0.047<\/p>\n\n\n\n

This small project is, of course, not sufficient to statistically determine the best results, but it shows that the results depend much more on the type of point generation than on the presence of anomalies in the Google Street View images.<\/p>\n\n\n\n

Comparison with aerial images.<\/h1>\n\n\n\n

We used the plugin to process both aerial photos with the same algorithm:<\/p>\n\n\n\n

\"\"<\/a><\/figure>\n\n\n\n

The images show the original image and the image with the pixels determined as \u00ab\u00a0green\u00a0\u00bb.<\/p>\n\n\n\n

The GVI of the North image is calculated as 0.1, and that of the South image as 0.07.<\/p>\n\n\n\n

Both are significantly higher than the average GVI calculated on the Google Street View images: 0.1 versus 0.05 and 0.07 versus 0.047 respectively.<\/p>\n\n\n\n

Conclusions<\/h2>\n\n\n\n

This small test does not claim to determine whether the Green View Index plugin gives \u00ab\u00a0better\u00a0\u00bb results than the traditional method of calculating from aerial photos. It simply shows that the results are different. The task now is to understand why and what more they can bring to the knowledge of the urban environment.<\/p>\n","protected":false},"excerpt":{"rendered":"

The QGIS Green View Index plugin is a powerful tool to measure the vegetation cover of a given area. It is interesting to see what it brings compared to the classical processing based on vertical aerial…<\/p>\n","protected":false},"author":1,"featured_media":10421,"comment_status":"open","ping_status":"closed","sticky":false,"template":"","format":"standard","meta":{"give_campaign_id":0,"_bbp_topic_count":0,"_bbp_reply_count":0,"_bbp_total_topic_count":0,"_bbp_total_reply_count":0,"_bbp_voice_count":0,"_bbp_anonymous_reply_count":0,"_bbp_topic_count_hidden":0,"_bbp_reply_count_hidden":0,"_bbp_forum_subforum_count":0,"sfsi_plus_gutenberg_text_before_share":"","sfsi_plus_gutenberg_show_text_before_share":"","sfsi_plus_gutenberg_icon_type":"","sfsi_plus_gutenberg_icon_alignemt":"","sfsi_plus_gutenburg_max_per_row":"","_monsterinsights_skip_tracking":false,"_monsterinsights_sitenote_active":false,"_monsterinsights_sitenote_note":"","_monsterinsights_sitenote_category":0,"_jetpack_newsletter_access":"","_jetpack_dont_email_post_to_subs":false,"_jetpack_newsletter_tier_id":0,"_jetpack_memberships_contains_paywalled_content":false,"_jetpack_memberships_contains_paid_content":false,"footnotes":"","jetpack_post_was_ever_published":false},"categories":[1377],"tags":[1977,1962,1979,1964,1386],"class_list":["post-10419","post","type-post","status-publish","format-standard","has-post-thumbnail","hentry","category-qgis-en","tag-aerial","tag-green-view-index-en","tag-photography","tag-plugins-en","tag-qgis-en"],"aioseo_notices":[],"jetpack_featured_media_url":"https:\/\/i0.wp.com\/www.sigterritoires.fr\/wp-content\/uploads\/2023\/05\/a1-3-2-scaled.jpg?fit=2560%2C1053&ssl=1","jetpack_shortlink":"https:\/\/wp.me\/p6XU0A-2I3","jetpack_sharing_enabled":true,"_links":{"self":[{"href":"https:\/\/www.sigterritoires.fr\/index.php\/wp-json\/wp\/v2\/posts\/10419","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/www.sigterritoires.fr\/index.php\/wp-json\/wp\/v2\/posts"}],"about":[{"href":"https:\/\/www.sigterritoires.fr\/index.php\/wp-json\/wp\/v2\/types\/post"}],"author":[{"embeddable":true,"href":"https:\/\/www.sigterritoires.fr\/index.php\/wp-json\/wp\/v2\/users\/1"}],"replies":[{"embeddable":true,"href":"https:\/\/www.sigterritoires.fr\/index.php\/wp-json\/wp\/v2\/comments?post=10419"}],"version-history":[{"count":0,"href":"https:\/\/www.sigterritoires.fr\/index.php\/wp-json\/wp\/v2\/posts\/10419\/revisions"}],"wp:featuredmedia":[{"embeddable":true,"href":"https:\/\/www.sigterritoires.fr\/index.php\/wp-json\/wp\/v2\/media\/10421"}],"wp:attachment":[{"href":"https:\/\/www.sigterritoires.fr\/index.php\/wp-json\/wp\/v2\/media?parent=10419"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/www.sigterritoires.fr\/index.php\/wp-json\/wp\/v2\/categories?post=10419"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/www.sigterritoires.fr\/index.php\/wp-json\/wp\/v2\/tags?post=10419"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}