Description
Hello, I'm using a RedEdge MX as a ground sensor and I would like to take into account also DLS for reflectance calibration.
Citing tutorial n. 3
"we assume a constant irradiance (over time) and a ratio of direct to diffuse irradiance of 6:1. These assumptions only hold approximately, and only for clear sky. Cloud cover can lead to dramatically different results"
and also
"If your intent is to simply use this code to correct your data, you can skip this section and access pose-corrected DLS data through the capture.dls_irradiance() method".
My questions are: does the function capture.dls_irradiance() give me a value that is not affected by the choice of the ratio of direct to diffuse irradiance? In other words, the values that this function gives me, can be used directly in the subsequent part of code, that I report her below:
panel_radiances = np.array(cap.panel_radiance())
irr_from_panel = math.pi * panel_radiances / panel_reflectance_by_band
dls_correction = irr_from_panel/dls_irradiances
cap.plot_undistorted_reflectance(dls_irradiances*dls_correction)
For the correcting calibration of reflectance of the other captures should I go with:
dls_irradiances=capture_to_analyze.dls_irradiance()
capture_to_analyze.undistorted_reflectance(dls_irradiances*dls_correction)
Moreover, every time I tried tutorial n.3 in this part of code
import numpy as np
import micasense.dls as dls
import os, glob
import micasense.capture as capture
images_path = os.path.join('.','data','0000SET','000')
image_names = glob.glob(os.path.join(images_path,'IMG_0000_*.tif'))
cap = capture.Capture.from_filelist(image_names)
# set panel corners manually if zbar is not installed
panelCorners = [[[809,613],[648,615],[646,454],[808,452]],
[[772,623],[613,625],[610,464],[770,462]],
[[771,651],[611,653],[610,492],[770,490]],
[[829,658],[668,659],[668,496],[829,496]],
[[807,632],[648,634],[645,473],[805,471]]]
cap.set_panelCorners(panelCorners)
panel_radiances = np.array(cap.panel_radiance())
gave me this error
Cell In[285], line 1
np.array(panelCap.panel_radiance())
File ~\anaconda3\envs\micasense\lib\site-packages\micasense\capture.py:396 in panel_radiance
mean, _, _, _ = p.radiance()
File ~\anaconda3\envs\micasense\lib\site-packages\micasense\panel.py:268 in radiance
return self.region_stats(radiance_img,
File ~\anaconda3\envs\micasense\lib\site-packages\micasense\panel.py:239 in region_stats
rev_panel_pts = np.fliplr(region) # skimage and opencv coords are reversed
File <__array_function__ internals>:200 in fliplr
File ~\anaconda3\envs\micasense\lib\site-packages\numpy\lib\twodim_base.py:98 in fliplr
raise ValueError("Input must be >= 2-d.")
ValueError: Input must be >= 2-d.
What am I doing wrong? Thanks