Nội suy dữ liệu bằng cấp sử dụng QGIS?

Tôi muốn tạo bản đồ gió. Có một phần giới thiệu cho ArcGIS và cho QGIS . Vì tôi không có quyền truy cập vào Nhà phân tích không gian, tôi đã thử phiên bản QGIS. Tuy nhiên, tôi không biết cách nội suy dữ liệu độ.

Ví dụ, hướng gió có thể là 359 ° tại một điểm và 1 ° tại điểm thứ hai. Các giá trị giữa điểm này có lẽ phải là 359 °, 0 ° hoặc 1 °.

Vì phép nội suy không biết rằng đó là dữ liệu độ, có vẻ như các giá trị nội suy nằm trong khoảng 180 °.

Làm cách nào tôi có thể nội suy dữ liệu độ bằng cách sử dụng QGIS?

Nội suy sin và cosin của góc, sau đó chuyển đổi trở lại một góc với hàm atan. Các chức năng này có sẵn trong công cụ biểu thức của QGIS. Có một atan2(dy,dx)chức năng giống như chức năng trong RI sử dụng bên dưới ...

Đây là một hàm R để minh họa. Tôi đã sử dụng meanở đây để đưa ra phép nội suy:

 dinterp = function(d){
             r=d*pi/180
             (180/pi)*atan2(mean(sin(r)),mean(cos(r)))
           }

cung cấp cho:

> dinterp(c(10,20))
[1] 15                *as expected*
> dinterp(c(340,20))
[1] 3.384676e-15      *approximately zero*
> dinterp(c(0,180))   
[1] 90                *could be +/- 90
> dinterp(c(359,180))
[1] -90.5             *small adjustment makes it go to ~-90

Trong QGIS, giả sử bạn đang tạo một raster được nội suy từ một tập dữ liệu điểm, quy trình làm việc của bạn sẽ như sau:

• Tạo hai thuộc tính mới tại mỗi điểm sinwindvà coswind, từ hướng gió của bạn theo độ. Hãy chắc chắn rằng bạn chuyển đổi sang radian nếu cần thiết.
• Tạo hai raster nội suy, một từ sinwindvà một từ coswind.
• Sử dụng máy tính raster để chuyển đổi trở lại hướng bằng cách tính toán atan2các raster hình sin và cosin nội suy.
• Chuyển đổi sang độ nếu cần thiết.

Bạn có thể tạo một quy trình xử lý cho việc này - cũng có thể đã có một plugin. Lưu ý rằng mặc dù điều này cung cấp cho bạn một phép nội suy theo hướng gió, nhưng mọi đánh giá về độ không đảm bảo của phép nội suy đó là một chút khó khăn và đòi hỏi một số phương pháp thống kê vòng tròn.