বাইরে কি অন্ধকার? একটি সূর্যের মানচিত্র আঁকুন!

Question

আমাদের নিকটতম তারা, সূর্য বেশ ফিদগেটে। এটি বেড়ে ওঠার সময়গুলি আপনি কোথায় আছেন এবং শীতকালীন কিনা তা নির্ভর করে।

আমরা যদি আমাদের বেসমেন্টের আরাম না রেখে সূর্যের বাইরে বাইরে জ্বলজ্বল করতে পারি - তবেই আমাদের একটি আপ-টু-ডেট সূর্য মানচিত্র (ওরফে দিবালোকের মানচিত্র) প্রয়োজন need আপনি সেই প্রোগ্রামটি লিখেছেন যা কেবল এটিই উত্পন্ন করে!

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

আমি যদি বলছি তারপরেও যদি আপনার কাছে কোনও ধারণা না থাকে তবে এখানে উইকিপিডিয়া থেকে একটি উদাহরণ:

এটি একটি জনপ্রিয়তা প্রতিযোগিতা । আপনি নীচের কোনটি অর্জনের চেষ্টা করছেন তা আপনার উত্তরে লক্ষ্য করা উচিত (একাধিক পছন্দ সম্ভব):

• শুদ্ধি। মনে রাখবেন যে বিধিগুলি 'একটি আনুমানিক' বলে - আপনার অনুমান যত ভাল, এই বিভাগে আরও পয়েন্ট। আপনি ওল্ফ্রাম আলফা , সময় এবং তারিখ বা ডাই ডট এর বিপরীতে আপনার বাস্তবায়ন পরীক্ষা করতে পারেন ।

• কার্যকারিতার। উদাহরণস্বরূপ, মিথস্ক্রিয়া সম্পর্কে কী? নির্দিষ্ট অবস্থান চিহ্নিত করছেন? অন্যান্য গ্রহ ম্যাপিং?

• নন্দনতত্ব। মহাদেশগুলি অঙ্কন? বোনাস পয়েন্ট. টেক্সচার মহাদেশগুলি? বোনাস পয়েন্ট. 3 ডি পৃথিবীতে? মেঘের সাথে? তারার? সঠিক তারা? প্রচুর বোনাস পয়েন্ট। ইত্যাদি।

• অস্বাভাবিক, পুরানো বা কেবল সাধারণ ভুল প্রযুক্তি ব্যবহার করা। অবশ্যই, আপনি ম্যাথমেটিকায় এইটিকে বেত্রাঘাত করতে পারেন তবে আপনি কি ব্যবহারের কথা ভেবে দেখেছেন m4? এসকিউএল? বের? x86 সমাবেশ?

• মজাদার। একটি ডায়ম্যাক্সিয়ন-প্রক্ষেপণ মানচিত্র ব্যবহার করতে চান ? এগিয়ে যান!

• শর্ট কোড. এটি কোড গল্ফ এসই, সর্বোপরি।

আনন্দ কর!

Comments

@ পিটারটেলর সম্ভবত! আমি মনে করি আপনাকে এমনকি আলোকসজ্জাও করতে হবে না: আপনি একটি 3 ডি পৃথিবী আঁকতে এবং এটিকে ঘোরতে পারেন যাতে দিবালোক অংশ (এবং অন্য কিছুই নয়) দর্শকের মুখোমুখি হয়। এটি গ্রহের রাতের সময় অংশটি দেখায় না, তবে এটির প্রয়োজন নেই।

— Wander Nauta

একটি বিকল্প সমাধান হ'ল উইন্ডোজ ইনস্টল করা। (বুনিয়াদ, আমি বলতে চাচ্ছি।)

— খুঁতখুঁতে ossifrage

@ কিউইউআর আমি m4এই কাজের জন্য সপ্তাহের যে কোনও দিন x86 সমাবেশ নেব ...

— নন্দা

এছাড়াও: "সূর্য বেশ বেঁচে থাকে কারণ এটি বিভিন্ন সময়ে উত্থিত হয় এবং অস্ত যায়।" অবশ্যই সূর্যের দোষ: P

— qwr

@ কিউইউআর কখনও কি ভেবে দেখেছেন যে আমাদের যদি 10-ঘন্টা দিন, 10-দিনের সপ্তাহ, 10-সপ্তাহ মাস এবং 10-মাস বছর এবং t = 0 এবং t = 5 এর মধ্যে সূর্য থাকে তবে সময় / তারিখ গণনা কত সহজ হবে? তবে না, সূর্যকে বিভিন্ন জায়গায় বিভিন্ন সময়ে তার কদর্য চেহারা দেখাতে হবে এবং পৃথিবীর চারপাশে যেতে খুব বেশি সময় নিতে হবে। কোনও অভিযোগ বিভাগ নেই। বিরক্তিকর।

— নওতা

Answer 1

হাস্কেল - নিম্ন মানের কোড

আমি যখন এটি লিখেছিলাম তখন আমি অত্যন্ত ক্লান্ত হয়ে পড়েছিলাম।

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

import Data.Time.Clock
import Data.Time.Calendar
import Control.Applicative
import Data.Fixed
import Data.Maybe

earth :: [[Int]]
earth = [[256],[256],[256],[256],[64,1,1,2,1,5,14,16,152],[56,19,3,27,1,6,50,1,2,1,90],[53,6,1,11,2,36,26,1,2,1,16,2,1,1,2,1,24,4,66],[47,2,5,14,4,35,22,7,54,2,1,3,60],[38,1,2,2,3,1,6,1,2,1,2,7,6,1,1,33,24,3,3,1,56,2,60],[34,2,1,4,2,1,3,1,1,3,3,2,15,3,3,29,57,5,19,1,2,11,17,1,1,1,34],[40,3,10,2,1,8,16,27,54,3,18,19,18,1,36],[33,6,5,3,2,3,1,3,2,2,1,5,16,21,1,2,53,2,10,1,6,19,1,7,4,3,9,2,33],[32,4,1,7,1,2,3,2,1,1,3,11,14,23,53,2,10,3,1,4,2,33,7,7,29],[8,5,25,10,5,3,2,14,10,2,1,18,1,2,31,6,18,1,7,4,1,60,22],[5,18,2,12,3,5,1,3,2,2,1,3,4,2,3,8,11,18,30,13,9,2,7,3,2,72,1,6,8],[4,36,2,1,1,4,3,7,1,4,3,9,8,15,34,18,2,2,2,17,1,78,4],[4,1,1,27,3,1,1,24,6,3,1,1,1,3,6,13,13,1,20,15,1,4,1,104,1],[3,31,1,24,1,2,4,8,10,9,12,6,18,7,3,7,1,1,2,99,3,2,2],[7,50,2,2,2,1,2,1,3,2,1,2,10,7,15,1,20,7,2,111,7,1],[4,35,1,15,9,1,1,3,4,1,12,5,34,8,3,110,10],[4,9,1,2,1,37,12,6,16,3,34,8,3,96,5,6,13],[6,6,1,1,8,32,12,6,3,1,49,9,4,2,1,86,1,3,4,2,19],[9,2,1,1,11,31,11,11,40,1,8,1,2,4,5,83,12,3,20],[8,1,16,33,9,11,39,2,8,1,2,3,3,83,13,5,19],[28,33,5,12,40,2,7,3,6,62,1,19,13,5,20],[27,36,2,15,34,3,2,2,6,71,1,22,11,2,22],[30,21,1,11,2,16,33,3,1,4,2,72,1,24,1,1,9,1,23],[31,21,1,26,39,4,1,98,1,1,33],[31,42,7,1,40,100,1,1,33],[33,25,2,15,4,4,35,102,36],[33,23,2,1,2,14,8,1,36,27,1,9,1,61,3,1,33],[33,26,5,14,42,10,1,11,2,2,2,7,3,5,1,9,1,44,38],[33,26,1,2,1,9,2,1,45,7,1,2,2,9,8,6,2,6,1,53,4,2,33],[33,26,1,4,1,6,44,8,6,2,3,7,9,5,3,56,1,1,4,3,33],[33,37,45,8,7,2,3,6,2,4,3,6,4,53,43],[33,36,46,6,6,1,4,1,2,2,3,16,3,47,1,5,8,2,34],[34,34,46,7,11,1,3,2,2,16,3,45,6,2,8,1,35],[34,33,48,5,11,1,4,1,4,16,2,49,3,2,6,2,35],[35,32,54,8,17,60,5,2,4,4,35],[36,30,50,12,18,60,8,2,1,1,38],[38,27,50,15,16,61,6,2,41],[38,25,51,18,3,4,6,62,6,1,42],[39,1,1,17,2,3,51,93,49],[40,1,1,11,9,2,49,31,1,10,2,50,49],[40,1,2,9,10,2,48,33,1,10,2,49,49],[41,1,2,8,11,1,47,34,2,10,5,44,50],[42,1,2,7,58,36,1,11,2,1,8,36,51],[46,6,58,36,2,15,7,34,2,1,49],[46,6,12,2,43,38,2,14,7,2,1,12,1,15,55],[46,6,5,2,7,2,41,38,2,14,10,10,4,10,59],[47,6,3,3,10,3,38,37,3,12,11,8,6,9,2,1,57],[49,10,51,38,3,9,13,7,8,9,9,2,48],[51,7,51,40,2,7,15,6,9,1,1,8,8,2,48],[55,7,47,41,1,6,17,4,12,8,8,1,49],[57,5,47,42,1,2,20,4,13,8,9,1,47],[59,3,8,1,38,43,22,4,13,1,2,4,10,2,46],[60,2,6,5,38,41,1,4,18,3,17,3,10,2,46],[61,2,1,1,2,3,1,7,34,45,18,2,18,1,60],[63,1,2,13,33,44,22,1,12,1,16,3,45],[66,14,33,43,22,1,13,1,14,1,1,1,46],[66,18,30,4,1,1,5,30,34,1,2,2,9,3,50],[66,19,43,27,34,2,2,1,7,3,52],[65,20,43,26,36,2,1,2,5,5,51],[65,21,42,24,39,3,4,7,2,1,1,1,1,1,44],[56,1,7,23,41,16,1,6,41,2,4,6,7,1,44],[64,25,39,16,1,5,42,3,4,5,2,1,8,1,2,1,37],[64,29,35,22,43,3,1,1,2,3,2,1,1,1,2,1,1,2,1,7,6,1,27],[63,31,35,20,45,2,11,1,9,7,4,2,26],[64,32,34,19,67,1,2,6,1,2,28],[65,31,34,12,1,6,48,4,18,6,31],[65,31,34,19,54,2,1,2,2,1,10,2,2,1,30],[66,29,36,14,1,3,57,1,19,2,28],[66,29,36,14,1,4,63,1,42],[67,27,36,15,1,4,63,5,3,2,33],[67,26,37,20,5,2,53,2,1,4,4,2,33],[68,25,37,20,4,3,52,9,3,3,32],[70,23,36,20,3,4,53,11,1,4,31],[71,22,37,17,5,4,51,18,31],[71,22,37,16,7,3,50,20,30],[71,21,39,15,6,3,5,1,42,24,29],[71,20,40,15,6,3,47,26,28],[71,17,43,15,6,3,46,28,27],[71,16,45,13,8,1,48,27,27],[71,16,45,12,58,28,26],[71,16,45,12,58,28,26],[70,16,47,10,59,28,26],[70,15,49,9,60,27,26],[70,14,50,7,62,7,6,13,27],[70,13,51,6,63,6,8,1,1,9,28],[70,10,138,10,28],[69,12,139,7,29],[69,11,141,5,19,3,8],[69,8,167,3,9],[69,8,166,1,1,1,10],[70,5,149,2,16,2,12],[69,6,166,3,12],[68,6,166,2,14],[68,5,166,3,14],[68,6,182],[67,6,183],[68,4,184],[68,4,6,2,176],[69,4,183],[70,5,20,1,160],[256],[256],[256],[256],[256],[256],[78,1,1,1,109,1,65],[75,2,115,1,23,1,39],[72,3,80,1,1,5,20,42,32],[74,1,70,1,4,21,5,52,2,1,25],[67,1,2,2,1,4,64,28,4,62,21],[69,9,34,1,1,1,1,1,1,1,2,48,3,69,15],[50,1,5,1,16,5,34,130,14],[32,1,1,2,4,1,3,1,4,29,32,128,18],[20,1,1,54,32,128,20],[17,49,34,137,19],[9,1,2,54,20,4,6,143,17],[16,51,18,5,10,135,21],[11,1,4,54,25,140,21],[12,66,4,155,19],[12,231,13],[0,6,9,5,2,234],[0,256],[0,256]]
main = do
    header
    mapM_ line [0..299]
    where
        header = do
            putStrLn "P3"
            putStrLn "# Some PPM readers expect a comment here"
            putStrLn "400 300"
            putStrLn "2"
        line y = mapM_ (\x -> pixel x y >>= draw) [0..399]
            where
                draw (r, g, b) = putStrLn $ (show r) ++ " " ++ (show g) ++ " " ++ (show b)
                pixel x y = fromMaybe (return (1, 1, 1)) $
                    mapRegion (\x y -> (50, -x, y)) (x - 50) (y - 50)
                    <|> mapRegion (\x y -> (-x, -50, y)) (x - 150) (y - 50)
                    <|> mapRegion (\x y -> (-x, y, 50)) (x - 150) (y - 150)
                    <|> mapRegion (\x y -> (-50, y, -x)) (x - 250) (y - 150)
                    <|> mapRegion (\x y -> (y, 50, -x)) (x - 250) (y - 250)
                    <|> mapRegion (\x y -> (y, -x, -50)) (x - 350) (y - 250)
                    where
                        mapRegion f x y = if x >= -50 && y >= -50 && x < 50 && y < 50 then
                            Just $ fmap (worldMap . shade) getCurrentTime
                            else Nothing
                                where
                                    t (x, y, z) = (atan2 y z) / pi
                                    p (x, y, z) = asin (x / (sqrt $ x*x+y*y+z*z)) / pi * 2
                                    rotate o (x, y, z) = (x, y * cos o + z * sin o, z * cos o - y * sin o)
                                    tilt o (x, y, z) = (x * cos o - y * sin o, x * sin o + y * cos o, z)
                                    shade c = ((t $ rotate yearAngle $ tilt 0.366 $ rotate (dayAngle - yearAngle) $ f x y)) `mod'` 2 > 1
                                        where
                                            dayAngle = fromIntegral (fromEnum $ utctDayTime c) / 43200000000000000 * pi + pi / 2
                                            yearAngle = (fromIntegral $ toModifiedJulianDay $ utctDay c) / 182.624 * pi + 2.5311
                                    worldMap c = case (c, index (t $ f x y) (p $ f x y)) of
                                            (False, False) -> (0, 0, 0)
                                            (False, True) -> (0, 0, 1)
                                            (True, False) -> (2, 1, 0)
                                            (True, True) -> (0, 1, 2)
                                            where
                                                index x y = index' (earth !! (floor $ (y + 1) * 63)) (floor $ (x + 1) * 127) True
                                                    where
                                                        index' [] _ p = False
                                                        index' (x:d) n p
                                                            | n < x = p
                                                            | otherwise = index' d (n - x) (not p)

এটি ঠিক - ত্রিভুজাকার- whereকোড, নেস্টেড caseএস, অবৈধ আইও ব্যবহার।

Comments

এটি বাঁকা প্রতিভা একটি কাজ। একটি পরামর্শ, কাছাকাছি fromIntegral (fromEnum $ utctDayTime c)হয় (realToFrac $ utctDayTime c)। (আমি কেবল আমার উত্তর লেখার সময় এটি শিখেছি)

— বাজরঘ

আমি সারা দিন .gif দেখতে পারি।

— মাইক দ্য লাইয়ার

mnlip, আমি @mikeTheLiar এর সাথে একমত আপনার জবাবটি আপনার উত্তরের শীর্ষে এম্বেড করা উচিত, আপনি সমস্ত ভোট পাবেন।

— বাজরাঘ

আমি সম্পূর্ণরূপে ন্যাকামি জিআইএফ-এ ভোট দিয়েছি। এটি বিশ্বের দিকে দেখার জন্য কেবল একটি বাঁকানো উপায়।

— অ্যালেন গোল্ড

কখনও তাই সামান্য

— ট্রিপ্পি

Answer 2

হাস্কেল, 'কারণ এটি আছে' বিভাগে

আমি কৌতুহল ছিল তাই আমি একটি লিখেছি। সূত্রগুলি যথাযথভাবে সঠিক [1], তবে তারপরে আমি যথাযথ প্লেট কারি মানচিত্রের পরিবর্তে কিছু আস্কি শিল্প ব্যবহার করি কারণ এটি দেখতে সুন্দর লাগে (পিক্সেলগুলিকে আমি ল্যাটে / দীর্ঘভাবে রূপান্তর করি যেভাবে কেবল প্লেট কারির জন্য সঠিকভাবে কাজ করে)

import Data.Time
d=pi/180
tau=2*pi
m0=UTCTime(fromGregorian 2000 1 1)(secondsToDiffTime(12*60*60))
dark lat long now =
  let
    time=(realToFrac$diffUTCTime now m0)/(60*60*24)
    hour=(realToFrac$utctDayTime now)/(60*60)
    mnlong=280.460+0.9856474*time
    mnanom=(357.528+0.9856003*time)*d
    eclong=(mnlong+1.915*sin(mnanom)+0.020*sin(2*mnanom))*d
    oblqec=(23.439-0.0000004*time)*d
    ra=let num=cos(oblqec)*sin(eclong)
           den=cos(eclong) in
       if den<0 then atan(num/den)+pi else atan(num/den)
    dec=asin(sin(oblqec)*sin(eclong))
    gmst =6.697375+0.0657098242*time+hour
    lmst=(gmst*15*d)+long
    ha=(lmst-ra)
    el=asin(sin(dec)*sin(lat)+cos(dec)*cos(lat)*cos(ha))
  in
  el<=0

td x = fromIntegral x :: Double
keep="NSEW"++['0'..'9']
pixel p dk=if dk && p`notElem`keep then if p==' ' then '#' else '%' else p
showMap t= do
  let w=length(worldmap!!0)
      h=length worldmap
  putStrLn (worldmap!!0)
  putStrLn (worldmap!!1)
  mapM_(\y->do
           mapM_(\x->let
                    lat=(0.5-td y/td h)*pi
                    long=(0.5-td x/td w)*tau
                    in
                     putStr [pixel ((worldmap!!(y+2))!!x) (dark lat long t)]) [0..(w-1)]
           putStrLn "") [0..(h-4)]
  putStrLn (last worldmap)

main = do {t<-getCurrentTime; showMap t}

worldmap=[
 "180 150W  120W  90W   60W   30W  000   30E   60E   90E   120E  150E 180",
 "|    |     |     |     |     |    |     |     |     |     |     |     |",
 "+90N-+-----+-----+-----+-----+----+-----+-----+-----+-----+-----+-----+",
 "|          . _..::__:  ,-\"-\"._       |7       ,     _,.__             |",
 "|  _.___ _ _<_>`!(._`.`-.    /        _._     `_ ,_/  '  '-._.---.-.__|",
 "|.{     \" \" `-==,',._\\{  \\  / {)     / _ \">_,-' `                mt-2_|",
 "+ \\_.:--.       `._ )`^-. \"'      , [_/(                       __,/-' +",
 "|'\"'     \\         \"    _L       oD_,--'                )     /. (|   |",
 "|         |           ,'         _)_.\\\\._<> 6              _,' /  '   |",
 "|         `.         /          [_/_'` `\"(                <'}  )      |",
 "+30N       \\\\    .-. )          /   `-'\"..' `:._          _)  '       +",
 "|   `        \\  (  `(          /         `:\\  > \\  ,-^.  /' '         |",
 "|             `._,   \"\"        |           \\`'   \\|   ?_)  {\\         |",
 "|                `=.---.       `._._       ,'     \"`  |' ,- '.        |",
 "+000               |    `-._        |     /          `:`<_|h--._      +",
 "|                  (        >       .     | ,          `=.__.`-'\\     |",
 "|                   `.     /        |     |{|              ,-.,\\     .|",
 "|                    |   ,'          \\   / `'            ,\"     \\     |",
 "+30S                 |  /             |_'                |  __  /     +",
 "|                    | |                                 '-'  `-'   \\.|",
 "|                    |/                                        \"    / |",
 "|                    \\.                                            '  |",
 "+60S                                                                  +",
 "|                     ,/           ______._.--._ _..---.---------._   |",
 "|    ,-----\"-..?----_/ )      _,-'\"             \"                  (  |",
 "|.._(                  `-----'                                      `-|",
 "+90S-+-----+-----+-----+-----+----+-----+-----+-----+-----+-----+-----+",
 "Map 1998 Matthew Thomas. Freely usable as long as this line is included"]

বছরের আউটপুট উদাহরণস্বরূপ, বছরের আরও আকর্ষণীয় সময় থেকে (আমরা বিষুব কাছাকাছি, তাই ওয়ান্ডার নওতার আয়তক্ষেত্রাকার ব্লবগুলি মোটামুটি নির্ভুল :)) - এটি 16 জানুয়ারী 13:55:51 ইউটিসি 2014 এর জন্য:

180 150W  120W  90W   60W   30W  000   30E   60E   90E   120E  150E 180
|    |     |     |     |     |    |     |     |     |     |     |     |
%90N%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
%##########%#%%%%%%%%##%%%%%%%#######%7#######%#####%%%%%#############%
%##%%%%%#%#%%%%%%%%%%%%%%####%########%%%#####%%#%%%##%##%%%%%%%%%%%%%%
%%%#####%#%#%%%%%%%%%%%##%##%#%%#####%#%#%%%%%%#%################%%%2%%
%#%%%%%%%#######%%%#%%%%%#%%######, [_/(         ##############%%%%%%#%
%%%%#####%#########%####%%#####  oD_,--'            ####%#####%%#%%###%
%#########%###########%%#####    _)_.\\._<> 6        ######%%%#%##%###%
%#########%%#########%######    [_/_'` `"(             ###%%%##%######%
%30N#######%%####%%%#%#####     /   `-'"..' `:._       ###%%##%#######%
%###%########%##%##%%#####     /         `:\  > \  ,-^. #%%#%#########%
%#############%%%%###%%###     |           \`'   \|   ?_)##%%#########%
%################%%%%%%%#      `._._       ,'     "`  |' %%#%%########%
%000###############%####`-._        |     /          `:`<_%%%%%%######%
%##################%####    >       .     | ,          `=.%%%%%%%#####%
%###################%%#    /        |     |{|              %%%%%#####%%
%####################%#  ,'          \   / `'            ,"#####%#####%
%30S#################%  /             |_'                |  %%##%#####%
%####################% |                                 '-'##%%%###%%%
%####################|/                                      ##%####%#%
%####################\.                                       #####%##%
%60S################                                          ########%
%##################   ,/           ______._.--._ _..---.-------%%%%###%
%####%%%%%%%%%%%%%--_/ )      _,-'"             "                ##%##%
%%%%%###########       `-----'                                    ##%%%
%90S%%%%%%%%%----+-----+-----+----+-----+-----+-----+-----+-----+----%%
Map 1998 Matthew Thomas. Freely usable as long as this line is included

[1] 0 থেকে 360 এর মধ্যে ডিগ্রি রাখা, 0 থেকে 24 এর মধ্যে ঘন্টা এবং 0 এবং 2pi এর মধ্যে রেডিয়েন্স বাদে আপনি অন্য কোথাও এটির মতোই পাবেন। আমি মনে করি যে দিনগুলি আমরা স্লাইড বিধি ব্যবহার করেছি সেগুলি হোল্ডওভারগুলি; ট্রিগ ফাংশনগুলি এই সীমার বাইরে ঠিক কাজ করে ...

Comments

উজ্জ্বল! আমি পছন্দ করি যে আপনি এখনও 'অন্ধকার' এর মাধ্যমে মানচিত্রটি দেখতে পাচ্ছেন। এছাড়াও, গণিতটি শক্ত দেখায়। আপনি উদাহরণটির জন্য আপনি যে তারিখটি ব্যবহার করেছেন তা যুক্ত করতে পারেন, যাতে অন্য লোকেরা কী তাদের সমাধানগুলি আপনার সাথে তুলনা করতে পারে?

— নাউটা

আমি দেখতে পেলাম যে আপনি তারিখটি যুক্ত করেছেন, ধন্যবাদ!

— নটটা

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

— বাজরাঘ

Answer 3

।

বাশ, 882 * অক্ষর

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

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

রামের প্রবেশ এবং আমার প্রথম প্রবেশের মতো এটিও seasonতু পরিবর্তনের জন্য নয়, যার অর্থ এটি খুব সুনির্দিষ্ট নয়। তবে এটি আমার প্রথম প্রবেশের চেয়ে খাটো, সুন্দর এবং আরও সুনির্দিষ্ট - এবং রামের প্রবেশের বিপরীতে মানচিত্রের ডেটা এবং জিআইএফ অ্যানিমেশন তৈরির কোড অন্তর্ভুক্ত রয়েছে।

                               echo '
                    iVBO  Rw0KGgoAAAA       NS
              UhE  U g      AAAEgAAAA                     kAQMAAAAQFe4lAAAABlB
    MVEUAFFwAbxKgAD63 AAAA   AWJLR0                  QAiAUdSAAAAAlwSFlzAAALEwAACx
 MB AJqcGAAAAAd0SU1FB9  4DE  hUWI   op      Fp5MAAADDSURBVBhXrcYhTsNQGADgr3ShE4Qi
    h4BeYQFBgqAJN8Lh    +r                jBb rArIJHPobgAgkzgeSQkVHT7MWThAHzq44
           /j/jezy6jSH  M6fB           gd  9T Nbxdl99R4Q+XpdNRISj4dlFRCz
            oI11FxIpup4uIRDe5           fokp0Y2W25jQFDfrGNGsDNsoqBaGj34D2
             bA7TcAwnmRoDZM             5tLkePUJb6uIT2rEq7hKaUhUHCXWpv7Q
             PqEv1rsuoc7X              RbV Bn2d   kGTYKMQ3C7H8z2+wc/eMd S
              QW39v8kAAA               AA      SUVOR K5CYII='|base64 \
               -di>t;for                X in     {0..23};do R=$((90-(\
                $X*15)                )); echo "camera{location <0,
                 0,                   -5> angle 38 }    light_source{
                  <0,0,               -1000> rgb < 2,2,   2>} sphere
                    {<0              ,0,0> 1 pigment      {
                      /**/            image_map{\"t\"        map_type
                        1}}                rotate           <0,$R,0>
                        }">s               ;povray             +Is +H300\
                        +Of$X.png          +W400
                        mogrify            -fill                     white    \
                        -annotate           +0+10                    "$X:00" \
                         -gravity           south                    f$X.png
                         done;              convert                -delay     \
                         100                -loop                 0 $(ls f*  \
                         |sort               -V)                  ani.gif
                        exit;

বোনাস হিসাবে , এখানে এমন একটি জিআইএফ রয়েছে যা স্পেস-সেভিং 1-বিট টেক্সচারের পরিবর্তে নাসার ব্লু মার্বেল চিত্রটি ব্যবহার করে, অর্থাত কোনও আকারের সীমাবদ্ধতা ছাড়াই ফলাফলটি কেমন হবে: http://i.imgur.com/AnahEIu.gif

^{*: 882 টি অক্ষর সজ্জিত শ্বেত স্পেস গণনা করছে না, মোট 1872 টি অক্ষর।}

Comments

এগুলি সমস্ত স্বনির্ভর করার জন্য +1। এবং স্ব-রেফারেন্টাল কোড তৈরির জন্য যা নিজেই একটি বিশ্বের মানচিত্রের মতো লাগে। চমৎকার কাজ.

— রাম নরসিমহান

আহা! অনেকগুলো নিবন্ধন. পাগল না করার জন্য এখন আমি আমার অজুহাত হারাচ্ছি ...

— বাজরঘ

অসুস্থ! ভাল লাগবে।

— পান্ডুবিয়ার

হেই, আমি এখন বুঝতে পেরেছি যে আমি বিশ্বের মানচিত্রের জন্য উত্স হিসাবে কোডের বিন্যাসটি নিজেই ব্যবহার করতে পারতাম (স্পেস = সমুদ্র, অন্য সবকিছু = ভূমি সহ) এবং আসলে কম অক্ষরের সাথে আরও ভাল রেজোলিউশন পেতে পারি । ওহ ভাল ...

— Wander Nauta

এটি মাইনক্রাফ্টের মতো দেখাচ্ছে।

— কাজ ওল্ফ

Answer 4

আমি সংক্ষিপ্ত কোড বিভাগে আমার নিজের একটি এন্ট্রি নিয়ে প্রতিযোগিতাটি শুরু করার সিদ্ধান্ত নিয়েছি । এটি 923 টি অক্ষর দীর্ঘ, নিউলাইনগুলি গণনা করছে না।

সি: 923 টি অক্ষর

কোডটি এখানে:

i;j;w=160;f=40;t;b;p;s;e;k;d=86400;q=599;
char* m="M('+z EDz :!#\"!*!8S$[\"!$!#\"\")\"!3R)V$'!!()1M./!F)\"!!!!)'/GE5@\"\"!&%.3&,Y$D\"!!%$)5i\"\"\"F\"%&&6%!e'A#!#!!#&$5&!f&A'$*\"5&!c-#'3''8\"$!!#\"U'\"=5$'8#$$\"S(#=7!*5\"!\"#['!A@6#!^H=!#6bH;!!!\"6_!!I;<&!&\"!!$\"F\"!I8;&\"#\"$&#\"C#\"I7<%#!\"/\"BP5=$*,\"=#\"$!L4A%&\"\"G\"\"\"#M1@)*F\"%P/@,!N#!S(E;!@W'E=!!!<Y&D7!&!\"$7\\$D8!)$4_$C8!('&#&!!a&@9!&(%$&g$>9!$*#(%h\">:!!-\"(%&!b!$&5:!\"+\"(!!#$!!!c+5<-!'!'!#!e)5:.!(!&!\"\"e,:25!!!\"!\"\"h-;07#\"$h.9/:\"\"$!!#\"a17-;'!\"$!!\"$!X46,<\"%\"&$\\45,>#&!$$#!W45,C!!!'!\"!$!V26,H\"#!$!\"!\"!S17-#!A!!#\"!_07,\"#A&!\"`.7+#\"A*.!Q.7*$\">/^-9)$\"=0^*<)$!>1]*<(D1])>&E2\\)>&F&!)\\)@#G$%(\\'w%]'x#,\"P%z .\"P%z .!R$z -\"S$z b#z c#z d#z 3";
main(){
t=(time(0)%d*160)/d;
printf("P2\n%d 62\n5\n",w);
for(;i<q;i++){
for(j=m[i]-' ';j>0;j--){
p=k%w,s=(t-f),e=(t+f);
printf("%c ","1324"[b*2+((p>s&&p<e)||(p>s+w&&p<e+w)||(p>s-w&&p<e-w))]);
k++;
}
b=!b;
}
}

এখানে কিভাবে এটা কাজ করে:

বিশ্বের ^* একটি অশোধিত বিটম্যাপটি স্ট্রিং হিসাবে এনকোডযুক্ত রান-দৈর্ঘ্যের। স্ট্রিংয়ের প্রতিটি অক্ষর স্থল বা সমুদ্রের পিক্সেলগুলির একটি রানকে উপস্থাপন করে। দীর্ঘমেয়াদে সমুদ্রকে এক রান সমুদ্রের মধ্যে বিভক্ত করা হয়, তারপরে 0 টি ল্যান্ড পিক্সেল, তারপরে সমুদ্রের আর একটি রান, স্ট্রিংয়ের অক্ষর অক্ষরগুলি এড়ানোর জন্য। পাইথন স্ক্রিপ্ট আমি এই বিন্যাসে PBM ফাইল রূপান্তর করতে লিখেছেন এখানে ।

আমি তখন সময় () 1 জানুয়ারী 1 জানুয়ারী মধ্যরাত, গ্রীনউইচ মধ্যে কত সেকেন্ড পেরিয়ে গেছে তা জানতে ব্যবহার করি I আমি তথ্যের সাহায্যে মানচিত্রের হালকা অংশটি স্থিত করার জন্য আজ সেখানে কত সেকেন্ড কেটে গেছে তা জানতে মডুলো করেছি- বা-সেই অনুযায়ী কম (আমি আশা করি)।

সঠিকতা একটি রসিকতা। কোনও গণিত নেই। কোডটি ধরে নিয়েছে যে পৃথিবীটি একটি সিলিন্ডার (ব্লক-আকারের দিন / রাত), যে সূর্য সরাসরি নিরক্ষীয় অঞ্চলে (গ্রীষ্ম / শীতকালীন নয়) এর উপরে থাকে এবং আপনি বর্ণ ধূসর (কোন বর্ণের) পছন্দ করেন না।

প্লাস দিকে, আমি মহাদেশগুলি আঁকছি।

আউটপুটটি পোর্টেবল গ্রেম্যাপ (পিজিএম) ফর্ম্যাটে রয়েছে, যা পরে ইমেজম্যাগিক বা জিআইএমপি-র মাধ্যমে পিএনজিতে রূপান্তরিত হতে পারে।

এখানে একটি উদাহরণ আউটপুট, পিএনজিতে রূপান্তরিত ( বৃহত সংস্করণ ):

^{*: অ্যান্টার্কটিকা ব্যতীত পুরো বিশ্ব, তবে যিনি যাইহোক সেখানে বাস করেন ...}

Comments

সুন্দর আউটপুট, অন্ধকার অঞ্চলটি বাঁকা থাকলে আরও ভাল হবে

— qwr

হাঁ! আউটপুটটি বাঁকানো হওয়ার অর্থ আপনাকে কিছু ত্রিকোণমিতি করতে হবে যদিও এটি এটি আরও দীর্ঘ করবে। (অথবা আমি অনুমান করি যে আপনি এটি সঠিক-ইশ দেখানোর জন্য সঠিক কোণটি গোল করতে পারেন ...)

— নটটা

@ ওয়ান্ডারনাউটা এমন কিছু পেঙ্গুইনকে ভীষণভাবে হতাশ করবে যে তারা তাদের প্রতারণা করবে যে আপনি তাদের বিবেচনায়

— নেন নি

@ ওলিওয়েস্ট যদি আপনি আপনার পেঙ্গুইনের মতো পোলগুলিতে বাস করেন তবে আপনার এই প্রোগ্রামটির দরকার নেই - মধ্যরাতের সূর্য এবং সমস্ত।

— নটটা

Answer 5

হাস্কেল - এটি হ্যামার সময়।

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

চিত্রটিকে একটি উচ্চতর রেজোলিউশনে পরিবর্তিত করতে সম্পাদনা করা হয়েছে তবে (ইচ্ছাকৃতভাবে) মোটা অন্তর্নিহিত বিটম্যাপ, যাতে আপনি সর্পিলটির প্রভাব দেখতে পারেন। এই 5000 পয়েন্ট (থেকে নমুনা ব্যবহার ~ 260000 ), একটি 50x100 বিটম্যাপ সমতূল্য, কিন্তু খুঁটি চেয়ে বিষুবরেখা আরো রেজল্যুশন দেয়।

কোডটি ব্যবহার করতে, ghc দিয়ে সংকলন করুন, একটি anচ্ছিক সংখ্যাসূচক প্যারামিটার দিয়ে চালান যা ঘন্টা অফসেট; ফাইলগুলি 'আর্থ0.pgm', 'আর্থ1.pgm' এর মতো উত্পন্ন হয়।

import System.Environment
import Data.List (intercalate,unfoldr)
import qualified Data.Set as Set
import Data.List.Split
import Data.List
import Data.Maybe (catMaybes)
import qualified Data.Map as Map
import Data.Time
import Debug.Trace
d=pi/180
tau=2*pi
m0=UTCTime(fromGregorian 2000 1 1)(secondsToDiffTime(12*60*60))
dark::Double->Double->UTCTime->Bool
dark lat long now =
  let
    time=(realToFrac$diffUTCTime now m0)/(60*60*24)
    hour=(realToFrac$utctDayTime now)/(60*60)
    mnlong=280.460+0.9856474*time
    mnanom=(357.528+0.9856003*time)*d
    eclong=(mnlong+1.915*sin(mnanom)+0.020*sin(2*mnanom))*d
    oblqec=(23.439-0.0000004*time)*d
    ra=let num=cos(oblqec)*sin(eclong)
           den=cos(eclong) in
       if den<0 then atan(num/den)+pi else atan(num/den)
    dec=asin(sin(oblqec)*sin(eclong))
    gmst =6.697375+0.0657098242*time+hour
    lmst=(gmst*15*d)+long
    ha=(lmst-ra)
    el=asin(sin(dec)*sin(lat)+cos(dec)*cos(lat)*cos(ha))
  in
  el<=0
infill(open, known)= 
  if null open then known else infill gen
  where
    neighbours (x,y)=catMaybes $ map ((flip Map.lookup) known) [(x+1,y),(x-1,y),(x,y+1),(x,y-1),(x+1,y+1),(x-1,y+1),(x-1,y-1),(x-1,y-1)] 
    vote a= if null a then Nothing
             else Just ((sum a)`div`(length a))
    fill x (open',  known')=
      case vote (neighbours x) of
        Nothing->(x:open',known')
        Just c->(open',(x,c):known')
    gen=(\(o,k)->(o,Map.fromList k))$foldr fill ([], Map.toList known) open
mpoint (a,b)=case a of Nothing->Nothing;Just c->Just(c,b)
grid w h n g lut= map (\y->map (\x->if Set.member (x,y) g then 3 else case Map.lookup (x,y) lut of Nothing->7;Just c->c) [1..w]) [1..h]
unknowns w h lut=concatMap (\y->concatMap (\x->let z=1-(2*x//w-1)^2-(2*y//h-1)^2 in case Map.lookup (x,y) lut of Nothing->if z<0 then [] else [(x,y)];_->[]) [1..w]) [1..h]
main=do
  args <- getArgs
  let off = if null args then 0 else read(args!!0)
  actual <- getCurrentTime
  let now=((fromIntegral off)*60*60) `addUTCTime` actual
  let tod=realToFrac(utctDayTime now)/86400+0.4
  let s=5000
  let w=800
  let h=400
  let n=6
  -- pbm <- readFile "earth.pbm"
  -- let bits=ungrid s$parsepbm pbm
  let bits=[0,23,4,9,1,3,1,2,6,10,1,10,4,1,3,7,10,7,4,2,2,1,2,6,12,1,1,2,1,5,4,1,8,1,3,
            1,21,7,2,2,35,1,4,3,2,2,2,2,16,1,25,1,2,8,1,4,1,2,13,3,2,1,26,1,1,10,3,3,8,
            2,3,6,1,3,25,2,1,10,15,5,1,6,2,3,30,10,15,19,32,11,16,20,35,11,1,2,14,22,27,
            1,8,14,16,22,2,1,22,1,1,2,1,1,2,1,2,1,3,16,14,25,1,2,21,1,6,1,2,1,1,2,3,17,
            14,26,1,2,1,1,26,1,1,3,3,1,1,19,13,28,4,1,26,6,6,21,11,35,40,21,11,37,41,20,
            2,4,4,1,1,39,19,1,6,1,16,19,2,4,5,40,18,2,7,1,17,19,1,1,1,1,1,2,3,46,7,1,5,
            4,25,16,3,1,1,3,5,44,1,4,5,4,3,6,4,1,19,22,5,46,2,3,4,6,2,9,22,22,2,50,1,5,
            2,1,1,6,1,8,24,15,5,1,2,51,2,5,1,1,1,5,1,10,23,14,9,55,1,4,2,17,16,1,4,14,9,
            57,4,1,3,17,13,20,11,54,2,1,3,1,2,20,12,18,13,47,4,3,8,21,10,17,15,44,5,1,1,
            4,1,3,2,22,10,15,16,46,4,3,1,2,2,25,9,17,15,47,1,1,3,30,9,18,13,46,2,1,4,25,
            2,1,11,16,13,46,8,24,2,2,9,16,11,45,12,22,1,3,7,17,10,45,12,21,1,3,7,19,8,
            43,12,25,6,19,8,41,12,25,5,20,7,40,11,25,4,20,6,40,5,3,2,48,6,38,3,54,4,30,
            1,6,2,55,2,29,1,5,1,53,3,28,1,55,3,49,1,30,2,76,1,284,3,4,1,15,1,17,10,1,9,
            7,1,13,21,4,4,1,2,6,17,2,8,3,63]
  let t(phi,lambda)=unitsphere$rx (-pi/4)$rz (-tod*4*pi)$sphereunit(phi, lambda)
  let hmr=(fmap (\(x,y)->(floor((fl w)*(x+4)/8),floor((fl h)*(y+2)/4)))).hammer.t
  let g=graticule hmr n
  let lut = Map.fromList$ catMaybes $map mpoint$map (\((lat,long),bit)->(hmr(lat,long),bit*4+2-if dark lat long now then 2 else 0))  $zip (spiral s) (rld bits)
  -- let lut = Map.fromList$ catMaybes $map mpoint$map (\((lat,long),bit)->(hmr(lat,long),bit))$zip (spiral s) (rld bits)
  let lut' = infill ((unknowns w h lut), lut)
  let pgm = "P2\n"++((show w)++" "++(show h)++" 7\n")++(intercalate "\n" $ map (intercalate " ")$chunksOf 35 $ map show(concat$grid w h n g lut'))++"\n"
  writeFile ("earth"++(show off)++".pgm") pgm

fl=fromIntegral
spiral::Int->[(Double,Double)]
spiral n=map (\k-> let phi=acos(((2*(fl k))-1)/(fl n)-1) in rerange(pi/2-phi,sqrt((fl n)*pi)*phi)) [1..n]
rld::[Int]->[Int]
rld bits=concat$rld' (head bits) (tail bits)
  where
   rld' bit []=[]
   rld' bit (run:xs) = (replicate run bit):(rld' (case bit of 1->0;_->1) xs)
rle::[Int]->[Int]
rle bits=(head bits):(map length$group bits)
sample::Int->Int->Int->[(Int,Int)]
sample n w h = map (\(phi, theta)->((floor((fl h)*((phi-(pi/2))/pi)))`mod`h, (floor((fl w)*(theta-pi)/(tau)))`mod`w )) $ spiral n
ungrid::Int->[[Int]]->[Int]
ungrid n g = rle $ map (\(y, x)->(g!!y)!!x) (sample n w h)
  where w = length$head g
        h = length g
parsepbm::[Char]->[[Int]]
parsepbm pbm=
    let header = lines pbm
        format = head header
        [width, height] = map read$words (head$drop 1 header)
        rest = drop 2 header
        d = ((map read).concat.(map words)) rest
    in chunksOf width d
rerange(phi,lambda)
 | abs(phi)>pi = rerange(phi - signum(phi)*tau, lambda)
 | abs(phi)>pi/2 = rerange(phi-signum(phi)*pi, lambda+pi)
 | abs(lambda)>pi = rerange(phi, lambda - signum(lambda)*tau)
 | otherwise = (phi, lambda)
laea(phi,lambda)=if isInfinite(z) then Nothing else Just (z*cos(phi)*sin(lambda),z*sin(phi)) where z=4/sqrt(1+cos(phi)*cos(lambda))
hammer(phi,lambda)=case laea(phi, lambda/2) of Nothing->Nothing; Just(x,y)->Just (x, y/2)
bresenham :: (Int, Int)->(Int, Int)->[(Int, Int)]
bresenham p0@(x0,y0) p1@(x1,y1)
  | abs(dx)>50||abs(dy)>50=[]
  | x0>x1 = map h$ bresenham (h p0) (h p1)
  | y0>y1 = map v$ bresenham (v p0) (v p1)
  | (x1-x0) < (y1-y0) = map f$ bresenham (f p0) (f p1)
  | otherwise = unfoldr (\(x,y,d)->if x>x1 then Nothing else Just((x,y),(if 2*(d+dy)<dx then(x+1,y,d+dy)else(x+1,y+1,d+dy-dx)))) (x0,y0,0)
      where 
        h(x,y)=(-x,y)
        v(x,y)=(x,-y)
        f(x,y)=(y,x)
        dx=x1-x0
        dy=y1-y0
globe n k= 
  (concatMap (\m->map (meridian m) [k*(1-n)..k*(n-1)]) [k*(1-2*n),k*(2-2*n)..k*2*n])
  ++(concatMap (\p->map (parallel p) [k*(-2*n)..k*2*n]) [k*(1-n),k*(2-n)..k*(n-1)])
  where
  meridian m p=(radians(p,m),radians(p+1,m))
  parallel p m=(radians(p,m),radians(p,m+1))
  radians(p,m)=rerange((p//(k*n))*pi/2,(m//(k*n))*pi/2)
graticule f n=Set.fromList $ concatMap (\(a,b)->case (f a,f b) of (Nothing,_)->[];(_,Nothing)->[];(Just c,Just d)->bresenham c d) (globe n 4)
rx theta (x,y,z) = (x, y*(cos theta)-z*(sin theta), y*(sin theta)+z*(cos theta))
ry theta (x,y,z) = (z*(sin theta)+x*(cos theta), y, z*(cos theta)-x*(sin theta))
rz theta (x,y,z) = (x*(cos theta)-y*(sin theta), x*(sin theta)+y*(cos theta), z)
sphereunit (phi, theta) = (rz theta (ry (-phi) (1,0,0)))
unitsphere (x,y,z) = (asin z, atan2 y x)
x//y=(fromIntegral x)/(fromIntegral y)    

Comments

এটা পাগলামি. পরম উন্মাদনা। আমি এটা ভালোবাসি.

— নুটা

এইটির ডানদিকে রূপান্তর পেতে আমার এত দীর্ঘ সময় লেগেছে। একাকী সর্পিল বিটম্যাপটি হাসিখুশি।

— বাজরাঘ

পুনরায় "এটি যথেষ্ট পরিমাণে স্কেল আপ করে" - আপনি কি সহজেই উচ্চতর-রেজোলিউশন জিআইএফ তৈরি করতে পারবেন তার অর্থ?

— ঘোরাফেরা নওতা

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

— বাজরঘ

সেখানে আপনি যান - ক্রেজিটাউনের দুর্দান্ত বড় চিত্র

— বাজরঘ

Answer 6

সি, পিএনএম চিত্র ব্যবহার করে

দেরী উত্তর, সঠিকতা এবং নান্দনিকতার উপর ফোকাস । আউটপুট হল দুটি ইনপুট ইমেজ (day.pnm এবং night.pnm) এর সংমিশ্রণ, গোধূলি একটি স্ট্রাইপ সহ। আমি এখানে নাসার নীল মার্বেলের উপর ভিত্তি করে চিত্রগুলি ব্যবহার করছি।

কোডটি স্পষ্টতার জন্য আমার নিজের img.h ব্যবহার করে (কেবলমাত্র কঠোর নিয়ম মেনে চলার জন্য .c তে ভারব্যাটিম অন্তর্ভুক্ত করা হয়েছে তা কল্পনা করুন ...)। সেখানকার সমস্ত কিছুই সি ম্যাক্রোর মাধ্যমে প্রয়োগ করা হয়। অ্যানিমেশনগুলি একাধিক ফ্রেম থেকে চিত্রগ্রাহকদের রূপান্তর করে তৈরি করা হয় - প্রোগ্রামটি নিজেই কেবল স্থির চিত্রগুলি আউটপুট দেয়। কোডটি নীচে।

এখন: (13 আগস্ট, ~ 13: 00 সিইএসটি)

একদিন: (১ লা জানুয়ারি)

এক বছর: (12:00 ইউটিসি)

sun.c

  #include <math.h>
  #include <time.h>

  #include "img.h"

  #ifndef M_PI
  #define M_PI 3.14159265359
  #endif

  double deg2rad(double x) {return x / 180.0 * M_PI;}
  double rad2deg(double x) {return x * 180.0 / M_PI;}

  double  sind(double x) {return  sin(deg2rad(x));}
  double  cosd(double x) {return  cos(deg2rad(x));}
  double asind(double x) {return rad2deg(asin(x));}

  double elevation(double latitude, double longitude, int yday, int hour, int min, int sec)
  {
     double fd = (hour + (min + sec / 60.0) / 60.0) / 24.0;
     double fyd = 360.0 * (yday + fd) / 366.0;

     double m = fyd - 3.943;
     double ta = -1.914 * sind(m) + 2.468 * sind(2 * m + 205.6);
     double hourangle = (fd - 0.5) * 360.0 + longitude + ta;
     double decl = 0.396 - 22.913 * cosd(fyd) + 4.025 * sind(fyd) - 0.387 * cosd(2 * fyd) + 0.052 * sind(2 * fyd) - 0.155 * cosd(3 * fyd) + 0.085 * sind(3 * fyd);

     return asind(cosd(hourangle) * cosd(decl) * cosd(latitude) + sind(decl) * sind(latitude));
  }

  int main(int argc, char* argv[])
  {
     Image day, night, out;
     int x, y;
     time_t t = time(0);
     struct tm* utc = gmtime(&t);
     int yday = utc->tm_yday, hour = utc->tm_hour, min = utc->tm_min, sec = utc->tm_sec;

     imgLoad(day, "day.pnm");
     imgLoad(night, "night.pnm");
     imgLoad(out, "day.pnm");
     for(y = 0; y < day.height; ++y)
     {
        double latitude = 90.0 - 180.0 * (y + 0.5) / day.height;
        for(x = 0; x < day.width; ++x)
        {
           double longitude = -180.0 + 360.0 * (x + 0.5) / day.width;
           double elev = elevation(latitude, longitude, yday, hour, min, sec);
           double nf = elev > -0.8 ? 0.0 : elev > -6.0 ? 0.5 : 1.0;
           double df = 1.0 - nf;
           Color dc = imgGetColor(day, x, y);
           Color nc = imgGetColor(night, x, y);
           imgDotC3(out, x, y, df * dc.r + nf * nc.r, df * dc.g + nf * nc.g, df * dc.b + nf * nc.b);
        }
     }
     imgSave(out, "out.pnm");
  }

img.h

  #include <stdlib.h>
  #include <stdio.h>
  #include <string.h>

  typedef struct
  {
     unsigned char r;
     unsigned char g;
     unsigned char b;
  } Color;

  typedef struct
  {
     Color* data;
     int width;
     int height;
     Color c;
  } Image;

  #define imgCreate(img, w, h)           {\
                                            int length;\
                                            (img).width = (w);\
                                            (img).height = (h);\
                                            length = (img).width * (img).height * sizeof(Color);\
                                            (img).data = malloc(length);\
                                            memset((img).data, 0, length);\
                                            (img).c.r = (img).c.g = (img).c.b = 0;\
                                         }

  #define imgDestroy(img)                {\
                                            free((img).data);\
                                            (img).width = 0;\
                                            (img).height = 0;\
                                            (img).c.r = (img).c.g = (img).c.b = 0;\
                                         }

  #define imgSetColor(img, ur, ug, ub)   {\
                                            (img).c.r = (ur);\
                                            (img).c.g = (ug);\
                                            (img).c.b = (ub);\
                                         }

  #define imgDot(img, x, y)              {\
                                            (img).data[(int)(x) + (int)(y) * (img).width] = (img).c;\
                                         }

  #define imgDotC3(img, x, y, ur, ug, ub) {\
                                            (img).data[(int)(x) + (int)(y) * (img).width].r = (ur);\
                                            (img).data[(int)(x) + (int)(y) * (img).width].g = (ug);\
                                            (img).data[(int)(x) + (int)(y) * (img).width].b = (ub);\
                                         }

  #define imgDotC(img, x, y, c)          {\
                                            (img).data[(int)(x) + (int)(y) * (img).width] = (c);\
                                         }

  #define imgGetColor(img, x, y)         ((img).data[(int)(x) + (int)(y) * (img).width])

  #define imgLine(img, x, y, xx, yy)     {\
                                            int x0 = (x), y0 = (y), x1 = (xx), y1 = (yy);\
                                            int dx =  abs(x1 - x0), sx = x0 < x1 ? 1 : -1;\
                                            int dy = -abs(y1 - y0), sy = y0 < y1 ? 1 : -1;\
                                            int err = dx + dy, e2;\
                                            \
                                            for(;;)\
                                            {\
                                               imgDot((img), x0, y0);\
                                               if (x0 == x1 && y0 == y1) break;\
                                               e2 = 2 * err;\
                                               if (e2 >= dy) {err += dy; x0 += sx;}\
                                               if (e2 <= dx) {err += dx; y0 += sy;}\
                                            }\
                                         }

  #define imgSave(img, fname)            {\
                                            FILE* f = fopen((fname), "wb");\
                                            fprintf(f, "P6 %d %d 255\n", (img).width, (img).height);\
                                            fwrite((img).data, sizeof(Color), (img).width * (img).height, f);\
                                            fclose(f);\
                                         }

  #define imgLoad(img, fname)            {\
                                            FILE* f = fopen((fname), "rb");\
                                            char buffer[16];\
                                            int index = 0;\
                                            int field = 0;\
                                            int isP5 = 0;\
                                            unsigned char c = ' ';\
                                            while(field < 4)\
                                            {\
                                               do\
                                               {\
                                                  if(c == '#') while(c = fgetc(f), c != '\n');\
                                               } while(c = fgetc(f), isspace(c) || c == '#');\
                                               index = 0;\
                                               do\
                                               {\
                                                  buffer[index++] = c;\
                                               } while(c = fgetc(f), !isspace(c) && c != '#' && index < 16);\
                                               buffer[index] = 0;\
                                               switch(field)\
                                               {\
                                                  case 0:\
                                                     if (strcmp(buffer, "P5") == 0) isP5 = 1;\
                                                     else if (strcmp(buffer, "P6") == 0) isP5 = 0;\
                                                     else fprintf(stderr, "image format \"%s\" unsupported (not P5 or P6)\n", buffer), exit(1);\
                                                     break;\
                                                  case 1:\
                                                     (img).width = atoi(buffer);\
                                                     break;\
                                                  case 2:\
                                                     (img).height = atoi(buffer);\
                                                     break;\
                                                  case 3:\
                                                     index = atoi(buffer);\
                                                     if (index != 255) fprintf(stderr, "image format unsupported (not 255 values per channel)\n"), exit(1);\
                                                     break;\
                                               }\
                                               field++;\
                                            }\
                                            imgCreate((img), (img).width, (img).height);\
                                            if (isP5)\
                                            {\
                                               int length = (img).width * (img).height;\
                                               for(index = 0; index < length; ++index)\
                                               {\
                                                  (img).data[index].r = (img).data[index].g = (img).data[index].b = fgetc(f);\
                                               }\
                                            }\
                                            else\
                                            {\
                                               fread((img).data, sizeof(Color), (img).width * (img).height, f);\
                                            }\
                                            fclose(f);\
                                         }

Comments

খুব সুন্দর! আমি এটা পছন্দ করি.

— ঘুরে বেড়াও

Answer 7

আর: জিজিপ্লট 2 এবং মানচিত্রের অভিক্ষেপ ব্যবহার করা হচ্ছে

@ এমনিপের পোস্ট দ্বারা অনুপ্রাণিত হয়ে আমি আর এর ম্যাপ্রোজ প্যাকেজটি ব্যবহার করার চেষ্টা করার সিদ্ধান্ত নিয়েছি, যেখানে আমরা প্রজেকশন গণনা করার সময় উত্তর মেরু কোথায় থাকতে হবে তা নির্দিষ্ট করে বিশ্বকে আলোকিত করতে পারি।

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

কোডটির বেশিরভাগটি কেবল নান্দনিকতা। আমি যে অংশটি বের করতে পেরেছিলাম তা হল "দুপুর দ্রাঘিমাংশ" গণনা করা, এটি দ্রাঘিমাংশ মান যেখানে এটি কোনও GMT সময় দেওয়া দুপুর ছিল was

library(ggplot2);library(maps);library(ggmap)
world <- map_data("world")# a lat-long dataframe from the maps package
worldmap <- ggplot(world, aes(x=long, y=lat, group=group)) + 
  geom_path(color="orange") + 
  theme(panel.background= element_rect("black"),  
        axis.text.y=element_blank(),
        axis.ticks=element_blank(),
        axis.title.x=element_blank(),
        axis.title.y=element_blank(),
        panel.grid.major = element_line(colour="blue", size=0.75),
        panel.grid.minor = element_line(colour="blue")
  )  

#Create a function that takes in the current GMT time
print_3d_coordmap <- function (current_gmt_time) {
  curr_gmt_mins <- as.POSIXlt(current_gmt_time)$hour*60 + as.POSIXlt(current_gmt_time)$min
  noon_longitude <- 180 - (curr_gmt_mins * 360/1440)
  #centered at wherever longitude where it is Noon now on (lat:equator)  
  worldmap + coord_map("ortho", orientation=c(0, noon_longitude, 0))
}

#test it out
print_3d_coordmap(Sys.time() + 7*60*60) # my location is 7 hours behind UTC

এরপরে আমি 24 টি চিত্র উত্পন্ন করতে আর অ্যানিমেশন প্যাকেজটি ব্যবহার করেছি এবং সেগুলিকে একটি জিআইএফ এ সেলাই করেছি।

Comments

দুর্দান্ত লাগছে! যদিও এটি গ্রীষ্ম এবং শীতকে সঠিকভাবে পরিচালনা করে? আমি আর এটি খুব ভাল জানি না তবে মনে হচ্ছে আপনার নিরক্ষীয় অংশটি সর্বদা চিত্রের কেন্দ্রে থাকে।

— নওতা

হ্যাঁ আপনি সঠিক. আমি একমাত্র অক্ষাংশ হিসাবে নিরক্ষীয় অঞ্চল ব্যবহার করে একটি দ্রুত এবং নোংরা বাস্তবায়ন করেছি। (পরিবর্তে অ্যানিমেশনটিতে কেন্দ্রীভূত।) প্রকল্পটিতে আমি ব্যবহার করছি না এমন প্রচুর বৈশিষ্ট্য রয়েছে। যদি এমন কোনও রেফারেন্স থাকে যা মরসুমের সাথে অক্ষাংশকে কীভাবে আলাদা করতে পারে তা দেখায়, আমি এটি চেষ্টা করে খুশি হব।

— রাম নরসিমহান

আসলে এই আমি কোথায় থেকে আমার এন্ট্রিতে সূত্র অনূদিত - এখানে আর, রাম একটি রেফারেন্স এর stackoverflow.com/questions/8708048/...

— bazzargh

@ বাজজার্গ ধন্যবাদ! দেখে মনে হচ্ছে আমার করণীয় অজিমুথ কোণগুলি সম্পর্কে অনেক কিছু রয়েছে।

— রাম নরসিমহান

Answer 8

জাভাস্ক্রিপ্ট - মার্টিন ক্লেপ্পে ( http://aem1k.com/ )

আমি জোর দিয়ে বলতে চাই যে এটি আমার কাজ নয়, মার্টিন ক্লেপ্পের কাজ। আমি কেবল মনে করি এটি এতটা নিখুঁতভাবে ফিট করে যে এটি এখানে নিখোঁজ হওয়া উচিত নয়:

অনলাইন ডেমো (বা এটি কেবল কনসোলে আটকান)

eval(z='p="<"+"pre>"/*        ######## */;for(y in n="zw24l6k\
4e3t4jnt4qj24xh2 x/*    *############### */42kty24wrt413n243n\
9h243pdxt41csb yz/*  #################### */43iyb6k43pk7243nm\
r24".split(4)){/*     *#################*   */for(a in t=pars\
eInt(n[y],36)+/*          ###############*   */(e=x=r=[]))for\
(r=!r,i=0;t[a/*               ############*   */]>i;i+=.05)wi\
th(Math)x-= /*                #############    */.05,0<cos(o=\
new Date/1e3/*                #########*       */-x/PI)&&(e[~\
~(32*sin(o)*/*                     ####*       */sin(.5+y/7))\
+60] =-~ r);/*                         *###    */for(x=0;122>\
x;)p+="   *#"/*                        #####  */[e[x++]+e[x++\
]]||(S=("eval"/*                      *##### */+"(z=\'"+z.spl\
it(B = "\\\\")./*      ###*           ####  */join(B+B).split\
(Q="\'").join(B+Q/*                  ###* */)+Q+")//m1k")[x/2\
+61*y-1]).fontcolor/*               ##   */(/\\w/.test(S)&&"#\
03B");document.body.innerHTML=p+=B+"\\n"}setTimeout(z)')//m1k\

Comments

যদি এটি আপনার কাজ না হয় তবে আপনার উত্তরটি একটি সম্প্রদায়ের উইকিতে পরিণত করা উচিত।

— কাইল কানোজ

টিপটির জন্য ধন্যবাদ, আগে কখনও এই চেকবক্সটি খেয়াল করেনি। সম্পন্ন!

— ইনগো বার্ক