최근 수심 데이터를 다운로드 할 때 NGDC는 사용자 지정 지역 다운로드를 제공하지 않지만 글로벌 데이터를 다운로드 할 수있어 여전히 사용이 불편하다는 것을 알았습니다. 제 동료가 저에게 블로그를 보냈습니다. 그는 독일 Helmholtz Kiel 해양 연구 센터의 해양 지구과학 박사 과정 학생 인 Ren Yu 라는 저명한 저자 입니다. 다음 내용은 그의 기사 https : //www.zhihu 에서 재현 한 것입니다. com / column / c_184222208
이 기사의 두 번째 부분 인 netcdf 및 arc ASCII 코드를 변환하는 방법은 매우 실용적입니다.
대부분의 연구자들에게 개방형 데이터베이스 사용은 매우 중요한 기술입니다. 이 기사는 내가 연구 작업에 사용하는 몇 가지 뛰어난 오픈 액세스 지오 데이터베이스를 나열합니다.
- 해저 지형 기복 데이터
1.1 글로벌 고정밀 수심 데이터
JAMSTEC 기억 해 확인 "수심 측량 (MBES)를"
1.2 지역 고정밀 수심 데이터
AusSeabed Bathymetry 및 후방 산란 데이터 , 호주
1.3 15 arc-second 글로벌 지형 구호 데이터
GEBCO 2019 Gridded Bathymetry 데이터
2. 지진 데이터
2.1 지진 카탈로그
2.2 초점 메커니즘 솔루션
National Institute of Geophysics and Volcanology-Time Domain Moment Tensor Catalog , 이탈리아
BATS CMT 카탈로그 검색 , 대만
2.3 지구 모델
CRUST 1.0 : 1x1도에서 글로벌 지각 모델
3. 중력
4. 지자기
5. 지열
6. 화산
수심-지형 데이터 형식 Arc ASCII와 GMT netCDF 사이의 변환, https://zhuanlan.zhihu.com/p/138883556
AcrGIS는 지리 정보 시스템 공간 분석에 고유 한 장점이 있으므로 현재 주제를 연구하는 주요 도구가되었습니다. 그러나 통신 과정에서 데이터 처리 및 그리기를위한 주요 도구가 다르기 때문에 데이터 형식 변환 문제가 자주 발생합니다. 이 기사에서는 GMT 4 명령을 사용하여 지형 데이터 형식 Arc ASCII와 GMT가 인식하는 netCDF 간의 변환을 실현하는 방법을 소개합니다.
- Arc ASCII를 GMT netCDF 형식으로 변환
필자의 칼럼 기사 Ocean Earth Science Open Database 에는 유용한 데이터베이스가 나열되어 있습니다. 아래 그림과 같이 먼저 GMRT에서 Arc ASCII 형식의 고정밀 다중 빔 파일 MAR.asc를 선택하여 다운로드합니다.
그런 다음 xyz2grd 명령을 사용하여 Arc ASCII 형식의 데이터를 GMT가 인식하는 netCDF 형식으로 직접 변환합니다 .
xyz2grd MAR.asc -GMAR.grd -E -V
JAMSTEC과 같은 데이터베이스에서 다운로드 한 ASCII 파일 인 경우 인접 보간법을 사용하는 데 익숙합니다.
blockmean MAR.dat -R-45.6/-44.2/25.3/27 -I0.001 > MAR.xyz
nearneighbor MAR.xyz -V -R -I -S0.002 -N2 -GMAR.grd
출력 MAR.grd 파일은 GMT netCDF 형식입니다.
- GMT netCDF를 Arc ASCII 형식으로 변환
GMT netCDF 형식 데이터를 Arc ASCII 형식으로 변환 할 때는 먼저 netCDF 파일을 다시 샘플링하여 경도 및 위도 방향의 샘플링 속도가 정확히 동일하도록 한 다음 grd2xyz 명령을 사용하여 Arc ASCII 형식으로 변환해야합니다.
grdsample MAR.grd -I0.001= -GMAR.nc
grd2xyz MAR.nc -E > MAR_G.asc
여기서는 리샘플링 명령의 -I 옵션 을 "=" 와 함께 추가 하여 출력 MAR.nc 파일이 경도 및 위도 방향에서 동일한 샘플링 속도를 갖도록해야하며, 이는 데이터 형식을 변환하는 데 사용할 수 있습니다. grd2xyz 명령의.
최종 출력 MAR_G.asc 파일은 Arc ASCII 형식입니다. ArcMAP을 사용하여 언덕 그늘을 계산하면 아래 그림에서 고정밀 해저 지형을 볼 수 있습니다.
- 밟기 쉬운 구덩이 요약
grd2xyz 및 xyz2grd 명령이 Arc ASCII 형식 데이터를 처리 할 때 "-E"옵션을 추가해야합니다!
b.GMT netCDF를 Arc ASCII로 변환 할 때 첫 번째 단계에서 리샘플링 명령의 "-I"옵션 을 "="!
최종 진술 :이 기사에 포함 된 GMT 코드는 모두 GMT 4 버전이며 더 높은 GMT 버전에는 적용되지 않습니다!