Skip to content

Review plant_height script output #210

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Open
4 of 8 tasks
craig-willis opened this issue Dec 4, 2017 · 16 comments
Open
4 of 8 tasks

Review plant_height script output #210

craig-willis opened this issue Dec 4, 2017 · 16 comments
Assignees

Comments

@craig-willis
Copy link

craig-willis commented Dec 4, 2017

Review task for plant_height script and output:
https://github.com/terraref/extractors-3dscanner/blob/master/plant_height/laser_to_height.sh

We seem to have lost track of this review task in the extractor repo:
terraref/extractors-3dscanner#9 (comment)

Note:

  • This script is not an extractor
  • The script operates on the raw PLY data (in Level 1), not the merged LAS data

Completion criteria:

  • Code review complete (no extractor, so just script)
  • Sample output data generated
  • Sample output data verified
  • Changes for 2017 release implemented
  • Issues for 2018 created
  • Code tagged/versioned
  • Script executed
  • Final output verified/reviewed
@craig-willis
Copy link
Author

Per terraref/extractors-3dscanner#9 (comment):

The first end-to-end test ran over the weekend on ROGER for 2017-06-02. Please check that the
Level_1 and trait data are as expected.

Data is on ROGER:
/projects/arpae/terraref/sites/ua-mac/Level_1/canopy_height/

@dlebauer
Copy link
Member

dlebauer commented Dec 4, 2017

(from terraref/extractors-3dscanner#9 (comment))

  • in the Level 1 data, I see 359 unique dates under scanner3DTop and 63 unique dates under scanner3DTop_mergedlas
  • in BETYdb, there appear to be some very small values from May 30 through June 2 2017

image

@ZongyangLi
Copy link

@dlebauer
Is that means only some of data points have such a very small and negative values? PI643008 is a gene id?

@craig-willis
Copy link
Author

@dlebauer What is the relation to mergedlas? I think this script operates only on the PLY files

@dlebauer
Copy link
Member

dlebauer commented Dec 4, 2017

@ZongyangLi PI643008 is a cultivar (== genotype) id.

@craig-willis sorry I missed that.

@dlebauer
Copy link
Member

dlebauer commented Dec 4, 2017

@ZongyangLi It looks like ther are no longer any records < 0 but all of the records with the 'Scanner 3d ply data to height' method are <1.1 cm

@ZongyangLi
Copy link

@dlebauer
Last time when we discuss about unit for canopy height: terraref/computing-pipeline#210 (comment), I was thinking I should update the output unit to meter instead of cm.

So it's in meter now. I am sorry if I make another misunderstanding, should I change it for cm?

@dlebauer
Copy link
Member

dlebauer commented Dec 4, 2017

@ZongyangLi yes please change to cm. It is an option to change all of our height data to m but at this point we have been using cm, so the easiest at this point is to use cm, and if it would be useful to change to m, we can decide with inputs from other groups and convert all of the code and data at once.

@dlebauer
Copy link
Member

dlebauer commented Dec 5, 2017

I purged the old data

delete from traits 
   where method_id = (select id from methods where name = 'Scanner 3d ply data to height')

@craig-willis
Copy link
Author

craig-willis commented Dec 5, 2017

@dlebauer Should we run another sample data run, or start full processing? (i.e., can we check more boxes in the completion criteria)

@dlebauer
Copy link
Member

dlebauer commented Dec 5, 2017

@craig-willis yes, everything else looked good; should be ready to go

@craig-willis
Copy link
Author

@ZongyangLi
There seems to be a problem handling season 1 data

        # MAC Field Scanner Season 4 Range 5 Column 6 W
        plot_record = [int(s) for s in site.split() if s.isdigit()]
        if len(plot_record) != 3:
            return 0, 0, 0, 0, 0, 0

This won't work for the season 1 names, for example: "MAC Field Scanner Season 1 Field Plot 542 E".

Is this expected?

@ZongyangLi
Copy link

@craig-willis
That's right, it can not handle that. I didn't imagine that plot records have different style, shouldn't it be constant?

Besides, I don't think season 1 or season 2 should be involve, since the laser scanner's offset(#44) is not constant during the first 2 season.

@craig-willis
Copy link
Author

@dlebauer
Is there any additional action to be taken here? I assume you want me to move forward with running on the later seasons.

@dlebauer
Copy link
Member

dlebauer commented Dec 6, 2017 via email

@ZongyangLi
Copy link

I could not find such an issue now, maybe I am wrong. The offset is constant in all the field seasons. Then we don't need to worry about this.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

3 participants