Gdal-dev মেলিং তালিকা থেকে সরানো:
সোমবার, ২ সেপ্টেম্বর, ২০১৩ সন্ধ্যা :0:০৯ এ ডেভিড শান লিখেছেন:
হাই তালিকা, আমি বিতরণের জন্য একক নেটসিডিএফ ফাইল হিসাবে অভিন্ন প্রজেকশন / ব্যাপ্তি / রেজোলিউশন সহ জিটিফ রেস্টারদের টাইমসারিগুলি প্যাকেজ করার চেষ্টা করছি। আমি গত ঘন্টা অনলাইন ডকটির সাথে পরামর্শ করে এবং কোনও সাফল্য ছাড়াই gdal_translate, gdalbuildvrt এবং gdalwarp নিয়ে খেলছি।
বিদ্যমান gdal কমান্ড লাইন ইউটিলিটিগুলি ব্যবহার করে কি এটি করার সহজ উপায় আছে? আমি অনুভব করেছি যে নেটসিডিএফএফ পাইথন এপিআই ব্যবহার করে একটি কাস্টম সমাধান নিয়ে যাওয়ার আগে আমি জিজ্ঞাসা করব।
ধন্যবাদ। -David
মঙ্গলবার, 3 সেপ্টেম্বর, 2013 সকাল 10: 15 টা এতিয়েন ট্যুরিগনি লিখেছেন:
আপনি যা চান তা সম্ভবত gdal এর আওতার বাইরে। এটির জন্য কিছু চালাক মেটাডেটা পরিচালনার প্রয়োজন হবে যাতে gdal_translate এগুলিকে একক ফাইলে রাখে ...
আমি আপনাকে পরামর্শ দিচ্ছি যে আপনি তাদের সকলকে gdal_translate ব্যবহার করে নেটcdf এ রূপান্তর করুন এবং তারপরে অস্থায়ী মাত্রায় স্ট্যাক করার জন্য পাইথন-নেটসিডিএফ 4 (নিম্পি / স্কিপি থেকে নয়) ব্যবহার করুন।
মঙ্গলবার, 3 সেপ্টেম্বর, 2013, সকাল 7:55 এ, "সিগনেল, রিচার্ড" লিখেছেন:
ডেভিড, আপনি যদি জিআইএস স্ট্যাকেক্সচেঞ্জ গ্রুপ /gis// এ আপনার প্রশ্ন পোস্ট করেন তবে আমি একটি উদাহরণ কোড সরবরাহ করব যা সহায়ক হওয়া উচিত।
-Rich
====================
আপডেট 9/3/13 17:04 পিডিটি
আমার এক ইনপুট ডেটাসেটের জন্য এখানে gdalinfo আউটপুট দেওয়া হয়েছে:
gdalinfo 20120901T2024_align_x+22.19_y+3.68_z+14.97_warp.tif
Driver: GTiff/GeoTIFF
Files: 20120901T2024_align_x+22.19_y+3.68_z+14.97_warp.tif
Size is 10666, 13387
Coordinate System is:
PROJCS["unnamed",
GEOGCS["WGS 84",
DATUM["WGS_1984",
SPHEROID["WGS 84",6378137,298.257223563,
AUTHORITY["EPSG","7030"]],
AUTHORITY["EPSG","6326"]],
PRIMEM["Greenwich",0],
UNIT["degree",0.0174532925199433],
AUTHORITY["EPSG","4326"]],
PROJECTION["Polar_Stereographic"],
PARAMETER["latitude_of_origin",70],
PARAMETER["central_meridian",-45],
PARAMETER["scale_factor",1],
PARAMETER["false_easting",0],
PARAMETER["false_northing",0],
UNIT["metre",1,
AUTHORITY["EPSG","9001"]]]
Origin = (-211346.063781524338992,-2245136.291794800199568)
Pixel Size = (5.000000000000000,-5.000000000000000)
Metadata:
AREA_OR_POINT=Area
Image Structure Metadata:
COMPRESSION=LZW
INTERLEAVE=BAND
Corner Coordinates:
Upper Left ( -211346.064,-2245136.292) ( 50d22'39.70"W, 69d23'55.59"N)
Lower Left ( -211346.064,-2312071.292) ( 50d13'22.38"W, 68d48'10.75"N)
Upper Right ( -158016.064,-2245136.292) ( 49d 1'33.33"W, 69d26'16.42"N)
Lower Right ( -158016.064,-2312071.292) ( 48d54'35.06"W, 68d50'27.28"N)
Center ( -184681.064,-2278603.792) ( 49d38' 1.32"W, 69d 7'17.04"N)
Band 1 Block=256x256 Type=Float32, ColorInterp=Gray
NoData Value=-32767লুক এর প্রস্তাবিত পদ্ধতির অনুসরণ করা।
ভিআরটি প্রজন্ম ভাল কাজ করে:
gdalbuildvrt -separate newtest.vrt *warp.tif
<VRTDataset rasterXSize="10666" rasterYSize="13387">
<SRS>PROJCS["unnamed",GEOGCS["WGS 84",DATUM["WGS_1984",SPHEROID["WGS 84",6378137,298.257223563,AUTHORITY["EPSG","7030"]],AUTHORITY["EPSG","6326"]],PRIMEM["Greenwich",0],UNIT["degree",0.0174532925199433],AUTHORITY["EPSG","4326"]],PROJECTION["Polar_Stereographic"],PARAMETER["latitude_of_origin",70],PARAMETER["central_meridian",-45],PARAMETER["scale_factor",1],PARAMETER["false_easting",0],PARAMETER["false_northing",0],UNIT["metre",1,AUTHORITY["EPSG","9001"]]]</SRS>
<GeoTransform> -2.1134606378152434e+05, 5.0000000000000000e+00, 0.0000000000000000e+00, -2.2451362917948002e+06, 0.0000000000000000e+00, -5.0000000000000000e+00</GeoTransform>
<VRTRasterBand dataType="Float32" band="1">
<NoDataValue>-3.27670000000000E+04</NoDataValue>
<ComplexSource>
<SourceFilename relativeToVRT="1">20110619T2024_align_x+15.51_y+1.15_z+12.10_warp.tif</SourceFilename>
<SourceBand>1</SourceBand>
<SourceProperties RasterXSize="10666" RasterYSize="13387" DataType="Float32" BlockXSize="256" BlockYSize="256" />
<SrcRect xOff="0" yOff="0" xSize="10666" ySize="13387" />
<DstRect xOff="0" yOff="0" xSize="10666" ySize="13387" />
<NODATA>-32767</NODATA>
</ComplexSource>
</VRTRasterBand>
<VRTRasterBand dataType="Float32" band="2">
<NoDataValue>-3.27670000000000E+04</NoDataValue>
<ComplexSource>
<SourceFilename relativeToVRT="1">20110802T2024_align_x+16.33_y+2.14_z+12.02_warp.tif</SourceFilename>
<SourceBand>1</SourceBand>
<SourceProperties RasterXSize="10666" RasterYSize="13387" DataType="Float32" BlockXSize="256" BlockYSize="256" />
<SrcRect xOff="0" yOff="0" xSize="10666" ySize="13387" />
<DstRect xOff="0" yOff="0" xSize="10666" ySize="13387" />
<NODATA>-32767</NODATA>
</ComplexSource>
</VRTRasterBand>
...কিন্তু যখন আমি এনসি তে অনুবাদ করার চেষ্টা করি তখন আমি নিম্নলিখিত ত্রুটিটি পাই:
gdal_translate -of netcdf newtest.vrt newtest.nc
Input file size is 10666, 13387
Warning 1: Variable has 0 dimension(s) - not supported.
0...10...20...30...40...50ERROR 1: netcdf error #-62 : NetCDF: One or more variable sizes violate format constraints .
at (netcdfdataset.cpp,SetDefineMode,1574)
ERROR 1: netcdf error #-39 : NetCDF: Operation not allowed in define mode .
at (netcdfdataset.cpp,IWriteBlock,1435)
ERROR 1: netCDF scanline write failed: NetCDF: Operation not allowed in define mode
ERROR 1: An error occured while writing a dirty block
...ERROR 1: netcdf error #-39 : NetCDF: Operation not allowed in define mode .
at (netcdfdataset.cpp,IWriteBlock,1435)
ERROR 1: netCDF scanline write failed: NetCDF: Operation not allowed in define mode
ERROR 1: netcdf error #-62 : NetCDF: One or more variable sizes violate format constraints .
at (netcdfdataset.cpp,~netCDFDataset,1548)সুতরাং কাছাকাছি পরিদর্শন করার পরে, দেখা যাচ্ছে যে আমি যে পোলার স্টেরিওগ্রাফিক প্রজেকশনটি ব্যবহার করছি (ইপিএসজি: 3413) তাতে জিডাল অসন্তুষ্ট। নেটকডিএফডিটাসেট.সি.পি. এর 1570-1582 লাইনগুলি দেখুন:
আমার প্রজেকশনটিতে একটি অক্ষাংশ_অফ_রিগিন নির্দিষ্ট করা আছে তবে নেটসিডিএফ ড্রাইভারের দ্বারা প্রত্যাশার মতো কোনও মান সমান্তরাল নেই।