이 netCDF 래스터는 어떤 좌표계입니까?


10

netCDF 래스터 파일을 얻었지만 래스터가 작성된 좌표계의 이름을 가져올 메타 데이터를 얻을 수 없습니다. 래스터 자체에는 좌표계가 포함되어 있지 않습니다. 나는 그것이 WGS84 일 것이라고 생각했고, 처음 보았을 때처럼 보였지만 ArcMap에서 추가 조사를 통해 그것이 일반적인 시스템이 아니라는 것을 알았습니다. 표시 방법은 다음과 같습니다. 여기에 이미지 설명을 입력하십시오

주황색 래스터는 WGS84의 일반 래스터로 비교 목적으로 여기에 삽입했습니다. 자주색은 좌표계가 알려지지 않은 래스터입니다. 이것이 무엇인지 힌트가 있습니까?

일부 업데이트 : 여기에 netCDF의 래스터는 다음과 같습니다 https://www.dropbox.com/s/nottbl9yt6dwss6/sic_average_nclimate.nc?dl=0 나는 또한 이미지 제공 업체에서 일부 메타 데이터를 얻을 수 있었다 :

netcdf sic_average_nclimate {
dimensions:
    nlon = 361 ;
    nlat = 90 ;
    nseas = 4 ;
variables:
    float SIC_Change(nlat, nlon) ;
        SIC_Change:Title = "Gridded Multi-Model Ensemble Mean Annual Mean Change in Ice Concentration 21C-20C" ;
    float SIC_Season_Change(nseas, nlat, nlon) ;
        SIC_Season_Change:Title = "Gridded Multi-Model Ensemble Mean Seasonal Mean Change in Ice Concentration 21C-20C" ;
    float SIC_Change_STD(nlat, nlon) ;
        SIC_Change_STD:Title = "Gridded Multi-Model Standard Deviation of the Annual Mean Change in Ice Concentration 21C-20C" ;
    float SIC_Season_Change_STD(nseas, nlat, nlon) ;
        SIC_Season_Change_STD:Title = "Gridded Multi-Model Standard Deviation of the Seasonal Mean Change in Ice Concentration 21C-20C" ;
    float LAT(nlat) ;
        LAT:Title = "Latitude" ;
    float LON(nlon) ;
        LON:Title = "Longitude" ;

// global attributes:
        :Title = "Ice Concentration metrics for Model subset as in Figure 1 of NCLIMATE paper" ;

위도는 경계를 길게 나타내지 만 좌표계에 관한 정보는 없습니다.


1
남극 대륙은 보통 극지 입체 투영법에 있으므로 모든 것이 남극에서 북쪽으로 가고 있습니다. nsidc.org/data/polar_stereo/ps_grids.html
Mapperz

1
이상하게도 0-360처럼 보이지만 왜 북반구에 그려지고 있습니까?
mkennedy

1
NetCDF 헤더에보고 된 익스텐트가 위도 / 경도처럼 보이지만 위도 : 0 ~?, 론 : 0 ~ 360 인 경우 : -90 ~ (?-90) 및 론 : -180 ~ + 180?
mkennedy 2014

@Mapperz & mkennedy에서 Polar Stereographic Projection을 지정하려고 시도했지만 작동하지 않았습니다. 수직으로 아래쪽으로 미러링 할 수는 있지만 오른쪽 절반은 여전히 ​​180E 경도를 초과합니다. 메모장으로 netCDF를 열었지만 위도, 경도와 관련된 정보를 찾을 수 없습니다. 당신이 당신쪽에 그것을 열어 궁금하다면 질문에 래스터를 첨부하고 있습니다.
multigoodverse

해당 링크가 작동하지 않습니다. dropbox에 "소유자가 액세스 권한을 부여하지 않았습니다"라고 표시되어 있습니다.
BradHards 2014

답변:


6

다음 GDAL 명령을 사용하여 데이터를 올바른 위치로 이동할 수 있습니다.

gdal_translate -a_srs EPSG:4326 -a_ullr 0 0 360 -90 NETCDF:"sic_average_nclimate.nc":SIC_Change change360.tif
gdalwarp -t_srs WGS84 change360.tif change180.tif  -wo SOURCE_EXTRA=1000 --config CENTER_LONG 0

(Frank Warmerdam의 약간의 도움으로 180 자오선을 잘라 0360 에서 -180180으로 래스터를 재 투영하는 방법 )

입체 투영에서 결과는 다음과 같습니다.

여기에 이미지 설명을 입력하십시오


1

(이것은 답변이 아니라 의견이 되기에는 너무 긴 프로세스 제안입니다)

두 래스터에서 공통 점을 식별 할 수있는 것처럼 GeoReferencer 도구를 사용하여 QGIS에서이를 지리 참조 할 수 있습니다. 이렇게하면 GCP 포인트 세트와 변환 매개 변수가있는 GDAL 스크립트가 생성됩니다. 두 투영의 차이점이 무엇인지 설명 할 수 있으므로이 스크립트를이 포인트 샘플로 질문에 추가 할 수 있습니다.

최악의 시나리오 인 경우 재배치 된 래스터를 기준으로 대략적인 값을 계산하여 목표를 달성하는 데 도움이 될 수 있습니다.


+ 1, 예. 지리 참조는 마지막 글 머리 기호 일 수 있지만 제어점을 일치시키기가 어렵 기 때문에 프로세스에서 약간의 오류가 발생합니다. 또한 netCDF 래스터는 동일한 주요 자오선이 없으므로 분할해야합니다. 예를 들어 남극 대륙의 꼬리는 오른쪽에 있고 주황색 래스터는 왼쪽에 있습니다.
multigoodverse

0

왼쪽 하단 모서리는 (0,0) 좌표이지만 나에게 (0,0)은 귀하의 경우 오른쪽 상단 모서리의 좌표이며 원점이 -90 인 지리적 좌표계를 모릅니다.

ArcGIS에서 작업 할 때 지리 참조 툴바 를 사용 하여 데이터를 올바른 위치로 이동할 수 있습니다. 오른쪽 상단을 (0,0)으로, 왼쪽 하단을 (-360, -90)으로 이동하면 정확한 변환이 이루어 지지만 여전히 날짜 표시 줄에 문제가 있습니다.

여기에 이미지 설명을 입력하십시오

따라서이 링크 를 사용하면 사용자 지정 지리적 좌표계를 만드는 데 도움이됩니다. 요약하자면, 참조 시스템에서 자오선이 -180 인 사용자 정의 WGS 84를 생성 한 다음 경도 회전을 사용하여 사용자 정의 변환을 작성하십시오.

여기에 이미지 설명을 입력하십시오

이 새로운 시스템에서 지리 참조는 원점을 (-180, -90)으로 설정하는 것으로 구성됩니다. 번역이 필요하기 때문에 하나의 제어점으로 충분합니다. 왼쪽 하단 픽셀을 신중하게 선택하고 마우스 오른쪽 버튼을 클릭 한 후 정확한 XY 값을 추가하십시오. 그러면 다음과 같아야합니다

여기에 이미지 설명을 입력하십시오

마지막으로, 좌표계가 WGS 84를 기반으로한다는 확신은 없습니다. 그러나 지리적 좌표계가 확실하고 래스터의 해상도가 잘못된 데이텀 선택으로 인한 잠재적 오류보다 거칠다는 것은 확실합니다.

당사 사이트를 사용함과 동시에 당사의 쿠키 정책개인정보 보호정책을 읽고 이해하였음을 인정하는 것으로 간주합니다.
Licensed under cc by-sa 3.0 with attribution required.