আর এগুলিতে দৃষ্টি আকর্ষণীয় ঘনত্বের তাপের মানচিত্র তৈরি করা


27

যদিও আমি জানি যে আরে তাপের মানচিত্র তৈরির জন্য কয়েকটি ক্রিয়াকলাপ রয়েছে, সমস্যাটি হ'ল আমি দৃষ্টি আকর্ষণীয় মানচিত্র তৈরি করতে অক্ষম। উদাহরণস্বরূপ, নীচের চিত্রগুলি উত্তাপের মানচিত্রের ভাল উদাহরণ যা আমি এড়াতে চাই। প্রথমটিতে পরিষ্কারভাবে বিশদের অভাব রয়েছে, অন্যটি (একই পয়েন্টের ভিত্তিতে) দরকারী হওয়ার জন্য খুব বিশদ। স্পটস্যাট আর প্যাকেজে ঘনত্ব () ফাংশন দ্বারা উভয় প্লট তৈরি করা হয়েছে ।

আমি কীভাবে আমার প্লটগুলিতে আরও "প্রবাহ" পেতে পারি? আমি যে লক্ষ্যটি লক্ষ্য করছি তা আরও স্পষ্টতই বাণিজ্যিক স্প্যাটিয়ালকি ( স্ক্রিনশট ) সফ্টওয়্যার তৈরি করতে সক্ষম হওয়ার ফলাফল look

কোনও ইঙ্গিত, অ্যালগরিদম, প্যাকেজ বা কোডের লাইন যা আমাকে এই দিকে নিয়ে যেতে পারে?

এখানে চিত্র বর্ণনা লিখুন


আপনি কি চিত্র () ফাংশনটি দেখেছেন? বিরতির সংখ্যা বাড়লে রেজোলিউশন বাড়ায়।
ম্যাক্রো

উত্তর:


28

দুটি জিনিস রয়েছে যা প্লটের সাবলীলতাকে প্রভাবিত করবে, আপনার কর্নেল ঘনত্বের প্রাক্কলনের জন্য ব্যান্ডউইথ এবং আপনার প্লটটিতে রঙ নির্ধারিত বিরতি।

আমার অভিজ্ঞতায়, অনুসন্ধানী বিশ্লেষণের জন্য আমি দরকারী ব্যান্ড না পাওয়া পর্যন্ত কেবল ব্যান্ডউইথকে সামঞ্জস্য করি। নীচে বিক্ষোভ।

library(spatstat)
set.seed(3)
X <- rpoispp(10)
par(mfrow = c(2,2))
plot(density(X, 1))
plot(density(X, 0.1))
plot(density(X, 0.05))
plot(density(X, 0.01))

এখানে চিত্র বর্ণনা লিখুন

কেবলমাত্র ডিফল্ট রঙের স্কিম পরিবর্তন করা কোনওরকম সাহায্য করবে না বা পিক্সেলগুলির রেজোলিউশন পরিবর্তন করবে না (যদি কোনও ডিফল্ট রেজোলিউশন খুব সুনির্দিষ্ট হয় তবে আপনার রেজোলিউশন হ্রাস করতে হবে এবং পিক্সেলগুলি আরও বড় করা উচিত)। যদিও আপনি নান্দনিক উদ্দেশ্যে ডিফল্ট রঙের স্কিম পরিবর্তন করতে চাইতে পারেন তবে এটি অত্যন্ত বৈষম্যমূলক হওয়ার উদ্দেশ্যে।

রঙটি সহায়তা করতে আপনি যে কাজগুলি করতে পারেন তা হল স্কেল স্তরটি লগারিদমে পরিবর্তন করা (যদি আপনি খুব অসাধারণ প্রক্রিয়া করেন তবেই কেবল তা সাহায্য করবে ), রঙের প্যালেটটি নীচের প্রান্তে আরও আলাদা হতে হবে (রঙ র‌্যাম্পের স্পেসিফিকেশনের ক্ষেত্রে পক্ষপাতিত্ব) আর), বা অবিচ্ছিন্ন পরিবর্তে পৃথক বিন্দুতে কিংবদন্তি সামঞ্জস্য করুন।

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

Z <- density(X, 0.1)
logZ <- eval.im(log(Z))
bias_palette <- colorRampPalette(c("blue", "magenta", "red", "yellow", "white"), bias=2, space="Lab")
norm_palette <- colorRampPalette(c("white","red"))
par(mfrow = c(2,2))
plot(Z)
plot(logZ)
plot(Z, col=bias_palette(256))
plot(Z, col=norm_palette(5))

এখানে চিত্র বর্ণনা লিখুন


শেষ চিত্রটিতে রঙগুলি স্বচ্ছ করতে (যেখানে প্রথম রঙের বিনটি সাদা) কোনও একটি কেবল রঙের র‌্যাম্প তৈরি করতে পারে এবং তারপরে আরজিবি স্পেসিফিকেশনকে স্বচ্ছ রঙের সাথে প্রতিস্থাপন করতে পারে। উপরের মত একই ডেটা ব্যবহার করে নীচের উদাহরণ।

library(spatstat)
set.seed(3)
X <- rpoispp(10)
Z <- density(X, 0.1)
A <- rpoispp(100) #points other places than density


norm_palette <- colorRampPalette(c("white","red"))
pal_opaque <- norm_palette(5)
pal_trans <- norm_palette(5)
pal_trans[1] <- "#FFFFFF00" #was originally "#FFFFFF" 

par(mfrow = c(1,3))
plot(A, Main = "Opaque Density")
plot(Z, add=T, col = pal_opaque)
plot(A, Main = "Transparent Density")
plot(Z, add=T, col = pal_trans)


pal_trans2 <- paste(pal_opaque,"50",sep = "")
plot(A, Main = "All slightly transparent")
plot(Z, add=T, col = pal_trans2)

এখানে চিত্র বর্ণনা লিখুন


আপনাকে ধন্যবাদ, খুব স্পষ্টভাবে উত্তর প্রতিক্রিয়া। এটি "স্প্যাটস্যাট.অপশনস (এনপিক্সেল = সি (100,100))" কমান্ডের সাথে মিলিত হয়ে আপনি কৌতুকটি করেছেন। আমি যদি আপনাকে আরও কিছুটা বিরক্ত করতে পারি: আপনি কী জানেন যে কীভাবে আমি "ব্যাকগ্রাউন্ড" রঙটি স্বচ্ছ করতে পারি?
ফিগারো

1
@ ফিগারো, আমি রঙগুলি আধা স্বচ্ছ করার উদাহরণ দিয়ে আপডেট করেছি।
অ্যান্ডি ডব্লিউ

ধন্যবাদ. আরেকটা জিনিস. যদিও এটি কিছুটা সামান্য বিষয় হতে পারে: আমি কীভাবে এই হিটম্যাপগুলি শেপ ফাইলে রফতানি করব? অথবা আমি কি কেবল তাদের সাথে রাস্টার হিসাবে আচরণ করব?
ফিগারো

আমার কখনই তা করা হয়নি, তাই আমি অফহন্ড বলতে পারি না। স্ট্যাকওভারফ্লোতে আর গুরুদের পক্ষে সম্ভবত একটি ভাল প্রশ্ন।
অ্যান্ডি ডব্লিউ

হিটম্যাপটিকে রাস্টার বাদে অন্য কিছু হিসাবে উপস্থাপন করা যেতে পারে?
ফিগারো

2

আকিমা প্যাকেজ থেকে আপনি ইন্টারপ ফাংশন থেকে উপকৃত হতে পারেন। এটি প্রয়োজন হলে আপনার ম্যাট্রিক্সকে অন্য রেজোলিউশনে ইন্টারপোল্ট করতে দেবে। আপনার লিঙ্কযুক্ত উদাহরণের মতো কিছু তৈরি করার জন্য আপনাকে বেশ সূক্ষ্ম গ্রিডে বিভক্ত করতে হবে (যুক্তিটি Xo এবং yo এর দৈর্ঘ্য 1000 ডলার সহ)। এটি আপনাকে একটি নতুন ম্যাট্রিক্স দেবে যা আপনি চিত্র () দিয়ে প্লট করতে পারেন। আপনি যদি স্বচ্ছতা চান তবে এটি কিছু অতিরিক্ত কাজ নেবে। রঙ প্যালেট দিয়ে এটি করা সহজ নয়, সুতরাং আপনি প্রতিটি গ্রিডকে একটি নির্ধারিত রঙের সাথে বহুভুজ হিসাবে প্লট করতে পারেন।


স্বচ্ছতা তৈরি করার অভিজ্ঞতা কি আপনার রয়েছে? উদাহরণ আছে কি? এটি পরবর্তী পদক্ষেপ হবে।
ফিগারো

ইন্টারপ () এর জন্য প্রয়োজনীয় x, y, z মানগুলি কীভাবে তৈরি করব? উদাহরণস্বরূপ, ঘনত্ব () যদি 100x100 গ্রিড ব্যবহার করে তবে প্রতিটি পয়েন্টের জন্য আমার কাছে z মান রয়েছে তবে প্রয়োজনীয় ভেক্টর বিন্যাসে আমার কাছে এক্স & ওয়াই নেই (আমার কাছে কেবল অক্ষগুলির স্থানাঙ্ক রয়েছে)।
ফিগারো

স্থানাঙ্ক জোড়া পেতে একটি লুপ তৈরি করেছেন, তবে এটির জন্য কোনও ফাংশন আছে?
ফিগারো

2
দেখে মনে হচ্ছে চিত্র সর্বোপরি স্বচ্ছ রঙ গ্রহণ করতে পারে। দেখুন: menugget.blogspot.de/2012/04/… । আপনি যদি এখনও ইন্টারপ ব্যবহার করতে চান তবে আপনি আপনার জেড স্তরের জন্য আপনার xy স্থানাঙ্কের সমন্বয়গুলি তৈরি করতে পারেন: "গ্রিড <- এক্সপ্রেসড্রিড (x = x, y = y)"
মার্ক

1

আপনি ggplot2 সন্ধান করতে পারেন। দেখে মনে হচ্ছে যে আপনি যে প্যাকেজের চেষ্টা করেছেন তাতে দুর্দান্ত রঙের স্কিম বা "প্রবাহ" নেই - আরকলারব্রেয়ারটি একবার দেখুন। এমন একটি ব্লগ রয়েছে যেখানে এটি এই প্যাকেজটিকে একটি সাধারণ উদাহরণ সহকারে বাস্তবায়িত করেছে।

আপনার লিঙ্কযুক্ত উদাহরণে যেমন দেখানো হয়েছে আপনি ভৌগলিক ডেটা প্লট করার চেষ্টা করছেন কিনা তা আমি নিশ্চিত নই তবে আপনি যদি জানেন তবে গুগল "স্ট্যাটিক ম্যাপস এপিআই ভি 2 বিকাশকারী গাইড" অফার করে এবং আপনি গুগল এবং আরকে একটি প্যাকেজ, আরজিগল্যাপস নামে একত্রিত করতে পারেন called ।

শুভ কামনা তোমার গবেষণার সাথে।


দুঃখিত যদি আমি অস্পষ্ট ছিলাম তবে এখনই সমস্যাটি রঙ স্কিম নয়। তবে আমি আজকের আগে আরকলারব্রেয়ার চেষ্টা করেছিলাম: দুর্দান্ত অল্প প্যাকেজ! এটি স্থানিক ডেটা, আমি গুগল এপিআই ব্যবহার করতে চাই তা নিশ্চিত নই, তবে আমার কিছুটা চেহারা থাকতে পারে।
ফিগারো

0

আপনি কি রেজুলেশনটি ক্র্যাঙ্ক করার চেষ্টা করেছেন density? তর্ক dimyx=c(512, 512)বা উচ্চতর চেষ্টা করুন ।

আমাদের সাইট ব্যবহার করে, আপনি স্বীকার করেছেন যে আপনি আমাদের কুকি নীতি এবং গোপনীয়তা নীতিটি পড়েছেন এবং বুঝতে পেরেছেন ।
Licensed under cc by-sa 3.0 with attribution required.