구형 지구 (파이썬 /베이스 맵)에서 2D 벡터 필드의 발산

Question

나는 여러 지점에서 바람의 u (zonal) 및 v (자오선) 구성 요소의 2D 배열을 가지며 위치는 해당 2D 위도 및 경도 배열.

이 바람의 발산을 계산하고 싶습니다. divergence는 du/dx + du/dy라고 이해하지만, x와 y 좌표는 위도와 경도이므로 올바르게 출력되지 않습니다. 예를 들어 지구상의 균일 한 자오선 바람 (v = 상시 모든 곳)은 모든 색인에서 v 배열이 일정한 숫자 임에도 불구하고 투표에서 분기/수렴을 실제로 가져옵니다.

필자는 Basemap 인스턴스에 분기를 플로팅하는 것을 궁극적으로 생각하고 있습니다.

어떻게해야합니까?

Answer 1

발산 수식은 사용자가 발견 한 좌표계에 따라 다릅니다. 수식을 수정하는 것은 가치있는 연습이지만 (좌표 수식과 체인 규칙의 변경 사용) wikipedia에는 필요한 수식의 유용한 테이블이 있습니다.

구면 좌표계의 공식을 원한다면 구체의 표면에서 작업 중이므로 을 잡고 싶을 것입니다. 이것은 첫 번째 용어를 삭제하고 다른 두 용어를 확장하는 효과가 있습니다.