Tìm lát Sentinel cho tọa độ Long / Lat cụ thể

Tôi đang xây dựng một công cụ tải hình ảnh từ "Sentinel-2 trên AWS"

Tôi gặp khó khăn khi dịch tọa độ Long / Lat sang ô Sentinel cụ thể. Tôi biết rằng Sentinel 2 sử dụng MGRS và tôi đã cố gắng chuyển đổi Long / Lat thành MGRS và với độ chính xác thấp, tôi nhận được ID gạch S2A. Và nó hoạt động, nhưng không phải lúc nào.

Ví dụ: Long / Lat 34.665,31.625 phân giải thành "36R XA" trong MGRS, nhưng ô Sentinel 36RXA không tồn tại.

Điều gì sẽ là cách đúng để xác định ô Sentinel bằng tọa độ Long / Lat?

Phân vùng S2 trong MGRS đặc biệt khó khăn khi dữ liệu nằm giữa các vùng UTM và các ô xếp chồng lên nhau, xem hình bên dưới. Các ô từ trái sang phải là 20_M_QE, 20_M_RE, 21_M_SV và 21_M_TV:

Một tìm kiếm được thực hiện thông qua openearch cho khu vực trả về một lần mua S2 bị hỏng trong hai ô, 20_M_RE và 21_M_TV:

Vì vậy, nếu bạn quan tâm đến dữ liệu trong vùng gạch 21_M_SV được đại diện, bạn nên sử dụng 20_M_RE. Lưu ý rằng biểu diễn được sử dụng trong các ví dụ (shapefile thu được từ NGA) sự trùng lặp giữa các ô không được biểu diễn.

Cách an toàn nhất để có được các ô dường như là tìm kiếm AOI của bạn bằng cách sử dụng opensearch và sau đó lấy các ô được trả về hoặc sử dụng công cụ sentinel2-sử dụng như bugmenot123 đã đề cập.

Nếu bạn chọn ánh xạ Long / Lat thành các ô, tôi khuyên bạn nên sử dụng KML do ESA cung cấp làm định nghĩa lưới ốp lát:

https://sentinel.esa.int/web/sentinel/missions/sentinel-2/data-products

Một công cụ như vậy đã tồn tại. Nó được gọi là Sentinelsat và nguồn có sẵn trên GitHub . Nó cung cấp giao diện dòng lệnh và API Python. Nó hoạt động với Sentinel 1 và 2. Truy vấn không gian dựa trên đa giác chứ không phải điểm, nhưng nếu không thì đây chính xác là những gì bạn cần.

EDIT: 1) bạn có thể trả về ID sản phẩm (hoặc danh sách ID sản phẩm) bằng chức năng truy vấn. Và từ dòng lệnh, bạn có thể bỏ qua tùy chọn -d để không khởi chạy tải xuống (và sử dụng tùy chọn -f để tạo dấu chân của mỗi hình ảnh) 2) Liên quan đến Geojson, bạn có thể thấy bên dưới có thể là một wkt sử dụng trực tiếp thay thế. Như được đề cập bởi @T BreathumaBob, bạn có thể tạo một đa giác nhỏ dựa trên tọa độ lat / long của mình. Nhưng trên thực tế, tôi đã thử nghiệm nó với một điểm GeoJason và nó cũng hoạt động tốt. Vì vậy, không cần phải chuyển đổi điểm của bạn thành đa giác.

# search by polygon, time, and Hub query keywords
footprint = geojson_to_wkt(read_geojson('map.geojson'))
products = api.query(footprint,
                     date = ('20151219', date(2015, 12, 29)),
                     platformname = 'Sentinel-2',
                     cloudcoverpercentage = (0, 30))

mà có thể được điều chỉnh thành

# search by point, time, and Hub query keywords

products = api.query('POINT({0} {1})'.format(longitude,latitude),
                     date = ('20151219', date(2015, 12, 29)),
                     platformname = 'Sentinel-2',
                     cloudcoverpercentage = (0, 30))