plotRGB function with Pix4DFields multispectral orthomosaic

[Max Feldman]

Dear Developers and Users,

I am trying to follow along with the fantastic github tutorial you have prepared. I have performed a UAS flight and captured images using a Micasense RedEdgeMX sensor and generated an orthomosaic file using Pix4DFields.

Here is a link to the complete orthomosaic:

Link to orthomosaic

Here is a link to a pseudo RGB of the orthomosaic

Link to pseudo RGB of orthomosaic

When going through the tutorial I have no problem loading in the stack and the resulting object seems like the appropriate size and dimensions.

5024, 9805,  6

Unfortunately, when I try to visualize the orthomosaic using the plotRGB() function I get a binary image where the orthomosaic is black and the unmeasured space is white (see attached). 

Also when I attempt to crop the image using the fieldCrop() function I get the following error message:

[1] "6 layers available"

Error in .local(.Object, ...) : 

Previously I have no problem working through this tutorial with an orthomosaic generated by DroneDeploy (I believe) from a Parrot Sequoia+ sensor. 

What am I doing wrong? How can I fix this issue?

Thank you for any suggestions you might have.

Max Feldman

[fmatias]

Hi Max,

 

Thank you for asking.

 

I guess to use stretch="lin" with plotRGB() can help you.

 

The function fieldCrop() works well with your image in my computer. 

Could you share your script with us?

EX1 <- stack("test_orthomosaic.data.tif")
plot(EX1) # layer 6 it's not necessary
plotRGB(EX1, stretch="lin")
EX1.Crop <- fieldCrop(mosaic = EX1,plot = T, fast.plot = T) # "fast.plot = T" Fast plotting step

 

Let me know if you have any other question,

Filipe

[Max Feldman]

Happy to report that these suggestions have solved the issue. Thanks Filipe!