annotate src/laser/rasterize.clj @ 8:0f48db8d2a05

going to upgrade to clojure 1.2
author Robert McIntyre <rlm@mit.edu>
date Sat, 21 Aug 2010 03:21:22 -0400
parents 5e167f275a93
children f952052e37b7
rev   line source
rlm@0 1 (ns laser.rasterize)
rlm@0 2
rlm@0 3 (import '(java.io File))
rlm@0 4 (import '(org.apache.commons.io FileUtils))
rlm@0 5 (import '(javax.imageio ImageIO) )
rlm@0 6 (import '(javax.swing JFrame))
rlm@0 7 (import '(java.awt Color BorderLayout))
rlm@0 8 (import '(ij ImagePlus IJ))
rlm@0 9 (import '(java.lang Math))
rlm@2 10 (import '(java.awt Graphics2D Panel))
rlm@0 11 (import '(ij Macro))
rlm@0 12
rlm@0 13 (import '(java.io BufferedReader InputStreamReader))
rlm@0 14 (import '(java.awt.image BufferedImage))
rlm@0 15
rlm@0 16 (use 'clojure.contrib.str-utils)
rlm@8 17 (use 'clojure.contrib.seq-utils)
rlm@0 18 (use 'clojure.contrib.combinatorics)
rlm@8 19 (use 'clojure.contrib.duck-streams)
rlm@0 20
rlm@0 21 (use 'clojure.contrib.repl-utils)
rlm@0 22
rlm@0 23 (set! *print-length* 20)
rlm@0 24
rlm@0 25
rlm@0 26
rlm@1 27
rlm@1 28 (def feed 120)
rlm@1 29 (def dpi [500, 500])
rlm@8 30
rlm@8 31
rlm@1 32
rlm@1 33
rlm@5 34 (defn preserve-meta [f]
rlm@5 35 (fn [& x] (with-meta
rlm@5 36 (apply f x)
rlm@5 37 (meta (last x)))))
rlm@2 38
rlm@8 39 (defmulti frame-hash-multi class)
rlm@1 40
rlm@1 41
rlm@8 42 (defmethod frame-hash-multi ImagePlus
rlm@8 43 [image+]
rlm@0 44 (with-meta
rlm@0 45 (let [buf (.. image+ getBufferedImage)
rlm@0 46 color (.getColorModel buf)]
rlm@0 47 (apply hash-map
rlm@0 48 (interleave
rlm@0 49 (doall (for [x (range (.getWidth image+)) y (range (.getHeight image+))]
rlm@0 50 (vector x y)))
rlm@0 51 (doall (for [x (range (.getWidth image+)) y (range (.getHeight image+))]
rlm@0 52 (let [data (.getRGB buf x y)]
rlm@0 53 (hash-map :r (bit-shift-right (bit-and 0xff0000 data) 16)
rlm@0 54 :g (bit-shift-right (bit-and 0x00ff00 data) 8)
rlm@0 55 :b (bit-and 0x0000ff data))))))))
rlm@8 56 {:width (.getWidth image+) :height (.getHeight image+)}))
rlm@0 57
rlm@0 58
rlm@8 59 (defmethod frame-hash-multi String
rlm@8 60 [image-name]
rlm@8 61 (let [image+ (ImagePlus. image-name)]
rlm@8 62 (frame-hash-multi image+)))
rlm@5 63
rlm@5 64
rlm@8 65 (defn frame-hash
rlm@8 66 "yields a convienent representation for the pixles in an image.
rlm@8 67 Because of the size of the structvre generated, this must only be used
rlm@8 68 in a transient way so that java can do it's garbage collection."
rlm@8 69 [something]
rlm@8 70 (frame-hash-multi something))
rlm@8 71
rlm@8 72 ;(def frame-hash (preserve-meta frame-hash))
rlm@8 73
rlm@5 74
rlm@5 75
rlm@5 76
rlm@1 77 (def white {:r 255, :g 255, :b 255})
rlm@1 78 (def black {:r 0, :g 0, :b 0})
rlm@1 79
rlm@1 80 (def expt #(Math/pow %1 %2))
rlm@1 81
rlm@1 82 (defn rgb-euclidian
rlm@1 83 [{r1 :r g1 :g b1 :b} {r2 :r g2 :g b2 :b} ]
rlm@1 84 (expt (+ (expt (- r1 r2) 2)
rlm@1 85 (expt (- g1 g2) 2)
rlm@1 86 (expt (- b1 b2) 2)) 0.5))
rlm@1 87
rlm@1 88 (defn b&w
rlm@1 89 "turn everything strictly black or white"
rlm@1 90 [window]
rlm@1 91 (with-meta
rlm@1 92 (zipmap
rlm@1 93 (keys window)
rlm@1 94 (map (fn [rgb]
rlm@1 95 (if (> (rgb-euclidian rgb white) (rgb-euclidian rgb black))
rlm@1 96 black white))
rlm@1 97 (vals window))) (meta window)))
rlm@1 98
rlm@1 99
rlm@1 100
rlm@3 101 (defn raster-preamble []
rlm@3 102 (str-join \newline
rlm@3 103 ["M63 P0\nG61"
rlm@3 104 (str \F feed)
rlm@3 105 "M101"
rlm@8 106 "M3 S1\n"]))
rlm@3 107
rlm@4 108 (defn raster-epilogue []
rlm@4 109 (str-join \newline
rlm@8 110 ["M63 P0"
rlm@4 111 "M5"
rlm@8 112 "M2\n"]))
rlm@3 113
rlm@1 114
rlm@4 115 (defn raster-comment [string]
rlm@4 116 (str "(" (re-gsub #"[()]" "" string) ")"))
rlm@1 117
rlm@4 118 (defn filter-keys [fun m]
rlm@4 119 (select-keys m (filter fun (keys m))))
rlm@6 120
rlm@5 121 (def filter-keys (preserve-meta filter-keys))
rlm@2 122
rlm@4 123 (defn filter-vals [fun m]
rlm@5 124 (into {} (filter (comp fun val) m)))
rlm@6 125
rlm@5 126 (def filter-vals (preserve-meta filter-vals))
rlm@1 127
rlm@0 128 (defn frame-hash->bufferedImage
rlm@0 129 [frame-hash]
rlm@0 130 (let [data (meta frame-hash)
rlm@0 131 image (BufferedImage. (:width data) (:height data) BufferedImage/TYPE_INT_BGR)]
rlm@0 132
rlm@0 133 (doall (for [element frame-hash]
rlm@0 134 (let [coord (key element)
rlm@0 135 rgb (val element)
rlm@0 136 packed-RGB
rlm@0 137 (+ (bit-shift-left (:r rgb) 16)
rlm@0 138 (bit-shift-left (:g rgb) 8)
rlm@0 139 (:b rgb))]
rlm@0 140 (.setRGB image (first coord) (last coord) packed-RGB))))
rlm@0 141 image))
rlm@5 142
rlm@5 143 (defmulti display "Creates a JFrame and displays a buffered image" class)
rlm@0 144
rlm@5 145 (defn- makePanel [image] (proxy [Panel] [] (paint [g] (.drawImage g image 0 0 nil))))
rlm@0 146
rlm@6 147 (defn select-row [x window]
rlm@6 148 (filter-keys (comp (partial = x) first) window))
rlm@6 149
rlm@6 150
rlm@0 151
rlm@5 152 (defmethod display
rlm@5 153 BufferedImage [image]
rlm@5 154 (let [panel (makePanel image)
rlm@5 155 frame (JFrame. "Oh Yeah!")]
rlm@5 156 (.add frame panel)
rlm@5 157 (.pack frame)
rlm@5 158 (.setVisible frame true )
rlm@5 159 (.setSize frame(.getWidth image) (.getHeight image))))
rlm@5 160
rlm@5 161 (defmethod display
rlm@5 162 ImagePlus [image]
rlm@5 163 (display (.getBufferedImage image)))
rlm@5 164
rlm@5 165 (defmethod display
rlm@5 166 clojure.lang.PersistentHashMap [frame-hash]
rlm@5 167 (display (frame-hash->bufferedImage frame-hash)))
rlm@5 168
rlm@5 169 (defmethod display
rlm@5 170 clojure.lang.PersistentArrayMap [frame-hash]
rlm@5 171 (display (frame-hash->bufferedImage frame-hash)))
rlm@5 172
rlm@5 173
rlm@5 174
rlm@5 175
rlm@5 176
rlm@7 177 ;this is a sequence of rows
rlm@7 178
rlm@7 179 ;(defn span [row]
rlm@7 180 ; (let [sorted-row (sort #(< (first %1) (first %2)) row)]
rlm@7 181 ; (vector (first sorted-row) (last sorted-row))))
rlm@7 182
rlm@7 183
rlm@7 184 (defn row->gcode [[x-dpi y-dpi] row]
rlm@7 185 (let [[x1 y1] (first row)
rlm@7 186 [x2 y2] (last row)]
rlm@8 187 (str (format "G0 X%.3f Y%.3f\n"
rlm@7 188 (float (* x1 (/ x-dpi)))
rlm@7 189 (float (* y1 (/ y-dpi))))
rlm@7 190
rlm@7 191 (format "G1 X%.3f Y%.3f\n"
rlm@7 192 (float (* x2 (/ x-dpi)))
rlm@7 193 (float (* y2 (/ y-dpi)))))))
rlm@7 194
rlm@7 195 (defn gather-row [row]
rlm@7 196 (let [base [[(first (first row)) (first (first row))]]]
rlm@7 197 ; (println base)
rlm@7 198 (reduce
rlm@7 199 (fn colapse [collection new-n]
rlm@7 200
rlm@7 201 (let [collection (apply vector collection)
rlm@7 202 prevoius (last (last collection))
rlm@7 203 range-start (first (last collection))]
rlm@7 204 ; (println new-n)
rlm@7 205 ; (println prevoius)
rlm@7 206 ; (println range-start)
rlm@7 207 (if (<= new-n (+ prevoius 1))
rlm@7 208 (do ;(println "join")
rlm@7 209 ;(println (butlast collection))
rlm@7 210 (conj (apply vector (butlast collection)) (vector range-start new-n)))
rlm@7 211 (conj collection (vector new-n new-n)))))
rlm@7 212
rlm@7 213 base
rlm@7 214 (map first row))))
rlm@7 215
rlm@7 216
rlm@7 217
rlm@7 218 (defn row->gmask [[x-dpi y-dpi] forward? row]
rlm@7 219
rlm@7 220 (let [start (float (* (/ x-dpi) (first (first row))))]
rlm@7 221
rlm@7 222 (let [preamble (if forward?
rlm@7 223 (str
rlm@7 224 (format "0 0 0 %.3f\n" start)
rlm@7 225 (format "0 0 1 %.3f\n" start)))
rlm@7 226 body
rlm@7 227 (for [[x y] (gather-row row)]
rlm@7 228 (let [x (float (* x (/ x-dpi)))
rlm@7 229 y (float (* y (/ x-dpi)))]
rlm@7 230 (if forward?
rlm@7 231 (str (format "0 0 1 %.3f\n" x)
rlm@7 232 (format "0 1 1 %.3f\n" y))
rlm@7 233
rlm@7 234 (str (format "0 0 0 %.3f\n" x)
rlm@7 235 (format "0 1 0 %.3f\n" y)))))]
rlm@7 236
rlm@7 237 (str preamble (str-join "" body)))))
rlm@7 238
rlm@7 239
rlm@7 240
rlm@7 241 (defn make-rows [pic]
rlm@7 242
rlm@7 243 (map (partial sort #(< (first %1) (first %2)))
rlm@7 244 (partition-by last
rlm@7 245 (sort (fn [[x1 y1][x2 y2]] (> y2 y1))
rlm@7 246 (map first (filter-vals (partial = black) pic))))))
rlm@7 247
rlm@7 248
rlm@7 249
rlm@7 250 (defn generate-gmask [pic]
rlm@7 251
rlm@7 252 (str "1 0 0 0\n"
rlm@8 253 (str-join "" (map (fn [[index row]]
rlm@8 254 (row->gmask dpi (odd? index) row))
rlm@8 255 (indexed (make-rows pic))))))
rlm@8 256
rlm@8 257 ;; 1 0 0 0
rlm@8 258 ;; 0 0 1 2.881
rlm@8 259 ;; 0 0 0 2.881
rlm@8 260 ;; 0 1 0 2.863
rlm@8 261 ;; 0 0 0 2.769
rlm@8 262 ;; 0 1 0 2.751
rlm@8 263 ;; 0 0 0 2.729
rlm@8 264 ;; 0 1 0 2.617
rlm@8 265 ;; 0 0 0 2.593
rlm@8 266 ;; 0 1 0 2.561
rlm@8 267 ;; 0 0 0 2.463
rlm@8 268 ;; 0 1 0 2.445
rlm@8 269 ;; 0 0 0 2.385
rlm@8 270 ;; 0 1 0 2.317
rlm@8 271 ;; 0 0 0 2.253
rlm@8 272 ;; 0 1 0 2.233
rlm@8 273 ;; 0 0 0 2.177
rlm@8 274
rlm@7 275
rlm@7 276
rlm@7 277 (defn generate-gcode [pic]
rlm@8 278 (str (raster-preamble)
rlm@8 279
rlm@8 280
rlm@8 281 (str-join "" (map (partial row->gcode dpi) (make-rows pic)))
rlm@8 282 (raster-epilogue)))
rlm@8 283
rlm@7 284
rlm@7 285
rlm@8 286 (defn rotate [degrees #^ImagePlus image]
rlm@8 287 (.rotate (.getChannelProcessor image) degrees)
rlm@8 288 image)
rlm@7 289
rlm@8 290 (defn map-keys [f m]
rlm@8 291 (into {} (map (fn [[key val]] [(f key) val]) m)))
rlm@8 292
rlm@8 293
rlm@8 294
rlm@8 295 (defn invert-frame-hash [pic]
rlm@8 296 (map-keys (fn [[x y]] [x (- (:height (meta pic)) y 1)]) pic ))
rlm@8 297
rlm@8 298
rlm@8 299 (defn generate-files [pic]
rlm@8 300 (let [image (invert-frame-hash (b&w (frame-hash (rotate 180 (ImagePlus. pic)))))]
rlm@8 301 (spit "/home/r/kevin/out.ngc" (generate-gcode image))
rlm@8 302 (spit "/home/r/kevin/out.gmask" (generate-gmask image))
rlm@8 303 image))
rlm@8 304
rlm@8 305
rlm@8 306
rlm@8 307 ;;;; testing
rlm@8 308
rlm@8 309 (defn init []
rlm@8 310 (let [stuff
rlm@8 311
rlm@8 312 (bound-fn []
rlm@8 313
rlm@8 314 (do
rlm@8 315 (println "hi everyone")
rlm@8 316 (def img "/home/r/kevin/sing.png")
rlm@8 317 (def pic (frame-hash (let [image (ImagePlus. img)]
rlm@8 318 (.rotate (.getChannelProcessor image) 180)
rlm@8 319 image)))
rlm@8 320
rlm@8 321
rlm@8 322 (def test-image
rlm@8 323 (invert-frame-hash (b&w (frame-hash (rotate 180 (ImagePlus. img))))))
rlm@8 324
rlm@8 325 (defn test-gmask []
rlm@8 326 (println (str-join "" (take 170 (generate-gmask test-image)))))
rlm@8 327
rlm@8 328 (println "ALL variables initialized!")
rlm@8 329
rlm@8 330 ))]
rlm@8 331 (.start
rlm@8 332 (Thread. stuff))))
rlm@8 333
rlm@8 334
rlm@8 335
rlm@8 336 (defn thread-test []
rlm@8 337
rlm@8 338 (let [temp *out*]
rlm@8 339 (.start
rlm@8 340 (Thread.
rlm@8 341 (fn []
rlm@8 342 (with-bindings {#'clojure.core/*out* temp}
rlm@8 343 (Thread/sleep 5000)
rlm@8 344 (println "hi")))))))