আমার একসাথে প্রায় 550Gb টিআইভি চিত্রাবলী মোজাইক করা দরকার এবং আমি যে সফ্টওয়্যারটি চেষ্টা করেছি তা ব্যর্থ হতে পারে। অঞ্চলটি জোনে বিভক্ত হয়েছে যাতে ক্ষুদ্রতম প্রায় 200 টাইল হয়।
আমি R.৩০ গিগা হার্টজ ইন্টেল জিওন ই 31245, ডিল, 16 জিবি র্যাম, 64-বিট উইন 7 প্রফেশনাল এর ERDAS (ইমেজিন এবং ম্যাপার), আর্কিনএফও এবং গ্লোবাল ম্যাপারের সর্বশেষ সংস্করণগুলি ব্যবহার করেছি। মুল্টি-কোর (মোট 4), হাইপার-থ্রেডেড (8 মোট) মেশিন। আমার সিতে 700 গিগাবাইট ফ্রি এবং ডি রয়েছে 1.5 টিবি।
আমি গ্রাস ব্যবহার করার চেষ্টা করছি (এর আগে কখনও ছিল না) তবে i.image.mosaic কেবল 4 টি ফাইল হ্যান্ডেল করে বলে মনে হচ্ছে ... আমার কারও কারও কাছে 600 টাইল রয়েছে। অন্য কোন বিকল্প বা ওপেনসোর্স সফ্টওয়্যার চেষ্টা করার জন্য?
দুঃখিতটি যুক্ত করা উচিত যে আমরা একটি মোজাইক ডেটাসেট (বা অন্যান্য সফ্টওয়্যারের সমতুল্য) ব্যবহার করতে পারি না কারণ আমাদের নির্ধারিত নো-ডেটা অঞ্চলগুলির সাথে অঞ্চল তৈরি করতে হবে যেমন কোনও জিআইএস সফ্টওয়্যার এ সেগুলি খোলা যায় এবং নিম্ন রেজোলিউশন / পুরানোের সাথে একত্রিত করা যায় ডেটা যখন নতুন ডেটা অবিচ্ছিন্নভাবে উপস্থিত না থাকে data
কিছু মোসাইজ করা ফাইলগুলি বিভিন্ন সফটওয়্যারগুলিতে কীভাবে দেখায় তার একটি উদাহরণ। গ্লোবাল ম্যাপার / ইআরডিএএস ঠিক আছে তবে অর্কগিসে এটি সঠিক নয়।
--- পুরানো তথ্য ---
মোটামুটি আঁকার জন্য দুঃখিত। সুতরাং রঙিন অঞ্চলগুলিকে 5 টি অঞ্চল হিসাবে রাখলে বড় এওআইয়ের কোনও ডেটা অঞ্চল হ্রাস করা যাবে না।
আরকিসে কোডটি নিম্নরূপে রয়েছে (এটি একটি মডেল হিসাবে চালিত হয় এবং পাইথনে নয় কারণ আমি টিফলিস্ট ইনপুট নিতে পারি না)।
arcpy.MosaicToNewRaster_management(tifList+";" +mask,RootOutput,"Tile1.tif","PROJCS['GDA_1994_MGA_Zone_55',GEOGCS['GCS_GDA_1994',DATUM['D_GDA_1994',SPHEROID['GRS_1980',6378137.0,298.257222101]],PRIMEM['Greenwich',0.0],UNIT['Degree',0.0174532925199433]],PROJECTION['Transverse_Mercator'],PARAMETER['False_Easting',500000.0],PARAMETER['False_Northing',10000000.0],PARAMETER['Central_Meridian',147.0],PARAMETER['Scale_Factor',0.9996],PARAMETER['Latitude_Of_Origin',0.0],UNIT['Meter',1.0]]","16_BIT_UNSIGNED","0.5","3","MAXIMUM","#")
# Replace a layer/table view name with a path to a dataset (which can be a layer file) or create the layer/table view within the script
# The following inputs are layers or table views: "test2"
arcpy.CopyRaster_management(OutputFile,RootOutput+"Tile1b.tif","#","256","256","NONE","NONE","16_BIT_UNSIGNED")
যেখানে টিএসএললিস্টটি কোনও সিএসভি ফাইল থেকে পড়তে হবে তবে এটি অজগরটিতে কাজ করে না তাই এর পরিবর্তে আমি উপরের মডেলটিতে চালাচ্ছি ...
আমার ড্রাইভে আমার 1.5TB + ফ্রি স্পেস রয়েছে তবে প্রক্রিয়াটি 9999 ত্রুটির সাথে ক্র্যাশ হয়ে যায়।
এমনকি 100 টাইলস প্রক্রিয়া করবে? -আমাদের আরও জোনে ভাঙ্গার দিকে নজর দেওয়া উচিত?