일반 위도/경도 그리드에 대한 극좌표 입체 투영

Question

이것은 예제와 데이터로 훨씬 더 구체적이라는 점에서 이전 질문과 같습니다. 나는 netcdf가 here을 얻을 수있는 극지방의 투영법에서의 304 x 448 배열의 얼음 데이터를 가지고있다. 각 격자 상자의 해당 위도/경도 점은 file에서 찾을 수 있습니다. 나는이 데이터를 각 그리드 박스가 학위 인 정규 lat/lon 그리드 (180x360)로 변환하고 싶다. 이 일을하는 파이썬 방법이 있습니까? Basemap을 사용하는 방법 또는 CDO (기후 데이터 연산자)를 사용하는 방법?

미리 감사드립니다. 제임스

Answer 1

먼저,

ncks -A -v lat,lon NSIDC_LatLon.nc sample.nc 

는 그런 다음 예와 그 재격 자화하는 NCO의 ncremap을 사용할 수 있습니다, 예를 들어와 데이터 파일에 그리드를 추가하는 경우

ncremap -i sample.nc -d 1x1.nc -o out.nc 

불행하게도이 의지 입력 파일에 NCO를 처리하기 위해 방금 패치 한 비정상적인 순서로 차원이 포함되어 있으므로 NCO의 현재 스냅 샷이 필요합니다 (릴리스 된 버전이 아님). regridded 출력의 샘플은 here입니다.

OP 요청에 따라 파일을 재생성 한 다음 업로드했습니다 (here). 나는 1x1 크기의 데이터 집합 (ncremap -g 옵션)을 사용 설명서에 표시된 내가 처음 그와 함께 regridded, 1x1 크기의 그리드 파일을 만들었 그래서 주위에 앉아 있지 않았다 답장을

ncremap -i ~/NSIDC.nc -g ${DATA}/grids/180x360_SCRIP.20150901.nc -o ~/NSIDC_1x1.nc