Làm cách nào để tôi giải quyết lỗi Gdalwarp 'quá nhiều điểm không thể chuyển đổi' để ánh xạ lại Geostationary sang Lambert tuân thủ?

Tôi đang cố gắng ánh xạ lại từ Geostationary sang Lambert tuân thủ bằng gdalwarp. Dữ liệu đầu vào của tôi nằm trong netcdf và ở tọa độ địa lý (độ) và tôi muốn xuất dữ liệu được ánh xạ lại thành netcdf. Tôi đã tạo một tệp vrt tương ứng cho dữ liệu netcdf đầu vào. Gdalwarp sẽ xuất tệp netcdf, nhưng dữ liệu đầu ra là tất cả các số không và tôi nhận được lỗi sau:

Creating output file that is 5120P x 5120L.
Processing input file netcdf.vrt.
ERROR 1: Too many points (441 out of 441) failed to transform,
unable to compute output bounds.
Warning 1: Unable to compute source region for output window 0,0,5120,5120, skipping.
0...10...20...30...40...50...60...70...80...90...100 - done.

Tôi đã thử lệnh sau:

/usr/bin/gdalwarp -s_srs "+proj=geos +h=35785831 +lon_0=-75 +x_0=-0.151844 +y_0=0.151844 +a=6378140 +b=6356754.99999591 +units=degrees +no_defs" -t_srs "+proj=lcc +ellps=clrk66 +a=6378137 +b=6378137 +e=0.0818191910435 +lat_0=24.9999 +lon_0=-95 +lat_1=24.9999 +lat_ts=25.0001 +units=meters +no_defs" -te -1952976.3246 -828316.5944 3248431.6754 4373091.4056 -of netCDF -geoloc -overwrite -r bilinear -ts 5120 5120 netcdf.vrt out.nc

Gdalwarp có thể ánh xạ lại từ tọa độ địa lý để chiếu? Hay tôi cần dịch địa lý để chiếu trước? Ngoài ra, gdalwarp có thể đọc thông tin chiếu trực tiếp từ netcdf hay bạn CẦN viết vào .vrt trước?

Đây là những gì gdalinfo xuất ra từ tệp đầu vào: (đó là tệp GOES 13 từ LỚP)

Size is 512, 512
Coordinate System is `'
Metadata:
  NC_GLOBAL#Conventions=CF-1.4
  NC_GLOBAL#Satellite Sensor=G-13 IMG    
  NC_GLOBAL#Source=McIDAS Area File
Subdatasets:
  SUBDATASET_1_NAME=NETCDF:"goes13.2013.100.165517.BAND_04.nc":auditTrail
  SUBDATASET_1_DESC=[3x80] auditTrail (8-bit character)
  SUBDATASET_2_NAME=NETCDF:"goes13.2013.100.165517.BAND_04.nc":data
  SUBDATASET_2_DESC=[1x665x2036] data (32-bit floating-point)
  SUBDATASET_3_NAME=NETCDF:"goes13.2013.100.165517.BAND_04.nc":lat
  SUBDATASET_3_DESC=[665x2036] lat (32-bit floating-point)
  SUBDATASET_4_NAME=NETCDF:"goes13.2013.100.165517.BAND_04.nc":lon
  SUBDATASET_4_DESC=[665x2036] lon (32-bit floating-point)
Corner Coordinates:
Upper Left  (    0.0,    0.0)
Lower Left  (    0.0,  512.0)
Upper Right (  512.0,    0.0)
Lower Right (  512.0,  512.0)
Center      (  256.0,  256.0)

Và thông tin gdal bổ sung về biến dữ liệu:

Driver: netCDF/Network Common Data Format
Files: goes13.2013.100.174518.BAND_04.nc
Size is 2036, 665
Coordinate System is `'
Metadata:
  data#coordinates=lon lat
  data#long_name=0-255 Brightness Temperature
  data#type=VISR
  NC_GLOBAL#Conventions=CF-1.4
  NC_GLOBAL#Satellite Sensor=G-13 IMG    
  NC_GLOBAL#Source=McIDAS Area File
  NETCDF_DIM_EXTRA={time}
  NETCDF_DIM_time_DEF={1,4}
  NETCDF_DIM_time_VALUES=1365615900
  time#long_name=seconds since 1970-1-1 0:0:0
  time#units=seconds since 1970-1-1 0:0:0
Geolocation:
  LINE_OFFSET=0
  LINE_STEP=1
  PIXEL_OFFSET=0
  PIXEL_STEP=1
  SRS=GEOGCS["WGS 84",DATUM["WGS_1984",SPHEROID["WGS 84",6378137,298.257223563,AUTHORITY["EPSG","7030"]],TOWGS84[0,0,0,0,0,0,0],AUTHORITY["EPSG","6326"]],PRIMEM["Greenwich",0,AUTHORITY["EPSG","8901"]],UNIT["degree",0.0174532925199433,AUTHORITY["EPSG","9108"]],AUTHORITY["EPSG","4326"]]
  X_BAND=1
  X_DATASET=NETCDF:"goes13.2013.100.174518.BAND_04.nc":lon
  Y_BAND=1
  Y_DATASET=NETCDF:"goes13.2013.100.174518.BAND_04.nc":lat
Corner Coordinates:
Upper Left  (    0.0,    0.0)
Lower Left  (    0.0,  665.0)
Upper Right ( 2036.0,    0.0)
Lower Right ( 2036.0,  665.0)
Center      ( 1018.0,  332.5)
Band 1 Block=2036x1 Type=Float32, ColorInterp=Undefined
  NoData Value=9.96920996838686905e+36
  Metadata:
    coordinates=lon lat
    long_name=0-255 Brightness Temperature
    NETCDF_DIM_time=1365615900
    NETCDF_VARNAME=data
    type=VISR

Chúng tôi rất trân trọng bất kỳ sự giúp đỡ nào!

— Katie J
nguồn

Phép chiếu geos sẽ không sử dụng độ; thử mét. Bạn đang nhận các giá trị + x_0 / + y_0 ở đâu? Dựa trên gdalinfo, tôi không chắc chắn raster đầu vào được tham chiếu địa lý chút nào. Trong srs mục tiêu, bạn có + a = + b là một hình cầu, nhưng cũng đặt + e. Tuy nhiên, + ellps dành cho một ellipsoid hoàn toàn khác. Các giá trị + lat khác nhau đều có vẻ kỳ quặc. lat_ts là lat của thang đo thực nên song song tiêu chuẩn giống như lat_1.

— mkennedy

Cảm ơn bạn. Tôi sẽ thử sử dụng mét. Tôi đang nhận được x_0 và y_0 (tỷ lệ và độ lệch) từ định nghĩa GOES mặc dù đây không phải là các đầu vào bắt buộc cho + proj = geos, vì vậy tôi có thể thử loại bỏ chúng. Và cảm ơn vì đã chỉ ra phép cộng + e ellipsoid Các định nghĩa lat cho t_srs là dành cho AWIPS definiton của lambert (một kích thước đầu ra cụ thể). Tôi sẽ thêm những gì thông tin gdal đưa ra cho biến dữ liệu cụ thể vào bài đăng câu hỏi (quá dài để bình luận)

— Katie J

Định nghĩa AWIPS mà tôi đang đề cập đến được mô tả trên trang này: nws.noaa.gov/noaaport/html/icdtb48_2.html (cái đầu tiên là Lambert tôi đang cố gắng ánh xạ lại)

— Katie J

Hmmm, do đó, nó nói lat / lon WGS84, nhưng tọa độ góc được báo cáo liên quan đến tôi vì chúng chỉ là các giá trị ô thô. LCC là trường hợp tiếp tuyến - song song tiêu chuẩn / vĩ độ gốc đều ở mức 25N. Tôi chưa từng làm việc với bất kỳ dữ liệu nào trong số này, vì vậy tôi chỉ xem thông tin siêu dữ liệu.

— mkennedy

Hình ảnh không được tham chiếu địa lý, nhưng một nguồn srs là nguồn cung cấp. Một vài câu hỏi: * Bạn có thể chạy với CPL_DEBUG = GDAL_netCDF không? Vì vậy CPL_DEBUG = GDAL_netCDF / usr / bin / gdalwarp ... Tôi nghi ngờ rằng có thể có vấn đề với các mảng định vị địa lý. * Bạn có thể làm cho dữ liệu của bạn có sẵn?

Nếu nguồn dữ liệu của bạn có các giá trị ô kinh độ và vĩ độ dưới dạng các tệp con riêng biệt, hãy thử xây dựng tệp vrt theo cách thủ công để phân phối lại, như được mô tả trong

Không thể làm cong các tệp HDF5 và gdal để từ chối dữ liệu VSCMO VIIRS

Trong bước đầu tiên, hãy chuyển hướng sang EPSG: 4326 bằng cách sử dụng các tập hợp con, sau đó đến bất kỳ CRS nào bạn muốn.

— AndreJ
nguồn