আমি জিওডেশনারি থেকে ল্যামবার্ট কনফর্মাল জিডিওয়াল্প ব্যবহার করে পুনরায় তৈরি করার চেষ্টা করছি। আমার ইনপুট ডেটা নেটসিডিএফ-এ রয়েছে, এবং ভৌগলিক স্থানাঙ্কে (ডিগ্রি) রয়েছে এবং আমি বাকী ডেটা নেটসিডিএফ-এ আউটপুট দিতে চাই। আমি ইনপুট নেটসিডিএফ ডেটার জন্য একটি সম্পর্কিত ভিআরটি ফাইল তৈরি করেছি। জিডালওয়ার্প নেটসিডিএফ ফাইল আউটপুট দেবে, তবে আউটপুট ডেটা সমস্ত শূন্য এবং আমি নিম্নলিখিত ত্রুটিটি পেয়েছি:
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জিডালওয়ার্প কি ভৌগলিক স্থানাঙ্ক থেকে পুনর্নির্মাণের প্রস্তাব করা যেতে পারে? নাকি আমার আগে ভৌগলিকটিকে প্রজেক্টে অনুবাদ করতে হবে? এছাড়াও, gdalwarp সরাসরি নেটসিডিএফ থেকে প্রজেকশন তথ্য পড়তে পারেন বা আপনি .vrt লিখতে প্রথমে প্রয়োজন?
ইনপুট ফাইল থেকে gdalinfo আউটপুট কি তা এখানে: (এটি CLASS থেকে একটি 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 রয়েছে যা একটি গোলক, তবে সেটও + ই। তবে, + এলপস সম্পূর্ণ ভিন্ন উপবৃত্তাকার জন্য। বিভিন্ন + ল্যাট মানগুলি খুব অদ্ভুত বলে মনে হয়। ল্যাট_টিস সত্য স্কেলের ল্যাট তাই ল্যাট_1 এর মতো একটি আদর্শ সমান্তরাল।
—
মাকনেডি
ধন্যবাদ. আমি মিটার ব্যবহার করার চেষ্টা করব। আমি GOES সংজ্ঞা থেকে x_0 এবং y_0 (স্কেল এবং অফসেট) পাচ্ছি যদিও এগুলি + প্রোজ = জিওসের জন্য প্রয়োজনীয় ইনপুট নয়, তাই আমি সেগুলি চেষ্টা করে দেখতে পারি। এবং + ই এলিপসয়েড সংযোজনটি নির্দেশ করার জন্য ধন্যবাদ t_srs এর ল্যাট সংজ্ঞাগুলি ল্যাম্বার্টের একটি অ্যাডব্লিউপিএস নির্দিষ্টকরণের জন্য (একটি নির্দিষ্ট আউটপুট আকার)। প্রশ্ন পোস্টে নির্দিষ্ট ডেটা ভেরিয়েবলের জন্য জিডিএল তথ্যটি কী ছড়িয়ে দেয় আমি যুক্ত করব (মন্তব্যের জন্য খুব দীর্ঘ)
—
কেটি জে
আমি যে AWIP সংজ্ঞাটি উল্লেখ করছি তার এই পৃষ্ঠায় বর্ণিত হয়েছে: nws.noaa.gov/noaaport/html/icdtb48_2.html (প্রথমটি ল্যামবার্ট আমি পুনর্নির্মাণের চেষ্টা করছি)
—
কেটি জে
হুম, সুতরাং এটি ল্যাজ / লম্বা ডাব্লু জিএস ৮৪ বলেছে, তবে রিপোর্ট করা কোণটি আমাকে উদ্বেগ জানায় কারণ এগুলি কেবল কাঁচা ঘর মান। এলসিসি একটি স্পর্শকাতর ক্ষেত্রে - একক মানের সমান্তরাল / অক্ষাংশের উত্স সমস্ত 25N এ থাকে। আমি এই ডেটার কোনওটির সাথেই কাজ করি নি, তাই আমি কেবল মেটাডেটা তথ্য দিয়ে যাচ্ছি।
—
মেকনেডি
চিত্রটি জিওরফারেন্সযুক্ত নয়, তবে একটি উত্স srs সরবরাহযোগ্য। কয়েকটি প্রশ্ন: * আপনি কি CPL_DEBUG = GDAL_netCDF দিয়ে চালাতে পারবেন? সুতরাং সিপিএল_ডিইউইউজি = জিডিএল_নেটসিডিএফ / ইউএসআর / বিন / জিডালওয়ার্প ... আমার সন্দেহ হয় যে ভূ-অবস্থানের অ্যারেতে কোনও সমস্যা হতে পারে। * আপনি কি আপনার ডেটা উপলব্ধ করতে পারবেন?