img2grdgmt - 클라우드의 온라인

프로그램:

이름

img2grd - Mercator 또는 지리적 형식으로 img 파일의 하위 집합 추출

개요

img2grd 이미지 파일 grd 파일 지방 유형 [ ] [ [최소라트/최대라트] ] [ ] [ 분 ] [ ] [ 탐색
] [ [규모] ] [ [수평] ] [ 맥론 ] [ -n]

참고 : 옵션 플래그와 관련 인수 사이에는 공백이 허용되지 않습니다.

기술

img2grd img 형식 파일을 읽고 하위 집합을 추출한 다음 그리드 파일에 씁니다. 그만큼 -M
옵션은 img 파일의 Spherical Mercator 투영이
보존되거나 메르카토르 투영을 실행 취소하여 지리적 그리드를 작성해야 하는 경우. 만약에
지리적 그리드를 선택하면 정확한 리샘플링을 요청할 수도 있습니다. -R 주어진.

필요한 인수

이미지 파일
해양 중력 또는 해저 지형 필드와 같은 Mercator img 형식 파일
Sandwell과 Smith의 위성 고도계 데이터에서 추정. 사용자가 설정한 경우
환경 변수 $GMT_DATADIR다음, img2grd 찾기 위해 노력할 것입니다 이미지 파일 in
$GMT_DATADIR; 그렇지 않으면 열려고 할 것입니다 이미지 파일 직접.

-Ggrd 파일
grd 파일 출력 그리드 파일의 이름입니다.

-NS[단위]서쪽/동쪽/남쪽/북쪽[/지민/지맥스][NS]
서쪽, 동쪽, 남쪽및 북쪽 관심 영역을 지정하고 다음을 지정할 수 있습니다.
십진법 또는 [+-]dd:mm[:ss.xxx][W|E|S|N] 형식으로 표시됩니다. 추가 r 낮은 경우
w/e/s/n 대신 왼쪽 및 오른쪽 위 지도 좌표가 제공됩니다. 둘
속기 -RG 과 -RD 전역 도메인(경도 0/360 및 -180/+180)을 나타냅니다.
각각 위도에서 -90/+90). 또는 그리드 생성을 위해 다음을 제공하십시오.
R코드론/위도/nx/ny어디로 암호 L, C, R의 2자 조합(왼쪽,
중앙 또는 오른쪽) 및 T, M, B는 상단, 중간 또는 하단을 나타냅니다. 예: 왼쪽 하단의 경우 BL.
이것은 직사각형 영역의 어느 지점을 나타냅니다 경도/위도 좌표 참조
및 그리드 치수 nx 과 ny 그리드 간격을 통해 -I 생성하는 데 사용됩니다
해당 지역. 또는 기존 그리드 파일의 이름을 지정합니다.
그리고 -R 설정(및 해당하는 경우 격자 간격)이 격자에서 복사됩니다.
사용 -R단위 선택한 것과 호환되는 투영된(데카르트) 좌표를 예상합니다. -J
실제 직사각형 지리적 영역을 결정하기 위해 역으로 투영합니다. 을위한
투시도(-p), 선택적으로 / 추가지민/지맥스. 투시도의 경우
(-p), z 범위(지민, 지맥스)를 추가하여 세 번째 차원을 나타낼 수 있습니다. 이것
사용할 때만 수행해야 합니다. -Jz 옵션, 만 사용하는 경우가 아님 -p 옵션을 선택합니다.
후자의 경우 평면의 투시도가 그려지며 세 번째는 없습니다.
치수.

선택 사항 인수

-C 투영 중심을 기준으로 x 및 y 메르카토르 좌표를 설정합니다. [기본값은
그리드의 왼쪽 하단 모서리를 기준으로]. 필요하다 -M.

-NS[최소라트/최대라트]
확장된 위도 범위 -80.738/+80.738을 사용합니다. 또는 추가
최소라트/최대라트 입력 img 파일의 위도 범위로. [기본값은
-72.006/72.006]. 범위를 결정할 수 있으므로 일반적으로 필요하지 않습니다.
파일 크기 검사.

-E 다음과 같은 경우에 사용할 수 있습니다. -M 최종 그리드가 정확히 동일하도록 설정되지 않았습니다.
요청한 지역 -R. 기본적으로 최종 영역은
원래 메르카토르 지역이며 일반적으로
요청된 위도 범위, 또한 위도의 그리드 증분은
경도 증분을 일치시킵니다. 그러나 추가 리샘플링은 작은
보간 오류이며 출력 그리드가 다음과 일치해야 하는 경우에만 사용해야 합니다.
요청된 영역이고 x_inc = y_inc입니다. 이 경우에 의해 설정된 영역 -R 해야합니다
증분의 배수로 제공됩니다(.eg, -R0/45/45/72).

-I 가리키다 분 경도의 분 단위로 입력 이미지 픽셀의 너비로.
[기본값은 2.0입니다]. 에서 픽셀 크기를 결정할 수 있으므로 일반적으로 필요하지 않습니다.
크기 검사.

-M 구형 메르카토르 격자선을 출력합니다[기본값은 지리적 경도/위도 격자입니다]. 그만큼
img 파일의 구형 메르카토르 투영이 보존되므로 영역이 -R
사용자가 설정한 내용이 약간 수정되었습니다. 수정된 영역은 가장자리에 해당합니다.
픽셀 [또는 그룹 탐색 픽셀]. 그리드 파일 헤더는 x와
y축 길이는 이미지의 서쪽 및 남쪽 가장자리로부터의 거리를 나타냅니다.
사용자 기본 단위로 측정 -Jm1 및 조정 -R. 설정하여
디폴트 값 프로젝트_ 타원체 = Sphere, 사용자는 조정된 크기로 오버레이를 만들 수 있습니다. -R
일치하도록 합니다. 보다 사용 예 아래에. 조정된 -R 그리드에도 쓰여있다
나중에 찾을 수 있도록 헤더 설명입니다. 보다 -C 에 상대적인 좌표 설정
프로젝션 센터.

-N탐색 입력 img 픽셀의 값을 평균으로 탐색 by 탐색 사각형, 만들기
각 정사각형에 대해 하나의 출력 픽셀. 함께 사용하는 경우 -T3 평균을보고합니다
0과 1 사이의 제약. 함께 사용하는 경우 -T2 출력은 평균 데이터 값이 됩니다.
또는 평균 제약 조건이 > 0.5인지 여부에 따라 NaN입니다. 탐색 균등하게 나누어야 한다
imgfile의 크기를 픽셀 단위로 변환합니다. [기본 1 평균화하지 않음].

-NS[규모]
img 파일 값에 다음을 곱하십시오. 규모 그리드 파일에 저장하기 전에. [기본값은
1.0]. 최근 img 파일의 경우: img topo 파일은 (수정된) 미터에 저장됩니다. [-S1];
mGal*10의 자유 공기 중력 파일 [-SmGal을 얻기 위해 0.1]; 수직 처짐 파일
마이크로라디안*10 [-S마이크로 라디안을 얻기 위해 0.1], 수직 중력 기울기 파일
엇보스*50 [-SEotvos를 얻으려면 0.02, 또는 -SmGal/km를 얻으려면 0.002]). 아니오인 경우 규모 주어진다
단서를 찾기 위해 파일 이름을 조사하여 규모를 결정하려고 합니다.

-T유형 유형 제약 정보의 인코딩을 처리합니다. 유형 = 0은 없음을 나타냅니다.
이러한 정보는 img 파일에 인코딩됩니다(1995년 이전 버전의
중력 데이터) 모든 데이터를 가져옵니다. 유형 > 0은 제약 정보가
인코딩(img 파일의 1995년 이후(현재) 버전)
다음과 같이 그리드 파일을 생성합니다. -T1 모든 지점에서 데이터 값을 얻습니다. -T2 데이터를 가져옵니다
제한된 지점의 값 및 보간된 지점의 NaN; -T3 에서 1을 얻는다
제한된 점과 보간된 점에서 0 [기본값은 1].

-V[수평] (더 ~)
상세 수준 [c]을 선택합니다. 보는 데 도움이되므로 특히 여기를 권장합니다.
좌표가 조정되는 방법.

-W맥론
가리키다 맥론 입력 img 파일의 최대 경도 범위로. 버전
1995년부터 맥론 = 360.0, 일부 이전 파일은 맥론 = 390.0.
[기본값은 360.0입니다].

-^ or 다만 -
명령 구문에 대한 짧은 메시지를 인쇄한 다음 종료합니다(참고: Windows
그냥 사용 -).

-+ or 다만 +
설명을 포함하여 광범위한 사용(도움말) 메시지를 인쇄하십시오.
모듈별 옵션(GMT 공통 옵션 아님)을 선택한 다음 종료됩니다.

-? or 아니 인수
옵션 설명을 포함하여 전체 사용법(도움말) 메시지를 인쇄한 다음
출구.

--번역
GMT 버전을 인쇄하고 종료합니다.

--show-datadir
GMT 공유 디렉토리의 전체 경로를 인쇄하고 종료합니다.

지리적 사용 예

XNUMXD덴탈의 -M 출력 그리드가 지리적 위치에 있어야 하는 경우 옵션을 제외해야 합니다.
좌표. 지역에서 데이터를 추출하려면 -R-40/40/-70/-30부터 world_grav.img.7.2 과
지리적 좌표를 얻기 위해 재투영하면 시도해 볼 수 있습니다.

img2grd world_grav.img.16.1 -Gmerc_grav.nc -R-40/40/-70/-30 -V

img 파일의 위도 간격은 메르카토르 단위로 등거리이므로
결과 그리드는 지정된 그리드와 일치하지 않습니다. -R 정확하고 위도 간격은
경도 간격과 같습니다. 정확한 일치가 필요한 경우 -R 그리고 같은 간격
경도와 위도에서 -E 옵션 :

img2grd world_grav.img.16.1 -Gmerc_grav.nc -R-40/40/-70/-30 -E -V

메르카토르 사용 예

img 파일은 메르카토르 투영에 있으므로 지리적 그리드를 추출하면 안 됩니다.
메르카토르 지도를 만들 계획이라면 당신이 그렇게했다면 당신은 투영을 끝내고
그리드를 재투영하여 단파장 세부 정보를 잃습니다. 더 나은 사용 -M 그리드를 플롯
원하는 메르카토르 투영법과 동일한 축척의 선형 투영법을 사용합니다(GMT 참조).
예 29). 지역에서 데이터를 추출하려면 -R-40/40/-70/-30부터 world_grav.img.7.2, 운영

gmt img2grd -M world_grav.img.7.2 -Gmerc_grav.nc -R-40/40/-70/-30 -V

참고로 -V 옵션은 범위가 다음으로 조정되었음을 알려줍니다.
-R-40/40/-70.0004681551/-29.9945810754. 우리는 또한 사용할 수 있습니다 grdinfo 그리드 파일을 찾으려면
헤더는 해당 지역이 -R0/80/0/67.9666667 이것은 우리가 얻을 x,y의 범위입니다.
다음을 사용하여 구형 메르카토르 투영에서 -R-40/40/-70.0004681551/-29.9945810754 및
-Jm1. 따라서 ship.lonlatgrav를 가져와 merc_grav.nc를 샘플링하는 데 사용하려면 다음과 같이 할 수 있습니다.

GMT 세트 PROJ_ELLIPSOID 구

gmt mapproject -R-40/40/-70.0004681551/-29.9945810754 -Jm1i ship.lonlatgrav | \
gmt grdtrack -Gmerc_grav.nc | GMT 맵 프로젝트 \
-R-40/40/-70.0004681551/-29.9945810754 -Jm1i -I > 배.lonlatgravsat

위와 같은 방법으로 데이터를 투영 및 투영 해제하는 것이 좋습니다.
보간 단계가 하나만 있기 때문입니다(in grdtrack). 하나라면
먼저 그리드 파일을 lon,lat로 변환한 다음 샘플을 시도합니다. 두 가지가 있습니다.
보간 단계(변환 및 샘플링).

위의 그리드에서 lon,lat 그리드를 만들려면 다음을 사용할 수 있습니다.

gmt grdproject merc_grav.nc -R-40/40/-70.0004681551/-29.9945810754 -Jm1i -I -D2m -Ggrav.nc

어떤 경우에는 이것이 쉽지 않을 것입니다. -R 두 좌표계에서 정렬되지 않을 수 있습니다.
잘. 이런 일이 발생하면 다음을 사용할 수도 있습니다(사실 항상 사용하는 것이 더 나을 수 있음).

gmt grd2xyz merc_grav.nc | GMT 맵 프로젝트 \
-R-40/40/-70.0004681551/-29.994581075 -Jm1i -I | \
GMT 표면 -R-40/40/-70/70 -I2m -Ggrav.nc

위 지역의 메르카토르 지도를 만들려면 gmt.conf 값이 다음과 같다고 가정합니다.
PROJ_LENGTH_UNIT는 인치입니다. 그러면 위의 merc_grav.nc 파일이 -Jm1i
너비는 80인치입니다. 를 사용하여 8인치 너비의 지도를 만들 수 있습니다. -Jx모든 지도 프로그램의 0.1i
이 그리드에 적용됨(예: grdcontour, 이미지, grdview), 오버레이에 대해
lon,lat에서 작업(예: 싸이, PSCoast) 위의 조정된 값을 사용할 수 있습니다. -R 과 -Jm0.1 얻을
일치하는 두 시스템.

그러나 우리는 이것보다 더 똑똑할 수 있습니다. 입력 img 파일에 픽셀 2.0이 있음을 인식
분 폭(또는 grdinfo merc_grav.nc로 nx 및 ny 확인)
merc_grav.nc는 img 파일의 전체 해상도를 사용했으며 2400 x 2039 픽셀이 있습니다.
너비가 8인치인 경우 인치당 300픽셀입니다. 우리는 그렇게 많이 필요하지 않다고 결정하고
인치당 100픽셀이면 만족할 것이므로 데이터를 3 x 3으로 평균화하려고 합니다.
사각형. (등고선 플롯을 원하는 경우 훨씬 더 많은 데이터를 평균화하도록 선택할 것입니다.
(예: 6 x 6) 부드러운 윤곽선을 얻습니다.) 2039는 3으로 나눌 수 없으므로
다른 조정 -R 이 시간:

gmt img2grd -M world_grav.img.7.2 -Gmerc_grav_2.nc -R-40/40/-70/-30 -N3 -V

이번에는 조정된 영역이 -R-40/40/-70.023256525/-29.9368261101 및
출력은 800 x 601 픽셀로 우리에게 더 좋은 크기입니다. 이제 우리는 인공물을 만들 수 있습니다
이를 위한 조명 파일 사용 그라데이션:

GMT grdgradient merc_grav_2.nc -Gillum.nc -A0/270 -Ne0.6

"grav.cpt"라는 CPT 파일도 있으면 색상 음영 기복 지도를 만들 수 있습니다.
이렇게 :

gmt grdimage merc_grav_2.nc -Iillum.nc -Cgrav.cpt -Jx0.1i -K > map.ps
gmt psbasemap -R-40/40/-70.023256525/-29.9368261101 -Jm0.1i -Ba10 -O >> map.ps

img 파일에서 위도/경도의 제한된 데이터 값만 가져오려고 한다고 가정합니다.
좌표. 그런 다음 실행 img2grd 와 더불어 -T2 옵션, 사용 grd2xyz 값을 덤프하려면 파이프
grep -v NaN을 통해 NaN을 제거하고 파이프를 통해 지도 프로젝트 역으로
위와 같이 투영합니다.

onworks.net 서비스를 사용하여 온라인에서 img2grdgmt 사용