Hi,
I want to merge RGB values from raster image to the las pointcloud.
Currently I tried:
las_rgb <- merge_spatial(las, img[[1]], attribute = "r") las_rgb <- merge_spatial(las_rgb, img[[2]], attribute = "g") las_rgb <- merge_spatial(las_rgb, img[[3]], attribute = "b") las_rgb <- add_lasrgb(las_rgb, as.integer(las_rgb$r), as.integer(las_rgb$g), as.integer(las_rgb$b))
But I guess there is a better way, since the intermediate r value seems obsolete. How would I do it properly?
Also it seems strange that plot(las_rgb, color = "RGB) cant find RGB values when I do not run add_lasrgb (since I named them with lower letters), but filter_poi seems to recognize the lower letters as RGB values and renames them to capital RGB, which then results in two attributes for each R, G and B.
las_rgb <- merge_spatial(las, img[[1]], attribute = "r") las_rgb <- merge_spatial(las_rgb, img[[2]], attribute = "g") las_rgb <- merge_spatial(las_rgb, img[[3]], attribute = "b") plot(las_rgb, color = "RGB") # Error: No 'RGB' attributes found. las_rgb <- add_lasrgb(las_rgb, as.integer(las_rgb$r), as.integer(las_rgb$g), as.integer(las_rgb$b)) plot(las_rgb, color = "RGB") # Error in if (colmax > 255) 16 else 8 : # missing value where TRUE/FALSE needed las_rgb <- filter_poi(las_rgb, !is.na(R) & !is.na(G) & !is.na(B)) # Attribute 'r' renamed 'R' to match with default attribute names. # Attribute 'g' renamed 'G' to match with default attribute names. # Attribute 'b' renamed 'B' to match with default attribute names. str(las_rgb@data) # Classes ‘data.table’ and 'data.frame': 336191 obs. of 23 variables: # $ X : num 567275 567275 567274 567274 567274 ... # $ Y : num 5708885 5708885 5708885 5708885 5708885 ... # $ Z : num 0 0.998 0 1.313 -0.025 ... # $ gpstime : num 567074 567074 567074 567074 567074 ... # $ Intensity : int 176 20 136 26 149 117 12 12 113 23 ... # $ ReturnNumber : int 5 5 6 4 5 3 3 4 5 2 ... # $ NumberOfReturns : int 5 6 6 5 5 3 3 5 5 4 ... # $ ScanDirectionFlag: int 0 0 0 0 0 0 0 0 0 0 ... # $ EdgeOfFlightline : int 0 0 0 0 0 0 0 0 0 0 ... # $ Classification : int 2 13 2 13 20 20 13 13 20 13 ... # $ Synthetic_flag : logi FALSE FALSE FALSE FALSE FALSE FALSE ... # $ Keypoint_flag : logi FALSE FALSE FALSE FALSE FALSE FALSE ... # $ Withheld_flag : logi FALSE FALSE FALSE FALSE FALSE FALSE ... # $ ScanAngleRank : int 30 30 30 30 30 30 30 30 30 30 ... # $ UserData : int 50 47 51 53 53 49 122 61 55 54 ... # $ PointSourceID : int 150 150 150 150 150 150 150 150 150 150 ... # $ Zref : num 322 323 322 323 322 ... # $ R : int 104 104 81 108 88 84 108 92 68 81 ... # $ G : int 96 96 73 100 80 72 100 84 55 73 ... # $ B : int 85 85 60 87 67 60 87 71 46 60 ... # $ R : int 26728 26728 20817 27756 22616 21588 27756 23644 17476 20817 ... # $ G : int 24672 24672 18761 25700 20560 18504 25700 21588 14135 18761 ... # $ B : int 21845 21845 15420 22359 17219 15420 22359 18247 11822 15420 ... # - attr(*, ".internal.selfref")=<externalptr>
RetroSearch is an open source project built by @garambo | Open a GitHub Issue
Search and Browse the WWW like it's 1997 | Search results from DuckDuckGo
HTML:
3.2
| Encoding:
UTF-8
| Version:
0.7.4