Esri.com- তে হোস্ট করা একটি ডেটাসেটের জন্য এই তথ্য পৃষ্ঠাটি দেখছেন: http://fema-services2.esri.com/arcgis/rest/services/2012_Sandy/ ইমেজগ্যাট_এনএলটি / ম্যাপসার্ভার / ক্লায়েন্টস
আমি কীভাবে কাঁচা ডেটা (কাঁচা ল্যাট / এলএনজি স্থানাঙ্ক বা এসএইচপি) অ্যাক্সেস পেতে পারি তা জানার চেষ্টা করছি।
দেখে মনে হচ্ছে ডেটা উপলভ্য, এই ইন্টারফেস থেকে কীভাবে এটি টানা যায় আমি ঠিক তা বুঝতে পারি না।
উত্তর:
এই মানচিত্র পরিষেবাটি একটি আর্কজিআইএস সার্ভার ডায়নামিক মানচিত্র পরিষেবা যা সাধারণত কেবলমাত্র ডাব্লুএমএসের মতো চিত্র এবং নির্দিষ্ট ক্যোয়ারী ফলাফল দেয়। কিছু আর্কজিআইএস সার্ভার ইমেজ পরিষেবাদি ডেটা ডাউনলোডের জন্য অনুমতি দেয় তবে এটি এর মধ্যে একটি নয়।
ক্যোয়ারী অপারেশনের মাধ্যমে আপনি যে তথ্যটি সন্ধান করছেন তা পেতে পারেন তবে এই পরিষেবাটিতে সর্বাধিক রেকর্ড রিটার্ন সীমা 1000 হওয়ায় এবং একমাত্র ভি 14 স্তরটিতে 58000 এরও বেশি রেকর্ড রয়েছে বলে একাধিক ক্যোয়ারী লাগবে।
প্রশ্ন জমা দেওয়ার জন্য, স্তর সমাপ্তি বিন্দুতে যান এবং সমর্থিত ক্রিয়াকলাপগুলির জন্য নীচে স্ক্রোল করুন, তারপরে ক্যোয়ারি ক্লিক করুন । মোট রেকর্ডের সংখ্যা খুঁজতে, যেখানে ক্ষেত্রের মধ্যে 1 = 1 লিখুন এবং কেবল রিটার্ন কাউন্টের জন্য সত্য নির্বাচন করুন , 58919 রেকর্ডের মতো দেখাচ্ছে। এটি একবারে কেবল 1000 টি রেকর্ড ফিরিয়ে দেবে, সুতরাং আপনাকে সমস্ত রেকর্ড পেতে 59 টি ক্যোয়ারী জারি করতে হবে। OBJECTID দ্বারা ফিল্টার করা সবচেয়ে সহজ বলে মনে হচ্ছে, সুতরাং প্রথম কোয়েরিটি হ'ল OBJECTID <1000, আউট ফিল্ডস = *, রিটার্ন জ্যামিতি = সত্য, রিটার্ন গণনা = মিথ্যা এবং ফর্ম্যাটটি এইচটিএমএল । পরিষেবাটি সমস্ত বৈশিষ্ট্যগুলি 1- 999 এবং তার বৈশিষ্ট্যগুলি প্রদান করে। আপনি আউটপুট ফর্ম্যাটটিকে কেএমজেডে পরিবর্তন করতে এবং পরে পছন্দসই ফর্ম্যাটে রূপান্তর করতে পারবেন, পুনরাবৃত্তি করুন এবং আপনার সমস্ত ডেটা থাকবে।
আরেকটি পদ্ধতি যা কাজ করবে তা হ'ল প্রোগ্রামের মাধ্যমে বৈশিষ্ট্যটির সমাপ্তি বিন্দুটি ব্যবহার করে পৃথক বৈশিষ্ট্য সম্পর্কিত তথ্য গ্রহন করা। উপরে থেকে এইচটিএমএল ফর্ম্যাট প্রতিক্রিয়ায়, ওবিজেইসিটিড বৈশিষ্ট্যটি কেবলমাত্র সেই বৈশিষ্ট্যের জন্য বৈশিষ্ট্যের একটি হাইপার লিঙ্ক । আরআরএল ইউআরএল ওবিজেসিটিডের সাথে শেষ হয়, যাতে আপনি কেবল এটিকে বাড়িয়ে দিতে পারেন এবং পার্স করা আরও সহজ করার জন্য প্রতিটি প্রতিক্রিয়া জেএসএনে গ্রাহ্য করতে পারেন।
দ্রষ্টব্য - এই বৈশিষ্ট্যগুলি / লিঙ্কগুলির মধ্যে কিছু কেবল আরকজিআইএস সার্ভার 10.1 রিস্টের শেষ প্রান্তগুলিতে প্রযোজ্য।
আমি সর্বদা নিজেকে এমন পরিস্থিতিতে খুঁজে পাই যেখানে আমাকে কোনও মানচিত্র পরিষেবা থেকে একটি চিত্রকলে রূপান্তর করতে হবে। এখানে ব্যবহারের জন্য একটি খুব সহজ ইউটিলিটি যা কোনও পরিষেবা থেকে প্রতিটি বৈশিষ্ট্য রফতানি করবে এবং আপনার প্রয়োজন হলে এটি শেফফাইল এবং জিওজসন হিসাবে সংরক্ষণ করবে। আপনার নোড.জেএস থাকা বা ইনস্টল করতে হবে
https://github.com/tannerjt/AGStoShapefile
একবার আপনি আপনার সিস্টেমে ফাইলগুলি রেখে দিলে, কেবল ফোল্ডারে নেভিগেট করুন, পরিষেবাগুলিতে আপনার মানচিত্র পরিষেবা যুক্ত করুন t txt ফাইল এবং কমান্ড লাইন থেকে এটি চালান:
node AGStoSHP.txt services.txt ./output/
আপনার পরিষেবার জন্য একটি নাম যুক্ত করতে আপনার পরিষেবাদি.টি.এস.টি. ফাইলে একটি পাইপ (|) রাখার বিষয়টি নিশ্চিত করুন।
আশা করি এটি অন্যদের যাদের এই কার্যকারিতাটি প্রয়োজন তাদের সহায়তা করে।
কমান্ড লাইন এবং পাইথন প্যাকেজ পাইজ্রিডম্প ব্যবহার করে এক সাথে একটি আর্কজিআইএসএসইএস ম্যাপ সার্ভারে সঞ্চিত ডেটা ডাউনলোড করুন ।
উদাহরণ কমান্ড:
esri2geojson http://gis.naperville.il.us/arcgis/rest/services/OpenData/OpenDataMapService/MapServer/4 naperville_parking_lots_122417.geojson
আমাকে সম্প্রতি এটি করতে হয়েছিল এবং এটি এখন পর্যন্ত আমার সেরা চেষ্টা ছিল। আমি মূলত "objectid non in {}".format(ids)এমনটি করার চেষ্টা করছিলাম যেখানে আইডিগুলি সংগ্রহ করা আপত্তিযুক্তদের একটি প্রধান অংশ হবে তবে ইউআরএল কোনও ডেটা ফেরত দেবে না, ক্লজ স্ট্রিংটি কত দীর্ঘ হতে পারে তার একটি সীমা থাকতে হবে। এই কোডটির কয়েকটি হার্ড কোডড এবং আইডিগুলি এই স্ক্রিপ্টের চেয়ে অ-অনুক্রমিক হলে সম্ভবত সম্ভবত কাজ করবে না। তবে যাইহোক আমি আশা করি এটি গাইডকে সহায়তা করে
import os, arcpy, json, requests
arcpy.env.workspace=r'C:\path'
arcpy.env.overwriteOutput=True
def non_esri_extract(url,where,idlist):
dic={"where": where,"outFields": "*","returnGeometry": "true","f":"json"}
resp=requests.get(url, params=dic)
data=resp.json()
for i in data['features']:
idlist.append(int(i['attributes']['OBJECTID']))
maximum=max(idlist)
minimum=min(idlist)
return maximum,minimum
def esri_extract(url,e_w):
fields="*"
esri_param="?where={}&outFields={}&returnGeometry=true&f=json".format(e_w, fields)
fsURL=url+esri_param
fs = arcpy.FeatureSet()
fs.load(fsURL)
outname="interm"+str(x)
arcpy.CopyFeatures_management(fs, outname)
x=0
maximum=0
minimum=1
baseURL="http://gismaps.vita.virginia.gov/arcgis/rest/services/VA_Base_layers/VA_Parcels/FeatureServer/0/query"
while maximum!=minimum:
print "number of loops = {}".format(str(x))
if x==0:
ids=[]
maximum,minimum=non_esri_extract(baseURL,"LOCALITY = 'Franklin County'",ids)
esri_where="objectid >={} and objectid <={} and LOCALITY = 'Franklin County'".format(minimum,maximum)
esri_extract(baseURL,esri_where)
x+=1
else:
ids=[]
interm_where="objectid >={} and objectid <={} and LOCALITY = 'Franklin County'".format(maximum,maximum+999)
maximum,minimum=non_esri_extract(baseURL,interm_where,ids)
esri_where="objectid >={} and objectid <={} and LOCALITY = 'Franklin County'".format(minimum,maximum)
esri_extract(baseURL,esri_where)
x+=1
fcs = arcpy.ListFeatureClasses()
arcpy.Merge_management(fcs, "Merged")
return IDs onlyযা সমস্ত 58919 আইডি প্রদান করে।