मैं गॉस्टरवर्प त्रुटि को कैसे हल करूं 'कई बिंदुओं को बदलने में नाकाम रहे', जोस्ट को लैमर्ट कंफर्मल के लिए जियोस्टेशनरी निकालने में विफल रहा?


11

मैं जियोस्टेशनरी से लाम्बर्ट के अनुरूप रीमैप करने की कोशिश कर रहा हूं जो कि गैडलवप का उपयोग कर रहा है। मेरा इनपुट डेटा netcdf में है, और भौगोलिक निर्देशांक (डिग्री) में है और मैं रीमैप किए गए डेटा को netcdf में आउटपुट करना चाहूंगा। मैंने इनपुट netcdf डेटा के लिए संबंधित vrt फ़ाइल बनाई है। Gdalwarp netcdf फ़ाइल का उत्पादन करेगा, लेकिन आउटपुट डेटा सभी शून्य है और मुझे निम्न त्रुटि प्राप्त होती है:

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.nc

क्या भूगर्भिक भौगोलिक निर्देशांक से रीमैप को अनुमानित किया जा सकता है? या क्या मुझे पहले अनुमानित करने के लिए भौगोलिक अनुवाद करने की आवश्यकता है? इसके अलावा, क्या gwarwarp प्रोजेक्शन जानकारी को netcdf से सीधे पढ़ सकता है या क्या आप .vrt पर लिखना चाहते हैं?

यहाँ इनपुट फ़ाइल से gdalinfo आउटपुट क्या है: (यह कक्षा से 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)

और डेटा चर पर अतिरिक्त गदल जानकारी:

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_ts सही पैमाने का अक्षांश है इसलिए lat_1 की तरह एक मानक समानांतर है।
mkennedy

धन्यवाद। मैं मीटर का उपयोग करने की कोशिश करूंगा। मुझे GOES परिभाषा से x_0 और y_0 (स्केल और ऑफ़सेट) मिल रहे हैं, हालांकि ये + proj = geos के लिए आवश्यक इनपुट नहीं हैं, इसलिए मैं उन्हें बाहर निकालने की कोशिश कर सकता हूं। और ई + दीर्घवृत्त के अलावा इंगित करने के लिए धन्यवाद t_srs के लिए अव्यक्त परिभाषाएं लैम्बर्ट (एक विशिष्ट आउटपुट आकार) के AWIPS निश्चित के लिए हैं। मैं जोड़ दूंगा कि क्या गदल जानकारी विशिष्ट डेटा चर के लिए प्रश्न पोस्ट (टिप्पणी के लिए लंबे समय तक) के लिए बाहर निकलती है
केटी जे

मैं जिस AWIPS परिभाषा का उल्लेख कर रहा हूं, वह इस पृष्ठ पर वर्णित है: nws.noaa.gov/noaaport/html/icdtb48_2.html ( लॅमबर्ट मैं जिसे रीमैप करने का प्रयास कर रहा हूं)
केटी '

1
हम्म, इसलिए यह कहते हैं कि लैट / लोन WGS84, लेकिन रिपोर्ट किए गए कोने मुझे चिंता करते हैं क्योंकि वे सिर्फ कच्चे सेल मान हैं। एलसीसी एक स्पर्शरेखा मामला है - एकल मानक समानांतर / अक्षांश का मूल सभी 25N पर हैं। मैंने इस डेटा में से किसी के साथ काम नहीं किया है, इसलिए मैं सिर्फ मेटाडेटा जानकारी द्वारा जा रहा हूं।
mkennedy

छवि को जियोजित नहीं किया जाता है, लेकिन एक स्रोत srs आपूर्ति है। कुछ प्रश्न: * क्या आप CPL_DEBUG = GDAL_netCDF के साथ चल सकते हैं? तो CPL_DEBUG = GDAL_netCDF / usr / bin / gdalwarp ... मुझे संदेह है कि जियोलोकेशन सरणियों के साथ कोई समस्या हो सकती है। * क्या आप अपना डेटा उपलब्ध करा सकते हैं?

जवाबों:


1

यदि आपके डेटा स्रोत के पास अलग-अलग सबडेटासेट के रूप में देशांतर और अक्षांश सेल मान हैं, तो वर्णन के लिए मैन्युअल रूप से एक वीआरटी फ़ाइल बनाने की कोशिश करें, जैसा कि वर्णित है

VSCMO VIIRS डेटा की अस्वीकृति के लिए HDF5 फ़ाइलों और गाल्ड को ताना देने में असमर्थ

पहले चरण में, EPSG: 4326 को सबडेटासेट का उपयोग करते हुए, फिर आप जो भी सीआरएस चाहते हैं, उसका पुनःप्रयोग करें।

हमारी साइट का प्रयोग करके, आप स्वीकार करते हैं कि आपने हमारी Cookie Policy और निजता नीति को पढ़ और समझा लिया है।
Licensed under cc by-sa 3.0 with attribution required.