Hệ tọa độ nào là raster netCDF này?


10

Tôi đã nhận được tệp raster netCDF, nhưng tôi không thể lấy bất kỳ siêu dữ liệu nào để lấy tên của hệ tọa độ mà raster đã được xây dựng. Bản thân raster không có hệ thống tọa độ nào được nhúng. Tôi nghĩ rằng nó sẽ chỉ là một WGS84, và nó trông giống như vậy ngay từ cái nhìn đầu tiên, nhưng với sự điều tra sâu hơn trong ArcMap, tôi thấy rằng đó là một hệ thống không phổ biến. Đây là cách nó hiển thị: nhập mô tả hình ảnh ở đây

Raster màu cam là một raster bình thường trong WGS84 mà tôi đã chèn vào đây để so sánh. Màu tím là raster với hệ tọa độ không xác định. Bạn có bất kỳ manh mối gì đây có thể là?

Một số cập nhật: Đây là raster netCDF: https://www.dropbox.com/s/nottbl9yt6dwss6/sic_aenses_nclimate.nc?dl=0 Tôi cũng có thể nhận được một số siêu dữ liệu từ nhà cung cấp hình ảnh:

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" ;

Chúng hiển thị ranh giới lâu dài, nhưng dường như không có bất kỳ thông tin nào liên quan đến hệ tọa độ.


1
Nam Cực thường nằm trong Polar Stereographic Projection nên mọi thứ đều đi về phía bắc từ cực nam. nsidc.org/data/polar_st bd / ps_grids.html
Mapperz

1
Thật kỳ lạ, trông giống như 0-360, nhưng tại sao nó lại âm mưu ở bán cầu bắc?
mkennedy 16/12/14

1
Nếu các phạm vi được báo cáo trong tiêu đề NetCDF trông giống lat / lon, nhưng là lat: 0 đến ?, Lon: 0 đến 360. Bạn có thể hack chúng thành lat: -90 đến (? - 90) và lon: -180 đến + 180?
mkennedy

@Mapperz & mkennedy, tôi đã cố gắng gán cho nó một Phép chiếu lập thể Polar, nhưng nó không hoạt động. Tôi có thể có thể phản chiếu nó theo chiều dọc xuống dưới, nhưng nửa bên phải của nó sẽ vẫn vượt quá kinh độ 180E. Tôi đã mở netCDF bằng notepad, nhưng không thể tìm thấy bất kỳ thông tin nào liên quan đến lat, lâu. Tôi sẽ đính kèm raster trong câu hỏi nếu bạn tò mò muốn mở nó về phía bạn.
đa cấp

Liên kết đó không hoạt động - dropbox cho biết "chủ sở hữu chưa có quyền truy cập".
BradHards

Câu trả lời:


6

Bạn có thể dịch chuyển dữ liệu vào đúng vị trí bằng các lệnh GDAL sau:

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

(với một chút trợ giúp của Frank Warmerdam: Cách chuyển hướng raster từ 0 360 xuống -180 180 bằng cách cắt 180 kinh tuyến )

và kết quả trông như thế này trong một phép chiếu lập thể:

nhập mô tả hình ảnh ở đây


1

(Đây không phải là một câu trả lời, mà là một đề xuất quá trình quá dài để có thể nhận xét)

Bạn có thể tham chiếu địa lý trong QGIS bằng công cụ GeoReferencer khi bạn dường như có thể xác định các điểm chung trên cả hai trình quét. Điều này sẽ tạo ra một tập lệnh GDAL với một tập hợp các điểm GCP và các tham số chuyển đổi. Bạn có thể thêm tập lệnh này với một mẫu các điểm này vào câu hỏi của bạn, vì nó có thể giải thích sự khác biệt giữa hai phép chiếu là gì.

Trường hợp xấu nhất, bạn sẽ có ước tính sơ bộ về các giá trị dựa trên raster được định vị lại có thể giúp bạn đạt được mục tiêu của mình.


+ 1 và có, hội nghị địa lý có thể là viên đạn cuối cùng, nhưng quá trình này sẽ gây ra một số lỗi vì khó khớp với các điểm kiểm soát. Ngoài ra, raster netCDF sẽ phải được phân chia vì nó không có cùng kinh tuyến gốc. Ví dụ, đuôi của Nam Cực nằm bên phải, trong khi raster màu cam hiển thị nó ở bên trái.
đa cấp

0

Góc dưới bên trái là tọa độ (0,0), nhưng đối với tôi (0,0) là tọa độ của góc trên bên phải trong trường hợp của bạn và tôi không biết bất kỳ hệ tọa độ địa lý nào có gốc tọa độ là -90.

Khi bạn làm việc trong ArcGIS, bạn có thể sử dụng thanh công cụ hội nghị địa lý để di chuyển dữ liệu của mình tại đúng vị trí. Di chuyển phía trên bên phải sang (0,0) và phía dưới bên trái sang (-360, -90) và điều này sẽ mang lại một biến đổi chính xác nhưng bạn vẫn sẽ gặp vấn đề với dòng thời gian.

nhập mô tả hình ảnh ở đây

Do đó, bạn có thể thử liên kết này sẽ giúp bạn tạo một hệ tọa độ địa lý tùy chỉnh. Để tổng hợp, trong hệ thống tham chiếu của bạn, hãy tạo WGS 84 tùy chỉnh với kinh tuyến gốc ở -180, sau đó chuyển đổi tùy chỉnh bằng cách sử dụng xoay kinh độ.

nhập mô tả hình ảnh ở đây

Hội nghị địa lý trong hệ thống mới này sẽ bao gồm cài đặt gốc tại (-180, -90). Một điểm kiểm soát là đủ vì bạn chỉ cần một bản dịch. Chọn cẩn thận pixel phía dưới bên trái, nhấp chuột phải và thêm các giá trị XY chính xác). Sau đó, nó sẽ trông như thế này

nhập mô tả hình ảnh ở đây

Như một nhận xét cuối cùng, không có gì chắc chắn rằng hệ tọa độ dựa trên WGS 84. Nhưng chắc chắn là hệ tọa độ địa lý tính theo độ và độ phân giải của raster của bạn có thể thô hơn lỗi tiềm ẩn do lựa chọn dữ liệu sai.

Khi sử dụng trang web của chúng tôi, bạn xác nhận rằng bạn đã đọc và hiểu Chính sách cookieChính sách bảo mật của chúng tôi.
Licensed under cc by-sa 3.0 with attribution required.