ระบบพิกัดใดที่ netCDF แรสเตอร์นี้

ฉันได้รับไฟล์ rCD netCDF แต่ฉันไม่สามารถรับเมทาดาทาใด ๆ เพื่อรับชื่อของระบบพิกัดที่ raster ได้ถูกสร้างขึ้น แรสเตอร์นั้นไม่มีระบบพิกัดใด ๆ ฝังอยู่ ฉันคิดว่ามันน่าจะเป็น WGS84 และดูเหมือนว่าในตอนแรก แต่การสอบสวนเพิ่มเติมใน ArcMap ฉันเห็นว่ามันเป็นระบบที่ค่อนข้างไม่ธรรมดา นี่คือวิธีแสดง:

แรสเตอร์สีส้มเป็นแรสเตอร์ทั่วไปใน WGS84 ซึ่งฉันใส่ไว้ที่นี่เพื่อวัตถุประสงค์ในการเปรียบเทียบ สีม่วงคือแรสเตอร์ที่มีระบบพิกัดที่ไม่รู้จัก คุณมีเงื่อนงำอะไรบ้างที่อาจเป็นเช่นนี้?

การอัปเดตบางอย่าง: นี่คือ raster 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" ;

พวกเขาแสดงขอบเขตละติจูดยาว แต่เห็นได้ชัดว่าไม่มีข้อมูลใด ๆ เกี่ยวกับระบบพิกัด

คุณสามารถย้ายข้อมูลไปยังตำแหน่งที่ถูกต้องโดยใช้คำสั่ง 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: วิธีการปฏิเสธ raster จาก 0 360 ถึง -180 180 ด้วยการตัด 180 meridian )

และผลลัพธ์จะมีลักษณะเช่นนี้ในการประมาณการภาพสามมิติ:

(นี่ไม่ใช่คำตอบ แต่เป็นคำแนะนำกระบวนการที่ยาวเกินกว่าจะแสดงความคิดเห็นได้)

คุณสามารถอ้างอิงทางภูมิศาสตร์ใน QGIS โดยใช้เครื่องมือ GeoReferencer เนื่องจากคุณสามารถระบุจุดร่วมใน rasters ทั้งสองได้ สิ่งนี้จะสร้างสคริปต์ GDAL พร้อมชุดคะแนน GCP และพารามิเตอร์การแปลง คุณสามารถเพิ่มสคริปต์นี้พร้อมกับตัวอย่างคะแนนเหล่านี้ในคำถามของคุณเนื่องจากอาจอธิบายความแตกต่างระหว่างการฉายสองภาพ

สถานการณ์กรณีที่เลวร้ายที่สุดคุณจะต้องประเมินค่าคร่าวๆตามค่าของแรสเตอร์ที่ปรับตำแหน่งใหม่ซึ่งอาจช่วยให้คุณบรรลุเป้าหมาย

มุมล่างซ้ายอยู่ที่พิกัด (0,0) แต่สำหรับฉัน (0,0) คือพิกัดของมุมขวาบนในกรณีของคุณและฉันไม่รู้ระบบพิกัดทางภูมิศาสตร์ใด ๆ ที่มีต้นกำเนิดอยู่ที่ -90

ขณะที่คุณทำงานใน ArcGIS คุณสามารถใช้แถบเครื่องมือการอ้างอิงทางภูมิศาสตร์เพื่อย้ายข้อมูลของคุณในตำแหน่งที่ถูกต้อง ย้ายมุมขวาบนไปที่ (0,0) และซ้ายล่างถึง (-360, -90) และสิ่งนี้ควรให้การแปลงที่แน่นอน แต่คุณยังคงมีปัญหากับ dateline

คุณสามารถลองลิงค์นี้ซึ่งจะช่วยให้คุณสร้างระบบพิกัดทางภูมิศาสตร์ที่กำหนดเอง ในการสรุปในระบบการอ้างอิงของคุณให้สร้าง WGS 84 แบบกำหนดเองด้วยค่าสูงสุดเฉพาะที่ -180 จากนั้นทำการแปลงแบบกำหนดเองโดยใช้การหมุนลองจิจูด

การอ้างอิงทางภูมิศาสตร์ในระบบใหม่นี้จะประกอบด้วยการตั้งค่าจุดเริ่มต้นที่ (-180, -90) จุดควบคุมเดียวก็เพียงพอแล้วเพราะคุณแค่ต้องการแปล เลือกพิกเซลด้านล่างซ้ายอย่างระมัดระวังคลิกขวาแล้วเพิ่มค่า XY ที่แน่นอน) จากนั้นควรมีลักษณะเช่นนี้

ตามคำพูดสุดท้ายไม่มีความแน่นอนว่าระบบพิกัดขึ้นอยู่กับ WGS 84 แต่แน่นอนว่าระบบพิกัดทางภูมิศาสตร์เป็นองศาและความละเอียดของแรสเตอร์ของคุณน่าจะหยาบกว่าข้อผิดพลาดที่อาจเกิดขึ้นเนื่องจากการเลือกข้อมูลผิด