আমার কাছে পাইথন কোড রয়েছে যা নিম্নলিখিত কর্মপ্রবাহের মাধ্যমে পয়েন্ট শেফফাইলগুলি ডিজাইন করার জন্য ডিজাইন করা হয়েছে:
- পয়েন্টগুলি মার্জ করুন
- পয়েন্টগুলি একীভূত করুন, যেমন একে অপরের 1 মিটারের মধ্যে যে কোনও পয়েন্ট এক বিন্দুতে পরিণত হয়
- বৈশিষ্ট্য স্তর তৈরি করুন, যেখানে z <10 সহ পয়েন্টগুলি নির্বাচন করা হয়েছে
- বাফার পয়েন্টস
- বহুভুজ থেকে রাস্টার 1 মি রেজোলিউশন
- পুনর্গঠিত করুন, যেখানে 1 - 9 = 1; নোডাটা = 0
প্রতিটি শেফফাইলে প্রায় 250,000 থেকে 350,000 পয়েন্ট থাকে ~ 5x7 কিমি covering ইনপুট হিসাবে ব্যবহৃত পয়েন্ট ডেটা গাছের অবস্থানের প্রতিনিধিত্ব করে। প্রতিটি পয়েন্ট (অর্থাত্ গাছ) এর সাথে একটি "জেড" মান থাকে যা মুকুট ব্যাসার্ধকে উপস্থাপন করে এবং বাফার প্রক্রিয়াতে ব্যবহৃত হয়। আমার উদ্দেশ্য হ'ল চাঁদীর কভারটি বর্ণনা করে একটি রাস্টার উত্পাদন করতে একটি পৃথক প্রক্রিয়াতে চূড়ান্ত বাইনারি আউটপুট ব্যবহার করা।
আমি চারটি শেফফাইল দিয়ে একটি পরীক্ষা চালিয়েছি এবং এটি একটি 700 এমবি রাস্টার তৈরি করেছে এবং 35 মিনিট (আই 5 প্রসেসর এবং 8 জিবি র্যাম) নিয়েছে। আমার এই প্রক্রিয়াটি 3500 শেফফাইলে চালানোর দরকার হবে তা দেখে, প্রক্রিয়াটি সহজ করার জন্য আমি কোনও পরামর্শের (সংযুক্ত কোড দেখুন) প্রশংসা করব। সাধারণভাবে বলতে গেলে, জিওপ্রসেসিং বিগ ডেটা মোকাবেলার সেরা উপায় কী? আরও সুনির্দিষ্টভাবে বলা যায় যে কোড বা ওয়ার্কফ্লোতে কোনও টুইট রয়েছে যা দক্ষতা বাড়াতে সহায়তা করতে পারে?
সম্পাদনা করুন :
জিওপ্রসেসিং কাজের জন্য সময় (মোটের%):
- মার্জ = 7.6%
- একীভূত = 7.1%
- লির থেকে বৈশিষ্ট্য = 0
- বাফার = 8.8%
- পলি থেকে রাস্টার = 74.8%
- পুনর্গঠিত = 1.6%
# Import arcpy module
import arcpy
# Check out any necessary licenses
arcpy.CheckOutExtension("spatial")
# Script arguments
temp4 = arcpy.GetParameterAsText(0)
if temp4 == '#' or not temp4:
temp4 = "C:\\gdrive\\temp\\temp4" # provide a default value if unspecified
Reclassification = arcpy.GetParameterAsText(1)
if Reclassification == '#' or not Reclassification:
Reclassification = "1 9 1;NODATA 0" # provide a default value if unspecified
Multiple_Value = arcpy.GetParameterAsText(2)
if Multiple_Value == '#' or not Multiple_Value:
Multiple_Value = "C:\\t1.shp;C:\\t2.shp;C:\\t3.shp;C:\\t4.shp" # provide a default value if unspecified
# Local variables:
temp_shp = Multiple_Value
Output_Features = temp_shp
temp2_Layer = Output_Features
temp_Buffer = temp2_Layer
temp3 = temp_Buffer
# Process: Merge
arcpy.Merge_management(Multiple_Value, temp_shp, "x \"x\" true true false 19 Double 0 0 ,First,#,C:\\#########omitted to save space
# Process: Integrate
arcpy.Integrate_management("C:\\gdrive\\temp\\temp.shp #", "1 Meters")
# Process: Make Feature Layer
arcpy.MakeFeatureLayer_management(temp_shp, temp2_Layer, "z <10", "", "x x VISIBLE NONE;y y VISIBLE NONE;z z VISIBLE NONE;Buffer Buffer VISIBLE NONE")
# Process: Buffer
arcpy.Buffer_analysis(temp2_Layer, temp_Buffer, "z", "FULL", "ROUND", "NONE", "")
# Process: Polygon to Raster
arcpy.PolygonToRaster_conversion(temp_Buffer, "BUFF_DIST", temp3, "CELL_CENTER", "NONE", "1")
# Process: Reclassify
arcpy.gp.Reclassify_sa(temp3, "Value", Reclassification, temp4, "DATA")