QGISã§æš¹æšå°åž³ãäœãæ¹æ³ïŒèªç©ºåçãšLiDARãã1æ¬ãã€ç»é²
å ¬éæ¥
QGISã§ã®æš¹æšå°åž³ãšã¯ã1æ¬ã®æšã«ã€ã1ç¹ã眮ãããã€ã³ãã¬ã€ã€ã§ãåç¹ã«IDãæš¹çš®ãå¹¹ã®å€ªããç¶æ ãæçµç¹æ€æ¥ãšãã£ã屿§ã®è¡ãä»ãããã®ã§ããèªç©ºåçãLiDARãããã°ãäœçœ®ãæš¹å ãµã€ãºãæš¹é«ã¯æºäžã§åŸãããŸããæš¹çš®ãå¹¹ãç¶æ ã¯çŸå°èª¿æ»ã§èšé²ããŸãã
ãã®åœ¹å²åæ ããå€ãã®ã¬ã€ãã¯èª¬æããŸããã20 cmã®ãªã«ãœãã©ãããæš¹çš®ãè æœããå¹¹ãèªã¿åããã¢ãã«ã¯ãããŸãããã ããæãæ©ãå°åž³ã¯ãæ°ããã®ãæºäžã§ã倿ããã®ãçŸå°ã§è¡ããã®ã§ãã
çµè«
| åºçºç¹ | æ¹æ³ | æºäžã§åŸããããã® |
|---|---|---|
| èªæ²»äœãçºæ³šè ããã§ã«å°åž³ãå ¬éããŠãã | æ¹æ³1ïŒãªãŒãã³ãªæš¹æšããŒã¿ | 暹皮ãå€ãã¯å¹¹ã®å€ªãä»ãã®ãã€ã³ã |
| ãªã«ãœãã©ããšæ°çŸæ¬ã®æš | æ¹æ³2ïŒãã€ã³ããæåã§ããžã¿ã€ãº | èªåã§çœ®ãããã€ã³ã |
| LiDARã®è¡šé¢ã¢ãã«ãšå°åœ¢ã¢ãã« | æ¹æ³3ïŒæš¹å é«ã¢ãã«ããæš¹é ç¹ãæœåº | æš¹é«ä»ãã®ãã€ã³ã |
| ãªã«ãœãã©ããšæ°åæ¬ã®æš | æ¹æ³4ïŒã»ã°ã¡ã³ããŒã·ã§ã³ã¢ãã« | æšããšã«1ã€ã®æš¹å ããªãŽã³ãããããæš¹å ãµã€ãºä»ãã®ãã€ã³ã |
4ã€ãšãæåŸã¯åãã§ããã¬ã€ã€ãã¿ãã¬ããã§çŸå°ã«æã¡åºããç»åã§ã¯åãããªãæ å ±ã誰ãã远å ããŸãã
æš¹æšå°åž³ã«èšé²ããå 容
ãã€ã³ãã1ã€çœ®ãåã«ã屿§ããŒãã«ã決ããŠãã ãããçŸå°èª¿æ»ã®ããšã§å€æŽãããšãããäžåºŠçŸå°ã«è¡ãããšã«ãªããŸãã
| 屿§ | æºäžã§ | çŸå°ã§ |
|---|---|---|
| æš¹æšID | å¯ãèªåçæ | æšã«çªå·æãããã°ããã®çªå· |
| äœçœ® | å¯ | æš¹å ã§å¹¹ãé ããŠããå Žåã¯ä¿®æ£ |
| æš¹å ã®çŽåŸãŸãã¯é¢ç© | å¯ãç»åãã | ç¢ºèª |
| æš¹é« | å¯ãLiDARãããã° | LiDARããªããã°å®æž¬ |
| 暹皮 | äžå¯ | å¯ |
| å¹¹ã®çŽåŸãŸãã¯åšå²é· | äžå¯ | å¯ãåºæºãå®ããé«ãã§ |
| ç¶æ ãæ¬ é¥ã次åç¹æ€ | äžå¯ | å¯ |
ãã€ãã®Baumkatasterã¯ãå®å šé æ ®çŸ©åã«åºã¥ãç¹æ€ãèšé²ãããã®ã§ãå€ãã®ææè ã¯FLLã®æš¹æšç¹æ€ã¬ã€ãã©ã€ã³ïŒ2020幎çïŒãšGALKã®æš¹æšå°åž³ã¬ã€ãã«åŸã£ãŠããŸãããããã®ææžãçŸå°ã§ã®èšé²é ç®ãå®ããŠããã®ã§ãäžã®è¡šã§ã¯ãªãããŸããã¡ãããå§ããŠãã ããã
ãŠã§ãã®æ å ±ã§èª€ããããç¹
ãAIã¯ç»åã ãã§å°åž³ããã¹ãŠäœãããã AIãæ°ããŠæž¬ããã®ã¯æš¹å ã§ããæš¹å ã®äžã®å¹¹ã¯èŠãããRGBåçãã倧åã®æš¹çš®ãå€å¥ããããšããæãèœã¡ããã©ããã®å€æãã§ããŸãããé£è¡1åã§å®æããå°åž³ãã§ãããšçŽæããããŒãžã¯ãäœæ¥ã®ååãã売ã£ãŠããŸããã
ãæš¹å ãããã¯æš¹æšå°åž³ã§ãããã æš¹å ãããã¯é¢ã§è¡šãããã®ã§ãããã®å°åºã®äœå²ãæ¥é°ããã«çããŸããå°åž³ã¯ãã©ã®æšããã©ãã«ãã£ãŠãã©ããªç¶æ ããã«çããŸããå¿ èŠãªã®ã被èŠçãªããQGISã§æš¹å ãå°å³åããã¬ã€ããé©åã§ãã
å§ããåã«
- QGIS 3.22以éã
- ä»å¹Žã®æšãåã£ãŠããæ°ããã®ç»åãèã®ããææã®ç»åã¯æš¹å ãæãèŠããããèœèæã®ç»åã¯å¹¹ãšæåãããèŠããããªããŸãã
- ã¡ãŒãã«åäœã®æåœ±åº§æšç³»ãé¢ç©ãšçŽåŸãã¡ãŒãã«ã§åºãŸãã
ç¡æã®ç»åã¯ãå€ãã®äººãæã以äžã«å å®ããŠããŸãã
- ããŒã©ã³ãïŒ Geoportalã®å šåœãªã«ãœãã©ãã¯ã2020幎6æ24æ¥ããçšéãåããç¡æã§ããè§£å床ã¯å€ãã25 cmãéœåžéšã¯10 cm以äžã§ããç¡æã®ãPobieracz danych GUGiKããã©ã°ã€ã³ã§ãQGISã®äžããããŠã³ããŒãã§ããŸãã
- ãã€ãïŒ ããŒãã³ã»ãŽã¥ã«ãã³ãã«ã¯å·ãããŒããŒã¶ã¯ã»ã³å·ãªã©ãããã€ãã®å·ã20 cmã®ãªã«ãœãã©ãã§ããDOP20ããªãŒãã³ããŒã¿ãšããŠå ¬éããŠããŸããé£éŠæ©é¢ã®ãã€ãå šååããµãŒãã¹ã¯ææã§ãã
- ã€ã³ã°ã©ã³ãïŒ ç°å¢åºããã»ãŒå šåã®1 mã®LiDAR衚é¢ã¢ãã«ãšå°åœ¢ã¢ãã«ããOpen Government Licenceã§å ¬éããŠããŸãã
- ã¢ã¡ãªã«ïŒ NAIPã®èªç©ºåçã¯ãUSDAããç¡æã§å ¥æã§ããŸãã
æ¹æ³1ïŒãªãŒãã³ãªæš¹æšããŒã¿
ãã§ã«å°åž³ãéçšããŠå ¬éããŠããèªæ²»äœã¯å€ããããŸãããã«ãªã³ã¯è¡è·¯æš¹ãšå ¬åã®æšãå ¬éãããã«ãã ã³ãã¯èªãã®ããŒã¿ãšèªç©ºåçããäœã£ãå°åž³ãå ¬éããã±ã«ã³ã¯2025幎7æã«å°åž³ããªã³ã©ã€ã³ã«åºããããªã¯219,432æ¬ããã€ã³ãã§æäŸããŠããŸãããã¡ã€ã«ããµãŒãã¹ãQGISã§èªã¿èŸŒãã°ã暹皮ããã§ã«å ¥ã£ãç¶æ ã§å§ããããŸãã
OpenStreetMapã«ããnatural=treeã®ãã€ã³ããšããŠåã
ã®æšãèŒã£ãŠããŸããQGISã§OSMããŒã¿ãããŠã³ããŒãããã§ãQuickOSMã䜿ã£ãååŸæ¹æ³ã説æããŠããŸãã
åããŠããã±ãŒã¹ïŒ æšãåžã®ç®¡çã§ãå°åž³ãæ°ãããšãã
åããªãã±ãŒã¹ïŒ ç§æå°ããã£ã³ãã¹ãæ°ããéçºå°ãæ±ããšãããªãŒãã³ãªå°åž³ã¯å ¬å ±ã®æšãŸã§ã§ãã
æ¹æ³2ïŒãã€ã³ããæåã§ããžã¿ã€ãº
Layer > Create Layer > New GeoPackage LayerâŠã§ã¬ã€ã€ãäœãããžãªã¡ããªã¿ã€ããPointã«ããŠãäžã®è¡šã®é
ç®ã远å ããŸããç·šéãæå¹ã«ãããªã«ãœãã©ãäžã§åæš¹å ã®äžå¿ã1åãã€ã¯ãªãã¯ããŸãã
æºäžã®æ¹æ³ã§ã¯æãæ£ç¢ºã§ãããæãé ããªããŸãã1ç¹ãã€èªåã§å€æããã®ã§ãçãå£ã12æ¬ã®æšã«ãªãããšããæ¥ãã2ã€ã®æš¹å ã1æ¬ã«ãªãããšããããŸããã
åããŠããã±ãŒã¹ïŒ æ°çŸæ¬çšåºŠã®ãšãããŸãã¯éçºèš±å¯ç³è«ã®ããã«1æ¬1æ¬ãéèŠãªçŸå Žã®ãšãã
åããªãã±ãŒã¹ïŒ æ°ãæ°åã«åã¶ãšãã
æ¹æ³3ïŒæš¹å é«ã¢ãã«ããæš¹é ç¹ãæœåº
LiDARãããã°ãåæšã®é ç¹ã¯åšèŸºã§æãé«ããã¯ã»ã«ã«ãªãã®ã§ãåŠç¿ãªãã§é ç¹ãèŠã€ããããŸãããŸãã衚é¢ã¢ãã«ããå°åœ¢ãåŒããæš¹å é«ã¢ãã«ããæš¹å ã¬ã€ãã®æé ã§äœã£ãŠãã ããã2ã€ã®å ¥åããŒã¿ã¯DSMãšDTMã®è§£èª¬ã§èª¬æããŠããŸãã
屿æå€§å€ãæ±ãã
æš¹å é«ã¢ãã«ã«å¯ŸããŠr.neighborsïŒGRASSïŒãå®è¡ããæ¹æ³ãmaximumãè¿åã5ãã¯ã»ã«ã«ããŸãã1 mã®ãã¯ã»ã«ãªã5 mã®çªã§ãã
æš¹é ç¹ãæ®ã
Raster > Raster CalculatorâŠã§ã屿æå€§å€ãšçããã3 mããé«ããã¯ã»ã«ãæ®ããŸãã("chm@1" = "chm_max@1") AND ("chm@1" > 3)
ãã¯ã»ã«ããã€ã³ãã«å€æ
çµæãRaster pixels to pointsã§å€æããExtract by attributeã§å€ã1ã®ãã€ã³ããæ®ããŸãã
æš¹é«ã远å
æš¹å é«ã¢ãã«ã«å¯ŸããŠSample raster valuesã䜿ããåãã€ã³ãã«æš¹é«ã远å ããŸãã
éèŠãªèšå®ã¯çªã®ãµã€ãºã§ããå°ãããããšã1æ¬ã®å€§ããªãªãŒã¯ãã3ã€ã®é ç¹ãåºãŸãã倧ãããããšã2æ¬ã®å°ããªæšã1ã€ã«ãªããŸããçŸå Žã®å žåçãªæš¹å å¹ ã«åããããµã³ãã«ãç®ã§ç¢ºèªããŠãã ããã
åããŠããã±ãŒã¹ïŒ LiDARããã£ãŠãæš¹é«ãå¿ èŠãªãšãã
åããªãã±ãŒã¹ïŒ ç»åãããªããšãããŸãã¯LiDARãæ°ããæšããå€ããšãã
æ¹æ³4ïŒã»ã°ã¡ã³ããŒã·ã§ã³ã¢ãã«
ã»ã°ã¡ã³ããŒã·ã§ã³ã¢ãã«ã¯æš¹å ã®èŠãç®ãåŠç¿æžã¿ã§ãéåžžã®RGBãªã«ãœãã©ãããæšããšã«1ã€ã®ããªãŽã³ãè¿ããŸããåœç€Ÿã¯AI Segmentationãäœã£ãŠããã®ã§ããã®ç¯ã¯èªç€ŸããŒã«ã®èª¬æã§ãã
QGISã§éãããªã«ãœãã©ãã®äžã«ãŸãŒã³ãæããtreeãšå
¥åããŸããã¢ãã«ã¯åœç€Ÿã®ãµãŒããŒã§åãã®ã§ãã°ã©ãã£ãã¯ã¹ã«ãŒãã¯èŠããŸãããåæš¹å ã¯ããããžã§ã¯ãã®CRSã®ããªãŽã³ãšããŠãarea_m2ãšä¿¡é ŒåºŠã®scoreä»ãã§è¿ã£ãŠããŸãã

æš¹å ãå°åž³ã®ãã€ã³ãã«ããã«ã¯ã次ã®ããã«ããŸãã
æš¹å ããšã«1ç¹ãäœã
æš¹å ã¬ã€ã€ã«Centroidsãå®è¡ããŸããåãã€ã³ãã¯ãå
ã®æš¹å ã®å±æ§ãåŒãç¶ããŸãã
æš¹å ã®çŽåŸã远å
ãã£ãŒã«ãèšç®æ©ã§ãå°æ°åã®ãã£ãŒã«ããšããŠæš¹å ã®çŽåŸã远å ããŸãã2 * sqrt("area_m2" / pi())
æšã«çªå·ãæ¯ã
@row_numberãå
¥ããtree_idãã£ãŒã«ãã远å ããç¶ããŠå±æ§ããŒãã«ã§æ±ºããçŸå°çšã®ç©ºã®ãã£ãŒã«ãã远å ããŸãã
åããŠããã±ãŒã¹ïŒ ãã§ã«æã£ãŠããç»åãããéããå ¬åãåºã®æ°åæ¬ãæ±ããšãã
åããªãã±ãŒã¹ïŒ æãéããŠãããšããæš¹å ã®éã«ããéããªããšãtreeã®å®è¡ã¯å¹¹ããšã®ããªãŽã³ã§ã¯ãªã1ã€ã®å€§ããªåœ¢ãè¿ããŸãããã®å Žåã¯æ¹æ³3ã䜿ã£ãŠãã ãããç»åã¿ã€ã«ã¯åœç€Ÿã®ãµãŒããŒã«éãããŸãããã·ã³ã®å€ã«åºããªãå Žåã¯ãæ¹æ³2ã3ã䜿ã£ãŠãã ããã
çŸå°ã§ä»äžãã
ãã€ã³ãã¬ã€ã€ã¯ãQGISãããžã§ã¯ãåãã®ç¡æã¢ãã€ã«ã¢ããªQFieldïŒOPENGIS.ch補ïŒã§çŸå°ã«æã¡åºããŸãããªãã©ã€ã³ã§åããQGISã§èšå®ãããã©ãŒã ããã®ãŸãŸéããæšããšã«åçã远å ã§ããŸããåºçºåã«QGISã§æš¹çš®ãšç¶æ ã®ããããããŠã³ãèšå®ããŠããã°ã調æ»å¡å šå¡ãåãèšèã§èšå ¥ã§ããŸãã
æ³çãªèšé²ã䌎ãå®æç¹æ€ã«ã¯ãå°çšã®å°åž³ãœããããããŸãããã€ãã®Baumsichtã®ããã«QGISäžã§åããã®ããããŸããããã§äœã£ããã€ã³ãã¬ã€ã€ã¯ããã¡ãã«åã蟌ããŸãã
èŠããŠããããš
- æš¹æšå°åž³ã¯ãæš¹å ã®é¢ã§ã¯ãªãã屿§ä»ãã®1æ¬1ç¹ã®ãã€ã³ãã§ãã
- æºäžã§åŸãããã®ã¯äœçœ®ãšæš¹å ãµã€ãºãLiDARãããã°æš¹é«ã§ããæš¹çš®ãå¹¹ãç¶æ ã¯çŸå°ã§åŸãŸãã
- ãŸããªãŒãã³ãªå°åž³ãæ°çŸæ¬ãªãæåã®ããžã¿ã€ãºãæš¹é«ãå¿ èŠãªãLiDARã®æš¹é ç¹ãç»åäžã®æ°åæ¬ãªãã»ã°ã¡ã³ããŒã·ã§ã³ã¢ãã«ã䜿ããŸãã
- ã©ã®èªåææ³ããæã§æ°ãããµã³ãã«ã§ç¢ºèªããŠããåèšãä¿¡çšããŠãã ããã
ãããã質å
æš¹æšå°åž³ã¯çŸ©åã§ããïŒ
å°åž³ãã®ãã®ãæ±ããæ³åŸã¯ãéåžžãããŸãããéè·¯ãæ©éã«é£æ¥ããæšã®ææè ã«ã¯å®å šã®ããã®ç¹æ€çŸ©åããããå°åž³ã¯ãã®ç¹æ€ã宿œãã蚌æ ã«ãªããŸãããã€ãã§ã¯ããã®å®å šé æ ®çŸ©åããå€ãã®èªæ²»äœãBaumkatasterãæã€çç±ã§ãã
èªç©ºåçããAIã§æš¹çš®ãå€å¥ã§ããŸããïŒ
å€ãã®äººãæã£ãŠããRGBã®ãªã«ãœãã©ãããã¯ã確å®ã«ã¯å€å¥ã§ããŸãããæš¹çš®ã®å°å³åã«ã¯ãã«ãã¹ãã¯ãã«ãŸãã¯ãã€ããŒã¹ãã¯ãã«ã®ç»åãšå°åã®æåž«ããŒã¿ãå¿ èŠã§ãããã§ãçŸå°ã§ç¢ºèªããŸããåœç€Ÿã®ãã©ã°ã€ã³ãè¿ãã®ã¯æš¹çš®ã§ã¯ãªãæš¹å ã§ãã
èªåã®æš¹æšæ€åºã¯ã©ããããæ£ç¢ºã§ããïŒ
ããŒã«ãããæš¹å ã®ç¶æ ã«å·Šå³ãããŸããå€ç«ããè¡è·¯æš¹ãåºæšã¯ããŸãåºãŸãããæ¥ããæš¹å ã¯1ã€ã«åããããŸããèªåã®ç»åã§50æ¬ãæã§æ°ããŠã¬ã€ã€ãšæ¯ã¹ãã°ããã®çŸå Žã§ã®èª€å·®çãåãããŸãã
QGISã§ç¡æã§æš¹æšå°åž³ãäœããŸããïŒ
ã¯ãããªãŒãã³ãªå°åž³ãæåã®ããžã¿ã€ãºãLiDARã®æš¹é ç¹ã¯ãæé以å€ã®è²»çšãããããŸãããQFieldãç¡æã§ããã»ã°ã¡ã³ããŒã·ã§ã³ã¢ãã«ã¯ã倧éã®æ¬æ°ã§æéãç¯çŽã§ããç¡æã§å§ããããŸãã
é¢é£ããäœæ¥ãšããŠããªã«ãœãã©ããšã¯äœãã§ãããã§æ±ã£ããã¹ãŠã®æ¹æ³ã®å ã«ãªãç»åã説æããŠããŸããåãã©ã°ã€ã³ã®åäœæ¡ä»¶ã¯QGIS AIãã©ã°ã€ã³ã®ããŒãžã§ç¢ºèªã§ããŸãã
ã§ã®æš¹æšå°åž³ãšã¯ã1æ¬ã®æšã«ã€ã1ç¹ã眮ãããã€ã³ãã¬ã€ã€ã§ãåç¹ã«IDãæš¹çš®ãå¹¹ã®å€ªããç¶æ
ãæçµç¹æ€æ¥ãšãã£ã屿§ã®è¡ãä»ãããã®ã§ããèªç©ºåçãLiDARãããã°ãäœçœ®ãæš¹å ãµã€ãºãæš¹é«ã¯æºäžã§åŸãããŸããæš¹çš®ãå¹¹ãç¶æ
ã¯çŸå°èª¿æ»ã§èšé²ããŸãã

