terraref / reference-data

Coordination of Data Products and Standards for TERRA reference data
https://terraref.org
BSD 3-Clause "New" or "Revised" License
9 stars 2 forks source link

Test scans for input on VNIR settings #248

Open NewcombMaria opened 6 years ago

NewcombMaria commented 6 years ago

Today @jdemieville-ua and I ran multiple VNIR test scans across the short axis of the field with targets deployed at a range of exposure settings, all at 0.04 m/sec velocity (y axis speed). The majority of the test runs were in triggered mode, one in untriggered mode. Exposures ranged from 15 to 35 milliseconds. Conditions were mostly cloudy during our test scan, but we had to proceed without waiting for a sunny day. We can repeat the tests again on day without clouds after planting. Photo copied below with the primary targets labeled. A couple potted sorghum plants were included.

@remotesensinglab would you be able to advise us on settings?
Let me know if you need additional information.

image

Paheding commented 6 years ago

@NewcombMaria and @jdemieville-ua Thanks for running the test scans. Were that VNIR data stored in Globus? If yes, could you please let me know the folder names? So that we can look into these specific data and find out the proper settings.

jdemieville-ua commented 6 years ago

@Paheding : It doesn't look like they've made it into Globus yet. I'll keep an eye on it and update you with a list of directories when they're available.

max-zilla commented 6 years ago

This issue is why they haven't been transferred via Globus yet: https://github.com/terraref/computing-pipeline/issues/492

We have a ticket open with Globus, the endpoint is reporting an expired certificate even though JD looked on the cache server in AZ and it reports the cert won't expire until 2021. until this is resolved, transfers can't occur.

Paheding commented 6 years ago

@jdemieville-ua and @max-zilla Thanks for the info. I look forward to seeing the images.

remotesensinglab commented 6 years ago

All, Patrick and I evaluated the data. We could run the quality check pipeline code on these as well to see if we can warn the quality issues automatically. But wanted to provide a quick note as we prepare for data collection. The report is available here. HSI_evaluation.pdf

I think the 20 ms exposure time with 50 frame rate is the best among the settings: (1) Reflectance at NIR is under or equal to 0.55 which is expected for plants (2) reflectance for 95% panel is where it is supposed to be

THe remaining issue is there is still distortion on the data. Pixels is not square. It is not a major and can be fixed either by changing the frame rate (reducing it) or increasing speed a little bit.

Vasit

max-zilla commented 6 years ago

@Paheding here are the datasets we transferred for VNIR from 8/18 and 8/22:

screen shot 2018-08-28 at 8 34 09 am

screen shot 2018-08-28 at 8 34 22 am

...this matches the files and sizes on the gantry side. For the 22nd the raw data goes from 30 GB on the first image to 20 MB on the last, but that matches the source.

remotesensinglab commented 6 years ago

Ok. We will review these ASAP and report back to the team. thanks

Paheding commented 6 years ago

@max-zilla Thanks for letting us know the status. The above-mentioned evaluation report was based on some of those datasets, we will do the evaluation on the rest that has been transferred completely.

rjstrand commented 6 years ago

With respect to triggered acquisition, we have no real frame rate. Frames are triggered based on position. Position is based on the motion of the camera box. We can change the trigger rate by either changing the velocity of the camera box, or the division rate in the frequency divider on the trigger circuit.

The Headwall frame period setting has no impact on the triggered acquisition. Nevertheless, it must be set to a value that is larger than the exposure time.

I will verity with Headwall, but I believe that exposure time cannot simply be set to a value of 1/trigger rate. There has to be some "processing" time between triggers. This is why we ran the VNIR with with a 35 ms exposure at a 50 ms trigger rate.

The y-axis encoder signal is currently divide down so that the cameras trigger at a rate of approximately 1 scan per 0.95 mm (As I remember, currently dividing primary pulse train by 65). We can modify this by changing the divisor on the frequency divider.

remotesensinglab commented 6 years ago

Hi Bob, Do not change anything until we evaluate other datasets. It seems there are more data than we evaluated that have been transferred to the server recently.

We will look at all systematically and get back to the team next day or two.

Thanks

Vasit

Paheding commented 6 years ago

Hi all,

We have run the image quality evaluation algorithm on the available VNIR test scans, and conducted some analysis on the spectral profiles of plant and white pannel. All results are attached here (HSI_evaluation.pdf). For convenience, a short summary is as follows:

image

image

Remaining issue: The geometrical distortion on the pixels, see the shapes of plant and color board in the following image:

image

dlebauer commented 6 years ago

Regarding the geometrical distortion - is this file on globus where I can take a look? If you are referring to the fact that the round objects appear oval, isn't this just a reflection of the fact that the size of the pixel in the left-right direction is determined by the camera and the scan direction (up) is controlled by the speed of the camera? i.e. is this an artifact of trying to fit rectangular information into a square raster image, and if we have correct geospatial reference does this distortion affect the interpretation of the results?

It isn't clear how the plants are distorted, but we have always had waviness that I've assumed is due to wind blowing the plants back and forth as the scan line is pushed forward.

Other teams have mentioned that geometric calibration is a uniquely challenging task, but that it may not directly impact many of the key features we are trying to measure.

What impacts will this distortion have on biological inference?

dlebauer commented 6 years ago

a few minor details on the soil mask data products ... the way the files are currently generated, the metadata is incomplete (if you compare the output of ncdump -h file.nc on a L1 vnir file compared to this, you will see that each variable in the L1 file has a number of attributes including name, long name, units, description, notes, etc. Also, the original dimension 'x' and 'y' have been replaced with 'x_dim' and 'y_dim' which is confusing, and in the soil mask, 'x', 'y', 'latitude', and 'logitude' are all empty.

It would be great to have a soil mask data product, but as it stands this would need some work to clean this up. NCO has a lot of utilities that should make it 'easy' to copy x,y,latitude,longitude from the input file to the output file, delete any obsolete variables, and add appropriate metadata (perhaps the python library does as well).

NewcombMaria commented 6 years ago

Thanks everyone for helpful input. @Paheding and @remotesensinglab the report and analysis of VNIR results are great. Bob Strand also passed on info received from Headwall: in triggered mode exposure must be less than the trigger period, otherwise the trigger that occurs during the exposure is ignored/skipped. This applies to VNIR and SWIR.

@jdemieville-ua wrote scan scripts with the VNIR and SWIR set in triggered mode at exposures 20ms for both. We've been running short hyperspectral scans starting mid-afternoon on 9/1 and 9/2 and today to collect data on soil surface moisture in proximity of the soil moisture probes. The data are also useful for continued tests of settings. At this time of the season the sun angle allows partial sun in afternoon images, but there's also unavoidable partial shade. Seedling emergence occurred on 8/31 so there are seedlings present in the images.

Paheding commented 6 years ago

@dlebauer The files are located at Globus-- /ua-mac/raw_data/VNIR/2018-08-22/. Typically, the distoration can be due to mirror scans, velocity variance, weather conditions, etc. I am not quite sure if the wind blowing occurs at that date that caused the distoration. Based on the experiences from UAV based hypersepecral image collection, the pixel distoration can be corrected by properly adjusting some parameters such as framerate, flight height, etc.

Paheding commented 6 years ago

@dlebauer The soil mask extractor (https://github.com/terraref/extractors-sunshade) will process a netCDF data file and export a binary mask of non-vegeation area from hyperspectral imagery, and generates background (typically soil) removed hyperspectral imagery with the associated NDVI value as well as soil fraction area. The georeference points are also added to the meta data. These outputs were dicussed and suggested by team. We can definitely add more variables if they are useful.

Paheding commented 6 years ago

@NewcombMaria We glad to hear the analytic report is helpful. And thanks again for collecting more data. Especially the data containing seedling emergence and lighting condition would help to robust analysis on the spatial and spectral content. We will provide more analytic results on the coming data and share with the team.

SeanHartling commented 6 years ago

All,

Just wanted to share some test data we gathered to examine the exposure settings for hyperspectral collection. The data below was collected from our UAV platform with Headwall Nano VNIR hyperspectral camera. Despite our best timing attempts, there were intermittent clouds during the collect. The two left scenes are from the same flight path. A-1 is the start of the mission, then A-2 is in the same flight line (arrows indicate direction of travel) and B is next flight line in the opposite direction. You can see the effect of lighting conditions and exposure between the three scenes. The right two scenes demonstrate the effect of exposure and lighting conditions on reflectance. All the reflectance values on the reflectance tarp are higher on scene B compared to scene A. This should be something to consider when dealing with a fixed exposure rate. We will continue testing to figure out the best approach for this phenomenon. headwall_exposure_reflectance_example

Paheding commented 6 years ago

@dlebauer Althought the correction on the pixel distoration is not current priority, we have quick-checked several VNIR data that were previously captured, It looks like the spatial resolution is better than recent one. A sample image is shown below:

Image tag: 2018-02-01__15-05-34-160 image

However, we need to check the further upcoming data in case of the pixel distortation is due to the wind blowing or any other weather conditions, as you suggested. Another possible reason would be some limitations of this new VNIR sensor compared to the previous one.

NewcombMaria commented 6 years ago

Thanks Patrick and Vasit for input on exposure settings. Now that the sun angle is lower and there is more sunlight, we are wondering if the 20 ms exposure setting is best for bright sunlight, or if 15 ms exposure for VNIR and SWIR may be better when the image is in full sun.

We collected data on 2018-09-22 to compare the 20 ms and 15 ms exposure settings.
SWIR 15 ms: /ua-mac/raw_data/SWIR/2018-09-22/2018-09-22__13-21-35-977/ SWIR 20 ms: /ua-mac/raw_data/SWIR/2018-09-22/2018-09-22__13-34-50-596/

VNIR 15 ms: /ua-mac/raw_data/VNIR/2018-09-22/2018-09-22__13-21-55-992/ VNIR 20 ms: /ua-mac/raw_data/VNIR/2018-09-22/2018-09-22__13-35-10-611/

@remotesensinglab would your group be able to evaluate and let us know which exposure is better now that the sun angle allows more sunlight under the camera box?

Paheding commented 6 years ago

@NewcombMaria Thanks for conducting experiments. We will evaluate those data and provide the feedback.

Paheding commented 5 years ago

@NewcombMaria We have conducted some experments. The results suggest that those two settings do not make too much difference. However, 20 ms may be slightly better for both VNIR and SWIR cases. See the graphs below in details.

image

image