gdalwarp를 사용하여 Geostationary에서 Lambert Conformal로 다시 매핑하려고합니다. 내 입력 데이터는 netcdf에 있으며 지리적 좌표 (도)이며 다시 매핑 된 데이터를 netcdf에 출력하고 싶습니다. 입력 netcdf 데이터에 해당하는 vrt 파일을 작성했습니다. Gdalwarp는 netcdf 파일을 출력하지만 출력 데이터는 모두 0이며 다음과 같은 오류가 발생합니다.
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.
다음 명령을 시도했습니다.
/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.ncgdalwarp가 지리적 좌표에서 투영으로 다시 매핑 될 수 있습니까? 아니면 지리적 위치를 먼저 투영 된 것으로 변환해야합니까? 또한, gdalwarp는 netcdf에서 직접 투영 정보를 읽을 수 있습니까? 아니면 먼저 .vrt에 기록해야합니까?
입력 파일에서 gdalinfo가 출력하는 내용은 다음과 같습니다 (CLASS의 GOES 13 파일).
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)
데이터 변수에 대한 추가 gdal 정보 :
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
어떤 도움이라도 대단히 감사하겠습니다!
1
지리 투영은 각도를 사용하지 않습니다. 미터를 사용해보십시오. + x_0 / + y_0 값은 어디에서 얻습니까? gdalinfo를 기반으로 입력 래스터가 지리 참조인지 확실하지 않습니다. 대상 srs에는 구형 인 + a = + b가 있지만 + e도 설정되어 있습니다. 그러나 + ellps는 완전히 다른 타원체를위한 것입니다. 다양한 + lat 값도 모두 이상하게 보입니다. lat_ts는 실제 스케일의 위도이므로 lat_1과 마찬가지로 표준 병렬입니다.
—
mkennedy
감사합니다. 나는 미터를 사용해보십시오. GOES 정의에서 x_0 및 y_0 (스케일 및 오프셋)을 얻습니다. 이것은 + proj = geos에 필요한 입력이 아니므로 제거 할 수 있습니다. + e 타원체 추가를 지적 해 주셔서 감사합니다. t_srs에 대한 위도 정의는 lambert (특정 출력 크기)의 AWIPS 정의에 대한 것입니다. 특정 데이터 변수에 대한 gdal 정보가 질문 게시물에 너무 오래 나오는 내용을 추가합니다 (댓글이 너무 깁니다)
—
Katie J
내가 언급 한 AWIPS 정의는이 페이지에 설명되어 있습니다 : nws.noaa.gov/noaaport/html/icdtb48_2.html (첫 번째는 내가 다시 매핑하려는 Lambert입니다)
—
Katie J
흠, 그것은 위도 / 경도 WGS84라고 말하지만,보고 된 모서리 좌표는 원시 셀 값이기 때문에 관심이 있습니다. LCC는 탄젠트 한 경우입니다. 단일 표준 병렬 / 위도는 모두 25N입니다. 나는이 데이터로 작업하지 않았으므로 메타 데이터 정보 만 사용합니다.
—
mkennedy