গ্রাসে প্যান-তীক্ষ্ণ ল্যান্ডস্যাট চিত্রগুলির মিথ্যা-সংমিশ্রণের রঙ-ব্যালেন্স কীভাবে করবেন?

Question

যখনই আমি ব্যবহার ঘাস মধ্যে কিছু ল্যান্ডস্যাট ইমেজ থেকে প্যান-তীক্ষ্ণ করুন composites চেষ্টা i.pansharpen, i.fusion.broveyবা আইএইচএস ধার পদ্ধতি, আউটপুট নিম্নোক্ত কোন একটি বা বৈশিষ্ট্য সব হবে:

• সংশ্লেষের রঙটি আন-শার্পেনযুক্ত সংমিশ্রণের তুলনায় আলাদা রঙে রয়েছে
• উজ্জ্বলতা স্তর গণ্ডগোল হয়
• পুরো সংমিশ্রণটি সমস্ত সাদা / সবুজ কালো হয়ে গেছে (প্রাক-প্রক্রিয়াযুক্ত চিত্রগুলি বায়ুমণ্ডলের উপরের প্রতিচ্ছবি বা পৃষ্ঠের প্রতিবিম্ব সংশোধনগুলিতে ব্যবহার করার সময় i.landsat.toar)

আমি নিম্নলিখিত সমস্ত চেষ্টা করেছি ; তবে রঙ / উজ্জ্বলতা একইরকম থেকে গেছে বা আরও খারাপ হয়েছে:

• i.landsat.rgbপ্যান-শার্পিং প্রক্রিয়া করার আগে এবং পরে প্রয়োগ করা হয়েছে
• ভিতরে -fবা -pপতাকা নিয়ে খেলেছেi.landsat.rgb
• r.colorsবর্ণের টেবিলটি ধূসর / গ্রে -255 / গ্রেই.ইকিতে সম্পাদনা করার চেষ্টা করা হয়েছে
• i.pansharpenসমস্ত ব্রোভে / আইএইচএস / পিসিএ পদ্ধতি ব্যবহার করার চেষ্টা করা হয়েছে
• নীল চ্যানেলটির ভারসাম্য বজায় রাখতে -lপতাকাটি খেলেনi.pansharpen

ঘাস জিআইএস ম্যানুয়াল কিভাবে প্যান স্পষ্টতার এবং রঙ-ব্যালেন্সিং সম্পাদন করতে উপর ব্যাখ্যা হয়নি, কিন্তু আমি চিন্তা করতে পারে না কিভাবে সমবর্তী কর্মপ্রবাহ উভয় প্রসেস একত্রিত করতে হবে। আমি সন্দেহ করেছিলাম যে এটি গ্রাসে রঙিন টেবিল, রঙ-হিস্টোগ্রাম ইত্যাদি সম্পর্কে আমার কম বোঝার কারণে ..

সুতরাং, কেউ আমাকে ব্যাখ্যা করতে পারেন - গ্রাসে চিত্র-প্রসেসিংয়ের পরে ল্যান্ডস্যাট চিত্রগুলির সাথে ডিল করার সময় আপনি কীভাবে রঙ-ব্যালেন্সিং সমস্যাগুলি মোকাবেলা করবেন? আপনি কি আমার সাথে আপনার প্রিয় কর্মপ্রবাহ / পদ্ধতিগুলি ভাগ করতে পারেন?

যে কোনও প্রতিক্রিয়ার জন্য অনেক ধন্যবাদ!

Answer 1

সংক্ষিপ্ত বিবরণ

প্যান-শার্পিংয়ের পরে গ্রহণযোগ্য রঙ-ভারসাম্যযুক্ত যৌগিক চিত্র পাওয়ার জন্য গ্রাস-জিআইএস সংস্করণ 7-এর মধ্যে একটি কার্যকারী পদ্ধতি

1. ইনপুট ডেটা 8-বিটের ভিতরে রয়েছে কিনা তা পরীক্ষা করুন [0, 255]
2. যদি ডেটা ভিতরে থাকে [0, 255] তারপর প্যান-শার্পিংয়ে এগিয়ে যান ( i.pansharpen)
3. যদি তথ্য [0, 255] এর ভিতরে না থাকে তবে তাদের এই ব্যাপ্তিতে পুনরুদ্ধার করুন ( r.rescale)
4. বৈশিষ্ট্যযুক্ত যেকোন পদ্ধতির সাথে প্যান-শার্প করুন (ব্রোভে, আইএইচএস, পিসিএ)
5. রঙ-ভারসাম্য স্বয়ংক্রিয়ভাবে i.landsat.rgbমডিউলটি ব্যবহার করে বা আগ্রহের ব্যান্ডগুলির রঙের টেবিলগুলিতে ম্যানুয়ালি সমন্বয় করে

বিশদ এবং উদাহরণ নির্দেশাবলী

প্যান-শার্পিং / ফিউশন

গ্রাস তে একটি উত্সর্গীকৃত প্যান-শার্পিং মডিউল ধারণ করে, i.pansharpenযা ব্রোভে রূপান্তরকরণ , শাস্ত্রীয় আইএইচএস পদ্ধতি এবং একটি যা পিসিএ ভিত্তিক একটি তীক্ষ্ণ করার জন্য তিনটি কৌশল যুক্ত করে ।

i.pansharpenএকটি ইনপুট হিসাবে 8-বিট রাস্টার মানচিত্রের সাথে সূক্ষ্ম কাজ করে। যদি প্রক্রিয়া করার ডেটা এই ব্যাপ্তির বাইরে থাকে তবে এর বাইরে [0, 255], গ্রাসের r.rescaleমডিউল ব্যবহার করে এগুলি এই পরিসরে ফিট করার জন্য তাদের পুনরুদ্ধার করা যেতে পারে ।

11-বিট বর্ণালী ব্যান্ডগুলির একটি সেট দেওয়া (উদাহরণস্বরূপ নীল, সবুজ, লাল, এনআইআর এবং প্যান) এর মধ্যে রয়েছে [0, 2047], উদাহরণস্বরূপ নীল ব্যান্ডটি জিজ্ঞাসা করলে ফিরে আসবে

r.info Blue_DNs -r
min=0
max=2047

এর মধ্যে বিস্তৃত নীল ব্যান্ডটিকে পুনরায় উদ্ধার করা [0, 255]

r.rescale in=Blue_DNs out=Blue_DNs_255 from=0,2047 to=0,255

একই পদক্ষেপটি মাল্টি-বর্ণালী ব্যান্ডের বাকি অংশ এবং আগ্রহের প্যানক্রোমেটিক ব্যান্ড উভয়ের ক্ষেত্রেই প্রযোজ্য।

গ্রাআরএসএসের সাথে কাজ করার সময় যথারীতি আগ্রহের অঞ্চলটি নির্ধারণ করা প্রয়োজন, অর্থাত ব্যান্ড (গুলি) এর সীমাটি বা অন্য কোনওটির সাথে মেলে। রেজোলিউশনটি মডিউল দ্বারা এই বিশেষ ক্ষেত্রে যত্ন নেওয়া হয় এবং ফলস্বরূপ প্যান-তীক্ষ্ণ রাস্টার মানচিত্রগুলি প্যানক্রোমেটিক ব্যান্ডের মতো একই উচ্চ (এআর) রেজোলিউশনের হবে।g.regionrast=Blue_DNs_255

আইএইচএস-ভিত্তিক প্যান-শার্পানিং ক্রিয়াকলাপের উদাহরণ কমান্ডটি দেখতে পেল

i.pansharpen pan=Pan_DNs_255 ms1=Blue_DNs_255 ms2=Green_DNs_255 ms3=Red_DNs_255 output=sharptest255 sharpen=ihs

রঙের ব্যালেন্সিং

প্রক্রিয়া সমাপ্তির পরে, মডিউল আউটপুট

...
The following pan-sharpened output maps have been generated:
sharptest255_red
sharptest255_green
sharptest255_blue

To visualize output, run: g.region -p rast=sharptest255.red
d.rgb r=sharptest255_red g=sharptest255_green b=sharptest255_blue

সাধারণত প্যান- শার্পিংয়ের পরে রঙগুলি পুনরায় ভারসাম্য করার জন্য পর্যাপ্ত হওয়া উচিত উদাহরণস্বরূপ i.landsat.rgb মডিউল বা একটি আরজিবি চিত্র রচনা করবে এমন তিনটি ব্যান্ডের প্রত্যেকটির ম্যানুয়াল সামঞ্জস্য।

স্ক্রীনশট

... যোগ করা হবে

Comments

আমি জানতাম আরও ভাল উপায় থাকতে হবে! এখন আমি নিখরচায় i.sharpen মডিউলটি ব্যবহার করতে পারি। R.rescale মডিউলটি নির্দেশ করার জন্য ধন্যবাদ। অসাধারণ কাজ নিকোস!

হাজিক, আমি নিশ্চিত নই যে 11-বিট ডেটা সেটটি 8-বিটে রূপান্তর করার সময় এবং "কতটা" সূক্ষ্ম বিবরণ হারিয়ে গেছে। উদাহরণস্বরূপ, কুইকবার্ড চিত্রাবলী 11-বিট সেন্সর। এগুলি 8-বিট এবং 16-বিট উভয় ফর্ম্যাটে উপলব্ধ। কী করা উচিত তা সিদ্ধান্ত নেওয়ার উপর নির্ভর করে ব্যবহারকারী। এটি অবশ্যই i.pansharpenসমস্ত ধরণের ফর্ম্যাটগুলি পরিচালনা করার জন্য সুন্দর হবে। দয়া করে কোনও সম্পর্কিত "টিকিট" দেখুন: টিকিট # 2048: i.pansharpen 8-বিটের চিত্রের মধ্যে সীমাবদ্ধ । অন্যদিকে, আমি কেবল স্টাফগুলি বুঝতে পারি না এবং এইভাবে, আইপ্যানশার্পেন সঠিকভাবে ব্যবহার করতে সক্ষম না হই ...?

— নিকস আলেকজান্দ্রিস

Answer 2

আমি উচ্চ-নীচে অনুসন্ধান করেছি এবং আমি মনে করি যে আমি আমার সমস্যার মূল খুঁজে পেয়েছি। আমি বিশ্বাস করি এখনই তাদের জন্য সমাধান পেয়েছি - তবে এটি কিছুটা অগোছালো। আমি নিশ্চিত যে এগুলি সমাধানের আরও ভাল উপায় আছে। আপনি যদি সহজ উপায় জানেন তবে ভাগ করুন!

সমস্যার মূল:

1. এর আউটপুট i.landsat.toar ভাসমান বিন্দুতে । আমি বুঝতে পেরেছি যে আমি যখন কোনও প্যান-শার্পিং পদ্ধতিতে ফ্লোটিং পয়েন্ট রাস্টারগুলি ব্যবহার করি তখন রঙগুলি গণ্ডগোল হয়। এই অ্যালগরিদম কোনওভাবে মূল পূর্ণসংখ্যার আকারে রেস্টারদের পছন্দ করে।
2. প্যান-শার্পিং মডিউলগুলি যেমন i.pansharpen এবং i.fusion.brovey মডিউলগুলি রঙগুলিকে মেস আপ করবে । তারা সেই মডিউলগুলিতে যে অ্যালগোরিদমগুলি ব্যবহার করেছিল তা আমি যথেষ্টভাবে উপলব্ধি করতে পারি নি - তবে কোনওভাবে রঙ-টেবিলগুলি প্রভাবিত হবে এবং ফলস্বরূপ প্যান-তীক্ষ্ণ চিত্রগুলি নষ্ট করে দেবে।

সমাধান:

1. আউটপুটটি i.landsat.toar ফ্ল্যাট থেকে ইনট এ রূপান্তর করে ব্যবহার করেr.recode ।
2. আইএইচএস প্যান-শার্পিং পদ্ধতিতে ( i.rgb.hisএবং i.his.rgb) ইনপুট হিসাবে রেস্টারদের ব্যবহার করুন । আমি ব্যবহার থেকে দূরে রাখতে i.pansharpenএবং i.fusion.brovey।

ওয়ার্কফ্লো:

1. প্রতিটি রাস্টার ব্যান্ডের সাথে ডিএন ন্যূনতম এবং সর্বাধিক মানগুলি পেতে প্রক্রিয়া করা হয়েছে যা পতাকাটির r.info সাথে ব্যবহার করুন । উদাহরণ স্বরূপ:-r i.landsat.toar

   > r.info -r BAND1
   min=0.01
   max=0.370064120902708
   

   যেমনটি আমরা দেখতে পাচ্ছি, মানগুলি 0-1 এর মধ্যে হয় যা মূলগুলির চেয়ে বেশ আলাদা (যা 0-255 এর মধ্যে থাকে)। এটি ব্যাখ্যা করে যে কেন প্যান-শার্পিং থেকে আউটপুট ফাঁকা পরিণত হয়েছিল, কারণ ব্যবহৃত পরিসরের মানটি খুব কম (1 এর নীচে)।

2. ব্যবহার করে সেই রাস্টার ব্যান্ডটি রূপান্তর করুন r.recode । পদক্ষেপ 1 থেকে প্রাপ্ত সর্বনিম্ন এবং সর্বাধিক মানগুলি 0-255 এর নতুন পরিসরে রূপান্তর করতে ব্যবহার করুন। কোড স্নিপেট একটি উদাহরণ:

   r.recode input=BAND1 output=NEWBAND1 rules=- << EOF
   0.01:0.370064120902708:0:255
   EOF
   

   আমরা এর সাথে নতুন রূপান্তরিত মানগুলি পরীক্ষা করতে পারি r.info:

   > r.info -r NEWBAND1
   min=1
   max=254
   

   মানগুলি 0-255- এ রয়েছে: এখন এটি প্যান-শার্পিং প্রক্রিয়াটির জন্য ব্যবহারযোগ্য।

3. রূপান্তরিত ব্যান্ডের সাথে ধূসর-স্কেল রঙের টেবিলটি প্রয়োগ করুন r.colors ।

   r.colors NEWBAND1 color=grey
   

   এখনও অবধি, ধূসর রঙের টেবিলটি ব্যবহার করে আমি সেরা ফলাফল পেয়েছি - প্যান-তীক্ষ্ণ মিশ্রণগুলি মূল সংমিশ্রণের সাথে ঘনিষ্ঠভাবে মিলেছে। অন্যান্য বিকল্পগুলি হল এর সাথে পতাকা color=grey.eqব্যবহার করে বা ধূসর রঙের টেবিলকে সমান -eকরা color=grey। অথবা আমরা এর i.landsat.rgbপরিবর্তে মডিউলটি ব্যবহার করতে পারি r.colors..

4. প্যান রাস্টার (ব্যান্ড 8) সহ কমপোজিট হিসাবে ব্যবহারের আমাদের মনস্থ হওয়া অন্যান্য রাস্টার ব্যান্ডগুলির সাথে 1-3 বারের পুনরাবৃত্তি করুন । স্ক্রিপ্ট ব্যবহার এখানে প্রশংসিত হবে।

5. তারপরে আইএইচএস প্যান-শার্পিং পদ্ধতিতে ইনপুট হিসাবে প্রক্রিয়াজাত রাস্টারগুলি ব্যবহার করুন । উদাহরণস্বরূপ, ব্যান্ড 7,4,2 এর সংমিশ্রণটি তৈরি করার সময়:

   i.rgb.his r=NEWBAND7 g=NEWBAND4 b=NEWBAND2 hue=HUE int=INT sat=SAT
   

   এটি 3 স্তর আউটপুট দেবে: হিউ স্তর HUE, একটি তীব্রতা স্তর INTএবং একটি স্যাচুরেশন স্তর SAT। আমরা তখন তীব্রতা স্তর প্রতিস্থাপন করবে INTপ্যান রাস্টার ব্যান্ডের সাথে NEWBAND8এ i.his.rgb:

   i.his.rgb hue=HUE sat=SAT int=NEWBAND8 r=COMP742_red g=COMP742_green b=COMP742_blue
   

   COMP742_red , COMP742_green , COMP742_blue এর লাল চ্যানেলগুলির ফলাফলের পরে d.rgbবা ব্যবহার করে একত্রিত করা যেতে পারে r.composite...

আগে এবং তার পরে নমুনা:

প্যান-শার্প করার আগে:

প্যান-শার্পিংয়ের পরে:

এই জাতীয় ছোট চিত্রগুলি থেকে দেখার সময়, তীক্ষ্ণতর পার্থক্যগুলি বলা শক্ত । তবে কী গুরুত্বপূর্ণ তা হ'ল প্যান-তীক্ষ্ণ চিত্রটির মূলটি মূল থেকে মিশ্রিতের সাথে মিলে। কার্যোদ্ধার!

অন্যান্য নোট:

• r.recode তাপ ব্যান্ডগুলি (ব্যান্ড 6) করবেন না । i.landsat.toarক্যালভিনের তাপমাত্রার মানগুলিতে এই তাপ ব্যান্ডগুলি আউটপুট দেয় (ডিএন মানগুলির সাথে কিছুই করার নয়)। রাখুন r.recodeস্বাভাবিক বহু-ভুতুড়ে এবং প্যান ব্যান্ড (ব্যান্ড 1-5,7,8) এ রুটিন।
• তাহলে আমরা কখনোই এমনকি ব্যবহার i.landsat.toar , কিন্তু এর ফলে composites সত্যই কোনও গোলোযোগ চেহারা , এটি সাধারণত কারণ আগে ও প্যান ধার প্রক্রিয়া পরে রঙ-টেবিল mismatching এর আছে। আমি r.colors RASTER color=greyপ্যান-শার্পিংয়ের আগে আসল রাস্টার ব্যান্ডগুলিতে এবং রঙগুলির ঘনিষ্ঠভাবে ম্যাচিং নিশ্চিত করতে প্যান-শার্পিংয়ের পরে ফলাফল চ্যানেলগুলিতে প্রয়োগ করেছি ।
• ভুল রঙের সংমিশ্রিত সমস্যার একটি সাধারণ ক্ষেত্রে : মূল রাস্টার ব্যান্ডগুলি থাকে color=grey255, আউটপুট i.landsat.rgbইন থাকে color=grey.eq। এতে অবাক হওয়ার মতোই না যে তারা উভয়ই আলাদা দেখাচ্ছে!
• ব্যবহারের জন্য ল্যান্ডসেট চিত্রগুলি প্রসেস করা সত্যিই একটি সময়-করের ক্রিয়াকলাপ হতে পারে । সবগুলি প্রক্রিয়াজাত হওয়ার জন্য অপেক্ষা করার সময় আরও কিছু করার দরকার হয় বা আপনি যখন থাকবেন তখন কমপক্ষে কিছু হাস্যকর পরিমাণে কফি এবং ভাল সংগীত পান;)

আশা করি এটি কারও উপকারে আসবে: ভুল কী তা খুঁজে পেতে আমাকে কয়েক দিন সময় লেগেছে ..

Comments

আপনার গবেষণা এবং প্রচেষ্টার জন্য ধন্যবাদ, আমার এসএনএন থেকে নতুন জিআরএসএস 7 এও এই সমস্যাটি ছিল। এখন পানশার্পেনড ইমেজের রঙগুলি ঠিক আছে। =)

— ভ্লাদিমির

হ্যাঁ, আমি গ্রাসের কোন সংস্করণটি চালিয়ে যাচ্ছি তা উল্লেখ করতে ভুলে গেছি - মনে করিয়ে দেওয়ার জন্য @ ভ্লাদিমিরনাওমোভকে ধন্যবাদ! আমি গ্রাস 7 এসভিএন ব্যবহার করছি; উপলব্ধি করা উচিত ছিল যে এই সমস্যাটি এমন কিছু হতে পারে যা

এটিতে ফিরে: আপনি r.mapcalcএর সংহত int()ফাংশনটি ব্যবহার করে ভাসমান পয়েন্ট রাস্টার মানচিত্রকে পূর্ণসংখ্যার রাস্টার মানচিত্রে রূপান্তর করতে পারেন । আমি মনে করি r.recodeআপনার অ্যাপ্লিকেশন প্রসঙ্গে ব্যবহার করা উচিত নয়।

— নিকস আলেকজান্দ্রিস

আমি মনে করি সামগ্রিক "সমস্যা" সম্পর্কিত i.pansharpenএটি হ'ল এটি 8-বিট রাস্টার মানচিত্রকে ইনপুট হিসাবে প্রত্যাশা করে, অর্থাত ইনপুট মানগুলি 0 থেকে 255 এর মধ্যে হওয়া উচিত

— নিকোস আলেকজান্দ্রিস

এছাড়াও, দ্রষ্টব্য যে i.pansharpenকেবল ঘাস 7_ ট্র্যাঙ্কে রয়েছে। এটি বিকাশ সংস্করণ ... :-)

— নিকস আলেকজান্দ্রিস