Skip to content
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

no data found for atl08 in ROI #436

Closed
linxiong100 opened this issue Oct 3, 2024 · 2 comments
Closed

no data found for atl08 in ROI #436

linxiong100 opened this issue Oct 3, 2024 · 2 comments

Comments

@linxiong100
Copy link

I got no error and no data for atl08. But there should be 55 matching granules in this region.

I am running the following testing code in Amazon forest:

from sliderule import icesat2
icesat2.init("slideruleearth.io")
grand_mesa = [ {"lon": -60, "lat": -3},
               {"lon": -59.8, "lat": -3},
               {"lon": -59.8, "lat": -2.8},
               {"lon": -60, "lat": -2.8},
               {"lon": -60, "lat": -3} ]
parms = { "poly": grand_mesa}
icesat2.atl08p(parms)
@jpswinski
Copy link
Member

@linxiong100 the parms currently default to land ice for the icesat2 module. I realize this is not ideal, and in the next release we are changing it so that the default is dynamic (meaning the server side code chooses the best surface type for you).

But until then, you will need to explicitly provide the surface reference type in your parameters by adding "srt":icesat2.SRT_DYNAMIC.

So your code should look like:

from sliderule import icesat2
icesat2.init("slideruleearth.io")
grand_mesa = [ {"lon": -60, "lat": -3},
               {"lon": -59.8, "lat": -3},
               {"lon": -59.8, "lat": -2.8},
               {"lon": -60, "lat": -2.8},
               {"lon": -60, "lat": -3} ]
parms = { "poly": grand_mesa, "srt": icesat2.SRT_DYNAMIC}
icesat2.atl08p(parms)

When I run the code above, I get the following:

                                                                canopy_h_metrics  ph_count  ...   h_te_median                    geometry
time                                                                                        ...                                          
2018-11-02 18:38:42.853801984  (2.0147335490085164e+34, 9.183269356106256e-41...       909  ...  4.674732e-42  POINT (-59.80225 -2.80870)
2018-11-02 18:38:42.856002048  (0.0, 9.183269356106256e-41, -1.34316132028482...      1043  ...  1.163078e-43  POINT (-59.80225 -2.80856)
2018-11-02 18:38:42.858902016  (0.0, 9.183269356106256e-41, -1.34316132028482...       967  ...  1.541428e-43  POINT (-59.80229 -2.80838)
2018-11-02 18:38:42.862001920  (0.0, 0.0, -1.3431613202848254e-17, -3.1765583...      1019  ...  4.674732e-42  POINT (-59.80230 -2.80818)
2018-11-02 18:38:42.864602112  (0.0, 9.183269356106256e-41, -1.34316132028482...      1067  ...  1.653532e-43  POINT (-59.80233 -2.80801)
...                                                                          ...       ...  ...           ...                         ...
2024-07-22 14:52:59.217272576  (0.0, 0.0, 6.228695307664391e+31, -2.352258838...       148  ...  4.971527e-41  POINT (-59.86527 -2.80091)
2024-07-22 14:52:59.220172544  (0.0, 0.0, 6.228695307664391e+31, -2.352258838...       137  ...  4.972508e-41  POINT (-59.86527 -2.80072)
2024-07-22 14:52:59.222922624  (0.0, 0.0, 6.228695307664391e+31, -2.352258838...       124  ...  4.973489e-41  POINT (-59.86530 -2.80055)
2024-07-22 14:52:59.226172672  (0.0, 0.0, 6.228695307664391e+31, -2.352258838...       129  ...  4.975310e-41  POINT (-59.86533 -2.80034)
2024-07-22 14:52:59.228572672  (0.0, 0.0, 6.228695307664391e+31, -2.352258838...       139  ...  4.976291e-41  POINT (-59.86534 -2.80019)

[72261 rows x 20 columns]

@linxiong100
Copy link
Author

oh, got it, need to set the surface type! thank you!

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

2 participants