Tôi cần tập hợp dữ liệu từ NetCDF bằng cách sử dụng một shapefile cụ thể. Dữ liệu là nhiệt độ mặt nước biển và màu đại dương ở độ phân giải 1/4 độ. Tôi có 4 đa giác mô tả về Hoa Kỳ. Thềm lục địa Đông Bắc thềm hệ sinh thái biển lớn và đó là các thành phần phụ mà tôi cần sử dụng để trích xuất dữ liệu. Tôi đang làm việc với các tệp tổng hợp hàng tháng từ 1982-2014, vì vậy thói quen trích xuất dữ liệu này cần được tự động hóa. Các tệp đã được đặt vào lưới khu vực làm việc gần đúng của [35, 45, -80, -60].
Trước đây, chúng tôi đã chuyển đổi các tệp dữ liệu HDF5 sang các trình quét trong R và xử lý chúng theo cách này, nhưng phương pháp này thực sự không hiệu quả và tôi chắc chắn có một giải pháp tốt hơn trong python bằng các tệp NetCDF hiện tại.
Cho đến nay tôi đã sử dụng GDAL và Fiona để đọc các tệp hình dạng và NetCDF4 để tải các tệp dữ liệu. Tôi không chắc chắn làm thế nào để đi về tập hợp dữ liệu. Tôi đã tìm thấy cái này:
GDAL cho Python: trích xuất tên miền phụ từ tệp NetCDF?
Nhưng tôi không có ý tưởng rõ ràng nhất về cách đặt tập tin netcdf bằng cách sử dụng bất cứ thứ gì ngoài hộp giới hạn đơn giản, điều mà các đa giác này chắc chắn là không.
Điểm trong các thói quen đa giác có thể sẽ mất nhiều thời gian để hoạt động, nhưng có lẽ tôi có thể tập hợp dữ liệu bằng cách sử dụng hộp giới hạn nhỏ hơn được xoay để phù hợp với các hình dạng này như một điểm bắt đầu bẩm sinh và sau đó thực hiện tìm kiếm theo điểm:
Đặt lại tệp netCDF curvilinear (đầu ra mô hình ROMS) bằng cách sử dụng hộp giới hạn lon / lat.
Có ý kiến gì không?
CHỈNH SỬA 1 :
Tôi vừa xem qua gói OpenClimateGIS có vẻ như nó hoàn toàn phù hợp với hóa đơn ... Tôi sẽ đi cùng với điều này để xem liệu tôi có thể làm cho nó hoạt động không: http://ncpp.github.io/ocgis/examples. html # nâng cao